Numerical study for cohesive zone model in delamination analysis based on higher-order B-spline functions

Abstract

A numerical aspect of the implementation cohesive zone model for delamination analysis is presented in this work by employing interface elements and higher-order B-spline functions. The stress oscillation is addressed, and Newton–Cotes integration scheme is considered as a good candidate to overcome this phenomenon. In order to track the nonlinear equilibrium path, a general arc-length constraint named dissipation-based arc-length method is applied. A numerical example is presented to perform the ability of the formulation in predicting the delamination behavior of studied structures.