Abstract: Composite laminated cylindrical shell is a common load-carrying structure. Delaminations may arise from its fabrication, transport and service, which will influence the load-carrying capacity of the cylindrical shell. So, it is necessary to model it correctly to investigate the load-carrying capacity of the delaminated cylindrical shell. In this paper, first, an axial-compressive cylinder shell including arbitrary delaminations in longitudinal and thickness and spanning the entire circumference was divided into multiple sublaminated shells whose axial and circumferential displacements were simulated with a cubic function of the thickness coordinate and a trigonometric series of circumferential coordinate. Then, the variational principle was applied to obtain the governing equations and determining conditions. Finally, the governing equations and determining conditions were described as state-space forms and state-space scheme was used to solve the problem.