Abstract In this article, the analytical buckling load coefficient formula for a cylinder with circumferential thickness variation subjected to varying external pressure is established for the first time by developing a quadratic perturbation method. Based on the presented formula, some specific examples that consider either circumferential shell thickness variation or varying lateral pressure are studied and compared with those in the published literature. First, classical cosine thickness variation is analyzed. The computed results of buckling load coefficients indicate that most of differences between the proposed formula and finite element (FE) analyses in previous literature are less than 5%. Second, wind-type pressure load is discussed, and the maximum difference of predicted buckling load coefficients between our formula and results in published article is 2.3%. Third, both circumferential thickness variation and linear liquid pressure are analyzed by the proposed analytical formula and validated by the Galerkin method when both thickness and load variations are small. Therefore, accuracy and reliability of the established formula are validated. For the purpose of engineering application, buckling of a circular shell with local circumferential thickness variation due to weld shrinkage under wind-type pressure is analytically investigated. The effects of thickness variation amplitude, load parameter, half angle of thickness variation, and shell dimensions on buckling load coefficients are studied in detail.