In this study, an efficient method is developed to investigate the compressive large deflection behavior of unsymmetric composite laminates with multiple through-thewidth delaminations. The analytical method is based on the first order shear deformation theory (FSDT) and its formulation is developed on the basis of the Rayleigh-Ritz approximation technique by the implementation of the simple and complete polynomial series. The method can handle both local deflection of the delaminated sublaminate and global deflection of the whole plate. Also, the contact among sublaminates is investigated. The three-dimensional finite element analysis is performed by using ANSYS5.4 general purpose commercial software, and the results are compared with those obtained by the analytical model.