4. Conclusion
A new method is proposed for modeling high cycle fatigue delamination in composite materials. The presented approach is based on a cohesive zone model and uses the level set method for tracking the crack front evolution in 3D problems. The integration of the traction-separation curve at integration points is used to calculate the energy release rate accurately. This method is superior to VCCT since it can be used for bi-material interfaces and no remeshing is required for arbitrary crack growth. Smoothing is proposed to reduce the oscillations of the calculated energy release rates at integration points. Since the velocity at the nodes is required by the level set method, the Newton–Cotes integration scheme is used. Unlike previous fatigue models based on the cohesive zone method, the presented approach does not introduce new material parameters or require calculation of an effective length. The benchmarks show the ability of the method to simulate fatigue-driven delamination and track the front evolution in 3D.