This paper presents the numerical part of the research program on concrete-filled steel columns. Nonlinear, three dimensional FE analysis of axial compression, was conducted using the finite element program ABAQUS. The numerical results were validated through comparison with experimental data in terms of ultimate loading and deformation modes. Modeling related problems such as the definition of boundary conditions, imperfections, concrete-steel interaction, material representation and others are investigated using a comprehensive parametric study. The developed FE models will be used for an enhanced interpretation of experiments and for the predictive study of cases not included in the experimental testing.