Predicting the holistic force-displacement relation of the periodontal ligament: in-vitro experiments and finite element analysis
© Chang et al.; licensee BioMed Central Ltd. 2014
Received: 3 April 2014
Accepted: 16 July 2014
Published: 30 July 2014
The biomechanical property of the periodontal ligament (PDL) is important in orthodontics and prosthodontics. The objective of this study was to evaluate the feasibility of measuring the biomechanical behavior of the periodontal ligament using micro-computed tomography (micro-CT).
A custom-made apparatus measured the force and displacement of a porcine PDL specimen within the micro-CT environment. Synchronized computed tomography (CT) images were used to obtain the deformation and displacement of the entire specimen and to reconstruct the three-dimensional mesh model. To match the experimental results, finite element analysis was then applied to simulate the biomechanical response of the PDL. The mechanical model of the PDL was assumed as the hyperelastic material in this study.
The volume variations of the tooth and the alveolar bone were less than 1%, which implies that tooth displacement was caused mostly by displacement of the PDL. Only translational displacement was observed with each load step because the transformation matrix acquired from the CT image registration was identical. The force-displacement curve revealed the nonlinear behavior of the PDL. There was a high correlation between the experimental displacement results and the simulation displacement results. The numerical results (based on the assumption that the PDL is the hyperelastic material) showed good agreement with the experimental results.
Nondestructive measurements by micro-CT obtained the biomechanical behavior of the PDL. Using the hyperelastic characteristic as the constitutive model can properly predict the force-displacement relation of the PDL after loading. This study provided a feasible approach for measuring the biomechanical behavior of the PDL for further dental application.
The periodontal ligament (PDL) has a determinative role in dental biomechanics. The PDL connects the tooth and alveolar bone, absorbs occlusal impact, and forms and resorbs the alveolar bone for tooth movement. Hence, the biomechanical property of the PDL is important in orthodontics and prosthodontics. Orthodontic treatment is a reiterative and time-consuming process during which a tooth gradually moves to its prospective position. The physiological mechanism of tooth movement is primarily because of the response of the PDL while an orthodontic force is applied .
The PDL is clearly important, but its biomechanical behavior remains unclear. Toms et al.  performed a shear test by using a transverse thin specimen of the premolar tooth. They revealed the anisotropic behavior of the PDL. Dorow et al.  designed an apparatus to test the response of the porcine anterior tooth under different strain rates. Their results showed that the nonlinear hysteresis behavior varied with the strain rate. Natali et al.  designed a laser-optical system to measure the minipig premolar tooth. They obtained the nonlinear force-displacement relation of the PDL. In 2009, Tohill et al.  acquired the relaxation characteristic of the premolar porcine tooth. However, these studies only demonstrated the simplex mechanical feature of the PDL.
To date, the PDL cannot be harvested by a nondestructive approach for evaluating the stress–strain relation. Analytical models have consequently become a popular approach for estimating the biomechanical behavior of the PDL. There are three types of governing equation for these models: (1) linear elastic model, (2) viscoelastic model, and (3) hyperelastic model. The linear elastic behavior of the PDL has been widely introduced in previous literature reports, but the Young’s modulus values range from 0.07 MPa to 1750 MPa . However, some studies indicate that the elastic property of the PDL is nonlinear [7–9]. Previous analytical models have assumed viscoelastic properties for the PDL such as creep, stress relaxation, and hysteresis [10–13]. However, most PDL researchers prefer the linear elastic behavior because of its simplicity. Because the force-displacement relation of the PDL is insufficient, an energy-based equation was developed and named the “hyperelastic model”. To predict tooth mobility or movement, many researchers have tried to establish a hyperelastic model of the PDL [4, 13–15]. Choosing hyperelastic models and parameters is arduous and may directly influence the simulation results. The dilemma of predicting PDL behavior has limited the advancement of dental biomechanics.
Micro-computed tomography has recently been widely applied for dental research in the microscale such as dental morphology , dental biology , dental biomechanics [18, 19], and dental material science . This method can detect the precise morphology and can inspect the details of the microstructure or microvariations (e.g., movement or deformation) in biological tissues. Therefore, we aimed to measure the biomechanical behavior of the PDL using micro-CT.
Custom device for installation into the micro-CT scanner
Finite element analysis
in which C 10 , C 01 , and C 11 are the material constants and d is the material incompressibility parameter. and are the first and second deviatoric strain invariant. The initial estimation of these material constants (C 10 = 0.04 MPa; C 01 = 0.02 MPa; C 11 = 0.04 MPa; and d = 0.02), to which previous researchers have referred , was established to obtain a proper combination for fitting the experiment results.
The boundary and loading conditions were followed by the experiments. The surface nodes of the perforated hole in the alveolar bone were fixed in every direction and the loadings, according to the experimental results, were applied to the surface nodes of the perforated hole in the tooth. The all interfaces between tissues were continuous to ensure the loading could be properly transmitted.
The variations in the volume of the alveolar bone and the tooth at each load step (relative to zero force and displacement)
Alveolar Bone Volume (mm 3)
Tooth Volume (mm 3)
In this study, micro CT was demonstrated to be a useful tool for measuring the biomechanical behavior of the periodontal ligament. This approach can be used to determine the response of the PDL according to the entire displacement with external stimuli applied to the tooth. Many studies have investigated the biomechanical behavior of the PDL by using specimens from animals such as rats, rabbits, and pigs. [4, 7, 21, 22]. Animal specimens are used because intact specimens of the tooth and adjacent alveolar bone from human cadavers are difficult to harvest. This study preserved the intact porcine PDL and adjacent bone for the investigation of PDL biomechanics. Furthermore, using the custom-made apparatus in the micro-CT environment made it possible to simultaneously record the displacement and the corresponding force of the tooth and the PDL. This setup may create an in vitro environment that is biomechanically compatible with the loading experienced by the tooth and PDL during occlusion.
Although it is easier to apply a simple linear elastic model, this assumption was not suitable for the soft tissue, which was well-known as a nonlinear material. The general behavior of the soft tissue was divided into three regions: toe, linear and failure regions. When the loading of the soft tissue increased initially, the stiffness of the toe region was low; the linear region was identified after the toe region and had a higher stiffness than the toe region; then the soft tissue failed eventually. It is not easy to clarify the transition zone between each region. Furthermore, the viscoelastic model was also not suitable in this study because the experimental results did not reveal the time-dependent characteristic of the PDL. Therefore, the hyperelastic model was more favorable to fit the nonlinear force-displacement behavior of the PDL.
The results of previous studies represented the local characteristic of the PDL instead of the whole tooth structure; hence, the biomechanical behavior would greatly differ among these studies. Therefore, finding a suitable constitutive model of the PDL is an important issue in dental biomechanics, especially in the development of prosthodontics and orthodontics. However, the parameters reported in this study resulted in greater error when a small force was applied. In addition, each tooth might have its own constitutive equation according to the morphology; therefore, the parameters may still need further validation.
From the sectional CT images, we observed that the transversal thickness of the PDL differed, based on the location of the tooth. The transversal thickness of the PDL was greater around the alveolar crest and the apical region than in other regions. Previous studies have also found similar results [23, 24]. In addition, most specimens were sliced transversely into a thin plate for performing a shear test , which could only represent the local biomechanical behavior of the PDL. Therefore, this study demonstrated the global biomechanical behavior of the PDL and may have greater differences in comparison to other studies. Furthermore, because multiple roots of a tooth could disturb the biomechanical performance of the PDL during occlusal pressure, this study applied a tension force to highlight the independent effect of the PDL. The resistance force increased dramatically with the displacement after passing the toe region, based on the experimental results (Figure 6). This would explain the protective mechanism that prevents the tooth from excessive extrusion.
Testing in the micro-CT environment has several advantages compared to the previous experiments. First, it is a nondestructive modality for performing biomechanical testing. Further biological or physiological analysis can be performed because the PDL specimen remains intact. Second, micro or local observations can be accomplished because micro-CT has a higher resolution than traditional modalities such as CT, magnetic resonance imaging, or ultrasound images. Analysis of the intact PDL or observation of internal morphology is available by 3D model reconstruction because no permanent damage occurs within the physiological loading range.
Bone remodeling is an important process in orthodontic biomechanics. Studies report a direct correlation of stress–strain fields in the PDL with alveolar bone resorption [25, 26]. While PDL can be a direct effector, the response of bone to mechanical load could be a delayed effect that is necessary for tooth translation. Instead of a time-dependent (i.e., viscoelastic) characteristic of the PDL, our results showed good agreement between the numerical simulation and the experiment. This implied that the hyperelastic constitutive model is practical for investigating the biomechanical response of the PDL. The mechanobiological interaction of the PDL and alveolar bone nevertheless remains indeterminable.
Our results demonstrated that the PDL is suitable for simulating hyperelastic material; however, some limitations still need to be depicted. Because of the scattering effect of the micro-CT, the material selected for the custom-made apparatus should be taken in account to eliminate its influence in image registration. The experimental results and the hyperelastic model can well interpret the steady state of the PDL after loading because the entire scanning process requires more than 30 minutes. More samples under various loading modes are required for a complete evaluation of the biomechanical property of the PDL. The detailed tooth structures were not modelled in this study, such as enamel, dentine and pulp. The Young’s modulus and stiffness of the anatomical layer of enamel and dentine was much greater than those of the PDL. In addition, the volume of pulp took a very small part of the entire tooth. Furthermore, from the results of micro-CT reconstruction, the tooth and bone was approximately undeformed. Therefore, the lack of these layers in FE model would not affect the results much. This study excluded the behavior of the PDL at a high strain rate. For this study, only one specimen of a porcine mandibular premolar tooth was harvested. However, finding a suitable constant set of the hyperelastic model was a time-consuming task. Besides, the individual difference was existed in every sample harvested from the porcine. Although more specimens may have enhanced the reliability of the experimental and numerical results, this study successfully established a feasible approach to predict the force-displacement relation of the PDL in a holistic observation.
Micro-computed tomography, which performs nondestructive measurements, is capable of obtaining the biomechanical behavior of the PDL. Using the hyperelastic characteristic as the constitutive model may properly predict the force-displacement relation of the PDL after loading.
We would like to acknowledge the technical support of medical image reconstruction from School of dentistry, China Medical University, Taiwan. This study was partly supported by the National Science Council, Taiwan (grant number: NSC 102-2221-E-214-012).
- Proffit WR, Fields HW, Sarver DM: Contemporary Orthodontics. 5th edition. St. Louis, MO: Elsevier/Mosby; 2013.Google Scholar
- Toms SR, Lemons JE, Bartolucci AA, Eberhardt AW: Nonlinear stress–strain behavior of periodontal ligament under orthodontic loading. Am J Orthod Dentofacial Orthop 2002, 122: 174–179.View ArticleGoogle Scholar
- Dorow C, Krstin N, Sander F-G: Determination of the mechanical properties of the periodontal ligament in a uniaxial tensional experiment. J Orofac Orthop 2003, 64: 100–107.View ArticleGoogle Scholar
- Natali AN, Carniel EL, Pavan PG, Bourauel C, Ziegler A, Keilig L: Experimental-numerical analysis of minipig’s multi-rooted teeth. J Biomech 2007, 40: 1701–1708.View ArticleGoogle Scholar
- Tohill R, Hien M, McGuinness N, Chung L, Reuben RL: Measurement Of The Short-Term Viscoelastic Properties Of The Periodontal Ligament Using Stress Relaxation. In 4th European Conference of the International Federation for Medical and Biological Engineering; 23–27 November, 2008. Edited by: Sloten J, Verdonck P, Nyssen M, Haueisen J. Belgium: Springer; 2009:1467–1470.Google Scholar
- Rees JS, Jacobsen PH: Elastic modulus of the periodontal ligament. Biomaterials 1997, 18: 995–999.View ArticleGoogle Scholar
- Kawarizadeh A, Bourauel C, Jager A: Experimental and numerical determination of initial tooth mobility and material properties of the periodontal ligament in rat molar specimens. Eur J Orthod 2003, 25: 569–578.View ArticleGoogle Scholar
- Poppe M, Bourauel C, Jager A: Determination of the elasticity parameters of the human periodontal ligament and the location of the center of resistance of single-rooted teeth a study of autopsy specimens and their conversion into finite element models. J Orofac Orthop 2002, 63: 358–370.View ArticleGoogle Scholar
- Ziegler A, Keilig L, Kawarizadeh A, Jager A, Bourauel C: Numerical simulation of the biomechanical behaviour of multi-rooted teeth. Eur J Orthod 2005, 27: 333–339.View ArticleGoogle Scholar
- Komatsu K, Shibata T, Shimada A, Viidik A, Chiba M: Age-related and regional differences in the stress–strain and stress–relaxation behaviours of the rat incisor periodontal ligament. J Biomech 2004, 37: 1097–1106.View ArticleGoogle Scholar
- Natali AN, Pavan PG, Carniel EL, Dorow C: Viscoelastic response of the periodontal ligament: an experimental–numerical analysis. Connect Tissue Res 2004, 45: 222–230.View ArticleGoogle Scholar
- van Driel WD, van Leeuwen EJ, Von den Hoff JW, Maltha JC, Kuijpers-Jagtman AM: Time-dependent mechanical behaviour of the periodontal ligament. Proc Inst Mech Eng [H] 2000, 214: 497–504.View ArticleGoogle Scholar
- Qian L, Todo M, Morita Y, Matsushita Y, Koyano K: Deformation analysis of the periodontium considering the viscoelasticity of the periodontal ligament. Dent Mater 2009, 25: 1285–1292.View ArticleGoogle Scholar
- Natali AN, Pavan PG, Scarpa C: Numerical analysis of tooth mobility: formulation of a non-linear constitutive law for the periodontal ligament. Dent Mater 2004, 20: 623–629.View ArticleGoogle Scholar
- Wood SA, Strait DS, Dumont ER, Ross CF, Grosse IR: The effects of modeling simplifications on craniofacial finite element models: the alveoli (tooth sockets) and periodontal ligaments. J Biomech 2011, 44: 1831–1838.View ArticleGoogle Scholar
- Verma P, Love RM: A micro CT study of the mesiobuccal root canal morphology of the maxillary first molar tooth. Int Endod J 2011, 44: 210–217.View ArticleGoogle Scholar
- Naveh GRS, Shahar R, Brumfeld V, Weiner S: Tooth movements are guided by specific contact areas between the tooth root and the jaw bone: a dynamic 3D microCT study of the rat molar. J Struct Biol 2012, 177: 477–483.View ArticleGoogle Scholar
- Chen G, Fan W, Mishra S, El-Atem A, Schuetz MA, Xiao Y: Tooth fracture risk analysis based on a new finite element dental structure models using micro-CT data. Comput Biol Med 2012, 42: 957–963.View ArticleGoogle Scholar
- Lin JD, Özcoban H, Greene JP, Jang AT, Djomehri SI, Fahey KP, Hunter LL, Schneider GA, Ho SP: Biomechanics of a bone–periodontal ligament–tooth fibrous joint. J Biomech 2013, 46: 443–449.View ArticleGoogle Scholar
- Cho E, Sadr A, Inai N, Tagami J: Evaluation of resin composite polymerization by three dimensional micro-CT imaging and nanoindentation. Dent Mater 2011, 27: 1070–1078.View ArticleGoogle Scholar
- Dorow C, Krstin N, Sander F-G: Experiments to determine the material properties of the periodontal ligament. J Orofac Orthop 2002, 63: 94–104.View ArticleGoogle Scholar
- Komatsu K, Chiba M: Synchronous recording of load-deformation behaviour and polarized light-microscopic images of the rabbit incisor periodontal ligament during tensile loading. Arch Oral Biol 2001, 46: 929–937.View ArticleGoogle Scholar
- Hurng JM, Kurylo MP, Marshall GW, Webb SM, Ryder MI, Ho SP: Discontinuities in the human bone-PDL-cementum complex. Biomaterials 2011, 32: 7106–7117.View ArticleGoogle Scholar
- Toms SR, Dakin GJ, Lemons JE, Eberhardt AW: Quasi-linear viscoelastic behavior of the human periodontal ligament. J Biomech 2002, 35: 1411–1415.View ArticleGoogle Scholar
- Katona TR, Paydar NH, Akay HU, Roberts WE: Stress analysis of bone modelling response to rat molar orthodontics. J Biomech 1995, 28: 27–38.Google Scholar
- Kawarizadeh A, Bourauel C, Zhang D, Götz W, Jäger A: Correlation of stress and strain profiles and the distribution of osteoclastic cells induced by orthodontic loading in rat. Eur J Oral Sci 2004, 112: 140–147.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.