Purpose Orthodontic implant has been used for temporary anchorage for orthodontic treatment.The stress states of loaded orthodontic implant–bone interfaces is related to success of orthodontic implant. Stress distribution were analyzed by finite element method(FEM), and determined the difference in the stress distribution for different thread design and different implant placement angle to identify risk factors for the loosening of orthodontic mini-implants. Methods The Taguchi method was employed to find out the significance of each thread design factor to affect bone stress. Using the computer program(ANSYS) to build up three-dimensional finite element models alveolar bone(cortical bone 2 mm thick) and orthodontic mini-implants(Ti-6Al-4V) corresponding to the Taguchi L9 array. Horizontal force(150 CN) was applied ; stimulation time was 1 second. We calculated the most equivalent stresses(Von Mises) on the bone elements and evaluated stress distribution according to implantation direction (0, 15, 30, 45, 60, 75 degree) and three different kind of orthodontic implant : self-design orthodontic implant corresponding to the Taguchi L9 array result, AbsoanchorR (Dentos, Daegu, Korea), LOMASR (Mondeal Medical System GmbH,Tuttlingen, Germany). Results Equivalent stress on the cortical bone ranged from 1.30 to 4.69 MPa and on the cancellous bone anged from 0.10 to 0.17 MPa depending on the implantation direction. The stress distribution on the surrounding bone of implant placement angle 0 degree was even and small. Bone surrounding stress is small in self-design orthodontic implant corresponding to the Taguchi L9 array result. Conclusions This report suggests information for clinicians to understand the stress distribution on alveolar bone with different placement direction. Implant thread design is also one the factors to facilitate implant stability in clinical treatmentz