In this paper, analytical and numerical homogenization methods are proposed to effectively simulate the macroscopic characteristics of a pultruded composite lamina. A continuum damage model was implemented via user material subroutine to model fiber failure, while the Mohr-Coulomb plastic criterion is employed to model matrix damage. In order to simulate the damage of the fiber-matrix interface, the relationship between traction and displacement is established. The proposed theoretical and numerical models were verified by tensile, compressive,...