Development of a cohesive zone model for three-dimensional simulation of joint de-bonding/delamination under mixed-mode I/II fatigue loading

In this work, the cohesive zone model (CZM) developed by some of the authors to simulate the propagation of fatigue defects in two dimensions is extended in order to simulate the propagation of defects in 3D. The paper aims to discuss this issue.

The procedure has been implemented in the finite element (FE) solver (Abaqus) by programming the appropriate software-embedded subroutines. Part of the procedure is devoted to the calculation of the rate of energy release per unit, G, necessary to know the growth of the defect.

The model was tested on different joint geometries, with different load conditions (pure mode I, mode II pure, mixed mode I/II) and the results of the analysis were compared with analytical solutions or virtual crack closure technique (VCCT).

The possibility to simulate the growth of a crack without any re-meshing requirements and the relatively easy possibility to manipulate the constitutive law of the cohesive elements makes the CZM attractive also for the fatigue crack growth simulation. However, differently from VCCT, three-dimensional fatigue de-bonding/delamination with CZM is not yet state-of-art in FE softwares.