A constitutive model is constructed to consider the resin matrix post-yield softening and progressive hardening behaviors. A user-defined material mechanical behavior (UMAT) subroutine is created

then the non-linear three-dimensional finite element analysis on the tensile processes of multi-fiber composites is conducted. The approximate 45° shear bands emanating from the matrix crack tip are found

being coincided with the experimental observations. The shear stress on the adjacent intact fiber/matrix interface is strongly influenced by the shear band and thus the stress concentration factor (SCF) changes obviously in the adjacent fibers. The distinct stress redistribution in the adjacent intact fibers implies the significant effect of the shear bands on the progressive fiber fracture initiation. As the inter-fiber spacing increases

the peak value of the SCF in the adjacent intact fiber decreases

whereas the overload zone becomes wider. The research has provided a helpful tool to evaluate the failure of fiber composites and optimize the composite performance through the proper selection of resin matrix properties and fiber volume fraction.