Continuum Shape Sensitivity analysis of a mode-I fracture in functionally graded materials

This paper presents a new method for conducting a continuum shape sensitivity analysis of a crack in an isotropic, linear-elastic, functionally graded material. This method involves the material derivative concept from continuum mechanics, domain integral representation of the J-integral and direct differentiation. Unlike virtual crack extension techniques, no mesh perturbation is needed to calculate the sensitivity of stress-intensity factors. Since the governing variational equation is differentiated prior to the process of discretization, the resulting sensitivity equations are independent of approximate numerical techniques, such as the meshless method, finite element method, boundary element method, or others. In addition, since the J-integral is represented by domain integration, only the first-order sensitivity of the displacement field is needed. Several numerical examples are presented to calculate the first-order derivative of the J-integral, using the proposed method. Numerical results obtained using the proposed method are compared with the reference solutions obtained from finite-difference methods for the structural and crack geometries considered in this study.