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.