This paper introduces a numerical model to simulate bond-slip behavior in composite beam bridges. Based on a linear bond stress-slip relation along the interface, the slip behavior is implemented into a finite element formulation. An adoption of the numerical model makes it possible to consider the slip behavior even in a beam element defined by two nodes only. The slip behavior can be determined from the force equilibrium, a constant curvature distribution and the piecewise linear distribution of the bending moment. Finally, Correlation studies between numerical results and experimental values were conducted to verify the validity of the proposed model. (C) 2010 Elsevier Ltd. All rights reserved.