Efforts are ongoing to characterise a comprehensive resistance function for cylindrical shells under uniform bending, a ubiquitous structural system that finds application in load-bearing circular hollow sections, tubes, piles, pipelines, wind turbine support towers, chimneys and silos. A recent computational study by Rotter et al. (2014) demonstrated that nonlinear buckling of perfect elastic cylinders under bending is governed by four length-dependent domains-'short', 'medium', 'transitional' and 'long'-depending on the relative influence of end boundary conditions and cross-section ovalisation. The study additionally transformed its resistance predictions into compact algebraic relationships for use as design equations within the recently-developed framework of Reference Resistance Design (RRD). This paper extends on the above to present a detailed computational investigation into the imperfection sensitivity of thin elastic cylindrical shells across the most important length domains, using automation to carry out the vast number of necessary finite element analyses. Geometric imperfections in three forms-the classical linear buckling eigenmode, an imposed cross-section ovalisation, and a realistic manufacturing 'weld depression' defectare applied to demonstrate that imperfection sensitivity is strongly length-dependent but significantly less severe than for the closely-related load case of cylinders under uniform axial compression. The axisymmetric weld depression almost always controls as the most deleterious imperfection. The data is processed computationally to offer an accurate yet conservative lower-bound algebraic design characterisation of imperfection sensitivity for use within the RRD framework. The outcomes are relevant to researchers and designers of large metal shells under bending and will appeal to computational enthusiasts who are encouraged to adopt the automation methodology described herein to explore other structural systems.