Implementation of the extended finite element method for coupled dynamic thermoelastic fracture of a functionally graded cracked layer

The effect of thermal shock on a cracked functionally graded material (FGM) layer is considered using the extended finite element method. Classical coupled thermoelastic equations are used in the calculations. The coupled dynamical system of equations obtained from the extended finite element discretization is solved by the Newmark method in the time domain. Micromechanical models for conventional composites are used to estimate the material properties of functionally graded layer. The interaction integral is then employed to calculate the dynamic thermal stress intensity factors (SIFs) at each time step. The effects of initial crack angle and volume fraction profiles of FGMs on SIFs are studied. Also crack propagation phenomenon is investigated in this paper. We have used MATLAB software to do the different stages of simulation from mesh generation to numerical computation of SIFs. Some numerical examples are implemented to investigate the validity and accuracy of attained results.