Abstract The objective of this paper is to develop a constrained shell Finite Element Method (cFEM) based on a force approach for elastic buckling analysis of thin-walled members. The new cFEM is able to separate the general deformation of thin-walled members into the three fundamental deformation mode classes, namely Global (G), Distortional (D), and Local (L), to enable the modal decomposition and identification. In this paper, four force-based mechanical criteria are defined to separate these mode classes. These mechanical criteria are implemented with the general shell finite element formulation without any special treatment of the element formulation. The constraint matrices of the G, D, and L mode classes are then constructed. Numerical examples are presented to demonstrate the capabilities of the new cFEM in modal decomposition and identification. In particular, the modal decomposition and identification results are compared with the cFSM solutions. Applicability of the new cFEM to other shell FE formulation and different loading conditions are illustrated. All these numerical examples demonstrate the potential of the developed cFEM in taking advantages of the modeling capability of the existing shell FE method.