Study on Stress and Displacement of Axisymmetric Circular Loess Tunnel Surrounding Rock Based on Joint Strength

: The development of an effective evaluation method suitable for loess-tunnel excavation is necessary to avoid the collapse accidents caused by tunnel excavation and any secondary disasters. Although the Fenner formulas and the modiﬁed Fenner formulas are widely used in tunnel engineering, a defect still exists in these formulas because the Mohr–Coulomb (M–C) criterion exaggerates the tensile strength of the surrounding rock of the loess tunnel. A newly modiﬁed Fenner formula was derived based on joint strength to overcome this deﬁciency. First, the expressions of stress and the radius of the plastic zone of the surrounding rock of the loess tunnel and the expressions of radial displacement were derived based on the stress-equilibrium equation of the axisymmetric plane and the joint strength. Then, the difference in the modiﬁed Fenner formulas based on the two kinds of strength criteria for the loess tunnel were compared. The results showed that the radius of the plastic zone and the radial displacement of the loess tunnel determined by the modiﬁed Fenner formula based on joint strength were larger than those determined by the modiﬁed Fenner formula based on M–C strength. However, the plastic stress of the plastic zone determined by the modiﬁed Fenner formula based on joint strength was smaller. The comparative analysis reveals that the modiﬁed Fenner formula based on joint strength can evaluate the stress and plastic-displacement ﬁeld of the surrounding rock of a loess tunnel more reasonably.


Introduction
Loess is a multiphase porous medium with a special structural behavior.With the rapid development of the transportation infrastructure in China over the last few decades, many new tunnels have been constructed or are being constructed through regions with challenging geological conditions.Hundreds of tunnels with large cross-sections have been constructed in the ground, such as in the loess area in China [1,2].Several engineering facilities have been produced, such as highways, dams, bridges, and tunnels, to expand the living space of mankind.Jefferson et al. updated a series of research topics proposed by G.A. Mavlyanov concerning engineering geology and discussed key concepts [3].A collapse accident caused by tunnel excavation can cause ground collapse, which may harm the ground-water systems and ecological environments.Based on a probabilistic background, Oggeri et al. proposed a global retrospective analysis method which can define the unknown parameters of the rock and soil in the construction stage, so as to obtain more accurate results [4].The stress field in the original stratum can change during tunnel excavation, which can lead to stress redistribution and the formation of a plastic zone or even a fracture zone in the rock surrounding the tunnel [5].The stress and displacement field in the surrounding rock of the tunnel can be used to evaluate the stability of the surrounding rock and the design of the lining.Thus, accurately determining the radius of the plastic zone based on different strength criteria has great engineering importance in optimising the design of underground support structures and ensuring the construction safety of tunnels.
The Mohr-Coulomb (M-C) strength criterion is commonly used in geotechnical materials and still plays an important role at present.The damaged zone of rock masses in some deep circular tunnels was analysed in view of the strain-softening-based plastic zone using the M-C strength criterion [6].However, the Drucker-Prager (D-P) strength criterion was put forward to overcome the corner defect in M-C strength in π plane; it is widely used in numerical analysis of geotechnical engineering.Multiphase-field modelling of crack propagation in geological materials and porous media was analysed based on D-P strength criterion [7].The king of theoretical solutions in tunnel-excavation processes was deduced through the cavity-expansion theory based on the D-P strength criterion [8].The Hoek-Brown (H-B) strength criterion was developed in 1980 to reasonably consider the tensile and compressive strengths of rock mass [9], and has been widely used due to its adequacy and convenience in rock-behaviour predictions and applications in rock engineering [10].Over the past few decades, the analysis of tunnel mechanical properties has been mainly based on the M-C, D-P and H-B strength criteria to determine the plastic zone and its elastic-plastic stress for the evaluation of the stress environment and displacement field in the surrounding rock during tunnel excavation.
The Fenner formula was proposed on the basis of the M-C strength criterion [11], and the expressions of stress and the radius of the plastic zone in a circular tunnel under the condition of an infinite homogeneous stratum were derived.However, the radial stress σ R at the interface between the elastic and plastic zones in the surrounding rock of the tunnel was determined.Fenner adopted a simplified formula by neglecting the cohesion c in the place.On the basis of the M-C strength criterion, Kastner's formula of plastic stress and the radius of the plastic zone in a circular tunnel were derived for a circular tunnel under the condition of infinite homogeneous stratum by using the unconfined compressive strength σ c and the internal friction angle ϕ as parameters [12].The cohesion c at the interface between the elastic and plastic zones in the surrounding rock of a tunnel could be expressed by the unconfined compressive strength σ c and internal-friction angle ϕ.Therefore, the cohesion c at the interface between the elastic and plastic zones was appropriately considered in the Kastner formula.By comparing the difference between the Fenner and Kastner formulas, a modified Fenner formula of elastic-plastic stress and the plastic radius in a circular tunnel, which considered the formula at the interface between the elastic and plastic zones, was derived on the basis of M-C strength criterion [13].As proven, the modified Fenner formula based on the M-C strength criterion was equivalent to the Kastner formula.Fenner, Kastner and modified Fenner formulas are widely used in tunnel engineering and are mainly used to determine the plastic zone of a tunnel and to analyse the radial displacement of a tunnel.
The stress and displacement of surrounding rock of the tunnel were solved by Park, and on the basis of M-C criterion and the large strain similarity solution, the effect of small and large strain solutions on displacement was studied [14].Ren et al. studied the effect of the hydraulic pressure on the direction of the principal stress in plastic zone of surrounding rock of the circular rock tunnel; then, they analysed the relationship between the hydraulic pressure of tunnel wall and the plastic zone of surrounding rock of the tunnel based on Fenner formula [15].On the basis of the M-C strength criterion, Li et al. deduced the analytical solutions of stress displacement and the range of the plastic zone in the elastic, elastoplastic and rheological stages [16].Liu et al. presented a semi-analytical solution for the tunnel excavation problem based on an extended D-P model [17].Some elasticplastic studies were conducted on the basis of the H-B strength criterion to better analyse the elastic-plastic mechanical properties of a rock tunnel.For the rock tunnel, Pan et al. considered the structural plane of the surrounding rock, performed an elastic-plastic analysis of the surrounding rock of a circular tunnel and obtained the Fenner formula based on the H-B strength criterion [18].Zou et al. studied the theoretical solutions of stress, displacement and the plastic radius of a circular tunnel based on the H-B strength criterion and the non-associated flow law of elastic-brittle-plastic rock [19].Wang et al. believed that the intermediate principal stress was an important factor to determine the plastic zone; thus, the H-B strength criterion was modified by considering the intermediate principal stress, and the theoretical formula of the plastic zone radius of a rock tunnel was proposed [20].Ranjbarnia et al. developed a new analytical-numerical procedure based on the H-B strength criterion to provide the stresses and strains around a circular tunnel in rock masses exhibiting different stress-strain behaviours; a simple relationship was also presented to determine the closure of a circular tunnel in terms of the rock mass strength and tunnel depth [21].Sharan compared the applicability of the M-C and H-B strength criteria in a rock tunnel, analysed the elastic-brittle-plastic behaviour of tunnel and evaluated the stress difference of the surrounding rock in a tunnel between the expansive and non-expansive rocks using a numerical method [22].On the basis of the difference in the analytical solutions of displacement of the surrounding rock between the M-C and H-B strength criteria, Park et al. used the non-associated flow rule to compare the displacement difference of the surrounding rock in a tunnel between soft and hard rocks with the effect of the angle of dilatancy [23].Serrano et al. deduced a mathematical expression to determine the radius of the plastic zone based on M-C and H-B strength criteria; then, they compared the convergence of circular tunnels in elastoplastic rock masses [24].
The abovementioned considerable literature has promoted research on the elasticplastic mechanical properties of tunnels.However, two problems still need to be solved: one is that the defect still exists in the Fenner formula and in the modified Fenner formula derived based on M-C strength criterion because this criterion exaggerates the tensile strength of surrounding rock of the tunnel; the other is that the H-B strength criterion is unsuitable for the analysis of a tunnel in earth, especially a loess tunnel, given that it is mainly based on the empirical criterion of rock failure, which is put forward by the curve-failure envelope of a large number of rock tests.Xie et al. studied the effects of matrix suction and net confining pressure on the deformation and strength characteristics of compacted loess by carrying out the triaxial consolidation shear test of unsaturated loess, and the results showed that under certain matrix suction conditions, the shear strength of compacted loess increased with the increase of net confining pressure, and under constant net confining pressure, the shear strength increased with the increase of matrix suction [25].When the strength of structural loess is higher, it not only has higher shear strength but also has relatively higher tensile strength, and tensile strength considerably affects the mechanical properties of structural loess [26].Therefore, the influence of tensile strength of structural loess should be considered when the stress and displacement of the surrounding rock of structural loess tunnel need to be studied.Li et al. put forward the joint strength formula of structural loess to reasonably consider the characteristics of tensile strength of structural loess; it overcame the shortcoming that M-C strength exaggerated the tensile strength of loess [27].Through the formula and calculation of loess earth pressures based on joint strength, it can be verified that the joint strength formula can not only can reflect the characteristics of tensile strength of loess but also can study the characteristics of the shear strength of the loess [28].
Therefore, accurately determining the radius of plastic zone of loess tunnel based on joint strength criteria is an issue to be solved in loess tunnel.A newly modified Fenner formula for the loess tunnel will be derived on the basis of the joint strength theory that can consider the tensile and shear properties of loess.Then, the change in stress environment in the surrounding rock of a loess tunnel under the action of the average stratum stress should be studied.Lastly, the differences in the stress distribution, the radius of plastic zone and the radial displacement will be analysed and compared on the basis of the joint strength and M-C strength.

The Basis of Derivation 2.1. Basic Hypothesis for Derivation
The deep axisymmetric tunnel may be described as the following plane-strain problem in polar coordinates (r, θ): a tunnel of circular cross-section with radius a is excavated in a homogeneous isotropic rock mass of infinite extent, which is subjected to a uniform in situ stress σ s in all directions, and a radial support pressure σ a is applied on the tunnel wall (Figure 1).ppl.Sci.2023, 13, x FOR PEER REVIEW

Basic Hypothesis for Derivation
The deep axisymmetric tunnel may be described as the following plan lem in polar coordinates (r, θ): a tunnel of circular cross-section with radius in a homogeneous isotropic rock mass of infinite extent, which is subjected in situ stress σs in all directions, and a radial support pressure σa is applied wall (Figure 1).In Figure 1, a uniform in situ stress σs = γh (γ is weight per unit volum central depth of tunnel), r is vector diameter in the surrounding rock, R is ra zone, σr is radial stress, R σ is radial stress at the interface between the elas zones and θ σ is tangential stress.
Therefore, the stress equilibrium equation of the axisymmetric plane p lar coordinates can be expressed as [5].

Joint Strength of Loess
The loess area is generally a semi-arid area in the middle tempera groundwater level is deep, the water content of the loess is low, and the li the natural loess is generally 0.19 [29].Natural loess is in a hard plastic st hardness and high strength.In the process of tunnel excavation, the surrou the loess shows a certain degree of self-stability.After tunnel excavation, th surrounding rock is redistributed.In the tensile shear-stress zone, the tens shear stress of a loess with high strength do not exceed the limit.
The M-C strength line of loess can be determined according to conven shear strength test and the three failure stress circles, as shown in Figure 2 normal stress and τ is shear stress.By determining the point of interse strength line and horizontal axis, the value of tensile strength correspon strength can be obtained.In Figure 1, a uniform in situ stress σ s = γh (γ is weight per unit volume and h is the central depth of tunnel), r is vector diameter in the surrounding rock, R is radius of plastic zone, σ r is radial stress, σ R is radial stress at the interface between the elastic and plastic zones and σ θ is tangential stress.
Therefore, the stress equilibrium equation of the axisymmetric plane problem in polar coordinates can be expressed as [5].

Joint Strength of Loess
The loess area is generally a semi-arid area in the middle temperate zone.The groundwater level is deep, the water content of the loess is low, and the liquid index of the natural loess is generally 0.19 [29].Natural loess is in a hard plastic state, with high hardness and high strength.In the process of tunnel excavation, the surrounding rock of the loess shows a certain degree of self-stability.After tunnel excavation, the stress of the surrounding rock is redistributed.In the tensile shear-stress zone, the tensile stress and shear stress of a loess with high strength do not exceed the limit.
The M-C strength line of loess can be determined according to conventional triaxial shear strength test and the three failure stress circles, as shown in Figure 2, where σ is normal stress and τ is shear stress.By determining the point of intersection of M-C strength line and horizontal axis, the value of tensile strength corresponding to M-C strength can be obtained.
In order to unite the tensile strength and shear strength of loess, the strength that can take the tensile characteristics and shear characteristics into consideration is called the joint strength, and the corresponding failure strength curve is expressed in the form of hyperbola.Here "joint" means a kind of united consideration.In order to unite the tensile strength and shear strength of loess, the str take the tensile characteristics and shear characteristics into consideratio joint strength, and the corresponding failure strength curve is expressed hyperbola.Here "joint" means a kind of united consideration.
To fit a hyperbolic curve of the joint strength line of structural loess on Li et al. assumed that the horizontal axis intercept of the curve represe strength σt, and that the M-C strength envelope was regarded as the asymp the joint strength formula corresponding to the dotted line, which can desc and shear strengths of structural loess comprehensively and reasonably, w the basis of the tensile and conventional triaxial shear strength tests.

Equivalent Equation in the Form of Principal Stress of the Joint Strength
The ultimate equilibrium state of the joint strength on the σ-τ plane is ure 3, where σ1 is major principal stress and σ3 is minor principal stress.To fit a hyperbolic curve of the joint strength line of structural loess on the σ-τ plane, Li et al. assumed that the horizontal axis intercept of the curve represents the tensile strength σ t , and that the M-C strength envelope was regarded as the asymptote [26].Then, the joint strength formula corresponding to the dotted line, which can describe the tensile and shear strengths of structural loess comprehensively and reasonably, was derived on the basis of the tensile and conventional triaxial shear strength tests.
Figure 2 shows that the tensile strength determined by the M-C strength is more than twice that of the tensile test.The results show that the M-C strength criterion obviously overestimates the tensile strength of the loess, whilst the joint strength overcomes the shortcoming of the M-C strength exaggerating the tensile strength of the loess.Therefore, the joint strength can reasonably evaluate the stress and displacement of the surrounding rock of the loess tunnel.

Equivalent Equation in the Form of Principal Stress of the Joint Strength
The ultimate equilibrium state of the joint strength on the σ-τ plane is shown in Figure 3, where σ 1 is major principal stress and σ 3 is minor principal stress.In order to unite the tensile strength and shear strength of loess, the strength take the tensile characteristics and shear characteristics into consideration is ca joint strength, and the corresponding failure strength curve is expressed in the hyperbola.Here "joint" means a kind of united consideration.
To fit a hyperbolic curve of the joint strength line of structural loess on the σ-Li et al. assumed that the horizontal axis intercept of the curve represents the strength σt, and that the M-C strength envelope was regarded as the asymptote [26 the joint strength formula corresponding to the dotted line, which can describe th and shear strengths of structural loess comprehensively and reasonably, was der the basis of the tensile and conventional triaxial shear strength tests.

(
) ( ) Figure 2 shows that the tensile strength determined by the M-C strength is m twice that of the tensile test.The results show that the M-C strength criterion ob overestimates the tensile strength of the loess, whilst the joint strength overco shortcoming of the M-C strength exaggerating the tensile strength of the loess.Th the joint strength can reasonably evaluate the stress and displacement of the surro rock of the loess tunnel.

Equivalent Equation in the Form of Principal Stress of the Joint Strength
The ultimate equilibrium state of the joint strength on the σ-τ plane is show ure 3, where σ1 is major principal stress and σ3 is minor principal stress.The stress Circle O 1 is tangent to the joint strength line, with the assumption that the coordinates of tangent Point F are (σ F , τ F ), because tangent Point F is on the joint strength line.Thus, Or In Formulas (3) and ( 4), The tangent slope k F at the differentiable Point F of Formula ( 4) is Then, the slope of the normal Line O 1 F at Point F is k N : According to the coordinates of Point F and Formula ( 7), the equation of straight Line The intersection of Formula (8) and the σ coordinate axis is the central coordinate of the Circle O The distance between Points O 1 and F is the radius of Circle O 1 .Thus, the equation of Circle O 1 is According to the half relation between the abscissa of centre of Circle O 1 and the value of (σ 1 + σ 3 ), we have Introducing Formula (4) into Formula (9) yields With the simultaneous solution of Formulas ( 10) and ( 11) and the elimination of σ F , the equivalent equation in the form of principal stress of the joint strength can be obtained as

Derivation of the New Modified Fenner Formula from Joint Strength
The influence parameters are shown in Figure 1.The surrounding rock of the loess adopts peak strength, which is expressed by the shear strength index.

Formula of Stress in Plastic Zone of Surrounding Rock of the Tunnel Based on Joint Strength
The surrounding rock of the plastic zone (a ≤ r ≤ R) is in the state of ultimate equilibrium.If the compressive stress is positive, then the major principal stress σ 1 is tangential stress σ θp and the minor principal stress σ 3 is the radial stress σ rp in the plastic zone.Introducing them into Formula (12) yields where σ rp is radial stress in the plastic zone and σ θp is tangential stress in the plastic zone.
To simplify the calculation, we let Introducing Formula ( 14) into Formula (13) yields With simultaneous solution of Formulas ( 14) and ( 15), the formulas of σ rp and σ θp based on joint strength can be obtained with η as a parameter.
The derivative can be obtained by Formula (16).
Introducing Formulas ( 16) and ( 17) into Formula (1) yields After separating the variables in Formula (18), Formula ( 19) can be obtained by the integral.ln where C is constant of integration.When r equals a, σ rp = σ a , η = η 0 .Introducing them into Formulas ( 16) and ( 19) yields By introducing Formula (21) into Formula (19), the equation between η and r can be obtained.
Formulas ( 16) and (22) show that the expressions of the tangential stress σ θp , the radial stress σ rp in the plastic zone and the vector diameter r of the calculated point regard η as a variable parameter based on the joint strength.By introducing the vector diameter r into Formula ( 22), the value of η the calculated point can be obtained.Introducing η into Formula ( 16) can determine the corresponding tangential stress σ θp and radial stress σ rp in the plastic zone.
3.1.2.Formula of Radius of Plastic Zone in Surrounding Rock of the Tunnel Based on Joint Strength According to Lomé's solution and boundary condition, the elastic stress [5] in the surrounding rock of the tunnel can be obtained from Formula (23).
where σ R is radial stress at the interface between the elastic and plastic zones and R is radius of plastic zone.The condition of stress continuity is satisfied at the interface between the elastic and plastic zones.When r equals R, Formula (24) holds.
where σ re is radial stress in the elastic zone and σ θe is tangential stress in the elastic zone.At the interface between the elastic and plastic zones, r equals R. Thus, Formula (23) becomes Formula (25).
At the interface between the elastic and plastic zones, r equals R and η equals η R .Accordingly, Formula (16) becomes Formula (26).
Introducing Formulas ( 24) and (25) into Formula (26) yields At the interface between the elastic and plastic zones, r equals R and η equals η R .Thus, plastic radius R based on joint strength can be derived from Formulas ( 22) and ( 28), as shown in Formula (29).

.2. Formula of Displacement in Surrounding Rock of the Tunnel Based on Joint Strength
At the interface between the elastic and plastic zones, the displacement of the surrounding rock of the tunnel can be determined from Formula (33) when r equals R and η equals η R .
Introducing Formulas ( 29) and (31) into Formula (33) can derive the displacement at the interface between the elastic and plastic zones based on joint strength.
With the assumption that the plastic zone conforms to the volume-compatibility condition in the deformation process, Formula (35) must hold.
After expanding Formula (35) and omitting higher-order terms of u 2 R and u 2 a , radial displacement of tunnel wall based on joint strength can be approximatively derived, as shown in Formula (36).
Introducing Formulas ( 29) and (34) into Formula (36) can obtain the radial displacement of the tunnel wall based on joint strength, as shown in Formula (37).

Results and Discussion
According to the abovementioned results, which include Formulas ( 16), ( 29), (32), (34) and (37), a series of analytical solutions comprises the modified Fenner formula based on joint strength.In this section, the difference between the modified Fenner formula based on joint strength and the modified Fenner formula based on M-C strength [7] for the unlined loess tunnel will be analysed and discussed.
For the unlined loess tunnel, the radius of the loess tunnel is 2.0 m, weight per unit volume γ is 15 kN/m 3 , cohesion c is 60 kPa, internal friction angle ϕ is 25 • , tensile strength σ t is 28 kPa, elastic modular E is 72 MPa and Poisson's ratio µ is 0.35.

Comparison of the Stresses in Surrounding Rock of Loess Tunnel
When the buried depth of unlined loess tunnel is 50 m, the uniform in situ stress σ s is 750 kPa.Thus, the distribution of stresses in the surrounding rock of the unlined loess tunnel of the modified Fenner formula based on M-C strength and the newly modified Fenner formula based on joint strength can be calculated and is shown in Figures 4 and 5.The tangential stresses of the modified Fenner formula based on M-C strength a the newly modified Fenner formula based on joint strength increase firstly and then d crease when the distance is far from the axis of the loess tunnel.In the meantime, t   The tangential stresses of the modified Fenner formula based on M-C stren the newly modified Fenner formula based on joint strength increase firstly and t crease when the distance is far from the axis of the loess tunnel.In the meanti The radial stresses of the modified Fenner formula based on M-C strength and the newly modified Fenner formula based on joint strength increase when the distance is far from the axis of the loess tunnel.However, the radial stress of surrounding rock of loess tunnel determined by modified Fenner formula based on joint strength is less than that determined by M-C strength.This is consistent with the conclusion in the literature [28] that the M-C strength criterion overestimates the tensile strength of loess, resulting in a deviation in the calculation results.
The tangential stresses of the modified Fenner formula based on M-C strength and the newly modified Fenner formula based on joint strength increase firstly and then decrease when the distance is far from the axis of the loess tunnel.In the meantime, the maximum tangential stresses occur at the interface of the elastic-plastic zone.
The tangential stress in the plastic zone of the surrounding rock determined by the modified Fenner formula based on joint strength is obviously smaller than that determined by the M-C strength, but the phenomenon is reversed in the elastic zone.The largest difference in tangential stress in the two kinds of modified Fenner formulas is observed at the elastic-plastic interface.

Comparison of the Radius of Plastic Zone and Radial Displacement
The radii of the plastic zone in the surrounding rock of the unlined loess tunnel of the modified Fenner formula based on M-C strength and the newly modified Fenner formula based on joint strength can be calculated is shown in Figure 6.The radii of the plastic zone calculated by the modified Fenner formula bas strength and the modified Fenner formula based on M-C strength increased w crease in tunnel depth.However, the radius of the plastic zone determined by fied Fenner formula based on joint strength was larger than that based on the la the same depth.
The radial displacement at the interface of elastic-plastic zone in the su rock of the unlined loess tunnel for the modified Fenner formula based on Mand the newly modified Fenner formula based on joint strength can be calcula shown in Figure 7.In the meantime, the radial displacement of the tunnel wall d by the two strength theories is shown in Figure 8.The radii of the plastic zone calculated by the modified Fenner formula based on joint strength and the modified Fenner formula based on M-C strength increased with the increase in tunnel depth.However, the radius of the plastic zone determined by the modified Fenner formula based on joint strength was larger than that based on the latter under the same depth.
The radial displacement at the interface of elastic-plastic zone in the surrounding rock of the unlined loess tunnel for the modified Fenner formula based on M-C strength and the newly modified Fenner formula based on joint strength can be calculated and is shown in Figure 7.In the meantime, the radial displacement of the tunnel wall determined by the two strength theories is shown in Figure 8.
The radial displacement at the interface of elastic-plastic zone in the surro rock of the unlined loess tunnel for the modified Fenner formula based on M-C s and the newly modified Fenner formula based on joint strength can be calculated shown in Figure 7.In the meantime, the radial displacement of the tunnel wall dete by the two strength theories is shown in Figure 8.The radial displacements determined by the two strength theories increase w increase in buried depth.However, the radial displacement determined by t strength theory is slightly larger than that determined by M-C strength.In a w M-C strength overestimates the tensile strength of structural loess, which lead smaller radius of the plastic zone and the smaller displacement of the surroundi of the unlined loess tunnel calculated by the modified Fenner formula based o strength.The joint strength can be used to reasonably evaluate the tensile and shear st of structural loess.Thus, the modified Fenner formula based on joint strength practical in the analysis of the stress environment and displacement field of sthe u ing rock of a loess tunnel.Therefore, a reasonable evaluation of the tunnel excava pact not only can eliminate potential collapse accidents in advance but also ca damage caused by ground subsidence to groundwater systems and ecological e ments.

Conclusions
A newly modified Fenner formula for a loess tunnel was derived on the bas joint-strength theory that can consider the tensile and shear properties of loess.T change in the stress environment in the surrounding rock of a loess tunnel unde tion of the average stratum stress was studied.Lastly, the differences in the stres bution, the radius of the plastic zone and the radial displacement were analysed a pared on the basis of the joint strength and M-C strength.The main conclusion follows: The radial displacements determined by the two strength theories increase with the increase in buried depth.However, the radial displacement determined by the joint strength theory is slightly larger than that determined by M-C strength.In a word, the M-C strength overestimates the tensile strength of structural loess, which leads to the smaller radius of the plastic zone and the smaller displacement of the surrounding rock of the unlined loess tunnel calculated by the modified Fenner formula based on M-C strength.
The joint strength can be used to reasonably evaluate the tensile and shear strengths of structural loess.Thus, the modified Fenner formula based on joint strength is more practical in the analysis of the stress environment and displacement field of sthe urrounding rock of a loess tunnel.Therefore, a reasonable evaluation of the tunnel excavation impact not only can eliminate potential collapse accidents in advance but also can avoid damage caused by ground subsidence to groundwater systems and ecological environments.

Conclusions
A newly modified Fenner formula for a loess tunnel was derived on the basis of the joint-strength theory that can consider the tensile and shear properties of loess.Then, the

Figure 1 .
Figure 1.Stress component of surrounding rock in plane-strain analysis.

Figure 1 .
Figure 1.Stress component of surrounding rock in plane-strain analysis.

Figure 2 .
Figure 2. M-C strength and the joint strength.

Figure 2
Figure 2 shows that the tensile strength determined by the M-C streng twice that of the tensile test.The results show that the M-C strength crite overestimates the tensile strength of the loess, whilst the joint strength shortcoming of the M-C strength exaggerating the tensile strength of the lo the joint strength can reasonably evaluate the stress and displacement of th rock of the loess tunnel.

Figure 2 .
Figure 2. M-C strength and the joint strength.

Figure 2 .
Figure 2. M-C strength and the joint strength.

Figure 3 .Figure 3 .
Figure 3. Ultimate equilibrium state in the joint strength.

Figure 4 .
Figure 4. Comparison of radial stresses in surrounding rock of loess tunnel.

Figure 5 .
Figure 5.Comparison of tangential stresses in surrounding rock of loess tunnel.

Figure 4 .
Figure 4. Comparison of radial stresses in surrounding rock of loess tunnel.

Figure 4 .
Figure 4. Comparison of radial stresses in surrounding rock of loess tunnel.

Figure 5 .
Figure 5.Comparison of tangential stresses in surrounding rock of loess tunnel.

Figure 5 .
Figure 5.Comparison of tangential stresses in surrounding rock of loess tunnel.

Figure 6 .
Figure 6.Relationship between radius of plastic zone and depth of loess tunnel.

Figure 6 .
Figure 6.Relationship between radius of plastic zone and depth of loess tunnel.

Figure 7 .
Figure 7. Relationship between radial displacement at elastic-plastic interface and depth o

Figure 7 .Figure 8 .
Figure 7. Relationship between radial displacement at elastic-plastic interface and depth of tunnel.

Figure 8 .
Figure 8. Relationship between radial displacement of the tunnel wall and depth of tunnel.
. Sci. 2023, 13, 6836 9 of 14 3.2.Formula of Elastic Stress and Displacement in Elastic Zone of Surrounding Rock of the Tunnel Based on Joint Strength 3.2.1.Formula of Elastic Stress in Surrounding Rock of the Tunnel Based on Joint Strength