Calculation method for brittle fracture of functional gradient materials

Combined the phase field model with the wavelet dummy node-virtual crack closure technique (WDN-VCCT) used for stress intensity factors (SIFs) calculation. Calculated the node displacement using an improved brittle fracture phase field method, established the relationship between node displacement and node force using WDN-VCCT, and calculated the SIFs at the crack tip. The correctness and accuracy of the proposed method were verified through functional gradient material (FGM) tensile experiment. The influence of crack inclination angle, crack position and gradient index on mechanical response, and SIFs values at the crack tip was discussed. This study provides important computational tools for estimating the service life of FGMs and important references for structural optimization design methods.