The main objective of this paper is to propose a novel method that provides an opportunity to evaluate an orthodontic process at early phase of the treatment. This was accomplished by finding out a correlation between the applied orthodontic force and thermal variations in the tooth structure. To this end, geometry of the human tooth surrounded by the connective soft tissue called the periodontal ligament and the bone was constructed by employing dental CT scan images of a specific case. The periodontal ligament was modeled by finite strain viscoelastic model through a nonlinear stress-strain relation (hyperelasticity) and nonlinear stress-time relation (viscoelasticity). The tooth structure was loaded by a lateral force with 15 different quantities applied to 20 different locations, along the midedge of the tooth crown. The resultant compressive stress in the periodontal ligament was considered as the cause of elevated cell activity that was modeled by a transient heat flux in the thermal analysis. The heat flux value was estimated by conducting an experiment on a pair of rats. The numerical results showed that by applying an orthodontic force to the tooth structure, a significant temperature rise was observed. By measuring the temperature rise, the orthodontic process can be evaluated. KEYWORDS finite element method, mechanical stress, orthodontic loading, thermal analysis, viscohyperelastic behavior