The paper deals with determination of flexural resistance buckling curves for welded I-section steel members made of high strength steel (S 690). In the paper the previously proposed BF analytical model is used for approximation of FEM results obtained using moderately large deformation shell theory and ABAQUS/Standard software. Final formulation of flexural resistance buckling curves is possible through the use of the Merchant-Rankine-Murzewski approach adopted extensively in the authors’ previous papers. For nonlinear optimization, which is needed for analytical model parameters determination, the Wolfram Mathematica package is used. Obtained results for S 690 steel are presented against the results for S 355 steel.