Rock bolting plays an important role in different geo-engineering applications and its numerical modelling is crucial for the analysis and design of rock structures. Continuum modelling simulation of bolt-reinforced rock masses requires specific techniques to properly model the reinforcement system and its interaction with the rock mass, which often exhibits a nonlinear softening/brittle response. In this context, strain localization might occur, which, in turn, may affect numerical convergence and the quality of results. This paper presents some advanced numerical techniques implemented in PLAXIS to overcome the abovementioned challenges. Firstly, a regularization technique is implemented for an extended version of the Hoek-Brown failure criterion with strain softening. Secondly, the formulation of the structural bolt element interacting with the rock mass is developed. Finally, the robustness and accuracy of these techniques are discussed via a numerical example of a typical underground mining excavation problem.