Direct estimation of generalized stress intensity factors using a three-scale concurrent multigrid X-FEM

A concurrent multigrid method is devised for the direct estimation of stress intensity factors (SIF) and higher order coefficients of the elastic crack tip asymptotic field. The proposed method bridges three characteristic length scales that can be present in fracture mechanics: the structure, the crack and the singularity at the crack tip. For each of them, a relevant model is proposed. First, a truncated analytical reduced order model based on the Williams’ expansion is used to describe the singularity at the tip. Then, it is coupled with a standard X-FEM model which is known to be suitable for the scale of the crack. A multigrid solver finally bridges the scale of the crack to that of the structure for which a standard FE model is often accurate enough. Dedicated coupling algorithms are presented and the effects of their parameters are discussed. The efficiency and accuracy of this new approach is exemplifyed using three benchmarks.