Upper Bound Solution of Surrounding Rock Pressure for Deep Cavity Using Nonlinear Hoek – Brown Failure Criterion

(e linear failure criterion is generally adopted in the stability problem of geotechnical engineering, whereas the experiments have indicated that there is a nonlinear relationship between the maximum and minimum principal stresses in weak surrounding rock. According to the characteristics of weak rock, the failure mechanism of deep cavity is constructed by combining the nonlinear Hoek–Brown failure criterion and the upper bound theorem of limit analysis. (e upper bound solution of the surrounding rock pressure was deduced using the tangent method. (e results show that the surrounding rock pressure of deep cavity is affected by surrounding rock grade, cavity depth, and section size. Especially, the influence of the disturbance factor is quite obvious. Upper bound solution based on nonlinear Hoek–Brown failure criterion can fully take the influence of parameters on surrounding rock pressure. (erefore, this method is more scientific than the linear failure criterion to calculate the surrounding rock pressure.


Introduction
e stability problem of deep cavity has been highlighted greatly and continuously in the field of the underground engineering, whereas the problem that needs to be solved urgently is how to accurately obtain the surrounding rock pressure when the deep cavern is destroyed [1,2].Currently, the analytic methods for solving the surrounding rock pressure mainly include the limit equilibrium method and the limit analysis method, both of which have been illustrated widely in literature.However, the former does not consider the constitutive relationship of rock, while the latter does by adopting the orthogonal flow rule.erefore, the upper bound theory is more rigorous than the limit equilibrium method.Especially, there is no need to study the whole process of rock's elastic-plastic deformation, but directly to pay attention to the ultimate failure state of rock.us, its calculation process is simple and has been regarded as one of the most effective approaches to solve the surrounding rock pressure problem in the deep cavity [3][4][5].
In recent years, some scholars have adopted the limit analysis method to study the surrounding rock pressure of deep cavity.According to the results of the model test, Davis et al. [6] proposed four failure mechanisms for underground cavity, which included the arch collapse mode and vault collapse mode bordered wall, and obtained the upper bound solution of the surrounding rock pressure by the limit analysis method.Takemura et al. [7] constructed the failure mechanism of the deep cavity using the centrifugal model test.e failure model is composed of 5 rigid blocks, and the limit analysis method is used to solve the surrounding rock pressure of the underground cavity under the anisotropic soil conditions.e results were compared with the measured values, and the validity of the results has been verified.Yang and Yang [8] proposed the destruction mechanism of underground chambers composed of n blocks.According to the requirement of calculation accuracy, the upper bound solution of surrounding rock pressure under different number of blocks was obtained.rough comparison, the destruction characteristics of deep cavity can be better illustrated with more rigid blocks, and accordingly, the surrounding rock pressure is more accurate.However, the above research studies are applied to shallow buried cavity where the Mohr-Coulomb failure criterion is used and more suitable for soil mass.
For deep cavity, Atkinson and Potts [9] established the failure mode of circular which were based on the results of model tests, obtained the upper bound solution of the surrounding rock pressure using the limit analysis method, and the correctness of the results was verified.And the curve failure mechanism of the deep cavity based on the limit analysis method was proposed by Fraldi and Guarracino [10][11][12], which the Hoek-Brown failure criterion and the variational method were applied to solve the analytic expression of the curve and got the surrounding rock pressure under the limit state of the vault collapse.Yang and Yang [8] combined the limit analysis theory with the variational principle, deduced the expression of the surrounding rock pressure of deep rectangle cavity and specified the corresponding failure range.In addition, Huang et al. [13] took the influence of pore water based on the previous research studies into account and got a reasonable result.e limit analysis method was used to calculate the surrounding rock pressure of deep buried circular cavity.Compared with the existing research results, the correctness of the calculation results was verified.Qin et al. [2] considered the multilayer rock mass and used the upper bound method to study the influence of the variation of the parameters of the Hoek-Brown failure criterion on the surrounding rock pressure of the deep buried rectangular cavity.e vault of deep cavity was only considered and assumed to obey a curve function f(x) in the abovementioned studies, and yet, the side wall was ignored in the analysis.However, it could not accord with the conditions in reality.erefore, by referring to the previous research results, the multiblock and logarithmic spiral curves are used to construct the failure mechanism of deep cavity in this paper.Based on the upper bound theorem of limit analysis and Hoek-Brown failure criterion, the vault and side wall damage state is investigated to determine the rock pressure when the deep cavity collapses and provide the scientific basis for the design in engineering.

Hoek-Brown Nonlinear Failure Criterion and Tangent Method
In 1980, Brown proposed the initial Hoek-Brown failure criterion on the basis of Hoek's experimental research studies on rock plastic behaviors.Subsequently, Hoek made a further revision to the failure criterion and put forward the modified generalized Hoek-Brown criterion.It is expressed as follows [14,15]: where σ 1 ′ and σ 3 ′ are the maximum and minimum effective principal stress of rock mass at critical failure, respectively; σ ci is the uniaxial compressive strength of rock; α is the physical index relating to the integrity of rock; m b and S indicate the dimensionless parameters related to geotechnical materials, which characterizes the integrity of the rock.e above parameters can be calculated through the GSI index; the specific process is as follows [16][17][18]: + e −GSI/15 − e −20/3  . (2) For the upper limit solution of the deep cavity under the Hoek-Brown failure criterion, it is necessary to use the tangent method to obtain the counterforce of the support.Among them, the correlations between c t and φ t can be written as follows [19,20]: In the above formula, the nonlinear shear strength index φ t can be optimized by the least energy principle when the upper bound is calculated.e c t is obtained by Equation (3) after the determination of φ t .

Building Destruction Mechanism
During the excavation process of the deep cavity, the collapse easily happens due to the untimely support or weak strength.Based on the existing research results [8,10] and the requirement of the upper bound theorem of limit analysis, the failure mode of deep cavity was structured, as shown in Figure 1. e width of the chamber is l, and the height is h.e wedge ABGOG 1 B 1 at the top of the cavity is collapsing vertically and downwards at a rate of v 0 .e BGC and B 1 G 1 C 1 segments, respectively, rotate failure around G and G 1 points.
e triangles CGD, DGE, and EGF and and E 1 G 1 F 1 would undergo translational failure.e surrounding rock pressure of vault and side wall are q and e, respectively, and e � Kq.
e velocity vector diagram corresponding to this failure mechanism is shown in Figure 2.
Calculation of speed and velocity discontinuity's length: Based on the analysis of the velocity field in the shear zone of the circular arc radiation, the following can be obtained: From Figure 2 and geometric triangle relationship can be obtained: 2 Advances in Civil Engineering e recurrence relation of various velocity discontinuities is as follows:

Calculation of Surrounding Rock Pressure
According to the upper bound theorem of limit analysis, the external power and internal energy dissipation rate can be obtained when the deep cavity collapses.If the external power and internal energy dissipation rate are equal, the surrounding rock pressure at the collapse of the deep cavity can be obtained.

Power of Gravity.
Area of each rigid block: e gravitational power generated by each rigid slider is as follows.Quadrilateral rigid block ABGO: Sector region GBC: e power calculation of circular shear zone BGC is shown in Figure 3.
Among them, Triangular rigid block GCD: Among them, Triangular rigid block GDE: Among them, Triangular rigid block GEF: Among them, Power of gravity: 4.1.2.Power of Supporting Force.
Among them, e � Kq, e total external power is equal to the sum of the deadweight power of the surrounding rock and the power of the supporting counterforce.at is, Among them, Energy dissipation of circular shear surface BC and shear zone GBC: Energy dissipation along discontinuity line CD: Energy dissipation along discontinuity line GD: Energy dissipation along discontinuity line DE: Energy dissipation along discontinuity line GE: Energy dissipation along discontinuity line EF: Total internal energy dissipation: (30)

Calculation of Supporting Counterforce.
According to the principle of virtual power, e expression of surrounding rock pressure is Above all of the formulas, h is the height of the cavity, measured in meter.l represents the width of the cavity, measured in meter.q is the vertical supporting pressure, and e is the lateral supporting pressures of which the unit is kPa.φ is the internal friction angle, measured in degree.c is the unit weight of soil, measured in kg/m 3 .c is the cohesion, measured in kPa.K is the lateral pressure coefficient.V 0 , V 1 , V 2 , V 3 , and V 4 are the speed in the velocity field.α 1 , α 2 , α 3 , and α 4 are the geometrical variables which determine the shape of the failure mechanism whose unit is degree.W is the external power, W exit is the external total power, W soil is the weight power of surrounding rock, W T is the power of support of the antiforce, and the unit is watt.W 1 , W 2 , W 3 , W 4 , and W 5 are variables.D is the dissipation power of internal energy, and f 0 , f 1 , f 2 , . .., f 14 are variables.
e constraint condition of Equation ( 32) is Invoking the command in the Matlab software optimization toolbox, it will generate a set of angle data under the constraint condition (33).According to the formula of the surrounding rock pressure, an upper bound solution q is obtained; then, by adjusting the variable parameter values, numerous q values are obtained, in which the maximum q is the upper bound solution of the surrounding rock pressure.

Result Analysis
Under the nonlinear Hoek-Brown failure criterion, the upper bound theorem is used to analyze the influence of related parameters on the surrounding rock pressure of deep cavity.According to the failure mechanism and method mentioned above, the calculation results are as follows.

Influence of c and K on Surrounding Rock Pressure.
To investigate the influence of c and K on the surrounding rock pressure, the parameters are assumed as follows: Cavity width l � 10 m, height h � 10 m, geological strength index GSI � 20, rock mass constant m i � 20, uniaxial compressive strength σ ci � 400 kPa, disturbance factor D � 0, and c is the soil gravity, which had taken the values of 16, 18, 20, 22, 24, and 26 kN/m 3 , respectively, in the experiment.Similarly, the coefficient K was taken 0.4, 0.6, 0.8, 1.0, 1.2, and 1.4, respectively.According to the failure mechanism constructed in this paper which combined translational destruction and rotational destruction together, the upper bound solution of the surrounding rock pressure obtained by the limit analysis method is shown in Tables 1 and 2 and Figures 4 and 5. e tables and figures present that the surrounding rock pressures q and e increase with the rising of soil gravity c, and the pressure of surrounding rock q decreases, the surrounding rock pressure e increases with the rising of coefficient K.It shows that the buried depth has certain influence on the surrounding rock pressure.
e graph is drawn from the above table, as follows.

Influence of Chambers Section Size on Surrounding Rock
Pressure.Likewise, to study the influence of chambers section size on the surrounding rock pressure, the parameters are soil gravity c � 20 kN/m 3 , coefficient K � 1.0, geological strength index GSI � 20, rock mass constant m i � 20, uniaxial compressive strength σ ci � 400 kPa, disturbance factor D � 0; the value of chamber width l and height h was taken as 5, 6, 7, 8, 9, and 10 m, respectively, in the experiment to conduct further analysis.According to the failure mechanism constructed in this paper which combined translational destruction and rotational destruction together, the upper bound solution of the surrounding rock pressure obtained by the limit analysis method is shown in Tables 3 and 4 and Figures 6 and 7. e tables and figures present that the surrounding rock pressures q and e increase with the rising of the chamber width l and height h, and the effect is obvious.Based on the research and above illustration, it is obvious that the chambers section size has a great influence on the pressure of surrounding rock.

Influence of GSI and m i on Surrounding Rock Pressure.
To study the influence of GSI and m i on the surrounding rock pressure, the parameters are soil gravity c � 20 kN/m 3 , coefficient K � 1.0, cavity width l � 10 m, height h � 10 m, uniaxial compressive strength σ ci � 400 kPa, disturbance factor D � 0, and the value of geological strength index GSI was taken as 10, 15, 20, 25, 30, and 35, respectively; rock mass constant value m i was taken as 5, 10, 15, 20, 25, and 30, respectively.According to the failure mechanism constructed in this paper which combined translational destruction and rotational destruction, the upper bound solution of the surrounding rock pressure obtained by the limit analysis method is shown in Tables 5 and 6 and Figures 8 and 9. e graph is drawn from the above table, as follows.When coefficient K � 1.0, the surrounding rock pressure e � q. e tables and figures show that the surrounding rock pressure q and e decrease with the increasing of the geological strength index GSI and rock mass constant m i , and the influence is apparent.It could be concluded that the

6
Advances in Civil Engineering quality of surrounding would have a significant and direct influence on the pressure of surrounding rock.

Influence of σ ci and D on Surrounding Rock Pressure.
To study the influence of σ ci and D on the surrounding rock pressure, the parameters are soil gravity c � 20 kN/m 3 , coefficient K � 1.0, cavity width l � 10 m, height h � 10 m, geological strength index GSI � 20, rock mass constant m i � 20, uniaxial compressive strength had taken the value of 200, 400, 600, 800, 1000, and 1200 kPa, and the disturbance factor D had taken the value of 0, 0.2, 0.4, 0.6, 0.8, and 1.0, respectively, to conduct further analysis.According to the failure mechanism constructed in this paper which combined translational destruction and rotational      Advances in Civil Engineering destruction, the upper bound solution of the surrounding rock pressure obtained by the limit analysis method is shown in Tables 7 and 8 and Figures 10 and 11.
e graph is drawn from the above table, as follows.When coefficient K � 1.0, the surrounding rock pressure e � q. e tables and figures show that the surrounding rock pressure q and e decrease with the increasing of the uniaxial compressive strength σ ci , while the surrounding rock pressure q and e increase with the rising of disturbance factor D. And the effect is obvious, which reveals that the disturbance factor would influence the pressure of surrounding rock directly.

Conclusion
(1) e nonlinear Hoek-Brown failure criterion is applied to the upper bound of limit analysis by the       Advances in Civil Engineering tangent method to construct the circular arc failure mechanism of deep cavity.According to the principle of virtual power, the analytic expression of surrounding rock pressure could be deduced, and the upper bound solution of the surrounding rock pressure of deep cavity could be obtained.(2) e surrounding rock pressures increase obviously with the rising of soil gravity c, tunnel width l, and tunnel height h.Besides, the pressure of surrounding rock q decreases with the rising of coefficient K, whereas the change rule of the surrounding rock pressure e is quite the opposite.(3) e surrounding rock pressures decrease significantly with the increase of GSI, m i and σ ci , while the surrounding rock pressure increase greatly with the rising of D.

Figure 3 :
Figure 3: Schematic diagram of the power calculation of circular shear zone BGC.

2 .
Internal Energy Dissipation Power.Energy dissipation along discontinuity line AB:

Figure 4 :
Figure 4: Effect of c and K on the surrounding rock pressure q.

Figure 7 :Figure 5 :
Figure 7: Effect of l and h on the surrounding rock pressure e.

Figure 6 :
Figure 6: Effect of l and h on the surrounding rock pressure q.

Figure 10 :
Figure 10: Effect of σ ci and D on the surrounding rock pressure q.

Figure 8 :
Figure 8: Effect of GSI and m i on the surrounding rock pressure q.

Figure 9 :
Figure 9: Effect of GSI and m i on the surrounding rock pressure e.

Figure 11 :
Figure 11: Effect of σ ci and D on surrounding rock pressure e.

Table 2 :
Upper bound solution of the surrounding rock pressure e at different c and K.

TABLE 1 :
Upper bound solution of the surrounding rock pressure q at different c and K.

Table 3 :
Upper bound solution of the surrounding rock pressure q at different l and h.

Table 4 :
Upper bound solution of the surrounding rock pressure e at different l and h.

Table 5 :
Upper bound solution of the surrounding rock pressure q at different GSI and m i .

Table 7 :
Upper bound solution of the surrounding rock pressure q at different σ ci and D.

Table 8 :
Upper bound solution of the surrounding rock pressure e at different σ ci and D.

Table 6 :
Upper bound solution of the surrounding rock pressure e at different GSI and m i .GSI e (kPa)m i � 5 m i � 10 m i � 15 m i � 20 m i � 25 m i � 30