Composite cylindrical shells and panels are widely used in aerospace structures. These are often subjected to defects and damage from both in-service and manufacturing events.Delamination is the most important of these defects. This paper deals with the computational modelling of delamination in isotropic and laminated composite cylindrical shells. The use of three-dimensional finite elements for predicting the delamination buckling of these structures is computationally expensive. Here combined double-layer and single-layer of shell elements are employed to study the effect of delamination on the global load-carrying capacity of such systems under axial compressive load. It is shown that through-the-thickness delamination can be modelled and analysed effectively without requiring a great deal of computing time and memory. A parametric study is carried out to study the influence of the delamination size, orientation and through-the-width position of a series of laminated cylinders. The effect of material properties is also investigated. Some of the results are compared with the corresponding analytical results. It is shown that ignoring the contact between the delaminated layers can result in wrong estimations of the critical buckling loads in cylindrical shells under compressive load.