Research on the Deformation Mechanism of Mining Roadway Stratiform Surrounding Rock with Nonuniform Stress Field

Aiming at the failure problems of mining roadway stratiform surrounding rock inﬂuenced by asymmetrical high stress, ﬁrstly, the rock instability conditions of diﬀerent structural plane angles were summarized based on the mechanical model of rock joints and instability criterion. Secondly, stresses, plastic zone, and displacement of surrounding rock with diﬀerent structural plane angles were calculated and analyzed with the theoretical calculation method or numerical simulation method, respectively. The distribution characteristics of plastic zone of surrounding rock under diﬀerent inﬂuence factors are further studied, and the results show that the lateral pressure coeﬃcient can change the size and shape of plastic zone at the same time. The damage mechanism was analyzed through the above study. Research studies show that the shear and tension failure area of plastic zone were directly determined by the maximum shear stress and the minimum principal stress. The eﬀect of structure planes on stress has two sides, one to destroy the continuity and the other to strengthen discreteness and laddering nature. When the angle between maximum shear stress and structure plane increases, the extended mode and distribution pattern of plastic zone change substantially as well. The #15 haulage roadway is taken as engineering background, and the above research’s results are veriﬁed by the comprehensive analysis. The research results can enrich the butterﬂy plastic zone theory and provide the basis for roadway supporting design.


Introduction
In the process of coal mining, the surrounding rock of the roadway is almost all layered rock, and the deformation and failure of the surrounding rock of the roadway have always been paid attention to by many experts and scholars [1,2]. Layered surrounding rocks include rock and bedding plane, whose physical and mechanical properties have an important impact on the deformation and failure of surrounding rocks; in particular, the bedding has a more significant impact on the stability of surrounding rocks [3,4]. At present, domestic and foreign experts believe that the deformation and failure of roadway surrounding rocks are closely related to the plastic zone, and the morphology and size of the plastic zone directly determine the failure form and degree of roadway surrounding rocks [5,6]. erefore, it is of positive guiding significance to clarify the plastic zone distribution characteristics of surrounding rock in layered roadway for the study of surrounding rock deformation mechanism, support scheme design, roof accident treatment, and many other engineering practices [7,8].
e research on the plastic zone of surrounding rock started from the mechanical analysis of circular roadway by domestic and foreign experts and scholars, of which abundant theoretical calculation results obtained during this period laid a solid foundation for the follow-up research. For example, through the method of mechanical analysis, people derived the famous formula of "Fenner-Kanerite" surrounding rock plastic zone of roadway [9,10], and the influence laws of lateral pressure, intermediate principal stress, internal friction angle of surrounding rock, cohesion, and other factors on the distribution characteristics of surrounding rock plastic zone were clarified [11][12][13][14]. However, due to the limitation of theoretical calculation, it is impossible to make a direct research on the relationship between the plastic zone and deformation and failure of roadway surrounding rock.
With the progress of science and technology, the visualization numerical simulation software can be used to conduct in-depth research on the plastic zone of roadway surrounding rock, and the research results can be combined with the actual engineering situation to some extent. e homogenate of the plastic zone of surrounding rock of roadway has been studied systematically by the numerical simulation method [15][16][17][18][19][20][21][22][23]. In the course of these studies, the plastic zone morphology of surrounding rock with different section shapes, the general process of "plastic point ⟶ plastic ring ⟶ local distortion ⟶ non-uniform expansion ⟶ malignant expansion" of the plastic zone development, and the influence characteristics of temperature on the plastic zone of surrounding rock are obtained. According to the research results of this period, the mechanism of surrounding rock caving and rock burst in the roadway based on the plastic zone were established, which laid a solid theoretical foundation for the study on the stability control technology of surrounding rock in this type of roadway.
At present, the study of the plastic zone of the surrounding rock mainly focuses on three aspects: ① theoretical calculation method of the radius of the plastic zone; ② the influence of cohesion, internal friction angle, section shape, span, excavation, temperature, and other factors on the plastic zone; and ③ the relationship between the shape of the plastic zone and the deformation of the surrounding rock. However, there are relatively few studies on the distribution of plastic zones in layered roadways. erefore, the influence of the lateral pressure coefficient, the shear strength of the bedding surface, and the tensile strength on the plastic zone of the surrounding rock of the layered roadway was studied in this paper. In addition, the distribution law and deformation characteristics of the plastic zone of the surrounding rock of layered roadways under different inclination angles are analyzed.

eoretical Analysis of the Failure of Layered Surrounding
Rock.
e surrounding rock of layered roadway has many vertical or horizontal bedding planes, showing obvious anisotropy and heterogeneity. In this case, there is a single bedding surface, a single group of bedding surfaces, and multiple groups of bedding surfaces, as shown in Figure 1. Figure 2 shows the mechanical model when a rock mass is cut by a single bedding plane AB. In such case, there are two situations that can cause rock mass destruction: one is that the rock part is damaged by tension or compression shear, and the other is the failure of the bedding surface under tension or compression.
According to the structural plane mechanics model shown in Figure 2, the following formula can be obtained: where σ and τ represent the normal stress and shear stress on the structural surface, respectively. σ 1 and σ 3 indicate the maximum and minimum principal stress, respectively. β signifies the angle between the bedding surface and the maximum principal stress direction. e strength curve of the structural surface obeys the Coulomb criterion, and the shear stress of the structural surface can be expressed as where C j and b j represent the cohesion and internal friction angle of the structural surface, respectively. Equation (3) can be obtained according to equations (1) and (2).
e angle between the bedding surface and the maximum principal stress direction can be summarized as When the principal stress of a certain point in the rock satisfies formula (1), it means that the point in the rock is in a state of limit equilibrium. When the principal stress on the bedding plane AB satisfies formula (2), it indicates that the principal stress on the bedding plane AB is in a state of limit equilibrium. It can be seen from formula (3) that when the bedding surface enters the plastic state, σ 1 � f(β). erefore, the strength of the surrounding rock mass is not a constant but changes with the inclination of the bedding surface.
When β � 90 ∘ or β � b j and σ 1 ⟶ ∞, it means that the rock is impossible to break along the bedding plane at this moment. However, it is also impossible for σ 1 to tend to infinity. Under this condition, it will be destroyed along a  certain bedding plane in the rock. en the strength of the rock is the strength of the rock mass. From the perspective of Mohrʼs circle (Figure 2(b)), line A is the strength envelope of the rock, line B is the strength envelope of the bedding plane, and M 2 represents the stress state of a certain structural plane. If the stress state of a certain section of the rock is exactly at the position where line A and the Mohr circle are tangent (point M 0 ), then the rock will enter plastic failure along the section; otherwise, it will not enter the plastic failure state. e point M 1 is exactly on the strength envelope of the bedding plane, which means that the rock mass will enter a shear plastic state along the bedding plane. e point M 2 is below the strength envelope B of the bedding plane, which means that the rock mass will not enter a shear plastic state along the bedding plane. It can be seen that the failure condition of the bedding surface is that the strength limit curve of the bedding surface intersects the Mohr circle, which is obviously different from the failure situation, where the Mohr circle and the strength limit curve are tangent to the general homogeneous body.
According to equation (3), the uniaxial or triaxial compressive strength of the rock mass is controlled by the angle between the loading direction and the bedding plane. If there is a group of bedding planes with the same properties in the same kind of rock, when the maximum principal stress is perpendicular to the bedding plane, the strength of the rock mass has nothing to do with the bedding plane. At this time, the rock mass strength is the rock strength. When β satisfies formula (4), the rock mass will be destroyed along the bedding plane, and the strength of the rock mass at this time is the strength of the bedding plane. When the maximum principal stress is parallel to the bedding plane, the rock mass will be destroyed due to the lateral expansion of the bedding plane. At this time, the strength of the rock mass will be between the above two conditions. It can be seen from the above analysis that the size of the inclination of the bedding plane has a significant impact on the state (shear or tension) of the surrounding rock of the roadway. e size and shape of the plastic zone of the surrounding rock determine the failure mode and degree of the surrounding rock of the roadway. e dimensions of the roadway section are length × width � 4 × 3 m. e number of grid divisions is between 45,000 and 60,000. e number of layers is between 30 and 35. 336 monitoring points were arranged at equal intervals (1 m) in the horizontal and vertical directions along the center of each roadway section. e upper boundary of the model is a stress boundary, and the remaining boundaries are displacement boundaries. Set the simulation depth to 600 m and the lateral pressure coefficient to 0.3. e interface command in the FLAC3D program was used to establish the structural plane, as shown in Figure 3.

Numerical Simulation of Deformation and
Because the number of mining roadways occupies the absolute advantage of all roadways and the mining roadways are mostly coal roadways, the mechanical parameters of coal seam No. 5 and the mechanical parameters of the original closed structural surface of coal seam No. 5 were assigned to the rock and structural surface characteristics in the model, respectively. e two parts of mechanical parameters come from the laboratory of the Institute of Coal Science and Technology, the Engineering Quality Inspection Center of the Yangtze River Scientific Research Institute of the Ministry of Water Resources, and the mechanical parameters of adjacent mines are given in Table 1. Shock and Vibration

Roadway Surrounding Rock Stress.
According to the failure criteria and the above analysis, the shear and tensile failure of the surrounding rock mass of the roadway are closely related to the principal stress and shear stress. erefore, to clarify the generation mechanism and distribution characteristics of the plastic zone, it is necessary to clarify the distribution law of maximum shear stress and maximum principal stress. According to the knowledge of rock and soil mechanics, the maximum shear stress affects the shear failure of the surrounding rock, and the maximum principal stress can determine whether there is tensile stress in the surrounding rock (the smaller the value is, the easier it is to produce tensile stress). Because the tensile strength of the surrounding rock mass of the roadway is generally much lower than the compressive strength, it is very important to grasp the maximum principal stress of the surrounding rock mass of the roadway. It is worth noting that the requirements for numerical symbols in the software are opposite to the actual ones, so the minimum principal stress analysis of the software calculation results should be used in the analysis.

Maximum Shear Stress.
Analyzing the calculation results of the seven models, it can be seen that the low stress  part of the maximum shear stress (less than 2 MPa) is closely distributed in a certain area above and below the roof and floor of the roadway, and the closer to the edge of the roof and floor, the lower the stress value. e value of this range is difficult to make the rock enter the state of the plastic zone, so the maximum shear stress of the two sides should be analyzed. From Figure 4(a), in the case of no structural surface, the high stress part (above 6.5 MPa) of the maximum shear stress is distributed in a butterfly shape vertically and symmetrically on both sides of the roadway. e 8.5-9 MPa part is distributed about 3 m from the two sides of the roadway. e "butterfly leaf" part reaches the farthest 15 meters from the top corner of the roadway. Figure 4(b) shows that, in the case of horizontal stratification, the overall stress value distribution interval changes a little. e butterfly shape of the maximum shear stress and the edges of different stress values appear sharpened. e "butterfly body" part of the low-stress area of the roof and floor of the roadway is broken into two parts, and the lower part close to the floor is disordered. Figure 4(c) shows that the butterfly shape of the maximum shear stress and the edge sharpening undulations of different stress values are more obvious. Moreover, there are high stress concentration areas of about 8 MPa at the two bottom corners of the roadway. Since the angle between the maximum principal stress and the structural plane is 60°, which satisfies β � (π/4) + (b j /2) � 60.5°, when the inclination angle of the structural plane is 30°, Figure 4(d) shows that the butterfly shape of maximum shear stress sharpens serious water splash. A high stress concentration zone of 7-8 MPa appears in the four corners of the roadway. Figure 4(e) shows that the high stress interval of the maximum shear stress is smaller than that of Figure 4 ere are 6-7 MPa high stress concentration areas at the two corners of the roadway roof. Figure 4(f ) shows that the butterfly shape of the maximum shear stress has changed significantly. e "butterfly leaf" on the lower left is missing a high-stress area, and the stress values in the upper left and lower right parts of the highstress area have sharp edges. ere are 6-7 MPa high-stress concentration areas in the lower left and upper right corners of the roadway. Figure 4(g) shows that the butterfly shape of the maximum shear stress and the edges of different stress values are also sharpened. However, the most notable feature is that the high-stress "butterfly leaf" part has a longer and narrower oblique direction and is obviously symmetrically distributed on both sides of the roadway.

Minimum Principal Stress.
From the simulation results, the high stress part of the minimum principal stress (greater than 18 MPa) is distributed symmetrically in a fan shape at a distance of 2 m away from the two sides of the roadway, and the closer to the fan shaft, the higher the stress value. It is difficult for the value in this area to have tensile stress to make the rock enter a tensile plastic state. erefore, the minimum value of the minimum principal stress of the roof and floor should be mainly analyzed.
It can be seen from Figure 5(a) that the low-stress part of the minimum principal stress is approximately triangular and symmetrically distributes within 7 m of the roof and floor. e stress within 1.2 m of the roof is below 3 MPa, and the farther away from the roof, the lower the stress value. Figure 5(b) shows that there is an obvious ladder phenomenon between the layers, which may be due to the delamination in the low-stress part. e stress is below 3 MPa in the 2.2 m range. Figure 5(c) displays that the lowstress area below 3 MPa is distributed at the upper right corner of the roof and the lower left corner of the floor with a local tip (1 m). e stress within 2.5 m is 6 MPa. Figure 5(d) shows that the high-stress parts of the two sides have higher dispersion.
e phenomenon of water splashing between layers is more prominent. e low-stress area (within 5 MPa) is concentrated within 1.5 m of the delamination tip. Figure 5(e) shows that the high-stress parts of the two sides are further reduced and staggered up and down along the inclined layer. e low-stress area (within 5 MPa) of the roof and floor is located in the middle of the roof and floor in the form of a single peak. A small area of low stress appears in the middle of the left side and the bottom of the right side along the delamination tip. Figure 5(f ) shows that the shape of the minimum principal stress has changed significantly, and the high-stress regions of both sides are sharply reduced. e low-stress area (within 5 MPa) is mainly distributed within 2 m of the roof, and the delamination tip of the left shoulder corner is within 1.5 m of the floor. Figure 5(g) shows that the high-stress part of the two sides is almost consistent with the case without structure plane. However, the stress value is more intensive. A narrow and long lowstress distribution area with the same width as the roadway appears on the roof and floor, and the low-stress area within 2 MPa is about 1.5 m. e maximum shear stress and minimum principal stress are closely related to the plastic zone of the surrounding rock of the roadway. e damage form and scale of the surrounding rock of the roadway are directly determined by the plastic zone. erefore, to grasp the deformation and failure mechanism of the surrounding rock of a layered roadway, it is necessary to analyze the relationship between the plastic zone and the deformation of the surrounding rock.

Plastic Zone of Surrounding Rock of Layered Roadway.
e distribution characteristics of the plastic zone of the surrounding rock under different conditions are shown in Figure 6. e causes and distribution characteristics of the plastic zone shape of each model are analyzed here. Figure 6(b) shows that the boundary line is in the form of a butterfly, and the boundary line matches well with the maximum shear stress of 7 MPa and the minimum principal stress of 1.2 MPa. e two sides of the roadway are mostly shear failure due to the high-stress part of the shear stress. Figure 6(c) shows that the boundary line shape, shear failure area, and tension failure area are consistent with Figure 6(b). e reason is that the maximum principal stress is exactly perpendicular to the structural plane at this time, and the Shock and Vibration strength of the rock mass is the rock strength. Figure 6(d) shows that the boundary line shape still maintains the butterfly shape, but the "butterfly leaf" becomes wider, and the maximum extent of the tension failure zone of the two sides reaches 3 m. Figure 6(e) shows that the butterfly shape is still relatively clear, but the continuity area of the plastic zone is further reduced. e tensile failure distribution of each inclined sublayer has obvious dispersion, which is closely related to the splash phenomenon of the maximum shear stress. e essence is that the angle between the maximum principal stress and the structural plane just satisfies formula (4), which is the most unfavorable angle. e strength of the rock mass is equal to the strength of the structural plane at this time, so large-scale damage occurs along the structural plane. Figure 6(f) shows that the butterfly shape has become smaller and more incomplete; in particular, the tip of the lower right corner has disappeared. e area of the plastic zone along the parallel direction of the structural plane is obviously larger than that of the vertical structural plane. is is because the maximum shear stress range is further reduced and the stress above 7 MPa in the lower right corner disappears. Figure 6(g) shows that the butterfly shape becomes smaller and severely incomplete; in particular, the lower left corner completely disappears. e main reason for this situation is that the 60°i nclination angle changes the maximum shear stress distribution obviously, and the stress in the lower right corner is reduced to a level below 5 MPa. Figure 6(h) shows that the butterfly shape becomes larger and very clear. e "butterfly leaf" became narrower and longer and reached a maximum of 14 m. e reason is that the maximum principal stress is parallel to the structural plane, and the plastic failure is more serious under the action of shear stress.
In short, the shear failure plastic zone on both sides of the roadway ledge changes with the change of the maximum shear stress. e boundary line is consistent with the shear stress value that reaches the critical shear failure strength. e tensile failure plastic zone of the roof and floor of the roadway changes with the change of the minimum principal stress, and its boundary line is consistent with the boundary line of the 2 MPa shear stress value.  6 Shock and Vibration

Deformation and Failure of Surrounding Rock of Layered
Roadway. e vertical and horizontal displacement cloud diagrams of the surrounding rock of the roadway are shown in Figures 7 and 8, respectively. From Figures 7(b) and 7(c), in the two cases of unstructured surface and 0°horizontal layering, the size and shape of the plastic zone are almost the same. However, the range of vertical displacement exceeding 9 cm is larger for the latter. Figures 6(b) and 7(c) also show a similar situation, which shows that when the shape and range of the plastic zone are the same, the more the surrounding rock is cut in layers, the larger the deformation is.
From Figures 6(c)-6(e) and 7(c)-7(e), the range of the continuous plastic zone ("butterfly" shape) gradually decreases with the increase of the inclination angle when the structural plane inclination angle is 0°∼30°. Under the structural plane with 30°, a large range of discrete plastic zones appear in the figure. In addition, the high displacement area of the vertical displacement of the surrounding rock gradually decreases with the decrease of the plastic zone.   Shock and Vibration    show that the plastic zone of the surrounding rock of the roadway lacks the shear failure zone in the lower left corner when the structural plane dip is 60°, which makes this part extremely irregular, so that the vertical displacement at the left and lower left corners is larger than that at the corresponding right side. Figures 6(g) and 8(g) show that this situation also causes the horizontal displacement to form a high-level displacement zone in this area.
Figures 6(h) and 7(h) show that when the structural plane dip is 90°, the "butterfly" shape of the plastic zone of the surrounding rock becomes larger and the "butterfly leaf" becomes thinner and longer. In this case, the distribution range and size of the vertical displacement are between the unstructured surface and the horizontal structure surface. Figures 6(h) and 8(h) show three obvious changes in the high displacement area of horizontal displacement: ① e range of high-level displacement above 3 cm within 3 m of the left and right sides increases significantly. ② Compared with the case where the structural plane inclination angle is 15°∼60°, the high-level displacement area 3 m away from the two sides is obviously larger. ③ e high-level displacement area with four inclination angles reappears in a butterfly shape and the range is increased.

Distribution Characteristics of Plastic Zone of Surrounding Rock of Layered Roadway under Different Influencing Factors
e plastic zone of surrounding rock can intuitively reflect the deformation and failure of surrounding rock. e main factors that affect the distribution characteristics of plastic zone of surrounding rock in layered roadway include stress environment, mechanical properties of rock mass, and bedding plane. In order to clarify the influence laws of different influence factors on the deformation and failure of surrounding rock, it is necessary to clarify the distribution characteristics of plastic zone of surrounding rock under different influence factors. erefore, the dip angle model of 15°bedding plane was taken as an example, and the single factor research method was adopted to carry out numerical simulation study one by one to analyze the deformation and failure characteristics of roadway surrounding rock under different influencing factors.

Lateral Pressure Coefficient.
In order to clarify the influence of stress environment on the distribution characteristics of plastic zone in layered roadway surrounding rock, the lateral pressure coefficients λ � 0.3, 0.6, 1.0, 1.5, 2.0, and 2.5 were set, respectively, for numerical simulation calculation. e calculation results are summarized as shown in Figure 9.
When λ � 0.3, the plastic zone as a whole is symmetrically distributed around the roadway in the form of "X," and  the diagonals are, respectively, inclined upward and downward, which is mainly due to the great nonuniformity of horizontal stress and vertical stress. When λ � 0.6, the overall range of the plastic zone decreases and the "X" shaped tip degenerates, which is mainly because the tip of the plastic zone decreases obviously with the increase of horizontal stress. When λ � 1.0, the shape of the plastic zone evolves into a regular ellipse, and the plastic zone area of roof and floor began to increase, indicating that the deformation and failure of roof and floor began to become serious with the increase of horizontal stress. When λ � 1.5, the plastic zone morphology changes obviously, and the plastic zone area begins to expand significantly on the roof and floor, which is the result of the enhanced action of horizontal stress on the roof and floor. When λ � 2.0 and 2.5, the plastic zone of roadway roof and floor expands greatly under the action of high horizontal stress, and the effect of vertical stress on the two sides is relatively weak. It can be seen that the change of lateral pressure coefficient directly changes the geometry of the plastic zone.

Cohesion and Internal Friction Angle of Rock Strata.
In order to clarify the influence of mechanical properties of rock strata on the distribution characteristics of plastic zone of roadway surrounding rock, the internal friction angles of 43°, 44°, 46°, and 47°and the cohesion forces of 2.81 MPa, 3.01 MPa, 3.21 MPa, and 3.41 MPa were, respectively, set for simulation calculation, and the plastic zones under various conditions were sorted out as shown in Figures 10 and 11. It can be seen from Figure 10 that the cohesion of rock mass has a significant influence on the size of the plastic zone of surrounding rock, but it has no obvious influence on the distribution form of the plastic zone of surrounding rock. When the cohesion of rock mass increases from 2.81 MPa to 3.41 MPa with a gradient of 0.2 MPa, the plastic zone of surrounding rock remains unchanged and the range is gradually reduced, which indicates that the plastic zone can be controlled by increasing the cohesion of surrounding rock of roadway, but its shape cannot be changed.
As can be seen from Figure 11, the internal friction angle of rock stratum has a significant influence on the size of the plastic zone of surrounding rock but has no significant influence on the distribution of the plastic zone of surrounding rock. When the internal friction angle of rock stratum increases from 44°to 46°with a gradient of 1°, the plastic zone of surrounding rock remains unchanged and the range decreases gradually, which indicates that the plastic zone size can be controlled by increasing the internal friction angle of rock stratum, but the shape of the plastic zone cannot be changed.

Internal Friction Angle and Cohesion of Bedding Surface.
In order to reveal the influence laws of friction angle and cohesion of bedding on the distribution of plastic zone of surrounding rock, four conditions of internal friction angle, 15°, 31°, 45°, and 60°, and four conditions of cohesion, 0.1 MPa, 0.6 MPa, 1.1 MPa, and 1.6 MPa, were, respectively, set for simulation calculation. e plastic zones under various conditions are sorted out as shown in Figures 12 and 13. Figure 12 shows that the plastic zone of surrounding rock is closely related to the internal friction angle of the bedding plane. When the internal friction angle of the bedding plane is 15°, the two sharp angles of the plastic zone show obvious dispersion phenomenon. When the internal friction angle of bedding plane increases to 31°, the dispersion phenomenon of two sharp angles in the plastic zone disappears. As the internal friction angle continues to increase, the change in the shape and size of the plastic zone of surrounding rock of roadway is relatively small. is indicates that increasing the friction angle in the roadway bedding plane can reduce the phenomenon of overbedding development in the plastic zone, but it cannot change its shape and size. Figure 13 shows that the plastic zone of surrounding rock presents a slight increase trend with the increase of the cohesion of bedding plane. e increase was relatively modest and no obvious change occurs when it increases to a certain extent. For example, when the cohesiveness of bedding plane is 0.1 MPa, the maximum boundary of plastic zone of surrounding rock is 6.21 m; when the cohesiveness of bedding plane increases to 1.6 MPa, the maximum boundary of plastic zone of surrounding rock is 7.12 m, which only increases by 1 m.
is phenomenon does not mean that increasing the cohesion of surrounding rock strata is increasing the deformation and failure of the roadway. On the contrary, increasing the cohesion of surrounding rock strata strengthens the homogeneity of surrounding rock of two layers and makes the plastic zone of surrounding rock of the roadway expand like in homogeneous body.

Normal Stiffness and Shear Stiffness of Bedding Surface.
Normal stiffness and the tangential stiffness are the key mechanics properties of bedding surfaces, which have important influence on stability of rock mass. In order to show the influence of these two factors on the distribution of plastic zone of surrounding rock directly, four cases of normal stiffness, 5e9 MPa·cm −1 , 10e9 MPa·cm −1 , 15e9 MPa·cm −1 , and 20e9 MPa·cm −1 , and four cases of tangent stiffness, 0.5e9 MPa·cm −1 , 2.5e9 MPa·cm −1 , 4.5e9 MPa·cm −1 , and 6.5e9 MPa·cm −1 , were, respectively, set for simulation calculation.
As can be seen from Figure 14, when the normal stiffness of the plastic zone of surrounding rock of the roadway is 5e9 MPa·cm −1 , the upper and lower left butterflies reach 10.81 m and 7.81 m, respectively, and their full development forms an asymmetrical distribution state with the right. When the normal stiffness is increased to 11e9 MPa·cm −1 , the asymmetry decreases obviously. With the further reduction of the normal stiffness, the butterfly morphology gradually becomes symmetrical, and the reduction of the left plastic zone becomes smaller. When the normal stiffness is increased to 20e9 MPa·cm −1 , the plastic zone morphology changes a little. It shows that the too low normal stiffness of bedding plane can cause the asymmetric distribution of the plastic zone of surrounding rock. Figure 15 shows the influence rule of tangential stiffness on the plastic zone of surrounding rock. When the tangential stiffness is 0.5e9 MPa·cm −1 , the plastic zone of surrounding rock of roadway has obvious bedding development phenomenon. When the tangential stiffness is 2.5e9 MPa·cm −1 , the distribution of the penetration in the plastic zone of surrounding rock disappears. With the further increase of tangential stiffness, the plastic zone of surrounding rock becomes symmetrical. When the tangential stiffness increases to a certain extent, the shape and size of the plastic zone do not change significantly.

Deformation Characteristics of the Roadway.
e surrounding rock deformation of the roadway is very serious due to the influence of high ground stress and structural surface.
is lane is a typical multilayered roadway controlled by inclined structural planes. e deformation and failure characteristics of the roadway are as follows: (1) Serious internal squeezing and deformation of the shoulders of the roadway roof. As shown in Figure 18(a), the shoulder corners of the roof are so severely squeezed that the original rectangular shape becomes an isosceles trapezoid shape. In addition, the failure of the bolt (cable) is prominent, and the roof sinks up to 1.5 m. From Figure 18(b), the layered distribution of the roof is very obvious. (2) e two corners of the roadway floor heave and deform seriously. e surrounding rock near the bottom corner of roadway floor heaves laterally, which makes the original regular right angle become a straight line. e floor heave is serious in the whole roadway, and the most serious floor heave section is 1.2 m. e total deformation of roof and floor is more than 2.5 m.
Borehole imager was used to observe the damage of surrounding rock. Five boreholes with 11 m in depth were arranged on the roof and the two sides of the roadway for observing the damage of surrounding rock. e borehole of the floor was not constructed due to the difficulty of construction and observation. e typical borehole fracture images of surrounding rock in different depths are shown in Figure 19.
In order to visually display the distribution of cracks in different parts of the surrounding rock, the distribution of fractures was redrawn according to Figure 19, and the result was shown in Figure 20.
From Figure 20, the crack development area can be divided into severe crack area, moderate crack area, and mild crack area along the roadway surface to deep. e outline connecting lines of these three areas are approximately saddleshaped. e maximum values of the severe crack area at the two sides, the shoulder corners, and the middle of the roof are 1.7 m, 3.4 m, and 1.3 m, respectively. e maximum values of the mild crack area at the two sides, the shoulder corners, and the middle of the roof are 6.5 m, 10.2 m, and 5.7 m, respectively.
is shows that the destruction depth of the surrounding rock is not the same in different distribution intervals, which directly leads to an increase in the intensity of the deformation of the surrounding rock in different parts. e contour line is consistent with the contour line of the plastic zone when the side pressure ratio is 1 : 3. In particular, the positions of the shoulders of the roadway are like butterfly leaves. e following phenomena exist in the roadway surrounding rock failure zone: shoulder corner > two sides > roof.

Deformation Characteristics of the Roadway.
e failure mechanism of the roadway can be further clarified by analyzing the plastic zone of the surrounding rock of the roadway. In order to clarify the deformation and failure  According to the measured mine pressure data and experience provided by the miner, the stress applied at the top of the model is 2500 × 700 m � 17.5 MPa, and the lateral pressure coefficient λ � 0.34, as shown in Figure 21.
According to the data in the rock mechanics manual and the mechanical parameters of the roof and floor provided by the miner, and combined with numerical simulation inversion, the mechanical parameters of the roadway roof and floor were obtained as shown in Table 2. Mechanical parameters of structural surfaces are consistent with those in Table 1.
According to the original support design drawing provided by the mining party, the numerical simulation support scheme design was carried out. e roof is supported by four rock bolts and three anchor cables, and the left sidewall and the right sidewall are supported by five rock bolts, respectively. e size of the rock bolt is φ × L � 22 mm × 2400 mm, the size of the anchor cable is φ × L � 17.8 mm × 8300 mm, and the preload forces of the bolt and anchor cable are 30 kN e rock bolt interval and row spacing are 800 mm and 700 mm, respectively, and the anchor cable interval and row spacing are 1500 mm and 2000 mm, respectively. e details are shown in Figure 22. Figure 23 shows the failure mode of the plastic zone in the track roadway of coal seam No. 15. e butterfly leaf part of the plastic zone (shear failure) of the roadway develops obliquely at the four corners of the roadway due to the influence of the layered joint surface and the nonuniform stress field. Tension failure is dominant in the positive direction of the roof and floor, which is distributed in a regular triangle pattern. e depth of the plastic zone in the four corners of the surrounding rock is the largest and is accompanied by obvious bedding penetration phenomenon.
With the expansion of the plastic zone, the surrounding rock of the roadway begins to show four corners bulging, roof sinking, floor bulging, and so forth. According to the mechanism of the expansion of butterfly-shaped plastic zone and the failure of surrounding rock, it can be preliminarily considered that the failure of surrounding rock of the track roadway is mainly caused by layered joint surfaces and nonuniform stress field. e displacement nephogram of #15 haulage roadway is shown in Figure 24; large deformation of surrounding rock is the most obvious feature. e maximum displacements of the roof and bottom plate are 0.45 m and 1.13 m, respectively, and the deformation of the left and right sides is about 0.48 m. ere is a positive correlation between the deformation of surrounding rock and the distribution of plastic zone. It can be seen that the deformation and failure of the surrounding rock of the roadway are very serious under the support of the original scheme, which has a good matching degree with the field situation in Figure 18. e field measurement data of nearly 100 days were organized into deformation curves of surrounding rock, as shown in Figure 25. It can be seen from Figure 25 that the deformation of surrounding rock began to increase after 20 days, and the increase of deformation began to decrease after 60 days, but it still kept increasing. After 100 days, the deformation of surrounding rock reaches the maximum value during the monitoring period, and the maximum displacements of floor, roof, left wall, and right wall are 1.35 m, 0.71 m, 0.63 m, and 0.51 m, respectively.
e numerical simulation results are consistent with the field measurement results, indicating that the deformation and failure of surrounding rock are closely related to the plastic zone distribution of surrounding rock.     Figure 24: Numerical simulation of displacement nephogram of #15 haulage roadway.

Conclusions
It is of positive guiding significance to clarify the plastic zone distribution characteristics of surrounding rock in layered roadway for the study of surrounding rock deformation mechanism, support scheme design, roof accident treatment, and many other engineering practices. Some conclusions are summarized as follows: (1) Based on the mechanical model and failure criterion of rock mass single structural plane, the failure conditions of rock mass are theoretically analyzed under different structural plane inclination angles. e uniaxial or triaxial compressive strength of the rock mass is controlled by the direction of loading and the inclination of the structural plane under a single structural plane. e strength of the rock mass is controlled by the angle between the maximum principal stress and the structural plane under a set of structural planes with the same properties.
(2) e shape of the area of maximum shear stress and minimum principal stress directly determines the shape of shear and tension failure in the plastic zone. e structure plane destroys the continuity of stress and strengthens the discreteness and laddering nature of stress.
(3) e horizontal structural plane has almost no effect on the plastic zone of the surrounding rock of the roadway. During the gradual increase of the structural plane inclination angle from 0°to 90°, the overall shape of the plastic zone undergoes a process of decrease, increase, decrease, and increase again. When the inclination angle of the structural plane is 30°, the structural plane of the plastic zone is distributed in star-point shape and the overall dispersion is the strongest. When the inclination angle of the structural plane is 90°, the butterfly shape of the plastic zone is thinner and longer. e deformation of the surrounding rock of the roadway is affected by the shape and size of the plastic zone. (4) Among many factors affecting the plastic zone of surrounding rock, only the lateral pressure coefficient can change the shape and size of the plastic zone of surrounding rock at the same time, the cohesion force and internal friction angle of rock strata can only change the size of plastic zone of surrounding rock, and the mechanical properties of bedding surface can change the size and dispersion degree of the plastic zone of surrounding rock. (5) Different dip angles of bedding plane have obvious influence on the deformation and failure of surrounding rock, which is the joint action of dip angles of bedding surfaces, stress state, and mechanical properties of rock strata. When the angle β between the maximum principal stress and bedding surfaces is closer to π/4 + b j /2 (b j is the internal friction angle of bedding surface), the deformation and failure of surrounding rock are more influenced by the bedding plane. When the angle β between the maximum principal stress and bedding plane is closer to π/2, the deformation and failure of surrounding rock are less influenced by the bedding plane. Data Availability e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e author declares that there are no conflicts of interest.

Authors' Contributions
All the research contents of this paper were independently completed by Yongxin Sun.