Abstract

Tunnels are generally designed for a sustained usage of 80 to 100 years, during which the safety of tunnel structures must be guaranteed. A common supporting form utilized in contemporary tunnel engineering is composite lining. To derive applicable parameters of the supporting form and therefore ensure the long-term safety of the tunnel structure, it is imperative to determine the extra acting force exerted onto the composite lining by the creep of the rock surrounding the tunnel and to calculate the stress-strain characteristics of composite lining. In the current study, this paper proposes an approach termed surrounding reinforcement, which is based on the homogenization method. Specifically, this paper defined the bolt force as the internal force of the surrounding rock, analyzed their viscoelastic-plastic properties using the unified strength theory, and derived an equation for calculating the stress-strain relationship of the composite lining. To further validate the method in tunnel structures, this paper applied the derived equation to a representative instance. The results of this paper show that the initial support force has also increased during the creep process of the surrounding rock, indicating that engineers should pay close attention to the coordination between the strength of initial support and the secondary lining and thus ensure an optimal distribution of the pressure from the surrounding rock when designing composite lining tunnel within weak strata. This paper proposes that the initial support not only would guarantee the tunnel safety during the construction stage but also could cooperate with the secondary lining to brace the stress caused by the creep, ensuring that the supporting structure stays stable across the whole period of tunnel operation. This paper provides an alternative to previous methods that is more comprehensive, with simpler calculations, and more applicable to the composite lining supporting design within weak strata.

1. Introduction

As a permanent underground building, tunnels need to ensure the stability of surrounding rock during excavation and need to be able to resist additional forces caused by the long-term creep of the surrounding rock during the tunnel operation [1]. Composite lining structure is the most common form of support in present-day tunnel engineering. Without feasible designing parameters for the secondary lining [2, 3], weak strata that have already shown conspicuous creep characteristics will eventually fail due to the excessive deformation pressure [48]. As shown in Figure 1, common forms of failed designs include tunnel roof subsidence, tunnel lining destruction, tunnel lateral wall destruction, and the destruction of tunnel support steel arch. Therefore, it is an urgent problem to analyze the stress characteristics of composite lining under the creep condition of the surrounding rock from the mechanical point of view and design the secondary lining reasonably [9].

At present, the secondary lining support mechanism considering the creep is mainly studied by means of numerical simulation [1012], theoretical analysis, and field data feedback [13]. Xu et al. [14], combined with the uniaxial compression creep test, proposed a numerical method based on the particle discrete element method to describe the typical creep behavior of phyllite. The failure process of the secondary tunnel lining in the phyllite stratum was evaluated by the numerical model method. The results showed that the failure process of the secondary lining was gradual, with four different stages. Huang et al. [15], based on the concept of disturbance state, proposed a new creep constitutive model for soft rocks. It played an essential role in the study of time-dependent damage mechanisms governing the deformation behavior of surrounding rocks. Wu et al. [16] divided the whole deformation of lining into two parts for processing: the deformation of highly elements and that of the segmental lining. The viscoelastic analytical solutions of the lining pressure and tunnel displacement at the two deformation stages were deduced respectively. Paraskevopoulou and Diederichs [17], based on mechanical analysis, proposed a convergence-constraint method considering the time effect. Through the above method, the trend of the tunnel lining stress over time could be predicted. Xu et al. [18] investigated the failure process of the secondary lining using the similarity model test and the particle-discrete-element numerical approach. The results showed that the results of the similar model test and numerical simulation were consistent in the order of crack occurring and the final failure pattern. Guan et al. [19], based on the classic Burger-MC rheological model, proposed a Burger-deterioration rheological model composed of a Kelvin unit, a Maxwell unit, and a deteriorating MC unit in series. It could be used to explain the mechanical problems for the delayed deformation after the tunnel was put into use. Wang et al. [20] proposed a high-efficiency and high-precision segmented optimization method. The authors aimed to optimize the minimum reserved deformation layer thickness and to evaluate the optimization by the plastic strain of the secondary lining after 100 years of creep occurrence in the surrounding rock. Chu et al. [21] deduced the expressions of the lining resistance, the lining displacement, and the displacement of surrounding rock in the context of the rheological properties of the surrounding rock. The method for determining the timing was studied in depth, and the values obtained from actual measurements were compared to the theoretical values of the surrounding rock vault subsidence to verify the theoretical accuracy. Wang et al. [22] used the finite element program ANSYS to establish a calculation model and analyzed the force law of the composite lining under different supporting times and different stiffness conditions for the secondary lining. The authors proposed that when designing the composite lining, the initial support should bear all basic loads, while the secondary lining should serve as a safety reserve, and consequently the primary support and the secondary lining could share their specific loads. Cristescu et al. [23] believed that the rock followed the law of non-Newtonian fluid and introduced a nonlinear creep model to analyze the tunnel creep law. Subsequently, they investigated the influence of the surrounding rock creep on the force deformation of the lining based on this model.

Although previous studies described above have done substantial research on the force of the composite lining in the creep process of tunnel surrounding rock, there are several limitations:(1)At present, there is no systematic theoretical understanding of the load transfer law of tunnel lining structure that considers creep action of tunnel surrounding rock, which leads to the lack of a reliable theoretical basis for composite lining design under soft rock conditions.(2)There is an insufficient study on the rationality of the initial support and secondary support load distribution in the creep process of the tunnel surrounding rock.

To address these limitations, the research ideas of this article are as follows.

Based on the unified strength theory [2427], the bolt supporting force is regarded as an internal force applied to the surrounding reinforcement, and its mechanical model is established to derive the expression of the tunnel stress and deformation under the composite lining effect. The generalized Kelvin model is used to analyze the creep characteristics of the surrounding rock [2831], and the law of stress law of the initial support and secondary lining in the creep process of the surrounding rock is obtained. A more reasonable design concept of composite lining is proposed to ensure the stability of the supporting structure and the surrounding rock during the tunnel operation.

2. Establishment of the Analytical Model

It is assumed that the radius of the circular tunnel is , in situ stress is , the supporting force applied by lining is , the supporting force applied by bolt is , and the radius of the plastic zone is . These include shotcrete, steel arch, and other supporting forces provided by the initial support and forces provided by the secondary lining. The bolt and the rock and soil in its anchorage range constitute the reinforcement body of the surrounding rock. The elastic-plastic stress of the tunnel should be solved as the plane strain problem. The basic assumptions of the mechanical model are listed as follows: (1) the tunnel section is equivalent to a circular one with infinite length. (2) The original rock stress is hydrostatic pressure, that is, the lateral pressure coefficient is 1 [32]. (3) The surrounding rock is homogeneous, isotropic, and incompressible material. (4) The surrounding rock reinforcement body is a homogeneous medium. The supporting force of bolt in the surrounding rock reinforcement body is uniformly loaded, and the radius of action is . The stress changes linearly with the depth of the surrounding rock. (5) The influence of creep of the supporting system on surrounding rocks has been ignored. (6) Compressive stress is defined as negative and tensile stress as positive. The calculation model is shown in Figure 2.

3. Viscoelastic-Plastic Analysis of Tunnel Surrounding Rock

3.1. Creep Analysis in Elastic Zone

The generalized Kelvin model is widely used in rock mass engineering because it can describe the elastic hysteresis of rock mass and reflect instantaneous deformation characteristics [3338]. Considering the time effect of stress redistribution in the plastic zone during its formation, this paper analyzed the creep characteristics of rock by using the generalized Kelvin model [39, 40]. The constitutive model of the generalized Kelvin model is

In equation (1): are parameters for the generalized Kelvin model; is a time-dependent differential operator, ; is the stress tensor; is a deviatoric tensor of strain, (x and y represent the horizontal and vertical directions in the plane coordinate system of Figure 2 and z represents the direction along the tunnel); represents mean strain; and represents mean stresses.

In the viscoelastic region, for a point of radius , is an arbitrary positive constant. The stress state only varies with the radius of the plastic region, but not with time. According to the elastic solution of the plane symmetry problem, equation (1) is deformed into

In equation (2), is the long-term shear modulus of surrounding rock, ; is slack time, ; represents radial stress; represents radial strain; represents hoop stress; and represents hoop strain.

Assuming that the volume strain of rock mass is zero during creep, where represents the radial displacement of surrounding rock.

When t = 0, there iswhere represents the displacement of surrounding rock when ; is the unified cohesion of surrounding rock in the elastic zone; is the unified internal friction angle of surrounding rock in the elastic zone.

Thus, the stress and displacement in the viscoelastic region can be obtained:

In equation (5),

In equations (5) and (6), is the radial stress in the viscoelastic region; is the hoop stress in the viscoelastic region; is the radial displacement of the viscoelastic region; and is the differential value between circumferential stress and radial stress of surrounding rock.

3.2. Creep Analysis in Plastic Zone

To study the load transfer law of supporting structure and surrounding rock more accurately and consider the influence of intermediate principal stress, this paper conducted the elastic-plastic analysis of tunnel surrounding rock by using unified strength theory [41]. Namely,where is the unified internal friction angle of geotechnical materials under different states, and is the unified cohesion of geotechnical materials under different states, where indicates the state of surrounding rock in elastic zone, fractured zone, and surrounding rock reinforcement body, respectively.

When ,where c and are the cohesive force and internal friction angle of surrounding rock, respectively. b is the influence parameter of intermediate principal stress, representing the coefficient of the influence degree of intermediate principal shear stress and normal stress on the corresponding plane on material failure, and the value range is .

When ,where is the cohesion of the rock mass in the loosened area, and is the internal friction angle of the rock mass in the loosened area.

When ,where is the cohesion of surrounding rock reinforcement body; can be taken; is the internal friction angle of surrounding rock reinforcement body; andwhere is the shear strength of the anchor; is the cross-sectional area of the anchor; and are the row spacing between the anchors.

If the bolt supporting force is regarded as a physical force, the equilibrium equation of the axisymmetric plane strain mechanics model in the surrounding rock reinforcement body is as follows:

According to the boundary bar when , there is , obtained by the simultaneous equations (7), (10), and (12). The stress distribution of surrounding rock reinforcement can be

The force acting on the bolt in the upper form is

where is the anchor force caused by the deformation of surrounding rock; is the pretension force of anchor; and the other symbols have the same meaning as before.

Among them,where is the elastic modulus of the bolt and are the displacements of the inside and outside ends of the bolt, respectively.

Obtained from equation (11), the stress conditions at the external boundary of surrounding rock reinforcement can be . Repeating the above analysis process and using equations (7) and (8), this paper can obtain the stress distribution of surrounding rock in the plastic zone of external reinforcement of surrounding rock:

On the elastic-plastic interface, () satisfying

Then, the radius of the plastic zone is

Substituting equation (18) into equation (16), obtainable stress on the elastic-plastic boundary is

Based on the assumption that the volume of the plastic zone is invariant, the radial displacement expression of the plastic zone is as follows:where is the elastic modulus of surrounding rock; is the shear modulus of surrounding rock; and is Poisson’s ratio of surrounding rock. The other symbols have the same meaning.

Transforming equation (20) with Laplace transformation, the following equations are obtained:

Then, the displacement of surrounding rock is obtained by Laplace reversal:

The supporting force provided by initial support iswhere is the displacement release of surrounding rock when bolt support is applied considering the spatial effect of tunnel excavation, which can be calculated by 60% of the elastic displacement of surrounding rock or 20% to 30% of the total displacement; is the stiffness coefficient of initial support; andwhere is the elastic modulus of shotcrete and is the thickness of shotcrete.

The supporting force provided by secondary lining iswhere is the radial displacement of surrounding rock before secondary lining construction. Assuming that the elastic-plastic deformation of surrounding rock is complete and the creep deformation of lining is not taken into account, the following equation can be obtained at this time:where is the stiffness coefficient of lining, andwhere is the shear modulus of secondary lining; is Poisson’s ratio of secondary lining; and , are the internal and external radius of the secondary lining.

Then, the stress distribution of surrounding rock reinforcement is

From equation (28), the force at the interface between reinforcement body and plastic zone () can be obtained at t time, and the radius expression of the plastic zone can be obtained:

4. Stress Form Analysis of Composite Lining

When using composite lining for tunnel support, it is generally required that the initial support should fully ensure the stability of surrounding rock during tunnel excavation, while the secondary lining mainly bears the additional stress caused by rock creep during service and provides appropriate safety reserve [42]. However, during the operation of the tunnel, the maximum radius of the surrounding rock fractured zone should not exceed the scope of the surrounding rock reinforcement body, to avoid the instability of the tunnel caused by the failure of bolts, which makes the weight load of the surrounding rock in the fractured zone attached to the secondary lining. Secondly, the initial support and the secondary lining during the creep process of the surrounding rock should be guaranteed to bear the additional stress of the surrounding rock creep, to ensure the surrounding rock and supporting structure of the tunnel during the operation are stable.

4.1. Stress Calculation of Support Structure

Assuming that the design life of the tunnel is T, substituting t = T into the theoretical deduction mentioned above, the radius of the plastic zone and displacement of surrounding rock and radial displacement of the external boundary of surrounding rock reinforcement during operation can be calculated. Then, the stress of the anchor bolt can be solved:

The force acting on shotcrete is as follows:

The force acting on the secondary lining is as follows:

4.2. Stability Checking Calculation of Supporting System

Under the creep action of surrounding rock during tunnel operation, the supporting force of the bolt gradually increases to a certain value. Suppose the supporting force is greater than the design value of the supporting strength of the bolt. In that case, the supportive effect will be lost due to the tension failure of the bolt and the instability and collapse of the surrounding rock in the fractured area. At this time, the secondary lining bears the additional force produced by rock creep and needs to resist the self-weight force of collapsed rock mass. It has adverse effects on tunnel stability. Therefore, the creep effect of surrounding rock should be considered, and the strength of the tunnel support structure should be checked.

4.2.1. Stability Check of Bolt Support System considering Creep of Surrounding Rock

To ensure the stability of the bolt support system [43], it should satisfywhere is the tension stress of a single bolt and is the design value of the tension strength of the bolt.

If the stress of the bolt does not satisfy equation (33) after checking calculation, the supporting effect of bolt should be strengthened or the strength of secondary lining should be increased to share the supporting force of the bolt.

4.2.2. Calculation of Lining Ultimate Bearing Capacity

The plastic failure of the inner edge under confining pressure is taken as the calculation condition of the limit supporting force of secondary lining. The inner and outer radius of the lining are and , respectively. acts on the outer edge of the secondary lining. The elastic stress solution of the secondary lining obtained from Lame’s solution is as follows:where represents the force on the outer edge of the secondary lining.

To obtain the ultimate bearing capacity of lining and make lining just enter the plastic state, that is, is an elastic-plastic interface at this time, according to the Mohr–Coulomb criterion, it can be concluded thatwhere is the ultimate bearing capacity of secondary lining; is the cohesion of secondary lining; is the angle of internal friction of secondary lining.

Then, the secondary lining supporting force derived from the theory in this paper should satisfywhere represents the maximum force on the secondary lining during the creep process.

Similarly, the ultimate bearing capacity of shotcrete can be obtained as follows:where is the ultimate bearing capacity of shotcrete; is the cohesive force of shotcrete; and is the internal friction angle of shotcrete.

The shotcrete supporting force should be satisfied:where is the maximum force on shotcrete during the creep process.

Equations (36) and (38) are used to evaluate the stability of lining. If the conditions mentioned above are not satisfied, measures such as increasing lining thickness or increasing lining label should be taken to ensure the stability of tunnel lining in operation.

The stress of composite lining support structure is calculated by equations (30)–(32) and checked by equations (33), (36), and (38). Finally, the stress characteristics under specific engineering geological conditions are determined, which provides a reference for the design of support structure and stability evaluation.

5. Engineering Example Analysis

The expressway from Shiyan city in Hubei province to Tianshui city in Gansu province is a horizontal link in the national expressway network design. It is an important component of Shaanxi’s expressway network that “connects east and west, connects north and south, covers the whole province, and reaches all neighbors.” The Mingyazi tunnel is located in Shaanxi province, which is the only specialty tunnel across the whole line. The length of the right line of the tunnel is 4980 m, and the length of the left line is 4944 m. It belongs to the superlong tunnel and is the control project of Shiyan–Tianshui Expressway construction. The site of the tunnel is shown in Figure 3 [44, 45].

This paper applied the theoretical results of this paper to analyze the force of supporting structure.

The equivalent radius of the tunnel is  m, the initial amplitude of Earth stress is P = 2.3 MPa, and the specific engineering geological parameters are shown in Table 1.

5.1. Force Law of Initial Support

As shown in Figure 4, given that other supporting conditions remain unchanged, the initial support should be designed in accordance with the principle that the joint force can guarantee the safety of the construction stage with the surrounding rock. The diameter of the bolt is calculated to be 18 mm. the spacing between rows is 1.5 m × 1.5 m, and the shotcrete thickness is 160 mm. That is to say, it meets the requirements of the above principles. At this time, the ultimate bearing capacity of shotcrete . With the creep development of surrounding rock, shotcrete stress increases gradually (as shown in Figure 3). When , the spray concrete has already reached the stress limit, after which the shotcrete will be destabilized and ultimately destroyed, leading to declining stability of surrounding rock reinforcement. This is not conducive to the long-term stability of the tunnel. Therefore, the thickness of shotcrete should be increased. Take , and according to the service life of 100 years to calculate and analyze (), it is concluded that the maximum force of shotcrete is , and at this time, the ultimate bearing capacity of concrete is . That is, shotcrete can provide support for the surrounding rock of the tunnel without damage during the whole operation period of the tunnel.

Figure 5 shows the variation law of bolt radial stress with time under different support conditions. When adopting bolt support scheme 1, bolt parameters take diameter , spacing between rows for initial bolt support. Although the stability of surrounding rock can be guaranteed during tunnel excavation, the stress of the anchor gradually increases when the surrounding rock enters the creep process. This eventually leads to the failure of anchor rod when the stress exceeds its ultimate bearing capacity, which has adverse effects on the stability of surrounding rock reinforcement. When adopting bolt support scheme 2 (bolt parameters take diameter , spacing between rows ) for initial bolt support, the stress level of bolt increases, indicating that the supporting effect of bolt increases. With the creep of surrounding rock, the maximum stress of bolt is , but the ultimate bearing capacity provided by bolts in this scheme is . Therefore, during the whole operation period of the tunnel, the anchor bolt can resist the additional stress caused by creep of surrounding rock without fracture failure.

5.2. Stress Law of Secondary Lining

As shown in Figure 6, with the creep of surrounding rock, the radial force of secondary lining increases, but the change rate decreases gradually. Other conditions remain unchanged. When the thickness of the secondary lining increases from 300 mm to 500 mm, the radial force of the secondary lining increases at the same time. The load-sharing ratio of the secondary lining (the ratio of secondary lining support force to initial support surrounding rock pressure) K increases with the creep of surrounding rock, indicating that the supporting effect of secondary lining increases gradually during the creep process to bear additional stress caused by the creep of surrounding rock. The ultimate load-sharing ratio of the two kinds of secondary lining thickness is 25.6% and 36.8%, indicating that increasing the thickness of secondary lining reduces the proportion of stress in surrounding rock and initial support and makes the stress of surrounding rock concentrate on the secondary lining.

As shown in Figure 7, when the initial support strength is relatively conservative (initial support scheme 1), the stress level and load-sharing ratio of secondary lining decrease during the creep process of surrounding rock. This is because the elastic modulus of surrounding rock reinforcement increases with the increase of initial support strength. Its load-bearing capacity increases with the rise in initial support strength. The process of coordinated deformation of the supporting structure and surrounding rock can share a larger proportion. Additional stress caused by the creep of surrounding rock results in the decrease of the load-sharing ratio of secondary lining [46].

6. Discussions

Based on the homogenization method, the concept of surrounding rock plus solid is proposed, and the stress variation law of primary support and secondary lining in the process of surrounding rock creep is obtained. The theoretical model adopted in this paper does not consider the influence of tunnel section shape, lateral pressure coefficient, and other factors on lining stress. The above problems will be further analyzed in the follow-up research.

7. Conclusions

(1)Based on the unified strength theory, this paper derived an equation that characterizes the relationship between stress and deformation by performing a viscoelastic-plastic analysis on the surrounding rock reinforcement body. In the model of this paper, both the bolt and the soil within its anchorage range are considered as the surrounding rock reinforcement body. The bolt force is regarded as the internal force that acts on the surrounding rock reinforcement body.(2)The example in this paper analysis reveals that, during the creep process of surrounding rock, the load-sharing ratio of secondary lining gradually increases, whereas its change rate gradually decreases. At the same time, the load-sharing ratio of secondary lining decreases with more initial support, while the opposite happens with more secondary lining support. Therefore, when designing composite lining, attention should be paid to the coordinated allocation of initial support and secondary lining support strength, optimizing their respective supporting role in the creep process of surrounding rock.(3)This paper proposes that when designing the composite lining in weak strata, the initial support should not only guarantee the tunnel safety during the construction stage but also cooperate with the secondary lining to brace the additional stress caused by the creep, ensuring that the supporting structure stays stable across the whole period of tunnel operation.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported in part by the China Postdoctoral Science Foundation (2019M663936XB) and the Scientific Research Program Funded by Shaanxi Provincial Education Department (19JK0399).