---

In this paper, the eXtended Finite Element Method is implemented to model the effect of the mechanical and thermal shocks on a cracked 2D orthotropic media. The uncoupled thermoelasticity equations are considered. Isoparametric four-node and eight-node rectangular elements are used to discrete governing equations. The dynamical stress intensity factors are computed by the interaction integral method. The Newmark and the Crank–Nicolson time integration schemes are used to numerical solve the spatial-discretized elastodynamic and thermal equations, respectively. A MATLAB code is developed to carry out all stages of the calculations from mesh generation to post-processing. Several elastic and thermoelastic numerical examples are implemented, to check the accuracy of the results and to investigate the effect of the orthotropic direction on the stress intensity factors.

---

In this paper, the extended Finite Element Method (XFEM) is implemented to compute the Stress Intensity Factors (SIFs) for rectangular media subjected to a hygrothermal loading. In governing hygrothermoelasticity equations, the cross coupled of temperature and moisture fields and temperature-dependent diffusion in some cases are considered. Furthermore, an interaction integral for hygrothermal loading is developed to compute the stress intensity factors. The non uniform mesh of isoparametric eight-nod rectangular element is used in XFEM to decrease the absolute error in SIFs computations. In order to numerical results validation, the SIF of mode I is obtained analytically. The coupled governing equations are firstly decoupled in terms of new variables and then solved by the separation of variable method. According to the results, the moisture concentration gradient has a significant effect on the SIFs so should be considered in the model. Up to reaching temperature to its steady state, the cross coupled of temperature and moisture synchronies their time variation which affects on the time variation of SIF. At early time of thermal shock, the SIF for shorter cracks is not necessarily lesser than the longer ones. Also, the mode I SIF for longer and inclined cracks is smaller. On the other hand, considering the moisture concentration as a temperature function increases the time required to reach the moisture steady state.