The Biomechanical Function of Periodontal Ligament Fibres in Orthodontic Tooth Movement

Orthodontic tooth movement occurs as a result of resorption and formation of the alveolar bone due to an applied load, but the stimulus responsible for triggering orthodontic tooth movement remains the subject of debate. It has been suggested that the periodontal ligament (PDL) plays a key role. However, the mechanical function of the PDL in orthodontic tooth movement is not well understood as most mechanical models of the PDL to date have ignored the fibrous structure of the PDL. In this study we use finite element (FE) analysis to investigate the strains in the alveolar bone due to occlusal and orthodontic loads when PDL is modelled as a fibrous structure as compared to modelling PDL as a layer of solid material. The results show that the tension-only nature of the fibres essentially suspends the tooth in the tooth socket and their inclusion in FE models makes a significant difference to both the magnitude and distribution of strains produced in the surrounding bone. The results indicate that the PDL fibres have a very important role in load transfer between the teeth and alveolar bone and should be considered in FE studies investigating the biomechanics of orthodontic tooth movement.     

A visco-elastic model for the prediction of orthodontic tooth movement

This study presents a biomechanical model of orthodontic tooth movement. Although such models have already been presented in the literature, most of them incorporate computationally expensive finite elements (FE) methods to determine the strain distribution in the periodontal ligament (PDL). In contrast, the biomechanical model presented in this work avoids the use of FE methods. The elastic deformation of the PDL is modelled using an analytical approach, which does not require setting up a 3D model of the tooth. The duration of the lag phase is estimated using the calculated hydrostatic stresses, and bone remodelling is predicted by modelling the alveolar bone as a viscous material. To evaluate the model, some typically used motion patterns were simulated and a sensitivity analysis was carried out on the parameters. Results show that despite some shortcomings, the model is able to describe commonly used motion patterns in orthodontic tooth movement, in both single- and multi-rooted teeth.     

The effects of the periodontal ligament on mandibular stiffness: a study combining finite element analysis and geometric morphometrics

It is generally accepted that the periodontal ligament (PDL) plays a crucial role in transferring occlusal forces from the teeth to the alveolar bone. Studies using finite element analysis (FEA) have helped to better understand this role and show that the stresses and strains in the alveolar bone are influenced by whether and how PDL is included in FE models. However, when the overall distribution of stresses and strains in crania and mandibles are of interest, PDL is often not included in FE models, although little is known about how this affects the results. Here we study the effect of representing PDL as a layer of solid material with isotropic homogeneous properties in an FE model of a human mandible using a novel application of geometric morphometrics. The results show that the modelling of the PDL affects the deformation and thus strain magnitudes not only of the alveolar bone around the biting tooth, but that the whole mandible deforms differently under load. As a result, the strain in the mandibular corpus is significantly increased when PDL is included, while the strain in the bone beneath the biting tooth is reduced. These results indicate the importance of the PDL in FE studies. Thus we recommend that the PDL should be included in FE models of the masticatory apparatus, with tests to assess the sensitivity of the results to changes in the Young's modulus of the PDL material.     

Semi-automatic generation of finite element meshes fo dental preparations]

The mechanical properties and elastic behaviour of periodontal tissue are a decisive factor in understanding initial tooth mobility and bone remodelling processes in orthodontics. An experimental set-up was designed to precisely determine a tooth's elastic response to different loading conditions. Segments of pig's maxilla bearing separated molars were used, and their mechanical response to loading was recorded. Subsequently, finite element analysis (FEA) was performed on the basis of the experimental data. The combination of experimental and numerical methods was used to determine the material properties of the periodontal ligament (PDL). The geometries of the preparations were reconstructed and FE meshes generated semi-automatically with the aid of the special computer program, CAGOG (Computer Aided Generator for Orthodontic Geometries) to optimally match the experimental geometry. Nonlinear material parameters were determined for the PDL and verified by comparing experimental and numerical results obtained in other specimens with an error of about 10%. This good correlation indicates that the selected method of mesh generation is appropriate for creating realistic FE models that can be compared with experimental results.     

Investigation of effective intrusion and extrusion force for maxillary canine using finite element analysis

Abstract

Orthodontic tooth movement is mainly regulated by the biomechanical responses of loaded periodontal ligament (PDL). We investigated the effective intervals of orthodontic force in pure maxillary canine intrusion and extrusion referring to PDL hydrostatic stress and logarithmic strain. Finite element analysis (FEA) models, including a maxillary canine, PDL and alveolar bone, were constructed based on computed tomography (CT) images of a patient. The material properties of alveolar bone were non-uniformly defined using HU values of CT images; PDL was assumed to be a hyperelastic–viscoelastic material. The compressive stress and tensile stress ranging from 0.47 to 12.8?kPa and 18.8 to 51.2?kPa, respectively, were identified as effective for tooth movement; a strain 0.24% was identified as the lower limit of effective strain. The stress/strain distributions within PDL were acquired in canine intrusion and extrusion using FEA; root apex was the main force-bearing area in intrusion–extrusion movements and was more prone to resorption. Owing to the distinction of PDL biomechanical responses to compression and tension, the effective interval of orthodontic force was substantially lower in canine intrusion (80–90?g) than in canine extrusion (230–260?g). A larger magnitude of force remained applicable in canine extrusion. This study revised and complemented orthodontic biomechanical behaviours of tooth movement with intrusive–extrusive force and could further help optimize orthodontic treatment.     

The effects of modeling simplifications on craniofacial finite element models: the alveoli (tooth sockets) and periodontal ligaments

Several finite element models of a primate cranium were used to investigate the biomechanical effects of the tooth sockets and the material behavior of the periodontal ligament (PDL) on stress and strain patterns associated with feeding. For examining the effect of tooth sockets, the unloaded sockets were modeled as devoid of teeth and PDL, filled with teeth and PDLs, or simply filled with cortical bone. The third premolar on the left side of the cranium was loaded and the PDL was treated as an isotropic, linear elastic material using published values for Young's modulus and Poisson's ratio. The remaining models, along with one of the socket models, were used to determine the effect of the PDL's material behavior on stress and strain distributions under static premolar biting and dynamic tooth loading conditions. Two models (one static and the other dynamic) treated the PDL as cortical bone. The other two models treated it as a ligament with isotropic, linear elastic material properties. Two models treated the PDL as a ligament with hyperelastic properties, and the other two as a ligament with viscoelastic properties. Both behaviors were defined using published stress-strain data obtained from in vitro experiments on porcine ligament specimens. Von Mises stress and strain contour plots indicate that the effects of the sockets and PDL material behavior are local. Results from this study suggest that modeling the sockets and the PDL in finite element analyses of skulls is project dependent and can be ignored if values of stress and strain within the alveolar region are not required.     

Mechanical responses of the periodontal ligament based on an exponential hyperelastic model: a combined experimental and finite element method

The V–W exponential hyperelastic model is adopted to describe the instantaneous elastic response of the periodontal ligament (PDL). The general theoretical framework of constitutive modeling is described based on nonlinear continuum mechanics, and the elasticity tensor used to develop UMAT subroutine is formulated. Nanoindentation experiment is performed to characterize mechanical properties of an adult pig PDL specimen. Then the experiment is simulated by using the finite element (FE) analysis. Meanwhile, the optimized material parameters are identified by the inverse FE method. The good agreement between the simulated results and experimental data demonstrates that the V–W model is capable of describing the mechanical behavior of the PDL. Therefore, the model and its implementation into FE code are validated. By using the model, we simulate the tooth movement under orthodontic loading to predict the mechanical responses of the PDL. The results show that local concentrations of stress and strain in the PDL are found.     

Experimental model of tooth mobility in the human "in vivo"]

The aim of the present study was to investigate experimentally the mechanical properties of tooth deflection under external loading. These properties have a significant impact on tooth movement during orthodontic treatment. The stresses and strains caused by tooth movement influence bone remodelling, which is the basis of orthodontic treatment. The movement of a tooth as a direct reaction to the forces acting on it is termed "initial" movement. It is nonlinear and has a clearly time-dependent component. While the initial tooth movement represents the totality of the reaction mechanisms of all the tissues of the tooth unit, it is determined primarily by the mechanical properties of the periodontal ligament (PDL). The PDL is the softest tissue of the tooth unit and is therefore subject to the largest deformations when forces act on the crown of the tooth. The objective of orthodontic treatment is to achieve as precise and rapid tooth movement as possible, without provoking such undesired effects as bone and root resorption. To enable the implementation of an optimal orthodontic force system that meets these requirements, a thorough knowledge of the biomechanics of tooth movement is a must.     

Finite element simulation of the behavior of the periodontal ligament: A validated nonlinear contact model

Due to its significance in tooth movement, the stress/deformation field of periodontium and the alveolar bone remodeling process, periodontal ligament (PDL) cannot be excluded from the studies investigating dental biomechanics regarding its excessive deformability. Therefore, many analytical and numerical researches are carried out to simulate its response and to create a constitutive model via experiments intending to discover the material properties of PDL. The aim of this study is to formulate a user specified contact model that can be used in conjunction with finite element (FE) software and reflects PDL’s influence on neighboring structures based on the currently available information, without requiring an actual volumetric finite element mesh of ligament. The results show good agreement with available experimental tooth mobility data. Smooth stress fields are obtained on the tooth root and alveolar bone, which is a significant aspect in bone-remodeling studies. The advantage of simulating PDL as a contact model at the interface of tooth root and the alveolar process instead of a solid-meshed FE model with poor geometric morphology and/or very dense mesh is expected to save pre/post-processing workforce, to increase the accuracy and to contribute to the smoothness of interface stress distributions.     

Determination of orthotropic bone elastic constants using FEA and modal analysis

Finite element models have been widely employed in an effort to quantify the stress and strain distribution around implanted prostheses and to explore the influence of these distributions on their long-term stability. In order to provide meaningful predictions, such models must contain an appropriate reflection of mechanical properties. Detailed geometrical and density information is now readily available from CT scanning. However, despite the use of phantoms, a method of determining mechanical properties (or elastic constants) from bone density has yet to be made available in a usable form.In this study, a cadaveric bone was CT scanned and its natural frequencies were measured using modal analysis. Using the geometry obtained from the CT scan data, a finite element mesh was created with the distribution of density established by matching the mass of the FE bone model with the mass of the cadaveric bone. The maximum values of the orthotropic elastic constants were then established by matching the predictions from FE modal analyses to the experimental natural frequencies, giving a maximum error of 7.8% over 4 modes of vibration. Finally, the elastic constants of the bone derived from the analyses were compared with those measured using ultrasound techniques. This produced a difference of <1% for both the maximum density and axial Young's Modulus. This study has thereby produced an orthotropic finite element model of a human femur. More importantly, however, is the implication that it is possible to create a valid FE model by simply comparing the FE results with the measured resonant frequency of the CT scanned bone.     

Mechanical Response of Bone under Short-term Loading of a Dental Implant with an Internal Layer Simulating the Nonlinear Behaviour of the Periodontal Ligament

We consider a non-standard design for a fixed dental implant, incorporating a soft layer which simulates the presence of the periodontal ligament (PDL). Instead of being aimed at causing an a priori defined stress/strain field within the surrounding bone, upon loading, such a design simply tries to better reproduce the natural tooth–PDL configuration. To do this, the mechanical properties of the internal layer match those of the PDL, determined experimentally to be strongly nonlinear. Three-dimensional finite element analyses show that the presence of such a layer produces (i) a prosthesis mobility very similar to that of a healthy tooth, for several loading conditions, and (ii) a stress/strain distribution substantially different from that arising, upon loading, around a conventional implant. The lack of knowledge of the real mechanical fields existing, under loading, in the bone around a healthy tooth makes it very difficult to state that the stress distribution produced by the modified implant is “better” than that produced by the standard one. Nevertheless, the comparison of the results obtained here, with those of previous refined analyses of the tooth–PDL–bone system, indicates that the modified implant tends to produce a stress distribution in the bone, upon loading, closer to “natural” than that given by the standard one, within the limits imposed by the presence of threads coupling the implant with the bone.     

Comparison of micro-level and continuum-level voxel models of the proximal femur

Continuum-level finite element (FE) models became standard computational tools for the evaluation of bone mechanical behavior from in vivo computed tomography scans. Such scans do not account for the anisotropy of the bone. Instead, local mechanical properties in the continuum-level FE models are assumed isotropic and are derived from bone density, using statistical relationships. Micro-FE models, on the other hand, incorporate the anisotropic structure in detail. This study aimed to quantify the effects of assumed isotropy, by comparing continuum-level voxel models of a healthy and a severely osteoporotic proximal femur with recently analyzed micro-FE models of the same bones. The micro-model element size was coarsened to generate continuum FE models with two different element sizes (0.64 and 3.04 mm) and two different density–modulus relationships found in the literature for wet and ash density. All FE models were subjected to the same boundary conditions that simulated a fall to the side, and the stress and strain distributions, model stiffness and yield load were compared. The results indicated that the stress and strain distributions could be reproduced well with the continuum models. The smallest differences between the continuum-level model and micro-level model predictions of the stiffness and yield load were obtained with the coarsest element size. Better results were obtained for both continuum-element sizes when isotropic moduli were based on ash density rather than wet density.     

Effect of variable periodontal ligament thickness and its non-linear material properties on the location of a tooth’s centre of resistance

In orthodontics, the 3D translational and rotational movement of a tooth is determined by the force–moment system applied and the location of the tooth’s centre of resistance (CR). Because of the practical constraints of in-vivo experiments, the finite element (FE) method is commonly used to determine the CR. The objective of this study was to investigate the geometric model details required for accurate CR determination, and the effect of material non-linearity of the periodontal ligament (PDL). A FE model of a human lower canine derived from a high-resolution µCT scan (voxel size: 50 µm) was investigated by applying four different modelling approaches to the PDL. These comprised linear and non-linear material models, each with uniform and realistic PDL thickness. The CR locations determined for the four model configurations were in the range 37.2–45.3% (alveolar margin: 0%; root apex: 100%). We observed that a non-linear material model introduces load-dependent results that are dominated by the PDL regions under tension. Load variation within the range used in clinical orthodontic practice resulted in CR variations below 0.3%. Furthermore, the individualized realistic PDL geometry shifted the CR towards the alveolar margin by 2.3% and 2.8% on average for the linear and non-linear material models, respectively. We concluded that for conventional clinical therapy and the generation of representative reference data, the least sophisticated modelling approach with linear material behaviour and uniform PDL thickness appears sufficiently accurate. Research applications that require more precise treatment monitoring and planning may, however, benefit from the more accurate results obtained from the non-linear constitutive law and individualized realistic PDL geometry.     

Comparison of strain measurement in the mouse forearm using subject-specific finite element models,strain gaging,and digital image correlation

Mechanical loading in bone leads to the activation of bone-forming pathways that are most likely associated with a minimum strain threshold being experienced by the osteocyte. To investigate the correlation between cellular response and mechanical stimuli, researchers must develop accurate ways to measure/compute strain both externally on the bone surface and internally at the osteocyte level. This study investigates the use of finite element (FE) models to compute bone surface strains on the mouse forearm. Strains from three FE models were compared to data collected experimentally through strain gaging and digital image correlation (DIC). Each FE model was assigned subject-specific bone properties and consisted of one-dimensional springs representing the interosseous membrane. After three-point bending was performed on the ulnae and radii, moment of inertia was determined from microCT analysis of the bone region between the supports and then used along with standard beam analyses to calculate the Young’s modulus. Non-contact strain measurements from DIC were determined to be more suitable for validating numerical results than experimental data obtained through conventional strain gaging. When comparing strain responses in the three ulnae, we observed a 3–14% difference between numerical and DIC strains while the strain gage values were 37–56% lower than numerical values. This study demonstrates a computational approach for capturing bone surface strains in the mouse forearm. Ultimately, strains from these macroscale models can be used as inputs for microscale and nanoscale FE models designed to analyze strains directly in the osteocyte lacunae.     

Influence of the change in stem length on the load transfer and bone remodelling for a cemented resurfaced femur

The effect of a short-stem femoral resurfacing component on load transfer and potential failure mechanisms has rarely been studied. The stem length has been reduced by approximately 50% as compared to the current long-stem design. Using 3-D FE models of natural and resurfaced femurs, the study is aimed at investigating the influence of a short-stem resurfacing component on load transfer and bone remodelling. Applied loading conditions include normal walking and stair climbing. The mechanical role of the stem along with implant–cement and stem–bone contact conditions was observed to be crucial. Shortening the stem length to half of the current length (long-stem) led to several favourable effects, even though the stress distributions in the implant and the cement were similar in both the cases. The short-stem implant led not only to a more physiological stress distribution but also to bone apposition (increase of 20–70% bone density) in the superior resurfaced head, when the stem–bone contact prevailed. This also led to a reduction in strain concentration in the cancellous bone around the femoral neck–component junction. The normalised peak strain in this region was lower for the short-stem design as compared to that of the long-stem one, thereby reducing the initial risk of neck fracture. The effect of strain shielding (50–75% reduction) was restricted to a small bone volume underlying the cement, which was approximately half of that of the long-stem design. Consequently, bone resorption was considerably less for the short-stem design. The short-stem design offers better prospects than the long-stem resurfacing component.     

Experimentally tested computer modeling of stress fractures in rats

The objective of this study was to develop a finite-element (FE) modeling methodology for studying the etiology of a stress fracture (SF). Several variants of three-dimensional FE models of a rat hindlimb, which differed in length or stiffness of tissues, enabling the analyses of mechanical strains and stress in the tibia, were created. We compared the occurrence of SFs in an animal model to validate locations of peak strains/stresses in the FE models. Four Sprague-Dawley male rats, age ~7 wk, were subjected to mechanical cyclic loads of 1.2 Hz and ~6 N, which were delivered to their hindlimb for 30 min, 3 times/wk, up to 12 wk, by using a specially designed apparatus. The results showed that 1) FE modeling predicted the maximal strains/stresses (~220,0 με and ~29 MPa, respectively) between the mid- and proximal thirds of the tibia; 2) in a longer shin, greater and more inhomogeneous tensile strains/stresses were evident, at the same location; 3) anatomical variants in shin length influenced the strain/stress distributions to a greater extent with respect to changes in mechanical properties of tissues; and 4) bone stiffness was more dominant than muscle stiffness in affecting the strain/stress distributions. In the animal study, 35,000 loading cycles were associated with the formation of a SF. The location of the identified SF in the rat limb verified the FE model. We find the suggested model a valuable tool in studying various aspects of SFs.     

Quantification of a rat tail vertebra model for trabecular bone adaptation studies

A feedback controlled loading apparatus for the rat tail vertebra was developed to deliver precise mechanical loads to the eighth caudal vertebra (C8) via pins inserted into adjacent vertebrae. Cortical bone strains were recorded using strain gages while subjecting the C8 in four cadaveric rats to mechanical loads ranging from 25 to 100 N at 1 Hz with a sinusoidal waveform. Finite element (FE) models, based on micro computed tomography, were constructed for all four C8 for calculations of cortical and trabecular bone tissue strains. The cortical bone strains predicted by FE models agreed with strain gage measurements, thus validating the FE models. The average measured cortical bone strain during 25-100 N loading was between 298 +/- 105 and 1210 +/- 297 microstrain (muepsilon). The models predicted average trabecular bone tissue strains ranging between 135 +/- 35 and 538 +/- 138 mu epsilon in the proximal region, 77 +/- 23-307 +/- 91 muepsilon in the central region, and 155 +/- 36-621 +/- 143 muepsilon in the distal region for 25-100 N loading range. Although these average strains were compressive, it is also interesting that the trabecular bone tissue strain can range from compressive to tensile strains (-1994 to 380 mu epsilon for a 100 N load). With this novel approach that combines an animal model with computational techniques, it could be possible to establish a quantitative relationship between the microscopic stress/strain environment in trabecular bone tissue, and the biosynthetic response and gene expression of bone cells, thereby study bone adaptation.     

Finite element and experimental cortex strains of the intact and implanted tibia

Finite Element (FE) models for the simulation of intact and implanted bone find their main purpose in accurately reproducing the associated mechanical behavior. FE models can be used for preclinical testing of joint replacement implants, where some biomechanical aspects are difficult, if not possible, to simulate and investigate in vitro. To predict mechanical failure or damage, the models should accurately predict stresses and strains. Commercially available synthetic femur models have been extensively used to validate finite element models, but despite the vast literature available on the characteristics of synthetic tibia, numerical and experimental validation of the intact and implant assemblies of tibia are very limited or lacking. In the current study, four FE models of synthetic tibia, intact and reconstructed, were compared against experimental bone strain data, and an overall agreement within 10% between experimental and FE strains was obtained. Finite element and experimental (strain gauge) models of intact and implanted synthetic tibia were validated based on the comparison of cortex bone strains. The study also includes the analysis carried out on standard tibial components with cemented and noncemented stems of the P.F.C Sigma Modular Knee System. The overall agreement within 10% previously established was achieved, indicating that FE models could be successfully validated. The obtained results include a statistical analysis where the root-mean-square-error values were always <10%. FE models can successfully reproduce bone strains under most relevant acting loads upon the condylar surface of the tibia. Moreover, FE models, once properly validated, can be used for preclinical testing of tibial knee replacement, including misalignment of the implants in the proximal tibia after surgery, simulation of long-term failure according to the damage accumulation failure scenario, and other related biomechanical aspects.     

Simulation of bone remodelling in orthodontic treatment

Orthodontic treatments not only displace irregular teeth but also induce responses in surrounding bone tissues. Bone remodelling is regarded as the regulatory mechanism triggered by mechanical loading. This study was aimed at investigating the effect of orthodontic loading on both tooth movement and neighbouring bone density distribution. A set of computational algorithms incorporating both external and internal remodelling mechanisms was implemented into a patient-specific 3D finite element (FE) model to investigate and analyse orthodontic treatment under four typical modes of orthodontic loading. The consequence of orthodontic treatment was reproduced numerically by using this FE-based technique. The results indicated that the diverse modes of orthodontic loading would result in different magnitudes of tooth movement and particular morphology of bone density distribution. It is illuminated that the newly developed algorithms may replicate the clinical situation more closely compared with the previous proposed method.     

Oscillatory shear loading of bovine periodontal ligament--a methodological study

This study examined the stress response of bovine periodontal ligament (PDL) under sinusoidal straining. The principle of the test consisted in subjecting transverse tooth, PDL and bone sections of known geometries to controlled oscillatory force application. The samples were secured to the actuator by support plates fabricated using a laser sintering technique to fit their contours to the tooth and the alveolar bone. The actuator was attached to the root slices located in the specimen's center. Hence the machine was able to push or pull the root relative to its surrounding alveolar bone. After determining an optimal distraction amplitude, the samples were cyclically loaded first in ramps and then in sinusoidal oscillations at frequencies ranging from 0.2 to 5 Hz. In the present study the following observations were made: (1) Imaging and the laser sintering technique can be used successfully to fabricate custom-made support plates for cross-sectional root-PDL-bone sections using a laser sintering technique, (2) the load-response curves were symmetric in the apical and the coronal directions, (3) both the stress response versus phase angle and the stress response versus. strain curves tended to "straighten" with increasing frequency, and (4) the phase lag between applied strain and resulting stress was small and did not differ in the intrusive and the extrusive directions. As no mechanical or time-dependent anisotropy was demonstrable in the intrusive and extrusive directions, such results may considerably simplify the development of constitutive laws for the PDL.