Abstract
This work presents a 3D progressive damage model based on Puck’s failure theory and linear damage evolution in fiber-reinforced plastic (FRP) laminates. It includes shear nonlinearity, in situ strengths, equivalent stress–strain, and mixed-mode fracture energy, and is implemented in abaqus/explicitTM through VUMAT subroutine. Various test cases were performed to validate the model and demonstrate its applications. The shear nonlinearity test shows that transverse compression retards matrix microcracking while transverse tension accelerates it. The open hole tension (OHT) test of laminates shows that delamination initiates around the holes and free edges, spreads the most, and propagates in different directions at different interfaces. Later, interfiber damage in 45 deg or −45 deg plies initiates and spreads at a slight inclination to the tip of the hole. Finally, fiber damage in 0 deg plies initiates at the tip of the hole, spreads the least, and propagates perpendicular to the loading direction. The ply-blocked laminates show around 30% higher strength and fracture strain than non-ply-blocked laminate due to delay in damage propagation, and are less sensitive to the hole size. Accordingly, their OHT strength reduces by 14.3% as opposed to 21.14% in the non-ply-blocked laminates, when the hole size increases from 6 to 9 mm. The damage location, magnitude, and propagation were corroborated with experimental findings in the literature.