Abstract: Objective：To analyze the stress distribution in the periodontal ligament (PDL) under different loading conditions at the stage of space closure by 3D finite element model of customized lingual appliances. Methods: The 3D finite element model was used in ANSYS 11.0 to analyze the stress distribution in the PDL under the following loading conditions: (1) buccal sliding mechanics (0.75 N,1.00 N,1.50 N), (2) palatal sliding mechanics (0.75 N,1.00 N,1.50 N), (3) palatalbuccal combined sliding mechanics (buccal 1.00 N + palatal 0.50 N, buccal 0.75 N + palatal 0.75 N, buccal 0.50 N+ palatal 1.00 N). The maximum principal stress, minimum principal stress and von Mises stress were evaluated. Results: (1) buccal sliding mechanics（0.75 N,1.00 N,1.50 N）: maximum principal stress: at the initial of loading, maximum principal stress, which was the compressed stress, distributed in labial PDL of cervix of lateral incisor, and palatal distal PDL of cervix of canine. With increasing loa-ding, the magnitude and range of the stress was increased. Minimum principal stress: at the initial of loading, minimum principal stress which was tonsil stress, distributed in palatal PDL of cervix of lateral incisor and mesial PDL of cervix of canine. With increasing loading, the magnitude and range of minimum principal stress was increased. The area of minimum principal stress appeared in distal and mesial PDL of cervix of central incisor. von Mises stress：it distributed in labial and palatal PDL of cervix of la-teral incisor and distal PDL of cervix of canine initially. With increasing loading, the magnitude and range of stress was increased towards the direction of root. Finally, there was stress concentration area at mesial PDL of cervix of canine. (2) palatal sliding mechanics（0.75 N,1.00 N,1.50 N）: maximum principal stress: at the initial of loading, maximum principal stress which was the compressed stress, distributed in palatal and distal PDL of cervix of canine, and distalbuccal and palatal PDL of cervix of late-ral incisor. With increasing loading, the magnitude and range of the stress was increased. Minimum principal stress: at the initial of loading, minimum principal stress which was tonsil stress, distributed in distalinterproximal PDL of cervix of lateral incisor and mesial-interproximal PDL of cervix of canine. With increasing loading, the magnitude and range of the stress was increased.von Mises stress: von Mises stress distributed in palatal and interproximal PDL of cervix of canine. With increasing loading, the magnitude and range of stress was increased. Finally, von Mises stress distributing area appeared at distalpalatal PDL of cervix of canine. (3) palatalbuccal combined sliding mechanics: maximum principal stress: maximum principal stress still distributed in distal-palatal PDL of cervix of canine. Minimum principal stress: minimum principal stress distributed in palatal PDL of cervix of lateral incisor when buccal force was more than palatal force. As palatal force increased, the stress concentrating area transferred to mesial PDL of cervix of canine.von Mises stress: it was lower and more well-distributed in palatal-buccal combined sliding mechanics than palatal or buccal sliding mechanics. Conclusion: Using buccal sliding mechanics，stress majorly distributed in PDL of lateral incisor and canine, and magnitude and range of stress increased with the increase of loading; Using palatal sliding mechanics，stress majorly distributed in PDL of canine, and magnitude and range of stress increased with the increase of loading; With palatal-buccal combined sliding mechanics, the maximum principal stress distributed in the distal PDL of canine. Minimum principal stress distributed in palatal PDL of cervix of lateral incisor when buccal force was more than palatal force. As palatal force was increasing, the minimum principal stress distributing area shifted to mesial PDL of cervix of canine. When using 1.00 N buccal force and 0.50 N palatal force, the von Mises stress distributed uniformly in PDL and minimal stress appeared.

Key words: Lingual orthodontics, Stress of periodontal ligament, 3 Dimensional finite element analysis