The main aim of the study is an assessment of models suitability for steel beams made of thin-walled cold-formed sigma profiles with respect to different numerical descriptions used in buckling analysis. The analyses are carried out for the sigma profile beam with the height of 140 mm and the span of 2.20 m. The Finite Element (FE) numerical models are developed in the Abaqus program. The boundary conditions are introduced in the formof the so-called fork support with the use of displacement limitations. The beams are discretized using S4R shell finite elements with S4R linear and S8R quadratic shape functions. Local and global instability behaviour is investigated using linear buckling analysis and the models are verified by the comparison with theoretical critical bending moment obtained from the analytical formulae based on the Vlasow beam theory of the thin-walled elements. In addition, the engineering analysis of buckling is carried out for a simple shell (plate) model of the separated cross-section flange wall using the Boundary Element Method (BEM). Special attention was paid to critical bending moment calculated on the basis of the Vlasov beam theory, which does not take into account the loss of local stability or contour deformation. Numerical shell FE models are investigated, which enable a multimodal buckling analysis taking into account interactive buckling. The eigenvalues and shape of first three buckling modes for selected numerical models are calculated but the values of critical bending moments are identified basing on the eigenvalue obtained for the first buckling mode.