A Numerical Study on the Influence of Strain Rate in Finite-Discrete Element Simulation of the Perforation Behaviour of Woven Composites

Predicting the perforation limit of composite laminates is an important design aspect and is a complex task due to the multi-mode failure mechanism and complex material constitutive behaviour required. This requires high-fidelity numerical models for a better understanding of the physics of the perforation event. This work presents a numerical study on the perforation behaviour of a satin-weave S2-glass/epoxy composite subjected to low-velocity impact. A novel strain-rate-dependent finite-discrete element model (FDEM) is presented and validated by comparison with experimental data for impacts at several energies higher and lower than their perforation limit. The strain rate sensitivity was included in the model by developing a novel user-defined material model, which had a rate-dependent bilinear traction separation cohesive behaviour, implemented using a VUSDFLD subroutine in Abaqus/Explicit. The capability of the model in predicting the perforation limit of the composite was investigated by developing rate-sensitive and insensitive models. The results showed that taking the strain rate into account leads to more accurate predictions of the perforation limit and damage morphology of the laminate subjected to impacts at different energies. The experimental penetration threshold of 89 J was estimated as 79 J by the strain-rate-sensitive models, which was more accurate compared to 52 J predicted by the strain-rate-insensitive model. Additionally, the coupling between interlaminar and intralaminar failure modes in the models led to a more accurate prediction of the delamination area when considering the rate sensitivity.

Automated summary

This study presents a novel numerical model for predicting the perforation limit and damage morphology of composite laminates by considering strain rate sensitivity. The results demonstrate that considering strain rate leads to more accurate predictions of the perforation limit and the coupling between interlaminar and intralaminar failure modes improves the prediction of the delamination area.