The buckling of axially‐loaded, imperfect, circular, cylindrical, thin shells is studied. A solution methodology, based on two different shell theory approximations to the kinematic relations, is developed and described. The two approximations are those of Sanders and Donnell. Critical loads (limit point loads) are obtained for two shell geometries, one isotropic and one laminated. Several parametric studies are performed in order to establish the accuracy of the Donnell approximation. This is accomplished by comparison to the results obtained by the Sanders kinematic relations, which are considered to be very accurate for thin shells. The parametric studies include the effect of cylinder to radius ration on the critical loads, as well as the effect of the radius to thickness ratio. Moreover, the imperfection sensitivity of the two geometries is investigated by both shell theory approximations.