In this paper, a numerical model for the analysis of bond-slip occurring in concrete filled steel tube (CFT) columns is introduced. Unlike the classical bond-link or bond-zone element using double nodes, the introduced model considers the bond-slip effect without taking double nodes by incorporation of the equivalent steel stiffness. Upon constructing the equations system on the basis of the force equilibrium and the compatibility condition at each nodal point, the deformations in a steel tube are determined, followed by evaluation of the bond-slip. Moreover, as a part of solving the equations system to evaluate the slip behavior, the mechanical properties for steel and bond-slip have been changed and updated through an iteration procedure. Finally, the validity of the introduced numerical model is verified by comparing the experimental data with the analytical results for CFT columns subjected to axial force and bending moment.