This paper investigates the nonlinear analysis of concrete-filled steel composite columns subjected to axial loading to predict the ultimate load capacity and behaviour of the columns. Finite element software LUSAS is used to conduct the nonlinear analyses. The accuracy of the finite element modelling is verified by comparing the result with the corresponding experimental result reported by other researchers. Nonlinear analyses are done to study and develop different shapes and number of cold-formed steel sheeting stiffeners with various thicknesses of cold-formed steel sheets. Effects of the parameters on the ultimate axial load capacity and ductility of the concrete-filled steel composite columns are examined. Effects of variables such as concrete compressive strength $f_c$ and cold-formed steel sheet yield stress $f_{yp}$ on the ultimate axial load capacity of the columns are also investigated. The results are shown in the form of axial load-normalized axial shortening plots. It is concluded from the study that the ultimate axial load capacity and behaviour of the concrete-filled steel composite columns can be accurately predicted by the proposed finite element modelling. Results in this study demonstrate that the ultimate axial load capacity and ductility of the columns are affected with various thicknesses of steel sheets and different shapes and number of stiffeners. Also, compressive strength $f_c$ of the concrete and yield stress $f_{yp}$ of the cold-formed steel sheet influence the performance of the columns significantly.