Elastoplastic Analysis for Circular Tunnel Based on Modified Lade Criterion considering Strain Softening and Dilatancy

Deng, Yansheng; Wang, Jianxiu; Zou, Baoping

doi:https://doi.org/10.1155/2021/4872558

Advances in Civil Engineering

On this page

Abstract Introduction Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 4872558 | https://doi.org/10.1155/2021/4872558

Elastoplastic Analysis for Circular Tunnel Based on Modified Lade Criterion considering Strain Softening and Dilatancy

Yansheng Deng,^1,2Jianxiu Wang,^2,3,4and Baoping Zou¹

Academic Editor: Yang Wu

Received23 Jun 2021

Revised20 Sept 2021

Accepted12 Oct 2021

Published25 Oct 2021

Abstract

Modified Lade criterion can not only describe the strength properties of many kinds of rocks well but also has simple and practical parameters. Although the elastoplastic solution of circular tunnel has been extensively investigated, the method based on modified Lade criterion considering the effect of the intermediate principal stress, strain-softening behavior, and dilatancy has not yet been studied. In this paper, a new numerical procedure based on modified Lade criterion is proposed to calculate the elastoplastic solutions for surrounding rock of the circular tunnel. The comparisons of stress, displacement, and plastic zone radius are carried out between the presented method and published literatures under axisymmetric and nonaxisymmetric original in situ stress field. Finally, a series of parametric analyses are executed and discussed. It can be concluded that the lateral pressure coefficient, λ, influences both the size of plastic zone and the development direction. The plastic zone radius shows a negative power function change with increasing critical deviatoric plastic strain and increases slightly with increasing dilation angle, ψ.

1. Introduction

Elastoplastic analysis of circular tunnels is one of the most fundamental problems in the fields of tunneling, mining, and geotechnical engineering. After tunnel excavation, the original state of stress equilibrium is disrupted, which causes redistribution of the stress of surrounding rock, elastoplastic deformation, or even failure at the same time. To ensure the safety of tunnel excavation, the stress and displacement solutions around a tunnel are required in the design stage of a tunnel. Based on Mohr–Coulomb (M-C) strength criterion, Fenner [1] first derived the analytical solution of circular tunnel surrounding rock, which was then modified by Kastner considering the nonhydrostatic stress field in 1949 [2]. For a long time, M-C and Hoek–Brown (H-B) failure criteria have been employed to carry out the elastoplastic analysis for the circular tunnel. Ogawa and Lo [3] investigated the influences of dilatancy and failure criteria (M-C and H-B criteria) on stress and displacement for a circular tunnel. Brown et al. [4] presented two new elastoplastic solutions for circular axisymmetric tunnel problem with H-B criterion. Sharan [5–7] took various rock mass behaviors into account to derive elastoplastic solutions with H-B and generalized H-B criterion. Alonso et al. [8] presented a one-dimensional numerical solution for circular tunnels with strain-softening surrounding rock. Lee and Pietruszczak [9] proposed a new numerical method accommodating both the M-C and H-B criteria to execute the elastoplastic analysis for circular opening with strain-softening rock mass. Based on the idea of the simplification of strain-softening process, Wang et al. [10] treated the strain-softening process as a series of stress drop and plastic flow with increasing plastic strain and proposed a new approach to calculate the elastoplastic solutions with M-C and H-B criteria. According to similar ideas, Zhang et al. [11] assumed the surrounding rock properties were uniform in a pretty small region and presented a multistep brittle-plastic model. Cui et al. [12] introduced a variable critical plastic softening parameter to carry the elastoplastic analysis of a circular opening in strain-softening rock mass.

Although various exact solutions or numerical solutions have been proposed in the above works, the influence of intermediate principal stress was ignored in the strength criteria, which has significant effect on the mechanical behaviors of surrounding rock [13–15]. Therefore, the yield criteria considering intermediate principal stress were gradually adopted to investigate the elastoplastic solutions of tunnel surrounding rock and achieved significant results, such as Drucker–Prager (D-P) criterion [16, 17], the unified strength theory [18–22], and Lade–Duncan (L-D) criterion [23]. However, the assumption of axisymmetric in situ stress field is still introduced to calculate the elastoplastic solution of rock mass, which is not in line with the actual stress state under most circumstances [24]. In addition, the parameters of these strength criteria are difficult to be obtained. Modified Lade criterion (MLC) was proposed by Ewy [25] based on original Lade’s criterion (LC) [26, 27], which is a simplified version of LC and is compelled to agree with M-C strength criterion. The MLC can not only describe the impact of intermediate principal stress on mechanical properties of surrounding rock but also have simple parameters. Its parameters S_a and η can be calculated as a function of the well-known c and φ from the M-C criterion. In addition, based on the comparison result of various strength criteria for different types of rock limited by polyaxial experimental data, the MLC shows the better agreement with the test data [28, 29]. Therefore, in this study, a new numerical procedure based on the MLC is proposed to carry out the elastoplastic analysis of surrounding rock under nonaxisymmetric original stress filed.

2. Problem Description

2.1. Basic Assumption

Due to the effects of geological structure and depth, it is unreasonable to assume that the in situ stress of tunnel surrounding rock is an axisymmetric stress field. Therefore, the elastoplastic analysis for the circular tunnel is based on the following assumptions:(1)The initial vertical in situ stress is P₀, and the horizontal in situ stress is λP₀, where λ is the lateral pressure coefficient(2)The circular tunnel is deep-buried and infinite, and the plane strain assumption is satisfied(3)Rock mass shows initial elasticity, isotropy, and strain-softening behavior obeying MLC

2.2. Solutions in the Elastic Zone

According to the basic assumption, the principle of stress superposition of structural mechanics can be employed to carry the mechanical analysis in the elastic zone [30, 31]. Therefore, the original nonaxisymmetric in situ stress field (Model A) can be decomposed into the sum of uniform in situ stress field (Model B) with λP₀ and vertical stress field (Model C) with (1-λ)P₀ (Figure 1).

2.2.1. Stress Solutions

The radial and tangential stress (σ_Br and σ_Bθ) in the elastic zone of Model B with uniform stress λP₀ are given aswhere P_i is the support pressure, λ is lateral pressure coefficient, R₀ is the radius of the circular tunnel, and r is radial distance from any point to the center of the tunnel.

The radial, tangential, and shear stress (σ_Cr, σ_Cθ, and τ_Crθ) in the elastic zone of Model C with vertical stress (1-λ)P₀ are given aswhere θ represents the angle from the horizontal direction.

From equations (1) and (2), the total stress components in the elastic zone of the Model A can be obtained as

2.2.2. Displacement Solutions

According to Hooke’s law, the radial displacement (u) of any point in the surrounding rock can be expressed aswhereand E is Young’s modulus of surrounding rock and is Poisson’s ratio.

Substituting equation (5) into equation (4), there is

2.3. Solutions in the Plastic Zone

2.3.1. Description of the MLC

Due to no limitation of the experimental data and the assumption of failure envelope shape, the LC can reflect well the strength properties of most of the rocks tested under various stress conditions [32]. However, there is no direct relationship between the parameters from LC and c and φ from the most commonly used M-C criterion, which limits the application of LC at some extent. For the purpose of more applicable, the parameters of MLC are determined directly by c and φ from M-C criterion. Therefore, the MLC not only can take the effect of intermediate principal stress into account but also is easier to be applied in engineering, which is expressed aswhere and are the appropriate first and third stress tensor invariants, σ₁, σ₂, and σ₃ are the principal stresses, S_a and η are material constants, which can be determined by the cohesion, c, and internal friction angle, φ, from the M-C criterion, and P_p represents the pore pressure.

In order to take the dilatancy into account at the plastic stage, the relationship of the principal stress and dilation angle is employed as follows [33]:where ψ is the dilation angle.

Substituting equation (9) into equation (7), the MLC can be reformulated as

2.3.2. Plastic Potential Function and Softening Parameter

The plastic potential function can be written aswhere γ_p represents the deviatoric plastic strain and k(γ_p) indicates the coefficient of dilation.

According to the plastic flow law, the relation between radial and tangential plastic strain increments can be expressed as

Considering the strength parameters as a function of deviatoric plastic strain γ_p, it is reasonable to assume the bilinear functions of γ_p can reflect the change of strength parameters aswhere ω indicates one of the c, φ, and ψ, ω^p and ω^r are the peak value and residual value of any one of the strength parameters (c, φ, and ψ), and represents the critical deviatoric plastic strain.

2.3.3. Stress and Displacement Solution in the Plastic Zone

The stress in the plastic zone is related to the properties of the tunnel and surrounding rock, not to the original in situ stress field. When the stress in the plastic zone of circular tunnel is in equilibrium, it should satisfy the axisymmetric static equilibrium equation as

Based on equation (3), the stress components of the boundary between elastic and plastic zone (r = R_p) can be given aswhere R_p is the radius of plastic area and σ_Rp represents the radial stress of the elastic-plastic boundary.

Substituting equation (16) into equation (10), there iswhere

According to equation (17), σ_Rp can be obtained numerically by using appropriate algorithm such as Cardano formula. Then, σ_θb can be calculated by equation (16).

As can be seen from equation (3), the total stress components (σ_Ar, σ_Aθ, τ_Aθr) in the elastic zone is function of θ, which means that the stress components (σ_rb, σ_Rp, and σ_θb) of the boundary between elastic and plastic zone (r = R_p) and the radius of plastic area (R_p) are also function of θ and is not a constant. Actually, equations (16) and (17) can also explain this. Therefore, in order to calculate the stress and displacement in the plastic zone, the plastic area can be divided into n (n ≥ 500) concentric annuli with the same radial stress under different θ and σ_rb, as shown in Figure 2 [9].

The normalized radius, ρ_i, is used to represent the inner radius of i^th annulus and can be expressed as

The normalized radii of the elastic-plastic boundary can be depicted by ρ₀, and there is

Therefore, the radial stress increment of i^th annulus can be given by

By substituting equation (10) into equation (15), there iswhere

Based on equation (15), the tangential stress σ_θ(i) of i^th annulus can be obtained by

To the plane strain problem, the elastic strain increment of i^th annulus can be written by Hooke’s law as follows:

The total strains in the plastic zone can be resolved into elastic and plastic strain, and there iswhile the compatibility equation must be satisfied as

Based on equations (26) and (27), the compatibility equation can be rewritten as

Combining with equations (13), (25), (24), and (28), can be obtained by the approximating differential formation, and there iswhere

Based on equation (26), the total strain of i^th annulus in the plastic zone can be expressed as

So, the radial displacement u_(i) of i^th annulus can be calculated from the normalized radial displacement, and there is

Because R_p is not a constant, different R_p can be obtained based on different θ with equations (19)–(30). For example, when θ = 0°(marked as θ_|θ ₌ _0°), R_p can be calculated by the proposed method and marked as R_p|θ ₌ _0°. In order to calculate R_p_|θ ₌ _0°, the plastic area can be divided into n (n ≥ 500) concentric annuli with the same radial stress (Δσ_r|θ=0°), and then, R_p_|θ ₌ _0° can be obtained by equations (19)–(30), as shown in Figure 3(a). Then, R_p|θ ₌ _1°, R_p|θ ₌ _2°, …, R_p|θ ₌ _90° can be calculated with the same method. The plastic zones were plotted in Figures 3(b)–3(d) to make it easy to understand when θ = 20°, θ = 40°, and θ = 90°. The flowchart for the sequence of calculations is shown in Figure 4.

(a)

(b)

(c)

(d)

3. Verification Examples

3.1. Verification under Axisymmetric Stress Field

In order to verify the proposed numerical procedure based on MLC, σ₂ = σ₃ and ψ = −90° in equation (9) are assumed to ignore the effect of intermediate stress principal, and then, the results from this study with λ = 1 can be compared with the M-C criterion. The dataset was taken from [9] as input data: R₀ = 5.0 m, P₀ = 3 MPa, P_i = 0 MPa, E = 10 GPa, , φ_p = 30°, φ_r = 26°, c_p = 0.5 MPa, and c_r = 0.2 MPa.

Figure 5 presents the comparisons of stress and displacement from degenerated MLC and M-C criterion. It can be seen that radial stress, tangential stress, and radial displacement calculated from the proposed method are in good agreement with the published literature.

(a)

(b)

3.2. Verification under Nonaxisymmetric Stress Field

Because the original in situ stress field is not ideal uniform stress field, it is significant to consider the effects of lateral pressure coefficient, λ. To verify the results under nonaxisymmetric stress field, the parameters used from [31] are R₀ = 2.965 m, P₀ = 10 MPa, P_i = 0 MPa, E = 20 GPa, , φ_p = 30°, φ_r = 30°, c_p = 0.3 MPa, c_r = 0.3 MPa, and λ = 1.0, 1.5, and 2.0, respectively. σ₂ = σ₃ and ψ = −90° are still used in equation (9) to ignore the effect of intermediate principal.

The results obtained by this study and published literatures are shown in Table 1. From Table 1, both the results from this study are slightly smaller than [31, 34] under different λ. When λ = 1.5 and 2.0, the results of this study with θ = 90° are a litter bigger than [35]. In general, the results from proposed numerical procedure are basically consistent with the literatures.

4. Parametric Analysis and Discussion

4.1. Influence of Lateral Pressure Coefficient, λ

The lateral pressure coefficient, λ, is used to represent the inhomogeneity of original in situ stress field. λ < 1 means the vertical stress is greater than the horizontal stress, while λ > 1 is on the contrary. λ = 1 indicates the original in situ stress is uniform stress field. In order to investigate the effect of λ on the radius of plastic zone, the datasets from Section 3.2 was taken as input data with λ = 0.25, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0, 2.5, and 3.0, respectively. The plastic zone and relationship between λ and R_mp/R₀ are shown in Figure 6, where R_mp represents the maximum plastic radius.

From Figure 6, when λ < 1, R_mp/R₀ decreases gradually with the increase of λ. The plastic zone is mainly located on both sides of the circular tunnel with λ = 0.25. With increasing λ, the shape of plastic zone is sliding into an ellipse with the major axis in the horizontal direction. When λ = 1.0 with the uniform stress field, a circular plastic zone is developed. When λ > 1, R_mp/R₀ presents positive logarithm relevant to λ. The elliptical plastic zone appears again, while the major axis turns into the vertical direction. The plastic zone develops mainly along the top and bottom when λ = 3.0. Therefore, λ affects not only the size of plastic zone but also the development direction, which should be reasonably taken into account in design and construction of the tunnel.

4.2. Influence of Critical Deviatoric Plastic Strain,

represents the start point of residual behavior of surrounding rock, which can indicate the softening degree to a certain extent. The surrounding rock enters the residual state quickly with no softening behavior when is very small. To study the influence of on the radius of the plastic zone, the datasets from Section 3.2 was taken as input data with λ = 1.0, and , respectively. Relationship between R_p/R₀ and is plotted in Figure 7.

Figure 7 shows that, with the increase of , R_p/R₀ decreases by means of a negative power function. Especially, R_p/R₀ drops from 8.73 to 3.2 rapidly when raises from 0.005 to 0.05. Then, with the increase of continuously, the rate of decrease gradually reduces. In this way, it is pretty significant to determine the correct for the calculation of the plastic zone.

4.3. Influence of the Dilation Angle, ψ

To study the influences of the dilation angle, ψ, on the plastic zone, the datasets from Section 3.2 was still taken as input data with λ = 1.0, and ψ = 5°, 10°, 15°, 20°, 25°, and 30°, respectively. Relationship between R_p/R₀, radial displacement, and different ψ is shown in Figure 8.

(a)

(b)

From Figure 8(a), with the increase of ψ, R_p/R₀ is approximately linear rise, whereas the slope of regression equation is 0.004, which means the gradient of increase is very small. When ψ increases from 5° to 30°, R_p/R₀ is just 4.9% more. In this way, ψ has no significant effect on the radius of plastic zone. Figure 8(b) plots the relationship between radial displacement and ψ, where the radial displacement increases exponentially with increasing ψ. R_p/R₀ increases more than 3 times when ψ raises from 5° to 30°. That means the dilation angle has an important role to affect the radial displacement, which are coincident with [3]. Hence, the correctness of the presented method is also verified to a certain extent.

5. Conclusions

Based on the relationship between the principal stresses and ψ, the MLC was first introduced to consider the effects of the intermediate principal stress and dilatancy on the elastoplastic behaviors of surrounding rock under the nonaxisymmetric stress field. A new numerical procedure was then proposed and verified. At last, a series of parametric analyses were performed. The following conclusions can be obtained:(1)Due to the convenience of parameters, the presented method can be efficiently employed to comprehensively carry out the elastoplastic analysis for the circular tunnel under both axisymmetric and nonaxisymmetric stress field.(2)The λ influences not only the size of plastic zone but also the developed direction. When λ < 1, the plastic zone grows along the horizontal direction, while develops vertically when λ < 1. Therefore, it is pretty significant to consider the inhomogeneity of the original in situ stress field.(3)It is important to determine accurately the value of based on experiments, which has great influence on the plastic zone radius, especially when it is smaller. In addition, the effect of ψ on the plastic zone radius is much greater than that on radial displacement.

Data Availability

The data used to verify the proposed numerical procedure are on pages 594–597 of [9] (https://www.sciencedirect.com/science/article/abs/pii/S0886779807001125).

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work was sponsored by National Key R&D Program of China (2017YFC0806000) and Zhejiang University of Science and Technology (Project no. 2021QN032).

References

R. Fenner, “Untersuchungen zur erkenntnis des Gebirgsdruckes,” Gluckauf, vol. 74, no. 681–695, pp. 705–715, 1938.
View at: Google Scholar
H. Kastner, “Uber den Echten Gebiergsdruck beim bau Tiefliengender, Tunnel,” Osterreich Bauzeitschrift, Vienna, Austria, vol. 10, 1949.
View at: Google Scholar
T. Ogawa and K. Y. Lo, “Effects of dilatancy and yield criteria on displacements around tunnels,” Canadian Geotechnical Journal, vol. 24, no. 1, pp. 100–113, 1987.
View at: Publisher Site | Google Scholar
E. T. Brown, J. W. Bray, B. Ladanyi, and E. Hoek, “Ground response curves for rock tunnels,” Journal of Geotechnical Engineering, vol. 109, no. 1, pp. 15–39, 1983.
View at: Publisher Site | Google Scholar
S. K. Sharan, “Elastic-brittle-plastic analysis of circular openings in Hoek-Brown media,” International Journal of Rock Mechanics and Mining Sciences, vol. 40, no. 6, pp. 817–824, 2003.
View at: Publisher Site | Google Scholar
S. K. Sharan, “Exact and approximate solutions for displacements around circular openings in elastic-brittle-plastic Hoek-Brown rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 42, no. 4, pp. 542–549, 2005.
View at: Publisher Site | Google Scholar
S. K. Sharan, “Analytical solutions for stresses and displacements around a circular opening in a generalized Hoek-Brown rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 1, pp. 78–85, 2008.
View at: Publisher Site | Google Scholar
E. Alonso, L. R. Alejano, F. Varas, G. Fdez-Manin, and C. Carranza-Torres, “Ground response curves for rock masses exhibiting strain-softening behaviour,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 27, no. 13, pp. 1153–1185, 2003.
View at: Publisher Site | Google Scholar
Y.-K. Lee and S. Pietruszczak, “A new numerical procedure for elasto-plastic analysis of a circular opening excavated in a strain-softening rock mass,” Tunnelling and Underground Space Technology, vol. 23, no. 5, pp. 588–599, 2008.
View at: Publisher Site | Google Scholar
S. Wang, X. Yin, H. Tang, and X. Ge, “A new approach for analyzing circular tunnel in strain-softening rock masses,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 1, pp. 170–178, 2010.
View at: Publisher Site | Google Scholar
Q. Zhang, B.-S. Jiang, S.-l. Wang, X.-r. Ge, and H.-q. Zhang, “Elasto-plastic analysis of a circular opening in strain-softening rock mass,” International Journal of Rock Mechanics and Mining Sciences, vol. 50, pp. 38–46, 2012.
View at: Publisher Site | Google Scholar
L. Cui, J.-J. Zheng, R.-J. Zhang, and Y.-K. Dong, “Elasto-plastic analysis of a circular opening in rock mass with confining stress-dependent strain-softening behaviour,” Tunnelling and Underground Space Technology, vol. 50, pp. 94–108, 2015.
View at: Publisher Site | Google Scholar
S. A. F. Murrell, “The effect of triaxial stress systems on the strength of rocks at atmospheric temperatures,” Geophysical Journal International, vol. 10, no. 3, pp. 231–281, 1965.
View at: Publisher Site | Google Scholar
J. Handin, H. C. Heard, and J. N. Magouirk, “Effects of the intermediate principal stress on the failure of limestone, dolomite, and glass at different temperatures and strain rates,” Journal of Geophysical Research, vol. 72, no. 2, pp. 611–640, 1967.
View at: Publisher Site | Google Scholar
C. Xu and C. Xia, “A new large strain approach for predicting tunnel deformation in strain-softening rock mass based on the generalized Zhang-Zhu strength criterion,” International Journal of Rock Mechanics and Mining Sciences, vol. 143, Article ID 104786, 2021.
View at: Publisher Site | Google Scholar
X. B. Zhang, G. M. Zhao, and X. R. Meng, “Elastoplastic analysis of surrounding rock on circular roadway based on Drucker-Prager yield criterion,” Journal of China Coal Society, vol. 38, no. S1, pp. 30–37, 2013.
View at: Google Scholar
X. Yang, H. Yuan, J. Wu, and S. Li, “Elastoplastic analysis of circular tunnel based on Drucker-Prager criterion,” Advances in Civil Engineering, vol. 2018, Article ID 5149789, 8 pages, 2018.
View at: Publisher Site | Google Scholar
H. X. Rong and Y. M. Hong, “New tri-shear failure criterion and its application in elasto-plastic analysis for the wall rock of tunnel,” Journal of China Coal Society, vol. 28, no. 4, pp. 389–393, 2003.
View at: Google Scholar
K. H. Ceng, H. Y. Ju, G. J. Cheng, and C. G. Zhang, “Elastic-plastic analytical solutions for surrounding rocks of tunnels and its applications,” Journal of China Coal Society, vol. 36, no. 5, pp. 752–755, 2011.
View at: Google Scholar
L. W. Chen, J. B. Peng, F. Wen, and S. Ping, “Analysis of surrounding rock mass plastic zone of round tunnel under non-uniform stress field based on the unified strength theory,” Journal of China Coal Society, vol. 32, no. 1, pp. 20–23, 2007.
View at: Google Scholar
Q. Zhang, H. Y. Wang, S. L. Wang, X. R. Ge, and G. Shao, “Deterioration elasto-plastic analysis of cracked surrounding rocks based on unified strength theory,” Journal of China Coal Society, vol. 35, no. 3, pp. 381–386, 2010.
View at: Google Scholar
L. Chen, X. Mao, Y. Chen, M. Li, Y. Hao, and D. Liu, “A new unified solution for circular tunnel based on a four-stage constitutive model considering the intermediate principal stress,” Advances in Civil Engineering, vol. 2018, Article ID 7912062, 14 pages, 2018.
View at: Publisher Site | Google Scholar
M. Zhu, Y. Yang, F. Gao, and J. Liu, “Analytical solution of tunnel surrounding rock for stress and displacement based on Lade-Duncan criterion,” Advances in Civil Engineering, vol. 2018, Article ID 5363658, 7 pages, 2018.
View at: Publisher Site | Google Scholar
E. T. Brown and E. Hoek, “Trends in relationships between measured in-situ stresses and depth,” International Journal of Rock Mechanics and Mining Science & Geomechanics Abstracts, vol. 15, no. 4, pp. 211–215, 1978.
View at: Publisher Site | Google Scholar
R. T. Ewy, “Wellbore-stability predictions by use of a modified Lade criterion,” SPE Drilling and Completion, vol. 14, no. 2, pp. 85–91, 1999.
View at: Publisher Site | Google Scholar
P. V. Lade, “Three-parameter failure criterion for concrete,” Journal of the Engineering Mechanics Division, vol. 108, no. 5, pp. 850–863, 1982.
View at: Publisher Site | Google Scholar
M. K. Kim and P. V. Lade, “Modelling rock strength in three dimensions,” International Journal of Rock Mechanics and Mining Science & Geomechanics Abstracts, vol. 21, no. 1, pp. 21–33, 1984.
View at: Publisher Site | Google Scholar
B. Bahrami, S. Mohsenpour, M. A. Miri, and R. Mirhaseli, “Quantitative comparison of fifteen rock failure criteria constrained by polyaxial test data,” Journal of Petroleum Science and Engineering, vol. 159, pp. 564–580, 2017.
View at: Publisher Site | Google Scholar
F. Feng, X. Li, K. Du, D. Li, J. Rostami, and S. Wang, “Comprehensive evaluation of strength criteria for granite, marble, and sandstone based on polyaxial experimental tests,” International Journal of Geomechanics, vol. 20, no. 2, Article ID 04019155, 2020.
View at: Publisher Site | Google Scholar
C. M. Zhu, Structural Mechanics, vol. 1, Higher Education Press, Beijing, China, 2004.
G. Z. Sun, Fundamentals of Rock Mass Mechanics, The Science Publishing Company, Beijing, China, 1983.
P. V. Lade, “Rock strength criteria: the theories and the evidence,” in Comprehensive Rock Engineering, J. Hudson and E. T. Brown, Eds., pp. 225–284, Elsevier, Amsterdam, Netherland, 1993.
View at: Google Scholar
Z. Q. Liu and D. M. Yu, “Elastoplastic stress and displacement analytical solutions to deep-buried circular tunnel considering intermediate principal stress and dilatancy,” Engineering Mechanics, vol. 29, no. 8, pp. 289–296, 2012.
View at: Publisher Site | Google Scholar
X. H. Cai and Y. P. Cai, Structural Stress Calculation for Hydraulic Pressure Tunnel, China Water & Power Press, Beijing, China, 2004.
J. S. Sun and W. B. Lu, “Analytical elastoplastic solutions to supporting rock masses of circular tunnels under asymmetric load,” Rock and Soil Mechanics, vol. 28, no. S1, pp. 327–332, 2007.
View at: Google Scholar

Copyright

Copyright © 2021 Yansheng Deng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies