Glued-in rod is a type of effective connector for timber structures. In the last decades, considerable researches have been conducted on the bond behavior between timber and glued-in rod experimentally and theoretically. This paper presents the numerical simulation on the bond behavior of glulam joints with glued-in steel rod. Glulam is modelled as the orthotropic elasto-plastic material. Adhesive and steel rods are modelled as isotropic materials. Non-linear spring element was employed to simulate the local bond stress-slip relations at different locations along the bonded length. The numerical results indicate that non-linear spring element can be used to simulate the bond behavior of the glued-in steel rod glulam joints effectively. Moreover, The numerical results obtained by considering the local bond relations at different locations is in a better agreement with the experimental results than that obtained based on the bond stress-slip relation at loaded end.