Nowadays, conventional materials have been progressively replaced by composite materials in a wide variety of applications. Particularly, fiber reinforced composite laminates are widely used. The appropriate design of elements made of this type of material requires the use of constitutive models capable of estimating their stiffness and strength. A general constitutive model for fiber reinforced laminated composites is presented in this paper. The model is obtained as a generalization of classical mixture theory taking into account the relations among the strains and stresses in the components and the composite in principal symmetry directions of the material. The constitutive equations for the laminated composite result from the combination of lamina constitutive equations that also result from the combination of fibers and matrix. It is assumed that each one of the components are orthotropic and elastoplastic. Basic assumptions of the proposed model and the resulting equations are first presented in the paper. The numerical algorithm developed for the implementation in a three-dimensional (3D) finite element nonlinear program is also described. The paper is completed with application examples and comparison with experimental results. The comparison shows the capacity of the proposed model for the simulation of stiffness and strength of different composite laminates.