This paper deals with an introduction of a numerical model to simulate bond-slip behavior in composite beams. The slip behavior is implemented into a finite element formulation and an adoption of the introduced numerical model makes it possible to take into account the slip behavior even in a beam element whose configuration is defined by both end nodes only. Correlation studies for bond-slip effect between numerical results were conducted, and a cyclic loaddisplacement relation of a composite beam is then evaluated to verify the validity of the proposed model.