This paper presents an analytical formula for the second-order normal stresses in flanges of the imperfect I beam-columns with sinusoidal web corrugation. To deduce the analytical expression, the beam is separated into three parts. The flanges are modelled as thin-walled beam-columns and the corrugated web is modelled as a thin equivalent orthotropic plate with uniform thickness. The indispensable elastic critical loads are derived from the total potential strain energy of the imperfect beam-column using an energy method. To assess the accuracy of the proposed analytical formula, full shell finite element models were created in Abaqus software. The result of the benchmark example showed that the derived expression is accurate and applicable for geometrically nonlinear analysis of beam-columns with sinusoidal web corrugation. Based on the proposed analytical expression, the Ayrton-Perry type solutions for corrugated web columns and beams are also presented in this paper.