A new numerical model based on the spline finite strip method is presented here for the analysis of buckling of built-up columns with and without end stay plates. The channels are modelled with spline finite strips while the connecting elements are represented by a 3D beam finite element, for which the stiffness matrix is modified in order to ensure complete compatibility with the strips. This numerical model has the advantage to give all possible failure modes of built-up columns for different boundary conditions. The end stay plates are also taken into account in this method. To validate the model a comparative study was carried out. First, a general procedure was chosen and adopted. For each numerical analysis, the lowest buckling loads and modes were calculated. The basic or "pure" buckling modes were identified and their critical loads were compared with solutions obtained using analytical methods and/or other numerical methods. The results showed that the proposed numerical model can be used in practice to study the elastic buckling of built-up columns. This model is considered accurate and efficient for the local buckling of short columns and global buckling for slender columns.