This study proposes to address the stability of roadway surrounding rock under three-dimensional stress field based on the Drucker-Prager criterion and the theory of elastic-plastic mechanics. Specifically, we deduced the boundary equation of plastic zone of roadway surrounding rock under three-dimensional stress field by introducing the axial stress solution, and analyzed the shape and size of plastic zone of surrounding rock under different stress field conditions. We compared the effect of five Drucker-Prager yield criteria on the shape and size of plastic zone under variations of cohesion, internal friction angle, Poisson's ratio and roadway radius in roadway surrounding rock. We then analyzed its implications for engineering practice based on the evolution characteristics. Results show that ① Variations in horizontal stress exerts significant impact on the plastic zone morphology of roadway surrounding rock: too high or too low horizontal stress will lead to a butterfly-shaped expansion in the roadway surrounding rock. ② Variations in axial pressure has little effect on the plastic zone shape of roadway surrounding rock, yet exerts major influence on the plastic zone size of surrounding rock: under the same horizontal stress, different axial stress will lead to a bowl-shaped development in the plastic zone size of surrounding rock. With variations in the cohesion, internal friction angle and roadway radius, we observed consistency in the size variation of plastic zone in roadway surrounding rock under different Drucker-Prager criteria. Therefore, the choice of different Drucker-Prager criteria will affect the calculation results of plastic zone in surrounding rock of roadway.

The morphology of the plastic zone of roadway surrounding rock has an important influence on the failure mode and degree of roadway. In order to explore the evolution of plastic zone morphology under three-dimensional stress field, this paper derives the axial stress expression based on elastic mechanics, and determines the approximate solution method of 3D plastic zone under 3D strength criterion according to the idea of solving the boundary equation of butterfly plastic zone. By determining the surrounding rock stress loading scheme through equal spherical stress p and equal deviatoric stress q with different Lode angles θ_σ, the morphological evolution of the plastic zone under different 3D strength criteria is studied in depth, and the low sensitivity of the criterion for butterfly failure is demonstrated. Based on the butterfly failure theory, the asymmetric deformation failure mechanism and control technology of 160206 return roadway in Yangchangwan are analyzed. The results show that: 1) Under the same p, q and different θ_σ stress loading conditions, the morphology of the plastic zone under the five strength criteria shows the evolution patterns of round, oval and butterfly shapes, and the morphology of the plastic zone of surrounding rock is basically consistent for each strength criterion under the same θ_σ. 2) Under the loading scheme with same stress state and different stress directions, the plastic zone morphology of surrounding rock varies greatly. The shape of the plastic zone is largely determined by the horizontal lateral pressure ratio. The axial lateral pressure has a greater influence on the size of the plastic zone, but less influence on the shape of the plastic zone. 3) Under the influence of superimposed mining, the roof of 160206 return roadway presents asymmetric large deformation and failure. Based on the support idea of butterfly plastic zone, the collaborative support technology of ' asymmetric anchor cable + advanced unit support + borehole pressure relief ' has been applied, and good supporting effect has been achieved.