Numerical simulation is a useful tool in investigating the loading performance of rock bolts. The cable structural elements (cableSELs) in FLAC3D are commonly adopted to simulate rock bolts to solve geotechnical issues. In this study, the bonding performance of the interface between the rock bolt and the grout material was simulated with a two-stage shearing coupling model. Furthermore, the FISH language was used to incorporate this two-stage shear coupling model into FLAC3D to modify the current cableSELs. Comparison was performed between numerical and experimental results to confirm that the numerical approach can properly simulate the loading performance of rock bolts. Based on the modified cableSELs, the influence of the bolt diameter on the performance of rock bolts and the shear stress propagation along the interface between the bolt and the grout were studied. The simulation results indicated that the load transfer capacity of rock bolts rose with the rock bolt diameter apparently. With the bolt diameter increasing, the performance of the rock bolting system was likely to change from the ductile behaviour to the brittle behaviour. Moreover, after the rock bolt was loaded, the position where the maximum shear stress occurred was variable. Specifically, with the continuous loading, it shifted from the rock bolt loaded end to the other end.