Objective: To construct a three-dimensional (3D) finite element model and analyze the biomechanical characteristics during mandibular second molar mesialization using customized lingual appliances and mini-implant. Methods: One adult student volunteer from The First Affiliated Hospital of Zhengzhou University with lower left first molar extraction was selected. After CT scanning, Mimics, Geomagic, ANSYS were employed to develop a 3D finite element model including customized lingual brackets, stainless steel lingual arch wire, buccal buckles, lower dentition, periodontal ligaments and alveolar bone. Four different loading methods (1. the force of 1.470 N loaded at the lingual side only; 2, 3, 4. the forces of 0.490, 0.735 and 0.980 N loaded at both buccal and lingual sides, respectively.) were included. The initial displacements of the lower second molar and stress distribution in the periodontal ligaments were analyzed. Results: More uniform stress distributions in the periodontal membrane of mandibular left second molar were found when the mesial force were loaded at both buccal and lingual sides than the force loaded at lingual side only. More inclination of the second molar was observed when force loaded at both sides than at lingual side only. With the values of the force increased, the mesial inclination potential of the second molar, the initial movement of the second molar and the stress distribution in the periodontal membrane increased. Under the force of 0.735 N at both buccal and lingual sides, mesially inclined initial displacements of the mandibuar second molar were observed and the stress values of von Mises was in the suitable area. Conclusions Less second molar rotation was found when the force loaded at both buccal and lingual sides than loaded at lingual side only. However, force loaded at both sides would increase the measial inclination potential of the second molar.