A stabilized finite element method for delamination analysis of composites using cohesive elements

Gourab Ghosh,Ravindra Duddu,Chandrasekhar Annavarapu

DOI: https://doi.org/10.48550/arXiv.2008.09015

2020-08-20

Abstract:We demonstrate the ability of a stabilized finite element method, inspired by the weighted Nitsche approach, to alleviate spurious traction oscillations at interlaminar interfaces in multi-ply multi-directional composite laminates. In contrast with the standard (penalty-like) method, the stabilized method allows the use of arbitrarily large values of cohesive stiffness and obviates the need for engineering approaches to estimate minimum cohesive stiffness necessary for accurate delamination analysis. This is achieved by defining a weighted interface traction in the stabilized method, which allows a gradual transition from penalty-like method for soft elastic contact to Nitsche-like method for rigid contact. We conducted several simulation studies involving constant strain patch tests and benchmark delamination tests under mode-I, mode-II and mixed-mode loadings. Our results show clear evidence of traction oscillations with the standard method with structured and perturbed finite element meshes, and that the stabilized method alleviates these oscillations, thus illustrating its robustness.

Computational Engineering, Finance, and Science,Numerical Analysis

What problem does this paper attempt to address?

This paper attempts to solve the problems of numerical instability and traction oscillation in the traditional finite element method (FEM) when delamination occurs in composite laminates under multi - directional loads. Specifically, the paper focuses on how to alleviate these problems through a stabilized finite element method based on the weighted Nitsche method, especially when using high cohesive stiffness values. The traditional penalty function method may encounter numerical instability and the need for artificial estimation of the minimum cohesive stiffness when dealing with such problems, while the proposed stabilized method can avoid these shortcomings, allowing the use of arbitrarily large cohesive stiffness values and without the need for engineering estimation of the minimum cohesive stiffness.

By defining the weighted interface traction, the paper realizes a smooth transition from the penalty function method of soft elastic contact to the Nitsche method of rigid contact in the stabilized method. The researchers conducted a number of simulation experiments, including constant - strain patch tests and benchmark delamination tests in mode I, mode II and mixed modes, to verify the effectiveness and robustness of the proposed method. The results show that, compared with the standard method, the stabilized method can significantly reduce traction oscillation, thus improving the accuracy of the analysis.