In this paper, a method for extracting stress intensity factors (SIFs) in orthotropic thermoelasticity fracture by the extended finite element method (XFEM) and interaction integral method is present. The proposed method is utilized in linear elastic crack problems. The numerical results of the SIFs are presented and compared with those obtained using boundary element method (BEM). The good accordance among these two methods proves the applicability of the proposed approach and conforms its capability of efficiently extracting thermoelasticity fracture parameters in orthotropic material.
