Next Article in Journal
Study on the Aging Behavior of Asphalt Binder Exposed to Different Environmental Factors
Next Article in Special Issue
Inversion of Surrounding Red-Bed Soft Rock Mechanical Parameters Based on the PSO-XGBoost Algorithm for Tunnelling Operation
Previous Article in Journal
Algorithm Research on Detail and Contrast Enhancement of High Dynamic Infrared Images
Previous Article in Special Issue
Simulations of Deformation Failure Process and Refinement of Reinforcement Scheme in the Dabenliu Quarry Slope Using Discontinuous Deformation Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calculation Method of the Design Thickness of a Frozen Wall with Its Inner Edge Radially Incompletely Unloaded

1
State Key Laboratory for Geomechanics and Deep Underground Engineering, School of Mechanics and Civil Engineering, China University of Mining and Technology, Xuzhou 221116, China
2
Capital Construction Office, Jiangsu Normal University, Xuzhou 221116, China
3
YunLong Lake Laboratory of Deep Underground Science and Engineering, Xuzhou 221116, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(23), 12650; https://doi.org/10.3390/app132312650
Submission received: 3 November 2023 / Revised: 21 November 2023 / Accepted: 23 November 2023 / Published: 24 November 2023
(This article belongs to the Special Issue Advances in Failure Mechanism and Numerical Methods for Geomaterials)

Abstract

:
The technology for freezing shaft sinking is widely used for shafts to pass through deep, unstable alluvia with the continuous exploitation of mineral resources. Due to the technique using the sectional excavation and shaft lining construction adopted in deep alluvia, the radial stress at the inner edge of a frozen wall is incompletely unloaded. In this paper, a mechanical model was established for a frozen wall with its inner edge radially incompletely unloaded. A parameter, α, expressing the degree of being unloaded was introduced, and then a new method of designing and calculating the thickness of the frozen wall was proposed. The range of parameter α was estimated based on the frozen wall–shaft lining interaction forces from field data from a given project. The results indicate that the range of α can be chosen to be from 0.05 to 0.15 in deep alluvia. The design thickness of the frozen wall can be reduced by at least 5% for the frozen wall with the inner edge radially incompletely unloaded. The design thickness is significantly influenced by the strength and elastic modulus of the frozen soil and the elastic modulus of the surrounding unfrozen alluvium. The design and calculation method of frozen wall thickness can provide new ideas for guiding the design of frozen walls in deep alluvia.

1. Introduction

The technology for freezing shaft sinking is a special construction technology that utilizes artificial ground freezing methods to form a closed frozen geotechnical curtain (frozen wall) around a supporting structure (a shaft) by temporarily freezing the water in strata to improve the bearing capacity and water-sealing performance of the strata and protect the construction safety of shafts [1,2,3]. It has been gradually popularized and applied due to its outstanding advantages of strong adaptability, high reliability, and good technical and social benefits, and it has become the main construction technology for underground structures in deep water-bearing unstable soil layers and water-rich rock strata [4,5,6,7]. The depth that the shaft is sunk through an alluvium is a landmark indicator to measure the level of technology for freezing shaft sinking [8]. For example, the auxiliary shaft of Wanfu Coal Mine was sunk through alluvium to a depth of 754.96 m (world record) using freezing sinking technology [9]. More than 1100 vertical shafts have been constructed in China using this technology.
The thickness of the frozen wall is closely related to the safety and cost of the freezing shaft sinking project [10,11]. For example, the thickness of the frozen wall in the auxiliary shaft of Wanfu Coal Mine reaches 12.5 m, with an average temperature of −23 °C, which requires a cooling capacity of up to 15,407 kW. Ground pressure and ground temperature increase with the depth that the shaft is sunk through the alluvium, which leads to a significant increase in the technical difficulty of freezing shaft sinking. Therefore, how to ensure the safety of freezing shaft sinking while taking into account the economy can provide new ideas for guiding the design thickness of frozen walls in deep alluvia.
Theories for designing the frozen wall thickness have developed with the widespread use of freezing shaft sinking in alluvia. The earliest formula [12] can be traced back to the Lame formula based on the elastic single-layer, thick-walled cylinder model to guide the freezing shaft sinking process in topsoil at a depth of approximately 120 m or less. As the Lame formula was demonstrated to be too conservative, the Domke formula [13] was presented based on an elastoplastic single-layer, thick-walled cylinder model with the limitation that the plastic zone outer radius is smaller than the radius at which freezing pipes are set. This formula was widely used in practical freezing shaft sinking projects at depths of approximately 120 m to 300 m. Considering the frozen soil as a kind of friction material, the Klein formula [14] improved the Domke formula by replacing the Tresca yield criterion with the Mohr–Coulomb yield criterion. Removing Domke’s limitation on the outer radius of the frozen wall plastic zone, Sanger [15] proposed a fully plastic single-layer, thick-walled cylinder model formula. Recently, there has been a surge in the utilization of unloading models to enhance the analysis of frozen walls [16,17,18,19,20]. Because the Domke formula has large errors when the modulus ratio of the frozen wall to its surrounding unfrozen alluvium is less than 10, Yang [17] started with the geostress unloading mechanism and presented a formula that comprehensively considered the excavation unloading effect, the interaction between the frozen wall and the surrounding unfrozen alluvium, and Domke’s limitation on the outer radius of the frozen wall plastic zone. This formula was successfully applied in an alluvium depth range of 400 m to 700 m. The results showed that the radial stress at the outer boundary of the frozen wall following excavation was evidently lower than the initial geostress, leading to a decrease in the thickness of the frozen wall by 10% to 36% in comparison to the results derived from the Domke formula.
The geotechnical support characteristic curves suggest that the stability of the underground structure is highly influenced by the support effect. The movement of the inner boundary of the frozen wall is closely linked to the stability of the shaft wall. The existing unloading models do not take into account the effect of displacement and stress unloading at the inner edge of the frozen wall on the design thickness of the frozen wall.
To reduce the displacement of frozen walls during excavations, a technique using short-sectional excavation and shaft lining construction is adopted [21]; it causes interaction between the frozen wall and the shaft, and the radial stress in the inner edge of the frozen wall is incompletely unloaded. With the increasing demand for underground resources, the depth that the shaft is sunk through the alluvium increases, and the excavation height–radius ratio for each section is generally reduced to 1/3~2/3. The characteristics of the inner edge of the frozen wall being radially incompletely unloaded become more obvious, and the deformation of the frozen wall is subject to greater constraints. Therefore, the design process should consider the attributes of being radially incompletely unloaded at the inner boundary of the frozen wall.
In this paper, based on the construction technique using short-sectional excavation and shaft lining, a mechanical model was established for a frozen wall with its inner edge radially incompletely unloaded, comprehensively considering the in situ stress and the frozen wall–shaft lining interaction. Analytical solutions were derived. A parameter α expressing the degree of being unloaded was introduced, and then a new method of designing and calculating values for an elastoplastic frozen wall thickness was developed based on α and the solutions. The influences of the mechanical parameters of the shaft lining, frozen wall, and surrounding unfrozen alluvium on the design thickness of the frozen wall were also discussed. The range of parameter α was estimated based on the frozen wall–shaft lining interaction forces from field data from a given project. A concrete application case is presented at the end of this paper.

2. Mechanical Model and Solution

The process for the method of frozen wall thickness design in this paper is shown in Figure 1.
The schematic for the technique using short-sectional excavation and shaft lining is illustrated in Figure 2. Following the lining of the inner edge of the frozen wall, excavation proceeds downward, causing the working face to move downward, thus creating a new exposed section of the frozen wall. The frozen wall at the depth of the exposed section is selected as the research object. The frozen soil below the exposed section and the shaft wall above it jointly restrict the wall’s deformation. With the gradual sinking of the working face, the exposed section of the frozen wall continues to deform and begins to interact with the shaft. During the interaction between the research object and the shaft, the constraint on the object from the working face gradually decreases until it disappears completely, while the research object–shaft interaction force synchronously increases until it stabilizes. The constraint on the exposed section is approximated to be transformed into the freezing wall–shaft interaction. For the reasons above, on the basis of Yang’s model for the frozen wall where the inner edge is radially incompletely unloaded, with some basic assumptions added, a mechanic model for the frozen wall where the inner edge is radially completely unloaded is ultimately created in this paper.

2.1. Mechanical Model

The mechanical model is shown in Figure 3, where N, I, II, and III represent the shaft lining, the frozen wall plastic zone, the frozen wall elastic zone, and the surrounding unfrozen alluvium, respectively; r a , r b , and r c are the outer radius of the shaft lining, the excavation, and the outer radii of the frozen wall, respectively; and t is the design thickness of the frozen wall.
The following assumptions are made in the model.
(1) The shaft lining, the frozen wall, and the surrounding unfrozen alluvium are axisymmetric, infinitely long, and thick-walled cylinders, and the displacement of the surrounding unfrozen alluvium at r → ∞ is zero.
(2) The shaft lining and the surrounding unfrozen alluvium are homogeneous, isotropic, linear elastic materials; the frozen wall is an ideal elastoplastic material, and the stresses in the plastic zone satisfy the Mohr–Coulomb yield criterion.
(3) Before and after the formation of the frozen wall, the in situ stresses remain constant.
σ r 0 = σ θ 0 = p 0
where p 0 is the initial horizontal ground stress, σ represents the stress, and the subscripts r and θ represent the radial and circumferential directions, respectively. The positive and negative signs of the stresses follow the conventions of elasticity theory.
(4) At the instant that the sectional excavation is completed, the inner edge of the frozen wall is displaced in the radial direction and then blocked by the outer edge of the shaft, causing the interaction force between the frozen wall and the shaft to be p .

2.2. Initial, Boundary, and Continuity Conditions

The initial states of the frozen wall and the surrounding unfrozen alluvium prior to the excavation are
σ r = σ θ = p 0
u = 0
where u denotes the displacement, with a negative notation, following the sign convention of displacement in elasticity theory.
At the interface between the elastic and plastic zones of the frozen wall, the displacement and stress continuity conditions are [22]
u I r = r p = u I I r = r p
σ r I r = r p = σ r I I r = r p = p p
σ θ I r = r p = σ θ I I r = r p
At the interface between the frozen wall and the surrounding unfrozen alluvium, the displacement and stress continuity conditions are
u I I r = r c = u I I I r = r c
σ r I I r = r c = σ r I I I r = r c
At the infinity of the surrounding unfrozen alluvium, the stress and displacement boundary conditions are
u I I I r = 0
σ r I I I r = σ θ I I I r = p 0

2.3. Displacement and Stress of the Frozen Wall Elastic Zone

To facilitate the solution, the mechanical model shown in Figure 3 is split into two basic models as shown in Figure 4 and Figure 5.
The analytical solutions of stress and displacement for the elastic zone of the frozen wall are equal to the solutions for the unloading model shown in Figure 3 superimposed over the initial stress and displacement [23,24].
The unloading amount of the radial stress at the inner edge of the frozen wall elastic zone is
Δ p p = p p p 0
According to the elastic thick-walled cylinder theory [25,26] and combining Equations (7), (10) and (11), the stress and the displacement of the unloading model can be obtained. Based on the aforementioned superposition, analytical solutions of the stress and displacement for the frozen wall elastic zone and the surrounding unfrozen alluvium are obtained as follows:
σ r I I = Δ p p ζ ζ r p 2 r 2 r p 2 r c 2 + r p 2 r 2 r p 2 r c 2 1 + p 0
σ θ I I = Δ p p ζ + ζ r p 2 r 2 r p 2 r c 2 r p 2 r 2 r p 2 r c 2 1 + p 0
u I I = r Δ p p 1 + μ I I E I I 1 2 μ I I 1 ζ 1 r p 2 r c 2 1 2 μ I I + r p 2 r 2
σ r I I I = Δ p p ζ r c 2 r 2 + p 0
σ θ I I I = Δ p p ζ r c 2 r 2 + p 0
u I I I = Δ p p ζ 1 + μ I I I E I I I r c 2 r
where ζ = 2 2 + M r c 2 r p 2 1 ; M = E I I E I I I 1 1 + μ I I I 1 μ I I 2 + 2 + μ I I I μ I I 1 μ I I 2 ; E and μ are the elastic modulus and the Poisson’s ratio, respectively; and superscripts II and III represent the frozen wall elastic zone and the surrounding unfrozen alluvium, respectively.
From Equations (13) and (14), the solutions for the circumferential stress and the displacement at the inner edge of the frozen wall elastic zone are
σ θ I I r = r p = Δ p p 2 ζ r p 2 r c 2 1 r p 2 r c 2 1 + p 0
u I I r = r p = r p Δ p p 1 + μ I I E I I 1 2 μ I I 2 1 ζ 1 μ I I 1 r p 2 r c 2

2.4. Displacement and Stress of the Frozen Wall Plastic Zone

From the schematic representation of ground–support interaction [27,28] and the intersection of the ground reaction and the support characteristic curves, the radial stress and the displacement at the outer edge of the shaft lining and the radial stress and the displacement at the inner edge of the frozen wall are determined to be
σ r N   r = r a = p
σ r I   r = r b = p
σ r N   r = r b = p k
u I r = r b = r a r b + u N r = r a
where superscript N represents the shaft lining, k is the stiffness of the shaft’s outer edge, and p represents the interaction force between the frozen wall and the shaft lining.
The Mohr–Coulomb yield criterion describes the relationship between the first and third principal stresses, and the axial stress is generally regarded as the intermediate principal stress [29,30,31]. The relationship between the radial and circumferential stresses in the frozen wall plastic zone is [32,33,34]
σ θ I = A σ r I σ c
where A = 1 + sin φ 1 sin φ , σ c = 2 c tan 4 5 ° φ 2 . φ , c , and σ c are the internal friction angle, cohesion, and uniaxial compressive strength of the frozen soils, respectively. Superscript I represents the frozen wall plastic zone.
The static equilibrium differential equation [35] for the frozen wall plastic zone (see Figure 5) is
d σ r d r + σ r σ θ r = 0
Substituting Equations (24) and (25) and combining them with Equation (21), solutions of the radial and circumferential stresses for the plastic zone of the frozen wall are
σ r I   = σ c A 1 r r b A 1 1 + p r r b A 1
σ θ I   = A σ c A 1 r r b A 1 1 + A p r r b A 1 σ c
From Equations (5) and (26),
p p = σ r I r = r p = σ c A 1 r p r b A 1 1 + p r p r b A 1
From Equation (27), the solution of the circumferential stress at the outer edge of the frozen wall plastic zone is
σ θ I r = r p = A σ c A 1 r p r b A 1 1 + A p r p r b A 1 σ c
Substituting Equations (18) and (29) into Equation (9), the relationship among rb, rc, and rp is obtained:
r p r b A 1 = p 0 p A 1 σ c 1 A 1 σ c p A 1 σ c 1 M 1 2 r p 2 r c 2 + A + 1 2 A 1
If r b = r a and p = 0 , the mechanical model in this paper (Figure 3) degenerates to that presented by Yang [17] for the frozen wall where the inner edge is completely unloaded in a radial manner; thus, Equation (30) degenerates to
r p r a A 1 = 1 A 1 p 0 σ c 1 M 1 2 r p 2 r c 2 + A + 1 2 A 1
Correspondingly, Equation (31) is the same as the formula presented by Yang [17].
Geomaterials do not follow the associated flow rule, and their plastic potential function φ is generally expressed as [36]
ϕ = β σ r σ θ
where β = 1 + sin ψ 1 sin ψ , and ψ is the dilation angle of frozen soil, 0 ≤ ψ φ .
Ignoring the contribution of the elastic deformation and according to plastic potential theory, the plastic strains of the frozen wall satisfy the following relationship:
ε r I + β ε θ I = 0
where ε r I and ε θ I are the radial and circumferential strains, respectively.
The geometric equations for the axisymmetric plane strain are
ε r = d u d r
ε θ = u r
The substitution of Equations (34) and (35) into Equation (33) results in
d u d r + β u r = 0
Combining Equation (36) with Equations (4) and (19), the analytical solution of the displacement for the frozen wall plastic zone is
u I = r p E I I 1 + μ I I r p r β σ c A 1 r p r b A 1 1 + p r p r b A 1 p 0 2 M 1 1 2 μ I I r c 2 r p 2 2 M 1 + r c 2 r p 2

3. Design for the Frozen Wall Thickness

3.1. Calculation Method and Formula

For the frozen wall with its inner edge radially incompletely unloaded, the deformation is constrained and the displacement of the inner edge is limited to a part of that with the inner edge radially completely unloaded. To express the degree of being unloaded, parameter α is introduced using the ratio of the displacement of the inner edge of the frozen wall with its inner edge radially incompletely unloaded to that of the inner edge of the same frozen wall with its inner edge radially completely unloaded; see Equation (38):
α = 1 u I r = r b U I r = r b
where α decreases with the degree of being unloaded, and α = 0 represents that the inner edge of the frozen wall is radially completely unloaded. U I r = r b is the displacement of the inner edge of the same frozen wall with its inner edge radially completely unloaded.
Substituting rp = Rp and p = 0 into Equation (37), U I r = r b is expressed as
U I r = r b = R p E I I 1 + μ I I R p r b β σ c A 1 R p r b A 1 1 p 0 2 M 1 1 2 μ I I r c 2 R p 2 2 M 1 + r c 2 R p 2
where Rp is the outer radius of the plastic zone of the same frozen wall with its inner edge radially completely unloaded.
Substituting rp = Rp and p = 0 into Equation (30) Rp can be solved using
R p r b A 1 = p 0 σ c + 1 A 1 1 M 1 2 R p 2 r c 2 + A + 1 2 A 1
Following Domke’s limitation imposed on the outer radius of the frozen wall plastic zone [13], r p = r b r c is substituted into Equations (30) and (37). Combining Equation (38) with Equations (38)–(40), the formulas for the design thickness of the frozen wall with its inner edge radially incompletely unloaded are determined as follows:
r c r b A 1 2 = p 0 p A 1 σ c 1 A 1 σ c p A 1 σ c 1 M 1 2 r b r c + A + 1 2 A 1 r c 1 + β 2 r b 1 β 2 1 + μ I I E I I σ c A 1 + p r c r b A 1 2 + σ c A 1 p 0 2 M 1 1 2 μ I I r c r b 2 M 1 + r c r b = p k + r a r b = R p 1 α 1 + μ I I E I I R p r b β σ c A 1 R p r b A 1 1 p 0 2 M 1 1 2 μ I I r c 2 R p 2 2 M 1 + r c 2 R p 2 R p r b A 1 = 1 A 1 p 0 σ c 1 M 1 2 R p 2 r c 2 + A + 1 2 A 1 t = r c r b
ra has a known value, and parameter α in Equation (41) will be discussed later. The mechanical parameters of the shaft lining, the frozen wall, and the surrounding unfrozen alluvium are substituted into Equation (41) to obtain t ,  p , r b , r c , and R p .
If the inner edge is radially completely unloaded, α = 0, p = 0, R p = r p , and u I r = r b = U I r = r b ; thus, Equation (41) degenerates to
r c r b A 1 2 = p 0 p A 1 σ c 1 A 1 σ c p A 1 σ c 1 M 1 2 r b r c + A + 1 2 A 1 t = r c r b
Correspondingly, Equation (42) is the same as the formula presented by Yang [17] for the frozen wall where the inner edge is completely unloaded in a radial manner.

3.2. Influences of the Parameters on the Design Thickness of the Frozen Wall

To facilitate the analysis of the influences of mechanical parameters, σc, EII, EIII, k, and t are made dimensionless, as shown in Equation (43):
t ~ = t r a ,   σ ~ c = σ c p 0 , E ~ I I = E I I p 0 , E ~ I I I = E I I I p 0 , k ~ = k r a p 0
Given typical engineering parameters, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 show that t ~ decreases significantly with α.
Figure 6 shows that t ~ decreases significantly with σ ~ c . The larger α is, the slower t ~ decreases with σ ~ c . α has a distinct effect on t ~ while σ ~ c 1 , and it has a slight effect while σ ~ c 1.2 . The difference between the curves decreases with σ ~ c . Figure 6 indicates that the smaller the radial unloading of the frozen wall, the thinner the design thickness, especially when the strength of the wall is relatively low.
Figure 7 and Figure 8 show that t ~ decreases significantly with E ~ I I I and increases significantly with E ~ I I . The higher the elastic modulus of the frozen wall, or the lower the elastic modulus of the surrounding unfrozen alluvium, the better for the stability of the frozen wall.
Figure 9 shows that t ~ decreases merely by 0.15% when k ~ varies from the typical value to infinity. Figure 10, Figure 11 and Figure 12 show that t ~ increases slightly with ψ, μII, and μIII, and the maximum increases are 3.5%, −2%, and 7.2%, while α = 0.3, respectively. With the change in these four variables, the curves of the design thickness of the frozen wall are maintained at their respective relatively stable levels.
Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 reveal that the design thickness of the frozen wall varies significantly in sensitivity to various influencing factors. Further, by means of single-factor analysis of variance in mathematics, quantitative data on the sensitivity of freezing wall design thickness to various factors by different values of α are shown in Table 1.
Of the many factors, identifying the dominant ones is crucial, as they have a significant impact on design, construction, and similar simulation tests related to the frozen wall; therefore, they require special attention.
As seen in Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 and Table 1, the sensitivity of the design thickness of frozen walls to the strength and elastic modulus of the frozen soil, and the elastic modulus of the surrounding unfrozen alluvium, change significantly as the value of α increases. When α is kept constant, if the frozen wall design thickness varies significantly with a factor, the variance associated with this factor in Table 1 is also relatively large, indicating that the frozen wall design thickness has a high sensitivity to this factor. Table 1 shows that the mean variance of frozen wall strength at different values of α clearly exceeds that of other factors. The ranking of these mean variances indicates that the strength and elastic modulus of the frozen soil and the elastic modulus of the surrounding unfrozen alluvium are the main influencing factors. The dilation angle and Poisson’s ratio of frozen soil, Poisson’s ratio of the surrounding unfrozen alluvium, and the stiffness of the shaft’s outer edge are secondary influencing factors.

4. Case Example

The background of this case example is a shaft to be constructed in North China with an alluvium depth exceeding 800 m. Equation (41) is used to calculate the design thickness t of the frozen wall at depths of 800 m, 900 m, 1000 m, and 1100 m.
In this case example [37], the average temperature of the frozen wall is −35 °C. The elastic modulus, Poisson’s ratio, internal friction angle, dilation angle, and uniaxial compressive strength of the frozen wall are EII = 400 MPa, μII = 0.2, φ = 6°, ψ = 0°, and σc = 12 MPa, respectively. The elastic modulus and Poisson’s ratio of the surrounding unfrozen alluvium are EIII = 150 MPa and μIII = 0.2, respectively. The outer radius of the shaft lining ra is 5 m, and the stiffness of the outer edge of the shaft lining k is 2000 MPa/m.
The initial horizontal ground pressure is determined by [37,38]
p 0 = 0.013   h
where p0 is the initial horizontal ground stress (MPa) and h is the depth (m).
Figure 13 shows the curves of t and p with h and different values of α . Based on Figure 13 and Table 2, the analyses are as follows:
(1) When α is constant, t increases nonlinearly and rapidly with h , and p decreases with h . When h is constant, t decreases with α, and p increases with α.
(2) With the parameter α kept constant, the differences between the frozen thickness–depth curves for both radially incomplete unloading and radial complete unloading conditions enlarge with the increase in depth. When the t curve with α = 0 is used as a baseline, t is approximately 95% ( α = 0.05), 90% ( α = 0.1), 85% ( α = 0.15), and 80% (α = 0.2) of the baseline at depths ranging from 800 m to 1100 m. When α increases by 0.05, t decreases by approximately 0.33 m ( h = 800 m), 0.42 m ( h = 900 m), 0.54 m ( h = 1000 m), and 0.68 m ( h = 1100 m), indicating that the effect of α on t increases with h .
(3) Young’s formula is derived based on the model for the frozen wall where the inner edge is completely unloaded in a radial manner, which should be a special case for α = 0 in this paper. Compared with results from Yang’s formula, the results when α = 0 in this paper should be the same as those when using Yang’s formula. As can be seen from Table 2, there is no difference between the results of the two methods, which verifies the effectiveness and accuracy of the method proposed in this paper.
The above analyses show that α has an important effect on the design thickness of frozen walls, and the effect intensifies as the depth increases. α is a construction parameter that depends on the excavation height–radius ratio and the extent of the freezing of the soils. The smaller the excavation height–radius ratio, the larger the value of α should be. Based on field data, α can be obtained using the results of the interaction force p between the frozen wall and the shaft lining. Based on to the field data of the freezing shaft sinking of Wanfu Mine [39], the results of p from the auxiliary shaft at the alluvium depth of 697.5 m are shown in Figure 11, which is of great reference. Prior to the construction of a section shaft lining, earth pressure cell sensors were affixed as evenly as possible to the inner edge of the frozen wall. After the construction of the section shaft lining, the interaction force p between the frozen wall and the shaft lining increases at a decreasing rate with the progress of construction, and finally tends towards steadiness, as shown in Figure 14. Due to the complexity of the construction, the duration curves of p in Figure 14 do not overlap ideally with each other as expected, which is a common problem in field surveys [40]. When α = 0.05, the difference between p in Figure 13 at the depth of 800 m and the max steady value (sensor NO. 6) in Figure 14 is no more than 0.4 Mpa, and when α = 0.15, the difference between p in Figure 11 at the depth of 800 m and the min steady value (sensor NO. 2) in Figure 14 is less than 0.3 Mpa, and thus the upper limit and lower limit of α are approximately 0.15 and 0.05, respectively.
According to the calculation results shown in Table 2, when α is chosen as the lower limit (α = 0.05), the calculated design thickness of the frozen wall is approximately 95% of that of the frozen wall with the inner edge radially completely unloaded. When α is chosen as the upper limit (α = 0.15), the calculated design thickness of the frozen wall is approximately 85% of that of frozen wall with the inner edge radially completely unloaded. These results indicate that the design thickness of the frozen wall in deep alluvia can be reduced by at least 5% for the frozen wall with the inner edge radially incompletely unloaded.

5. Conclusions

According to the construction technique using sectional excavation and shaft lining, a mechanical model was established for a frozen wall with its inner edge radially incompletely unloaded, comprehensively considering the in situ stress and the frozen wall–shaft lining interaction. Analytical solutions were derived. A new method of designing and calculating an elastoplastic frozen wall thickness was developed on the solutions. The main conclusions are as follows:
(1) The increase in the design thickness with the depth, based on the new model, is smaller than that based on the mechanical model of radial complete unloading at the inner edge of the frozen wall.
(2) The design thickness of the frozen wall is significantly influenced by the strength and elastic modulus of the frozen soil and the elastic modulus of the surrounding unfrozen alluvium. The dilation angle and Poisson’s ratio of frozen soil, Poisson’s ratio of the surrounding unfrozen alluvium, and the stiffness of the shaft’s outer edge are secondary influencing factors.
(3) The degree of being radially unloaded has an important effect on the design thickness of the frozen wall, and the effect becomes more pronounced with increasing depth. Using the new method of designing and calculating the frozen wall thickness, the design thickness in the case example of this paper can be reduced by at least 5%. This paper provides theoretical support for guiding the design of frozen wall thickness in deep alluvia.

Author Contributions

Conceptualization, T.H., W.Y. and C.H.; Methodology, T.H. and Z.Y.; Formal analysis, W.Y. and C.H.; Data curation, T.H. and C.H.; Writing—original draft, C.H.; review and editing, T.H., Z.Y. and W.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, grant number 2016YFC0600904, the National Natural Science Foundation of China, grant number 41171072, and the Foundation Research Project of Xuzhou, grant number KC22061.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to a degree of secrecy surrounding the National Key Research and Development Program of China.

Acknowledgments

We thank the anonymous reviewers for their professional suggestions in improving this submission. We thank Haipeng Li and Xin Huang for the valuable advice and support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cai, H.; Yang, L.; Pang, C.; Li, M.; Lu, C.; Hong, R. Model Test Study on Natural Thawing Temperature Field of Artificial Ground Frozen Wall. Sustainability 2023, 15, 3186. [Google Scholar] [CrossRef]
  2. Harris, J.S. Ground Freezing in Practice; Thomas Telford: London, UK, 1995. [Google Scholar]
  3. Hu, X.; Wu, Y.; Li, X. A Field Study on the Freezing Characteristics of Freeze-Sealing Pipe Roof Used in Ultra-Shallow Buried Tunnel. Appl. Sci. 2019, 9, 1532. [Google Scholar] [CrossRef]
  4. Andersland, O.B.; Ladanyi, B. Frozen Ground Engineering; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2004. [Google Scholar]
  5. Dorn, E.; Kaledin, O. Modern and innovative shaft sinking and construction technology with examples from current projects/Moderne und innovative Schachtbautechnik am Beispiel aktueller Abteufprojekte. Geomech. Tunn. 2013, 6, 574–581. [Google Scholar] [CrossRef]
  6. Schmall, P.C.; Braun, B. Ground Freezing—A Viable and Versatile Construction Technique. In Proceedings of the 13th International Conference on Cold Regions Engineering, Orono, ME, USA, 23–26 July 2006; 2006; pp. 1–11. [Google Scholar]
  7. Zhou, Y.; Zhou, G. Intermittent freezing mode to reduce frost heave in freezing soils—Experiments and mechanism analysis. Can. Geotech. J. 2012, 49, 686–693. [Google Scholar] [CrossRef]
  8. Xie, H.; Zhang, R.; Zhang, Z.; Gao, M.; Li, C. Reflections and explorations on deep earth science and deep earth engineering technology. J. China Coal Soc. 2023, 1–21. [Google Scholar]
  9. Wang, Y.S.; Yang, R. Monitoring and analysis of the stress and deformation of shaft lining and the influence of freezing tube fracture in deep topsoil. Cold Reg. Sci. Technol. 2022, 193, 103420. [Google Scholar] [CrossRef]
  10. Zhou, F.C.; Zhou, P.; Li, J.Y.; Ge, T.C.; Lin, J.Y.; Wang, Z.J. Key Parameters Design Method of AGF Method for Metro Connecting Passage in Water-Rich Coastal Area. KSCE J. Civ. Eng. 2022, 26, 5301–5317. [Google Scholar] [CrossRef]
  11. Alzoubi, M.A.; Zueter, A.; Nie-Rouquette, A.; Sasmito, A.P. Freezing on demand: A new concept for mine safety and energy savings in wet underground mines. Int. J. Min. Sci. Technol. 2019, 29, 621–627. [Google Scholar] [CrossRef]
  12. Andersland, O.B.; Ladanyi, B. An Introduction to Frozen Ground Engineering; Springer Science & Business Media: Dordrecht, The Netherlands, 2013. [Google Scholar]
  13. Domke, O. Über die Beanspruchung der Frostmauer beim Schachtabteufen nach dem Gefrierverfahren. Glückauf 1915, 51, 1129–1135. [Google Scholar]
  14. Klein, J. Present State of Freeze Shaft Design in Mining. In Developments in Geotechnical Engineering; Farmer, I.W., Ed.; Elsevier: Amsterdam, The Netherlands, 1981; Volume 32, pp. 147–153. [Google Scholar]
  15. Sanger, F.J.; Sayles, F.H. Thermal and rheological computations for artificially frozen ground construction. Eng. Geol. 1979, 13, 311–337. [Google Scholar] [CrossRef]
  16. Yang, W.; Du, Z.; Yang, Z.; Bo, D. Plastic design theory of frozen soil wall based on interaction between frozen soil wall and surrounding rock. Chin. J. Geotech. Eng. 2013, 35, 1857–1862. [Google Scholar]
  17. Yang, W.; Yang, Z.; Bo, D. Elastic-plastic design theory of frozen soil wall based on interaction between frozen wall and surrounding rock. Chin. J. Geotech. Eng. 2013, 35, 175–180. [Google Scholar]
  18. Zhang, W.; Wang, B.; Wang, Y. Elastic Analysis of Nonhomogeneous Frozen Wall under Nonaxisymmetric Ground Stress Field and in State of Unloading. Adv. Mater. Sci. Eng. 2018, 2018, 2391431. [Google Scholar] [CrossRef]
  19. Hu, X.D. Determination of Load on Frozen Soil Wall in Unloaded State. J. Tongji Univ. (Nat. Sci.) 2002, 1, 6–10. [Google Scholar]
  20. Bo, Z.; Weihao, Y.; Cheng, J. An elastoplastic design method for deep frozen wall considering shaft-flankdisplacement convergence. J. Min. Saf. Eng. 2023, 40, 112–119. [Google Scholar]
  21. Zhang, Y.; Sego, D.C.; Morgenstern, N.R. Deformation of artificially frozen shafts during excavation. In Proceedings of the Seventh International Symposium on Ground Freezing (Ground Freezing 94), Nancy, France, 24–28 October 1994. [Google Scholar]
  22. Yu, H.S. Cavity Expansion Methods in Geomechanics; Springer: Dordrecht, The Netherlands, 2000. [Google Scholar]
  23. Yang, R.; Wang, Q.; Yang, L. Closed-form elastic solution for irregular frozen wall of inclined shaft considering the interaction with ground. Int. J. Rock Mech. Min. Sci. 2017, 100, 62–72. [Google Scholar] [CrossRef]
  24. Lu, A.; Zhang, N.; Xu, Y.; Cui, P. Stress-displacement solution for a lined vertical borehole due to non-axisymmetric in situ stresses. Int. J. Rock Mech. Min. Sci. 2013, 57, 64–74. [Google Scholar] [CrossRef]
  25. Brady, B.H.; Brown, E.T. Rock Mechanics for Underground Mining; Springer: Dordrecht, The Netherlands, 2004. [Google Scholar]
  26. Vincenzo, V. Circular Cylinders and Pressure Vessels; Springer International Publishing: Cham, Switzerland, 2014. [Google Scholar]
  27. Carranza-Torres, C.; Fairhurst, C. The elasto-plastic response of underground excavations in rock masses that satisfy the Hoek-Brown failure criterion. Int. J. Rock Mech. Min. Sci. 1999, 36, 777–809. [Google Scholar] [CrossRef]
  28. Carranza-Torres, C.; Fairhurst, C. Application of the Convergence-Confinement method of tunnel design to rock masses that satisfy the Hoek-Brown failure criterion. Tunn. Undergr. Space Technol. 2000, 15, 187–213. [Google Scholar] [CrossRef]
  29. Liu, C. Distribution Laws of in-Situ Stress in Deep Underground Coal Mines. Procedia Eng. 2011, 26, 909–917. [Google Scholar] [CrossRef]
  30. Reed, M.B. The influence of out-of-plane stress on a plane strain problem in rock mechanics. Int. J. Numer. Anal. Methods Geomech. 1988, 12, 173–181. [Google Scholar] [CrossRef]
  31. Cai, M.; Chen, C.; Peng, H.; Ji, H.; Qiao, L.; Tan, Z. In-situ stress measurement by hydraulic fracturing technique in deep position of Wanfu coal mine. Yan Shi Li Xue Yu Gong Cheng Xue Bao 2006, 25, 1069–1074. [Google Scholar]
  32. Vrakas, A. A finite strain solution for the elastoplastic ground response curve in tunnelling: Rocks with non-linear failure envelopes. Int. J. Numer. Anal. Methods Geomech. 2017, 41, 1077–1090. [Google Scholar] [CrossRef]
  33. Yu, H.S.; Rowe, R.K. Plasticity solutions for soil behaviour around contracting cavities and tunnels. Int. J. Numer. Anal. Methods Geomech. 1999, 23, 1245–1279. [Google Scholar] [CrossRef]
  34. Peng, S.; Xu, Y.; Cao, G.; Pei, L. Research on the Elastoplastic Theory and Evolution Law of Plastic Zone Contours of Horizontal Frozen Walls under Nonuniform Loads. Appl. Sci. 2023, 13, 9398. [Google Scholar] [CrossRef]
  35. Vrakas, A.; Anagnostou, G. A simple equation for obtaining finite strain solutions from small strain analyses of tunnels with very large convergences. Géotechnique 2015, 65, 936–944. [Google Scholar] [CrossRef]
  36. Vrakas, A.; Anagnostou, G. A finite strain closed-form solution for the elastoplastic ground response curve in tunnelling. Int. J. Numer. Anal. Methods Geomech. 2014, 38, 1131–1148. [Google Scholar] [CrossRef]
  37. Zhang, B.; Yang, W.; Wang, B. Plastic Design Theory of Frozen Wall Thickness in an Ultradeep Soil Layer Considering Large Deformation Characteristics. Math. Probl. Eng. 2018, 2018, 8513413. [Google Scholar] [CrossRef]
  38. Weng, J.J. Special Construction Engineering of Mine Shaft and Drift; Coal Industry Press: Beijing, China, 1991. [Google Scholar]
  39. Yang, Z. Study on the Elastoplastic Design Theory of Icewall in Thick Soil Based on the Interaction with Surrounding Stratum. Ph.D. Thesis, China University of Mining and Technology, Xuzhou, China, 2018. [Google Scholar]
  40. Wang, T.; Zhang, H.; Zhou, G.; Liu, X. Evaluation of variability characteristics of mechanical parameters of artificially frozen clay in deep alluvium. Cold Reg. Sci. Technol. 2020, 171, 102978. [Google Scholar] [CrossRef]
Figure 1. Flowchart for the method of frozen wall thickness design.
Figure 1. Flowchart for the method of frozen wall thickness design.
Applsci 13 12650 g001
Figure 2. Schematic for the technique of short-sectional excavation and shaft lining.
Figure 2. Schematic for the technique of short-sectional excavation and shaft lining.
Applsci 13 12650 g002
Figure 3. Mechanical model for a frozen wall where the inner edge is partially unloaded in a radial manner.
Figure 3. Mechanical model for a frozen wall where the inner edge is partially unloaded in a radial manner.
Applsci 13 12650 g003
Figure 4. Unloading model.
Figure 4. Unloading model.
Applsci 13 12650 g004
Figure 5. Mechanical model for the plastic zone of a frozen wall where the inner edge is partially unloaded in a radial manner.
Figure 5. Mechanical model for the plastic zone of a frozen wall where the inner edge is partially unloaded in a radial manner.
Applsci 13 12650 g005
Figure 6. Influence of σ ~ c on t ~ .
Figure 6. Influence of σ ~ c on t ~ .
Applsci 13 12650 g006
Figure 7. Influence of E ~ I I I on t ~ .
Figure 7. Influence of E ~ I I I on t ~ .
Applsci 13 12650 g007
Figure 8. Influence of E ~ I I on t ~ .
Figure 8. Influence of E ~ I I on t ~ .
Applsci 13 12650 g008
Figure 9. Influence of k ~ on t ~ .
Figure 9. Influence of k ~ on t ~ .
Applsci 13 12650 g009
Figure 10. Influence of ψ on t ~ .
Figure 10. Influence of ψ on t ~ .
Applsci 13 12650 g010
Figure 11. Influence of μII on t ~ .
Figure 11. Influence of μII on t ~ .
Applsci 13 12650 g011
Figure 12. Influence of μIII on t ~ .
Figure 12. Influence of μIII on t ~ .
Applsci 13 12650 g012
Figure 13. Calculation results of case example.
Figure 13. Calculation results of case example.
Applsci 13 12650 g013
Figure 14. Duration curves of p from field data of auxiliary shaft sinking of Wanfu Mine.
Figure 14. Duration curves of p from field data of auxiliary shaft sinking of Wanfu Mine.
Applsci 13 12650 g014
Table 1. Results of single-factor analysis of variance.
Table 1. Results of single-factor analysis of variance.
FactorαAverage
00.10.20.3
σ ~ c 3.99413.21042.52811.95562.9221
E ~ I I I 0.03150.03930.0470.05840.044
E ~ I I 0.00170.00300.00450.00640.0039
μIII0.00140.00170.00200.00240.001
ψ0.0000.00020.0000.00060.0003
μII0.00014.0 × 10−63.9 × 10−50.00020.0001
k ~ 9.3 × 10−102.6 × 10−61.3 × 10−62.2 × 10−61.5 × 10−6
Table 2. Computational results of frozen wall thickness t (m) with different values of α.
Table 2. Computational results of frozen wall thickness t (m) with different values of α.
Formulah/m
80090010001100
Yang’s formula [17]6.858.9011.3214.21
α 06.858.9011.3214.21
0.056.508.4610.7513.51
0.16.158.0010.2012.79
0.155.837.599.6512.11
0.25.487.159.1011.43
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, C.; Yang, Z.; Han, T.; Yang, W. Calculation Method of the Design Thickness of a Frozen Wall with Its Inner Edge Radially Incompletely Unloaded. Appl. Sci. 2023, 13, 12650. https://doi.org/10.3390/app132312650

AMA Style

Hu C, Yang Z, Han T, Yang W. Calculation Method of the Design Thickness of a Frozen Wall with Its Inner Edge Radially Incompletely Unloaded. Applied Sciences. 2023; 13(23):12650. https://doi.org/10.3390/app132312650

Chicago/Turabian Style

Hu, Chenchen, Zhijiang Yang, Tao Han, and Weihao Yang. 2023. "Calculation Method of the Design Thickness of a Frozen Wall with Its Inner Edge Radially Incompletely Unloaded" Applied Sciences 13, no. 23: 12650. https://doi.org/10.3390/app132312650

APA Style

Hu, C., Yang, Z., Han, T., & Yang, W. (2023). Calculation Method of the Design Thickness of a Frozen Wall with Its Inner Edge Radially Incompletely Unloaded. Applied Sciences, 13(23), 12650. https://doi.org/10.3390/app132312650

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

Article Metrics

Back to TopTop