A numerical model to simulate bond-slip behavior of composite beam bridges is introduced in this paper. Assuming a linear bond stress-slip relation along the interface between the slab and girder, the slip behavior is implemented into a finite element formulation. Adopting the introduced model, the slip behavior can be taken account even in a beam element which is composed of both end nodes only. Governing equation of the slip behavior, based on the linear partial interaction theory, can be determined from the force equilibrium and a constant curvature distribution across the section of a composite beam. Since the governing equation for the slip behavior requires the moment values at both end nodes, the piecewise linear distribution of the constant bending moment in an element is assumed. Analysis results by the model are compared with numerical results and experimental values, and load-displacement relations of composite beams were then evaluated to verify the validity of the proposed model.