Composite cylindrical shells belong to typical basic structure in engineering, whose buckling and postbuckling behaviours determine the stability and safety of overall structures under complex service loadings. In this paper, a nonlinear mechanical model for anisotropic laminated composite cylindrical shells with general distributed initial imperfection is derived to analyse the buckling and postbuckling behaviours under torsional loads, in which the general distributed imperfections are represented by a generic two-dimensional imperfection function reconstructed by using measured data. Based on a higher order shear deformation shell theory, the governing equations are established with consideration of von Ka<acute accent>rma<acute accent>n-Donnell-type of kinematic nonlinearity and the coupling effects of extension/twist, extension/flexural and flexural/twist. By using a two-step perturbation technique, a new set of explicit solutions is obtained to determine the buckling loads and postbuckling equilibrium paths of the torsional composite cylindrical shells. Furthermore, the internal quantitative relationship between the shear stress with torsional characteristics along with an associate compressive stress of anisotropic laminated shell is for the first time obtained. Numerical analyses are performed to validate the accuracy and effectiveness of the present explicit solutions, and obtain the postbuckling response of different types of imperfections, anisotropic laminated cylindrical shells with different values of shell parameters.