Abstract

The selection of materials has a considerable role on long-term stability of implants. The materials having high resistance to fatigue are required for dental implant applications since these implants are subjected to cyclic loads during chewing. This study evaluates the performance of different types of materials (AISI 316L stainless steel, alumina and its porous state, CoCr alloys, yttrium-stabilized zirconia (YSZ), zirconia-toughened alumina (ZTA), and cp Ti with the nanotubular TiO2 surface) by finite element analysis (FEA) under real cyclic biting loads and researches the optimum material for implant applications. For the analysis, the implant design generated by our group was utilized. The mechanical behavior and the life of the implant under biting loads were estimated based on the material and surface properties. According to the condition based on ISO 14801, the FEA results showed that the equivalent von Mises stress values were in the range of 226.95 MPa and 239.05 MPa. The penetration analysis was also performed, and the calculated penetration of the models onto the bone structure ranged between 0.0037389 mm and 0.013626 mm. L-605 CoCr alloy-assigned implant model showed the least penetration, while cp Ti with the nanotubular TiO2 surface led to the most one. However, the difference was about 0.01 mm, and it may not be evaluated as a distinct difference. As the final numerical evaluation item, the fatigue life was executed, and the results were achieved in the range of 4 × 105 and 1 × 109 cycles. These results indicated that different materials showed good performance for each evaluation component, but considering the overall mechanical performance and the treatment process (implant adsorption) by means of surface properties, cp Ti with the nanotubular TiO2 surface material was evaluated as the suitable one, and it may also be implied that it displayed enough performance in the designed dental implant model.

1. Introduction

The optimal dental implants require the integration of material, mechanical, biologic, and economic factors. Especially, the stability of the implant under cyclic biting loads, the integrity of the implant to the bone tissue, and stress transformation from the implant to the bone tissue are the most important concerns in the success of implantation. In order to make considerable progress in this technology, it is essential to determine the reasons underlying the failure of implants [14].

The long-term mechanical reliability of dental implants is one of the most important concerns, so the cyclic life occurs as a failure problem with them. In this type of failure, fracture mechanism time is assumed as the prime factor to be considered since the damage accumulates over time. When the deformation reaches a critical level, catastrophic failure takes place in dental implants [5]. Thus, the fatigue can be implied as the main problem affecting the cyclic life. Stress concentrations, which are due to the threads relatively in the structure of the implant and a significant degradation of the material surface containing flaws, determine the fatigue behavior and the life of the implant material [6]. For this reason, material selection was paid attention as to provide enough fatigue life in this study.

The fixation is a frequently observed problem in dental implants. The integration mechanism of the implant is based on the placement and then the quantity and quality of cortical and trabecular bones, which affect the primary stability and so the fixation of the implant [4]. Herein, the osseointegration, which consists of anchorage of metallic implants to the related bone, makes the implant-bone interaction site the most important issue. Regarding the interaction, the loss of bone density, which is occurred due to the significant reduction of the dimensions of the mandible (i.e., induced by aging) and the massive atrophic changes, becomes the crucial point [2, 7]. The reduction in size and density of the bone impairs the placement of implants [1]. To solve this problem, bone density and the other factors, for example, applied loads, the length and diameter of the implant, its geometry and surface, the bone-implant surface, and the quality and quantity of the surrounding bone are considered and researched for obtaining the optimization.

Bone resorption is another serious complication in the case of integration in dental implant applications. By modifying the surface properties of implants, bone healing may be accelerated and stimulation of bone growth can be provided. For example, the rougher surfaces stimulate differentiation, growth, and attachment of bone cells and increase mineralization [8]. Therefore, the term of roughness hereof becomes important, so certain implants are known as they are manufactured with rough surfaces purposely [9]. In recent years, the surface properties of the implants are searched to develop for increasing the adsorption to the bone structure. Thus, focusing on the porous materials and their comparison becomes one of the basis of this study.

The other main problem can be marked as inconsistency of the implant with the surrounding bone tissue. In the study of Elias et al., it is implied that rigidity inconsistency between the bone and the implant materials causes stress shielding [10]. This problem is remarked as a high difference between the moduli of elasticity of the bone and implant that induces insufficient transfer of stress to the bone due to the high modulus of the implant. To overcome the problem, they pointed that the materials with lower elastic modulus, which provides closer value to the one of the bone tissues, can be utilized, since they have better stress distribution at the implant-bone interface and lead to less bone atrophy. Additionally, considering the porous materials and adjusting porosity to match the elastic modulus of the human bone may be stated as another technique in the literature for annihilating the problem [1113].

In order to sort out all these mentioned problems, the implants were investigated in view of the material in this study. Therefore, different types of biomaterials were assigned as a dental material in the implant model developed by our group. As a dental material, AISI 316L stainless steel, alumina and its porous state, CoCr alloys, yttrium-stabilized zirconia (YSZ), zirconia-toughened alumina (ZTA), and cp Ti with the nanotubular TiO2 surface were tested. The stress distribution over the implant, bone-implant interaction, and the life and the stability of implants were investigated by simulations under real cyclic biting loads. The effects of materials on the performance of implants were also discussed, and the optimum material properties for dental implants under fatigue were determined.

2. Materials and Methods

2.1. CAD Model Preparation of the Implant and Abutment

First, a reference implant-abutment model was needed to set an FEA model, so considering the criterions, a design of the implant that was generated by SolidWorks 2012 (Dassault Systèmes, France) was used in this work. This developed implant model has a 2.0 mm height of the nonthreaded cylindrical neck and an 11.0 mm height of the threaded part. Its diameter is 4.5 mm and length is 13.0 mm. An implant model (AdVent) of a trademark that is Zimmer Dental, Inc. (CA, USA) was used as a base for the thread detail [14]. In Figure 1, technical drawing of the implant-abutment system is shown.

2.2. FEA Model Preparation of the Implant and Abutment

To assess the effect of different dental implant materials, the CAD data were prepared for the FEA, in which the condition was set according to the ISO 14801. The generated implant, the abutment, and also the implant holder were assembled in CAD software and then exported to the software ANSYS Workbench 16.2 (ANSYS, Inc., USA).

Mesh configuration was set by the method of “tetrahedrons” to generate 10-node tetrahedral elements. Body sizing option was used for whole bodies to generate finer meshes with the option of “sphere of influence.” The center of sphere was located at the critical region that is near to the second thread from the head of the implant, and the element size was taken as 0.6 mm.

Using the study of Ferreira et al., the frictional connection was defined with frictional contacts, in which the coefficient was taken as 0.3 [15]. The “Augmented Lagrange” option of contact solution formulation was set, and the “Adjust to Touch” option was also imposed for the interface treatment because of the assembled structure bodies that were in touch with each other.

The lateral surface nodes of the implant holder were fix supported for the boundary condition. In the literature, it was revealed that a 100 N load is in the range of normal bite forces. Therefore, in this study, 100 N was applied to the top surface of the abutment with an inclination of 30° according to the ISO 14801 standard [16, 17]. According to recommendation of the implant manufacturer, the tightening torque of 30 N·cm was applied to the abutment for demonstrating as possible as the real condition [14]. Figure 2 shows the brief FEA conditions.

2.3. Material Selection and Assignment to the Bodies of the Model

For the FE (finite element) solution of the model, the material properties had to be assigned to all of the bodies in the FE model properly, so the conventional and novel biomaterials were researched in detail according to the criteria as mentioned in Introduction of this study for the ability of being used in dental application. It is known that Ti alloys, AISI 316L stainless steel, CoCr alloys, alumina, zirconia-toughened alumina (ZTA), and yttrium-stabilized zirconia (YSZ) are still studied and asserted as the alternative materials for dental implantation [1821].

Herein, it would be better to make an elaboration on them and general comparisons between the applicable materials for supporting the reasons of the selection of them.

When the literature is searched, it can be observed that AISI 316L stainless steel is preferred and used in dental implant applications among the mentioned materials because of its relatively low cost, corrosion resistance, and ease of manufacturing [18]. Moreover, additional coating is also preferred for decreasing the ion accumulation in the interacted bone region due to the compound elements in the alloy.

The other alloy type is CoCr that can also be used for dental implant applications, and they are regarded as excellent in the corrosion resistance since they have adherent layer of chrome-based oxides on the surface and it leads to a passivating effect. Thus, the implant surface is less rough due to both before and after saliva exposure. The situation indicates that the active material degradation process is less than the other commercial alloys such as titanium [19]. However, manufacturing of this material may be a little bit difficult and take long time relatively, so its usage is limited.

On the other hand, usage of ceramic materials can also be asserted as rising up. The ceramic implants become popular and familiar with the patients and clinicians due to their significant wear resistance and superior biocompatibility relatively. One of the popular ceramic material structures is porous alumina that is known as biocompatible in the porous form and it permits tissue ingrowth that provides an important property of implant stabilization [21]. For that reason, the possible novel ceramic materials like porous alumina were investigated for availability of being used as dental implants in this study. In this context, yttria-stabilized tetragonal zirconia (YSZ) was also studied. It is known as another preferred one in dental implant applications due to its superior flexural strength that is in the interval of 800 and 1000 MPa, toughness, corrosion, and wear resistance [22, 23].

Moreover, the alumina material can also be strengthened and toughened by zirconia. Tetragonal-phase zirconia addition to alumina leads to higher strength and fracture toughness with little reduction in hardness and elastic modulus relatively [24, 25]. Kurtz et al. also imply that higher fracture toughness allows for the manufacture of thinner liners to reduce the risk of impingement and dislocation and improve the stability of the implants [25].

Besides these ceramic materials, when the literature is investigated, a common base material of Ti appears as either a metallic alloy or a ceramic structure that is derived by controlled oxidation. Metallic alloy states of Ti are known as already used in dental implants. They provide superior fatigue life and strength. However, they may have a possible toxic effect resulting from released vanadium and aluminum, so their commercially pure grades are preferred [26].

Furthermore, the implants are desired with having enough roughness for enhancing the integrity with trabecular bone tissue, so it can be implied that this demand may be met by the manufacturing stage or material properties. So, considering the material properties and their effects, the novel materials with the surfaces of nanoporous anodic alumina and nanotubular titanium oxide were also studied for comparison with the other materials that were mentioned previously. They were assigned as the materials with their mechanical properties, which are listed in Table 1, in the FEA solver. Additionally, the fixture was assigned as the trabecular bone to observe its behavior after the loading.

2.4. Fatigue Analysis Procedure

Besides the static analysis, a well dental implant must perform enough fatigue life [6]. Fatigue FEA is utilized generally for determination of the cyclic performance of materials and designs. Thus, the fatigue analysis of the implant model was run by using the related material properties in the fatigue tool of ANSYS Workbench 16.2. To perform the analysis, the related materials’ S-N curves data were imposed to the solver [37, 4043]. In this process, there are exemption data about the cyclic properties of the porous and tubular structures. Since they are generally evaluated by the properties of the main or base material, the S-N data were selected according to the bulk cp Ti and alumina materials for the structures of nanotubular and nanoporous. Additionally, fatigue life analysis was solved by the Goodman fatigue theory that considers the tensile strength values of the materials since it is preferred for the brittle materials especially [44].

3. Results and Discussion

According to the static FEA results, the equivalent stress and also the penetration onto the bone structure were calculated and discussed. Then, by using the equivalent stress results, the fatigue life calculations were executed and they were also discussed.

3.1. Equivalent (von Mises) Stress Values

First of all, the equivalent (von Mises) stress values and their concentrated regions are shown in Figure 3.

It can be obviously seen that the maximum von Mises stresses occurred on the first thread for all the models. Thus, the critical region can also be implied as the first knuckle. It is consistent with the results of the investigation of Wang et al. [45]. Hence, such implants may be predicted to be fractured or have a failure on the related thread region, so it is related directly to the design; however, the design issue is not an objective of this study. It is not a parameter in those FEA because the design of all the models is the same. When the quantitative evaluation is made, the calculated equivalent stress values are observed in the range of 226.95 MPa and 239.05 MPa. They are close to each other relatively since the equivalent stress is dependent on dimensions of the implant, loading types and magnitudes, elastic modulus, and Poisson’s ratio of materials.

The same implant dimensions and loading conditions. However, when the FEM theory is considered, it is known that results are generally evaluated in the elastic behavior for static analysis and the stress values are generated according to the basic strength calculations and theories. They are also in relation to deformation amount, and the couple terms compensate with each other. Herein, the results should be discussed in this approach, and the different but close values may be clarified with this compensation, so the deformation values may be implied as much different for the close stress-induced models. Additionally, when Poisson’s ratio values are considered, the model assigned with YSZ may be interpreted as tending to more transversal deformation than the other ones because its Poisson’s ratio is the highest one among the others. Thus, the first knuckle of this material’s model may probably have the most deformation, so the most equivalent stress value is observed on it.

3.2. Penetration onto the Bone

Separately from the deformation, the main constraints are the stress state and cyclic life and also the conformity of the rigidity of implant materials with the bone for this study. In order to provide better stress transfer from the implant to the bone and avoid the stress-shielding effect, the elastic modulus of the implant is desired to be as close as possible to that of the bone [46]. By this fact, the materials of cp Ti with the nanotubular TiO2 surface, YSZ, and AISI 316L stainless steel become prominent since their elastic modulus is relatively low. Also, the lowest elastic modulus value can be seen in nanotubular titania structure. Among the studied materials, this value is the closest one to the value of bone. Moreover, this material caused the highest von Mises stress value on the implant model because its low elastic modulus provides more deformation on the structure. Thus, this prevents the stress shielding and also provides smooth stress transfer instead of the severe stress transfer to the bone. At this stage, the proximateness of Young’s modulus becomes important. Here, it can be possible to mention about existing different methods to calculate Young’s modulus of cp Ti with the nanotubular TiO2 surface. One of them is calculation by an equation based on the Berkovich indentation technique consisting of contact stiffness, maximum load, and depth [47]. The other one is simpler but useful to have a close value. It is based on the rule of mixtures technique that is used for continuous fiber-reinforced composites generally [48]. It considers the structure as in the isostrain condition or the isostress condition, so the orientation determines Young’s modulus of the composite. Then, it can be implied that Young’s modulus of the cp Ti with the nanotubular TiO2 surface is variable according to the fraction of the nanotubes and the direction. When the adjustment of its tubular structure is considered, Young’s modulus can be adjusted. Thus, it may be thought that nanotubular titania tends to absorb the stress and deform much more relatively. At this point, studying the penetration amounts of the implant models onto the bone becomes appropriate. The bone penetration amounts are shown in Figure 4.

The highest amount of penetration onto the bone structure was calculated on the implant model, in which the material was assigned by cp Ti with nanotubular TiO2 surface, as 0.013626 mm. The lowest value was also calculated on the implant model, in which the material was assigned by L-605 CoCr alloy, as 0.0037389 mm. In view of microscale, the highest result can be assumed as a much penetration value. However, the difference between these extremum penetration amounts is 0.01 mm approximately; hence, just the penetration amount is solely not sufficient for the evaluation. Therefore, However, the surface roughness issue and thus the bone adsorption get importance at this point. Coelho et al. imply that different methods of implant surface engineering may lead to different and unique surface properties, which might affect the host-to-implant response, and state that the new implant surfaces should be tested and evaluated as new biomaterials [9]. By the principle of the bone healing and growth that can be accelerated and provided by rougher surfaces, the results seem that the materials with the structures of nanotubular and nanoporous may be applicable for manufacturing the implant. However, when their results are evaluated with their unique mechanical properties, the magnitude of the von Mises stress of the model with nanoporous anodic alumina is 237.56 MPa and exceeds its ultimate tensile strength value, and the nanoporous region will probably get damaged, so this material should be studied and developed much more. Thus, the nanotubular titania becomes prominent again according to another fact.

3.3. Fatigue Life Results

The life prediction is based on the term of the cycle number corresponding to the cyclic stress value. Yona et al. also imply that materials are sought to establish an “infinite” life in the S-N curve under corresponding loading condition for dental implants to clear the problem [6]. In Figure 5, the fatigue life prediction results are shown. The S-N curve data found in the literature are based on R = −1, so the solver is run by a loading ratio of R = −1.

Cycle results are in the interval of 4 × 105 and 1 × 109. For the evaluation of the results, the value of 5 × 106 cycles, which is recommended by the guideline of Dental Device Branch of FDA (U.S. Department of Health and Human Services-Food and Drug Administration) [49], can be taken as a reference for meeting the life cycle demand from implants. In this section, the implant modelled with cp Ti with the nanotubular TiO2 surface resulted in a life value of 1 × 107 cycles and CoCr alloy-assigned model resulted in a life of 1 × 109 cycles, so they can also be marked as meeting and even exceeding the reference life requirement. Remaining materials seem not enough for the life criteria for the related implant model.

On the other hand, the implants made of titanium alloys are known as they lead to failure of the treatment due to high cyclic loading that leads to resorption of the peri-implant bone, which increases the bending moments on implants and the metal fatigue, and so the fracture [36]. In contrast to that problem, the bone resorption is predicted to decrease by usage of the nanotubular TiO2 surface structure of the cp Ti. Contrary to the issue again, the bone growth and adsorption is probably increased; thus, the problem can be overcome. Briefly, the usage of cp Ti with the nanotubular TiO2 surface seems reasonable according to all of the FEA results.

3.4. Thread Number Effect on the Fixation of the Implant

Consequently, cp Ti with the nanotubular TiO2 surface was evaluated as the best performing material overall according to the conditions in this study. Then, an extra-analysis about fixation was conducted to understand the effect of the thread numbers on the implant structure. Since the fixation is another issue in the implant system, roughness surface may increase the friction and avoid the loosening of the implant screw. In this section, a method of generating an additional thread was applied to improve the fixation in the implant model. The two- and three-threaded implant designs having the same total pitch were developed, and then, they were analysed to observe the stress amount and the distribution state. The analysis results of the two- and three-threaded designs with the mechanical properties of cp Ti with the nanotubular TiO2 surface are shown in Figure 6.

These von Mises stress values are close to each other, and they are also reasonable; however, the result of the one-threaded design is 239.05 MPa. It is higher than the results of the multithreaded designs. The stress may be distributed on much more surface than the basic one; therefore, the stress values are resulted lower. On the other hand, the discontinuity due to the more grooves of the multithreads may also lead to notch effect and thus cause unpredictable failure or damage. Hence, this multithreaded design approach should be verified by tests in vitro. Design of sharp thread should be avoided, and the thread angles at the junction region of the implant body should also be reduced as possible. In this context, depth of thread should be developed for avoiding the interfering of the thread with the bone.

By means of this method, the friction can be increased and the potential failure of implant due to the implant screw loosening may be prevented. Moreover, different from the usage of porous or nanotubular structure, the roughness may also be increased by utilizing multiple threaded design or changing the thread section for getting the mechanical locking between the implant surface and the surrounding bone tissue.

4. Conclusion

The effects of biomaterials on the stress state and fatigue life of dental implants and bone interaction were discussed. Within this scope, various biomaterials, which are appropriate for implant application, were investigated under cyclic biting loads and the following results were obtained:(i)The mechanical properties of implant materials cause stress distribution on the implant. Especially, the elastic modulus affects the stress magnitude and distribution. Moreover, Poisson’s ratio is another factor that influences the terms, but it is in a mutual relation with the elastic modulus. They change the deformation and stress together.(ii)When the elastic modulus of an implant converges to the bone numerically, any applied load leads to stress and much more deformation on the implant. This fact provides smooth transmission of the stress from the implant to the bone and also decreases the transmitted energy magnitude. Therefore, the implant and the bone interaction becomes smooth. It may be evaluated that this case leads to more deformation in both implant and bone tissue. However, this fact may also be evaluated as less bone damaging, since the implant is not rigid much more and so it leads less pulverization of bone.(iii)However, the decrease of the elasticity leads to increase of the implant to bone penetration due to the increase in the deformation of the implant. Thus, it should also be recommended to research the medical interaction in this phenomenon.(iv)To provide more fixation, the friction must be increased. For this reason, the multiple threaded designs may also be a solution to obtain required mechanical locking between the implant surface and the surrounding bone.(v)The nanoporous/nanotubular TiO2 dental implant having multiple threaded designs may supply the needs of implantation technology. It has similar elastic modulus with bone tissue, and this provides smooth stress transformation from implant to bone. The nanoporous surface and the multiple threaded design increase the surface roughness and the friction. Hence, the fixation of the implant can be increased, and the bone growth on the surface of the implant can be supported. This implant can also display long-term stability under cyclic biting loads due to its novel design and its high strength.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.