Research on Failure Mechanism and Parameter Sensitivity of Zonal Disintegration in Deep Tunnel

With the increase of excavation depth, the zonal disintegration phenomenon appears in the deep rock mass, which is quite different from the failure mode of shallow tunnel. In order to analyse the failure mechanism of this phenomenon, an elastoplastic softening damage model was put forward based on the softening damage characteristics of deep rock mass. 1e constitutive equations, the equilibrium equations, and the failure criterion were deduced. 1e theoretical solutions of radial displacement and radial stresses and tangential stresses of deep surrounding rock mass were calculated. 1e distribution law of zonal disintegration in deep tunnel was obtained. 1e theoretical solutions presented an oscillating mode. 1e theoretical calculated widths of fracture zones were in good agreement with the in situ test data. Besides, the sensitivity of different parameters to fracture morphology was calculated and analysed. 1e results show that the relative loading strength has a controlling role in the zonal disintegration morphology, followed by the cohesion force and deformationmodulus, and the internal friction angle is the least.1is study reveals themorphological characteristics and influencing factors of zonal disintegration, which provides a basis for the prediction and support control of fracture modes.


Introduction
With the increasing demand for deep underground space utilization and various resource developments, significant nonlinear deformation and failure occurs in the surrounding rock of the tunnel.Many zonal disintegration phenomena [1][2][3] were found during the monitoring of deep rock mass engineering (Figure 1), and many experts and scholars have studied this phenomenon.
Different test methods were used to study the mechanism of zonal disintegration: the cement mortar was used as the model test material, and the zonal disintegration phenomenon was observed [4] (Figure 2).e failure characteristics of the zonal disintegration in the deep surrounding rock were analysed by the geomechanical model tests [5].Many experts have studied the fragmentation phenomenon by numerical simulation [6,7].At present, many theories about the rupture of surrounding rock have been proposed: the occurrence of rupture of surrounding rock is related to the magnitude of the geostress and the strength of the rock [8]; the surrounding rock of the deep roadway will be subject to zonal rupture under the influence of continuous tensile and brittle fracture [9][10][11].
e current research studies mainly focus on the formation mechanism of zonal disintegration, and less attention is paid to fracture morphology caused by intrinsic factors of surrounding rock.
e fracture morphology contains the number and the depth of the fracture zones and will affect the stability of the surrounding rock.So, the sensitivity of intrinsic factors determined to the zonal disintegration morphology should be sorted, and then its influence on the stability of surrounding rock can be analysed.
An elastoplastic model and failure criterion was put forward based on the damage property of rock mass.e morphology of zonal disintegration was obtained through theoretical calculation.
en, the sensitivity of different parameters to the fracture morphology was analysed.is study reveals the morphological characteristics and influencing factors of zonal disintegration, which provides a basis for the prediction and support control of fracture modes.

The Calculation of Elastoplastic Displacement Field and Stress Field in Deep Rock Mass
2.1.Elastoplastic Softening Damage Model.For rock masses, the usual yield criterion is the Mohr-Coulomb yield criterion (Figure 3).According to the yield criterion, the shear strength of the rock mass τ depends on the normal stress σ on the shear slip plane and the shear strength can be expressed as where c is cohesive force, φ is internal friction angle.
According to classical elastic theory, After substituting equation ( 2) into the yield criterion formula and the yield criterion formula can be written as where A � (1 + sin φ)/(1 − sin φ), σ c � 2c cos φ/(1 − sin φ), and the physical meaning of σ c is the uniaxial compressive strength of rock.e surrounding rock under high geostress conditions exhibits certain ductility.After exceeding the peak strength, the rock masses are presented the strain softening characteristics [12], as shown in Figure 4.
Figure 4 shows the stress-strain curves of deep rock masses under high stress conditions, where σ f is peak stress, ε f is peak strain, and ε u is ultimate strain.According to the postpeak strain softening mechanical properties of rock under deep high geostress conditions, the stress-strain relationship is simplified (Figure 5) [13].
In Figure 5, the OA segment is simplified as a linear elastic zone, in which E is the elastic modulus; the AC segment is the plastic bearing zone [14], in which the AB segment is the prepeak plastic hardening segment, in which E 0 is the deformation modulus; and the BC segment is the postpeak plastic damage softening segment, in which |M| is the softening modulus, ε 1 ′ is the deformation in the direction of the principal stress of the postpeak

2
Advances in Civil Engineering plastic softening section, and ε ″ 1 is the irreversible deformation produced by the prepeak plastic hardening section; the CD section is the plastic ow area and σ 0 is the residual stress strength.e total irreversible deformation throughout the plastic deformation process is written as e microcracks begin to propagate near the yield point, and the volumetric strain of rock increases as well (Figure 6).erefore, the damage threshold of rock is at the yield point [15].d is introduced to describe the surrounding rock damage, and the mechanic characteristics of damaged rock under the high geostress are described with the following damage evolution law [16]: e evolution law hypothesizes that there is no damage at the elastic segment, and the rock masses enter damage state when yielding to plasticity.
e relationship between the principal stresses in the plastic softening section is where β is postpeak Poisson ratio, β −ε 3 ′ /ε 1 ′ ; ξ is modulus parameters, and ξ |M/E 0 |.So, the irreversible deformation can be expressed as e strain at the elastic section is e total strain is written as  Advances in Civil Engineering e plastic constitutive formula can be expressed as follows:

Linear elastic section
where

Elastoplastic Stress Field of Surrounding Rock in Circular
Tunnel.Taking the deep circular horizontal tunnel as an example, the stress state of the rocks around the tunnel is analysed.e surrounding rock in the deep tunnel is abstracted as a thick-walled cylinder model.e pressure on the outer wall of the cylinder is P 0 , the radius is b, the excavation radius of the tunnel is a, and the radius of the plastic bearing zone in the surrounding rock is r 12 (Figure 7).
Due to the axial symmetry of the thick-walled cylinder model, it is more convenient to analyse it in the polar coordinate system (r, θ), and the displacement at any position in the surrounding rock is represented by u.Further analysis shows that there are two nonzero items ε r and ε θ in the strain tensor ε ij , so the simplified geometric equation is written as e constitutive equations ( 8) and (10) are expressed in cylindrical coordinates as follows: (i) e elastic zone (ii) e plastic bearing zone (1) Stress solution in plastic bearing zone Assuming that physical strength is neglected, the equivalent equation is written as After substituting equation ( 14) into equation ( 15), the equivalent equation can be expressed as follows:  4 Advances in Civil Engineering where Plastic zone displacement u p can be calculated as where After substituting equation ( 17) into the plastic constitutive equation, the stresses can be expressed as follows: where e inner boundary (r � a) of the plastic bearing zone satisfies the following boundary conditions: (2) Stress solution in elastic zone e equilibrium equation for the elastic zone can also be expressed as After substituting equation ( 13) into equation ( 21), the equivalent equation can be expressed as follows: Elastic zone displacement u e can be calculated as where C 1e and C 2e are undetermined coefficients.After substituting equation ( 23) into the elastic constitutive equation, the stresses can be expressed as follows: e outer boundary (r � b) of the plastic bearing zone satisfies the following boundary conditions: At the elastoplastic boundary (r � r 12 ), the stresses satisfied the yield criterion shown in equation ( 3), so the following boundary conditions are satisfied as After substituting the boundary conditions (20), (25), and (26) into equations ( 18) and (24), respectively, the values of the determined parameters C 1p , C 2p , C 1e , C 2e and the range of the plastic loading zone r 12 can be calculated, and then the elastoplastic stress field of entire circular tunnel can be calculated.

Determination of Fracture Zones and Nonfracture Zones.
As shown in Figure 8, the tangential strain ε θ � 0 and the stress σ θ � σ r in the elastic region away from the tunnel.In the plastic zone near the tunnel, the stresses satisfied equation (26) at the boundary between the plastic bearing zone and elastic zone.
e strength of the rock mass is determined by the bite force, bonding force, and friction of the rock mass.
In the plastic bearing zone, with the increase of deformation, the unevenness on the sliding surface is flattened, and the bite force gradually decreases, so the strength of the rock body gradually decreases to residual strength σ 0 at the boundary r � r 23 of the plastic flow zone.
e residual strength is determined by the bond strength and friction of the rock mass.At this boundary, the stress condition is satisfied as follows: When the elastoplastic stress field of the surrounding rock in the deep circular tunnel is obtained, the boundary conditions equation ( 27) can be used to determine the fracture zone and nonfracture zone of the deep rock mass.Assuming that the stresses of the deep rock mass in the plastic bearing zone satisfied equation ( 27), the rock mass enters the plastic flow zone and ruptured phenomenon occurs, resulting in the first fracture zone with a concomitant decrease in tangential stress and radial stress.e radial Advances in Civil Engineering displacement also increases.e stresses wave continues to propagate to the interior.When the redistribution stresses of the surrounding rock satisfied equation ( 27) again, the stresses are further released and a second fracture zone is obtained.is cycle continues until the redistribution stresses condition cannot meet the boundary.e fracture zone created by each stresses' release is generally a certain distance away from the new excavation boundary, and this zone forms a nonfracture zone.As the location of the fracture zone is far from the inner boundary of the tunnel, the stress field is gradually reduced, the destruction of the surrounding rock is gradually stopped, and the rock body gradually stabilizes, and thus, the phenomenon that the fracture zones and the nonfracture zones appear alternately, shown in Figure 7.
In the range of the plastic bearing zone, substituting equation (18) into equation (12), an expression for the radius r 23 of the plastic flow zone is obtained: rough the above solution, the radius r 23 of the plastic flow zone can be obtained, and then the distribution of the surrounding rock zonal disintegration in the circular tunnel can be obtained.

Selection of Calculation Parameters.
Taking the surrounding rock of tunnel with a near-kilometer level in Huainan Dingji Coal Mine as a zonal disintegration example, the equivalent radius of the tunnel is a � 2.5 m, the radius of the calculated area is b � 20 m, the geostress is P 0 � 34 MPa, the rock density is ρ � 2.70 g/cm 3 , the elastic modulus is E � 77.82 GPa, the deformation modulus is E 0 � 12.97 GPa, cohesion c � 9 MPa , and the internal friction angle is φ � 32 °.

Calculation Result Analysis.
Using the above calculation parameters, the radius of the plastic zone of the surrounding rock is calculated as ρ � 12.63 m; that is, the surrounding rock is in the plastic state within the range of 2.5 m ≤ ρ ≤ 12.63 m, and the surrounding rock is in the elastic state within the range of 12.63 m < ρ ≤ 20 m. Figure 8 is a calculated solution of radial displacement u, radial stress σ r , and tangential stress σ θ of tunnel surrounding rock.
Figure 9 shows that the radial displacement, radial stress, and tangential stress around the tunnel exhibit an oscillation attenuation law in which peaks and troughs appear alternately, and the amplitude of the oscillation decreases gradually with the distance from the tunnel wall.
e theoretical value of the fracture zone was calculated by equation ( 27), as shown in Table 1.Table 1 shows that the number of fracture zone is four and the fracture zone depth is 11.36 m.
e radius of each fracture zone is in good agreement with the in situ measured values [17], confirming the applicability of this calculation model to analyse the zonal disintegration phenomenon.With the help of borehole television, the zonal disintegration is observed which is alternation phenomenon of the fracture zone and nonfracture zone in reference [17].At the same time, the range and radius of the fracture zone were also measured.
From Figure 9 and Table 1, it can be seen that the stresses of the surrounding rock in the nonfracture zone are in the peak area, the stress value is relatively large, the displacement is in the trough area, and the displacement value is relatively small; otherwise, the stresses of the surrounding rock in the fracture zone are in the trough area, and the stress value is relatively small, the displacement is in the peak area, and the displacement value is relatively large, which is consistent with the actual situation.

Analysis of Sensitivity Parameters of Zonal Disintegration Morphology
Zonal disintegration, a special form of deep rock damage, may cause large-scale disaster.Current research studies are focused on exploring its formation mechanism but ignore to explore the production conditions and failure patterns of zonal disintegration and its related influencing factors.e problem is the basic mechanism to predict disintegrative range and take preventive measures such as grouting, anchoring, and so on.e analysis of sensitivity parameters of zonal disintegration morphology has important implications to reveal failure patterns of zonal disintegration and its influencing factors and to build prediction, prevention, and control measures for different evolution mode of zonal disintegration.

Sensitivity Analysis.
Sensitivity analysis is a method for analysing system stability and has been applied in geotechnical engineering [18].With a system, its system characteristics P are mainly determined by several parameters a � (a 1 , a 2 , . . ., a n ).In a certain reference parameter a * � (a * 1 , a * 2 , . . ., a * n ), the reference characteristic is P * .e trend of the system characteristics P deviating from the state 6 Advances in Civil Engineering of the reference characteristic P * is analysed due to the changes in these parameters within their respective possible ranges.is analysis method is called sensitivity analysis.Firstly, a system model should be established, which is a functional relationship between system characteristics and parameters P � f(a 1 , a 2 , . . ., a n ).
is kind of functional relationship, if possible, is expressed analytically, for more complex systems; it can also be represented by calculated methods (such as finite element models, etc.) or graphs.Establishing a system model that is close to the actual system is a crucial task for effectively performing parameter sensitivity analysis.
After setting up the system model, a reference parameter set must be given based on the specific issues to be discussed.For example, in order to analyse the susceptibility of the stability of an underground project to changes in its rock mechanics parameters, the recommended values of the rock mechanics parameters for the project site may be taken as the reference parameter set.After the reference parameter set is determined, sensitivity analysis can be performed on each parameter.When analysing the influence of parameter a k on the characteristics P, the rest of the parameters can be taken as the reference value and fixed, and within the range of possible changes, the system characteristics are expressed as If small changes in a k cause large changes in P, a k is a highly sensitive parameter; however, if large changes in a k cause slight changes in P, a k is a low sensitive parameter.
To compare the sensitivities of different parameters, the sensitivity factor is defined as where S(a k ) is the sensitivity of the parameter a k ; U * is the reference characteristic value corresponding to the reference parameters set; and U a k max and U a k min are the maximum and minimum characteristic values within the variation range of the f parameter a k , respectively.

Reference Parameters.
In order to analyse the sensitivity effects of various parameters on the number of fracture The radial displacement (mm) The radius of tunnel (m) The radius of tunnel (m) The radius of tunnel (m) The radial stress (MPa) The tangential stress (MPa)  e parameters change ranges are selected according to the parameter distribution rules indicated by the test results and shown in Table 2.

e Impact of Various Parameters on the Fracture Zone.
Based on the above-mentioned ranges of basic parameters and theoretical calculation, the effect of each parameter on the number and depth of the fracture zones in the surrounding rock are obtained (Figure 10).
With the increase of surrounding rock mechanical parameters (deformation modulus, cohesion, and internal friction angle), the depth of the fracture zone decreases approximately linearly, and the number of fracture zones is also trapezoidally decreased with the increase of mechanical parameters.
With the increase of the relative loading strength, the depth of the fracture zone firstly increases gently, and when the relative loading strength continues to increase, the trend of the increase in the depth of the fracture zone became acute, and it increases steadily in the later period.It shows that there is a critical value and when the relative loading strength is higher than this value, the depth of the fracture zone will increase rapidly.It can also be seen from the number of fracture zone that when the relative loading strength is small, only one layer fracture zone occurs in the wall of the tunnel; when the relative loading strength is greater than a certain amount, the zonal disintegration phenomenon occurs.

Comparison of Sensitivity of Each Parameter.
In order to further explain the influence of various parameters on the morphology of zonal disintegration, the sensitivity analysis histogram is made according to the sensitivity definition given in equation ( 30) and each parameter within the change range (Figure 11).
Figure 11 shows that the sensitivity of parameters affecting the number and depth of fracture zone is in the order relative loading strength > cohesion force > internal friction angle > deformation modulus.e sensitivity of the relative loading strength is the largest, which shows that the strength of the surrounding rock itself and the state of the geostress are the most sensitive to the zonal disintegration of the tunnel.
Relevant theoretical research studies [19] and the similar model tests [4,5] about zonal disintegration indicate that the relative loading strength is an important condition for the zonal disintegration in the deep rock mass.e relative loading strength sensitivity is significantly higher than the other four factors in Figure 10, indicating that the relative loading strength is the primary condition to determine whether the zonal disintegration occurs.From the calculation results in Figure 9, it can also be seen that only when the relative loading strength is greater than a certain amount, the surrounding rock will have a zonal fracture phenomenon.
During the excavation process of deep tunnels, different excavation methods lead to different unloading times during the adjustment and redistribution of geostress.e sensitivity of the unloading times is second only to the relative loading strength in the formation process and destruction law of zonal disintegration; unloading time plays an important role in surrounding rock of deep tunnel.
e cohesion force, deformation modulus, and internal friction angle of the rock mass are the basic mechanical parameters of the rock mass, and the sensitivity of the internal friction angle is the lowest.erefore, the accuracy of the cohesion force and deformation modulus must be taken seriously under the evaluation of the surrounding rock mechanics parameters.

Conclusion
(1) Based on the softening damage characteristics of deep rock masses, the displacement and stress distributions of tunnel surrounding rock are calculated using the elastoplastic softening damage model.e radial displacement, radial stress, and tangential stress of the surrounding rock alternately oscillate between peaks and troughs.According to the law of change, the width and number of fracture zones in the theoretical calculations are in good agreement with the in situ measured values, confirming the applicability of the model in analysing the zonal disintegration phenomenon.For deep rock mass under high geostress conditions, due to the oscillation-type changing rules of radial displacement and stress field in the surrounding rock of deep tunnel, the fracture zone and the nonfracture zone alternately appear, resulting in the occurrence of zonal disintegration.(2) e parameter sensitivity affecting the zonal disintegration morphology of surrounding rock in deep tunnel is in the order relative loading strength > cohesion force > internal friction angle > deformation modulus.e relative loading strength is the decisive condition for the zonal disintegration phenomenon.e ranking of the relevant mechanical parameters is of great significance for the prediction, prevention, and remediation measures of different zonal disintegration evolution patterns.
(3) is theoretical calculation method has a broader application prospect.en, it also indicates the   Advances in Civil Engineering direction of effort to improve the calculation accuracy of the calculated model through the parameter sensitivity analysis and provides a basis for the prediction and support control of different zonal disintegration modes.

12 Figure 8 :
Figure 8: Zonal disintegration of surrounding rock of tunnel after excavation.

Figure 9 :
Figure 9: Calculated solutions of displacement and stress of surrounding rock in tunnel.

Table 1 :
Comparison between theoretical calculation result and in situ test value.depth of fracture zones, the mechanical parameters of tunnel surrounding rock in Dingji Coal Mine are used as reference parameter values.e values are as follows: deformation modulus E 0 � 12.97 GPa, cohesion c � 9 MPa, internal friction angle φ � 32 ∘ , and relative loading strength κ � P z /σ c � 2, where κ is the ratio of tunnel axial stress to compressive strength of surrounding rock which is a parameter can reflect the geostress state.

Table 2 :
Basic values and variations of parameters.