Finite element model of load adaptive remodelling induced by orthodontic forces

Anchorage effects of ligation and direct occlusion in orthodontics: A finite element analysis

BACKGROUND AND OBJECTIVE

During orthodontic treatment, the figure-of-eight ligature and the physiological occlusion play an important role in providing anchorage effects. However, their effects on reaction forces of tooth and stress state in periodontal ligament (PDL) have not been quantitatively evaluated yet. In this study, we presented a finite element analysis process for simulating posterior molar ligature and direct occlusion during orthodontics in order to quantitatively assess their anchorage effects.

METHODS

A high precision 3D biomechanical model containing upper and lower teeth, PDL, brackets and archwire was generated from the images of computed tomographic scan and sophisticated modelling procedures. The orthodontic treatment of closing the extraction gap was simulated via the finite element method to evaluate the biomechanical response of the molars under the conditions with or without ligation. The simulations were divided into experimental and control groups. In the experimental group, orthodontic force of 1 N was first applied, then direct occlusal forces of 3 and 10 N were applied on each opposite tooth. While in the control group, occlusal forces were applied without orthodontic treatment. The tooth displacement, the stress state in the PDL and the directions of the resultant forces on each tooth were evaluated.

RESULTS

In the case of molars ligated, the maximum hydrostatic stress in the molars' PDL decreases by 60%. When an initial tooth displacement of several microns occurs in response to an orthodontic force, the direction of the occlusal force changes simultaneously. Even a moderate occlusal force (3 N per tooth) can almost completely offset the mesial forces on the maxillary teeth, thus to provide effective anchorage effect for the orthodontics.

CONCLUSIONS

The proposed method is effective for simulating ligation and direct occlusion. Figure-of-eight ligature can effectively disperse orthodontic forces on the posterior teeth, while a good original occlusal relationship provides considerable anchorage effects in orthodontics.

An in-silico approach to modeling orthodontic tooth movement using stimulus-induced external bone adaptation

Bone remodeling process has been used in orthodontics to treat malposition of teeth in patients by applying stimuli outside of usual everyday loads to promote tooth movement by affecting equilibrium state of the surrounding bone tissue. Accurate modeling of long term orthodontic tooth movement (OTM) is crucial in the field of dental biomechanical research since it allows to predict the behavior and interaction of bone-tooth environment in a non-destructive way, and helps to gain more insight on how exactly tooth motion progresses over time. Existence of such predictive tools might help to avoid the adverse effects of OTM on teeth and the surrounding tissues during this clinical procedure. In this study a new numerical approach to simulating long-term OTM is proposed, that involves external bone adaptation with strain energy density of the bone taken as the stimulus parameter and bone adaptation modeled by nodal movements at the bone-tooth interface using Abaqus UMESHMOTION subroutine. Contrary to conventional re-meshing algorithms, where the mesh of resorbed-apposed bone region is constantly updated and element deletion/creation is performed for each increment, the proposed method only moves nodes without changing the initial mesh topology. For this study, a 3D model of right central maxillary incisor tooth and its surrounding maxillary bone was used for the modeling of OTM for a duration of 1 week. Two test cases were performed and the results from induced tooth motion were investigated. Results indicate tooth movement values that were quite close to clinical values provided in the literature and this method is easily applicable to validate various postulates of OTM via adapting the stimulus-adaption rate relation and patient-specific planning of orthodontic patients as well.

Orthodontic Tooth Movement Studied by Finite Element Analysis: an Update. What Can We Learn from These Simulations?

PURPOSE OF REVIEW

To produce an updated overview of the use of finite element (FE) analysis for analyzing orthodontic tooth movement (OTM). Different levels of simulation complexity, including material properties and level of morphological representation of the alveolar complex, will be presented and evaluated, and the limitations will be discussed.

RECENT FINDINGS

Complex formulations of the PDL have been proposed, which might be able to correctly predict the behavior of the PDL both when chewing forces and orthodontic forces are simulated in FE models. The recent findings do not corroborate the simplified view of the classical OTM theories. The use of complex and biologically coherent FE models can help understanding the mechanisms leading to OTM as well as predicting the risk of root resorption related to specific force systems and magnitudes.

Numerical simulation of optimal range of rotational moment for the mandibular lateral incisor, canine and first premolar based on biomechanical responses of periodontal ligaments: a case study

OBJECTIVES

The objective of this study was to investigate the optimal range of rotational moment for the mandibular lateral incisor, canine and first premolar to determine tooth movements during orthodontic treatment using hydrostatic stress and logarithmic strain on the periodontal ligament (PDL) as indicators by numerical simulations.

MATERIAL AND METHODS

Teeth, PDL and alveolar bone numerical models were constructed as analytical objects based on computed tomography (CT) images. Teeth were assumed to be rigid bodies, and rotational moments ranging from 1.0 to 4.0 Nmm were exerted on the crowns. PDL was defined as a hyperelastic-viscoelastic material with a uniform thickness of 0.25 mm. The alveolar bone model was constructed using a non-uniform material with varied mechanical properties determined based on Hounsfield unit (HU) values calculated using CT images, and its bottom was fixed completely. The optimal range values of PDL compressive and tensile stress were set as 0.47-12.8 and 18.8-51.2 kPa, respectively, whereas that of PDL logarithmic strain was set as 0.15-0.3%.

RESULTS

The rotational tendency of PDL was around the long axis of teeth when loaded. The optimal range values of rotational moment for the mandibular lateral incisor, canine and first premolar were 2.2-2.3, 3.0-3.1 and 2.8-2.9 Nmm, respectively, referring to the biomechanical responses of loaded PDL. Primarily, the optimal range of rotational moment was quadratically dependent on the area of PDL internal surface (i.e. area of PDL internal surface was used to indicate PDL size), as described by the fitting formula.

CONCLUSIONS

Biomechanical responses of PDL can be used to estimate the optimal range of rotational moment for teeth. These rotational moments were not consistent for all teeth, as demonstrated by numerical simulations.

CLINICAL RELEVANCE

The quantitative relationship between the area of PDL internal surface and the optimal orthodontic moment can help orthodontists to determine a more reasonable moment and further optimise clinical treatment.

Dynamic measurement of orthodontic force using a tooth movement simulation system based on a wax model

BACKGROUND

Orthodontic force is often statically measured in general, and only the initial force derived from appliances can be assessed.

OBJECTIVE

We aimed to investigate a technological method for measuring dynamic force using tooth movement simulation.

METHODS

Tooth movement was simulated in a softened wax model. A canine tooth was selected for evaluation and divided into the crown and root. A force transducer was plugged in and fixed between the two parts for measuring force. Forces on this tooth were derived by ordinary nickel-titanium (Ni-Ti) wire, hyperelastic Ni-Ti wire, low-hysteresis (LH) Ti-Ni wire and self-made glass fibre-reinforced shape memory polyurethane (GFRSMPU) wire. These forces were measured after the tooth movement.

RESULTS

The canine tooth moved to the desired location, and only a 0.2 mm deviation remained. The changing trends and magnitudes of forces produced by the wires were consistent with the data reported by other studies. The tooth had a higher moving velocity with ordinary Ni-Ti wires in comparison to the other wires. Force attenuation for the GFRSMPU wire was the lowest (40.17%) at the end of the test, indicating that it provided light but continuous force.

CONCLUSIONS

Mimicked tooth movements and dynamic force measurements were successfully determined in tooth movement simulation. These findings could help with estimating treatment effects and optimising the treatment plan.

Finite element modeling of the periodontal ligament under a realistic kinetic loading of the jaw system

Purpose

The stresses and deformations in the periodontal ligament (PDL) under the realistic kinetic loading of the jaw system, i.e., chewing, are difficult to be determined numerically as the mechanical properties of the PDL is variably present in different finite element (FE) models. This study was aimed to conduct a dynamic finite element (FE) simulation to investigate the role of the PDL (PDL) material models in the induced stresses and deformations using a simplified patient-specific FE model of a human jaw system.

Methods

To do that, a realistic kinetic loading of chewing was applied to the incisor point, contralateral, and ipsilateral condyles, through the experimentally proven trajectory approach. Three different material models, including the elasto-plastic, hyperelastic, and viscoelastic, were assigned to the PDL, and the resulted stresses of the tooth FE model were computed and compared.

Results

The results revealed the highest von Mises stress of 620.14 kPa and the lowest deformation of 0.16 mm in the PDL when using the hyperelastic model. The concentration of the stress in the elastoplastic and viscoelastic models was in the mid-root and apex of the PDL, while for the hyperelastic model, it was concentrated in the cervical margin. The highest deformation in the PDL regardless of the employed material model was located in the caudal direction of the tooth. The viscoelastic PDL absorbed the transmitted energy from the dentine and led to lower stress in the cancellous bone compared to the elastoplastic and hyperelastic material models.

Conclusion

These results have implications not only for understanding the stresses and deformations in the PDL under chewing but also for providing comprehensive information for the medical and biomechanical experts in regard of the role of the material models being used to address the mechanical behavior of the PDL in other components of the tooth.