An Analytical Solution for the Frost Heaving Force considering the Freeze-Thaw Damage and Transversely Isotropic Characteristics of the Surrounding Rock in Cold-Region Tunnels

Zhang, Jiabing; Zhang, Xiaohu; Fu, Helin; Wu, Yimin; Huang, Zhen; Chen, Deyue

doi:https://doi.org/10.1155/2020/6654778

Advances in Civil Engineering

On this page

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

Special Issue

Advancements in Tunnel Construction

View this Special Issue

Research Article | Open Access

Volume 2020 | Article ID 6654778 | https://doi.org/10.1155/2020/6654778

An Analytical Solution for the Frost Heaving Force considering the Freeze-Thaw Damage and Transversely Isotropic Characteristics of the Surrounding Rock in Cold-Region Tunnels

Jiabing Zhang,^1,2Xiaohu Zhang,^1,2Helin Fu,^3,4Yimin Wu,^3,4Zhen Huang ,^1,2and Deyue Chen^1,2

Academic Editor: Ma Jianjun

Received02 Nov 2020

Revised17 Nov 2020

Accepted18 Nov 2020

Published03 Dec 2020

Abstract

Frost damage is a frequent occurrence in cold regions and can threaten the normal use and structural stability of tunnel engineering projects. To accurately determine the frost heaving force and effectively evaluate the frost damage in cold-region tunnels, an analytical solution for the frost heaving force considering the freeze-thaw (F-T) damage and transversely isotropic characteristics of surrounding rock is presented based on complex variable theory and the power series method. The calculation results indicate that the frost heaving force acts on the lining considering that the transversely isotropic characteristics of surrounding rock are significantly greater than those when assuming the surrounding rock is homogeneous isotropic media. This demonstrates that the transversely isotropic characteristics of surrounding rock have a considerable impact on the frost heaving force and should be considered. The frost heaving force continuously increases as the bedding angle increases from 0° to 90°, and the maximum frost heaving force in the Guanjiao tunnel (the rock mass bedding angle is 30°) of the Xining-Geermu Railway in China is approximately 1.04 MPa. In addition, the influence of F-T cycles on the frost heaving force in cold-region tunnels is investigated based on the analytical solution of the frost heaving force proposed in this paper. The frost heaving force acting on the lining decreases with an increasing number of F-T cycles due to the deterioration of the mechanical parameters of the surrounding rock.

1. Introduction

The frost heaving of rock is a critical problem that exists in many geotechnical engineering applications, such as tunnel and mining engineering in cold regions [1–3]. With the increasing investment in the construction of infrastructure in the Qinghai-Tibet Plateau in China, tunnel and underground project construction in cold regions has progressively increased [4, 5]. However, most tunnels during the operation period suffer from serious frost damage, as shown in Figure 1, thereby severely affecting the normal use and structural stability of tunnel engineering projects and even causing accidents [6, 7]. When the underground water of tunnel in cold region freezes into ice, the surrounding rock expands in volume, and the deformation pressure is frost heaving force, which is determined by the amount of frost heaving deformation, the stiffness of surrounding rock, and the stiffness of the lining structure. Under this action, frost damage phenomena such as lining cracking will occur [8]. Therefore, the effective calculation of the frost heaving force is significant for engineering applications, helping to provide technical support and guarantee structural safety and scientifically based construction in cold-region tunnels.

Many studies have been performed to obtain a reasonable frost heaving force in cold-region tunnels [9–12]. These studies demonstrated that the frost heaving force is primarily induced by the 9% volume expansion from the water-ice phase change in the surrounding rock [13–15]. In addition, a number of theoretical models have been developed in recent decades to predict the frost heaving forces in cold-region tunnels [16–20]. The frost heave of the surrounding rock is assumed to be isotropic in all of the above analytical solutions; however, very few publications are available where the frost heave of the surrounding rock is assumed to be transversely isotropic. In reality, the transversely isotropic characteristics of surrounding rock have a considerable impact on the frost heaving force [21].

In addition, many multifield coupled frost heave models have been proposed to effectively evaluate the frost damage in cold-region tunnels [22, 23]. Lai et al. [24] developed a hydrothermal-mechanical model to describe the freezing process of silty soil. Furthermore, Zhang et al. [25] advanced a hydrothermal-salt-mechanical model in freezing saline soil. However, to the best of the authors’ knowledge, F-T damage to rock masses is also a critical problem during tunnel engineering projects [26, 27], and the frost heaving force acting on the lining changes constantly after the recurrent action of freezing and thawing. Therefore, F-T damage to the surrounding rock should be considered in the frost heaving force calculation.

To accurately determine the frost heaving force in cold-region tunnels, an analytical solution for the frost heaving force acting on the lining that considers the F-T damage and the transversely isotropic characteristics of the surrounding rock has been first presented based on complex variable theory and the power series method. The mechanical parameters of surrounding rock with different F-T cycles were then determined by a series of experiments [28]. Finally, the effects of the number of F-T cycles on the frost heaving force in cold-region tunnels were discussed in detail based on the analytical solution of the frost heaving force proposed in this paper.

2. Analytical Solution for the Frost Heaving Force Acting on the Lining considering the F-T Damage and Transversely Isotropic Characteristics of the Surrounding Rock

2.1. Calculation Model and Basic Assumptions

The calculation model of the frost heaving force considering the F-T damage and transversely isotropic characteristics of the surrounding rock is shown in Figure 2(a). The calculation model includes three parts: the lining region (S₁), the freeze-thawed surrounding rock (S₂), and the unfrozen surrounding rock region (S₃). Among them, r₁ and r₂ are the inner and outer radii of the lining, respectively; r₃ is the radius of the freeze-thawed surrounding rock; β is the bedding angle of the surrounding rock; and σ_x^∞, σ_y^∞, and τ_xy^∞ are the stress components acting at infinity. To obtain the analytical solution of the frost heaving force in cold-region tunnels, the following assumptions are introduced: (1) the calculation model for the frost heaving force can be approximated as an infinite plane hole problem, which is not affected by the boundary effect, and the ground stress change gradient can be ignored; (2) the freeze-thawed surrounding rock and unfrozen surrounding rock are heterogeneous with transversely isotropic properties, while the lining is a homogeneous material with isotropic and elastic mechanical properties; (3) the calculation model for the frost heaving force can be regarded as a plane strain problem; and (4) the cross section of the tunnel is equivalent to a circle, and the frost heaving and F-T damage to the lining are not considered.

As shown in Figure 2, according to complex variable theory, the mapping function that transforms the outer region of the circular hole in the z-plane into the outer region of the unit circle in the ζ-plane can be given aswhere r₂ is the outer radius of the lining.

2.2. Mapping Function for the Elastic Mechanic Analysis of Transversely Isotropic Rocks

In anisotropic elastic mechanics, the cross section no longer remains flat after deformation due to the presence of positive shear and shear coupling [29]. To solve this problem, three mapping functions z = ω(ζ), z₁ = ω₁(ζ₁), and z₂ = ω₂(ζ₂) are introduced to transform the outer regions of the circular hole in the z-, z₁-, and z₂-planes into the outer regions of the unit circle in the ζ-, ζ₁-, and ζ₂-planes, respectively. The mathematical relationships between the three mapping functions z, z₁, and z₂ can be given aswhere , , , and .where the constants α₁, α₂, β₁, and β₂ can be determined according to the elastic constants of the anisotropic rock mass.

For the transversely isotropic plane strain problem, when the coordinate z-axis is perpendicular to one plane of elastic symmetry, the compatibility equation for Airy’s stress function without considering body forces can be given aswhere a_ij represents the elastic constants of the transversely isotropic material and , , , , , , , , , and ; E_x, E_y, and E_z represent the elastic moduli in the x-, y-, and z-axis directions, respectively; , , and represent Poisson’s ratios; and represents the shear modulus of the transversely isotropic material. For the transversely isotropic plane strain problem, as shown in Figure 2, the values of the elastic moduli, Poisson’s ratios, and shear moduli should satisfy the following relations: , , , , , and .

The solution to equation (4) is related to the roots of the following equation:where the four roots in equation (6) must be conjugate complex roots (i.e. , , , and ). Here, µ₁ and µ₂ are exactly the same as those expressed in equation (3). Therefore, the mapping functions z₁ and z₂ can be calculated by equation (2) if the complex roots µ₁ and µ₂ are known.

Substituting equation (1) into equation (2), we can obtain the following equations:where is an analytical function of ζ. However, its conjugate, , is not an analytical function of ζ [29]. Therefore, the mapping functions z₁ and z₂ in equation (7) are not analytical functions of ζ. To determine the relationships between ζ₁, ζ₂, and ζ and to determine ω₁(ζ₁) and ω₂(ζ₂), we assume that ζ₁ = ζ₂ = ζ = σ = e^iα along the unit circle, which corresponds to the boundary of the hole. Substituting this assumption into equation (7) and replacing σ with ζ₁ and ζ₂ in the expressions ω₁(ζ₁) and ω₂(ζ₂), respectively, we can obtain the following relationships between ζ₁, ζ₂, and ζ:

Therefore, ζ₁ and ζ₂ can be determined by equation (8) after the value of ζ is determined, and then the three mapping functions ω(ζ), ω₁(ζ₁), and ω₂(ζ₂) can be obtained.

2.3. Stress and Displacement Components of the Surrounding Rock and Lining

The complex function method is used to determine the stress and displacement components of the surrounding rock and lining, and the two complex potential functions Φ₁(z₁) and Ψ₁(z₂) in the unfrozen surrounding rock region (S₃) can be expanded by the Laurent series, as follows:where the complex constant , , and a_1k, a_2k, b_1k, and b_2k are undetermined coefficients. , , and can be determined according to the stress components acting at infinity (i.e., σ_x^∞, σ_y^∞, and ), as follows:

The two complex potential functions Φ₂(z₁) and Ψ₂(z₂) in the freeze-thawed surrounding rock region (S₂) can be expanded by the Laurent series, as follows:where the complex constants , , , and and c_1k, c_2k, d_1k, d_2k, e_1k, e_2k, f_1k, and f_2k are undetermined coefficients.

Similarly, two complex potential functions Φ₃(z) and Ψ₃(z) in the lining region (S₁) can be expanded by the Laurent series, as follows:where the complex constants , , , and and , , h_1k, h_2k, m_1k, m_2k, n_1k, and n_2k are undetermined coefficients.

When the lining is a homogeneous isotropic material, the stress components at any point on the z-plane can be given by the complex potential functions, as follows:where Re is the real part of a complex number and Im is the imaginary part of a complex number.

The freeze-thawed surrounding rock and unfrozen surrounding rock are heterogeneous with transversely isotropic materials; hence, the stress components at any point on the z-plane in the freeze-thawed and unfrozen surrounding rock regions can be determined by the complex potential functions, as follows:

Assuming that the volume of the lining and unfrozen surrounding rock are not affected by the change in the surrounding rock temperature, the displacement components of any point on the z-plane in the lining and unfrozen surrounding rock region can be determined by the complex potential functions, as follows:where G is the shear modulus of the lining. Since the longitudinal length of the tunnel is much larger than the diameter of the section, this problem can be regarded as a plane strain problem, and then . , , , and .

Since the volume of the freeze-thawed surrounding rock region is affected by temperature changes, the displacement components at any point on the z-plane in the freeze-thawed surrounding rock region cannot satisfy the expression of a complex variable function. After mathematical deduction, the displacement equation considering the transversely isotropic characteristics of surrounding rock can be obtained as follows:where ε represents the linear expansion coefficient of surrounding rock in the freeze-thawed surrounding rock region on the isotropic plane y′oz′, represents Poisson’s ratio perpendicular to the isotropic plane, and θ represents the volume change rate of surrounding rock in the freeze-thawed surrounding rock region.

2.4. Boundary Conditions and Solution of Frost Heaving Force

In terms of the stress and displacement solutions of elastic mechanics, boundary conditions are usually divided into two categories: displacement boundary conditions and stress boundary conditions. The complex potential function must satisfy the continuous boundary conditions of displacement and stress at interfaces ρ = ρ₂ and ρ = ρ₃ and the stress boundary conditions at interface ρ = ρ₁. Therefore, when frost heaving of rock occurs, the displacement boundary conditions and the stress boundary conditions at interfaces ρ = ρ₁, ρ = ρ₂, and ρ = ρ₃ are as follows.

2.4.1. Stress Boundary Conditions

2.4.2. Displacement Boundary Conditions

It should be noted that the expressions and in equation (4) represent the formation displacement component caused by the initial ground stress of the surrounding rock, and thus these expressions should be removed when calculating the formation displacement caused by tunnel excavation in equation (19b). In addition, it is assumed that there is no relative slip between the contact surfaces in equations (17) and (18). From the mapping function (1), we can obtain the following:

Based on the power series solution method proposed by Chen [30] and through the combination of the stress boundary conditions and displacement boundary conditions, infinite linear equations can be obtained. The coefficients in the complex potential functions can be obtained by solving the infinite linear equations, and the stress distributions of the surrounding rock and lining can be obtained by equations (13) and (14), respectively. Then, the stress components in the rectangular coordinate system are transformed into the polar coordinate system by equation (21) to obtain the radial stress on the outer boundary of the lining. Finally, the analytical solution of the frost heaving force in cold-region tunnels, which considers the transversely isotropic characteristics of surrounding rock, can be obtained by analysing the radial stress on the outer boundary of the lining before and after frost heaving.

3. A Case Study of the Frost Heaving Force considering the F-T Damage and Transversely Isotropic Characteristics of the Surrounding Rock

3.1. Determination of the Calculation Model

To study the frost heaving force considering the F-T damage and transversely isotropic characteristics of the surrounding rock, section DyK304 + 370 of the Guanjiao tunnel of the Xining-Geermu Railway in China is taken as a case study. The Guanjiao tunnel is 32.605 km in length and 200 m in depth. The weight of the surrounding rock is 22.0 kN/m³, the bedding angle is 30°, the lateral pressure coefficient is 1.0, and the volume change rate of the surrounding rock under the condition of frost heaving is 9%.

In the above theoretical analysis of the frost heaving force, the tunnel section is treated as a circular section, while the actual excavated section of the Guanjiao tunnel is a curved wall section. At present, a noncircular section is usually treated as a circular section by using mathematical transformation. The main treatment methods for circular sections are as follows:where r is the radius of the tunnel after equal circle treatment, b is the section span of the original tunnel, and h is the section height of the original tunnel. The radius of the Guanjiao tunnel after equal circular treatment is approximately 6.27 m according to equation (22). Combined with the initial support thickness of 30 cm and the maximum freezing depth of 2.99 m in the Guanjiao tunnel, the radii r₁ = 5.97 m, r₂ = 6.27 m, and r₃ = 9.26 m can be calculated for the z-plane in Section 2.1. Furthermore, the radii ρ₁ = 5.97 m, ρ₂ = 6.27 m, and ρ₃ = 9.26 m can be obtained by equation (1) for the ζ-plane.

3.2. Determination of the Mechanical Parameters of Surrounding Rock

Before the calculation of the frost heaving force acting on the lining considering the F-T damage and transversely isotropic characteristics of the surrounding rock, the basic mechanical parameters of the surrounding rock should be determined. In general, the determination of the mechanical parameters of the surrounding rock is relatively complicated, and the methods used are mainly test methods, engineering analogy methods, inverse analysis methods, etc. In this study, a series of uniaxial compression experiments with different numbers of F-T cycles (the temperature control curve for each F-T cycle is shown in Figure 3) was conducted to obtain the elastic parameters of surrounding rock [28]. The basic elastic parameter test results for different numbers of F-T cycles are listed in Table 1.where E₁ and μ₁ are the elastic modulus and Poisson’s ratio in the direction parallel to the bedding surface, respectively, and E₂, μ₂, and are the elastic modulus, Poisson’s ratio, and shear modulus in the direction perpendicular to the bedding surface, respectively.

Combined with the elastic parameter test results in Table 1, the elastic constants of the surrounding rock, a′_ij, can be obtained for different numbers of F-T cycles in the local coordinate system, as shown in Table 2. Then, the elastic parameters of the surrounding rock are substituted into equations (4) and (5), and the four complex roots of the compatible equation in the local coordinate system under different numbers of F-T cycles can be determined, as shown in Table 3. Finally, the four complex roots of the compatible equation in the global coordinate system under different numbers of F-T cycles can be obtained according to the principle of coordinate transformation (as shown in equation (23)), and the constants α₁, α₂, β₁, and β₂ in equation (3) can be determined.

3.3. Determination of the Mechanical Parameters of the Lining

The initial support considers the action of only the I-beam support and C25 shotcrete. Among them, the steel arch support parameters are equivalently converted by improving the elastic modulus of the shotcrete. The specific conversion equation is as follows:where E is the elastic modulus of the converted shotcrete, E_c is the elastic modulus of the original shotcrete, E_s is the elastic modulus of the steel arch, A_s is the sectional area of the steel arch, h is the thickness of the original shotcrete, and d is the spacing between steel arches. According to equation (24), the tunnel’s lining elasticity modulus is 34.27 GPa, Poisson’s ratio is 0.25, and the shear modulus is 13.71 GPa.

3.4. Calculation and Analysis of the Frost Heaving Force in the Guanjiao Tunnel

With the zero freezing and thawing cycles as an example, the analytical solution of the stress distribution in a transversely isotropic rock mass tunnel can first be obtained according to the boundary conditions in equations (18) and (19). Then, the stress components in the rectangular coordinate system are transformed to a polar coordinate system by equation (21), and the radial stress on the outer boundary of the lining can be obtained. Finally, the frost heaving force considering the transversely isotropic frost heave of surrounding rock can be obtained by comparing and analysing the calculated results for analytical unfrozen heaving and analytical frozen heaving, as shown in Figure 4. The following conclusions can be obtained:(1)Figure 4(a) is the calculation model of the frost heaving force acting on the lining assuming that the surrounding rock is homogeneous isotropic continuous media. The radial stress on the outer boundary of the lining considering the transversely isotropic frost heave of the surrounding rock is significantly greater than that when assuming the surrounding rock is homogeneous isotropic media, and then it gradually increases with increasing bedding orientation of the surrounding rock.(2)Based on the complex variable theory and power series solution method, the frost heaving force considering the transversely isotropic frost heave of the surrounding rock is calculated to be between 0.21 MPa and 3.59 MPa (0 F-T cycles), and the frost heaving force continuously increases as the bedding angle increases from 0° to 90°. When the bedding orientation is 90°, the frost heaving force reaches a maximum. In addition, the rock mass bedding angle of the Guanjiao tunnel is 30°, and the maximum frost heaving force acting on the lining is approximately 1.04 MPa.(3)The positions of the maximum frost heaving force acting on the lining differ because the mechanical properties of transversely isotropic rocks exhibit distinct anisotropy and directional dependence. When the surrounding rocks are not layered, the frost heaving force acting on the lining is distributed symmetrically. The frost heaving force is the highest at the arch waist (θ = 0° or 180°) of the tunnel and the lowest at the arch top (θ = 90°) or arch bottom (θ = 270°) of the tunnel (as shown in Figure 4(a)). When the rock mass bedding angle is 0°–15°, the maximum frost heaving force occurs near the arch waist of the tunnel (θ = 0° or 180°) (as shown in Figures 4(b)-4(c)). When the rock mass bedding angle is 45°–60°, the maximum frost heaving force occurs near the right arch shoulder (θ = 40°) or the left arch foot (θ = 220°) (as shown in Figures 4(e)-4(f)). When the rock mass bedding angle is 75°–90°, the maximum frost heaving force occurs near the arch top (θ = 90°) or the arch bottom (θ = 270°) (as shown in Figures 4(g)-4(h)).(4)Figure 4 shows that the frost heaving force acting on the lining considering the transversely isotropic characteristics of surrounding rock is significantly greater than that assuming that the surrounding rock is homogeneous isotropic continuous media. It can be concluded that the anisotropy of the rock mass has a significant influence on the frost heaving force. Additionally, the lining structures of transversely isotropic rock mass tunnels are more prone to frost damage, which makes the maintenance of tunnels very difficult and poses a serious threat to the structure and operational safety of the tunnels. Therefore, the transversely isotropic characteristics of surrounding rock have a considerable impact on the frost heaving force and should be considered.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5 shows the circumferential stress on the outer boundary of the lining for different bedding angles of the surrounding rock. The circumferential stress on the outer boundary of the lining considering the transversely isotropic characteristics of the surrounding rock is significantly smaller than that when assuming that the surrounding rock is homogeneous isotropic media, and then it gradually decreases with increasing bedding orientation of the surrounding rock. Additionally, the circumferential stress on the outer boundary of the lining significantly increases because of the influence of frost heaving on the rock. When the surrounding rock is not layered, the maximum increment of the outer boundary circumferential stress of the lining is approximately 3.48 MPa (as shown in Figure 5(a)). Moreover, the increment of the outer boundary circumferential stress of the lining decreases with increasing bedding orientation of the surrounding rock (as shown in Figures 5(b)–5(h)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

4. Effect of F-T Cycles on the Frost Heaving Force in Cold-Region Tunnels

To comprehensively understand the effect of F-T cycles on the frost heaving force in cold-region tunnels, different mechanical parameters were used for the surrounding rock after different numbers of F-T cycles were applied. The maximum frost heaving force acting on the lining after F-T cycles can be obtained based on the analytical solution for the frost heaving force proposed in this paper. As revealed in Figure 6, the maximum frost heaving force acting on the lining decreases with an increasing number of F-T cycles when the bedding angle of the surrounding rock is constant. For example, when the bedding angle is 30°, the maximum frost heaving forces acting on the lining are 1.04 MPa, 0.98 MPa, 0.94 MPa, 0.91 MPa, 0.89 MPa, and 0.85 MPa for 0, 5, 10, 15, 20, and 30 F-T cycles, respectively. The mechanism for this phenomenon is that the deformation of rock is more likely to occur after the action of external forces due to the rock deterioration gradually increasing as the number of F-T cycles increases [13]. This cushions the lining structure’s stress and results in a gradual reduction in the frost heaving force in cold-region tunnels.

(a)

(b)

(c)

(d)

(e)

(f)

5. Discussion

The seasonal frozen and permafrost regions are widely distributed throughout China, accounting for 75.8% of the total land area. Along with the implementation of the strategy of “the Belt and Road,” a large number of tunnels are under construction in cold regions. The frost heaving force is one of the most critical challenges that can cause serious frost damage in cold-region tunnels. Therefore, a reasonable calculation of the frost heaving force that considers the F-T damage and the transversely isotropic characteristics of the surrounding rock is expected to provide references for the design and construction of cold-region tunnels.

At present, the commonly used calculation models of the frost heaving force are mainly divided into three categories. The first model assumes that the surrounding rock forms a freezing-thawing circle and the pore water or fissure water in the surrounding rock of the freezing-thawing circle will freeze and expand under a negative temperature condition, thus generating a frost heaving force. The second model suggests that there is a water accumulation space between the initial support and the secondary lining of the tunnel and the pore or fissure water in the surrounding rock will continuously supply to the water accumulation space. Under a negative temperature condition, the water will freeze and expand in volume, and the existence of the primary support and secondary lining restricts such volume expansion, thus forming the frost heaving force. The third model holds that the surrounding rock of any tunnel all has a weathering layer and the water in the weathering layer will freeze and expand under the condition of negative temperature, thus causing the frost heaving force. However, due to the complexity of the tunnel and underground engineering system, the existing calculation models of the frost heaving force are all studied with isotropic surrounding rocks as a precondition. In fact, as the tunnel construction gradually moves to cold regions, the surrounding rock of the tunnel will be in a state of alternating positive and negative temperatures for a long time, leading to a stronger anisotropy of the rock mass. In this case, the anisotropy of the rock mass cannot be ignored [4]. Therefore, an analytical solution for the frost heaving force in cold-region tunnels considering the transversely isotropic characteristics of the surrounding rock, which is presented in this paper, may improve the design techniques and safe operation of a transversely isotropic rock mass tunnel in cold regions.

In addition, the tunnel section is treated as an equal circle in the calculation of the frost heaving force in cold-region tunnels, and the model corresponding to the actual tunnel section is not used for the calculation. Due to the difference in the models, these research results are somewhat different from a realistic situation, so it is still necessary to further study the problem of the frost heaving force considering an actual section of tunnel in future studies.

6. Summary and Conclusions

To accurately determine the frost heaving force acting on the lining in cold-region tunnels, an elastic analytical solution of the frost heaving force in cold-region tunnels that considers the F-T damage and transversely isotropic characteristics of the surrounding rock has been developed. The following conclusions can be drawn based on the limited research work reported in this paper:(1)The frost heaving force acting on the lining continuously increases as the bedding angle increases from 0° to 90° when the number of F-T cycles is constant. The rock mass bedding angle of the Guanjiao tunnel is 30°, and the maximum frost heaving force acting on the lining is approximately 1.04 MPa.(2)The frost heaving force acting on the lining considering the transversely isotropic characteristics of surrounding rock is significantly greater than that when assuming the surrounding rock is homogeneous isotropic continuous media. Therefore, the transversely isotropic characteristics of surrounding rock have a considerable impact on the frost heaving force and should be considered.(3)The positions of the maximum frost heaving force acting on the lining differ because the mechanical properties of transversely isotropic rocks exhibit distinct anisotropy. Additionally, the frost heaving force acting on the lining decreases with an increasing number of F-T cycles due to the deterioration of the mechanical parameters of the surrounding rock.

Data Availability

All data are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was financially supported by the National Natural Science Foundation of China (grant nos. 51978668 and 51578550), the National Science Foundation for Distinguished Young Scholars of China (grant no. 51508181), and the Guangxi University Young and Middle-Aged Teachers’ Basic Scientific Research Ability Improvement Project (2020ky01011).

References

Y. Lai, W. Hui, W. Ziwang, L. Songyu, and D. Xuejun, “Analytical viscoelastic solution for frost force in cold-region tunnels,” Cold Regions Science and Technology, vol. 31, no. 3, pp. 227–234, 2000.
View at: Publisher Site | Google Scholar
P. Wang and G. Zhou, “Frost-heaving pressure in geotechnical engineering materials during freezing process,” International Journal of Mining Science and Technology, vol. 28, no. 2, pp. 287–296, 2018.
View at: Publisher Site | Google Scholar
X. Zhang, Y. Lai, W. Yu, and Y. Wu, “Forecast analysis for the re-frozen of Feng Huoshan permafrost tunnel on Qing-Zang railway,” Tunnelling and Underground Space Technology, vol. 19, no. 1, pp. 45–56, 2004.
View at: Publisher Site | Google Scholar
J. Zhang, H. Fu, Z. Huang, Y. Wu, W. Chen, and Y. Shi, “Experimental study on the tensile strength and failure characteristics of transversely isotropic rocks after freeze-thaw cycles,” Cold Regions Science and Technology, vol. 163, pp. 68–77, 2019.
View at: Publisher Site | Google Scholar
X. Zhao, X. H. Yang, H. W. Zhang, H. P. Lai, and X. Y. Wang, “An analytical solution for frost heave force by the multifactor of coupled heat and moisture transfer in cold–region tunnels,” Cold Regions Science and Technology, vol. 175, Article ID 103077, 2020.
View at: Publisher Site | Google Scholar
Y. M. Wu, P. Xu, W. B. Li, Z. J. Wang, Z. Y. Cai, and S. Shao, “Distribution rules and key features for the lining surface temperature of road tunnels in cold regions,” Cold Regions Science and Technology, vol. 172, Article ID 102979, 2020.
View at: Publisher Site | Google Scholar
C. Xia, Z. Lv, Q. Li, J. Huang, and X. Bai, “Transversely isotropic frost heave of saturated rock under unidirectional freezing condition and induced frost heaving force in cold region tunnels,” Cold Regions Science and Technology, vol. 152, pp. 48–58, 2018.
View at: Publisher Site | Google Scholar
C. C. Xia, Z. T. Lv, and Y. S. Wang, “Research progress and consideration on calculating method of frost heaving force of tunnels in cold region,” China Journal of Highway and Transport, vol. 33, no. 5, pp. 35–43, 2020.
View at: Google Scholar
S. B. Huang, Y. H. Ye, X. Z. Cui, A. P. Cheng, and G. F. Liu, “Theoretical and experimental study of the frost heaving characteristics of the saturated sandstone under low temperature,” Cold Regions Science Technology, vol. 174, Article ID 103036, 2020.
View at: Publisher Site | Google Scholar
H. Liu, X. Yuan, and T. Xie, “A damage model for frost heaving pressure in circular rock tunnel under freezing-thawing cycles,” Tunnelling and Underground Space Technology, vol. 83, pp. 401–408, 2019.
View at: Publisher Site | Google Scholar
Z. Lv, C. Xia, and Q. Li, “Experimental and numerical study on frost heave of saturated rock under uniform freezing conditions,” Journal of Geophysics and Engineering, vol. 15, no. 2, pp. 593–612, 2018.
View at: Publisher Site | Google Scholar
P. Yang, J.-m. Ke, J. G. Wang, Y. K. Chow, and F.-b. Zhu, “Numerical simulation of frost heave with coupled water freezing, temperature and stress fields in tunnel excavation,” Computers and Geotechnics, vol. 33, no. 6-7, pp. 330–340, 2006.
View at: Publisher Site | Google Scholar
Q. Liu, S. Huang, Y. Kang, and X. Liu, “A prediction model for uniaxial compressive strength of deteriorated rocks due to freeze-thaw,” Cold Regions Science and Technology, vol. 120, pp. 96–107, 2015.
View at: Publisher Site | Google Scholar
X. Tan, W. Chen, H. Liu, L. Wang, W. Ma, and A. H. C. Chan, “A unified model for frost heave pressure in the rock with a penny-shaped fracture during freezing,” Cold Regions Science and Technology, vol. 153, pp. 1–9, 2018.
View at: Publisher Site | Google Scholar
Z. Zhou, W. Ma, S. Zhang, Y. Mu, and G. Li, “Effect of freeze-thaw cycles in mechanical behaviors of frozen loess,” Cold Regions Science and Technology, vol. 146, pp. 9–18, 2018.
View at: Publisher Site | Google Scholar
Q. Feng, B.-S. Jiang, Q. Zhang, and L.-P. Wang, “Analytical elasto-plastic solution for stress and deformation of surrounding rock in cold region tunnels,” Cold Regions Science and Technology, vol. 108, pp. 59–68, 2014.
View at: Publisher Site | Google Scholar
G. Y. Gao, Q. S. Chen, Q. S. Zhang, and G. Q. Chen, “Analytical elasto-plastic solution for stress and plastic zone of surrounding rock in cold region tunnels,” Cold Regions Science and Technology, vol. 72, pp. 50–57, 2012.
View at: Publisher Site | Google Scholar
J. H. Huang, C. C. Xia, C. L. Han, and M. Huang, “The analytical solution of frost heave force acting on cold–region tunnel liner considering the anisotropy frost heaving of the surrounding rock,” Chinese Journal of Rock Mechanics and Engineering, vol. 34, no. S2, pp. 3764–3774, 2015.
View at: Google Scholar
I.-M. Lee and S.-W. Nam, “Effect of tunnel advance rate on seepage forces acting on the underwater tunnel face,” Tunnelling and Underground Space Technology, vol. 19, no. 3, pp. 273–281, 2004.
View at: Publisher Site | Google Scholar
Z. Lv, C. Xia, Q. Li, and Z. Si, “Empirical frost heave model for saturated rock under uniform and unidirectional freezing conditions,” Rock Mechanics and Rock Engineering, vol. 52, no. 3, pp. 955–963, 2019.
View at: Publisher Site | Google Scholar
Z. Lv, C. Xia, Y. Wang, and J. Luo, “Analytical elasto-plastic solution of frost heaving force in cold region tunnels considering transversely isotropic frost heave of surrounding rock,” Cold Regions Science and Technology, vol. 163, pp. 87–97, 2019b.
View at: Publisher Site | Google Scholar
Y. Kang, Q. Liu, and S. Huang, “A fully coupled thermo-hydro-mechanical model for rock mass under freezing/thawing condition,” Cold Regions Science and Technology, vol. 95, pp. 19–26, 2013.
View at: Publisher Site | Google Scholar
E. L. Liu and Y. M. Lai, “Thermo–poromechanics–based viscoplastic damage constitutive model for saturated frozen soil,” International Journal of Plasticity, vol. 128, Article ID 102683, 2020.
View at: Publisher Site | Google Scholar
Y. Lai, W. Pei, M. Zhang, and J. Zhou, “Study on theory model of hydro-thermal-mechanical interaction process in saturated freezing silty soil,” International Journal of Heat and Mass Transfer, vol. 78, pp. 805–819, 2014.
View at: Publisher Site | Google Scholar
J. Zhang, Y. M. Lai, J. F. Li, and Y. H. Zhao, “Study on the influence of hydro–thermal–salt–mechanical interaction in saturated frozen sulfate saline soil based on crystallization kinetics,” International Journal of Heat and Mass Transfer, vol. 146, Article ID 118868, 2020.
View at: Publisher Site | Google Scholar
H. Fu, J. Zhang, Z. Huang, Y. Shi, and W. Chen, “A statistical model for predicting the triaxial compressive strength of transversely isotropic rocks subjected to freeze-thaw cycling,” Cold Regions Science and Technology, vol. 145, pp. 237–248, 2018.
View at: Publisher Site | Google Scholar
C. M. Grossi, P. Brimblecombe, and I. Harris, “Predicting long term freeze-thaw risks on Europe built heritage and archaeological sites in a changing climate,” Science of the Total Environment, vol. 377, no. 2-3, pp. 273–281, 2007.
View at: Publisher Site | Google Scholar
H. L. Fu, J. B. Zhang, Z. Huang, H. W. Huang, and Y. Shi, “Experimental study on elastic parameters and uniaxial compressive strength of slate under freeze-thaw cycles,” Rock Soil Mechanics, vol. 38, no. 8, pp. 2203–2213, 2017.
View at: Google Scholar
A. Z. Lu, N. Zhang, X. L. Zhang, D. H. Lu, and W. S. Li, “Analytic method of stress analysis for an orthotropic rock mass with an arbitrary shaped tunnel,” International Journal of Geomechanics, vol. 15, no. 4, Article ID 04014068, 2015.
View at: Publisher Site | Google Scholar
Y. Z. Chen, “Multiple circular hole problem for an elastic half-plane,” Computers & Structures, vol. 52, no. 6, pp. 1277–1281, 1994.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Jiabing Zhang 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

Views

469

Downloads

1059

Citations