Finite Element Analysis of the Effect of Dental Implants on Jaw Bone under Mechanical and Thermal Loading Conditions

Dental implants have been studied over the years to replace missing teeth. One of the conditions for the success of implants is their stability and resistance under the applied forces and minimal tension in the surrounding bone. )e purpose of this dissertation is numerical and three-dimensional analysis of jaws with implants under mechanical and thermal loading by the finite element method. For this purpose, implant simulations (including ceramic crown, titanium root, and jaw bone) under dynamic and thermal load have been performed in Abacus software. In this simulation, it is considered that the jawbone is composed of two areas, one area is the superficial bone tissue (cortical) and the other part is the spongy tissue. Implants are usually made of different metals or ceramics with a bone-like structure that are compatible with body tissues. Implants are currently made of titanium metal. )erefore, titanium metal has been used for modeling implants in this dissertation. )e implant crown is also considered as a ceramic material. In the simulation, the effect of stresses imposed by the implant on the jawbone is performed. In this simulation, mechanical force is applied to the upper part of the implant and force enters the jawbone through the implant, which causes tension at the junction of the implant to the jawbone. To investigate the effect of thermal loads, different temperature conditions are considered by considering the decrease in temperature and increase in temperature on the tooth surface and its effect on the implant and the jaw bone. After validation and ensuring the accuracy of the modeling, it has been observed that, with increasing mechanical load, the stresses created in all parts of the ceramic coating, titanium implants, and jawbone have increased. It is also observed that the stress created in the titanium implant due to the application of negative heat flux was about twice as much as the stress created due to the application of positive heat flux.


Introduction
Given the importance of the human body and the importance and impact of the type and properties of implants in the jaw, as well as its impact on human body function, choosing the best parameters both in terms of mechanical properties and biocompatibility (compatibility with the human body) is essential. Also, the function of dental implants and its effect on the jaw bone is affected by mechanical and thermal loads, which is also important and necessary to study this issue. erefore, by modeling and analyzing the effect of dental implants on the jaw bone in terms of mechanical and thermal loading, as well as considering different properties for the implant, it is possible to have a detailed study on its performance, which is very important [1][2][3].
We may lose one or more of our teeth for a variety of reasons. Not having one or more teeth can cause problems such as loss of beauty.
ere are several alternatives to missing teeth, among which dental implants are the best choice. However, various reasons, such as the inability to withstand the forces and torques and the intolerance of high temperatures, can lead to failure in the implantation process. A dental implant is a titanium screw that replaces a missing tooth by being inserted into the jawbone. Dental implants can be one-piece or two-piece, with the crown and root joined by a screw. Dental implants are an alternative to shaven-tooth bridge prostheses, and the primary goal of introducing this therapy (dental implant) is to ensure that no teeth will be shaved or damaged in the future owing to the bridge base.
Oguz et al. [4] studied the static, dynamic, and fatigue behavior of dental implants using the finite element method. In this study, a dynamic load is applied to the occlusal surface for 5 seconds and the fatigue life is calculated based on Goodman, Soderbergh, and Gerber criteria. ey found that von Mises maximum stress in dynamic loading is greater than static loading. Kong et al. [5] investigated the effect of thread change in maximum stress on bone and implants through finite element analysis. e results showed that the thread pitch plays an important role in the strength of the implant under axial load. Djebbae et al. and Tian et al. [6,7] obtained the stress distribution in dentures by the finite element method. e amount of stress, especially at the junction of the implant bone, was examined. ey found that force and direction of loading have a large effect on the amount of stress. Huang et al. and Dorogoy et al. [8,9] used 11 different finite element models to study the stress distribution and slip of implants. e contact of the surfaces is of frictional type. e results show that, in fastloading implants, especially in off-axis loads, the stress on the bone around the implant is very high. As the number of threads and the contact surface of the bone and implant increase, the stress distribution and slip of the implant and bone surface decrease. Guan et al. and Yazdi [10,11] performed dynamic modeling and simulation of the dental implant placement process using the finite element method. In this study, by examining the effect of different placement depths on spongy bone and dense bone, it was found that increasing the placement depth increases the amount of stress in dense bone. e authors of [12][13][14] presented finite element analysis of thermal implants exposed to heat. In this study, by observing the results of thermal stresses, it can be seen that thermal stresses have small values, and it is also due to small temperature changes in the complex. e highest stress is related to the abutment, which is also caused by the stiffening torque. e authors in [15,16] reviewed the mechanical design requirements of dental implants and found that compressive loading stabilizes the implant in the jawbone and tensile loading of the implant loosens and shear loading causes the implant to fail. Creating rough surfaces also improves the adhesion between the implant and the bone. Using the finite element method, the authors of [17,18] studied the influence of the tooth and jaw bone-implant contact model on maximum implant stress. e largest stress can occur in the neck area of the implant and in the connection points with the dense bone under oblique loading, according to the stress distribution. e authors of [19,20] used the finite element method to undertake static and dynamic analysis of dental implants that Abacus specialist software was used in this investigation. Modeling geometry, specifying material characteristics, boundary conditions, contact conditions, loading, and elementing for bone-implant simulation were all done in Abacus program. ere will be two analyses, each with a distinct force. e ceramic cap is subjected to a force of 100 Newtons at a 45-degree angle during static analysis. All modeling phases are inserted into dynamic analysis with a force of 400 Newtons in 0.01 seconds, identical to static analysis. In static analysis, the highest von Mises stress was 82.5 MPa, but in dynamic analysis, it was 770.2 MPa. In contrast to dynamic analysis, the maximum stress in static analysis is lower. Dynamic analysis allows for a smaller maximum displacement than static analysis. e authors of [21,22] used finite element simulation to evaluate the biomechanical behavior of mini-implants under real-world working settings. Stress analysis of two distinct quality of the D2 and D3 jawbone around three types of mini-implants was performed using the finite element method in this study due to the relevance of placing short implants in confined spaces between the edentulous region. ree varieties of Osteocare, Dio, and Dentis mini-implants were included in this experimental investigation. e highest component of the abutment was loaded with a vertical force of 100 N and a lateral force of 100 N at a 45-degree angle.
ABAQUS software was used to analyze stress levels in the mini-implant and surrounding bone. e findings revealed that the amount of von Mises stress in D3 bone for all implants is higher than D2 bone and that the level of stress in the cortical bone is higher than the spongy bone. Furthermore, all the systems studied had identical stress distributions in the cortical bone. Additionally, the initial implant thread caused the most stress in the implant's neck, but the Osteocare mini-implant created less tension in the bone. e topology optimization technique in dental implants was examined using Abacus software by [23][24][25]. is article discusses how to use Abacus software to optimize dental implant architecture. One of the types of optimization in Abacus software is topology optimization. Topological optimization is a mathematical method for determining the best material distribution shape for a structure in a given space. Optimization minimizes the weight of the material, resulting in lower assembly costs and time. e authors of [26][27][28] used the finite element approach to simulate the mechanical performance of dental implants comprised of memory materials. A sample of existing dental implants, as well as a portion of the jawbone, were modeled in Abacus software and statically assessed for this purpose. Finally, the implant's level of stress was compared to that of other implants made of standard materials. e results show that when the memory implant was utilized, the stress values in the implant were lower than when the nickel-titanium implant was used. e amount of stress transferred from the nitinol implant to the nickeltitanium implant to the jawbone was also fewer in the regions which were subjected to increased force and stress. e authors of [29,30] used the finite element approach to investigate the effect of geometric and mechanical features on the stress distribution of the dental implant system. e findings demonstrate that the angle and step of fastening are critical in boosting implant stability and minimizing bone stress.
Over the years, dental implants have been researched as a way to replace missing teeth. One of the requirements for implant success is their stability and resistance to applied stresses, as well as minimum tension in the surrounding bone.
e goal of this work is to use the finite element method to perform numerical and three-dimensional analysis of jaws with implants under mechanical and thermal loading. Implant simulation (containing ceramic crown, titanium root, and jaw bone) is performed in Abacus program under dynamic and thermal load for this purpose. e jawbone is divided into two sections in this simulation: the superficial bone tissue (cortical) and the spongy tissue. Implants are typically made of metals or ceramics that have a bone-like structure and are compatible with bodily tissues. Currently, titanium metal is used to make implants [31]. As a result, in this article, titanium metal is used to simulate the implant. A ceramic substance is also used to make the implant crown. e effect of implant-induced strains on the jawbone is simulated in this simulation. Mechanical force is given to the top section of the implant and force is applied to the jawbone through the implant in this simulation, resulting in strains and concentrations of stresses at the implant-jawbone junction. Different temperature settings are investigated to investigate the effect of thermal stresses, including temperature decreases and increases, on the tooth surface, as well as their effect on the implant and the jaw bone. High temperature tolerance can arise when drilling the jawbone or when drinking hot liquids, according to the findings of this study. As a result, the goal of this research is to use the finite element method to conduct a numerical and three-dimensional analysis of jaws with implants under mechanical and thermal loads. us, we derive the stress distribution, strain, and displacement in the implant and jawbone by modeling the set of implants and jawbone while considering mechanical and thermal loads, and we investigate the effect of various parameters on them.

Modelling Software.
e study approach involves modeling the geometric model of the jaw bone in Solidwork's software in two states: cortical and spongy. Solid-Works software is also used to construct the geometric model of the implant set and its cover, which is then inserted in the bone. e boundary and force properties and conditions are determined after importing the designed geometric models into the Abaqus analytical software, and the model is then evaluated and examined after meshing. Mechanical and thermal loads are applied to the model in this finite element study, and the heat load due to drinking hot and cold drinks is applied to the whole surface of the ceramic coating in the form of heat flux per unit area. e size of the element converges to improve the accuracy of the simulation findings. e outcomes of modeling this research are also compared to the results of valid articles for validation. e jaw and the implant are presumed to be in perfect contact in this study and that there is no slippage between them. In fact, the bone and the titanium implant are joined, to use medical terminology. e step is utilized directly or implicitly in the analysis, which is dynamic. e simulation approach is offered after the stages of modeling and specifying the required problem in this program are accomplished step by step in different modules of Abacus software.

Create a Geometric Model.
e first stage in issue software modeling is to develop a geometric model of the problem's pieces. e geometry of the ceramic crown, titanium root, and jawbone, as depicted in Figures 1-3, is initially modeled for this purpose. ese parts' modeling is three-dimensional and adaptable.

Definition of Material Properties.
e crown of the implant is constructed of ceramic, while the root is composed of titanium Ti6Al4V metal, and the jawbone is divided into two groups: cortical and spongy. As a result, mechanical properties such as density, Young's modulus, Poisson's ratio, and plasticity qualities, as well as thermal properties such as specific heat, are determined for each of the materials mentioned. e Drucker-Prager plastic model is combined with the ductile damage model to create the ceramic material for the implant crown. e metal root of the implant is likewise subjected to Johnson Cook's damage model. e jawbone has also benefited from Johnson Cook's elastic and plastic characteristics. Table 1 shows the mechanical properties of implant components.   Table 2 shows the mechanical properties of bone components. Table 3 shows the thermal properties of the model components.

Assembly of Parts.
Different pieces of the model, including as the ceramic crown, titanium root, and jawbone, are positioned adjacent to each other and their relative positions are defined for this reason. Figure 4 depicts the model parts being put together.

Define the Type of Analysis.
e analysis is classified as linked temp-displacement. is problem clearly considers the type of solver.
is solution examines the governing equation system using the element diagonal mass matrix and the law of explicit integration. Figure 5 depicts the software's definition of this type of analysis for massive transformations.

Define Boundary Conditions and Loadings.
e jawbone's lateral surfaces are considered joint support in this article. e upper surface of the ceramic crown has also been subjected to two forms of loading: mechanical and thermal.
e entire upper surface of the ceramic crown is paired with a reference point on the same surface in order to establish mechanical and thermal loading, and then, the intended loads are applied to this point. On the upper surface of the crown, a compressive force is exerted. ree distinct values for the amount of compressive force have been evaluated in various analyses to investigate the effect of the amount of force on the behavior of the dental implant. A heat flux is delivered to the upper surface of the ceramic crown to provide a thermal burden.
ree distinct values for the amount of heat flow have been investigated in various analyses to investigate the effect of heat flux on the behavior of dental implants. As a result, Abacus software is used to define the above boundary conditions and loads. e articular support on the lateral sides of the jawbone is defined in Figure 6. Figure 7 shows the application of mechanical load as a compressive force to the upper surface of the crown. e application of heat flux is seen in Figure 8. It is worth noting that the starting temperature and surrounding environment are both 27 degrees Celsius.

Meshing.
e mesh of the model parts in the form of hexagonal meshes with an element size of 0.2 mm is taken into account. Temperature-displacement coupling elements, explicit library, and quadratic geometric order are among the elements in its family, with a total of 1,986,000 elements. In Section 3, we will look into mesh independence for this amount of elements. e model set correlation is depicted in Figure 9.

Check the Independence of the Mesh.
e simulation results for numerous different element sizes were explored to ensure that the results from the elements and meshes were independent. For varying sizes of model elements, Table 4 indicates the maximum output stress values in the ceramic coating.  As can be observed, the results are nearly consistent from element sizes equal to 0.2 mm onwards, indicating that the analysis findings are reliable and that the solution has attained convergence with a high degree of accuracy. e findings are also shown to be independent of the number and size of elements and meshes when using this element size. As a result, various findings have been derived using this element's size.

Validation.
A comparison was done between the results of the current work and the results connected to Niroumand and Jafari in order to validate and evaluate the accuracy of the modeling and numerical solution method. Table 5 shows the maximal von Mises stress of the current investigation and compares it to the results of robust and extensive research for this goal. ese values are provided for the model's various components. e stresses in the table are measured in megapascals.
Because there is a small discrepancy between the results of this study and the results of Niroumand and Jafari's publication, it is concluded that the current modeling is accurate and legitimate.

Stress Results.
For three distinct loads, Figure 10 depicts the stress distribution in the ceramic coating area (2500 N, 5000N, and7500).   Mathematical Problems in Engineering e stress distribution in the titanium implant region for three distinct loads is shown in Figure 11 (2500 N, 5000 N, and 7500 N). Figure 12 depicts the distribution of stress in the bone under three different loads (2500 N, 5000N, and 7500 N). Figure 13 shows the strain distribution in the ceramic coating section for three different loads (2500 N, 5000N, and 7500). e strain distribution in the titanium implant part for three distinct loads is shown in Figure 14 (2500 N, 5000 N, and 7500 N). Figure 15 shows the strain distribution in the bone for three different loads (2500 N, 5000 N, and 7500 N). Figure 16 shows the deformation distribution in the ceramic coating section for three different loads (2500 N, 5000 N, and 7500).     Mathematical Problems in Engineering                 3.6. Investigating the Effect of Heat Flux. Table 6 shows the temperature values in different parts of the model, due to the application of negative and positive heat flux. Table 7 shows the maximum stress values caused by negative and positive heat flux in various regions of the model.

Conclusion
Over the years, dental implants have been researched as a way to replace missing teeth. One of the requirements for implant success is their stability and resistance to applied stresses, as well as minimum tension in the surrounding    e goal of this work is to use the finite element method to perform numerical and three-dimensional analysis of jaws with implants under mechanical and thermal loading. Implant simulations (containing ceramic crowns, titanium roots, and jawbone) under dynamic and thermal load were done in Abacus program for this purpose, and the following results were found after validation and assuring the reliability of modeling: (1) e strains induced in all components of the ceramic covering, titanium implants, and jawbone have increased as the mechanical load has grown. It is worth noting that the stress levels shown are related to the von Mises result stress, which may be calculated using the equation below: (2) e highest levels of stress were seen in the titanium implant at varied loads, whereas the lowest quantities were observed in the jawbone due to the stress distribution.
(4) It was discovered that, under different loads, the maximum values of strain occurred in the titanium implant and the lowest amount occurred in the ceramic covering due to the strain distribution.
(5) As the mechanical stress on the ceramic coating, titanium implants, and jawbone increases, deformations in all portions of the ceramic coating, titanium implants, and jawbone rise. It is worth noting that the values of the presented deformations are related to the result's deformation.
(6) e maximum amount of distortion occurred in the ceramic covering under varied weights, while the lowest amount occurred in the jawbone, owing to the deformation distribution. (7) It has been noticed that, due to the effects of positive and negative heat flux, the temperature of the high surface of the ceramic coating that is exposed to the current has increased to 48.5°C when positive heat flux is applied and has reduced to 23.48°C when negative heat flux is applied. Furthermore, due to positive heat flow, the temperature of the titanium implant increased to 42.73°C and reduced to 30.58°C due to negative heat flux. (8) When positive heat flux was applied to the jaw bone, the temperature of the spongy part of the bone rose to 39.30°C, and when negative heat flux was applied, the temperature of the spongy part of the bone rose to 34.65°C. Due to the positive heat flux, the temperature of the cortical section of the bone increased to 41.90°C and reduced to 31.58°C. (9) In thermal analysis of the implant and bone set, the maximum temperature transferred to the bone tissue through the implant is investigated. How much heat has reached the bone tissue and the contact surface of the implant with the bone? With some proteins denaturing at 42°C and above, as well as exposing the bone to 47°C for 1 minute, rising bone temperature can have fatal consequences. As a result of the findings of this study, it may be regulated that conditions of detrimental effects on bone tissue do not exist in patients with dental implants. (10) When comparing the maximum stress values in different regions of the model, it was discovered that the stress caused by negative and positive heat flux in the ceramic covering, sponge bone, and cortical bone is about equivalent. While the stress induced in the titanium implant by negative heat flow was roughly twice that created by positive heat flux. is is due to the temperature gradient in the two circumstances being different. In the negative flux mode, the temperature gradient in the implant was around 7°, while in the positive flux mode, it was about 5°.

Suggestions.
In this study, finite element analysis of the effect of dental implants on the jawbone under mechanical and thermal loading conditions and its results are presented. e following are suggestions for research in this area: (1) Considering different contact models in modeling the placement of implants on the jawbone (2) Examining different types of implants with different geometry and materials and comparing their performance (3) Optimizing the performance of dental implants according to various geometric and physical parameters Data Availability e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.