Numerical Simulation of Mandible Bone Remodeling under Tooth Loading: A Parametric Study

• Published in: Scientific reports Vol. 9; no. 1; pp. 14887 - 12
• Main Authors: Su, Kangning, Yuan, Li, Yang, Jie, Du, Jing
• Format: Journal Article
• Language: English
• Published: London Nature Publishing Group UK 17.10.2019; Nature Publishing Group
• Subjects: 631/57/2266; 639/166/985; Algorithms; Biomechanical Phenomena; Bite Force; Biting; Bone density; Bone Density - physiology; Bone remodeling; Bone Remodeling - physiology; Boundary conditions; Chewing; Computer Simulation; Cone-Beam Computed Tomography; Convergence; Dental Models; Finite Element Analysis; Humanities and Social Sciences; Humans; Loading; Mandible; Mandible - diagnostic imaging; Mandible - physiology; Mastication - physiology; Mathematical models; multidisciplinary; Parametric statistics; Science; Science (multidisciplinary); Stress, Mechanical; Teeth; Tooth - diagnostic imaging; Tooth - physiology
• ISSN: 2045-2322; 2045-2322
• DOI: 10.1038/s41598-019-51429-w

Bone adapts to the change of mechanical stimulus by bone remodeling activities. A number of numerical algorithms have been developed to model the adaptive bone remodeling under mechanical loads for orthopedic and dental applications. This paper examines the effects of several model parameters on the computed apparent bone density in mandible under normal chewing and biting forces. The density change rate was based on the strain energy density per unit mass. The algorithms used in this study containing an equilibrium zone (lazy zone) and saturated values of density change rate provides certain stability to result in convergence without discontinuous checkerboard patterns. The parametric study shows that when different boundary conditions were applied, the bone density distributions at convergence were very different, except in the vicinity of the applied loads. Compared with the effects of boundary conditions, the models are less sensitive to the choice of initial density values. Several models starting from different initial density values resulted in similar but not exactly the same bone density distribution at convergence. The results also show that higher reference value of mechanical stimulus resulted in lower average bone density at convergence. Moreover, the width of equilibrium zone did not substantially affect the average density at convergence. However, with increasing width, the areas with the highest and the lowest bone density areas were all reduced. The limitations of the models and challenges for future work were discussed for the better agreement between the computed results and the in vivo data.