 Research
 Open Access
 Published:
Establishment, FEM analysis and experimental validation of tooth movement prediction model of orthodontic archwire Tloop
BMC Oral Health volumeÂ 22, ArticleÂ number:Â 406 (2022)
Abstract
Background
The Tloop has been used clinically to close gap between teeth. And it is a typical orthodontic archwire bending method. However, the design of the Tloop parameters for different patients is based on the clinical experience of the dentists. The variation in dentists' clinical experience is the main reason for inadequate orthodontic treatment, even high incidence of postoperative complications.
Methods
Firstly, the tooth movement prediction model is established based on the analysis of the Tloop structure and the waxy model dynamic resistance. As well as the reverse reconstruction of the complete maxillary 3D model based on the patient CBCT images, the oral biomechanical FEM analysis is completed. A maxillary waxy dental model is manufactured to realize the waterbath measurement experiment in vitro mimicking the oral bioenvironment. Thus, the calculated, simulation and experimental data are obtained, as well as obtaining a cloud of total deformation from the simulation analysis.
Results
The growth trend of the 11 sets of simulation data is the same as that of the experimental data. And all of them show that the tooth displacement is positively correlated with the crosssectional size of the archwire, and the clearance distance. As well as the higher Young's modulus of the archwire material, the greater the tooth displacement. And the effect of archwire parameters on tooth displacement derived from simulation and experimental data is consistent with the prediction model. The experimental and calculated data are also compared and analyzed, and the two kinds of data are basically consistent in terms of growth trends and fluctuations, with deviation rates ranging from 2.17 â€‰to â€‰10.00%.
Conclusions
This study shows that the accuracy and reliability of the tooth movement prediction model can be verified through the comparative analysis and deviation calculation of the obtained calculated, simulation and experimental data, which can assist dentists to safely and efficiently perform orthodontic treatment on patients. And the FEM analysis can achieve predictability of orthodontic treatment results.
Background
In recent years, oral malocclusion has shown a high prevalence and is one of the three major oral diseases [1]. It affects oral function, craniofacial development, appearance and interpersonal relationships, and even causes psychological disorders in minors [2, 3, 4]. Orthodontic techniques can effectively treat malocclusion. During orthodontic treatment, due to the phenomenon of orthodontic crowding, some teeth are selectively extracted and a corresponding extraction gap is formed. For this reason, the closure of the extraction gap is a very important aspect [5], which is directly related to the smooth implementation of the entire orthodontic plan [6, 7, 8].
Tloop is widely used to close gap between teeth due to the structural characteristics that can generate horizontal, vertical and torsional force [9, 10, 11, 12, 13]. During the process of closing gap, the teeth are subjected to both the orthodontic force generated by the Tloop and the biological resistance of the periodontal tissuesÂ [14, 15, 16]. In turn, the orthodontic force generated by the archwire are inextricably linked to the clearance distance, crosssectional size and material properties of the Tloop [17, 18].
In traditional orthodontic diagnosis, depending on the treatment stage and specific orthodontic purpose, dentists select archwire parameters empirically [19, 20, 21, 22]. This leads to unpredictable tooth movement, and can easily cause irreparable damage to the patients and reduce treatment efficiency [13]. During the treatment, if the tooth movement is too small, the orthodontic effect is not obvious, and if the tooth movement is too large, it can increase the patient's dental pain, leading to periodontal tissue damage, tooth loosening, and even serious root and alveolar bone resorption problems [23, 24, 25]. Therefore, dentists will find the optimal position by repeatedly adjusting the parameters of the orthodontic archwire. However, repeated adjustments can reduce the mechanical properties of the orthodontic archwire. And it also increases the clinical treatment time and causes inconvenience to both the dentist and the patient. Therefore, by revealing the relationship between archwire parameters, orthodontic force and tooth movement, it enables the dentist to develop an individualized orthodontic treatment plan by taking into account the patient's malocclusion, the mechanical behavior of archwire and the tooth movement pattern. This further improves the efficiency and safety of fixed orthodontic treatment and reduces the patient's pain.
In this paper, based on the dynamic resistance model of the waxy model and the dynamic orthodontic force model of the Tloop, a tooth movement prediction model is established to realize the tooth movement prediction analysis. The Tloop parameters and tooth movement distance are quantified. And FEM analysis and experimental validation are both conducted. Experimental data are collected using the designed tooth movement measurement method. And the accuracy of the tooth movement prediction model is demonstrated by data comparison and deviation analysis.
Therefore, the tooth movement prediction model established in this paper can assist dentists in selecting the appropriate Tloop parameters for the movement distance and angle of the tooth. Thus, it provides a theoretical basis for the development of individualized orthodontic treatment plans for patients. At the same time, it reduces the clinical treatment time of patients and dentists, improves orthodontic efficiency, as well as achieves predictability of orthodontic treatment effect.
Methods
Dynamic resistance analysis of waxy model
In order to restore the relative tooth movement under realistic conditions, a nonrigid waxy model is chosen to simulate the biological properties of the periodontal environment. Nonrigid waxy models are widely used in dental clinical experiments and training because the viscous fluid properties vary with temperature, such as the Typodon orthodontic training model [26, 27]. As well as Liu et al. experimentally verified the feasibility of the waxy model for simulating tooth movement at a certain temperature [28]. Meanwhile, according to the theory of viscous fluid dynamics, during orthodontic treatment, the impedance force of tooth in waxy model is composed of resistance and inertia force [29]. Through experimental verification, the optimal waterbath temperature for tooth movement in the waxy model is determined to be 75Â Â°C [30, 31]. At the same time, an expression for the variation of the waxy dental density with time is established. According to the expression, the resistance equation and the inertia force equation, the dynamic resistance model of the waxy dental model in the process of simulated tooth movement is established [32, 33], as shown in Fig.Â 1.
Modeling of dynamic orthodontic force in Tloop
During orthodontic treatment, the tooth movement is influenced by a combination of the periodontal tissue's biological resistance and the orthodontic force generated by the archwire. Besides, orthodontic force changes dynamically as treatment progresses [32]. Therefore, before establishing a tooth movement prediction model, a dynamic orthodontic force model needs to be established.
Analysis of Tloop working as follows: the dentist pulls the horizontal arm of Tloop to create a clearance, causes elastic deformation of the vertical arm and the arc portion, and then fixes the horizontal arm on the bracket. The tooth is pulled to a predetermined position using the restoring force of the Tloop itself, thus completing the closure of the gap between teeth [19, 20]. At the same time, Tloop with symmetrical structure is used, that is, the Alpha/Beta structure on both sides of Tloop is equal, so that the orthodontic force released by the horizontal arm can be equal and opposite [34]. Therefore, only one side of Tloop needs to be analyzed. The Tloop deformation analysis is shown in Fig.Â 2.
The orthodontic force generated by the Tloop is composed of the vertical arm and the arc portion of the deformation superimposed. So the vertical arm and the arc portion are modeled mechanically separately. And then the total orthodontic force of the Tloop is calculated by superposition. Firstly, the vertical arm is subjected to force analysis and mechanical modeling. As shown in Fig.Â 2, the bending radius of the arc portion is R, the total height is h, the length of the vertical arm is y. By pulling the horizontal arm of the Tloop creates a clearance distance b. And the component of the force applied to the horizontal arm that deforms the vertical arm is p.
The approximate differential equation for the vertical arm deflection curve is as follows:
where E is the elastic modulus of material, M(x) is the bending moment of the vertical arm at x distance. I_{z} is the rotational inertia of the archwire zaxis. The rotational inertia of a round archwire is I_{z}â€‰=â€‰Ï€D^{4}/64, and D is the diameter of the round archwire. For a rectangular archwire I_{z}â€‰=â€‰c_{1}c_{2}^{3}/12. The c_{1} is length of the side perpendicular to the zaxis on the rectangular archwire crosssection. The rotation angle equation Î¸(x) and the deflection equation Î½(x) for the vertical arm can be obtained by integrating Eq.Â (1) once and twice, respectively.
where C_{0} and D_{0} are integration constants that can be determined by the boundary conditions at special points. The expression for the bending moment of the vertical arm is:
Substituting Eqs. (1) and (3) into Eq.Â (2) for integration, obtaining:
Using the special points with known deflection and rotation angle, the integration constants in Eq.Â (4) can be determined. Since the longitudinal symmetry plane exists at the intersection of horizontal arm of the arc portion and vertical arm, and the symmetry plane is subject to external force, the axis of the deformed curved beam is still located within the longitudinal symmetry plane. The intersection point deformation belongs to the plane bending deformation of the curved beam. Therefore, the arc at the intersection can be equated to a curved beam with radius, curvature and section diameter. To take the radius microelement dÎ± cross section, as shown in Fig.Â 3.
After deformation, the differential equation of the deflection curve at the intersection of horizontal arm of the arc portion and vertical arm is as follows:
where u is the displacement of the curved beam in the xdirection intersection cross section. The bending moment condition at the intersection is known to be \(M_{0} = \left. M \right_{x = 0} =  Py\). \(I_{\omega }\) is the rotational inertia of the curved beam's cross section on the Ï‰ axis. The vertical arm is consistent with the bending form of the curved beam at the intersection, so \(I_{\omega } = I_{z}\). Also it is defined the arc length equation dsâ€‰=â€‰RdÎ±, and obtained the differential equation of the deflection curve at the intersection point after substitution, as shown in Eq.Â (6):
The nonassociative differential equation with constant coefficients for the deflection curve after deformation of the curved beam at the intersection point is:
The boundary condition at the intersection is \(u{}_{{\alpha { = }\pi /2}} = 0\), \(\left. {\frac{{{\text{d}}u}}{{{\text{d}}\alpha }}} \right_{\alpha = \pi /2} = 0\). Therefore, it can be determined that A_{0}â€‰=â€‰0, \(B_{0} { = }\frac{{M_{0} R^{2} }}{{EI_{\omega } }}\), and bring them into Eq.Â (7) to obtain Eq.Â (8):
Also, the corner equation of the curved beam can be obtained as:
By introducing of boundary conditions \(\left. v \right_{x = 0} = \left. u \right_{\alpha = 0} =  \frac{{M_{0} R^{2} }}{{EI_{\omega } }}\), \(\left. \theta \right_{x = 0} = \left. \beta \right_{\alpha = 0} = \frac{{M_{0} R^{{}} }}{{EI_{\omega } }}\), to obtain \(C_{0} = \frac{{M_{0} R^{{}} }}{{EI_{\omega } }}\), \(D_{0} =  \frac{{M_{0} R^{2} }}{{EI_{\omega } }}\).
Bringing \(C_{0}\) and \(D_{0}\) into Eqs. (4) and (10), obtaining:
According to the established deflection differential equation, the corresponding boundary conditions are replaced: the maximum deflection value m, and the maximum rotation angle at xâ€‰=â€‰y. The external force that deforms the vertical arm is then obtained as:
Further according to the principle of action and reaction force, the force P that deforms the vertical arm is equal to the counter force of the orthodontic force F_{1} applied to the tooth during the orthodontic procedure. Thus, it can be obtained that:
As mentioned earlier, the arc portion and the deformation of the vertical arm together constitute the orthodontic force generated by the Tloop during correction. After the analysis of the vertical arm, a mechanical modeling.
analysis of the arc portion is performed. Similar to the analysis of the vertical arm, the horizontal arm of the arc portion as the main deformation part, need to establish the corresponding deformation equation, as shown in Fig.Â 3.
During the deformation of the Tloop, the position of the connection between the horizontal arm of the arc portion and the vertical arm is constantly changing. It is not easy to measure the distance that the horizontal arm of the arc portion moves. Therefore, it is necessary to simplify the movement of the intersection symmetry center during the Tloop deformation. The symmetry center of the intersection before deformation is overlapped, and the difference between the theoretical length and the initial length of the vertical arm after deformation is calculated. This difference is the bending deformation s of the horizontal arm of the arc portion, which can be obtained from the following equation:
where y is the length of the vertical arm in the natural state and also the bending deflection of the vertical arm. Similar to Eq.Â (1), the deflection curve differential equation for the horizontal arm of the arc portion is:
The corner Eq.Â (15) and the deflection Eq.Â (16) for the horizontal arm of the arc portion can be determined by integrating Eq.Â (14) once and twice, respectively:
where C_{1} and D_{1} are integration constants, and the bending moment expression for the horizontal arm of the arc portion is:
where G is the component that deforms the horizontal arm of the arc portion, w is the total length of the horizontal arm of the arc section, and R is the bending radius of the arc portion. Substituting Eqs. (14) and (17) into Eqs. (15) and (16) for integration, the following equation is obtained:
The lateral arc on one side can be equated to a curved beam of radian Ï€/4, and the radian dÎ² is taken to be infinitesimal. This is the same process as solving for the vertical arm boundary condition, so the boundary condition can be given as:
The boundary conditions of the arc portion are \(\left. u \right_{\beta = \pi } = 0\), \(\left. {\frac{{{\text{d}}u}}{{{\text{d}}\beta }}} \right_{\beta = \pi } = 0\), and the general solutions obtained are \(A_{1} = \frac{{M_{0} R^{2} }}{{EI_{\omega } }}\), \(B_{1} = 0\). Then the deflection equation of the arc portion is:
The corner equation of the arc portion is:
Because the boundary conditions are \(\left. v \right_{{l{ = 0}}} = \left. u \right_{\beta = 0} = 0\) and \(\left. \theta \right_{l = 0} = \left. \varepsilon \right_{\beta = 0} = 0\), it is able to obtain C_{1}â€‰=â€‰0 and D_{1}â€‰=â€‰0. Substituting C_{1}â€‰=â€‰0 and D_{1}â€‰=â€‰0 into Eqs. (18) and (19) obtains Eqs. (23) and (24):
The corresponding boundary conditions are substituted as follows: when lâ€‰=â€‰wR, the horizontal arm of the arc portion reaches the maximum deflection at its end, which is vâ€‰=â€‰s, approximated by the same Eq.Â (13). Further can get the bending force of the horizontal arm of the arc portion, as follows:
From further mechanical analysis, it is obtained that the component force along the horizontal arm should be equal to the orthodontic force F_{2} released by the arc portion. The following equation is obtained:
From the mechanical analysis, the restoring force from the deformation of the vertical arm and the arc portion together constitute the orthodontic force F_{0} released by the horizontal arm. That is, the following equation:
Establishment of tooth movement prediction model
After the previous calculations, the dynamic resistance model for the tooth movement simulated in the waxy model can be expressed by Eq.Â (28):
where f is the dynamic resistance to tooth movement in the waxy model. Bringing in the orthodontic force F_{0} released by the Tloop, it is obtained that:
In physics, the relationship between the displacement r and the velocity vector v, Eq.Â (30), can be determined by integrating:
In Newton's Second Law, the relationship between the acceleration a_{0}, the mass m_{0} and the force F can be expressed by Eq.Â (32):
The acceleration a_{0} of the tooth moving in the waxy model, subject to orthodontic force F_{0}, inertial force f_{L} and resistance f_{D}, can be expressed by Eq.Â (32):
By substituting Eq.Â (33) into (31), the tooth movement prediction model can be obtained. The following equation can be obtained:
where S_{m} is the amount of tooth movement under the action of the Tloop in the waxy model, m_{0} is the mass of the tooth, and t_{1} is the time of tooth movement.
Tooth movement FEM analysis
In this study, a 3D model of the complete maxillary bone, dental row and periodontal ligament is reconstructed in reverse by using CBCT images (EWOOVATECH Co., Ltd., Korea, scan thickness 0.5Â mm, 301 scanned images) based on the patient's skull. And according to the clinical orthodontic treatment process, a 3D model with brackets and Tloop is installed on the reconstructed dental. At last, the orthodontic biomechanical FEM analysis is performed for the complete assembly model. The flowchart of the reverse 3D reconstruction model is shown in Fig.Â 4. And the simulation cloud is shown in Fig.Â 5. In the simulation process, the archwire is adopted of stainless steel archwire and Australian archwire, which are commonly used in clinical orthodontics. In addition, the material properties of the brackets and maxillary components are shown in Table 1 [34, 35, 36].
The simulation parameters shown in Fig.Â 5 are: (1) The archwire is a 0.014inch round Australian archwire. (2) The material of the bracket is stainless steel. (3) The top of the maxilla is fixedly constrained. (4) The contact between the archwire and the brackets, between the brackets and the teeth is bound. (5) The archwire is meshed with 0.1Â mm size and the others are meshed with 1.0Â mm, so that the deformation of the archwire can be calculated more accurately. And the total number of nodes and cells of the mesh is 235,827 and 138,031 respectively. In addition, the displacement load is defined as follows: set 10 analysis steps so that the clearance distance grows linearly with 0.1 of the step length, and take the clearance distance as 0.3â€“1.2Â mm for FEM analysis. And the maximum displacement of NO. 13 tooth is 1.4575Â mm at the cusp and the minimum displacement is 0.2512Â mm at the root, which meets the clinical orthodontic requirements [37, 38].
Tooth movement experimental measurement method
In this study, the waxy model is used to simulate the periodontal tissue structure. Firstly, alginate is used as the impression material, and the separated teeth are placed into the alginate negative mold, then the molten waxy base is poured, left to cool and removed. Thus, the preparation of waxy dental model for experimental use is completed. Due to the viscoelastic properties at the corresponding temperature, the waxy model under the action of external force will produce a corresponding elastic deformation [28]. Finally, the teeth on the waxy model will complete the simulation of the orthodontic process with the change of time. The experimental system composition and the waxy dental model are shown in Fig.Â 6. Also detailed waxy model preparation process and modeling methods have been described in the study of Huang et al. [32].
The position of each maxillary teeth is marked using the FDI (FÃ©dÃ©ration Dentaire Internationale) marking method, and maxillary teeth No. 13 and No. 15 are selected as the teeth for the experimental study. Before clinical orthodontic treatment, a total of four teeth in the maxilla and mandible are usually extracted to allow enough space for remaining teeth to move. In the experiment, tooth No.14 is extracted and the two horizontal arms of the Tloop are ligated into the brackets of teeth No. 13 and No. 15. The highest points of the marginal crests of the two teeth are marked, as shown in Fig.Â 6. The distance between the two marker points is measured before the waterbath experiment, and then the waxy model is immersed in a constant temperature waterbath at 75Â Â°C for 2Â min. After taking out, the distance between the two marker points is measured again with a vernier caliper [32]. The difference of the distance between the two marker points before and after the waterbath is sought as the experimental value of the tooth movement distance under the action of Tloop.
Results
In clinical treatment, the magnitude of orthodontic force is changed by adjusting the characteristic parameters of the Tloop. And the magnitude of orthodontic force has an important influence on the tooth movement, the periodontium disease, tooth root resorption, and the alveolar bone reconstruction [39]. Therefore, in this study, the accuracy and reliability of the tooth movement prediction model is verified by comparing and analyzing the calculated, simulation and experimental data. In the FEM analysis, 11 kinds of Tloop 3D models with different parameters are established based on clinical orthodontic treatment. The loading process of the FEM analysis is the same as described before, with the aim of investigating the effects of different clearance distances, crosssectional and materials on tooth movement. To be easily represented, the orthodontic archwires are codenamed as follows: IS1616, IIS1622, IIIS1625, IVS0014, VS0016, VIS0018, VIIS0020, VIIIA0014, IXA0016, XA0018, XIA0020. Where S is stainless steel archwire, A is Australian archwire, 1616 indicates 0.016â€‰Ã—â€‰0.016 inch rectangular archwire, and 0014 indicates 0.014 inch round archwire. Similarly, the 11 kinds of Tloop bending parameters used in the experimental are the same as those of the 3D model, and 10 measurement points are set at 0.1Â mm intervals to collect data during the loading measurement stage of the experiment. The calculated, simulation and experimental data are shown in Table 2. As well, the total strain clouds for the biomechanical FEM analysis of tooth movement under the action of 11 kinds of Tloop are shown in Fig.Â 7.
In this paper, the simulation data and experimental data are compared and analyzed, so that the accuracy and reliability of simulation analysis and experimental measurement can be verified with each other.
FigureÂ 8 shows the effect of rectangular archwire with different crosssectional widths (Codenamed as I, II, and III) on the variation of tooth displacement with Tloop clearance distance. Both simulation and experimental results shows that the tooth displacement is positively correlated with the crosssectional width for the same clearance distance. Similarly, the tooth displacement is positively correlated with the clearance distance for the same crosssectional width. In addition, the effect of crosssectional width on tooth displacement increases as the clearance distance increases.
Comparing different round crosssectional diameter (Codenamed as IV, V, VI, and VII), the effect of tooth displacement with the variation of Tloop clearance distance is shown in Fig.Â 9. Both simulation and experimental results show that the tooth displacement is positively correlated with the crosssectional diameter of the Tloop for the same clearance distance, which is similar to the conclusion for rectangular archwire. From the simulation and experimental results, although the crosssectional diameter of the archwire affects the tooth displacement, its effect is weaker than that of the clearance distance of the Tloop.
In clinical orthodontic treatment, stainless steel archwire and Australian archwire are mainly used. Therefore, two sets of results based on material are divided. The data of four identical archwire types in the two sets are compared to investigate the effect of different archwire materials on tooth displacement [5, 32]. The stainless steel archwires are codenamed as IV, V, VI, and VII, and the Australian archwires are codenamed as VIII, IX, X, and XI. As shown in Fig.Â 10, both the simulation and experimental data show that the effect of tooth displacement with the same crosssectional dimensions is more pronounced with increasing clearance distance using stainless steel archwire than Australian archwire. In addition, this conclusion became more evident as the diameter of the archwire increased.
While mutually verifying the accuracy and reliability of the simulation analysis and experimental measurements, this paper analyses the calculated and experimental data, and performs deviation calculations to verify the accuracy and validity of the tooth movement prediction model. The calculated data is shown in Table 2. And in Figs. 11, 12, 13, the experimental and calculated data, and deviation rates are shown by solid, dashed and dotted lines, respectively.
A comparison of experimental and calculated data of tooth displacement versus clearance distance for rectangular stainless steel archwire of different crosssectional widths (Codenamed as I, II and III) is shown in Fig.Â 11. The obtained calculated and experimental data have the same trend and the deviation varies from 2.17 â€‰to â€‰8.86%. When other conditions are the same, the trend of the curves shows that the calculated data of tooth displacement are positively correlated with the crosssectional width, and the effect of the crosssectional width on tooth displacement increases with the increase of the clearance distance.
For the round archwire, the comparison of experimental and calculated data of tooth displacement for different crosssectional sizes of stainless steel archwire (Codenmed as IV, V, VI, VII) and Australian archwire (codenamed as VIII, IX, X, XI) at different clearance distances, and deviation rates are shown in Figs. 12 and 13. The trends of the calculated and experimental data are consistent in different archwire materials. The deviation rate of calculated and experimental data for stainless steel archwire ranges from 2.36 â€‰to â€‰10.00%, while the deviation rate of Australian archwire ranges from 2.47 â€‰to â€‰9.38%. For the two kinds of orthodontic archwire with different materials, the larger the crosssectional size of the archwire, the greater the tooth displacement effect produces when the clearance distance is the same. When the crosssectional size are all the same, the tooth displacement effect produced by the round orthodontic archwire is positively correlated with the clearance distance.
Discussion
In order to reveal the tooth movement characteristics under different force during clinical orthodontic procedure, researchers have studied it through experimental measurements, FEM analysis and theoretical derivation. For example, researchers have established FEMs of tooth movement under different orthodontic force [40, 41, 42, 43]. During the simulation, the magnitude and direction of the orthodontic force is artificially set without reference, which is not consistent with the clinical reality. In addition, periodontal tissues cannot be effectively differentiated by materials, so corresponding experimental validation should be performed. However, experimental materials are mostly rigid models, such as plaster models or metal models [44], which lack biological properties. In addition, the orthodontic process is a dynamic one, and the orthodontic force will become weaker with the reduction of archwire deformation. Therefore, the real orthodontic effect cannot be effectively reproduced and the prediction function is not yet perfect. In contrast, the Tloop is a common orthodontic arch for clinical treatment of such oral malocclusion as closing the gap between teeth. Although the relationship between archwire shape and orthodontic force for Tloop has been studied to some extent [45, 46, 47, 48], but the archwire bending parameter and orthodontic force relationship, orthodontic force and tooth movement relationship have not been established. Therefore, the current studies still cannot assist dentists to effectively develop individualized Tloop parameter designs for different patients, as well as to accurately predict orthodontic treatment effect.
Thus, this paper firstly establishes a tooth movement prediction model based on the Tloop analysis and the waxy model dynamic resistance. Subsequently, to verify the accuracy of the prediction model, a reverse reconstruction model of the maxillary and dentition based on the patient's CBCT images is used to achieve a biomechanical FEM analysis of the tooth movement. Finally, a waxy dental model is established, and a waterbath measurement experiment mimicking the biological environment of the oral cavity is realized. Simulation analysis and experimental measurements are compared to verify each other's accuracy and reliability. And as shown in Figs. 8 and 9, the conclusion about the influence of different archwire parameters on tooth movement derived from the simulation analysis and experimental measurements is consistent with the established tooth movement prediction model. This further illustrates the accuracy of the simulation analysis and experimental measurements.
Afterwards, the prediction model and experimental measurements are compared and analyzed, and the deviation is calculated, and the total deviation range between the calculated and experimental data is obtained: 2.17 â€‰~â€‰10.00%, which is within the acceptable range and is approved by the dentists of Peking University School of Dentistry. In addition, it can be determined from the trends shown in Figs. 11, 12 and 13 that the movement of the target tooth is positively correlated with both square and round archwire in terms of clearance distance and crosssectional parameters. And the tooth displacement in the stainless steel archwire move faster than those in the Australian archwire as the clearance distance increases. The tooth displacement of the Australian archwire is milder for the same Tloop parameters. These trends and correlations are corrected with the existing academic studies in this field.
Conclusions
In this paper, based on the analysis of the Tloop structure and the waxy model dynamic resistance, a tooth movement prediction model is established. The relationship between Tloop parameters, orthodontic force and tooth movement are revealed and quantified. The validity and reliability of the simulation analysis and experimental measurements, as well as the accuracy of the tooth movement prediction model, are verified by comparing and analyzing the simulation analysis and experimental measurements with each other, as well as calculating and analyzing the deviation of the experimental and calculated data.
Therefore, it can be affirmed that the established tooth movement prediction model has a high degree of confidence and can assist dentists in selecting the appropriate archwire parameters when bending Tloop to obtain the desired tooth movement distance. At the same time, the tooth movement prediction model provides dentists with an effective theoretical guide to avoid inefficient bending of archwire and improve the treatment effect of malocclusion. Moreover, the process and results of the FEM analysis enable predictive orthodontic treatment effects, thus laying the foundation for the realization of digital orthodontic treatment.
In future work, other detailed factors that may affect tooth movement during orthodontic treatment will be further explored. Such as the relative static friction and relative sliding friction between the orthodontic archwire and brackets, and how to choose the right archwire parameters to achieve the best orthodontic effects and speed up the treatment process.
Availability of data and materials
All relevant data are within the present paper; however, the ethics committee strictly forbids sharing the CBCT files of patients and control subjects that they contain some identifying information, according to its regulations on the data access. Therefore, CBCT data of the samples will not be shared. Also the reverse reconstructed 3D models in this paper are belong to the authors and are therefore available only upon request, after approval by all the authors. And if required, please contact the corresponding author of this paper (Correspondence: jiangjingang@hrbust.edu.cn).
Abbreviations
 FEM:

Finite element method
 3D:

Three dimensional
 FDI:

FÃ©dÃ©ration dentaire internationale
References
Jiang JG, Huang ZY, Qian W, et al. Registration technology of augmented reality in oral medicine: a review. IEEE Access. 2019;7:53566â€“84.
He L, Zhang JT. Progress in research on opening and closing malformation. Electr J Gen Stomatol. 2019;6(27):11â€“2.
Yao H, Zhao EY, Qin JY. Research progress on the prevalence of malocclusion. China Primary Health Care. 2003;17(10):74.
Zeng J, Zheng DH, Wang XH, et al. Effects of malocclusion on the psychological status of adults. J Shandong Univ. 2014;52(6):94â€“7.
Jiang JG, Ma XF, Zhang YD, et al. Prediction model and examination of open vertical loop orthodontic force. Arab J Sci Eng. 2019;44(2):1489â€“99.
Liu CJ, You WJ, Chen QF. Comparison of the effects of two different ligation methods on the anterior and posterior arch planes when the Tloop is closed in the active wing treatment technique. Prim Med Forum. 2017;28:3877â€“8.
Wang S, Qin K. Comparison of two methods of closing the extraction space using straight wire appliance. Chin J Pract Stomatol. 2011;4(6):359â€“61.
Wu HP. Clinical study on technical activities wing correction closing tooth extraction gap of. Chin J Aesthet Med. 2015;24(17):42â€“5.
Chen YM, Li YR, Zhang C. Comparison the torque of upper incisors using two closing loops in the treatment of classII division 1 malocclusion. Hebei Medicine. 2007;13(6):655â€“7.
Huang CY. Treatment of 21 cases of adult Anâ€™s II~1 malocclusion with high and low Tloop. Guangxi Med. 2010;32(3):293â€“4.
Zhang YL. Observation of clinic effect of tbending on the correction of incisors crossbite during mixed dentition period. J Kunming Med Univ. 2008;30(2):206â€“7.
Jiang JG, Ma XF, Zuo SH, et al. Digital expression and interactive adjustment method of personalized orthodontic archwire for robotic bending. J Adv Mech Design Syst Manuf. 2019;13(2):18â€“59.
Deng H, Xia ZY, Weng SK, et al. Motion planning and control of a robotic system for orthodontic archwire bending. IEEE/RSJ Int Conf Intell Robots Syst. 2015;2015:3729â€“34.
Teruko TY, Kiyo S, Goudarzi F, et al. Synergistic acceleration of experimental tooth movement by supplementary highfrequency vibration applied with a static force in rats. Sci Rep. 2017;7(1):13969.
Tochigi K, Oda S, Arai K. Influences of archwire size and ligation method on the force magnitude delivered by nickeltitanium alloy archwires in a simulation of mandibular right lateral incisor linguoversion. Dent Mater J. 2015;34(3):388â€“93.
Ayumi SM, Takehito O, Mikihito H, et al. Osteocyte regulation of orthodontic forcemediated tooth movement via RANKL expression. Sci Rep. 2017;7(1):8753.
Pascal G, Laurence J. Stiffness in bending of a superelastic NiTi orthodontic wire as a function of crosssectional dimension. Angle Orthod. 2004;74(5):691â€“6.
Cai YQ. Finite element analysis of archwire parameters and activation forces on the M/F ratio of vertical, L and Tloops. BMC Oral Health. 2020;20(1):70.
Arora V, Sharma R, Chowdhary S. Comparative evaluation of treatment effects between two fixed functional appliances for correction of Class II malocclusion: a singlecenter, randomized controlled trial. Angle Orthod. 2018;88(3):259â€“66.
Zhang JY. Orthodontic treatment of dentition sparse and repair of small teeth in 8 cases. Guangxi Med. 2009;31(2):291â€“2.
Jiang JG, Ma XF, Zhang YD, et al. Study on threedimensional digital expression and robot bending method of orthodontic archwire. Appl Bion Biomech. 2018;2018:1â€“10.
Duan YF, Xiong L, Feng L, et al. Retrospective study of 100 recurrent patients after orthodontics. Stomatology. 2005;25(3):156â€“7.
Li XG, Liu SC, Xiao YY, et al. Different orthodontic forces on the root development of rat molars. Chin J Conserv Dent. 2018;28(8):450â€“4.
Ya M, Eri K, Justin F, et al. Excessive orthodontic force induces odontoclastogenesis in root cementum as evidenced by the expression of RANK/RANK. Int J OralMed Sci. 2011;8(1):38â€“46.
Wu JL, Liu YF, Wei P, et al. A biomechanical case study on the optimal orthodontic force on the maxillary canine tooth based on finite element analysis. J Zhejiang Univ Sci B. 2018;19(7):535â€“46.
Gu M, Li JZ, Liang R, et al. The Typodont investigation of the effect on orthodontic intrusion of upper anterior teeth. J Clin Stomatol. 2004;20:435â€“7.
Nunez DW, Taleghani M, Wathen WF, et al. Typodont versus live patient: predicting dental studentsâ€™ clinical performance. J Dent Educ. 2012;76:407â€“13.
Liu YF, Zhang Q, Zhou D, et al. Mechanical properties and resistance characteristic of orthodontic wax models. J Funct Mater Dev. 2015;21(6):181â€“7.
Zheng M, Du W. Phase behavior, conformations, thermodynamic properties, and molecular motion of multicomponent paraffin waxes: a Raman spectroscopy study. Vib Spectrosc. 2006;40(2):219â€“24.
Higa RH, Henriques JFC, Janson G, et al. Force level of small diameter nickeltitanium orthodontic wires ligated with different methods. Prog Orthod. 2017;18(1):21â€“4.
Lin XP. Basic biomechanical principles of orthodontic treatment. Anal Biomech Clin Orthodontics. 2012;5(1):8â€“9.
Jiang JG, Huang ZY, Ma XF, et al. Establishment and experiment of utility archwire dynamic orthodontic moment prediction model. IEEE Trans Biomed Eng. 2019;67(7):58â€“68.
Jiang JG, Chen HJ, Huang ZY, et al. Orthodontic force prediction model of Tloop closing spring based on dynamic resistance model. Proc Inst Mech Eng [H]. 2020;234(7):095441192094343.
Laviana A. The effect of different Î± and Î² gable on stainless steel and TMA segmental Tloop on the resulting force â€“ an in vitro study. IOP Conf Series: Mater Sci Eng. 2019;550: 012039.
Spyridon NP, Ludger K, Istabrak H, et al. Effect of material variation on the biomechanical behaviour of orthodontic fixed appliances: a finite element analysis. The Eur J Orthodontics. 2015;38(3):300â€“7.
Zhang XF, Wang C, Xia X, et al. Threedimensional finite element biomechanical analysis of lingual inward maxillary anterior teeth with the fine wire arch technique. West China J Stomatol. 2015;33(3):291â€“5.
Ren Y, Maltha JC, Hof MAV, et al. Optimum force magnitude for orthodontic tooth movement: a mathematic model. Am J Orthod Dentofac Orthop. 2004;125(1):71â€“7.
Guan LY, Wu LP. Photoelastic study and clinical application of orthodontic force and tooth movement. J Jilin Univ Technol. 1997;27(2):5.
Yao S, Zhang LJ, Wang GP, et al. Comparison of mechanical properties of different orthodontic wires and their clinical significance. J Mod Denty. 2000;14(2):92â€“3.
Li QH, Wang YR, Xiong GP. Construction of 3D finite element model using the customized lingual bracket system during retracting maxillary anterior teeth. J Clin Stomatol. 2016;32(8):462â€“4.
Wu J, Liu Y, Wei P, et al. Initial orthodontic force analysis of shape memory polymer archwire based on finite element simulation. Chin J Biomed Eng. 2016;35(2):202â€“10.
Cai Y, Yang X, He B, et al. Finite element method analysis of the periodontal ligament in mandibular canine movement with transparent tooth correction treatment. BMC Oral Health. 2015;15(1):106â€“17.
Viecilli RF, Budiman A, Burstone CJ. Axes of resistance for tooth movement: does the center of resistance exist in 3dimensional space. Am J Orthod Dentofac Orthop. 2013;43(2):163â€“72.
Fuck LM, Drescher D. Force systems in the initial phase of orthodontic treatmenta comparison of different leveling arch wires. J Orofac OrthopFortschritte Der Kieferorthop. 2006;67(1):6â€“18.
Jiang JG, Ma XF, Han YS, et al. Experimentation and simulation of second sequential loop orthodontic moment prediction modeling. IEEE Access. 2018;6:56258â€“68.
Ding R, Man YN, Wu PL, et al. Establishment of a threedimensional finite element model of transmission straight wire appliance. Chin J Tissue Eng Res. 2013;17(33):5942â€“8.
Kinzinger G, SyrÃ©e C, Fritz U. Molar distalization with different pendulum appliances: in vitro registration of orthodontic forces and moments in the initial phase. J Orofac Orthop. 2004;65(5):389â€“409.
Lapatki BG, Bartholomeyczik J, Ruther P, et al. Smart bracket for multidimensional force and moment measurement. J Dental Res. 2007;86(1):73â€“8.
Acknowledgements
This study was supported by Heilongjiang Provincial Natural Science Foundation of China (Grant No. LH2021E081), the Fundamental Research Foundation for Universities of Heilongjiang Province (Grant No. LGYC2018JQ016).
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
All authors made substantial contributions to the present study. JGJ, LY, XFM and YFG contributed to conception and design, acquisition of data, analysis and interpretation of data; they were, moreover, involved in writing and editing the manuscript. In details, LY reconstructed a 3D model of the maxilla and dental row by reverse CBCT images, and performed simulation analysis of the FEM; XFM realized the measurement experiment of tooth movement by a waxy model in a waterbath environment; YFG analyzed and interpreted all data: together, they were the major contributors in preparing and writing the manuscript. JGJ, YDZ and YL contributed to the conception of the study and revised the manuscript. All authors have read and approved the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The study design was approved by the Medical Ethics Committee of the 2nd Affiliated Hospital of Harbin Medical University. And all methods and procedures in this paper were conducted in accordance with relevant guidelines and regulations, as well as the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Furthermore, this paper does not contain any studies performed by any of the authors on animals. It is noteworthy that according to the ethical standards and the approval of the Medical Ethics Committee of the 2nd Affiliated Hospital of Harbin Medical University, all participants provided verbal informed consent.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Jiang, J., Yao, L., Zhang, Y. et al. Establishment, FEM analysis and experimental validation of tooth movement prediction model of orthodontic archwire Tloop. BMC Oral Health 22, 406 (2022). https://doi.org/10.1186/s12903022024309
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12903022024309
Keywords
 Tooth movement prediction model
 Tloop
 Finite elements method (FEM) analysis
 Orthodontic force measurement