In this study, a progressive damage model was developed for the mechanical response and damage evolution of carbon fiber stitched composite laminates under low-velocity impact (LVI). The three-dimensional Hashin and Hou failure criteria were used to identify fiber and matrix damage. The cohesive zone model was adopted to simulate the delamination damage, combined with the linear degradation discounting of the equivalent displacement method to characterize the stiffness degradation of the material, and the corresponding user material subroutine VUMAT was coded. The finite element analysis of the LVI of stitched composite laminates under different energies was finished in Abaqus/Explicit. Furthermore, the simulation predictions matched well with the results of the experimental tests. Based on this, composite laminates' mechanical response and damage forms with different thicknesses and stitch densities were analyzed. The findings show that the main damages of composite laminates were matrix tensile damage and delamination. The stitching process could improve the impact tolerance of composite laminates, inhibiting delamination and reducing the area of the delamination damage. The higher the density of the stitching, the more noticeable its inhibition would be. The thickness of the laminate also had a more significant effect on the damage to the laminate. Thin plates were more prone to matrix tensile damage due to their lower flexural rigidity, whereas thick plates were more susceptible to delamination because of their higher flexural rigidity.
Keywords: low-velocity impact; mechanical response; numerical simulation; stitching composite laminate; thickness.