Experimental Studies and Analytical Approaches on Non-limited Earth Pressures of Soilbag Reinforced Retaining Structures with Surcharge Loads

1. Introduction

Geosynthetic materials, such as geotextiles, geomembranes, and geogrids, are of-ten made from recycled or recyclable materials. Their use in civil engineering projects helps in reducing the dependency on natural resources and minimizes environmental footprints [1,2,3,4,5]. Soilbags are a form of application of geosynthetic materials in civil engineering [6]. Traditionally, due to their limited durability, soilbags are primarily employed in constructing temporary embankments during floods and for building temporary structures in post-disaster scenarios [7]. However, following aging tests conducted by Matsuoka and Liu [8], it was asserted that soilbags exhibit enhanced durability when manufactured with anti-aging agents in the woven bags. This advancement has led to the gradual adoption of these enhanced soilbags in permanent infrastructure projects. Consequently, they are now extensively utilized in reinforcing foundations [9,10,11], building retaining structures [12,13,14,15,16], stabilizing slopes [17,18], and constructing small dams [19].

The soilbag reinforced retaining structure is an important form in which soilbags applied in engineering, and several application cases and the well performance of this new type of structure have been reported [14]. To establish a well-supported design guideline for scientifically engineering these structures, Fan et al. [13] explored the distribution of earth pressures behind soilbag reinforced retaining walls through laboratory experiments. The stability of such walls is intricately linked to the interlayer friction of the soilbags. Addressing this, Liu et al. [20] conducted comprehensive laboratory tests to examine the impact of the materials filled in the bags, the layered arrangement of the soilbags, and the influence of water on this interlayer friction. Fan et al. [13] then developed a method for analyzing the sliding stability of the wall, based on the equilibrium of earth pressure and the interface resistance of soilbags. Additionally, Wang et al. [16,21] investigated the effects of soilbag arrangements and tail length on the displacement and earth pressure of the wall. They proposed a method for estimating the wall’s safety using the upper-bound approach. These studies collectively contribute valuable insights for designing soilbag reinforced retaining structures.

In a large number of retaining structure projects, additional surcharge loads, such as those from vehicles and buildings, are often present on the surface of the backfill behind the structures. These loads negatively impact the stability of the retaining structure [22,23,24]. It is therefore essential to determine the earth pressure exerted on soilbag reinforced retaining walls when subjected to such surcharge loads. However, existing research has primarily focused on soilbag reinforced retaining walls under limit states [13,16]. In practical engineering situations, this state of limit equilibrium is often not reached, resulting in a non-limit state for the wall. This divergence is attributed to factors such as the restricted movement of soilbags and the complex behavior of soil. As a result, there is a need for a more sophisticated approach to accurately estimate the active earth pressure on soilbag reinforced retaining walls under surcharge loads in a non-limited state, to better mirror the actual conditions faced by these structures.

In this study, physical model tests were performed to investigate the earth pressure exerted on the soilbag reinforced retaining wall subjected to surcharge loads under the non-limited state. To predict the non-limited lateral earth pressure, an analytical method grounded in the force equilibrium of differential elements was formulated. The validity and effectiveness of this analytical approach were subsequently corroborated by the results derived from these model tests.

2. Experimental Investigations

2.1. Model Test

Fan et al. [13] investigated the failure mode and earth pressure distribution of the soilbags reinforced retaining wall subjected to surcharge loads in the limited state through a model test depicted in Figure 1. During this test, a notable phenomenon was observed when the surcharge stress reached 8.7 kPa: a sliding surface emerged in the backfill behind the wall, which aligned nearly as a straight line. This sliding surface formed an angle of 60 ° with the horizontal plane, closely resembling the theoretical angle of 45 degrees plus half the internal friction angle (φ/2), as illustrated in Figure 1b. Within the structure of the soilbag reinforced retaining wall, a unique stepped sliding surface was identified, which penetrated through three layers of soilbags. This ladder-like sliding surface formation is attributed to the interlayer insertion of soilbags [5]. Additionally, the upper portion of this stepped sliding surface exhibited horizontal movement. The critical height of the sliding wall (H_crit) in this scenario was measured at 0.95 m.

This present research employs the same experimental framework but focuses on assessing the earth pressure on the wall subjected to surcharge loads under non-limited conditions. In the model tests, a soilbag reinforced retaining wall measuring 80cm in length, 40 cm in width, and 120 cm in height, were constructed from two sizes of staggered soilbags, one set measuring 20 cm × 20 cm × 5 cm and the other 20 cm × 10 cm × 5 cm. The soilbags were made of woven bags with a weight of 150 g/m². Details on the main performance parameters of the woven bags are provided in Table 1. The bag is filled with river sand, which was also used as the backfill soil for the wall, with a density of 1.76 g/cm³. The maximum and minimum dry densities of the sand are 1.66 g/cm³ and 1.89 g/cm³, respectively, with an unevenness coefficient of C_u=2.0, a curvature coefficient of C_c=1.28, and an internal friction angle of 35.4°. It belongs to medium sand. A 1 m wide section of backfills lies behind the wall, with a 0.6 m wide loading plate placed at top near the wall. A centrally positioned hydraulic jack on this plate applies a consistent vertical load to the backfill, increasing the vertical load by 2 kPa each time until the wall fails. Lateral earth pressure on the soilbag reinforced retaining wall during this process are monitored using soil pressure cells, as depicted in Figure 1a.

2.2. Test Result

Figure 2 shows the variation in lateral earth pressure distribution along the wall height at different stages of surcharge stress, specifically at 2 kPa, 4 kPa, 6 kPa, and 8 kPa. Notably, during these stages, the retaining wall remains undamaged. The distribution of lateral earth pressure on the wall is characterized as non-linear. As the surcharge stress increases, there is a corresponding increase in the lateral earth pressure exerted on the wall.

3. Analytical Approaches

The assessment of earth pressure is a crucial component in studies focusing on retaining structures, as it significantly influences both the stability and the cost of retaining walls. Within this realm, the Rankine and Coulomb earth pressure theories are recognized as foundational approaches for estimating earth pressure. However, a limitation shared by these theories is their assumption of soil remaining static, disregarding any wall displacement. Through comprehensive model experiments, Tazaghi [25] demonstrated that the active earth pressure displays a nonlinear pattern when retaining walls experience either translation or rotation, thereby rendering traditional earth pressure theories inadequate. Consequently, numerous academics have developed alternative approaches for calculating earth pressure under varying displacement scenarios of retaining walls. Widely adopted methods include the differential element method [24] and the finite element numerical method [26,27,28]. The differential element method, specifically, is formulated based on conditions involving a retaining wall with a vertical back, sandy soil fill, and a horizontal surface behind the wall. This method segments the sliding wedge behind the wall into horizontal soil elements, presuming a uniform vertical stress distribution across these elements. Subsequently, the earth pressure calculation solution is derived by considering both the static and moment equilibrium of the soil elements, revealing a nonlinear soil pressure distribution. This method has been effectively applied in calculating active earth pressure for various types of retaining walls, yielding favorable results [29,30,31,32].

Originally, the differential element method was utilized predominantly for calculating the active earth pressure of retaining walls in a limit state. To extend its application to non-limit states, several scholars have refined this approach. Notably, Tang and Chen [33] contributed by developing theoretical formulations linking wall displacement to the friction angle in non-limit states. Subsequent methodologies introduced by other researchers have largely followed the foundational principles set by Tang and Chen, with variations in computing the effective friction angle under non-limit state scenarios [34].

In the conducted model test, it was observed that the active earth pressure exerted on the soilbag reinforced retaining wall under non-limit conditions displayed a nonlinear behavior. Furthermore, the conditions under which the wall operated were consistent with the prerequisites for applying the differential element method. As a result, this method was effectively employed to compute the non-linear earth pressures that were recorded in the test.

3.1. Basic Equation for Non-Limited Active Earth Pressure

In the context of analyzing a sliding wedge as an isolated entity, as depicted in Figure 3a, a differential flat element with a thickness of dz is selected from the wedge, positioned at a depth z beneath the backfill surface. In conditions where surcharge stress is absent (q=0), the pressure exerted by the backfill due to filling is characterized as static earth pressure. Conversely, when surcharge stress, which does not compromise the integrity of the retaining wall (q=q_int), is applied to the backfill, the soil’s immediate response is hindered, leading to a swift increase in lateral earth pressure, as demonstrated in Figure 3b. This escalation in lateral earth pressure consequently prompts compression in the neighboring soilbag, resulting in augmented compression deformation of the soilbag. Concurrently, there is an observable increase in the swelling deformation along with a decrease in the backfill’s lateral earth pressure. Equilibrium is attained once the lateral earth pressure equilibrates with the lateral compressive stress applied on the soilbag. At this juncture, the lateral swelling deformation of the backfill matches the lateral compressive deformation of the soilbag, as shown in Figure 3b. The relationship established upon achieving equilibrium following the application of surcharge stress q_int is as follows:

$$\begin{array}{c}{p}_{\mathrm{x}-\mathrm{b}\mathrm{f}}={\sigma }_{\mathrm{x}-\mathrm{s}\mathrm{b}}\\ {\epsilon }_{\mathrm{x}-\mathrm{b}\mathrm{f}}{B}_{b\mathrm{f}}={\epsilon }_{\mathrm{x}-\mathrm{s}\mathrm{b}}{B}_{\mathrm{s}\mathrm{b}}\end{array}$$

(1)

where p_x-bf and ε_x-bf are the lateral earth pressure and strain of the backfill, σ_x-sb and ε_x-sb are the lateral compressed stress and strain of soilbag, B_bf is the width of differential flat element, and B_sb is the width of the soilbag reinforced retaining wall.

Building upon the findings of Wang [23], the lateral earth pressure exerted by the backfill can be expressed as:

$$p_{\mathrm{x}-\mathrm{b}\mathrm{f}}=K_{\mathrm{q}}q+K_{\mathsf{\gamma }}\gamma z$$

(2)

where K_q and K_γ are the surcharge pressure coefficient and the backfill pressure coefficient, respectively. These coefficients are inherently linked to the strain experienced by the differential flat element.

The lateral compressed stress of soilbags can be expressed as:

$${\sigma }_{\mathrm{x}-\mathrm{s}\mathrm{b}}=E_{\mathrm{x}-\mathrm{s}\mathrm{b}}{\epsilon }_{\mathrm{x}-\mathrm{s}\mathrm{b}}+K_0\gamma z$$

(3)

where E_x-sb is the elastic modulus of soilbag, and K₀ is the static earth pressure coefficient of the backfill.

By solving simultaneous Equations (1-3), the non-limited active earth pressure exerted on soilbags reinforced retaining walls subjected to surcharge loads can be determined.

3.2. Determination of K_q and K_γ

In limited-state scenarios, retaining walls exhibit complete sliding wedges, where the backfill’s internal friction angle (φ_cirt) and the friction angle at the backfill-wall interface (δ_crit) peak. In contrast, non-limited states yield only partial wedge development, resulting in lower, or intermediate, friction angles (φ_int for backfill internal friction and δ_int for backfill-wall friction). These intermediate values, lower than φ_cirt and δ_crit, indicate a less intense stress condition. Fan et al. [35] outline a procedure to calculate φ_int and δ_int,, suggesting these angles lie between their initial (φ₀ and δ₀, representing the starting internal and interface friction angles) and critical states. This approach helps in precisely gauging the frictional behavior of retaining walls under varying conditions.

The calculation equations provided by Fan et al. for φ_int and δ_int are as follows:

$$\begin{array}{c}{\tan \phi }_{\mathrm{i}\mathrm{n}\mathrm{t}}={\tan \phi }_0+\left({\tan \phi }_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-{\tan \phi }_0\right)\frac{s}{s_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}}(0\le s<s_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}})\\ {\tan \phi }_{\mathrm{i}\mathrm{n}\mathrm{t}}={\tan \phi }_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}(s\ge s_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}})\end{array}$$

(4)

$$\begin{array}{c}{\tan \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}={\tan \delta }_0+\left({\tan \delta }_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-{\tan \delta }_0\right)\frac{s}{s_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}}(0\le s<s_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}})\\ {\tan \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}={\tan \delta }_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}(s\ge s_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}})\end{array}$$

(5)

where s_crit is the critical lateral displacement necessary to reach the fully active state of the backfill. For walls with layered backfill, δ₀ is proposed to be half of φ_cirt (i.e., δ₀ = φ_cirt/2), and φ₀ can be derived from Equation (6):

$${(\frac{1}{{\cos \phi }_0}+\sqrt{{{\tan }^2\phi }_0+{\tan \phi }_0{\tan \phi }_0})}^2=\frac{1}{1-{\sin \phi }_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}}$$

(6)

Sherif et al. [36] provided an empirical equation for s_crit in retaining walls with varying densities of sand through model tests:

$$s_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}=H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}(7.0-0.13{\sin \phi }_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}})\left({10}^{-4}\right)$$

(7)

By solving Equations (4-7), the values of φ_int and δ_int can be determined.

From the equilibrium analysis of horizontal and vertical forces, as well as the moment equilibrium about the midpoint on the sliding surface of a differential flat element, a set of equations can be derived:

Equation (8) represents the balance of horizontal forces:

$$p_{\mathrm{x}-\mathrm{b}\mathrm{f}}{\cos \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}dz-R_{\mathrm{b}\mathrm{f}}\cos (90-\theta +{\phi }_{\mathrm{int}})dz/\sin \theta =0$$

(8)

Equation (9), addressing the vertical force equilibrium, are given by:

$$\begin{array}{c}{q}_{\mathrm{b}\mathrm{f}}\left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z\right)\cot \theta -\left(q_{\mathrm{b}\mathrm{f}}+{dq}_{\mathrm{b}\mathrm{f}}\right)\left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z-dz\right)\cot \theta +d{W}_{\mathrm{b}\mathrm{f}}-p_{\mathrm{x}-\mathrm{b}\mathrm{f}}{\sin \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}dz\\ -R_{\mathrm{b}\mathrm{f}}\mathrm{c}\mathrm{o}\mathrm{s}(90-\theta +{\phi }_{\mathrm{i}\mathrm{n}\mathrm{t}})dz/\mathrm{s}\mathrm{i}\mathrm{n}\theta =0\end{array}$$

(9)

Equation (10) pertains to the moment equilibrium:

$$\begin{array}{c}{p}_{\mathrm{x}-\mathrm{b}\mathrm{f}}\left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z-dz/2\right)\cot \theta {\sin \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}dz-q_{\mathrm{b}\mathrm{f}}\left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z\right)\cot \theta \left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z-dz\right)\cot \theta /2\\ +\left(q_{\mathrm{b}\mathrm{f}}+{dq}_{\mathrm{b}\mathrm{f}}\right)\left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z-dz\right)\left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z\right)\cot \theta /2\\ -\gamma \left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z-dz\right)\cot \theta \left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z\right)\cot \theta dz/2=0\end{array}$$

(10)

In those equations, θ is the sliding angle (θ=45°+ φ/2), q_bf is the force on the top of the element, R_bf is the normal reaction of the backfillsoil at rest, and dW is the weight of the element (dW_bf =γ(H_crit-z)cotθdz).

Equations (8) and (10) can be simplified as:

$$R_{\mathrm{b}\mathrm{f}}=p_{\mathrm{x}-\mathrm{b}\mathrm{f}}\frac{{\cos \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}\mathrm{s}\mathrm{i}\mathrm{n}\theta }{\mathrm{s}\mathrm{i}\mathrm{n}(\theta -{\phi }_{\mathrm{i}\mathrm{n}\mathrm{t}})}$$

(11)

$${dq}_{bf}=\gamma dz-2p_{x-bf}\tan \theta {\mathit{sin}\delta }_{int}dz/\left(H_{crit}-z\right)=0$$

(12)

Ignoring the second-order differential terms yields, Equation (9) can be expressed as:

$$\begin{array}{ll}{dq}_{\mathrm{b}\mathrm{f}}=\gamma dz+& q_{\mathrm{b}\mathrm{f}}dz/\left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z\right)\\ & -p_{\mathrm{x}-\mathrm{b}\mathrm{f}}\tan \theta {\mathit{sin}\delta }_{int}[1+{\cot \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}\cot \left(\theta -{\phi }_{\mathrm{i}\mathrm{n}\mathrm{t}}\right)]dz/\left(H_{\mathrm{c}\mathrm{r}\mathrm{i}\mathrm{t}}-z\right)\end{array}$$

(13)

Substituting Equation (12) into Equation (13) yields

$$q_{\mathrm{b}\mathrm{f}}=\tan \theta {\mathit{sin}\delta }_{int}\left({\cot \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}\cot \left(\theta -{\phi }_{\mathrm{i}\mathrm{n}\mathrm{t}}\right)-1\right)p_{\mathrm{x}-\mathrm{b}\mathrm{f}}$$

(14)

Substituting Equation (14) into Equation (12) yields:

$$\frac{{dq}_{bf}}{dz}=-\frac{2N{q}_{\mathrm{b}\mathrm{f}}}{\left(H_{crit}-z\right)}+\gamma$$

(15)

where N = 1/ (cotδ_intcot(θ-φ_int)-1).

By differentiation, the general solution of Equation (15) is

$$q_{\mathrm{b}\mathrm{f}}=C{\left(H_{crit}-z\right)}^{2N}+\frac{\gamma }{2N-1}\left(H_{crit}-z\right)$$

(16)

where C is a constant, which can be determined by the boundary condition (z = 0, q_bf = q,). C is determined as

$$C=\left(q-\frac{\gamma H_{crit}}{2N-1}\right){H_{crit}}^{-2N}$$

(17)

Thus, Equation (16) can also be expressed as:

$$q_{\mathrm{b}\mathrm{f}}=\left(q-\frac{\gamma H_{crit}}{2N-1}\right){\left(\frac{H_{crit}-z}{H_{crit}}\right)}^{-2N}+\frac{\gamma }{2N-1}\left(H_{crit}-z\right)$$

(18)

Substituting Equation (18) into Equation (14), p_x-bf can be obtained:

$$\begin{array}{ll}{p}_{\mathrm{x}-\mathrm{b}\mathrm{f}}=\frac{N}{\tan \theta {\mathit{sin}\delta }_{int}}& ({\frac{H_{crit}-z}{H_{crit}})}^{2N}q\\ & -\frac{N}{\left(2\mathrm{N}-1\right)\tan \theta {\mathit{sin}\delta }_{int}}\frac{H_{crit}-z}{z}{[1-(\frac{H_{crit}-z}{H_{crit}})}^{2N-1}]\gamma z\end{array}$$

(19)

Solving for Equations (19) and (2) yields:

$$K_{\mathrm{q}}=\frac{N}{\tan \theta {\mathit{sin}\delta }_{int}\left({\cot \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}\cot \left(\theta -{\phi }_{\mathrm{i}\mathrm{n}\mathrm{t}}\right)-1\right)}({\frac{H_{crit}-z}{H_{crit}})}^{2\left({\cot \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}\cot \left(\theta -{\phi }_{\mathrm{i}\mathrm{n}\mathrm{t}}\right)-1\right)}$$

(20)

$$K_{\mathsf{\gamma }}=\frac{H_{crit}-z}{z\left(3-{\cot \delta }_{\mathrm{int}}\cot \left(\theta -{\phi }_{\mathrm{int}}\right)\right)\tan \theta {\mathit{sin}\delta }_{int}}{[1-\left(\frac{H_{crit}-z}{H_{crit}}\right)}^{2{\cot \delta }_{\mathrm{i}\mathrm{n}\mathrm{t}}\cot \left(\theta -{\phi }_{\mathrm{i}\mathrm{n}\mathrm{t}}\right)-3}]$$

(21)

3.3. Determination of E_x-sb

The lateral compressive stress-strain characteristics of soilbags were investigated using a direct shear test apparatus, details of which are elaborated in the work of Fan et al. [5]. In this setup, two layers of soilbags were positioned on the apparatus’s base. The lower layer was anchored using two plates fixed to the base, ensuring stability during testing. Above the soilbags, a loading plate was placed, accompanied by a right-side plate. This arrangement was designed to apply normal stresses to the soilbags. A critical component of the setup was a displacement sensor attached to the right-side plate. This sensor’s role was to precisely track any lateral displacement that occurred during the testing process. For the application of lateral force, a controlled speed of 1 mm/min was maintained, and a load cell was incorporated to measure the lateral force exerted.

The results of these tests are depicted in Figure 4, which illustrates the lateral compressive stress-strain relationship of the soilbags under varying normal stresses. The data indicates a roughly linear relationship between lateral stress and strain across different levels of normal stress. This linear trend suggests a consistent response of the soilbags to lateral compression. Figure 5 presents the correlation between the lateral compressive modulus (E_x-sb) and the normal stress (σ_z-sb) of the soilbag. This relationship reveals that the E_x-sb value progressively increases as σ_z-sb increases, demonstrating a distinct exponential connection between these two parameters. This observed relationship can be quantitatively described by an exponential function, where the lateral stress-strain relationship of the soilbag is modeled as:

$${\sigma }_{\mathrm{x}-\mathrm{s}\mathrm{b}}=201{\mathrm{e}}^{0.15{\sigma }_{\mathrm{z}-\mathrm{s}\mathrm{b}}}{\epsilon }_{\mathrm{x}-\mathrm{s}\mathrm{b}}$$

(22)

Further simplification leads to the expression for the lateral compressive modulus (E_x-sb) as:

$$E_{\mathrm{x}-\mathrm{s}\mathrm{b}}=201{\mathrm{e}}^{0.15{\sigma }_{\mathrm{z}-\mathrm{s}\mathrm{b}}}$$

(23)

3.4. Validation

In the described model test, key parameters include the unit weight of the backfill (γ) at 17.6 kN/m³, the internal friction angle (φ_crit) at 35.4°, the frictional angle between the back of the soilbag and the backfill (δ_crit) at 28.1°, the width of the soilbag reinforced retaining wall (B_sb) at 0.4 m, and the critical height of the wall (H_crit) at 0.95 m. Utilizing the previously mentioned equations, the lateral displacement and lateral earth pressures of the backfill along the wall height were calculated under varying surcharge loads of 2 kPa, 4 kPa, 6 kPa, and 8 kPa.

The findings revealed that under a 2 kPa surcharge load, the calculated lateral displacement along the wall height remained below the critical lateral displacement (s_crit). This suggests that, in this scenario, the wall did not reach its limit state, as depicted in Figure 6a. However, as the surcharge load increased to 4 kPa, a notable change was observed. The lateral displacement of the soilbag at the upper part of the wall exceeded s_crit, indicating that this portion of the wall had reached the limit state, as shown in Figure 6b. Furthermore, the extent of the wall reaching the limit state increased progressively with the increment of the surcharge loads, as shown in Figure 6c,d.

Another key observation from Figure 6a to d is the close agreement between the calculated earth pressure and the experimental data. This alignment validates the accuracy of the calculations and the effectiveness of the equations used in predicting the behavior of the soilbag and the retaining wall under different loading conditions.

4. Conclusions

Model tests on the soilbag reinforced retaining wall subjected to surcharge loads under the non-limited state were carried out to analyze the non-limited earth pressures. The findings yielded several key insights:

• The lateral earth pressure distribution on the soilbag reinforced retaining wall subjected to surcharge loads under the non-limited state is non-linear. As the surcharge loads increased, there was a corresponding increase in the earth pressure.
• It was determined through experimentation that the lateral stress-strain behavior of the soilbags adheres to a linear model. Furthermore, an exponential relationship between the lateral compression modulus of the soilbags and the normal stress was established.
• The analytical solution, developed based on the force equilibrium of differential elements, proved to be highly effective in predicting the non-limited lateral earth pressure on the soilbag reinforced retaining wall under surcharge loads. The close agreement of this analytical solution with the experimental data highlights its potential as a reliable and accurate tool for predicting lateral earth pressure in similar scenarios.