Computational micromechanical modeling of transverse tensile damage behavior in unidirectional glass fiber-reinforced plastic composite plies: Ductile versus brittle fracture mechanics approach

A detailed micromechanical finite element analysis methodology is presented to predict the transverse tensile (fiber perpendicular) failure behavior of a unidirectional (UD) glass fiber-reinforced plastic composite ply. In order to understand the constituent-level stress–strain and damage behavior, finite element analysis is accomplished using representative volume element (RVE) that consists of random fiber distribution as observed in the microscopic image of an actual composite ply. For modeling the fiber/matrix interface failure behavior, cohesive zone module (cohesive surface/cohesive element) of Abaqus® is used. In order to capture the epoxy matrix stiffness and strength degradation, the following two different approaches are used: (i) initially, the linear Drucker–Prager plasticity model in combination with a ductile fracture criterion is used; (ii) later, a brittle failure approach such as the quadratic normal stress criterion within the framework of eXtended finite element method is used. From the detailed micromechanical analysis of the RVE, it is observed that the initial damage in the RVE occurs in the form of fiber/matrix interface decohesion. With increasing tensile load, interface crack propagates and creates a stress concentration region in the matrix material, adjacent to the crack tip. Further load application causes both interface crack tip and matrix stress concentration to move away from the load application direction. As soon as the interface crack tip reaches approximately 60° to 70° away from the load application direction, the conjunction of the matrix damage with the interface crack leads to the RVE final failure. The predicted average stress–strain curves from the above-mentioned two different epoxy matrix failure criterions (ductile and brittle) correlate very well with the experimental results, indicating that the brittle failure behavior of a UD fiber-reinforced plastic composite ply under transverse tensile load is mainly controlled by the fiber/matrix interface properties.