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, and shear test results. The outcomes of this study indicated that both theoretical, numerical prediction values agree well with experimental verification results confirming the validity of the proposed methodology in providing a reliable reference for structural design of pultruded fiber reinforced polymeric (FRP) composite structures.