Structural cohesive element for the modelling of delamination in composite laminates without the cohesive zone limit

Xiaopeng Ai,Boyang Chen,Christos Kassapoglou
2024-05-27
Abstract:Delamination is a critical mode of failure that occurs between plies in a composite laminate. The cohesive element, developed based on the cohesive zone model, is widely used for modeling delamination. However, standard cohesive elements suffer from a well-known limit on the mesh density-the element size must be much smaller than the cohesive zone size. This work develops a new set of elements for modelling composite plies and their interfaces in 3D. A triangular Kirchhoff-Love shell element is developed for orthotropic materials to model the plies. A structural cohesive element, conforming to the shell elements of the plies, is developed to model the interface delamination. The proposed method is verified and validated on the classical benchmark problems of Mode I, Mode II, and mixed-mode delamination of unidirectional laminates, as well as on the single-leg bending problem of a multi-directional laminate. All the results show that the element size in the proposed models can be ten times larger than that in the standard cohesive element models, with more than 90% reduction in CPU time, while retaining prediction accuracy. This would then allow more effective and efficient modeling of delamination in composites without worrying about the cohesive zone limit on the mesh density.
Computational Engineering, Finance, and Science,Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is that in the delamination modeling of composite laminates, the standard cohesive elements (CEs) are restricted by the size of the cohesive zone. Specifically, standard CEs require that the mesh size must be much smaller than the Cohesive Zone Length (CZL) in order to accurately predict delamination. However, this leads to problems of high computational cost and difficult meshing. The paper proposes a new structural cohesive element, aiming to overcome this limitation and allowing for efficient delamination modeling with larger element sizes while maintaining prediction accuracy. ### Main contributions 1. **New element development**: Developed a triangular Kirchhoff - Love shell element and a matching structural cohesive element for modeling composite laminates and their interfaces. 2. **Performance verification**: Verified the effectiveness and accuracy of the new method through classic Mode I, Mode II and mixed - mode delamination tests as well as the single - leg bending problem of multi - directional laminates. 3. **Efficiency improvement**: The results show that the element size of the new model can be ten times larger than that of the standard cohesive element model, the CPU time is reduced by more than 90%, while the prediction accuracy remains unchanged. ### Problems solved - **Cohesive zone size limitation**: Standard CEs require very fine meshes in delamination modeling, which leads to high computational costs and complex meshing. - **Computational efficiency**: The new method significantly reduces the computational time and improves the modeling efficiency. - **Applicability**: The new method is applicable to composite delamination modeling in three - dimensional space without artificially reducing the material strength or presetting cracks. ### Key technologies - **Triangular Kirchhoff - Love shell element**: Used for modeling composite laminates, with each node having three degrees of freedom (one out - of - plane displacement and two rotations). - **Structural cohesive element**: Matches the shell element and is used for modeling interface delamination. It places more integration points through an adaptive integration scheme to capture the stress gradient within the cohesive zone. ### Conclusion The new method proposed in the paper successfully overcomes the cohesive zone size limitation of standard cohesive elements in delamination modeling, allowing for efficient delamination modeling with larger element sizes while maintaining the prediction accuracy. This provides strong support for the reliable design of composite structures.