A progressive failure analysis procedure for composite laminates is completed in here. An anisotropic plastic constitutive model for fiber-reinforced composite material is implemented into computer program for a predictive analysis procedure of composite laminates. Also, in order to describe material behavior beyond the initial yield, the anisotropic work-hardening model and subsequent yield surface are implemented into a computer code, which is Predictive Analysis for Composite Structures (PACS). The accuracy and efficiency of the anisotropic plastic constitutive model and the computer program PACS are verified by solving a number of various fiber-reinforced composite laminates with and without geometric discontinuity. The comparisons of the numerical results to the experimental and other numerical results available in the literature indicate the validity and efficiency of the developed model.