Next Article in Journal
Regional Governance for Food System Transformations: Learning from the Pacific Island Region
Next Article in Special Issue
Assessment of the Blasting Efficiency of a Long and Large-Diameter Uncharged Hole Boring Method in Tunnel Blasting Using 3D Numerical Analysis
Previous Article in Journal
Vegetation Monitoring and Post-Fire Recovery: A Case Study in the Centre Inland of Portugal
Previous Article in Special Issue
A New Discrete Form of Hoek–Brown Criterion and Its Application to Limit Equilibrium Analysis of Rock Slope Stability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Active Earth Pressure for Narrow Unsaturated Backfills Considering Soil Arching Effect and Interlayer Shear Stress

1
School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
2
Henan Provincial Transportation Planning and Design Institute Co., Ltd., Zhengzhou 450015, China
3
School of Architecture and Civil Engineering, Xiamen University, Xiamen 361005, China
4
College of Civil Engineering, University of South China, Hengyang 421001, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(19), 12699; https://doi.org/10.3390/su141912699
Submission received: 13 September 2022 / Revised: 25 September 2022 / Accepted: 27 September 2022 / Published: 6 October 2022
(This article belongs to the Special Issue Advances in Rock Mechanics and Geotechnical Engineering)

Abstract

:
Affected by climatic conditions such as rainfall, evaporation and air temperature, most of the backfill soil is in an unsaturated state, and the decrease in matric suction leads to the failure of the retaining structure. In view of this, this study takes the vertical rigid retaining wall with narrow unsaturated backfill as the research object, assuming that the backfill behind the wall forms a circular soil arch and considering the interlayer shear stress; thus, the analytical solution of the active earth pressure of narrow unsaturated soil is derived based on the thin layer element method. The reliability of this method is verified by comparing with the experimental and existing theoretical results. A parameter analysis demonstrates that with the increase in the interface friction angle of the moving wall–soil, the average shear stress coefficient of zone I and zone II increases gradually, but with the increase in the interface friction angle of the fixed wall–soil, the average shear stress coefficient of zone I decreases; with the increase in effective internal friction angle and effective cohesion, the active earth pressure decreases and the tension crack depth increases; with the increase in the interface friction angle, the active earth pressure in the upper part of the retaining wall increases slightly, while the active earth pressure in the lower part decreases obviously; with the increase in matric suction, the active earth pressure first decreases rapidly and then increases gradually, and the tension crack depth first increases and then decreases, but the distribution pattern of the horizontal active earth pressure remains unchanged; the active earth pressure decreases with the decrease in the aspect ratio, and when the aspect ratio is smaller, the attenuation is more obvious; until the aspect ratio reaches a certain value, the active earth pressure is basically unchanged.

1. Introduction

The lateral pressure of the retaining wall is a key factor affecting the safety and stability of the retaining wall. An accurate estimation of the earth pressure behind a retaining wall is the premise of rational design. Coulomb’s and Rankine’s earth pressure theory have been used by engineering designers until now because of their simple form and clear concept. Both of them assume that the backfill behind the wall is semi-infinite and is fully saturated or dry. Furthermore, the earth pressure obtained by the two theories is linearly distributed along the depth. However, it is proven that the active earth pressure distribution behind the wall is nonlinear due to the rough wall–soil interfaces [1,2,3,4,5].
With the rapid development of urban construction, underground buildings are becoming more and more intensive, so that the backfills behind many retaining walls cannot meet the conditions of semi-infinite soil. The adjacent underground structures prevent the sliding surface in the backfill from fully developing to the backfill surface, resulting in a large difference between the earth pressure behind the wall and that in the semi-infinite soil case. Several studies have demonstrated that the active earth pressure acting on the retaining wall decreases with the narrowing of the backfill behind the wall [3,5,6,7,8,9,10,11,12]. In addition, some scholars have discussed the shape of the slip surface in the soil under the condition of narrow backfill. A simple but effective method is to assume that the linear sliding surface is truncated by the adjacent wall to form a trapezoidal sliding wedge [5,13]. Greco [14] believed that the sliding surface will reflect and form a multi fold line failure mode after it develops to the adjacent wall. Yang et al. [15] conducted a series of model tests and observed that the slip surfaces are curvilinear planes developed from the heel of the retaining wall to the crest of the backfill.
On the other hand, after the soil arch theory was put forward by Janssen [16], the relevant earth pressure research was carried out by many geotechnical scholars. Handy [17] defined the active soil arch as the trajectory of small principal stress, assumed that the trajectory of the soil arch between two parallel retaining walls is catenary, and established the expression of the earth pressure distribution behind the wall by the slice method. Paik and Salgado [18] assumed that the soil arch is a circular arc, and the sliding surface of the backfill is the Rankine sliding surface, but the deflection of the principal stress at the sliding surface is not considered. Goel [19] assumed the soil arch as a parabola and solved the theoretical expression of active earth pressure in the form of polar coordinates in the case of the linear failure mode and parabola failure mode. Many subsequent studies have proven that the soil arching effect cannot be ignored in the calculation of earth pressure [9,20,21,22].
Additionally, due to the complex natural geological environment, various retaining structures are inevitably backfilled by unsaturated soil. In this case, the traditional earth pressure theories treat backfills as completely saturated or dry soils. Nevertheless, many engineering cases demonstrate that economic and safe retaining structure design can only be carried out if the lateral pressure change caused by the change of matric suction in unsaturated backfill is reasonably considered [23,24,25]. The following work has been carried out by researchers: Pufahl et al. [26] obtained the Rankine earth pressure solution of unsaturated soil based on the strength theory of two stress variables [27]. Liang et al. [28] combined the strength theory of unsaturated soil with the Coulomb earth pressure theory, and obtained the unified solution of Coulomb earth pressure in unsaturated soil. Since then, the limit equilibrium method has been widely used to solve the earth pressure of unsaturated soil [24,29,30]. Moreover, the slip line theory [31,32] and the upper bound method of limit analysis [33,34,35] have also been extended to the calculation of earth pressure in unsaturated soils. However, most unsaturated earth pressure theories still assume that the backfill is in semi-infinite soil conditions when analyzing the finite fill behind the wall, which is seriously inconsistent with the actual situation, resulting in too conservative calculation results. In addition, two ideas are generally adopted for the theoretical derivation of unsaturated earth pressure at present. One is based on the traditional Coulomb earth pressure method, and the overall sliding wedge is used to solve the problem, but the distribution of earth pressure cannot be obtained. The other method is to use the horizontal thin-layer slicing method, but most theories do not consider the comprehensive influence of the soil arching effect and interlayer shear stress in the analysis, resulting in a large deviation between the predicted results and the actual situation.
Based on the above analysis, this study takes the narrow unsaturated backfill as the analysis object, assumes that the trajectory of small principal stress is an arc and the soil slip surface behind the wall is a Coulomb slip surface, and the Mohr stress circle is used to calculate the average vertical stress and shear stress of the horizontal interlayer; on this basis, the expression of active earth pressure, its resultant force and the height of the action point is given by the horizontal thin layer analysis method, and it is compared with the model test and other theoretical results in the relevant literature. At the same time, through parameter analysis, the effects of fill aspect ratio, cohesion, internal friction angle, matrix suction and wall–soil interface friction angle on active earth pressure, its resultant force and the height of the action point are studied.

2. Shear Strength for Unsaturated Soils

In order to more clearly describe the characteristic that the shear strength of unsaturated soil increases with the increase in matric suction, Fredlund et al. [27] introduced a new material variable φb related to matric suction, and proposed a shear strength theory of two independent stress state variables for calculating the shear strength of unsaturated soil. As shown in Figure 1, in the three-dimensional space of two stress variables and shear stress, the shear strength envelope is a plane. The expression of shear strength is given by:
τ f = σ u a tan φ + c + u a u w tan φ b
where σ is the total normal stress; ua and uw are the pore air pressure and the pore water pressure; (σua) represents the net normal stress on the failure plane; φ’ is the effective internal friction angle; c’ is the effective cohesion; (uauw) represents the matric suction; φb is matric suction angle, which indicates the rate of increase in the soil strength with respect to an increase in matric suction.
Subsequent studies [36,37] found that the ϕb is not a fixed value. With the change of matric suction, φb will take different values. Based on the microscopic analysis of unsaturated soil, Vanapalli et al. [38] suggested that the SWCC model be used to correlate the relationship between ϕb and volumetric water content of soil, which has the following form:
tan φ b = θ θ r θ s θ r tan φ
where θ is the volumetric water content; θr is the residual volumetric water content; and θs is the saturated volumetric water content. Using the SWCC model proposed by Van Genuchten [39], the matric suction angle corresponding to different matric suction can be obtained from the following formula:
tan φ b = tan φ 1 + α u a u w n 1 1 / n
where α and n are SWCC fitting parameters; parameter α is approximate to the inverse of air entry suction; parameter n reflects the pore size distribution.
Through the above analysis, when the value of matric suction is determined, the shear strength provided by the matric suction can be calculated. The intercept of the shear strength envelope in the corresponding net stress and shear stress plane coordinate system is the sum of the shear strength provided by the effective cohesion and matric suction, as shown in Figure 1. In the subsequent calculation, we will express the contribution of matric suction to shear strength through total cohesion:
c t = c + u a u w tan φ b

3. Lateral Pressure Coefficient and Shear Stress Coefficient

Terzaghi [40] defines the phenomenon of earth pressure transfer from the yield zone to the adjacent static zone as a soil arching effect, and the existence of the soil arching effect changes the distribution of the earth pressure. The basic idea of calculating the earth pressure by soil arching effect is to obtain the active lateral earth pressure coefficient by assuming the geometric shape of the principal stress arch, and then solve the average vertical pressure and earth pressure of the soil layer according to the differential element force balance equation. In the calculation of the lateral earth pressure, the circular earth arch is widely used and the earth pressure distribution is in good agreement with the test results [18,19,41]. However, few studies have considered the influence of adjacent underground buildings on the soil arch.
As shown in Figure 2, in the case of narrow backfill, assuming that the backs of two rigid walls are vertical, and the backfill surface is horizontal, it is considered that there is an overload q0 on the surface of the backfill, and the backfill of the wall is in an unsaturated state. Suppose the wall height is H, the soil weight is γ, the soil cohesion is c’, the soil friction angle is φ’, the friction angle between the retaining wall and the soil is δ1, and the friction between the existing building and the soil is angle is δ2. The active slip surface intersects with the adjacent underground buildings to form a trapezoidal soil wedge, the Coulomb slip surface is assumed to be the slip surface under this working condition, and the inclination angle β of the Coulomb-type slip surface passing through the wall heel is:
β = arctan tan 2 φ + tan φ tan φ + δ 1 + tan φ
Different lateral earth pressure coefficients and interlayer shear stress coefficients will be obtained when the soil arching effect is analyzed in the upper rectangular region (zone I) and the lower triangular region (zone II). In this section, the lateral earth pressure coefficient and the interlayer shear stress coefficient of each region are solved first, and the influence of the wall soil friction angle on them is discussed. Then, in the next section, the distribution of earth pressure, resultant force and their action points in each region are solved according to the horizontal thin layer element method.
Considering the total cohesion of unsaturated soils, the coordinate system σ is obtained by translating the coordinate system σoτ to the left by ctcotϕ’, as shown in Figure 3. Next, the lateral earth pressure coefficient and interlayer shear stress coefficient will be calculated by analyzing the soil arching effect in different areas.
Zone I:
From the stress Mohr’s circle in Figure 3, the angles θA and θB between the large principal stress and the horizontal direction at points A and B in Figure 2 can be expressed as:
θ A = π 2 1 2 arcsin sin δ 1 sin φ + δ 1 2
θ B = π 2 1 2 arcsin sin δ 2 sin φ + δ 2 2
where δ1 and δ2 are the interface friction angle of the moving wall and the fixed wall.
The width of the horizontal thin layer element in zone I can be calculated as:
b = R cos θ A + cos θ B
where R is the radius of circular soil arch.
In the new coordinate system, the horizontal stress, vertical stress and shear stress at any point in the soil arch can be expressed as:
σ x = σ 1 cos 2 θ + K a sin 2 θ
σ z = σ 1 sin 2 θ + K a cos 2 θ
τ = σ 1 1 K a sin θ cos θ
where σ1 is the major principal stress of the new coordinate system; Ka is the Rankine active earth pressure coefficient, which can be calculated by K a = 1 sin φ 1 + sin φ .
Handy [17] and Paik [18] defined the lateral active earth pressure coefficient as the ratio of the horizontal stress behind the retaining wall to the average value of the vertical stress on the soil arch trajectory. Similarly, we define the interlayer shear stress coefficient as the ratio of the average shear stress on the soil arch trajectory to the average vertical stress on the soil arch trajectory.
The average shear stress and the average vertical stress on the soil arch trajectory can be obtained by the integral method:
σ ¯ z 1 = θ A π θ B σ 1 sin 2 θ + K a cos 2 θ R sin θ d θ b = σ 1 1 + K a 1 cos 3 θ A + cos 3 θ B 3 cos θ A + cos θ B
τ ¯ 1 = θ A π θ B σ 1 1 K a sin θ cos θ R sin θ d θ b = σ 1 1 K a sin 3 θ B sin 3 θ A 3 cos θ A + cos θ B
The horizontal reaction force of the retaining wall in the new coordinate system can be expressed as:
σ x A = σ 1 cos 2 θ + A K a sin 2 θ A
Therefore, in the new coordinate system, the lateral active earth pressure coefficient and interlayer shear stress coefficient are:
K 1 = σ x A σ ¯ z 1 = 3 cos θ A + cos θ B cos 2 θ + A K a sin 2 θ A 3 cos θ A + cos θ B + K a 1 cos 3 θ A + cos 3 θ B
k 1 = τ ¯ 1 σ ¯ z 1 = 1 K a sin 3 θ B sin 3 θ A 3 cos θ A + cos θ B + K a 1 cos 3 θ A + cos 3 θ B
Zone II:
The deflection angle of point C has the same expression as that of point A. Considering the deflection of the principal stress at the sliding surface, θC and θD is obtained, as follows:
θ C = π 2 1 2 arcsin sin δ 1 sin φ + δ 1 2
θ D = π 4 + β φ 2
The width of the horizontal element in zone II and the radius of the soil arch conform to the following relationship:
b z = R cos θ C cos θ D
Referring to Equations (12)–(14), the average vertical stress, average shear stress and horizontal reaction force of the retaining wall in zone II can be expressed as:
σ ¯ z 2 = θ C θ D σ 1 sin 2 θ + K a cos 2 θ R sin θ d θ b z = σ 1 1 + K a 1 cos 3 θ C cos 3 θ D 3 cos θ C cos θ D
τ ¯ 2 = θ C θ D σ 1 1 K a sin θ cos θ R sin θ d θ b z = σ 1 1 K a sin 3 θ D sin 3 θ C 3 cos θ C cos θ D
σ x C = σ 1 cos 2 θ + C K a sin 2 θ C
Then, the lateral pressure coefficient and interlayer shear stress coefficient of zone II in the new coordinate system can be expressed as:
K 2 = σ x C σ ¯ z 2 = 3 cos θ C cos θ D cos 2 θ + C K a sin 2 θ C 3 cos θ C cos θ D + K a 1 cos 3 θ C cos 3 θ D
k 2 = τ ¯ 2 σ ¯ z 2 = 1 K a sin 3 θ D sin 3 θ C 3 cos θ C cos θ D + K a 1 cos 3 θ C cos 3 θ D
Figure 4 shows the relationship between the interlayer shear stress coefficients (k1 and k2) and the interface friction angle of the fixed wall–soil δ2 under different values of the interface friction angle of the moving wall–soil δ1. The curves in the figures are obtained from Equations (16) and (24). As shown in the Figure 4a, the increase of δ1, k1 increases obviously, but decreases with the increase of δ2. Moreover, from Figure 4b, k2 increases with the increase of δ1, but does not change with the increase of δ2. This is because the soil arch in zone II is not formed to the fixed wall. Additionally, when δ1 = δ2, k1 = 0. As δ2 continues to increase, k1 will take a negative value, and the direction of the interlayer shear stress will change.

4. Derivation of Active Earth Pressure Formula

For unsaturated backfills, the tension cracks on the backfill surface cannot be ignored due to matric suction. Combining Equation (9) and the transformation relationship between the two coordinate systems, at the critical depth (z = hc), the following relationship is given:
σ x A = σ 1 cos 2 θ + A K a sin 2 θ A c t cot φ = γ h c + q 0 + c t cot φ cos 2 θ + A K a sin 2 θ A c t cot φ = 0
Rewriting the above equation, the tension crack depth can be expressed as:
h c = 1 γ c t cot φ cos 2 θ + A K a sin 2 θ A c t cot φ q 0
In Equation (26), when hc < 0, take hc = 0.
Zone I:
Carry out force balance analysis for rectangular thin layer element, as shown in Figure 5, in which σx1 and σx2 are the horizontal reaction force of the retaining wall and fixed wall; τw1 and τw2 are the interface shear stress of the retaining wall and fixed wall; σ ¯ z and σ ¯ z + d σ ¯ z are the average vertical stress acting on the upper surface and lower surface; τ ¯ and τ ¯ + d τ ¯ are the average interlayer shear stress acting on the upper surface and lower surface; dW is the self-weight of thin layer element.
The force equilibrium equations in the vertical and horizontal directions are established as:
σ x 1 d z σ x 2 d z + τ ¯ + d τ ¯ b τ ¯ b = 0
τ w 1 d z + τ w 2 d z σ ¯ z b + σ ¯ z + d σ ¯ z b d W = 0
where
τ w 1 = σ x 1 tan δ 1 + c w t 1 τ w 2 = σ x 2 tan δ 1 + c w t 2 σ x 1 = K 1 σ ¯ z + K 1 1 c t cot φ τ ¯ = k 1 σ ¯ z + k 1 c t cot φ d W = γ b d z
The adhesive force of wall and soil can be considered as the same as the friction force, which is given by:
c w t 1 = c t tan δ 1 tan φ ,   c w t 2 = c t tan δ 2 tan φ
Neglecting the second-order terms, the first-order differential equation can be obtained by solving Equations (27) and (28):
d σ ¯ z d z + M 1 σ ¯ z + N 1 = 0
where
M 1 = K 1 tan δ 1 + tan δ 2 b k 1 tan δ 2 + 1 N 1 = K 1 c t cot φ tan δ 1 + tan δ 2 γ b b k 1 tan δ 2 + 1
By integrating Equation (30), the average vertical stress can be solved as:
σ ¯ z = C 1 e M 1 z N 1 M 1
Substituting the boundary conditions ( z = h c , σ ¯ z = q 0 + γ h c ) into Equation (33), C1 is obtained as follows:
C 1 = e M 1 h c q 0 + γ h c + N 1 M 1
Zone II:
As shown in Figure 6, according to the equilibrium conditions of stresses acting on the trapezoidal thin layer element, the two equations are expressed as:
σ x d z σ n d z sin β sin β + τ s d z cos β sin β + τ ¯ + d τ ¯ b z d z cot β τ ¯ b z = 0
τ w d z + τ s d z sin β sin β + σ n d z cos β sin β + σ ¯ z + d σ ¯ z b z d z cot β σ ¯ z b z d W = 0
in which
τ w = σ x tan δ 1 + c w t 1 τ s = σ n tan φ + c t σ x = K 2 σ ¯ z + K 2 1 c t cot φ τ ¯ = k 2 σ ¯ z + k 2 c t cot φ b z = H z cot β d W = γ b z d z 1 2 γ d z cot β d z
where σn and τs are the normal stress and shear stress on the slip surface.
Combining Equations (35) and (36), a differential equation is obtained as:
d σ ¯ z d z + M 2 H z σ ¯ z + N 1 H z + P 2 = 0
where
M 2 = K 2 tan π 2 β + φ + tan δ 1 cot β 1 + k 2 tan π 2 β + φ 1 N 2 = c t tan π 2 β + φ ( K 2 1 ) cot φ + cot β 1 k 2 cot φ + K 2 cot φ tan δ 1 + 1 cot β 1 + k 2 tan π 2 β + φ P 2 = γ 1 + k 2 tan π 2 β + φ
Integrating Equation (38) can be solved as
σ ¯ z = C 2 H z M 2 + P 2 1 M 2 H z N 2 M 2
By substituting the boundary value z = h , σ ¯ z = C 1 e M 1 h N 1 M 1 , C2 is obtained, as follows:
C 2 = C 1 e M 1 h N 1 M 1 P 2 1 M 2 H h + N 2 M 2 H h M 2
Therefore, the distribution of the horizontal active earth pressure is obtained as:
σ x = K 1 C 1 e M 1 z N 1 M 1 + K 1 1 c t cot φ h c z h K 2 C 2 H z M 2 + P 2 1 M 2 H z N 2 M 2 + K 2 1 c t cot φ h z H
The resultant force of the horizontal active earth pressure and the total overturning moment of the retaining wall can be expressed, respectively, as
E x = h c H σ x d z = h c h K 1 C 1 e M 1 z N 1 M 1 + K 1 1 c t cot φ d z + h H K 2 C 2 H z M 2 + P 2 1 M 2 H z N 2 M 2 + K 2 1 c t cot φ d z  
M = h c H σ x H z d z = h c h K 1 C 1 e M 1 z N 1 M 1 H z + K 1 1 c t cot φ H z d z + h H K 2 C 2 H z M 2 + P 2 1 M 2 H z N 2 M 2 H z + K 2 1 c t cot φ H z d z  
It should be noted that a special working condition may also occur. When h < hc, there will be no zone I. In this case, the boundary conditions ( z = h c , σ ¯ z = q 0 + γ h c ) are substituted into Equation (40), and C2 is obtained as follows:
C 2 = q 0 + γ h c P 2 1 M 2 H h c + N 2 M 2 H h c M 2
Then, the horizontal active earth pressure, resultant force of horizontal active earth pressure and the total overturning moment of the retaining wall can be expressed, respectively, as
σ x = 0 0 z h c K 2 C 2 H z M 2 + P 2 1 M 2 H z N 2 M 2 + K 2 1 c t cot φ h c z H
E x = h c H σ x d z = h c H K 2 C 2 H z M 2 + P 2 1 M 2 H z N 2 M 2 + K 2 1 c t cot φ d z  
M = h c H σ x H z d z = h c H K 2 C 2 H z M 2 + P 2 1 M 2 H z N 2 M 2 H z + K 2 1 c t cot φ H z d z  
The distance from the location of the resultant of the active earth pressure to the heel of the wall, denoted as Ha, can be easily derived as the ratio of M to Ex:
H a = M E x

5. Comparison and Verification

The novelty of this study is to propose a theoretical framework for calculating the active earth pressure of narrowed unsaturated backfill. To the knowledge of the authors, no similar experimental and theoretical methods have been developed so far, so the proposed method is only compared with the results in the absence of the matric suction. Frydman et al. [6] carried out centrifugal tests of limited sand active earth pressure for vertical retaining walls under different aspect ratio (b/H) conditions. In the test, b/H = 0.235, γ = 15.8 kN/m3, φ’ = 36°, δ1 = δ2 = 25°, and the backfill width b = 2 m. Additionally, Chen et al. [13] used the horizontal differential element method to calculate the σx/(γH) value of the narrow backfill behind the wall. Figure 7 shows the comparison results of the normalized active earth pressure σx/(γH) value obtained by the method in this paper and the centrifugal test [6] and Chen et al. [13]. It can be observed that the method in this paper is in good agreement with the results of the centrifugation test, which demonstrates the rationality and accuracy of the method in this paper. In addition, the horizontal differential element method adopted by Chen et al. [13] ignores the shear force between the adjacent elements; thus, the calculation result is larger than that of the method in this study.

6. Parametric Study

In this section, a parametric study was performed to analyze the effective internal friction angle φ’, surcharge pressure q0, interface friction angle δ, effective cohesion c’, matric suction (uauw), and width-height ratio of backfill b/H on the distribution of the horizontal active earth pressure. In the analysis, these parameters were fixed: γ = 18 kN/m3, α = 0.02 kPa−1, n = 3 and H = 10 m. Other parameters for analysis can be found in the analysis figures. In addition, the interface friction angle δ1 of the moving wall–soil and the interface friction angle δ2 of the fixed wall–soil are equal, that is, δ1 = δ2 = δ.
Figure 8 depicts the normalized horizontal active earth pressure (σx/γH) along the normalized height (z/H) for the various effective internal friction angle φ’. To investigate the effect of the effective internal friction angle on the horizontal active earth pressure against the retaining structure, different values of φ’ (i.e., φ’ = 20°, 25°, 30°, 35° and 40°) were used. From Figure 8, it is obvious that the horizontal earth pressure first increases and then decreases along the depth, showing a nonlinear “drum” distribution, reaching a peak above the wall bottom. With the effective internal friction angle increases, the active earth pressure acting on the retaining wall at each depth decreases significantly, and the area where the earth pressure is zero gradually increases, indicating that the depth of the tension crack also increases with the increase in the effective internal friction angle.
Figure 9 shows the normalized horizontal active earth pressure (σx/γH) along the normalized height (z/H) for different values of the surcharge pressure q0. Herein, the aim is to discuss the effect of the surcharge pressure on the horizontal active earth pressure against the retaining structure, and thus different constant values of q0 (i.e., q0 = 10 kPa, 20 kPa, 30 kPa, 40 kPa and 50 kPa) were used. Figure 9 shows that the horizontal active earth pressure increases significantly with the increase in surcharge pressure applied on the backfill surface, whereas the shape of the horizontal active earth pressure distribution remains unchanged. Moreover, the depth of the tension crack decreases with the increase in the surcharge pressure, and when values of q0 are large enough, the tension crack will not form.
Figure 10 illustrates the normalized horizontal active pressure (σx/γH) along the normalized height (z/H) for the various interface friction angle δ. To discuss the effect of interface friction angle on the horizontal active earth pressure against the retaining structure, different values of δ varying from 0.2φ’ to 0.8φ’ were used here. As shown in Figure 10, with the increase in the interface friction angle, the horizontal active earth pressure in the upper part of the retaining wall increases slightly, while the horizontal active earth pressure in the lower part decreases obviously. Furthermore, the action point of the horizontal earth pressure resultant force has a tendency to move upward with the increase in interface friction angle.
Figure 11 shows the normalized horizontal active earth pressure (σx/γH) along the normalized height (z/H) for the various values of soil effective cohesion c’. To analyze the impact of the backfill inclination angle on the horizontal active earth pressure against the retaining wall, different values of c’ (i.e., c’ = 5 kPa, 10 kPa, 15 kPa, 20 kPa and 25 kPa) were used. Figure 11 shows a pronounced decrease in the horizontal active earth pressure against the rigid retaining wall with the increasing value of c’, and the depth of the tension crack in soil increases correspondingly. From the earth pressure distribution curves, with the increase in the effective cohesion of the soil, the earth pressure distribution curves gradually shift to the left axis.
Figure 12 depicts the normalized horizontal active earth pressure distribution (σx/γH) along the normalized height (z/H) for various matric suction (uauw). To study the effect of the matric suction (uauw) on the horizontal active earth pressure against the retaining wall, various values of (uauw) varying from 0 kPa to 200 kPa were used. From Figure 11, with the increase in matric suction, the horizontal active earth pressure first decreases rapidly and then increases gradually, and the depth of the tension crack first increases obviously and then decreases, while the shape of the horizontal active earth pressure distribution remains unchanged. The reason is that the total cohesion increases first and then decreases with the increase in matric suction, and this is consistent with the conclusion given by Song [42]. the horizontal active earth pressure reaches its minimum value at the matric suction equaling 50 kPa; this is because the matric suction takes the air entry pressure (approximately equal to the reciprocal of parameter α), and the apparent cohesion provided by the matric suction is the largest.
The normalized horizontal active earth pressure distribution (σx/γH) along the normalized height (z/H) for different width-height ratio of backfill b/H is shown in Figure 13. In order to study the effect of the width of backfill on the horizontal active earth pressure against the retaining wall, various values of b/H (i.e., b/H = 0.2, 0.3, 0.4, 0.6, and 0.8) were employed. As shown in Figure 13, the magnitude of the horizontal active earth pressure decreases with the decrease in the b/H value, and the attenuation becomes clearer when the b/H value is smaller. When the width of backfill reaches a certain value, the horizontal active earth pressure remains unchanged with the increase in the b/H value. This is because the backfill with a sufficient width allows the slip surface to fully develop.

7. Conclusions

The variation of soil suction has a significant effect on the active earth pressure of narrow unsaturated backfill in retaining structures. An analytical framework was developed for evaluating the active earth pressures of narrow unsaturated backfills affected by arching when the interlayer shear stress was considered. By comparing the results of the proposed method with other theoretical and experimental results, the method’s effectiveness was verified. A detailed parametric study was carried out to discuss the effect of various factors on active earth pressure. The main conclusions include the following:
(1)
With the increase in the interface friction angle of the moving wall–soil, the average shear stress coefficient of zone I and zone II increases gradually, but with the increase in the interface friction angle of the fixed wall–soil, the average shear stress coefficient of zone I decreases, while the average shear stress coefficient of zone II remains unchanged. The horizontal earth pressure shows a nonlinear “drum” distribution along the depth, which increases first and then decreases, and reaches a peak above the wall bottom.
(2)
With the increase in effective internal friction angle and effective cohesion, the horizontal active earth pressure decreases significantly, while the tensile depth increases. In addition, the horizontal active earth pressure increases significantly with the increase in surcharge pressure, but the distribution shape of the horizontal active earth pressure remains unchanged. Additionally, when the interface friction angle increases, the horizontal active earth pressure at the upper part of the retaining wall slightly increases, while the horizontal active earth pressure at the lower part obviously decreases. In addition, the action point of the horizontal earth pressure force tends to move up with the increase in the friction angle.
(3)
With the increase in matric suction, the horizontal active earth pressure first decreases rapidly and then increases gradually, and the tension crack depth first increases obviously and then decreases, but the distribution pattern of the horizontal active earth pressure remains unchanged.
(4)
The horizontal active earth pressure decreases with the decrease in the aspect ratio. When the aspect ratio is smaller, the attenuation is more obvious. When the aspect ratio reaches a certain value, the horizontal active earth pressure is basically unchanged.

Author Contributions

Formal analysis, M.Y. and H.G.; Funding acquisition, M.Y. and B.D.; Investigation, M.Y.; Methodology, H.G.; Writing—original draft, D.F. and B.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Transportation Science and Technology project of Henan province (2019J-2-12, 2021J7) and Science Foundation for Youths of Hunan Province of China (2021JJ40460), which are gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tsagareli, Z.V. Experimental investigation of the pressure of a loose medium on retaining walls with a vertical back face and horizontal backfill surface. Soil Mech. Found. Eng. 1965, 2, 197–200. [Google Scholar] [CrossRef]
  2. Fang, Y.S.; Ishibashi, I. Static earth pressures with various wall movements. J. Geotech. Eng. 1986, 112, 317–333. [Google Scholar] [CrossRef]
  3. Take, W.A.; Valsangkar, A.J. Earth pressures on unyielding retaining walls of narrow backfill width. Can. Geotech. J. 2001, 38, 1220–1230. [Google Scholar] [CrossRef]
  4. O’Neal, T.S.; Hagerty, D.J. Earth pressures in confined cohesionless backfill against tall rigid walls—A case history. Can. Geotech. J. 2011, 48, 1188–1197. [Google Scholar] [CrossRef]
  5. Yang, M.; Tang, X. Rigid retaining walls with narrow cohesionless backfills under various wall movement modes. Int. J. Geomech. 2017, 17, 04017098. [Google Scholar] [CrossRef]
  6. Frydman, S.; Keissar, I. Earth pressure on retaining walls near rock faces. J. Geotech. Eng. 1987, 113, 586–599. [Google Scholar] [CrossRef]
  7. Tang, Y.; Chen, J. A computational method of active earth pressure from finite soil body. Math. Probl. Eng. 2018, 2018, 9892376. [Google Scholar] [CrossRef] [Green Version]
  8. Yang, M.; Deng, B. Simplified method for calculating the active earth pressure on retaining walls of narrow backfill width based on DEM analysis. Adv. Civ. Eng. 2019, 2019, 1507825. [Google Scholar] [CrossRef] [Green Version]
  9. Xie, M.; Zheng, J.; Zhang, R.; Cui, L.; Miao, C. Active earth pressure on rigid retaining walls built near rock faces. Int. J. Geomech. 2020, 20, 04020061. [Google Scholar] [CrossRef]
  10. Liu, H.; Kong, D. Active earth pressure of finite width soil considering intermediate principal stress and soil arching effects. Int. J. Geomech. 2022, 22, 04021294. [Google Scholar] [CrossRef]
  11. Wang, Z.; Liu, X.; Wang, W.; Tao, Z.; Li, S. Inclined layer method-based theoretical calculation of active earth pressure of a finite-width soil for a rotating-base retaining wall. Sustainability 2022, 14, 9772. [Google Scholar] [CrossRef]
  12. Hu, H.; Yang, F.; Tang, H.; Zeng, Y.; Zhou, J.; Gong, X. Field Study on Earth Pressure of Finite Soil Considering Soil Displacement. Appl. Sci. 2022, 12, 8059. [Google Scholar] [CrossRef]
  13. Chen, J.J.; Li, M.G.; Wang, J.H. Active earth pressure against rigid retaining walls subjected to confined cohesionless soil. Int. J. Geomech. 2017, 17, 06016041. [Google Scholar] [CrossRef]
  14. Greco, V. Active thrust on retaining walls of narrow backfill width. Comput. Geotech. 2013, 50, 66–78. [Google Scholar] [CrossRef]
  15. Yang, M.; Tang, X.; Wu, Z. Slip surface and active earth pressure of cohesionless narrow backfill behind rigid retaining walls under translation movement mode. Int. J. Geomech. 2020, 20, 04020115. [Google Scholar] [CrossRef]
  16. Janssen, H.A. Versuche uber getreidedruck in silozellen. Z. Ver. Dtsch. Ing. 1895, 39, 1045–1049. [Google Scholar]
  17. Handy, R.L. The arch in soil arching. J. Geotech. Eng. 1985, 111, 302–318. [Google Scholar] [CrossRef]
  18. Paik, K.H.; Salgado, R. Estimation of active earth pressure against rigid retaining walls considering arching effects. Geotechnique 2003, 53, 643–653. [Google Scholar] [CrossRef]
  19. Goel, S.; Patra, N.R. Effect of arching on active earth pressure for rigid retaining walls considering translation mode. Int. J. Geomech. 2008, 8, 123–133. [Google Scholar] [CrossRef]
  20. Nadukuru, S.S.; Michalowski, R.L. Arching in distribution of active load on retaining walls. J. Geotech. Geoenviron. Eng. 2012, 138, 575–584. [Google Scholar] [CrossRef]
  21. Khosravi, M.H.; Pipatpongsa, T.; Takemura, J. Theoretical analysis of earth pressure against rigid retaining walls under translation mode. Soils Found. 2016, 56, 664–675. [Google Scholar] [CrossRef]
  22. Patel, S.; Deb, K. Study of active earth pressure behind a vertical retaining wall subjected to rotation about the base. Int. J. Geomech. 2020, 20, 04020028. [Google Scholar] [CrossRef]
  23. Yoo, C.; Jung, H.Y. Case history of geosynthetic reinforced segmental retaining wall failure. J. Geotech. Geoenviron. Eng. 2006, 132, 1538–1548. [Google Scholar] [CrossRef]
  24. Vahedifard, F.; Leshchinsky, B.A.; Mortezaei, K.; Lu, N. Active earth pressures for unsaturated retaining structures. J. Geotech. Geoenviron. Eng. 2015, 141, 04015048. [Google Scholar] [CrossRef]
  25. Tang, H.; Wasowski, J.; Juang, C.H. Geohazards in the three Gorges Reservoir Area, China–Lessons learned from decades of research. Eng. Geol. 2019, 261, 105267. [Google Scholar] [CrossRef]
  26. Pufahl, D.E.; Fredlund, D.G.; Rahardjo, H. Lateral earth pressures in expansive clay soils. Can. Geotech. J. 1983, 20, 228–241. [Google Scholar] [CrossRef]
  27. Fredlund, D.G.; Morgenstern, N.R.; Widger, R.A. The shear strength of unsaturated soils. Can. Geotech. J. 1978, 15, 313–321. [Google Scholar] [CrossRef]
  28. Liang, W.B.; Zhao, J.H.; Li, Y.; Zhang, C.G.; Wang, S. Unified solution of Coulomb’s active earth pressure for unsaturated soils without crack. In Applied Mechanics and Materials; Trans Tech Publications Ltd: Bäch, Switzerland, 2012; Volume 170, pp. 755–761. [Google Scholar]
  29. Deng, B.; Yang, M. Analysis of passive earth pressure for unsaturated retaining structures. Int. J. Geomech. 2019, 19, 06019016. [Google Scholar] [CrossRef]
  30. Abdollahi, M.; Vahedifard, F.; Abed, M.; Leshchinsky, B.A. Effect of tension crack formation on active earth pressure encountered in unsaturated retaining wall backfills. J. Geotech. Geoenviron. Eng. 2021, 147, 06020028. [Google Scholar] [CrossRef]
  31. Vo, T.; Russell, A.R. Slip line theory applied to a retaining wall–unsaturated soil interaction problem. Comput. Geotech. 2014, 55, 416–428. [Google Scholar] [CrossRef]
  32. Vo, T.; Russell, A.R. Interaction between retaining walls and unsaturated soils in experiments and using slip line theory. J. Eng. Mech. 2017, 143, 04016120. [Google Scholar] [CrossRef]
  33. Li, Z.W.; Yang, X.L. Active earth pressure for soils with tension cracks under steady unsaturated flow conditions. Can. Geotech. J. 2018, 55, 1850–1859. [Google Scholar] [CrossRef]
  34. Li, Z.; Yang, X. Three–dimensional active earth pressure for retaining structures in soils subjected to steady unsaturated seepage effects. Acta Geotech. 2020, 15, 2017–2029. [Google Scholar] [CrossRef]
  35. Chehade, H.A.; Dias, D.; Sadek, M.; Jenck, O.; Chehade, F. Upper bound seismic limit analysis of geosynthetic–reinforced unsaturated soil walls. Geotext. Geomembr. 2020, 48, 419–430. [Google Scholar] [CrossRef]
  36. Escario, V.; Saez, J. The shear strength of partly saturated soils. Geotechnique 1986, 36, 453–456. [Google Scholar] [CrossRef]
  37. Zhang, L.L.; Fredlund, D.G.; Fredlund, M.D.; Wilson, G.W. Modeling the unsaturated soil zone in slope stability analysis. Can. Geotech. J. 2014, 51, 1384–1398. [Google Scholar] [CrossRef]
  38. Vanapalli, S.K.; Fredlund, D.G.; Pufahl, D.E.; Clifton, A.W. Model for the prediction of shear strength with respect to soil suction. Can. Geotech. J. 1996, 33, 379–392. [Google Scholar] [CrossRef]
  39. Van Genuchten, T.M. A closed–form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef] [Green Version]
  40. Terzaghi, K. Theoretical Soil Mechanics; John Wiley and Sons: Hoboken, NJ, USA, 1943. [Google Scholar]
  41. Rao, P.; Chen, Q.; Zhou, Y.; Nimbalkar, S.; Chiaro, G. Determination of active earth pressure on rigid retaining wall considering arching effect in cohesive backfill soil. Int. J. Geomech. 2016, 16, 04015082. [Google Scholar] [CrossRef]
  42. Song, Y.S. Suction stress in unsaturated sand at different relative densities. Eng. Geol. 2014, 176, 1–10. [Google Scholar] [CrossRef]
Figure 1. Extended Mohr-Coulomb failure envelope of unsaturated soils.
Figure 1. Extended Mohr-Coulomb failure envelope of unsaturated soils.
Sustainability 14 12699 g001
Figure 2. Narrow backfill behind retaining wall and its trajectory of minor principal stress.
Figure 2. Narrow backfill behind retaining wall and its trajectory of minor principal stress.
Sustainability 14 12699 g002
Figure 3. Stress state behind retaining wall. (a) Mohr stress circle; (b) strss state.
Figure 3. Stress state behind retaining wall. (a) Mohr stress circle; (b) strss state.
Sustainability 14 12699 g003
Figure 4. Relationship between interlayer shear stress coefficients and δ2: (a) k1; (b) k2.
Figure 4. Relationship between interlayer shear stress coefficients and δ2: (a) k1; (b) k2.
Sustainability 14 12699 g004
Figure 5. Horizontal thin-layer element model of finite soil of element I.
Figure 5. Horizontal thin-layer element model of finite soil of element I.
Sustainability 14 12699 g005
Figure 6. Horizontal thin-layer element model of finite soil of element II.
Figure 6. Horizontal thin-layer element model of finite soil of element II.
Sustainability 14 12699 g006
Figure 7. Comparison of centrifugal test and theoretical calculation results [6,13].
Figure 7. Comparison of centrifugal test and theoretical calculation results [6,13].
Sustainability 14 12699 g007
Figure 8. Values of σx/γH versus z/H for different values of the effective internal friction angle φ’.
Figure 8. Values of σx/γH versus z/H for different values of the effective internal friction angle φ’.
Sustainability 14 12699 g008
Figure 9. Values of σx/γH versus z/H for different values of the surcharge pressure q0.
Figure 9. Values of σx/γH versus z/H for different values of the surcharge pressure q0.
Sustainability 14 12699 g009
Figure 10. Values of σx/γH versus z/H for different values of the interface friction angle δ.
Figure 10. Values of σx/γH versus z/H for different values of the interface friction angle δ.
Sustainability 14 12699 g010
Figure 11. Values of σx/γH versus z/H for different values of the effective cohesion c’.
Figure 11. Values of σx/γH versus z/H for different values of the effective cohesion c’.
Sustainability 14 12699 g011
Figure 12. Values of σx/γH versus z/H for different values of the matric suction (uauw).
Figure 12. Values of σx/γH versus z/H for different values of the matric suction (uauw).
Sustainability 14 12699 g012
Figure 13. Values of σx/γH versus z/H for different values of the width–height ratio b/H.
Figure 13. Values of σx/γH versus z/H for different values of the width–height ratio b/H.
Sustainability 14 12699 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fu, D.; Yang, M.; Deng, B.; Gong, H. Estimation of Active Earth Pressure for Narrow Unsaturated Backfills Considering Soil Arching Effect and Interlayer Shear Stress. Sustainability 2022, 14, 12699. https://doi.org/10.3390/su141912699

AMA Style

Fu D, Yang M, Deng B, Gong H. Estimation of Active Earth Pressure for Narrow Unsaturated Backfills Considering Soil Arching Effect and Interlayer Shear Stress. Sustainability. 2022; 14(19):12699. https://doi.org/10.3390/su141912699

Chicago/Turabian Style

Fu, Daxi, Minghui Yang, Bo Deng, and Hutao Gong. 2022. "Estimation of Active Earth Pressure for Narrow Unsaturated Backfills Considering Soil Arching Effect and Interlayer Shear Stress" Sustainability 14, no. 19: 12699. https://doi.org/10.3390/su141912699

APA Style

Fu, D., Yang, M., Deng, B., & Gong, H. (2022). Estimation of Active Earth Pressure for Narrow Unsaturated Backfills Considering Soil Arching Effect and Interlayer Shear Stress. Sustainability, 14(19), 12699. https://doi.org/10.3390/su141912699

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