Acta Polytechnica CTU Proceedings doi:10.14311/APP.2017.13.0097 Acta Polytechnica CTU Proceedings 13:97–101, 2017 © Czech Technical University in Prague, 2017 available online at http://ojs.cvut.cz/ojs/index.php/app GEOMETRICAL OPTIMIZATION OF DENTAL IMPLANTS WITH REGARD TO OSSEOINTEGRATION Luboš Řehounek∗, František Denk, Aleš Jíra Czech Technical University in Prague, Faculty of Civil Engineering, Thákurova 7, 166 29 Prague 6, Czech Republic ∗ corresponding author: lubos.rehounek@fsv.cvut.cz Abstract. A newly developed reference dental implant specimen type was subjected to numerical simulations of osseointegration. The goal of these tests was to optimize the geometry of the implant so as to reduce local stress concentrations and provide a better flow of stress through the whole implant body. Conditions for osseointegration were considered when evaluating the anchoring system of the implant in regard to its placement in the human cancellous and cortical bone. Numerical simulations showed that stress concentrations occur mostly in the upper cylindrical part of the implant. By increasing the width of this cylindrical part, we were able to reduce the maximum values of von-Mises stress by 20 %. Keywords: Dental implant, optimization, trabecular, numerical, FEM. 1. Introduction The main shortcoming of commonly manufactured dental implants is their predisposition towards early failure [1–3]. This unwanted factor can be attributed to uneven distribution of stresses in the implant, peri- implant area and bone [4, 5]. Our main goal in de- veloping new implants was to reduce stress concen- trations occuring mainly in the upper cylindrical part of the implant, just below the abutment, while pre- serving the osseointegration benefits of the implant. The tested implant has been patented in the Czech Republic [6]. The numerical analysis describes the mechanical be- havior of the implant in regard to mechanical loads and the corresponding biomechanical stablity in the bone tissue. The tests focus on two aspects – mechanical behavior of the implant under a prescribed mechan- ical load and osseointegration of the implant inside the human mandible. The first part of the analysis fo- cuses on the osseointegration of the implant in human mandible. The aim of the second part of the tests was to determine the ultimate bearing capacity of the den- tal implant stem according to the ČSN EN ISO 14081 – “Stomatology – Implants – Dynamic Fatigue Testing of Intraosseous Dental Implants” standard. This pa- per will only focus on the first part of the tests (os- seointegration) and leave the second part (mechanical tests) for further studies, as some results are already available, but additional simulations are still needed. The analysis provides an overview of stress and de- formation relations of the tested implant stem in a pre- scribed environment in all of its parts. This is sim- ulated by embedding the implant stem in a block of human bone (Fig. 2) and subsequently evaluating the stress and deformation relations at surface areas of the stem and at the interfaces of the stem and bone tissue. Numerical simulations were done in the ANSYS APDL system by incorporating previously created 3D geometries (Fig. 1) into the environment of the soft- ware. The geometry of the implant stem was created in detail at the surface areas of the implant. The inter- nal part of the stem has been created with less detail (i.e.: without the screw-thread). 2. Numerical Testing of Osseointegration For the purpose on numerical analysis in ANSYS APDL, we created an environment where the den- tal implant is fixed into the bone tissue of human mandible. The implant consists of the intraosseous stem and the implant superstructure, which repre- sents the abutment and the crown. The connection of individual volumetric bodies is attained using rigid constraints generated by the volumetric “vglue” oper- ations. The mesh of the model consists of quadratic tetrahedrons type SOLID 187 generated in the “smart size” mode. The material of all components of the den- tal implant is the titanium alloy Ti-6Al-4V. The stem of the implant is fixed into a volumetric model of human mandible, situated according to real implantation procedures. The geometrical orientation of the ribs was arranged in the least favorable position (in regard to the results of mechanical tests). The ex- act geometry of the bone was generated by Boolean op- erations (subtracting the full body volume of the stem from the body of the bone). The connection between contact surfaces of the stem and the bone was at- tained by volumetric “vglue” operations. The numer- ical model of the mandible was then subsequently remeshed according to the new boundary conditions (Fig. 2). The front faces of the mandibles have symmetric 97 http://dx.doi.org/10.14311/APP.2017.13.0097 http://ojs.cvut.cz/ojs/index.php/app L. Řehounek, F. Denk, A. Jíra Acta Polytechnica CTU Proceedings Figure 1. 3D model of the designed implant 0001 (ANSYS APDL): whole 3D model (left), longitudinal section (middle) and FEM mesh (right). Figure 2. Numerical model of the mandible block which supports the intraosseous part of the implant stem – volumetric model (left), FEM mesh (middle) and total deformation isolines (right). boundary conditions and the lower strip of nodes is set to have zero displacement value. Physiological load is represented by a single slant force applied to the central point of the upper base of the super- structure. The force magnitude is F = 130 N. The di- rection of the force is defined by contribution of indi- vidual force components according to the total load imposed on a single tooth. The magnitude of normal force is 125 N and tangent force is 35 N (Fx = 25 N, Fy = 25 N, Fz = 125 N). 3. Optimization of the Cylindrical Part The following optimization will assume that the small strip of the cylindrical stem forming the upper struc- ture of the stem is always the area where values of stress concentrate. These extremities are caused Figure 3. Numerical analysis of the reference implant 0001: displacement – umax = 25 µm (left), equivalent von Mises stress – σeqv,max = 91.7 MPa (middle and right). 98 vol. 13/2017 Geometrical Optimization of Dental Implants Figure 4. Optimization of the geometrical shape in the upper cylindrical part: reference model 0001A (cyl. width 0.50 mm) (left), model 0001B (cyl. width 1.25 mm) (middle) and model 0001C (cyl. width 3.00 mm) (right). Figure 5. Numerical analysis of the reference im- plant 0001B: displacement – umax = 25 µm (left), equivalent von Mises stress – σeqv,max,1 = 97.1 MPa, σeqv,max,1 = 73.2 MPa (middle and right). by complex geometrical shapes and also a relatively small width of the cylindrical strip. The goal of the op- timization will be to effectively modify the geomet- rical shape of the implant in the area of fixation of the stem into cortical bone of the human mandible. Also, the overall stability of the intraosseous part must not be compromised. Two new variants have been tested – variant 0001B and 0001C. The modification represents increasing of the width of the cylindrical part. The width of the cylindrical part of the reference implant 0001A is 0.50 mm. The model 0001B had the width increased up to 1.25 mm and model 0001C up to 3.00 mm. Identical simulations were run on all of the men- tioned variants and the outcomes of the numerical analyses are shown in Fig. 3, Fig. 5 and Fig. 6. On the basis of equivalent von-Mises stress dis- Figure 6. Numerical analysis of the reference implant 0001C: displacement – umax = 25 µm (left), equiva- lent von Mises stress σeqv,max = 73.6 MPa (middle and right). tribution at the surface areas of the imlpant stem, we can observe that the widening of the cylindrical strip of the upper part without any other geometrical modification positively affects the stress distribution by lowering extremities. The variant 0001B (cylindrical part 1.25 mm) has a mean value of σeqv,avg = 61.6 MPa and the stress flow is distributed over a larger area. There is an ex- treme of σeqv,max = 73.2 MPa. There is also another extreme of σeqv,max = 97.1 MPa spread over a very small area situated in the edge of the cylindrical part. The system of vertical ribs is situated in both cortical and cancellous bone. The variant 0001C (cylindrical part 3.00 mm) has a maximum value of σeqv,max = 73.6 MPa, but the sys- tem of vertical ribs is situated now only in cancellous bone. 99 L. Řehounek, F. Denk, A. Jíra Acta Polytechnica CTU Proceedings Implant type osseointegrated in bone tissue of mandible Displacement of the implant and bone tissue Equivalent von-Mises stress Optimization influence ut [µm] σeqv,max [MPa] [%] Implant 0001A 25 91.70 100.00 Implant 0001B 25 73.20 79.80 Implant 0001C 25 73.60 80.00 Table 1. Overview of the results provided by geometrical modification. Reference implant 0001A and modified implants 0001B and 0001C. Figure 7. Two STL files showing the final geometrical solution of the implant (trabecular outer shell second from left) (left), the final implant model (middle) and homogenized layer of the outer trabecular structure for easier simulations (right). 4. Conclusions 4.1. Results The results of the numerical simulations show that by increasing the upper cylindrical portion of the stem, we can improve the distribution of stresses at the sur- face areas of the implant and also reduce extreme values of stress. The variant 0001B showed that by increasing the width of the cylindrical part by 0.75 mm (from 0.50 mm to 1.25 mm), we can reduce stresses by ap- proximately 20 % and achieve a more even distribution. There is, however, an extreme value of 97.1 MPa sit- uated at the edge of the cylindrical part over a very small area. This extreme is caused by discontinuities in geometry (sudden change of curvature), which will be accounted for by rounding the edge of the cylinder by a radius of at least 0.1 mm during the production process. Biomechanical stability is achieved by situat- ing the vertical ribs of the implant inside the cortical bone of the mandible. The variant 0001C (width from 0.75 mm to 3.00 mm, 2.50 mm increase) also showed a 20 % reduction in val- ues of stress. The cylindrical part of the implant does not permeate through the whole layer of cortical bone, so the distribution of stress is more even and there are no larger extremes in values of stress. This geomet- rical solution, however, might not provide sufficient biomechanical stability as the vertical ribs, which are supposed to anchor the implant in the bone, are only situated in cancellous bone. Considering the biomechanical implications of the presented results, we evaluated the variant 0001B with a 1.25 mm wide cylindrical strip as the most efficient solution, which reduces extreme values of stress as well as provides sufficient stability inside thevbone tissue. 4.2. Future Prospects The chosen variant 0001B has good mechanical stabil- ity and stress distribution and was chosen as the vari- ant for further development. With the shape solution of the implant stem complete, the next step in the re- search will be to create a numerical model of a com- posite implant (Fig. 8). This implant will comprise of both homogeneous and trabecular portions. Tra- becular structure has been discussed at length in one of other author’s works [7]. The implant will maintain its homogeneous stem. The outer shell with the stabilizing ribs will, how- ever, be comprised exclusively of the trabecular struc- ture. Further simulations will be made and multiple variants of different layer thickness analysed in order to determine which shape solution is the best for os- seointegration, mechanical stability and machanical bearing capacity. The approach towards modelling the trabecular structure in future research will be focused on sustain- able, repeatable analyses. We have already found this 100 vol. 13/2017 Geometrical Optimization of Dental Implants Figure 8. Whole implant body of the implant for further research. Implant has a homogeneous stem and an outer layer comprising of the trabecular struc- ture. The implant has been patented in Czech Repub- lic [6]. to only be possible by making a virtual homogenized material (Fig. 7) with specific curve-fitted properties determined by experiment [7]. This approach allows us to apply previously determined mechanical proper- ties of the trabecular structure onto a homogeneous material that is easy to analyze in numerical environ- ment. By running a few simple analyses, we found this to be the optimal approach as it is not very demand- ing in computation times and it maintains sufficient precision. By manipulating the shape of the outer trabecular structure, we will then be able to deter- mine which variant is optimal in regard to reduction of stress shielding and mechanical bearing capacity. Acknowledgements The financial support by the Faculty of Civil Engineering, Czech Technical University in Prague (SGS project No. SGS17/168/OHK1/3T/11) is gratefully acknowledged. References [1] L. Le Guéhennec, A. Soueidan, P. Layrolle, Y. Amouriq. Surface treatments of titanium dental implants for rapid osseointegration. Dental Materials 23(7):844–854, 2007. doi:10.1016/j.dental.2006.06.025. [2] R. Osman, M. Swain. A Critical Review of Dental Implant Materials with an Emphasis on Titanium versus Zirconia. Materials 8(3):932–958, 2015. doi:10.3390/ma8030932. [3] T. Jemt, M. Olsson, F. Renouard, et al. Early Implant Failures Related to Individual Surgeons: An Analysis Covering 11,074 Operations Performed during 28 Years. Clinical Implant Dentistry and Related Research 18(5):861–872, 2016. doi:10.1111/cid.12379. [4] M. I. Z. Ridzwan, S. Shuib, A. Y. Hassan, et al. Problem of stress shielding and improvement to the hip implant designs: A review, 2007. doi:10.3923/jms.2007.460.467. [5] M. Long, H. J. Rack. Titanium alloys in total joint replacement–a materials science perspective. Biomaterials 19(18):1621–1639, 1998. doi:10.1016/S0142-9612(97)00146-4. [6] F. Denk, F. Denk, A. Jíra. A dental implant shaft, Czech Republic, 306457, 2017. [7] L. Řehounek, A. Jíra. Experimental and numerical analyses of a 3D-printed titanium trabecular dental implant. Acta Polytechnica 57(3):218–228, 2017. doi:10.14311/AP.2017.57.0218. 101 http://dx.doi.org/10.1016/j.dental.2006.06.025 http://dx.doi.org/10.3390/ma8030932 http://dx.doi.org/10.1111/cid.12379 http://dx.doi.org/10.3923/jms.2007.460.467 http://dx.doi.org/10.1016/S0142-9612(97)00146-4 http://dx.doi.org/10.14311/AP.2017.57.0218 Acta Polytechnica CTU Proceedings 13:98–102, 2017 1 Introduction 2 Numerical Testing of Osseointegration 3 Optimization of the Cylindrical Part 4 Conclusions 4.1 Results 4.2 Future Prospects Acknowledgements References