INTRODUCTION
Lingual orthodontics is widely recognized and accepted by adult orthodontic patients because of its relatively beautiful features. However, compared with labial orthodontics, excessive lingual inclination of the maxillary anterior teeth and vertical and horizontal “arched effect” of dental arch easily appeared in the space closure stage during lingual orthodontics, because of large variability of anterior tooth lingual surface , small effect of bracket spacing on the rigidity and friction of the arch wire, and small bracket span between the orthodontic force point and the tooth impedance center (Lombardo et al., 2012; Long et al., 2013; Antoszewska-Smith et al., 2017).This has affected the orthodontics effect, and limited the clinical application of lingual orthodontics. In this study we investigated the biomechanics effects of three different wire types applied the retraction force and the labial intrusive force in newly designed Duet-slot lingual bracket featuring with occlusal and lingual slots by using finite element analysis (FEA). We examined the effect of initial tooth movement and stress distribution of periodontal ligament (PDL) and tooth root for providing the lingual technique with valuable information to achieve the best orthodontic result.
MATERIAL AND METHOD
Establishment of 3D finite element models of microimplants-lingual straight wire appliance during retracting maxillary anterior teeth. One normal cranial bone with normal tooth, complete dentition and symmetrical arch were selected and scanned with cone-beam CT (KAVO, German). Data with DICOM file were output. Two-dimensional CT scan results were introduced into Mimics 17.0 software. 3D information of maxilla and dentition was extracted by adjusting threshold so that the images of the maxilla and dentition were clearly visible. The maxillary models containing the entire dentition were obtained by the Boolean operation to simulate the perfectly aligned teeth after orthodontic treatment. Above 3D models were modified with software Unigraphics NX (version 8.5, Siemens Inc, Germany) to establish models of lingual brackets and arch wires. The geometric models of the Duet-slot lingual Bracket were based on standard sample of DLB from the School of Stomatology, Kyungpook National University, Daugu, South Korea (Fig. 1). In accordance with clinical conditions and properties of biomaterials, the above mentioned models were classified into cortical bone, cancellous bone, tooth, periodontal ligament, arch wire, traction hook, coil spring, micro-implant and bracket. The periodontal ligament (PDL) was simulated as a 0.2 mm thick layer around the roots, and the cortical bone was 0.5 mm thick (Kojima & Fukui, 2006; Kojima et al., 2007, 2012). Isotropic and linearly elastic behavior was assumed for all materials, and the mechanical properties were taken from previous studies (Pelsue et al., 2009; Kojima & Fukui, 2012; Xia et al., 2013; Facchini et al., 2017) (Table I).
Finite element analysis software ANSYSnsys Workbench 15 (ANSYS, USA) was used to generate tetrahedral 10-node element. It is assumed that the maxillary surface is a fixed end. The teeth-brackets, periodontal ligament-teeth, alveolar bone-periodontal ligament, and cortical bone-cancellous bone were connected with bonded; bracket-archwire was connected with No separation. Therefore, free axial rotation movement of the arch-wire in the brackets was allowed, while friction between the archwire and brackets along the axial direction was ignored. The established 3D finite element models of micro-implantslingual straight wire appliance are exhibited in Figure 2.
In the system studied, the x-axis was the midsagittal line of the dental arch on the occlusal view, The direction from right to the left of arch was defined as +X. The direction of Y was perpendicular to the X on the plane, and the direction from Labial surface to the lingual surface of the central incisors was defined as +Y. The direction of Z was perpendicular to the XY plane. The occlusal direction towards gingival direction referred to +Z (Fig. 2).
Loading conditions. The position of the microimplants(MIA) located the apical areas of the first and second molars for retraction the anterior teeth was assumed to be 8 mm from the alveolar crest (Fig. 3). The micro-implants situated the labial region between the two central incisors were assumed to apply vertical traction for clinic simulation the intrusion anterior from the threads fixed in the labial crown using light-cured resins (Fig. 4). The retraction hook(RH) was placed bilaterally between the lateral incisor and the canine at a height of 6 mm. Retraction force of 1.5N from the RH to MIA traction was applied, and the additional intrusion force (50 g) from the two incisors was combined to simulate the effect on labial crown torque.
The working wires for en-masse anterior retraction with MIA by sliding mechanics were divided into three groups: 0.016 X 0.022- in stainless steel(SS) rectangular wire (group A); 0.016 -in SS round wire (group B); 0.016 X 0.022in SS + 0.016-in SS (group C) . the initial displacement of the anterior dentition and the stress of each their root and periodontal ligament were analyzed.
RESULTADOS
Displacement of anterior dentition. In accordance with 3D finite element analysis models and three groups of loading conditions, initial displacement of anterior dentition was calculated and results were displayed in Figure 5. The initial displacement was the largest in the incisal edge of lateral incisors compared with the canine central incisors. The displacement gradually decreased from the incisal edge to the root apex. The displacement of anterior teeth of group C was more uniform than that of groups A and B (Table II).
Horizontal displacement of anterior dentition. As shown in Figure, anterior dentition of the three groups presented obvious expansion of dental arch in canine area because the crown and root of canine were both moved labially. Moreover, the amount of displacement on the cuspid of canine in the group B was approximately twice than that in the groups A and C. As shown in Table II, central incisors presented the tipping movement with crown moved medially and apex moved distally in the horizontal direction. The horizontal displacement of the lateral incisors was similar to that of the central incisors with more pronounced mesial tipping.
Table II Initial displacement of crown and root on X, Y, Z axis of the anterior teeth in three groups (mm).

Sagittal displacement of anterior dentition. The maximum of the sagittal displacement was at the incisal edge of the lateral incisors in all experimental groups. In particular, the lateral incisors in the group B presented obvious tipping movement with crown towards lingual direction. In the groups A and C, the central incisors moved close to the lingual bodily movement (Fig. 7). Moreover, the difference of crown and root displacement was smaller in the group C than in the group A. The canine teeth presented tipping movement of crown to distal direction and apex to mesial direction in the groups A and C, however, in group B the canine moved mesially both crown and root (Table III).These results indicate that canines have the tendency of mesial eversion. A pronounced horizontal “arched effect” was seen in the group B.

Fig. 6 Horizontal displacement of anterior teeth. U1 Right center incisor U2 Right lateral incisor U3 Right canine).
Table III Initial displacement of crown and root on Y axis of the anterior teeth in three groups (mm).

Vertical displacement of anterior dentition. The dental cervix and root apex of the lateral incisors demonstrated remarkable intrusion in all groups A, B and C. The incisal edge of lateral incisor was elongated. The maximum vertical displacement was much smaller in the groups A and C than that in the group B (Figs. 8 and 9). The central incisors in the three groups were extruded . Also , the canines in the groups A and C were extruded. However, the canine presented crown intrusion, and apex extrusion in the group B, Which indicated the canine obviously tent to crown lingual tipping movement.
Anterior tooth root stress. Figure 10 showed the distribution of S1 and S3 from the cervical margin to the apex of the maxillary anterior tooth root, The maximum (S1) stress was produced and mostly concentrated near the cervix on the lateral incisors and the canine. This suggested that dental cervix of maxillary anterior teeth on lingual side was a critical stress-strain region, which was consistent with the tendency of secondary tooth movement after bone reconstruction. In the group A, the stress of 1.1297 MPa was produced in the cervix of lateral incisors. In the group B, the maximum stress of 1.8762 MPa was produced in the lingual cervix of lateral incisors and canine teeth. In the group C, the stress of 1.0986 MPa was produced in the cervix of lateral incisors. It is thus evident that the root stress of anterior teeth was smallest in models in the group C; moreover, the stress distribution was uniform and reasonable, and the external root resorption was difficult to occur.
Stress distribution in the periodontal ligament of anterior teeth.Stress distribution of periodontal ligament was similar to the stress of tooth root in the three groups, and was significantly concentrated at cervix of lateral incisors and canines. Maximum stress of periodontal ligament was distributed at the cervix of lateral incisors on the lingual side in the group A (0.14059 MPa), in the group B which located the mesial cervix of lateral incisors (S1 0.25341 MPa), however, the cervix of lateral incisors was concentrated on the lingual side in the group C (0.13733 MPa). The stress of periodontal ligament was larger in the group B than in the groups A and C. The stress of periodontal ligament in the group C was smallest compared with the groups A and B (Fig. 11).
DISCUSSION
Lingual orthodontic system invented by Dr. Kurz and Dr. Fujita in 1970s, because of its highest esthetics, has developed rapidly in recent years.Compared with the labial orthodontocs, lingual orthodontics is difficult to control the maxillary incisors torque during the movement of retraction. The expression of the torque of the lingual brackets is usually achieved by the use of a personalized lingual resin base.However, full-size archwire can cause discomfort and difficulty in arch wire insertion.Therefore, after the dental arch is completely aligned and leveled, the 0.016 inch × 0.022 inch SS is commonly used to close the space . But, the torque provided by the bracket itself is not sufficient to provide the required torque in some clinical applications. In the case of tooth extraction, the “arched effect” easily appears in lingual orthodontics, Thus the micro-implant anchorage is inserted in palatal alveolar bone between the first molar and the second molar for en-masse space closure in lingual orthodontics. In some patients, however, retroclined anterior teeth with only retraction force can cause a torque lost of anterior teeth. Dr Liang analyzed initial displacement of the tooth and the stress and strain distribution of the periodontal ligament in 3D finite element models, and found that the lingual retraction force was more likely to cause anterior teeth crown lingual tipping movement (Liang et al., 2009). It is essential to achieve the better torque control towards root lingual movement. Thus, in this study we designed a labial vertical intrusive force in addition to the horizontal retraction force for space closure by palatal microimplant.
Duet Lingual Orthodontic Bracket (DLB) featuring two slots used in this study has been more applied in and outside China.The occlusal slot with veitical direction is convenient for insertion and removal of archwire, The lingual slot with horizontal insertion direction is closer to the CR compared to the occlusal slot, and facilitates the control of the torque and the mesiodistal tipping.So, two types arch wires can be simultaneously inserted in the two slot.This study focused on analysis of the biomechanical characteristics about the stress and movement of anterior teeth in 3D finite element models after inserting 0.016 inch × 0.022 inch stainless steel square wire (lingual slot), 0.016 inch stainless steel round wire (lingual slot), 0.016 inch × 0.022 inch stainless steel square wire (lingual slot) + anterior teeth 0.016 inch stainless steel round wire (occlusal slot).
With the rapid development of digital technology, the finite element analysis has become an effective research method of oral biomechanics (Geron et al., 2004; Cattaneo et al., 2015; Wang et al., 2017). 3D finite element method is an important method for studying the movement of orthodontic teeth (Murakami & Wakabayashi, 2014).The reaction of teeth to sliding mechanics can be analyzed by using finite element analysis (FEA) (Namburi et al., 2017). With 3-dimensional (3D) computer models, various conditions can be simulated by varying the simulation parameters. The initial reactions of the teeth, periodontal ligament (PDL), and alveolar bone can be evaluated qualitatively and quantitatively.
The torque control of the anterior teeth has always been the focus of lingual orthodontic design using microimplants during retracting maxillary anterior teeth. In this study , the results showed that it presented noticeable horizontal “arched effect”, and the initial displacement of anterior teeth was tipping movement (with lingual crown and labial apex movement) in group B with only 0.016 -in round wire, whereas the initial displacement of lateral incisors was almost translationin the groups A and C, also in the group C the anterior teeth appeared better bodily movement than in group A. The possible reasons for the biomechanical difference was the 0.016 inch × 0.022 inch SW can provide the better torque control of the anterior teeth than 0.016 inch RW. In addition, when inserting two wire in the occlusal and lingual slot, round wires in the occlusal slot can maintain the arch form and control the rotation of anterior teeth.
Although we designed the intrusion force vector during retraction, the central incisors still showed light extrusion in all groups. From mechanics, the intrusive force generates two orthodontic effects: intrusion of the maxillary incisor and an effect equal to lingual root torque. But, in orthodontic practice, a heavy intrusive force could cause root absorption more easily. Therefore, in Lingual orthodontic treatment, it is more important to control torque with light lingual retraction force and appropriately intrusive force to counteract the tendency of lingual crown tipping (Melsen, 2001). In the horizontal direction, initial displacement of mesioclination could be seen in the central incisors and lateral incisors in the three groups. However, Maximum initial displacement in the group B was twice as much as that in thegroups A and C. The stress was uneven in the group B. Canine teeth presented labial displacement and widened dental arch in the three groups, but tipping displacement was not obvious in the group C than in the group A.
Lee reported that the maximum stress that the periodontal ligament can sustain is 0.026 Mpa; if the strength exceeds the limit, the periodontal ligament will suffer permanent damage (Lee, 1965). In this study, the maximum stress of the periodontal ligament was 0.14059 MPa and 0.13733 MPa in the groups A and C, respectively, which was far less than 0.026 Mpa, and could not cause permanent damage to the periodontal ligament. Nevertheless, the maximum stress of the periodontal ligament was 0.25341 MPa in the group B, which was close to the critical value and had a risk of injury to the periodontal ligament. Moreover, the stress of tooth root was relatively large in the group B, suggesting the possibility of external resorption of tooth root. So, the group C with DW can better control the movement of anterior teeth in the space closure stage during lingual orthodontics.
The 3D-FEA analysis reveal that the group C namely double wire is highly efficient in torque control compared with the conventional single round wire or rectangular wire in lingual orthodontics. These results suggest that en-masse bodily movement of anterior teeth seems to be difficult though increase the vertical intrusion force by using MIA. However, these results are right only for the initial movement, which is produced by elastic deformation of the PDL. This is a limitation of this study. Long-term orthodontic movement might not be the same as the initial movement.