In recent years, an increasing interest in designing magnesium biomedical implants has been presented due to its biocompatibility, and great effort has been employed in characterizing it experimentally. However, its complex anisotropic behaviour, which is observed in rolled alloys, leads to a lack of reliable numerical simulation results concerning residual stress predictions. In this paper, a new model is proposed to focus on anisotropic material hardening behaviour in Mg base (in particular AZ31B as a representative alloy) materials, in which the particular stress cycle involved in Laser Shock Processing (LSP) treatments is considered. Numerical predictions in high extended coverage areas obtained by means of the implemented model are presented, showing that the realistic material’s complex anisotropic behaviour can be appropriately computed and—much more importantly—it shows a particular non-conventional behaviour regarding extended areas processing strategies.