Abstract

The small strain theory underestimates the self-bearing capacity of rock masses, especially for a soft rock tunnel under high geostress. To perform an efficient and accurate calculation and provide a reference for the stiffness design of a tunnel, the finite strain solution for a circular tunnel in Mohr–Coulomb strain-softening rock masses with a non-associated flow rule was derived as three sets of differential equations under the Lagrangian coordinate, which are in the residue region, the softening region, and the elastic region, respectively. Based on the bisection method, an iteration procedure for solving the finite strain solution was proposed to approximate the boundary condition at infinity, the values of two adjacent boundaries, and the initial values on the excavation boundary. This numerical procedure was verified by comparing with self-similar solutions, recursive solutions, and FLAC simulation results. In the calculation example, the relative error on boundaries can be decreased to less than 10−8 after only 10 times iteration and the time for each calculation is less than 15 s. Applying this procedure on the sensibility analysis and stiffness reliability design for the Zhongyi tunnel, a support stiffness of 4.3 MPa/m is recommended to guarantee a tunnel displacement lower than 0.5 m.

Highlights

  • With the rapid construction of transportation infrastructure, more and more tunnels have been built in soft rock under high geostress

  • Similar problems are encountered in the Muzhailing Railway Tunnel [3], Zhegu Mountain Tunnel [4], Huangjiazhai Tunnel [5], Minxian Tunnel [6] and have been predicted to appear during the construction of 43 soft rock tunnels with high geostress in the Sichuan–Tibet railway [7]

  • Scholars have been conducting considerable research on the mechanical calculation of tunnels excavated in strain-softening rock masses

Read more

Summary

A New Numerical Finite Strain Procedure for a Circular Tunnel

Citation: Chen, W.; Zhang, D.; Fang, Abstract: The small strain theory underestimates the self-bearing capacity of rock masses, especially for a soft rock tunnel under high geostress. To perform an efficient and accurate calculation and provide a reference for the stiffness design of a tunnel, the finite strain solution for a circular tunnel in Mohr–Coulomb strain-softening rock masses with a non-associated flow rule was derived as three sets of differential equations under the Lagrangian coordinate, which are in the residue region, the softening region, and the elastic region, respectively. Based on the bisection method, an iteration procedure for solving the finite strain solution was proposed to approximate the boundary condition at infinity, the values of two adjacent boundaries, and the initial values on the excavation boundary. This numerical procedure was verified by comparing with self-similar solutions, recursive solutions, and FLAC simulation results. In the calculation example, the relative error on boundaries can be decreased to less than 10−8 after only 10 times iteration and the time for each calculation is less than 15 s. Applying this procedure on the sensibility analysis and stiffness reliability design for the Zhongyi tunnel, a support stiffness of 4.3 MPa/m is recommended to guarantee a tunnel displacement lower than 0.5 m. Keywords: finite strain; circular tunnel; strain-softening rock masses; global sensibility analysis; stiffness design Application. Appl. Sci. 2022, 12, 2706. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Licensee MDPI, Basel, Switzerland. 4.0/).

Introduction
Model Description
Governing
Softening Parameters, Yield Criterion, and Evolutional Law
Solutions for Stress and Displacement of Rock Mass
Solution in the Residue Region
Solution in the Softening Region
Solution in the Elastic Region
Numerical Implementation Procedure
Numerical Procedures for Calculating Limit Support Pressure
Numerical Procedure for Elastic State
Numerical Procedure
Numerical Procedure for Elastic-Softening-Residue State
Verification for Finite
Procedure
Following thethe tening residue region listed in Appendix
Design
MPa/m, can that the tunnel the tunnel displacement u will not exceed
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call