Numerical Simulation of Cracking Failure andWeakening Law of Roadway Surrounding Rock under High Stress

In deep roadway mining, high water pressure causes rock mass cracking and weakens the overall strength, affecting the stability of underground metal mine mining roadways. 'erefore, using a numerical simulation method, we analyzed the strain softening characteristics of rocks after the inflexion point of elasticity and studied the strain distribution and the minimum support resistance under high-pressure groundwater conditions. Considering the groundwater pressure and effective porosity on the strain softening characteristics of the surrounding rocks, we investigated the critical groundwater pressure under which the surrounding rocks would remain stable. Actual engineering verification helps to obtain the supporting characteristic curves under different influencing factors. We found that water pressure and effective porosity are the significant factors that decide the development scope of the plastic zone. 'e more significant the increase of the plastic zone, the more notable the changes in the support curve. Moreover, the plastic zone is likely to occur when the hydraulic head is between 30 and 50m; when the hydraulic head exceeds 50m, it is likely to produce a relaxation zone.


Introduction
e stress behavior caused by groundwater pressure in roadway excavation of nonferrous metal mines is a major factor affecting the stability of surrounding rocks. e surrounding rocks of roadways excavated under high water pressure show obvious softening characteristics and reduced strength [1,2]. Because the surrounding rock of excavated roadways is subject to long-term physical and chemical reactions caused by pore water, its effective stress, rock stiffness, and stability usually decrease when it remains saturated with water [3,4]. Sun et al. carried out triaxial tests under water-saturated conditions to study carbonaceous shale [5]. ey found that the compressive strength of rocks changed as the confining pressure changed. Under low confining pressures, the compressive strength of rocks was affected by the time they saturated with water. Besides, the attenuation curve shows an exponential distribution. However, under high confining pressures, the compressive strength of rocks was not affected by the time they saturated with water [6,7].
Huang et al. obtained the permeability change law and cracking mechanism of deep rock mass under hydraulic pressure through experimental research on the permeability of deep roadway floor rock mass. Huang et al. used the true triaxial experimental system to study the influence of hydraulic fracture propagation and proposed a system of controlled fracture methods such as directional hydraulic slitting and fracturing [8,9]. Zhu et al. studied the mechanism of rock instability in mining roadway by establishing the model and obtained the law of the change after the stress disturbance of surrounding rock [2,10,11]. Jiang et al. analyzed the process of rock burst caused by the breakdown force of high water pressure by studying groundwater's accumulation effect in the mining process of deep roadway. ey put forward the corresponding prevention and control measures [12][13][14]. However, they did not analyze the stressstrain relationship caused by rock burst [15]. Most published studies have researched the elastic-plastic deformation of surrounding rocks in excavated roadways under water-rich conditions by combining theoretical research with numerical calculations. However, since these studies do not consider the effect of strength attenuation after the peak strength of surrounding rocks is reached, the results cannot be applied to actual experiments. Most softening constitutive models of surrounding rocks are linear softening models. e effect of groundwater seepage, especially the effect of high groundwater pressure, is often neglected in structural stress analyses of the surrounding rock mass. Since the water depth of high water pressure mining roadways is generally more than 30 m, and its initial ground stress is relatively low [16,17]; it is necessary to consider the mechanical characteristics of nonlinear strain softening of surrounding rocks in order to truly reflect the stress-strain relationship of deep surrounding rocks.
is is especially necessary for elucidating the stress-strain relationship of locally fractured rocks or soft rocks under high water pressure. erefore, based on the nonlinear softening model of surrounding rocks [18], the present study analyzed the elastic-plastic zone of surrounding rocks in deeply excavated roadways, discussed the various parameters that would influence the mechanical characteristics of surrounding rocks, and determined a relatively accurate law for displacement and stress distribution of the elastic-plastic zone of surrounding rocks, considering the seepage of high-pressure groundwater and different drainage conditions of the mining roadway.  [19][20][21], which can accurately reflect the relationship between peak intensity and strain. e model is described as follows:

Theoretical Analysis of the
where E is the elastic modulus; σ c is the peak strength expressed as σ c � 2ccosϕ/(1 − sinϕ), in which c is the cohesion of the rock and φ is the angle of internal friction; g is a test parameter and g � 7.264 × 10 − 4 σ c 2 ; the value of g differs before and after the peak strength, which is applicable to all rocks, and 0.25 is taken as the inflexion point for the value of g.
In order to test and verify the stress-strain changes and strain softening characteristics of rocks after the peak strength, this study analyzes the test results of four types of rocks as in Table 1.
When φ � 30 ∘ (average), 4 MPa, 5 MPa, 6 MPa, and 7 MPa are taken as values of σ c to illustrate the relationship between σ c and ε; the relationship curves are shown in Figure 1.
According to Figure 1, the stresses and strains of the four rocks show strain softening characteristics after reaching the peak strength. After reaching the peak strength, the stress decreases and the strain increases. e results show that the model has universal applicability to describing the strain softening characteristics of rocks after the peak strength is reached. When this model is applied to actual engineering problems, however, it is necessary to carry out fitting analyses with real rock test data to obtain the actual expressions for g and ε c .

Basic Assumptions for Mechanical Calculation.
(1) It is assumed that the section of the excavated roadway is circular. e area includes an elastic zone and a plastic zone. e radius of the roadway is set as a, the calculated radius of the plastic zone is R, and the calculated radius of the elastic zone is b. e calculation model is shown in Figure 2.
(2) To facilitate calculation, it is assumed that the rock strata above the calculated area of the roadway are not subject to any influence from mining, and the thickness of the overlying strata is h � b − a. (3) It is assumed that the surrounding rock is homogeneous and isotropic and has μ as its Poisson's ratio, φ as its effective porosity, the ratio of the volume of interconnected pores in the rock to the total volume of the rock, the characteristics of porous media, and a lateral pressure coefficient of 1 (λ � 1). e gravity stress is σ b , pore water pressure is p f , groundwater pressure is p 0 , and roadway support force is p a at the center of the roadway. (4) Because the roadway is axisymmetric and subject to uniform pressure, σ 1 � σ θ and σ 3 � σ r . If the excavated roadway is infinitely long, the strain problem can be treated as a plane problem, that is, ε z � ε 2 � 0. (5) It is assumed that groundwater seepage has a stable flow and follows Darcy's law of seepage and hydraulic movement principles [22,23]. (6) It is assumed that the calculation parameters do not change with the groundwater environment after rock burst, and the calculation parameters are invariants in the calculation process [24][25][26].

Distribution Law of Pore Water Pressure in Elastic-Plastic
Zone of Surrounding Rock. Affected by highly enriched groundwater, there is a pressure difference between the inside and the outside of the roadway, forming an elastic zone and a plastic zone. In addition, the groundwater will flow along the stratum fissures and cause seepage [27,28]. It is assumed that the seepage has stable flow and thus follows Darcy's law of seepage. According to the hydraulic movement principles of groundwater in rock strata, the pore water pressure in the elastic-plastic zone presents the following distribution law [29].
(1) In the plastic zone, a ≤ r ≤ R, and thus where p 1 is the pore water pressure on the wall of the roadway and p R represents the pore water pressure in the plastic zone.
(2) e flow velocity r � R, which satisfies the condition of continuous interstitial flow: where ρ is the coefficient ratio, ρ � K p /K e , K p is the permeability coefficient of surrounding rocks in the plastic zone, and K e is permeability coefficient of surrounding rocks in the elastic zone.

Analysis of Stress and Displacement of Surrounding Rock in Elastic
Zone. According to the mechanical equilibrium equation (30), the stress of the surrounding rock mass satisfies the following relation: where ϕ is the pore water pressure coefficient. By deriving (4) and combining it with (2), the following equation is obtained: Combining the plane physical equation and the axisymmetric geometric equation [19], it can be found that ε r � du/dr, ε θ � u/r, and c rθ � 0. e displacement equilibrium equation is described as follows:  Plastic zone where μ is the volume compression coefficient, and u is the radial displacement. By solving (6), the expression for the radial displacement u is obtained: where the two parameters C 1 and C 2 are integral constants. e above integral constants, C 1 and C 2 , can be determined by the displacement boundary conditions in the elastic zone, that is, e expression of stress distribution can be obtained as follows: where ] 1 � μ/1 − μ, and σ R is the radial stress (at the interface between elastic and plastic zones). e displacement calculated by using (9) is expressed as the absolute displacement. However, before the roadway is excavated, there is a certain initial displacement in the surrounding rock, and this initial displacement should be deducted in calculations for practical scenarios. us, the relative radial displacement after deduction is expressed as

Analysis of Stress and Displacement of Surrounding Rock in Plastic
Zone. In the plastic zone of surrounding rocks in mining roadways, the rock has different degrees of compressive deformation. To facilitate displacement analysis, it is assumed that the volume is incompressible in the plastic zone. erefore, u � A/r, where A is an integral constant. According to the assumption made in Section 2.2, when the mining roadway is infinitely long (ε z � 0), the equivalent effect ε i in the plastic zone is replaced by (12), and thus, When r � R, ε i � ε c , thus, e displacement in (12) and the equivalent strain in (14) Since the volume is incompressible in the plastic zone, , and the equivalent stress can be expressed as follows: It is assumed that the ultimate bearing capacity occurs in the complex stress state, and it is not unloaded during the loading process. Consequently, the proportional relationship of the stress components varies with the loading. erefore, the traditional constitutive equation of the plastic stage, i.e., (1), can be generalized and applied by using Henck's total strain theory in plastic mechanics [30,31]. e relationship between the equivalent stress and equivalent strain of the rock mass under complex stress states can be obtained as follows: In addition, under the condition of seepage, the equivalent stress follows the equilibrium equation [32,33], which satisfies the equilibrium equation expressed in By substituting (15) and (16) into (18), and combining the result with (17), the following equation is obtained: 4 Shock and Vibration where F 2 � ϕp R − p 1 /In(a/R), and C is an integral constant. According to the model established by [34,35], the value g takes 0.25 as the inflexion point: the value exceeds 0.25 before the strength peak, and it is lower than 0.25 after the peak.
According to the boundary conditions of the plastic zone, that is, σ p r | r�a � p a , the radial stress of the plastic zone is obtained as follows: In the same way, the relation between σ p θ and σ p r in (16) can be obtained as follows:

Determination of Radius R of Plastic Zone.
When r � R, the mining roadway is symmetrical, and therefore the stress in the elastic zone of surrounding rocks on the other side of the roadway is also in accordance with (16). us, (22) is obtained.
e stress in the plastic zone of the surrounding rocks on the other side of the mining roadway can be obtained by using (19), that is, According to the continuity condition of the elastic zone and plastic zones, that is, σ e R � σ p R , the radius R of the plastic zone can be obtained by numerical calculation.

Interaction between Roadway Surrounding Rock and
Roadway Support. After the roadway is excavated, due to the stress release of surrounding rocks and the effect of Shock and Vibration groundwater, the displacement of the tunnel wall is released, leading to a new stress equilibrium [36,37]. After adding support to the surrounding rocks, the groundwater, surrounding rocks, and lining interact with each other again. To reflect this interaction, lining is regarded as a thick-walled cylinder under a uniform external pressure p a . According to the Lame equation [35,38,39], the following is obtained: Here, r 0 is the inner radius, a is the outer radius, E ' is the elastic constant of the surrounding rock material, u ′ is the elastic constant of the lining material, u a represents the total displacement of the external tunnel wall after adding support, and Δu is the displacement of surrounding rocks released before adding support.

Example Verification and Discussions
is study uses geological data about surrounding rocks of mined roadways and lining in Beishan Lead-Zinc Mine, Guangxi, as calculation parameters. e lithology is limestone. e groundwater level elevation in the study area is 300.87-414.66 m. e groundwater is abundant, and the annual variation range is relatively large, generally around 40 m. Frequent and large water level fluctuations result in high water pressure and negative pressure effect on mining roadways, the surrounding rock, and the main factor of cracking. Table 2 indicates the mechanical parameters involved in this example.
When constructing groundwater seepage prevention and drainage systems in mining roadways, the safety of the roadway is always the top priority. However, it is impossible to avoid groundwater seepage completely. erefore, in order to compare the characteristic curves of surrounding rocks of roadways under different drainage conditions, this study considers the following two drainage conditions by referring to Beishan Lead-Zinc Mine data: p 1 � 0 (full drainage) and p 1 � 1.0p f (no drainage), where p 1 is the roadway water pressure and p f is the full water pressure.

Factors Influencing Radius of Plastic Zone after Rock
Burst.
(1) e mining depth of the roadway is 145 m, and the resistance under support is not considered. Under the two drainage conditions described above, the correlation curve between the groundwater pressure on surrounding rocks p 0 and the radius of the plastic zone R/a is shown in Figure 3. As shown in Figure 3, as the groundwater pressure rises when there is no drainage (i.e., in the water saturation state), the radius of the plastic zone R increases and the bearing capacity of surrounding rocks decreases gradually. When the bearing capacity of the surrounding rocks is completely lost, the curve tends to decline; under drainage conditions, the curve for the radius of the plastic zone R rises before leveling off, and the bearing capacity of surrounding rocks increases before reaching a relatively stable value. e curves under the two working conditions show that the stress of surrounding rocks in the upper part of the roadway will be significantly affected during excavation when there is a large volume of groundwater and changes in groundwater depth are significant. When the volume of groundwater is large, the radius of the plastic zone R is significantly affected by changes in groundwater depth. When the critical depth is reached, the stress of surrounding rocks will decrease rapidly and remain unsteady, and a slight disturbance of the surrounding rocks will lead to water inrush. (2) When the depth of groundwater in the upper part of the roadway is 30 m, there is no support in the roadway (i.e.p a � 0). us, the major factors affecting the radius of the plastic zone R are the depth of groundwater in the roadway and the effective porosity of the surrounding rock. e relationship curves are shown in Figure 4.
According to Figure 4, the radius of the plastic zone R increases as the groundwater depth and effective porosity increase, and the influence of the effective porosity grows as the degree of water drainage increases. e pore characteristics of rocks are different from those of general loose media, and the runoff of water in rocks is mainly layer fissure flow; thus, the effective porosity coefficient is used to modify the curve of effective pore stress change [33,39,40].

Correlation between Roadway Support and Surrounding
Rock. When the height of groundwater in the upper part of the roadway is 70 m and the radius of the excavated roadway is 4 m, the curve of surrounding rock deformation characteristics is drawn according to two different drainage conditions, as shown in Figure 5. Figure 5 shows that the surrounding rock support stress decreases and the displacement gradually increases under the two conditions; the characteristic curves gradually approach the x-axis while moving rightward across the graph. However, when the water is drained in the roadway, the stress is higher than that under the nondrainage condition. Besides, the change in displacement is smaller than that under the nondrainage condition, indicating the stability of surrounding rocks under the nondrainage condition is relatively weak.
Next, the roadways under the two working conditions in Figure 5 are supported by two different linings. e thicknesses of the linings are 0.2 m and 0.4 m, respectively, and the lining material is C 30 concrete. When the displacement of the supporting roadway wall is released, Δu � 3cm, the relation curve of the supporting characteristics is shown in Figure 6. Figure 6 shows that when the volume of water drainage remains the same, as the thickness of the lining increases, the support stiffness and resistance increase, and the   displacement of the roadway wall decreases. However, when the thickness of the lining continues to increase, the role of the supporting structure is not evident in further reducing the displacement of the roadway wall and improving the bearing capacity of surrounding rocks. erefore, when adding support to mined roadways, it is necessary to determine the proper thickness of the lining, and the critical thickness should be determined according to the volume of groundwater drainage. When the critical thickness is balanced with the volume of groundwater drainage, an economical and reasonable lining thickness can be obtained.

Law of Stress Distribution in Surrounding Rock.
When the depth of groundwater in the upper part of the roadway is 30 m, p a � 0, and p 1 � 0, the distribution of tangential stress and radial stress in the surrounding rock of the roadway after excavation presents a certain law of distribution, as described by the curves in Figure 7.
As Figure 7 shows, the ratio of original rock stress to radial stress presents different changing trends at the same point in the surrounding rock of the roadway under drainage and nondrainage conditions. Under the nondrainage condition, the tangential stress of surrounding rocks first increases and then decreases as the drainage volume grows. Under the drainage condition, the tangential stress increases gradually from the inner wall of the roadway to the outside and finally reaches a steady state. e two scenarios described above indicate that the peak stress of surrounding rocks decreases continuously in the expansion process of the plastic zone of surrounding rocks of roadways under the nondrainage condition; meanwhile, this peak stress increases continuously under the drainage condition.

Analysis and Determination of Minimum Support
Resistance of Surrounding Rock. As shear slips continue to occur in the plastic zone in the roadway mining process, the surrounding rock begins to loosen and collapse to form a relaxation zone [39,40], and the initial ground stress of the surrounding rock in the relaxation zone is higher than the normal tangential stress. e relation between the change in the loosening pressure and that in the deformation pressure of surrounding rocks can be reflected by the characteristic curve of surrounding rocks [41,42]. e intersection point of the two curves is where the minimum surrounding rock pressure p min is most likely to occur. It is essential to determine the value of p min in support design. p min can be where r s is the radius of the relaxation zone.
Where r � r s , the condition for continuity is σ p θ � σ b + p 0 . Where r � R, the radial stress σ e R � σ p R is the condition for continuity, and p a in 20, (23), and (24) are replaced by p min (in 25). Matlab can be used to solve for r s and R. When the groundwater depth in the upper part of the roadway exceeds 30 m, the correlation among r s , R, and p min under different thicknesses of overlying strata is shown in Figures 8 and 9.
According to Figures 8 and 9, when the height of the groundwater head in the upper part of the roadway is higher than 30 m, the surrounding rock in the locally soft or fractured section appears in the relaxation zone, and the minimum support resistance of the lining decreases notably. In practice, when the support resistance reaches the minimum value p min , the surrounding rock will undergo ultimate deformation, the pressure of the surrounding rock after loosening is borne by the support structure, and the interaction between the surrounding rock and the support structure is deficient. is negatively influences the system's overall stability and is a significant cause for many highpressure water mining accidents. erefore, this kind of situation should be avoided.

Conclusions
(1) According to the numerical simulation method, the distribution law of stress and displacement of surrounding rock of roadways with deep and highpressure water is studied. Considering the factors influencing the seepage of high-pressure water and different drainage conditions, we derived the analytical solution of the stress and displacement distribution of surrounding rock of the roadway under a complex stress state. (2) During groundwater seepage, when the high-pressure head pressure of groundwater reaches the critical pressure for the stability of the surrounding rock of the roadway, the equilibrium state of the surrounding rock of the mined roadway is destroyed. At this time, surrounding rock easily loses stability and collapse may occur under slight disturbance, causing water inrush in the mine. (3) e extent of plastic zone development is mainly affected by drainage conditions and effective porosity and increases with the increase in drainage and effective porosity. e supporting characteristic curve of fracturing range also changes with the increase in development of this zone. erefore, in the design of waterproof linings and drainage systems during excavation of the tunnel, the seepage flow and the effective porosity should be determined after identifying the impacts of different drainage volumes so as to minimize their impact on the support structure and the stability of the surrounding rocks.

Data Availability
e test data used to support the findings of this study are included within the article.

Conflicts of Interest
No conflicts of interest exit in the submission of this manuscript.