The cohesive zone model (CZM) has found a wide acceptance as a tool for the simulation of delamination in composites and debonding in bonded joints. Recently, fatigue-devoted implementations of CZM have been proposed. In earlier works, the authors have developed a model of the cohesive zone able to correctly simulate the propagation of fatigue defects in two-dimensional geometry. The procedure has been implemented in the finite element solver (Abaqus) by programming the appropriate software-embedded subroutines. Part of the procedure is devoted to the calculation of the strain energy release rate, G, necessary to know the growth of the defect. The model has been then extended to 3D cracks with quasi-straight crack front, where G could be evaluated by the contour-integral on parallel slices along the crack front. The aim of this work is to extend the cohesive zone model calculation of G using contour integral and the crack front evolution to 3D cracks with non-straight crack front.
展开▼