Next Article in Journal
CFD Modeling and Cold Physical Model Simulation on Single Molten Slag Ligament Disintegration into Droplets
Previous Article in Journal
Synthesis of Ye’elimite from Anthropogenic Waste
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Approach for Analyzing Circular Tunnels in Nonlinear Strain-Softening Rock Masses Considering Seepage Force

State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mines, Anhui University of Science and Technology, Huainan 232001, China
*
Author to whom correspondence should be addressed.
Minerals 2023, 13(2), 138; https://doi.org/10.3390/min13020138
Submission received: 11 December 2022 / Revised: 8 January 2023 / Accepted: 16 January 2023 / Published: 17 January 2023
(This article belongs to the Topic Support Theory and Technology of Geotechnical Engineering)

Abstract

:
Accurate calculation of the stresses and deformations of tunnels is of great importance for practical engineering applications. In this study, a three-region model for tunnels considering seepage force was established. A new nonlinear strain-softening model is proposed. Then, a unified solution for the stresses and deformations of tunnels is deduced. Through a series of discussions, the effects of seepage force, softening modulus coefficient of cohesion, and initial support resistance on the stress distribution, radii of the post-peak zone, and surface displacement around the tunnel are discussed. Results show that the tangential stresses are always larger than the radial stresses. As the distance from the tunnel center increases, the radial stress continues to increase, while the tangential stress first increases and then decreases. With the increases in seepage force, the radii of the post-peak zone and surface displacement all increase. With the increases in softening modulus coefficient of cohesion, the radii of the post-peak zone increase while the surface displacement decreases. Tunnels with a higher initial support resistance experience lower radii of the post-peak zone and surface displacement.

1. Introduction

With the rapid development of society and the economy, many huge infrastructure projects have been put into construction. As one of the most common structures, tunnels are widely used in water conservancy engineering and underground traffic engineering [1,2,3]. The stress state in the underground rock will change after the excavation of tunnels [4,5]. Once the changed stress exceeds the ultimate strength of the rock mass, the tunnel surrounding the rock would be destroyed, and three different regions (a damaged region, a strain-softening region, and an elastic region) would be formed, which would exert a significant impact on the stability of the tunnel [6,7,8,9,10]. Therefore, accurate stress analysis of the surrounding rock is very important to the stability control of the tunnel.
The appropriate constitutive models for different types of rock materials are the key to solving the stresses and deformation of tunnels. In past decades, the elastic–brittle-plastic constitutive model (Figure 1a), suitable for hard rock masses with Geological Strength Index (GSI) values above 75, and the elastic–perfectly plastic constitutive model (Figure 1b) for soft rock masses with GSI values below 25, were usually used to research this problem [11,12,13,14,15,16]. However, the strain-softening behavior of rock mass was ignored by the first two constitutive models. The third model, called the linear strain-softening constitutive model (Figure 1c), considering the linear deterioration of strength parameters in the strain-softening region, can stand for a wider range of rock masses with 25 < GSI < 75 [17,18,19,20]. Moreover, the elastic–perfectly plastic constitutive model (the slope of the strain-softening stage is equal to null) and the elastic–brittle-plastic constitutive model (the slope of the strain-softening stage is equal to infinity) are both special cases of the linear strain-softening constitutive model. In fact, the mechanical parameters of rock mass, such as cohesion, tend to decrease nonlinearly in the post-peak stage [21,22]. Thus, the nonlinear strain-softening constitutive model (Figure 1d) was presented and selected as the constitutive model for the rock mass in this study.
Groundwater is one of the important components of underground engineering structures. The seepage force caused by the flow of groundwater also has a non-negligible influence on the stress state around the tunnel [23,24,25]. Therefore, the seepage force is expected to be taken into consideration in the theoretical analysis of tunnels. In order to accurately analyze the stresses and deformations of tunnels in nonlinear strain-softening rock masses considering seepage force. In the present study, a mechanical model for tunnels in water-rich areas is established. A new nonlinear strain-softening model that considers the nonlinear degradation of cohesion is proposed. Based on the Mogi–Coulomb criterion, a unified solution for the stresses and displacements of the surrounding rock is determined. Through a series of case studies, the effects of seepage force, softening modulus coefficient of cohesion, and initial support resistance on the stress distribution, radii of the post-peak region, and surface displacement around the tunnel are discussed.

2. Model Establishment

As shown in Figure 2, assuming that there is a cylindrical tunnel in the water-rich area and the surrounding rock is a nonlinear strain-softening media, the original in-situ stress is σ0. After the tunnel excavation and stress redistribution, a crushed region, strain-softening region, and elastic region were formed in the surrounding rock around the tunnel. The radii of the tunnel and the three regions are sequentially recorded as R0, Rc, Rs, and Re. In addition, an initial support resistance (Pi) was uniformly applied to the excavation surface to maintain the stability of the tunnel.
We assume that the seepage force is evenly distributed outside the elastic region of the tunnel. Based on Darcy’s law, the differential equation for the seepage force is given by the following equation:
d 2 P w d r 2 + 1 r d P w d r = 0
where Pw is the seepage force anywhere around the tunnel.
Combining the boundary conditions, Pw = 0, r = R0 and Pw = P0, and r = Re, the seepage force can be solved as
P w = P 0 ln ( R 0 / r ) ln ( R 0 / R e )
Generally, the cohesion in the strain-softening region decreases nonlinearly, and the softening modulus of cohesion increases with the increase in strain, as shown in Figure 3. We assume that the softening modulus of cohesion is as follows:
E = α ε θ
where α is the softening modulus coefficient of cohesion, which reflects the brittleness of the rock mass. The rock mass is ideal elastic-strain-softening when α = 0, and the greater the softening coefficient, the more brittle the rock.

3. Mechanical Analysis

To solve the three-region model shown in Figure 1, the following basic formulas should be given as [26,27]:
1.
Equilibrium differential formula:
d σ r d r + σ r σ θ r + η d P w d r = 0
where σr and σθ are the radial and tangential stresses, respectively, and η is the seepage force coefficient.
2.
Geometric formula:
{ ε r = d u d r ε θ = u r
where εr and εθ are the radial and tangential strains, respectively.
3.
Constitutive formula:
{ ε r = 1 μ 2 E ( σ r μ 1 μ σ θ ) ε θ = 1 μ 2 E ( σ θ μ 1 μ σ r )
where E is the elastic modulus of rock mass, and μ is the Poisson’s ratio of rock mass.
4.
Mogi-Coulomb strength criterion:
σ θ = M σ r + N
where M = ( 3 + 2 sin φ ) / ( 3 2 sin φ ) and N = 4 c cos φ / ( 3 2 sin φ ) .

3.1. Stresses and Displacement in the Elastic Region

By substituting Equations (2), (5), and (6) into Equation (4), the following differential equation of the displacement in the elastic region can be deduced:
d 2 u e d r 2 + 1 r d u e d r u e r 2 = η F v r E
where v = ( 1 + μ ) ( 1 2 μ ) 1 μ and F = P 0 / ( ln ( R 0 / R e ) ) .
Solving Equation (8), we can obtain
u e = C 1 r + C 2 r + η F v r ln r 2 E
where C1 and C2 are undetermined constants.
By substituting Equation (9) into Equation (5), the following strains can be deduced:
{ ε r e = C 1 C 2 r 2 + η F v 2 E ( ln r + 1 ) ε θ e = C 1 + C 2 r 2 + η F v ln r 2 E
where C1 and C2 are undetermined constants.
By integrating Equations (6) and (10), we can obtain the radial and tangential stresses:
{ σ r e = E C 1 ( 1 + μ ) ( 1 2 μ ) E C 2 ( 1 + μ ) r 2 + η F ln r 2 ( 1 μ ) + η F 2 σ θ e = E C 1 ( 1 + μ ) ( 1 2 μ ) + E C 2 ( 1 + μ ) r 2 + η F ln r 2 ( 1 μ ) + μ η F 2 ( 1 μ )
With the stress boundary conditions σ r = σ r e s , r = Rs and σ r = σ 0 + P 0 at r = Re, C1 and C2 can be deduced as follows:
{ C 1 = ν ( 1 μ ) E ( σ 0 + P 0 ) + R s 2 R e 2 R s 2 v ( 1 μ ) E ( σ 0 + P 0 σ r e s ) + η F v 2 E R s 2 R e 2 R s 2 ln R s R e η F v ( 1 μ ) 2 E η F v 2 E ln R e C 2 = 1 + μ E R e 2 R s 2 R e 2 R s 2 ( σ 0 + P 0 σ r e p ) + η F ( 1 + μ ) 2 E ( 1 μ ) R e 2 R s 2 R e 2 R s 2 ln R s R e
By substituting Equation (12) into Equation (11), the stresses in the elastic region can be solved:
{ σ r e = σ 0 + P 0 + η F 2 ( 1 μ ) ln r R e + R s 2 R s 2 R e 2 ( R e 2 r 2 1 ) [ σ 0 + P 0 σ r e s + η F 2 ( 1 μ ) ln R s R e ] σ θ e = σ 0 + P 0 + η F 2 ln ( r / R e ) + 2 μ 1 1 μ R s 2 R s 2 R e 2 ( R e 2 r 2 + 1 ) [ σ 0 + P 0 σ r e s + η F 2 ( 1 μ ) ln R s R e ]
The radial and tangential stresses should satisfy Equation (7) at r = Rs, so we can obtain
σ r e s = 2 R e 2 ( σ 0 + P 0 ) + η F 2 ( 1 μ ) [ ( R s 2 R e 2 ) ( ln R s R e + 2 μ 1 ) ( R s 2 + R e 2 ) ln R s R e ] N e ( R s 2 R e 2 ) M ( R s 2 R e 2 ) R s 2 R e 2
By substituting Equation (11) into Equation (6), the radial and tangential strains in the elastic region can be deduced:
{ ε r e = C 1 C 2 r 2 + η F v 2 E ( ln r + 1 ) ε θ e = C 1 + C 2 r 2 + η F v ln r 2 E
Combining with Equations (5) and (15), the displacement in the elastic region can be derived:
u e = C 1 r + C 2 r + η F v r ln r 2 E

3.2. Stresses and Displacement in the Strain-Softening Region

Generally, the volume of roadway surrounding rock will expand after exaction. The relationship between the dilation coefficient and strains is shown in Figure 4.
According to the non-associated flow law, the relationship between radial and tangential strains in strain-softening and crushed regions can be denoted as
{ ε r s + β s ε θ s = 0 ε r c + β c ε θ c = 0
where βs and βc are the dilation coefficient in strain-softening and crushed regions, respectively. They can be obtained by
{ β s = 1 + sin ψ s 1 sin ψ s β c = 1 + sin ψ c 1 sin ψ c
where ψ s and ψ c are the dilation angle in strain-softening and crushed regions, respectively.
The total strains consist of elastic and strain-softening parts in the strain-softening region:
{ ε r = ε r s + ε r e s ε θ = ε θ s + ε θ e s
where ε r e s and ε θ e s are the radial and tangential strains on the interface between the elastic and strain-softening regions, respectively.
By integrating Equations (5), (17), and (19), the differential equation for displacement in the strain-softening region can be deduced:
d u s d r + β s u s r = ε r e s + β s ε θ e s
When u s = u e s and r = Rs, Equation (20) can be solved:
u s = [ u e s ε r e s R s β s + 1 ] ( R s r ) β s + ( ε r e s + β s ε θ e s ) r β s + 1
The strains in the strain-softening region can be deduced by substituting Equation (21) into Equation (5):
{ ε r s = β s ( ε r e s ε θ e s ) β s + 1 ( R s r ) β s + 1 + ε r e s + β s ε θ e s β s + 1 ε θ s = ε θ e s ε r e s β s + 1 ( R s r ) β s + 1 + ε r e s + β s ε θ e s β s + 1
As shown in Figure 2b, the cohesion in the nonlinear strain-softening region can be written as
c s = c 0 α ε θ s ( ε θ s ε θ e s ) = c 0 α [ ε θ e s ε r e s β s + 1 ( R s r ) β s + 1 + ε r e s + β s ε θ e s β s + 1 ] [ ε θ e s ε r e s β s + 1 ( R s r ) β s + 1 + ε r e s + β s ε θ e s β s + 1 ε θ e s ] = c 0 α ( ε θ e s ε r e s ) β s + 1 [ ε θ e s ε r e s β s + 1 ( R s r ) 2 ( β s + 1 ) + 2 ε r e s + ( β s 1 ) ε θ e s β s + 1 ( R s r ) β s + 1 ε r e s + β s ε θ e s β s + 1 ]
where c0 is the initial cohesion.
By substituting Equations (2), (7) and (23) into Equation (4), the differential equation for the radial stress in the strain-softening region can be deduced:
d σ r d r + ( 1 M ) σ r r + A ( R s r ) 2 ( β s + 1 ) + B ( R s r ) β s + 1 D r N e + η F r = 0
where
{ A = 4 cos φ α ( ε θ e s ε r e s ) 2 ( 3 2 sin φ ) ( β p + 1 ) 2 B = 4 cos φ α ( ε θ e s ε r e s ) ( β p ε θ e s ε θ e s + 2 ε r e s ) ( 3 2 sin φ ) ( β s + 1 ) 2 D = 4 cos φ α ( ε θ e s ε r e s ) ( ε r e s + β s ε θ e s ) ( 3 2 sin φ ) ( β s + 1 ) 2
When r = Rs and σ r = σ r e s , the stresses in the strain-softening region can be solved:
{ σ r s = ( σ r e s N e + η F + D 1 M ) ( R s r ) 1 M A 2 β s + M + 1 [ ( R s r ) 1 M ( R s r ) 2 ( β s + 1 ) ] B β s + M [ ( R s r ) 1 M ( R s r ) β s + 1 ] + N e + η F + D 1 M σ θ s = M { ( σ r e s N e + η F + D 1 M ) ( R s r ) 1 M A 2 β s + M + 1 [ ( R s r ) 1 M ( R s r ) 2 ( β s + 1 ) ] B β s + M [ ( R s r ) 1 M ( R s r ) β s + 1 ] + N e + η F + D 1 M } + N e A ( R s r ) 2 ( β s + 1 ) B ( R s r ) β s + 1 + D

3.3. Stresses and Displacement in the Crushed Region

The total strains consist of strain-softening and crushed parts in the crushed region and we can obtain
{ ε r = ε r c + ε r s c ε θ = ε θ c + ε θ s c
where ε r s d and ε θ s d are the radial and tangential strains at r = Rc, respectively.
By integrating Equations (5), (17), and (27), the differential equation for the displacement in the crushed region can be deduced:
d u c d r + β c u c r = ε r s c + β c ε θ s c
When r = Rc and u c = u s c , the displacement in the strain-softening region can be solved:
u c = [ u s c ε r s c R c β c + 1 ] ( R c r ) β c + ( ε r s c + β c ε θ s c ) r β c + 1
Combining Equations (5) and (29), the strains in the crushed region can be deduced:
{ ε r c = β c ( ε r s c ε θ s c ) β p + 1 ( R c r ) β s + 1 + ε r s c + β c ε θ s c β d + 1 ε θ c = ε θ s c ε r s c β c + 1 ( R c r ) β c + 1 + ε r s c + β c ε θ s c β c + 1
By substituting Equations (2) and (7) into Equation (4), the differential equation for radial stress in the crushed region can be obtained:
d σ r d r + ( 1 M ) σ r N d r η F r = 0
When r = R0 and σ r d = P i , the stresses in the crushed region can be derived:
{ σ r c = ( p i N c + η F 1 M ) ( R 0 r ) 1 M + N c + η F 1 M σ θ c = M ( p i N c + η F 1 M ) ( R 0 r ) 1 M + M N c + η F 1 M

3.4. Radii of Strain-Softening and Crushed Regions

When r = Rc and σ r s = σ r c , combining Equations (26) and (33), we can obtain
( σ r e s N e + η F 1 M ) ( R s R c ) 1 M + N e η F 1 M + 4 cos φ / ( 3 2 sin φ ) α ( ε θ e s ε r e s ) ( β s + 1 ) ( β s + M ) [ ( R s R c ) 1 β s ( R s R c ) 1 M ] + 4 cos φ / ( 3 2 sin φ ) α ( ε θ e s ε r e s ) ( β s + 1 ) ( 1 M ) [ 1 ( R s R c ) 1 M ] = ( p i N c + η F 1 M ) ( R 0 R c ) 1 M + N c + η F 1 M
Based on Equation (23), the cohesion at r = Rc can be deduced:
c c = c 0 α ( ε θ e s ε r e s ) β s + 1 [ ε θ e s ε r e s β s + 1 ( R s r ) 2 ( β s + 1 ) + 2 ε r e s + ( β s 1 ) ε θ e s β s + 1 ( R s r ) β s + 1 ε r e s + β s ε θ e s β s + 1 ]
Combining Equations (33) and (34), Rs and Rc can be obtained by an iterative method. The flowchart of the above calculation can be shown in Figure 5.

4. Special Cases Discussion

4.1. Kastner’s Formula

When α = 0, ψs = ψc = 0, c0 = cc, d η = 0, the current nonlinear strain-softening model degenerates for the elastic–perfectly plastic model, and the radius of the plastic region can be determined by the following equation:
R s = R 0 ( ( 2 σ 0 N e ) / ( 1 + M ) N e / ( 1 M ) P i N e / ( 1 M ) ) 1 / ( M 1 )
Equation (35) is Kastner’s formula [15].

4.2. Wilson’s Formula

When α = ∞, ψs = ψc = 0, and η = 0, the current nonlinear strain-softening model degenerates for the elastic–brittle-plastic model, and the radius of the plastic region can be determined by the following equation:
R s = R 0 ( ( 2 σ 0 N e ) / ( 1 + M ) N c / ( 1 M ) P i N c / ( 1 M ) ) 1 / ( M 1 )
Equation (36) is Wilson’s formula [16].

5. Parameter Sensitivity Analysis

As shown in Figure 6, taking the main return laneway in Qingdong Coal Mine in Anhui Province, China, as an example, the influence of some important factors on the stresses and deformations is discussed next. The geometrical and mechanical parameters are as follows: R0 = 3 m, σ0 = 18.9 MPa, Pi = 0.4 MPa, P0 = 4 MPa, η = 1, Re/R0 = 20, μ = 0.25, E = 2550 MPa, c0 = 4.5 MPa, cc = 1.1 MPa, φ = 26°, ψs = ψc = 10°, and α = 3×104 MPa.

5.1. Seepage Force

Figure 7 shows the effect of seepage force on the radii of strain-softening and crushed regions and surface displacement around the tunnel. It can be seen that the Rs, Rc, and u0 values all increase with the increasing P0. For example, as P0 increases from 3.0 MPa to 6.0 MPa, the Rs, Rc, and u0 increase by 1.61 m, 1.45 m, and 0.14 m, with an increment of 32.07%, 37.21%, and 91.98%, respectively. Therefore, the larger seepage force will aggravate the deformation of the tunnel.
The radial and tangential stresses in the tunnel surrounding rock with different seepage forces are shown in Figure 8. It can be seen that the radial stress is always less than the tangential stress. The radial stress at the tunnel wall is consistent with the initial support resistance, and it presents an increasing trend with the increase in radius. The tangential stress first increases in the broken and strain-softening regions and then decreases in the elastic region, and the peak value of tangential stress appears at the interface between elastic and strain-softening regions. The peak value of tangential stress increases with the increase in seepage force. For example, as the P0 value increases from 3.0 MPa to 6.0 MPa, the peak value of tangential stress increase by 6.88 MPa, with an increment of 16.95%.

5.2. Softening Modulus Coefficient of Cohesion

The effect of softening modulus coefficient of cohesion on the radii of strain-softening and crushed regions and surface displacement around the tunnel is shown in Figure 9. It can be seen that the radii of strain-softening and crushed regions both increase with the increase in softening modulus coefficient. However, the ratio of the radii of the above two regions is getting smaller with the increasing softening modulus coefficient. For example, as α increases from 3 × 104 MPa to 8.0 × 104 MPa, the ratio of Rs to Rc decreases from 1.27 to 1.13, with a reduction of 11.02%. The surface displacement also decreases with the increase in softening modulus coefficient. For example, the surface displacement is 0.19 m when the softening modulus coefficient is 3 × 104 MPa, and surface displacement reduces to 0.18 m when the softening modulus coefficient reaches 8 × 104 MPa.
The radial and tangential stresses in the tunnel surrounding rock with different softening modulus coefficients of cohesion are shown in Figure 10. It can be seen that the tangential stress is always larger than the radial stress. As the distance from the tunnel center increases, the radial stress continues to increase, and the tangential stress first increases in the crushed and strain-softening regions and then decreases in the elastic region, which is similar to the laws in Figure 8. However, the radial and tangential stresses in the crushed region are very close under different softening modulus coefficients of cohesion. Additionally, as the softening modulus coefficient increases, the peak value of tangential stress does not change, which means the sensitivity of softening modulus coefficient to stress distribution is quite low.

5.3. Initial Support Resistance

Figure 11 shows the effect of initial support resistance on the radii of strain-softening and crushed regions and surface displacement around the tunnel. It can be seen that the Rs, Rc, and u0 values all decrease with the increasing Pi. For example, as the Pi value increases from 0 to 1.0 MPa, the Rs, Rc, and u0 values decrease by 1.26 m, 0.99 m, and 0.09 m, with a reduction of 20.49%, 20.45%, and 38.51%, respectively. Therefore, the larger initial support resistance can control the deformation of the tunnel. As a result, some support measures, such as concrete lining and bolts, can be used to increase the initial support resistance and ensure the stability of the tunnel.
The radial and tangential stresses in the tunnel surrounding rock based on different initial support resistance are shown in Figure 12. It can be seen that with the increase in initial support resistance, the whole radial stress and the tangential stress in the crushed and strain-softening regions increases while the tangential stress in the elastic region decreases. However, the peak value of tangential stress remains constant as the initial support resistance changes.

6. Conclusions

Considering the influence of seepage force, a closed-form solution for the stress distribution and deformations around the tunnels in the water area is deduced based on the new nonlinear strain-softening model. Moreover, the effects of seepage force, softening modulus coefficient of cohesion, and initial support resistance on the stress distribution, radii of the post-peak zone, and surface displacement are also discussed. The conclusions can be summarized as follows:
(1)
Compared with other solutions, when the softening modulus coefficient of cohesion and seepage force coefficient are both zero, the current solution degenerates for the elastic–perfectly plastic solution by Kastner. When the softening modulus coefficient of cohesion is large enough, and the seepage force coefficient is zero, the current solution degenerates for the elastic–brittle-plastic solution by Wilson.
(2)
The tangential stress is always larger than the radial stress. As the distance from the tunnel center increases, the radial stress continues to increase, and the tangential stress first increases in the broken and strain-softening regions and then decreases in the elastic region. The peak value of tangential stress appears at the interface between elastic and strain-softening regions.
(3)
As the seepage force increases, the radii of strain-softening and crushed regions and surface displacement all increase. As the softening modulus coefficient of cohesion increases, the radii of strain-softening and crushed regions increase while the surface displacement decreases. Additionally, the radii of strain-softening and crushed regions and surface displacement all decrease with the increasing initial support resistance.
This study only analyzes the stresses and deformations of tunnels through theoretical methods. Further research on numerical simulation should be carried out in the future.

Author Contributions

Conceptualization, H.F. and L.W.; methodology, H.F.; validation, L.W. and S.L.; formal analysis, H.F.; resources, S.L.; data curation, S.L.; writing—original draft preparation, H.F.; writing—review and editing, L.W.; supervision, S.L.; funding acquisition, L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This study is supported by the Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology (No. 2021yjrc15), Major Special Projects of Science and Technology of Anhui Province (No. 202203a07020010), Collaborative Innovation Project of Anhui Province (GXXT-2020-055), and Opening Foundation of the State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mines (No. SKLMRDPC21KF09).

Data Availability Statement

The data used for conducting classifications are available from the corresponding author upon request.

Acknowledgments

The authors would like to thank the engineers in Qingdong Coal Mine for their engineering data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Detournay, E. Elastoplastic model of a deep tunnel for a rock with variable dilatancy. Rock Mech. Rock Eng. 1986, 19, 99–108. [Google Scholar] [CrossRef]
  2. González-Cao, J.; Varas, F.; Bastante, F.G.; Alejano, L.R. Ground reaction curves for circular excavations in non-homogeneous, axisymmetric strain-softening rock masses. J. Rock Mech. Geotech. Eng. 2013, 5, 431–442. [Google Scholar] [CrossRef] [Green Version]
  3. Luo, Y.; Gong, F.Q.; Zhu, C.Q. Experimental investigation on stress-induced failure in D-shaped hard rock tunnel under water-bearing and true triaxial compression conditions. Bull. Eng. Geol. Environ. 2022, 81, 76. [Google Scholar] [CrossRef]
  4. Liang, R.Z.; Wu, J.; Sun, L.W.; Shen, W.; Wu, W.B. Performances of adjacent metro structures due to regiond excavation of a large-scale basement in soft ground. Tunn. Undergr. Space Technol. 2021, 117, 104123. [Google Scholar] [CrossRef]
  5. Tonon, R.C.F. Closed-form solutions for a circular tunnel in elastic-brittle-plastic ground with the original and generalized Hoek–Brown failure criteria. Rock Mech. Rock Eng. 2011, 44, 169–178. [Google Scholar]
  6. Lv, A.; Masoumi, H.; Walsh, S.D.C.; Roshan, H. Elastic-Softening-Plasticity around a borehole: An analytical and experimental study. Rock Mech. Rock Eng. 2019, 52, 1149–1164. [Google Scholar] [CrossRef]
  7. Carranza-Torres, C. Elasto-plastic solution of tunnel problem using the generalized form of the Hoek-Brown failure criterion. Int. J. Rock Mech. Min. 2004, 41, 480–481. [Google Scholar] [CrossRef]
  8. Hao, Y.; Wu, Y.; Mao, X.; Li, P.; Zhang, L.; Tao, J. An elastoplastic softening damage model for hydraulic fracturing in soft coal seams. Adv. Civ. Eng. 2018, 2018, 2548217. [Google Scholar] [CrossRef] [Green Version]
  9. Park, K.H.; Tontavanich, B.; Lee, J.G. A simple procedure for ground response curve of circular tunnel in elastic-strain softening rock masses. Tunn. Undergr. Space Technol. 2008, 23, 151–159. [Google Scholar] [CrossRef]
  10. Pourhosseini, O.; Shabanimashcool, M. Development of an elasto-plastic constitutive model for intact rocks. Int. J. Rock Mech. Min. Sci. 2014, 66, 1–12. [Google Scholar] [CrossRef]
  11. Brown, E.; Bray, J.; Ladanyi, B.; Hoek, E. Ground response curves for rock tunnels. J. Geotech. Eng. 1983, 109, 15–39. [Google Scholar] [CrossRef]
  12. Masoudian, M.S.; Hashemi, M.A. Analytical solution of a circular opening in an axisymmetric elastic-brittle-plastic swelling rock. J. Nat. Gas Sci. Eng. 2016, 35, 483–496. [Google Scholar] [CrossRef]
  13. Sharan, S.K. Elastic-brittle-plastic analysis of circular openings in Hoek-Brown media. Int. J. Rock Mech. Min. 2003, 40, 817–824. [Google Scholar] [CrossRef]
  14. Zareifard, M.; Fahimifar, A. Analytical solutions for the stresses and deformations of deep tunnels in an elastic-brittle-plastic rock mass considering the damaged region. Tunn. Undergr. Space Technol. 2016, 58, 186–196. [Google Scholar] [CrossRef]
  15. Kastner, H. Statik des Tunnel-Und Stollenbaues; Springer: New York, NY, USA, 1962. [Google Scholar]
  16. Wilson, A.H. A method of estimating the closure and strength of lining required in drivages surrounded by a yield zone. Int. J. Rock Mech. Min. Geomech. Abstr. 1980, 17, 349–3550. [Google Scholar] [CrossRef]
  17. Ghorbani, A.; Hasanzadehshooiili, H.; Sadowski, L. Neural prediction of tunnels’ support pressure in elasto-plastic, strain-softening rock mass. App. Sci. 2018, 8, 841. [Google Scholar] [CrossRef] [Green Version]
  18. Han, J.X.; Li, S.C.; Li, S.C.; Yang, W.M. A procedure of strain-softening model for elasto-plastic analysis of a circular opening considering elasto-plastic coupling. Tunn. Undergr. Space Technol. 2013, 37, 128–134. [Google Scholar] [CrossRef]
  19. Li, Y.; Cao, S.; Fantuzzi, N.; Liu, Y. Elastoplastic analysis of a circular borehole in elastic-strain softening coal seams. Int. J. Rock Mech. Min. 2015, 80, 316–324. [Google Scholar] [CrossRef]
  20. Wang, S.L.; Yin, X.T.; Tang, H.; Ge, X. A new approach for analyzing circular tunnel in strain-softening rock masses. Int. J. Rock Mech. Min. 2010, 47, 170–178. [Google Scholar] [CrossRef]
  21. Jing, W.; Xue, W.; Yao, Z. Variation of the internal friction angle and cohesion of the plastic softening region rock in roadway surrounding rock. J. China Coal Soc. 2018, 43, 2203–2210. [Google Scholar]
  22. Lu, Y.; Wang, L.; Yang, F.; Li, Y.; Chen, H. Post-peak strain softening mechanical properties of weak rock. Chin. J. Rock Mech. Eng. 2010, 29, 640–648. [Google Scholar]
  23. Fahimifar, A.; Ghadami, H.; Ahmadvand, M. The influence of seepage and gravitational loads on elastoplastic solution of circular tunnels. Sci. Iran. 2014, 21, 1821–1832. [Google Scholar]
  24. Fahimifar, A.; Zareifard, M.R. A thoretical solution for analysis of tunnels below groundwater considering the hydraulic-mechanical coupling. Tunn. Undergr. Space Technol. 2009, 24, 634–646. [Google Scholar] [CrossRef]
  25. Huangfu, M.; Wang, M.S.; Tan, Z.S.; Wang, X.Y. Analytical solution for steady seepage into an underwater circular tunnel. Tunn. Undergr. Space Technol. 2010, 25, 391–396. [Google Scholar] [CrossRef]
  26. Fan, H.; Wang, L.G.; Wang, S.; Jiang, C.Y. A new unified solution for deep tunnels in water-rich areas considering seepage force. Geofluids 2021, 2021, 6696757. [Google Scholar] [CrossRef]
  27. Mogi, K. Flow and fracture of rocks under general triaxial compression. Appl. Math. Mech. 1981, 2, 635–651. [Google Scholar] [CrossRef]
Figure 1. Different constitutive models. (a) elastic–perfectly plastic; (b) elastic–brittle-plastic; (c) linear strain-softening; (d) nonlinear strain-softening.
Figure 1. Different constitutive models. (a) elastic–perfectly plastic; (b) elastic–brittle-plastic; (c) linear strain-softening; (d) nonlinear strain-softening.
Minerals 13 00138 g001
Figure 2. Mechanical models. (a) Three-region model of a tunnel cylindrical considering seepage force. (b) Nonlinear strain-softening model of surrounding rock.
Figure 2. Mechanical models. (a) Three-region model of a tunnel cylindrical considering seepage force. (b) Nonlinear strain-softening model of surrounding rock.
Minerals 13 00138 g002
Figure 3. Nonlinear strain-softening model of cohesion.
Figure 3. Nonlinear strain-softening model of cohesion.
Minerals 13 00138 g003
Figure 4. Dilation model of rock mass.
Figure 4. Dilation model of rock mass.
Minerals 13 00138 g004
Figure 5. Flowchart of the calculation.
Figure 5. Flowchart of the calculation.
Minerals 13 00138 g005
Figure 6. Location of the main return laneway.
Figure 6. Location of the main return laneway.
Minerals 13 00138 g006
Figure 7. Effect of seepage force on the radii of strain-softening and crushed regions and surface displacement. (a) Rs and Rc; (b) u0.
Figure 7. Effect of seepage force on the radii of strain-softening and crushed regions and surface displacement. (a) Rs and Rc; (b) u0.
Minerals 13 00138 g007
Figure 8. Effect of seepage force on the stress distribution in the tunnel surrounding rock. (a) σr; (b) σθ.
Figure 8. Effect of seepage force on the stress distribution in the tunnel surrounding rock. (a) σr; (b) σθ.
Minerals 13 00138 g008
Figure 9. Effect of softening modulus coefficient of cohesion on the radii of strain-softening and crushed regions and surface displacement. (a) Rs and Rc; (b) u0.
Figure 9. Effect of softening modulus coefficient of cohesion on the radii of strain-softening and crushed regions and surface displacement. (a) Rs and Rc; (b) u0.
Minerals 13 00138 g009
Figure 10. Effect of softening modulus coefficient of cohesion on the stress distribution in the tunnel surrounding rock. (a) σr; (b) σθ.
Figure 10. Effect of softening modulus coefficient of cohesion on the stress distribution in the tunnel surrounding rock. (a) σr; (b) σθ.
Minerals 13 00138 g010
Figure 11. Effect of initial support resistance on the radii of strain-softening and crushed regions and surface displacement. (a) Rs and Rc; (b) u0.
Figure 11. Effect of initial support resistance on the radii of strain-softening and crushed regions and surface displacement. (a) Rs and Rc; (b) u0.
Minerals 13 00138 g011
Figure 12. Effect of initial support resistance on the stress distribution in the tunnel surrounding rock. (a) σr; (b) σθ.
Figure 12. Effect of initial support resistance on the stress distribution in the tunnel surrounding rock. (a) σr; (b) σθ.
Minerals 13 00138 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fan, H.; Wang, L.; Li, S. A New Approach for Analyzing Circular Tunnels in Nonlinear Strain-Softening Rock Masses Considering Seepage Force. Minerals 2023, 13, 138. https://doi.org/10.3390/min13020138

AMA Style

Fan H, Wang L, Li S. A New Approach for Analyzing Circular Tunnels in Nonlinear Strain-Softening Rock Masses Considering Seepage Force. Minerals. 2023; 13(2):138. https://doi.org/10.3390/min13020138

Chicago/Turabian Style

Fan, Hao, Lei Wang, and Shaobo Li. 2023. "A New Approach for Analyzing Circular Tunnels in Nonlinear Strain-Softening Rock Masses Considering Seepage Force" Minerals 13, no. 2: 138. https://doi.org/10.3390/min13020138

APA Style

Fan, H., Wang, L., & Li, S. (2023). A New Approach for Analyzing Circular Tunnels in Nonlinear Strain-Softening Rock Masses Considering Seepage Force. Minerals, 13(2), 138. https://doi.org/10.3390/min13020138

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop