An elastoplastic anisotropic damage numerical model for rock is established by adopting Mohr-Coulomb criterion as the shearing yield condition

ultimate strain as the shearing damage criterion

and maximal tensile stress as the tension damage criterion. The staged damage evolution of underground stope is emphatically numerically studied in the construction process with supporting structure and without supporting structure; and the deformation field of underground stope is also numerically calculated on the supposing that rock as an elastic material. The numerical test reproduces the staged damage evolution of underground stope in the construction process; and it provides a means to optimize the design and construction process of underground stope

to effectively prevent and control some engineering accidents such as collapse.