Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 1598-6721(Print)
ISSN : 2288-0771(Online)
The Korean Society of Manufacturing Process Engineers Vol.17 No.6 pp.24-30
DOI : https://doi.org/10.14775/ksmpe.2018.17.6.024

Optimization of the Groove Depth of a Sealing-type Abutment for Implant Using a Genetic Algorithm

Hyeon-Yeol Lee*, Dae-Sun Hong*#
*School of Mechanical Engineering, Changwon National University.
Corresponding Author : dshong@changwon.ac.kr Tel: +82-55-213-3627, Fax: +82-55-263-5221
22/07/2018 31/07/2018 13/08/2018

Abstract


Dental implants are currently widely used as artificial teeth due to their good chewing performance and long life cycle. A dental implant consists of an abutment as the upper part and a fixture as the lower part. When chewing forces are repeatedly applied to a dental implant, gap at the interface surface between the abutment and the fixture is often occurred, and results in some deteriorations such as loosening of fastening screw, dental retraction and fixture fracture. To cope with such problems, a sealing-type abutment having a number of grooves along the conical-surface circumference was previously developed, and shows better sealing performance than the conventional one. This study carries out optimization of the groove shape by genetic algorithm(GA) as well as structural analysis in consideration of external chewing force and pretension between the abutment and the fixture. The overall optimization system consists of two subsystems; the one is the genetic algorithm with MATLAB, and the other is the structural analysis with ANSYS. Two subsystems transmit and receive the relevant data with each other throughout the optimization processes. The optimization result is then compared with that of the conventional one with respect to the contact pressure and the maximum stress. The result shows that the optimized model gives better sealing performance than the conventional sealing abutment.



유전자알고리즘을 이용한 임플란트용 실링어버트먼트의 홈 깊이 최적화에 관한 연구

이 현열*, 홍 대선*#
*창원대학교 기계공학부

초록


    © The Korean Society of Manufacturing Process Engineers. All rights reserved.

    This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

    1. Introduction

    Since the concept of osseointegration, which involves the direct contact between living bone and metal of an implant, was first reported by Branemark et al.,[1] dental implant technology has advanced rapidly on the back of experimental testing and clinical applications, and dental implants are now widely used as a reliable dental treatment.

    Fig. 1 shows an example of a dental implant consisting of a fixture, an abutment, and a prosthetic crown. The fixture acts a dental root implanted in the bone. At one end of the abutment is a column screwed into the fixture, which plays a role as an implant support[2]. The other end of the abutment is attached to the crown.

    After implantation in the oral cavity, a dental implant replaces the role of a tooth. As such, it is subject to a continuous dynamic load and various problems can occur because of this. In particular, foreign matter can come between the prosthetic materials due to loosening between the abutment and the fixture, which can significantly degrade the functions of the implant[3,4].

    A significant number of studies have been conducted to solve this problem, proposing modifications to the implant components, abutment, screw thread size and shape[5,6], chewing force, and pretension[7-9]. In particular, a sealing-type abutment has been developed to improve the contact performance between it and the fixture, characterized by the presence of a number of grooves around the circumference of the abutment (Fig. 2a)[10]. More recently, a complete sealing-type abutment[11] in which a gold ring is inserted into the top groove has been also developed (Fig. 2b).

    This study aimed to optimize the abutment groove depth (d) to further improve the sealing performance between the sealing abutment and the fixture, as shown in Fig. 2(a).

    The performance of the new design was compared with that of an existing product using structural analysis. A series of optimization simulations was conducted using the proposed design. The simulation results demonstrated that the optimized model reduced stress when an external force was applied compared to the existing product and thus improved sealing performance.

    2. Structural analysis of the abutment

    2.1 Structural analysis conditions

    Structural analysis was conducted on an existing sealing abutment prior to optimizing the groove dimensions. As shown in Fig. 2(a), the target product contained five grooves with a groove depth (d) of 0.05 mm, a groove height (h) of 0.1 mm, and a gap between the grooves (b) of 0.3 mm.

    The structural analysis followed the method presented in ISO14801,[12] which is applied to the fatigue testing of dental implant systems. ANSYS Classic[13] was used as the structural analysis tool. The abutment and the fixture consisted of titanium alloy Ti-6A1-4V. Table 1 presents the physical properties of this material[14].

    Two loading conditions were applied to the abutment: a chewing force, which acted as an external force, and pretension against the initial clamping force applied to the screw between the abutment and the fixture. Here, the chewing force was a 30° inclination force following the above ISO standard (250 N). The initial tightening torque between the abutment and the fixture was 0.35 Nm, and a pretension of 177.4 N was calculated using the formula[15] for screw clamping force and torque.

    If an external force is applied to the sealing abutment, inter-sliding occurs between the abutment and the fixture. In consideration of this contact, non-linear structural analysis should be conducted. The boundary conditions for the contact between the abutment and the fixture was set to "frictional", and the friction coefficient was set to 0.2. In ANSYS, the abutment was classified as the “contact surface” and the fixture as the “target surface,” while the outside of the fixture was classified as “fixed support.” Finally, because the implant was symmetrical, it was modeled as such, and the “symmetry” condition was applied to the cut section[13].

    2.2 Structural analysis results

    Finite element analysis was conducted on the existing sealing abutment, which was modeled prior to optimization based on the previously set boundary and loading conditions. The results are presented in Fig. 3. The stress distribution (Fig. 3a) indicates that maximum stress occurred in the uppermost area of the contact between the abutment and the fixture. The maximum stress was 241 MPa, and the safety factor was 4.0, which indicated that the strain occurred within the elastic section. The contact pressure distribution (Fig. 3b) revealed that the maximum contact pressure (355 MPa) occurred in the lowermost area of the contact between the abutment and the fixture.

    3. Optimization of the abutment grooves

    3.1 Groove shape and optimization variables

    A sealing abutment was fabricated by machining five grooves on the external side of a general abutment at regular distances. Fig. 4(a) shows the cross-section of the sealing abutment including the grooves, and Fig. 4(b) presents the tip of the cutting tool used to create the grooves. This study employed the tool in Fig. 4(b). The gap between the grooves (b) shown in Fig. 4(a) was fixed at 0.3 mm as in the existing product. The groove shape was modified by changing the groove depth d only as marked in the two figures. Groove depth d was set as the optimization variable while the number of grooves was set at five.

    3.2 Configuration of the optimization system

    A genetic algorithm[16] was used to optimize the shape of the grooves in the sealing abutment. In the genetic algorithm, the optimization variable (in this case, groove depth) is represented by a single entity consisting of multiple strings or bits, and a group is formed by a number of entities. This genetic algorithm allows non-linear complex problems to be solved through the evolution of the group[16]. The optimization algorithm was run using MATLAB to search for the optimal groove depth. MATLAB and ANSYS exchanged data in real time to perform the optimization task via the sharing of text files. Data during the search processes was delivered to ANSYS, which performed structural analysis based on the variable. The maximum stress and maximum contact pressure were calculated, and the data transferred back to MATLAB iteratively.

    3.3 Genetic algorithm

    To acquire an optimization solution using a genetic algorithm, a fitness function F(x) first needs to be defined. Generally, a fitness function consists of objective and penalty functions. In this study, the objective function O(x) was set to the following equation to ensure it was proportional to the contact pressure Pcontact between the abutment and the fixture.

    O ( x ) = 1 + P c o n t a c t / σ Y
    (1)

    Here, σ Y is 970 MPa, which is the yield point, and contact pressure Pcontact is the maximum contact pressure analyzed in ANSYS for the corresponding entity. As presented in the above equation, the objective function value is always larger than 1.

    The penalty function P(x) was set as follow considering the yield point σ Y and safety factor Sf.

    P ( x ) = { 1 ,  if  σ max > σ Y / S f 0 , o t h e r w i s e
    (2)

    Here, σ max is the maximum stress value obtained after analyzing the corresponding entity in ANSYS, and the safety factor Sf was set to 5. As presented above, a penalty was applied if the maximum stress of the abutment was larger than the allowable stress ( σ Y / S f ) .

    Next, the fitness function Fx was set as follows using the above two equations.

    F ( x ) = O ( x ) P ( x )
    (3)

    As shown in the equation, O(x) is always larger than 1, and the maximum value of P(x) is 1. Thus, F(x) was set to be a positive value that was always larger than 0.

    Table 2 presents the parameter values for the genetic algorithm used in the optimization. An individual consisted of genes that have 60 binary values. If a population size is too small, it will quickly converge to a local solution. If the population size is too large, the computation time will increase disproportionately compared to the increase in performance. Thus, this study set the population size to 8. In addition, the maximum number of generations, which establishes when the genetic algorithm will terminate, was set to 30 after identifying the convergence trend through simulations. The crossover probability (Pc) and mutation probability (Pm) were also set within generally known ranges (Pc: 80-95% and Pm: 0.01 -1%)[17].

    Fig. 5 shows the overall optimization process for the genetic algorithm and structural analysis.

    3.4 Optimization results and analysis

    As described above, the optimization was conducted by setting the groove depth d as the optimization variable (Fig. 4a) using the proposed genetic algorithm.

    Fig. 6 shows the change in the fitness function by generation. Groove height h was automatically determined according to depth d by the shape of the tool that produces the groove (Fig. 4b). As shown in Fig. 6, the fitness function tended to increase with the number of generations. Because no change was shown after 11 generations, this was regarded as the optimal solution.

    The groove depth d that corresponded to the optimum solution was found to be 0.031 mm, a reduction of 0.019 mm from the existing model (0.05 mm). Accordingly, the groove height was slightly educed as well, and the protrusion area in the abutment increased, thereby making the contact surface slightly larger.

    The results of the structural analysis for the optimal solution are shown in Fig. 7. The maximum stress of the optimization model was 170 MPa, a reduction of 71 MPa (29.5%) compared to the initial model (241 MPa). The maximum contact pressure was 870 MPa, an increase of 515 MPa (145.0%) compared to the initial model (355 MPa). The above results are summarized in Table 3.

    As presented in the summary, the maximum stress was reduced by approximately 30% compared to that of the existing model under the same xternal force while the contact pressure increased 1.5 times. This result indicates that the sealing performance was improved by the optimization model.

    4. Conclusions

    This study performed structural analysis and dimension optimization of five grooves on the outer surface of a sealing abutment using a genetic algorithm to improve the sealing performance abutment and the fixture in a dental implant.

    The number of grooves and the gap between the grooves were set to be the same as an existing product, while the groove depth was set as the optimization variable. Based on the results of the optimization process, the performance of the proposed design and that of the existing product were compared using structural analysis. The optimization results showed that the optimal groove depth was slightly smaller than that of the existing model, which reduced stress by approximately 30% when an external force was applied, while the maximum contact pressure, which is directly related to sealing performance, increased by 145%. These results indicate that the optimized sealing abutment model improved sealing performance over the existing model.

    Acknowledgments

    This research is financially supported by Changwon National University in 2017~2018.

    Figure

    KSMPE-17-24_F1.gif
    Constitution of a dental implant
    KSMPE-17-24_F2.gif
    Two types of sealing abutments recently developed
    KSMPE-17-24_F3.gif
    Results of structural analysis on the existing sealing abutment model
    KSMPE-17-24_F4.gif
    Configuration of the cutting tool and the grooves of the abutment
    KSMPE-17-24_F5.gif
    Flow chart of the proposed genetic algorithm
    KSMPE-17-24_F6.gif
    Change of the fitness function according to generation
    KSMPE-17-24_F7.gif
    Structural analysis on the optimized model

    Table

    Material property of the titanium alloy
    Parameters of the GA
    Comparison of the optimized model with the initial model

    Reference

    1. Brånemark, P. I., Breine, U., Adell,, R. Hansson, B. O., Lindstrom, J. and Olsson, A., “Intraosseous Anchorage of Dental Prostheses. I. Experimental Studies”, Scand J. Plast Reconstr Surg., Vol. 3, pp.81-100, 1969.
    2. Carlson B. and Carlsson G., "Prosthodontic Complications in Osseointegrated Dental Implant Treatment,” Int. J. Oral and Maxillofacial Implants, Vol. 9, No. 1, pp.90, 1994.
    3. Kim, K, N, Choi, K, H. and Oh, H, S., “The evaluation of Fitness on Implant Bolt Joints as Tightening Torques,” The Korean Society of Mechanical Engineers, Vol. 34, No. 11, pp.1786-1791, 2010.
    4. Kwon, Y. J., “Structural Design of a Dental implant (I) : Comparative Structural Analysis,” Journal of Computation Structural Engineering Institute of Korea, Vol. 25, No. 5 pp. 421-432, 2012.
    5. Siamos, G., Winkler S. and Boberick, K.G., “The Relationship between Implant Preload and Screw Loosening on Implant-Supported Prostheses,” Journal of Oral Implantology, Vol. 28, No. 2, pp.67-73, 2002.
    6. Haack, J. E., Sakaguchi, R. L., Sun, T. and Coffey J.P., "Elongation and Preload Stress in Dental Implant Abutment Screws," The International journal of oral & maxillofacial implants, Vol. 10, No. 5, pp. 529-536, 1995.
    7. Martin, W. C., Woody, R. D., Miller, B. H. and Miller, A.W., “Implant Abutment Screw Rotations and Preloads for Four Different screw Materials and Surfaces”, The Journal of Prosthetic Dentistry, Vol. 86, No. 1, pp.24-32, 2001.
    8. Lee, M. Y., Heo, S. J., Park, E. J. and Park, J. M., “Comparative Study on Stress Distribution around Internal Tapered Connection Implants according to Fit of Cement-and-Screw-Retained Prosthesis,” The Journal of Advanced Prosthodontics, Vol. 5, No. 3, pp.312-318, 2013.
    9. Kim, N., Choi, K. H. and Oh, H, S., “The Evaluation of Fitness on Implant Bolt Joints as Tightening Torques,” Proc. of The Korean Society of Mechanical Engineers of Autumn Conference, pp. 1786-1791, 2010.
    10. Kim, J. M. and Hong, D. S., "Development of a Sealing-Type Abutment for Implant and the Performance Evaluation via Structural Analysis," J. of the Korean Society for Precision Engineering, Vol. 33, No. 9. pp.769-775, 2016.
    11. Hong, D. S., Jeon, Y. J., Lee, H. Y. and Lee, S. Y., "Development and Performance Evaluation of a Sealing-Type Abutment Inserted with a Gold Ring," J. of the Korean Society of Manuf. Process Engineers, Vol. 15, No. 2. pp.97-103, 2016.
    12. ISO (International Organization for Standardization), ISO 14801 : Dentistry- Implants-Dynamics Fatigue Test for Endosseous Dental Implants, pp.1-16, 2016.
    13. ANSYS Inc., ANSYS Mechanical APDL Structural Analysis Guide, TEA SUNG S&E Inc, 2013.
    14. Boyer, R., Welsch, G. and Collings, E. W. eds, Materials Properties Handbook: Titanium Alloys, ASM International, Materials Park, 1994.
    15. Hong, J. P., Mechanical Design, BooksHill,pp.100-106, 2001.
    16. Goldberg, D. E., Genetic algorithms in search, optimization, and machine learning, Addison-Wesley, 1989.
    17. Jin, G. K., Genetic Algorithms and its Applications, Kyowoo Publication, pp. 113-194, 2002.