Abstarct: In this research, eXtended Finite Element Methode (XFEM) is used in computation of Stress Intensity Factors (SIFs) for an isotropic 2D finite domain with stationary edge crack subjected to thermal shock. Fully coupled displacement, temperature and moisture fields are considered in governing equations. Fourier’s and Fick’s law used for heat and moisture flux respectively. Also Soret flux is used for the caused moisture flux by heat. Interaction integral method is used for hygrothermal loadings to obtain stress intensity factors. Galerkin’s Method is used to discrete the governing equations and then the Newmark time integration scheme is used to solve them. In several numerical examples, the effect of hygrothermal shock on stress intensity factors as well as temperature and moisture distribution is studied. baxsed on the results, from the beginning of hygrothermal load untill the stress wave reaches to crack tip, the stress intensity factors increases and as the stress wave passes through the crack tip the stress intensity factor changes are small and then when the stress wave passes the crack tip, the stress intensity factors decrease with fast speed. Also the maximum of stress intensity factors depends on crack angle and when the angle of crack increases the stress intensity for mode I decreases and mode II increases. The Soret number also have effect on stress intensity factors beside the effect on moisture diffusion.