In order to study the nonlinear behaviors and interactions among the constituents for the composite material structure under the tensile load, multiscale damage model using generalized method of cells (GMC) and a lamina-level progressive damage model were established, respectively, for fiber reinforced composite laminates with a central hole, which were based on the thermodynamic Schapery Theory (ST) at either the micro-level or the lamina level. Once the nonlinear progressive degradation of the matrix material reached the lower limit value for the ST method, matrix failures naturally occurred, the failure of the fiber was determined by the maximum stress failure criterion. For the multiscale progressive damage model, the GMC model consisting of a fiber subcell and three matrix subcells was imposed at each integral point of FEM elements, and the three matrix subcells undergo independent damage evolution. The load versus displacement curves and failure modes of the open-hole laminates were predicted by using the two progressive failure models, and the results were compared with that obtained by the Hashin-Rotem progressive failure model and the experimental results. The results show that the ST based method can obtain the nonlinear progressive damage evolution states and failure states of the composite at both the lamina level and the multiscale level. Finally, the damage contours and failure paths obtained are also presented.

The failure types and ultimate loads for eight carbon-epoxy laminate specimens with a central circular hole subjected to tensile load were tested experimentally and simulated using two different progressive failure analysis (PFA) methodologies. The first model used a lamina level modeling based on the Hashin criterion and the Camanho stiffness degradation theory to predict the damage of the fiber and matrix. The second model implemented a micromechanical analysis technique coined the generalized method of cells (GMC), where the 3D Tsai–Hill failure criterion was used to govern matrix failure, and the fiber failure was dictated by the maximum stress criterion. The progressive failure methodology was implemented using the UMAT subroutine within the ABAQUS/implicit solver. Results of load versus displacement and failure types from the two different models were compared against experimental data for the open hole laminates subjected to tensile displacement load. The results obtained from the numerical simulation and experiments showed good agreement. Failure paths and accurate damage contours for the tested specimens were also predicted.