Abstract
Objective
A three-dimensional-printed individual titanium plate was applied for maxillary protraction to eliminate side effects and obtain the maximum skeletal effect. This study aimed to explore the stress distribution characteristics of sutures during maxillary protraction using individual titanium plates in various directions and locations.
Methods
A protraction force of 500 g per side was applied at forward and downward angles between 0° and 60° with respect to the Frankfort horizontal plane, after which the titanium plate was moved 2 and 4 mm upward and downward, respectively. Changes in sutures with multiple protraction directions and various miniplate heights were quantified to analyze their impact on the maxillofacial bone.
Results
Protraction angle of 0–30° with respect to the Frankfort horizontal plane exhibited a tendency for counterclockwise rotation in the maxilla. At a 40° protraction angle, translational motion was observed in the maxilla, whereas protraction angles of 50–60° tended to induce clockwise rotation in the maxilla. Enhanced protraction efficiency at the lower edge of the pyriform aperture was associated with increased height of individual titanium plates.
Class III skeletal malocclusion is a common clinical dental abnormality. Maxillary protraction promotes maxillary growth by simulating suture growth around the maxilla. Conventional facemask therapy inevitably leads to dental side effects such as maxillary incisor protrusion, molar distalization, and clockwise rotation of the mandible.1 Several researchers have attempted to use facemasks with temporary anchorage devices to eliminate side effects and obtain the maximum skeletal effect. Sar et al.,2 Cevidanes et al.,3 and Zhou et al.4 have applied titanium plates for maxillary protraction to achieve favorable maxillary protraction outcomes. However, traditional titanium plate surgery is associated with trauma and the risk of tooth root damage. A three-dimensional (3D)-printed individual titanium plate was designed based on the surface morphology of the individual bone tissue, cortex bone thickness, and root distribution location. The titanium plate perfectly matched the cortical bone of the implant area and was positioned away from the root of the tooth.5,6
The 3D finite element method can be employed to analyze the distribution pattern of stress and possible biomechanical changes after loading. However, limited finite element studies have been conducted on maxillary protraction using miniplates. The locations of the miniplates and directions of the orthopedic forces varied in each experiment. Despite these variations, no study has presented sufficient evidence to support this hypothesis. Additionally, the impact of miniplate height on maxillary protraction has not been previously reported. Therefore, individual titanium plates were designed in different areas based on root distribution and bone condition.
To explore the underlying mechanics, a finite element model of maxillary hypoplasia, including the suture, was established. The stress distribution characteristics under different traction directions and miniplate heights were analyzed to provide a theoretical basis for the clinical application of maxillary protraction using an individual titanium plate.
A 3D finite element model was developed using cone-beam computed tomography (CBCT) (slice thickness of 0.250 mm) scan data from an 11-year-old boy diagnosed with maxillary hypoplasia. The hand-wrist radiograph indicated a growth peak, identified by the appearance of the medial seed bone of the thumb and the epiphytic caps on the epiphyses of the middle and radius segments of the third section of the thumb. This study was approved by the Institutional Review Board (IRB) of the Beijing Stomatological Hospital, Capital Medical University (IRB number: CMUSH-IRB-KJ-PJ-2017-04F). The CBCT data were processed using Mimics software (version 15.01; Materialise, Leuven, Belgium), resulting in geometric models of each bone block in the stereolithography (STL) file format through threshold segmentation and 3D reconstruction. Following bone segmentation, adjacent sections were selected slice-by-slice to reconstruct the suture model. Three-dimensional STL files of the bone tissue were then imported into Geomagic software (Geomagic Studio 2013; Raindrop Corp., Wilmington, NC, USA) for additional processing and the generation of a solid model. The shapes of the titanium plate and screws in the anterior tooth area were designed using Creo software (version 3.0; Parametric Technology Corp., Needham, MA, USA) based on the bone surface morphology. Subsequently, the designs were imported into Geomagic to further optimization and assembly in the specified maxillary position. Considering the root distribution, the titanium plate was positioned 4 mm below the pyriform aperture and 10 mm above the alveolar ridge crest. The thickness of each titanium plate was 1 mm, with a hook located on the labial side of the alveolar crest. The miniscrew implants were simplified to a cylindrical shape (1.3 × 6 mm), and the cortical bone thickness at the implant site was approximately 1.5 mm. The Initial Graphics Exchange Specification files were imported into HyperMesh software (version 13.0; Altair Corp., Troy, MI, USA). The entire craniofacial complex model was assumed to be a continuous, homogeneous, isotropic, and linearly elastic material, with no division between cortical and cancellous bones. Referring to previous study, the material properties of the finite element model were defined, including those of cortical bone, suture line, titanium plate, and miniscrew (Table 1).7-9 Tetrahedral mesh division was utilized in the experiment. The finite element model comprised the cranial-maxillary complex and an individual titanium plate, consisting of 435,816 elements and 806,552 nodes (Figure 1). The model included nine defined sutures: frontomaxillary, bilateral zygomaticotemporal, zygomaticomaxillary, zygomaticofrontal, and pterygopalatine.
The finite element model was imported into the Abaqus software (version 6.13; Dassault System, Vélizy-Villacoublay, France). An encastre boundary condition was imposed on the nodes around the foramen magnum and frontal bone to constrain all translational and rotational degrees of freedom in the X-, Y-, and Z-axes of those nodes.
All contact points in this experiment were treated as absolute static contact points. Completely staticity was observed both between the bones and titanium plate, as well as between the titanium plate and miniscrews, indicating an absence of relative sliding.
A protraction force of 500 g per side was applied to the traction hook at forward and downward angles ranging from 0° to 60° relative to the Frankfort horizontal plane (Figure 2). The angle was increased in 10° increments. The von Mises stresses and the 1st principal stresses of the sutures were analyzed. The von Mises stresses provided a comprehensive overview of suture stress, whereas the 1st principal stress was utilized to determine the stress direction of the sutures. A linear graph and stress distribution nephogram were created based on the experimental data.
In the 3D finite element analyses focusing on maxillary protraction with individual titanium plates at varying titanium plate heights, the titanium plate was displaced 2 and 4 mm upward and 2 and 4 mm downward. The position of the traction hook in the mouth remained relatively fixed. Therefore, we redesigned the titanium plate shape and adjusted the length of the hook arm to maintain the height of the force position (Figure 3). A force of 500 g was applied to the titanium plate hooks. In the third step of our experiment, the traction direction was selected to induce translational motion in the maxilla. The von Mises stresses and 1st principal stresses of the sutures were analyzed, and a linear graph and stress distribution nephogram were generated based on the experimental data.
When a 500 g protraction force was applied at angles ranging from 0° to 60° relative to the Frankfort horizontal plane, the von Mises stresses in various sutures exhibited distinct patterns. The stresses in the pterygopalatine, zygomaticomaxillary, and zygomaticotemporal sutures decreased with an increasing protraction angle. For the frontomaxillary and zygomaticofrontal sutures, stresses initially decreased and subsequently increased, reaching a minimum value at a protraction angle of 40°. Among these, von Mises stresses were highest in the pterygopalatine and zygomaticotemporal sutures, whereas they were lowest in the zygomaticomaxillary suture (Figure 4).
The 1st principal stresses in the pterygopalatine, zygomaticotemporal, and zygomaticomaxillary sutures exhibited positive values with an increase in the angle, with the pterygopalatine suture experiencing the highest stress levels. Conversely, the stress in the frontomaxillary suture was negative within an angle of 0° to 30°, decreasing as the angle increased. Between 30° and 40°, the stress transition from negative to positive, subsequently increasing by 40° to 60°. The stress pattern in the zygomaticofrontal suture showed positivity at angles of 0° to 40°, decreasing with an increase in the angle. Between 40° and 50°, the stress changed from positive to negative, becoming negative from 50° to 60°. However, the absolute value of the stress continued to rise with increasing angle. At an angle of 40°, the 1st principal stress of each suture was positive and exhibited a uniform distribution (Figure 5).
The stress in each suture exhibited varying degrees of decrease with the downward movement of the titanium plate when a 500 g protraction force was applied at an angle of 40° relative to the Frankfort horizontal plane (Figure 6).
During the downward movement of the titanium plate, the 1st principal stress in each suture remained positive. Notably, the stresses in the pterygopalatine and zygomaticotemporal sutures decreased, whereas the changes in the stresses in the zygomaticomaxillary, frontomaxillary, and zygomaticofrontal sutures were minimal (Figure 7).
In a comparative analysis of maxillary protraction using a titanium plate in the zygomatic alveolar crest and lateral nasal wall areas, Lee and Baek8 observed greater stress values in the infrazygomatic crest. However, the titanium plate in the nasal sidewall area can reduce the counterclockwise rotation of the maxilla. Kim et al.10 revealed that maxillary protraction with the palatal plate resulted in a broader stress distribution, increased forward displacement, and a greater tendency for counterclockwise rotation than the miniplate at the infrazygomatic crest area and conventional tooth-borne appliances. For patients with Class III skeletal malocclusion, counterclockwise rotation of the maxilla is not conducive to improving the facial pattern. Therefore, in this study, titanium plates were placed in the anterior teeth. In traditional titanium plate surgery, a titanium plate is implanted into the lateral edge of the pyriform aperture to prevent damage to the root and permanent germ.11 However, this approach often results in severe trauma. The shape of the individual titanium plate and distribution of the titanium screws can be designed based on the CBCT data. Under the premise of avoiding damage to the root, an individual titanium plate can be implanted at the lower edge of the pyriform aperture to decrease surgical trauma and postoperative reactions. In this study, the range of vertical movement of the titanium plate was selected from the lower edge of the pyriform aperture to 6 mm above the alveolar crest, covering the clinical application range.
To determine the reference plane, previous studies, including the occlusal plane10,12 and Frankfort horizontal plane,13,14 were revisited. The occlusal plane, influenced by malocclusion severity and subject to significant variation among individuals, led us to choose the Frankfort horizontal plane as the reference.
Significant differences were observed in protraction force selection among prior studies. Lee and Baek,8 Yan et al.,15 and Zhang et al.16 utilized a force of 500 g per side to simulate maxillary protraction in their finite element models, whereas Tanne et al.,9 Miyasaka et al.,17 and Gautam et al.12 selected a force of 1,000 g. Zou and Lin18 reported that the von Mises stresses, maximum principal stresses, and displacement of the cranial-maxillary complex increased with higher protraction forces. However, orthopedic forces > 500 g per side resulted in an evident counterclockwise rotational tendency and narrowing of the dental arch in the maxilla, as observed in both clinical practice and previous studies. Based on this information, we chose a protraction force of 500 g for this study.
In the finite element analysis, each node was influenced by principal stresses in three directions, with the first principal stress having the largest absolute value among them. A negative value for the 1st principal stress indicates compressive stress in the sutures, whereas a positive value suggests tensile stress. Because the protraction force guides the growth direction of the maxilla, we focused on the first principal stress, representing the maximum traction force on the sutures used to analyze the growth direction of the maxilla. The frontomaxillary suture experienced compressive stress at a protraction angle of 0° to 30°, whereas the other sutures experienced tensile stress. Situated on the anterosuperior side of these sutures, the frontomaxillary suture showed a tendency for counterclockwise rotation in the maxilla. At a protraction angle of 40°, all sutures experienced tensile stress, indicating translational motion of the maxilla. Between a protraction angle of 50° and 60°, the zygomaticofrontal suture received compressive stress, whereas the other sutures received tensile stress. Positioned on the posterosuperior side of sutures, the zygomatic frontal suture indicated a tendency for clockwise rotation in the maxilla. The pterygopalatine, zygomaticotemporal, and zygomaticomaxillary sutures, distributed nearly longitudinally, generally experienced tensile stress under the horizontal component of maxillary protraction. For Class III patients with a deep overbite tendency, a protraction angle < 40° with respect to the Frankfort horizontal plane is recommended. Patients with a normal overbite may benefit from a protraction angle of approximately 40°, allowing translational movement of the maxilla along the translational direction. For Class III patients with an open bite tendency, an angle > 40° is suggested.
According to previous studies, Tanne et al.9 and Miyasaka et al.17 supported the idea that the optimal choice for maxillary protraction involves a forward and downward force angled at 30° with respect to the occlusal plane at the canines. The maxilla undergoes forward translation, and stress distribution is uniform. In this study, a titanium plate, chosen as the protraction device, was positioned at the lower edge of the pyriform aperture, causing its force point to move upward compared with that in the tooth anchorage, potentially influencing outcomes. Furthermore, prior studies did not construct a suture model and overlooked the significance of sutures in maxillary protraction. Yan et al.15 suggested translational motion of the cranial-maxillary complex when the force direction is approximately 20° relative to the occlusal plane during bone anchorage protraction, with bone anchorage at the infrazygomatic crest, which is not consistent with our findings. Qu et al.14 conducted maxillary protraction using an implant anchorage positioned between the lateral incisors and canines. They observed that the maxilla moved forward and rotated upward when the protraction angle was < 30°. At an angle of 40°, the suture in the nasion did not experience tensile or compressive stress, indicating translational motion of the maxilla. The maxilla moves forward and rotated downward when the angle exceeded 50°. The protraction loading position in this study was similar to that in Qu et al.,14 with both studies using the Frankfort horizontal plane as the reference plane. Therefore, our study results closely resemble those reported by Qu et al.14 This divergence might be attributed to the larger stress area of maxillary protraction with the titanium plate than the implant screws. Previous studies exhibited variations in maxillary protraction outcomes due to differences in finite element models, anchorage position, reference plane selections, and analysis indexes. In clinical applications, the direction of maxillary protraction should be determined based on the anchorage location and the vertical facial type of the patient.
In a 3D solid structure, the stress distribution is intricate and cannot be captured by a single directional stress value. The von Mises stress is often employed for heterogeneous objects with complex mechanical properties, including three principal stresses in directions equivalent to the same axis. The von Mises stresses of the pterygopalatine, zygomaticotemporal, and zygomaticomaxillary sutures decreased with an increasing protraction angle. The highest stress was observed when the force direction was approximately 0° with respect to the Frankfort horizontal plane. These results are consistent with those obtained by Yan et al.15 This consistency is attributed to the rough longitudinal distribution of the sutures. As the protraction angle increases, the component of the forward force gradually decreases, indicating that the protraction value should be appropriately increased with an increase in the protraction angle. The von Mises stresses on the frontomaxillary and zygomaticofrontal sutures initially decrease and then increase in accordance with the change in the 1st principal stress. In other words, the stress on the zygomaticomaxillary suture changes from negative to positive, whereas that on the zygomaticofrontal suture shows the opposite result.
Stress magnitude is influenced by several factors. First, stress magnitude is associated with the distribution of the sutures. The results showed higher stresses on the pterygopalatine and zygomaticotemporal sutures due to their vertical distribution and exposure to a greater forward component force. Secondly, stress is defined as the external force per unit area, correlating with the suture area. Kambara19 observed increased activity in the reconstruction of the zygomaticotemporal suture during maxillary protraction. The narrow shape of the suture led to stress concentration. Furthermore, stress is influenced by the distance from the protraction point and buffer during the force-transfer process. Therefore, stress magnitudes are the outcome of these combined factors.
Individual titanium plates were custom-designed based on bone thickness and root distribution, offering a range of options for implant height. The position of the traction hook in the mouth remained fixed, whereas the titanium plate was adjusted upward by 2 and 4 mm and downward by 2 and 4 mm. This study included various miniplate heights to encompass a broad spectrum of clinical scenarios. The analysis of stress distribution trends aimed to provide guidance for determining the optimal position of the titanium plate.
When the titanium plate moved downward, the von Mises stress of each suture decreased, albeit at varying amplitudes. This suggests an improvement in protraction efficiency with an increase in the titanium plate position. In clinical practice, when designing a titanium plate in a lower position, the force applied should be appropriately increased. The decrease in stress may be attributed to the increased distance between the application point and the sutures. Although this study designed the titanium plate with different heights, the position of the hooks remained constant. Therefore, factors such as hook length and angle variation could potentially contribute to these observed outcomes.
Throughout the movement of the titanium plate, the 1st principal stresses consistently exhibited positive values. This indicates that all sutures experienced tensile stress, resulting in translational motion of the maxilla. The findings suggest that the height of the titanium plate had an insignificant effect on maxillary rotation within the designed range. Therefore, in clinical applications, the choice of titanium plate position can be made based on considerations such as root distribution and bone mass.
This study had several limitations. The accuracy of the model was constrained by the thickness of the image layer. Additionally, the assumption of homogeneous, continuous, and isotropic bone tissue in the analysis deviates from the actual complex structure. Although finite element analysis offers insights for clinical applications, the influence of soft tissues such as muscles and ligaments on growth could lead to variations between experimental outcomes and real-world effects. The shape of an individual titanium plate can be customized based on various requirements, offering more choices for positioning. Further research is necessary to enhance guidance for the design of titanium plates.
Various protraction directions are suitable for patients with different vertical bone surfaces. A protraction angle of 0° to 30° relative to the Frankfort horizontal plane is recommended for Class III patients with a deep overbite tendency. For patients with a normal overbite, a protraction angle of 40° is deemed suitable. Class III patients with an open bite tendency may benefit from a protraction angle of 50° to 60°.
When the titanium plate is placed at a lower position, the protraction force can be increased accordingly. Within this range, the impact of the height of the titanium plate on maxillary rotation was found to be less significant than that of the degree of protraction.
Notes
AUTHOR CONTRIBUTIONS
Conceptualization: YB. Data curation: QC. Formal analysis: FW. Funding acquisition: YB. Investigation: FW. Methodology: FW, QC. Project administration: SL. Resources: QC, FW. Software: QC. Supervision: YB. Validation: QC. Visualization: SL. Writing–original draft: FW, QC. Writing–review & editing: SL, YB.
REFERENCES
1. Dietrich UC. 1970; Morphological variability of skeletal Class 3 relationships as revealed by cephalometric analysis. Rep Congr Eur Orthod Soc. 46:131–43. https://pubmed.ncbi.nlm.nih.gov/5287485/.
2. Sar C, Arman-Özçırpıcı A, Uçkan S, Yazıcı AC. 2011; Comparative evaluation of maxillary protraction with or without skeletal anchorage. Am J Orthod Dentofacial Orthop. 139:636–49. https://doi.org/10.1016/j.ajodo.2009.06.039. DOI: 10.1016/j.ajodo.2009.06.039.
3. Cevidanes L, Baccetti T, Franchi L, McNamara JA Jr, De Clerck H. 2010; Comparison of two protocols for maxillary protraction: bone anchors versus face mask with rapid maxillary expansion. Angle Orthod. 80:799–806. https://doi.org/10.2319/111709-651.1. DOI: 10.2319/111709-651.1. PMCID: PMC2930261.
4. Zhou YH, Ding P, Lin Y, Qiu L. 2007; Preliminary study on miniplate implant anchorage for maxillary protraction. Chin J Orthod. 14:102–5. Chinese. https://caod.oriprobe.com/articles/52754/Preliminary_study_on_miniplate_implant_anchorage_f.htm.
5. Liang S, Xie X, Wang F, Chang Q, Wang H, Bai Y. 2020; Maxillary protraction using customized mini-plates for anchorage in an adolescent girl with skeletal Class III malocclusion. Korean J Orthod. 50:346–55. https://doi.org/10.4041/kjod.2020.50.5.346. DOI: 10.4041/kjod.2020.50.5.346. PMID: 32938827. PMCID: PMC7500574.
6. Liang SR, Wang F, Zhou DQ, Chang S, Bai YX. 2017; Three-dimensional printed miniplate used for maxillary protraction. Zhonghua Kou Qiang Yi Xue Za Zhi. 52:753–5. Chinese. https://pubmed.ncbi.nlm.nih.gov/29275570/. DOI: 10.3760/cma.j.issn.1002-0098.2017.12.008. PMID: 29275570.
7. Yu HS, Baik HS, Sung SJ, Kim KD, Cho YS. 2007; Three-dimensional finite-element analysis of maxillary protraction with and without rapid palatal expansion. Eur J Orthod. 29:118–25. https://doi.org/10.1093/ejo/cjl057. DOI: 10.1093/ejo/cjl057. PMID: 17218719.
8. Lee NK, Baek SH. 2012; Stress and displacement between maxillary protraction with miniplates placed at the infrazygomatic crest and the lateral nasal wall: a 3-dimensional finite element analysis. Am J Orthod Dentofacial Orthop. 141:345–51. https://doi.org/10.1016/j.ajodo.2011.07.021. DOI: 10.1016/j.ajodo.2011.07.021. PMID: 22381495.
9. Tanne K, Hiraga J, Sakuda M. 1989; Effects of directions of maxillary protraction forces on biomechanical changes in craniofacial complex. Eur J Orthod. 11:382–91. https://doi.org/10.1093/oxfordjournals.ejo.a036010. DOI: 10.1093/oxfordjournals.ejo.a036010.
10. Kim KY, Bayome M, Park JH, Kim KB, Mo SS, Kook YA. 2015; Displacement and stress distribution of the maxillofacial complex during maxillary protraction with buccal versus palatal plates: finite element analysis. Eur J Orthod. 37:275–83. https://doi.org/10.1093/ejo/cju039. DOI: 10.1093/ejo/cju039.
11. Kircelli BH, Pektas ZO. 2008; Midfacial protraction with skeletally anchored face mask therapy: a novel approach and preliminary results. Am J Orthod Dentofacial Orthop. 133:440–9. https://doi.org/10.1016/j.ajodo.2007.06.011. DOI: 10.1016/j.ajodo.2007.06.011. PMID: 18331946.
12. Gautam P, Valiathan A, Adhikari R. 2009; Maxillary protraction with and without maxillary expansion: a finite element analysis of sutural stresses. Am J Orthod Dentofacial Orthop. 136:361–6. https://doi.org/10.1016/j.ajodo.2008.02.021. DOI: 10.1016/j.ajodo.2008.02.021. PMID: 19732670.
13. Lin QI, Jing-Hua T, Zhen-Jin Z, Yang Z. 2013; Study of stress distribution of craniofacial sutures following maxillary protraction with palatal expansion- a three-dimensional FEM study. Chin J Med Phys. 30:3933–5. http://zgyxwlxzz.paperopen.com/en/oa/DArticle.aspx?type=view&id=201301024.
14. Qu Y, Liu J, Yang S, Huang Y. 2016; 3-D finite element study on maxillary body protraction using implant anchorage. J Pract Stomatol. 32:58–62. https://caod.oriprobe.com/articles/48123795/3_D_finite_element_study_on_maxillary_body_protrac.htm.
15. Yan X, He W, Lin T, Liu J, Bai X, Yan G, et al. 2013; Three-dimensional finite element analysis of the craniomaxillary complex during maxillary protraction with bone anchorage vs conventional dental anchorage. Am J Orthod Dentofacial Orthop. 143:197–205. https://doi.org/10.1016/j.ajodo.2012.09.019. DOI: 10.1016/j.ajodo.2012.09.019. PMID: 23374926.
16. Zhang GH, Cai Z, Lu Q. 2000; Three-dimensional finite element method study on the maxillary protraction: the direction of the protraction. J Med Biomech. 15:208–11. http://www.mechanobiology.cn/yyswlxen/article/abstract/200004004?st=article_issue.
17. Miyasaka J, Tanne K, Tsutsumi S, Sakuda M. 1986; Finite element analysis of the biomechanical effects of orthopedic forces on the craniofacial skeleton. Construction of a 3-dimensional finite element model of the craniofacial skeleton. Osaka Daigaku Shigaku Zasshi. 31:393–402. Japanese. https://pubmed.ncbi.nlm.nih.gov/3469386/.
18. Zou M, Lin WJ. 2011; A three-dimensional finite element study on the influence of different protraction force on cranio-maxillary complex. Shanghai Kou Qiang Yi Xue. 20:88–92. Chinese. https://pubmed.ncbi.nlm.nih.gov/21451906/.
19. Kambara T. 1977; Dentofacial changes produced by extraoral forward force in the Macaca irus. Am J Orthod. 71:249–77. https://doi.org/10.1016/0002-9416(77)90187-7. DOI: 10.1016/0002-9416(77)90187-7. PMID: 402814.