
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39266606
72094
10.1038/s41598-024-72094-8
Article
Passive earth pressure on vertical rigid walls with negative wall friction coupling statically admissible stress field and soft computing
Shiau Jim 1
Nguyen Tan nguyentan@tdtu.edu.vn

23
Bui-Ngoc Tram 45
1 https://ror.org/04sjbnx57 grid.1048.d 0000 0004 0473 0844 School of Engineering, University of Southern Queensland, Toowoomba, QLD 4350 Australia
2 https://ror.org/01drq0835 grid.444812.f 0000 0004 5936 4802 Smart Computing in Civil Engineering Research Group, Ton Duc Thang University, Ho Chi Minh City, Vietnam
3 https://ror.org/01drq0835 grid.444812.f 0000 0004 5936 4802 Faculty of Civil Engineering, Ton Duc Thang University, Ho Chi Minh City, Vietnam
4 https://ror.org/02ryrf141 grid.444823.d 0000 0004 9337 4676 Laboratory for Computational Mechanics, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, Vietnam
5 https://ror.org/02ryrf141 grid.444823.d 0000 0004 9337 4676 Faculty of Mechanical - Electrical and Computer Engineering, School of Technology, Van Lang University, Ho Chi Minh City, Vietnam
12 9 2024
12 9 2024
2024
14 2132214 6 2024
3 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
It is well known that the roughness of a wall plays a crucial role in determining the passive earth pressure that is exerted on a rigid wall. While the effects of positive wall roughness have been extensively studied in the past few decades, the study of passive earth pressure with negative wall friction is rarely found in the literature. This study aims to provide a precise solution for negative friction walls under passive wall conditions. The research is initiated by adopting a radial stress field for the cohesionless backfill and employs the concept of stress self-similarity. The problem is then formulated in a way that a statically admissible stress field be developed throughout an analyzed domain using a two-step numerical framework. The framework involves the successful execution of numerical integration, which leads to the exploration of the statically admissible stress field in cohesionless backfills under negative wall friction. This, in turn, helps to shed light on the mechanism of load transfer in such situations so that reliable design charts and tables be provided for practical uses. The study continues with a soft computing model that leads to more robust and effective designs for earth-retaining structures under various negative wall frictions and sloping backfills.

Keywords

Negative wall friction
Admissible stress field
Passive earth pressure
Soft computing
Subject terms

Scientific data
Statistics
Civil engineering
http://dx.doi.org/10.13039/100007224 National Foundation for Science and Technology Development NCUD.02-2022.08 Nguyen Tan issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The presence of wall roughness can have a significant impact on the mechanism of load transfer in cohesionless backfill as well as the magnitude of lateral earth pressures. It is important to note that there are two distinct types of wall roughness. Depending on the relative movement between the soil in the backfill and the wall, positive wall friction occurs when the soil moves upward faster than the wall, while negative wall friction occurs when the wall moves upward faster than the soil1–4.

In active scenarios, friction is predominantly positive, serving as a stabilizing factor for retaining walls and structures4. However, certain conditions, such as excavation for repair work, can lead to situations where negative wall friction occurs, necessitating additional support or bracing against the backfill soil to ensure structural stability. In passive conditions, both positive and negative wall friction can influence the overall stability of the structure3. Understanding the influence of negative wall friction is particularly crucial as it significantly impacts the lateral earth pressure exerted on retaining structures. This phenomenon is especially important in the design of retaining walls and anchors, as it affects the overall load distribution and the potential for failure. For example, in conditions where the wall or anchor moves upward faster than the surrounding soil, the resultant increase in earth pressure may necessitate a more robust design to prevent structural failure. Additionally, the negative wall friction angle is a critical factor in determining the uplift capacity of anchors5, as it influences the relative movement between the anchor or wall and the surrounding soil. Typically analyzed using a curved rupture surface and the limit equilibrium method6, the direction of wall friction (positive or negative) must be determined based on anticipated ground movements under specific field conditions, which is essential for optimizing the design and ensuring the safety of geotechnical structures.

The introduction of wall friction transforms the rupture surface from a simple, planar configuration to a curved one, an observed phenomenon documented by Savage et al.1. Notably, the presence of negative wall–soil friction not only impacts the magnitude of the lateral earth thrust exerted on the wall but also results in substantial deviations in the failure surface, distinct from scenarios involving positive wall–soil friction3. Therefore, gaining a comprehensive understanding of the load transfer mechanisms associated with negative wall–soil friction holds exceptional significance within the realm of geotechnical engineering, as it can carry profound implications for the stability and design of retaining walls and structures.

In contrast to the well-established research with positive wall–soil friction7,8, published work on negative wall–soil friction is scarce. To date, limited attention has been given to the study of passive earth pressure applied on rigid retaining walls under negative wall–soil friction. One of the few works in relation to the effect of negative soil-wall friction is provided by Kérisel and Absi9, who presented a comprehensive table of passive earth pressure with negative soil-wall friction using the slip-line method. Kumar and Rao2 later conducted a study on the passive earth pressure acting on a wall by varying the wall–soil friction from negative to positive, using the limit equilibrium method by adopting Terzaghi’s failure surface10. Choudhury and Rao11 extended this to the seismic context, emphasizing the strong influence of seismic accelerations on passive resistance, especially with negative wall friction. On the other hand, the well-known Coulomb’s solution proposed by Poncelet12 incorporates both negative and positive wall–soil frictions3,13. It’s important to recognize that these methods focus on satisfying the kinematic condition, often at the expense of the equilibrium condition. Consequently, the passive load calculated by these traditional approaches generally exceeds the actual collapse load, as noted by James and Bransby14.

Very recently, Nguyen15 presented precise solutions for passive earth pressures under positive soil-wall friction using a fully plastic solution. The method is also referred to as the Incipient Failure Everywhere (IFE) solution or the rigid-plastic solution16,17. It has a long-established history in geotechnical engineering, and it adopts the assumption of critical limit equilibrium and the Mohr–Coulomb failure criterion18. The fully plastic solution was comprehensively investigated by Booker19 and has its roots in Rankine’s seminal paper20 on the study of slopes with traction-free surfaces and arbitrary inclination. Later, Rankine’s theory was combined with the method of characteristics by Sokolovskii21 for the infinite wedge theory. Nguyen’s solutions15, rooted in a statically admissible stress field, are highly valued for potentially providing an exact solution, as highlighted by Lancellotta22. These solutions incorporate the strengths of Sokolovski’s approach21, including mathematical elegance, a lower bound solution, and effectiveness in handling large soil deformations near walls, as observed by James and Bransby14.

The fully plastic solution, detailed by Nguyen et al.4, effectively addresses negative wall–soil friction in active earth pressure. Leveraging advanced soft computing techniques, this study significantly advances the estimation of active earth pressure in cohesionless backfill. Employing a Bayesian regularization backpropagation neural network, the model demonstrates remarkable alignment with measured values, presenting itself as a versatile tool for retaining wall design, especially in the presence of negative wall–soil friction, thereby pushing the boundaries of geotechnical engineering.

In the realm of passive earth pressure, drawing inspiration from the work of Nguyen and Shiau3, the study investigates the impact of negative wall–soil friction on various rigid wall geometries. Utilizing statically admissible stress fields, conservative solutions for the passive earth pressure coefficient are derived, providing insights into load transfer mechanisms within cohesionless backfill. The study identifies the stress state transition zone and elucidates the influence of wall inclination on passive earth pressure, contributing to the creation of comprehensive design charts for engineering practice. Consequently, a research gap is identified concerning the determination of statically admissible stress fields in cohesionless backfill with a vertical wall’s rear face and negative wall–soil friction under passive earth pressure conditions, posing a challenge in optimizing retaining wall design in engineering practice.

To bridge this research gap, this study explores the statically admissible stress field in the backfill of a rigid retaining wall with a vertical wall’s rear face under negative wall–soil friction, employing the fully plastic solution. Machine learning is integrated to gain deeper insights into the load transfer mechanism, proposing a user-friendly tool for practical retaining wall structure design. The primary objectives of this study are: (1) to determine the accurate passive earth pressure exerted on a rigid wall in the presence of negative soil-wall friction, (2) to investigate the statically admissible stress field using the fully plastic solution, and (3) to utilize machine learning techniques for additional insights into the physical behavior of the backfill material under diverse loading conditions. The outcomes of this research promise to contribute significantly to various fields within granular mechanics, deepening understanding of retaining wall structures, uplift anchor structures, and silo structures.

The theory

Equilibrium and yield criteria

The problem with negative wall–soil friction is presented with a potential curved sliding surface that is concave. Savage et al.1 experimentally demonstrated the existence of this failure surface (ref. Fig. 1a). Note that the stress components in the backfill must satisfy the equilibrium conditions in the polar coordinate system as presented in Eqs. (1) and (2) (ref. Fig. 1b).1 ∂σrr∂r+1r∂τrθ∂θ+σrr-σθθr=γcosθ

2 ∂τrθ∂r+1r∂σθθ∂θ+2τrθr=-γsinθ

where r is equal to z/cos(θ) (ref. Fig. 1b).3 σr-σθ2+4τrθ2=σr+σθ2sin2ϕ

Fig. 1 (a) Passive earth pressure with negative wall–soil friction; (b) Geometry and stress components along various planes of a rigid retaining wall subjected to negative wall friction.

Equation 3 shows a Coulomb-type yield criterion. Using this equation as well as the in-plane mean stress (p) and the relative angle ψ, the polar stress components are shown in Eq. (4).4 σrr=p(1+sinϕcos2ψ)σθθ=p(1-sinϕcos2ψ)τrθ=psinϕsin2ψ

where p = (σrr + σθθ)/2, i.e., p is a mean normal stress. Using the hypothesis of self-similarity proposed by Sokolovskii21, Eqs. (5) and (6) are presented.5 p=γrχθ

6 ψ=ψθ

where χ(θ) is a dimensionless mean normal stress function with the key parameter (θ) in the given radial stress field.

Substituting the polar stress components in Eq. (4) and the “self-similarity” equations in Eqs. (5) and (6) into Eqs. (1), (2), a new set of the form of equations for stress equilibrium is generated. Note that the two variables (χ) and (ψ) are involved in the equations as shown in Eqs. (7), (8).7 sinϕsin2ψdχdθ+2χsinϕcos2ψdψdθ+1+χ1+sinϕcos2ψ=cosθ

8 1-sinϕcos2ψdχdθ+2χsinϕsin2ψdψdθ+1+χsinϕsin2ψ=-sinθ

Further, a pair of Ordinary Differential Equations (ODEs) is presented in Eqs. (9) and (10) based on Eqs. (7),(8).9 dχdθ=χθsin(2ψ)-sin(2ψ+θ)cos2ψ-sinϕ

10 dΨdθ=dψdθ+1=sinϕcos(2ψ+θ)+χ(θ)cos2ϕ-cosθ2χθsinϕsinϕ-cos2ψ

For the processing of numerical solutions, it is necessary to transform the polar coordinate stress system into a Cartesian coordinate system. This is shown in Eqs. (11)–(15).11 x=rsinθ

12 z=rcosθ

13 σzz=σrr+σθθ2+σrr-σθθ2cos2θ-τrθsin2θ=p+qcos2Ψ

14 σxx=σrr+σθθ2-σrr-σθθ2cos2θ+τrθsin2θ=p-qcos2Ψ

15 τxz=σrr-σθθ2sin2θ+τrθcos2θ=qsin2Ψ

where16 Ψ=ψ+θ

Singular line in a radial stress field

A singular line is a characteristic line or a line of stress discontinuity representing wall–soil roughness and the respective geometry conditions3. The passive stress state of positive wall friction causes the major principal stress to rotate in a clockwise direction15. Conversely, in the case of negative wall friction, the direction of the major principal stress rotation is counterclockwise. As a result, the singular line in the backfill of negative wall friction differs from that of positive wall friction, and it exhibits distinct characteristics in comparison to backfill with positive wall friction.

The area between the singular line and the traction-free surface is referred to as Rankine’s region. Depending on the degree of negative wall friction, this region may either vanish or take up the entire backfill region if the traction-free surface coincides with the singular line (see Fig. 2a–c). In the Rankine region, the relative angle (ψ) on the traction-free slope can be determined by solving Eq. (10) with the condition that Ψ remains constant and the slope is traction-free (χ = 0)23. To prevent soil particles from sliding down the traction-free slope, the traction-free surface must maintain an angle that is less than the angle of repose for uncompacted cohesionless sand.Fig. 2 (a) A possible line of stress characteristics for a smooth wall; (b) The Rankine region encompassing the entire backfill for a negatively rough wall; (c) A possible line of stress discontinuity for a perfectly negative rough wall; (d) A Mohr–Coulomb diagram illustrating the stress relationship between two adjacent regions. It is important to note that within the Rankine regions, the direction of the major principal stress remains constant, and the stress components follow a strictly linear pattern. Conversely, outside the Rankine regions, both the direction of the major principal stress and the stress components exhibit nonlinearity.

Specifically, when the angle between the slope and the rear face of the wall is an acute angle, the backfill is always bisected by a line of stress discontinuity. On the other hand, if the angle between the slope and the rear face of the wall is obtuse, the backfill will be separated by a stress characteristic line for low levels of negative wall friction. However, for high levels of negative wall roughness, the line of stress characteristic will be replaced by a stress discontinuity line.

The inclusion of a stress discontinuity in the backfill is imperative as it would enable the establishment of a stress field that satisfies all relevant boundary conditions. The study of stress discontinuity in radial stress fields has been explored in numerous works, including Lee24, Shield25, Booker19, Savage18; thus, it will not be repeated here. To ensure stability and balance, the continuity of both normal and tangential stress components must be insisted across the stress discontinuity line. The stress discontinuity line, as a result of self-similarity, is a straight line that separates regions C1 and C2 (see Fig. 2d). The stress jump along the stress discontinuity line, represented by an angular coordinate θD, can be expressed mathematically by Eqs. (17),(18).17 tanΨ1-θDtanΨ2-θD=1-sinϕ1-sinϕ1+sinϕ1+sinϕ

18 sin2Ψ1-θDsin2Ψ1-θDsin2Ψ2-θDsin2Ψ2-θD=p2p2p1p1

Interestingly, it is to be noted that the stress discontinuity vanishes when the stress discontinuity line coincides with a line of stress characteristics.

Numerical procedures

Equations (9) and (10), which represent the governing equations, involve two primary variables: χ(θ) and ψ(θ). By applying the specified boundary conditions, a numerical solution was obtained for these ordinary differential equations using the initial-boundary value problem approach, where the boundary values for the relative angle (ψ) and the mean stress (χ) are specified at the same angular coordinate (θ). Alternatively, these equations can also be solved using the boundary-value problem method, where (χ) and (ψ) a are defined at different locations.

The boundary conditions for solving boundary-value problems are shown in Table 1. Note that either theoretical postulates or numerical methods can be utilized to convert them into initial-value problems.Table 1 Initial boundary conditions.

Planes	ψ	χ	
Traction-free surface	Nil	Zero	
On the rear face of wall	ψw	Nil	
Equation (29) is used to compute ψw.

Solutions to the governing equations (Eqs. (9) and (10)) were computed using two distinct approaches. For initial-value problems, a 5th order adaptive Runge–Kutta method was chosen to address instability issues in non-smooth domains. In contrast, boundary-value problems were tackled using a combination of the shooting method and the 5th order adaptive Runge–Kutta method, which transformed them into initial-value problems with an interpolation technique. The effectiveness of these numerical approaches in achieving precise solutions is a significant focus of this study.

Runge–Kutta method for an initial-value problem

Mathematically, a typical initial value problem can be stated as,19 y′=f(x,y),y(x=θo)=α

where y is a vector, and it consists of two variables (χ and ψ) that are independent of each other as shown in Eqs. (9) and (10), respectively. The global variable (x) is angular coordinate (θo) and α has two known values (α1 and α2) for the vector y at the same location, χθ=θ0=α1,ψθ=θ0=α2. The adaptive step size procedure will be adjusted based on the estimated error (Eq. 20) when the integration process reaches a singular line,20 h0=h1Δ0Δ1

where Δ0 is the chosen accuracy, and Δ1 is the error yielded by the step size h1.

A truncation error Δ(h) is estimated as below for the mth order formula,21 Δ(h)=ym+1(x+h)-ym∗(x+h)=∑i=16ci-ci∗ki

Note that the ym∗ and ym+1 are the 5th and 4th orders embedded formulas presented in Fehlberg 26. This is defined as follows,

The 5th order formula:22 ym+1(x+h)=ym(x)+∑i=16ciki+Oh6

The 4th order formula:23 ym+1∗(x+h)=ym(x)+∑i=16ci∗ki+Oh5

24 k1=hfx,y

25 ki=hfx+aih,y+∑j=0i-1bijkj,i=2,3,...,6

Noting that the set of coefficients ai, bi, ci, and c*i proposed by Cash and Karp27 provide better accuracy than Fehlberg’s original values28, hence, they are adopted in this study (ref. Table 2).Table 2 Cash-Karp coefficients for embedded Runge–Kutta method.

i	ai	bij	ci	ci∗	
1	1	–	–	–	–	–	37378	282527648	
2	15	15	–	–	–	–	0	0	
3	310	340	940	–	–	–	250621	1857548384	
4	35	310	-910	65	–	–	125594	1352555296	
5	1	-1154	52	-7027	3527	–	0	27714336	
6	78	163155296	175512	57513824	44275110592	2534096	5121771	14	

The method employed in this study involves adjusting the step size in the integration process based on the truncation error (Δ1) in comparison to a pre-defined error threshold (Δ0). If the truncation error is less than the threshold (Δ1 < Δ0), the step size will be reduced, but if it is greater than or equal to the threshold (Δ1 ≥ Δ0), the step size will be increased. This adaptive step size approach is shown in Fig. 3. The integration process requires setting both an initial step size and a numerical tolerance for each iteration. The numerical tolerance is consistently maintained at eTol = 1E−6 across all runs, The initial step size is initially set to an extremely small value, such as h = 1E−20. However, this value may be dynamically adjusted within a range, potentially as small as h = 1E−50, in response to numerical perturbations that arise when addressing singularities during the computation.Fig. 3 Demonstration of the shooting method for the case (λ = 0, ϕ = δ = 40°, β = − 30°).

Shooting method for the two-point boundary value problem

The initial value solution for the domain bounded by a characteristics line and the rigid wall may not be an accurate representation of the true one. This is because the initial-value solution fails to extend from Rankine’s region (bounded by the traction-free surface and characteristic line) to the non-Rankine region (bounded by the rear face of the wall and characteristic line). In other words, the solution does not necessarily link the numerical results (χ and ψ) between Rankine’s and non-Rankine’s zones, as reported by Nguyen and Pipatpongsa29. It follows that the problem can be defined as a two-point boundary value problem, as discussed in the works of Kiusalaas30.26 χ′=χ′χ,ψ,θ=dχdθψ′=ψ′(χ,ψ,θ)=dψdθ=dΨdθ-1,χθ=θ1=uψθ=θ1=ψ1,χθ=θ2=χ2ψθ=θ2=v

The numerical solution to the set of ordinary differential equations (Eqs. (9) and (10)) involves finding the two values u and v by integrating the equations between the initial point (θ1) and the final point (θ2). This is known as the two-point boundary value problem and is graphically shown in Fig. 3. The shooting method is adopted to obtain accurate numerical results for this problem.

A root-finding method based on Brent’s algorithm21 was adopted to solve for the unknown value of “u”. It is important to note that the solution of “χ” at θ = θ2 is dependent on the unknown value of “u” at θ = θ1, meaning that “u” is the root of Eq. (27):27 r(u)=f(u)-χθ=θ2=f(u)-χ2=0

herein, r(u) is the boundary residual between the specified boundary-value χθ=θ2=χ2 and, the computed boundary-value for χθ=θ1.

The numerical procedure, which involves the use of Brent’s algorithm, is executed in three consecutive steps. In the first step, a range of constraints for the root “u” of Eq. (27) is specified. It is important to note that the results for each problem can vary greatly, so the range of constraints for “u” was carefully determined for each problem. In the second step, Brent’s algorithm is executed iteratively to find “u” from Eq. (27), and the evaluation of f(u) is required for each iteration by solving an initial-value problem using the 5th order adaptive Runge–Kutta method. Finally, the numerical procedure is completed in the third step by solving the differential equations (Eqs. (9) and (10)) once more as an initial-value problem, until the root (u) is found.

Rankine’s and non-Rankine’s regions

A crucial part of the analysis involves identifying the Rankine region within the backfill. This identification requires considering various factors15, including the geometry of the backfill and the level of wall friction. To determine the direction of principal stress at the rear face of the wall, several assumptions related to wall roughness are made, which are outlined below.

For negative soil-wall friction1, it is assumed that the following condition holds at the rear face of the wall (θ = 0).28 τrθ=-σθtanδ

By substituting Eq. (4) into Eq. (28), we obtain the boundary condition (Eq. 29) for the “relative angle” of the major principal stress (ψ) at the wall rear face (θ = 0),29 ψw=12π+δ+arcsinsinδsinϕ

To begin an analysis, the first step is to specify the level of wall roughness for both the wall and backfill geometry (see Fig. 4). Consequently, Rankine regions for the backfill can be determined. When the whole backfill is within the Rankine region, the direction of the major principal stress is constant throughout the backfill (referred to as case 1). However, for higher degrees of wall roughness, a stress discontinuity line divides the backfill into Rankine’s and non-Rankine’s regions (case 2). Conversely, for lower degrees of wall roughness, a stress characteristic line divides the backfills (case 3). For a positive sloping backfill stabilized at the soil’s internal friction angle, the stress characteristic line intersects the traction-free surface, leading to the disappearance of the Rankine region from the backfill (case 4).Fig. 4 The effect of negative wall friction on the mechanism of load transfer.

Problems with a stress discontinuity line

The direction of the major principal stress at the traction-free surface is uncertain, as previously indicated in Table 1. Thus, the numerical procedure begins with solving Eq. (30) to obtain the value of ψ at the traction-free surface, as suggested by Booker19.30 cosθ-sinϕcos2ψθ+θ=0

Moreover, the direction of the major principal stress at the rear face of the wall is influenced by the wall–soil friction angle, causing it to tilt downward from the horizontal direction. Hence, Eq. (29) is employed to determine the relative angle of the major principal stress (ψw) with respect to the wall–soil friction angle (δ). The flowchart showing the numerical procedure for problems with a stress discontinuity line is presented in Fig. 5.Fig. 5 Problems with a stress discontinuity line ~ a flow chart showing the numerical procedure.

To handle the singularities at the rear face of the wall and the traction-free surface, a numerical perturbation technique is used. At the traction-free surface, the initial values for χ and ψ are slightly perturbed by the amount of 10–28 and 10–12, respectively. On the other hand, the integration process for this two-dimensional problem can reach the rear face of the wall without causing numerical stability issues.

Problems with a line of characteristics

The problem with a line of characteristics can be tackled by using the numerical procedures for solving the two-point boundary value problem under various soil backfill geometries and wall roughness.

For cases where the inclination angle of the traction-free surface is less than the internal friction angle, the procedure involves solving the initial value boundary problem in the region bounded by the characteristic line and the traction-free surface, and then shooting the pair of governing equations from the characteristic line to the wall rear face.

For cases where the soil internal friction angle equals the inclination angle of the traction-free surface, the numerical procedure is stepwise, involving an iterative process to find the value of the mean stress (χ) at the traction-free surface, and then extending the solution to the traction-free surface using 6th order polynomial fitting. The details of these procedures are explained in a flowchart shown in Fig. 6.Fig. 6 A flowchart showing the numerical procedure for problems with a characteristic line.

Results and discussion

Vertical wall with horizontal backfill

This study compares earth pressure coefficients Kp acting on rigid walls with negative wall friction for a vertical wall (λ = 0) with horizontal backfill (β = 0). The obtained theoretical results are compared with those published ones using the limit equilibrium method and the slip line method. Figure 7 shows that the results obtained in this study are in good agreement with those of Kumar and Rao2, whereas the slip line method proposed by Kérisel and Absi9 overestimated the passive earth pressure coefficient for rough walls with approximately δ/φ > 0.7.Fig. 7 Comparison of passive earth pressure coefficient with various wall friction angles between the present solution and those of Kerisel and Absi9 using slip-line method and Kumar and Rao2 using the limit equilibrium method.

The numerical results of mean stress (χ) and relative angle of the major principal stress (ψ) for various degrees of negative wall friction are shown in Fig. 8. Note the presence of negative wall friction is to induce a stress discontinuity in the backfill. The location of this line of stress discontinuity is further away from the wall as the negative wall friction increases, accompanied by a pronounced mean stress jump.Fig. 8 Results of the relative angle of the major principal stress (ψ) and mean stress (χ) in case of (ϕ = 30°, λ = 0, and β = 0) with various degrees of negative wall friction.

To the best of knowledge, the research conducted by Savage et al.1 provides an in-depth analysis of failure mechanisms in cohesionless soil behind rigid walls subjected to negative soil-wall friction. This seminal study combines both analytical and experimental methods, with the analytical results showing excellent alignment with experimental data, especially concerning the effects of negative wall–soil friction.

A critical aspect of this study is the use of a stress discontinuity surface to meet static equilibrium requirements, which is essential for addressing the boundary conditions associated with negative wall–soil friction. Savage et al.1 validated this approach by demonstrating a close match between theoretical predictions and experimental observations of failure surfaces. The analysis, as compared to the findings of Savage et al.1, is depicted in Fig. 9. This comparison shows a strong correlation, particularly in the location of the stress discontinuity surface, where the stress jump occurs. This key feature is crucial for establishing a statically admissible stress field and supports the validity of the proposed model.Fig. 9 Comparison between the present solution and Savage et al.1: (a) Relative angle of the major principal stress (ψ) and (b) Mean stress (χ), (c) the passive earth pressure coefficient (Kp) for the case of ϕ = 38°, λ = 0, and β = 0 with various degrees of negative wall friction.

Furthermore, the passive earth pressure coefficient (Kp) is compared with the results from Savage et al.1 (see Fig. 9c). The comparison reveals a significant level of agreement, further reinforcing the robustness of the theoretical model.

The stress distribution in the backfill is significantly affected by negative wall friction. Figure 10 shows that the vertical stress increases as the negative wall roughness increases, whereas the increase in wall friction reduces both the horizontal and shear stresses. Interestingly, the stress state changes from passive to active when the wall roughness reaches a perfectly rough state (δ = − φ). Also note that the wall friction only affects the stress distribution in the non-Rankine’s region bounded by the rear face of the wall and the line of stress discontinuity, whereas stress distribution in the Rankine region is independent of the magnitude of wall roughness.Fig. 10 Cartesian stress components at a height of interest in the backfill in case of (ϕ = 30°, λ = 0, and β = 0) with various degrees of negative wall friction: (a) normalized vertical stress σz/γz; (b) normalized horizontal stress σx/γz; (c) normalized shear stress σx/γz.

In the case of a perfectly smooth wall and horizontal backfill, the Rankine region occupies the entire backfill. However, for a rough wall, a line of stress discontinuity appears, dividing the backfill into two separate regions. The major principal stress direction is horizontal in the Rankine region, but in the negative wall backfill, the direction of the major principal stress tilts upward (see Fig. 11). This direction of the major principal stress aligns with the mechanism of negative wall friction, in which the wall moves faster than the soil.Fig. 11 Trajectories of the principal stresses in the backfill in case of (ϕ = 30°, λ = 0, and β = 0) with various degrees of negative wall friction (the stress discontinuity is shown in a red dashed line): (a) perfectly smooth wall (δ = 0); (b) rough wall (δ = − 15°); (c) perfectly rough wall (δ = − 30°).

Vertical wall with negative sloping backfills

It is worth noting that the earth pressure coefficient under a smooth wall is equivalent to the earth pressure coefficient at the central line of a symmetrical sand heap under a wedge shape. This has been extensively studied by Booker19 and Michalowski and Park31 using the fully plastic solution for the sand heap problem. Furthermore, the horizontal earth pressure at the center of the symmetrical sand heap is equivalent to the horizontal earth pressure acting on a vertical wall without wall roughness. Thus, their works are used to validate the obtained results in Fig. 12.Fig. 12 Comparison between the present method and those of Booker19 and Michalowski and Park31.

Figure 12 shows a comparison of the presented results with those Booker19 and Michalowski and Park31. Numerical results have shown that19’s numerical approximation of the IFE solution was not entirely convincing for small values of tan (θ°) (ref. Fig. 12a and b). By checking the weight balance condition, it is found that the error between the present method and the total weight of the symmetrical sand heap is only 0.00574%, indicating that the results presented by19 can be improved. The slight discrepancies between the results of the current study and those reported by Michalowski and Park31 could be attributed to the limited application of limit analysis for arch action.

The numerical results for the mean stress (χ) and relative angle of principal stress (ψ) are shown in Fig. 13. The degree of wall roughness affects the location of the stress discontinuity, with greater degrees of roughness causing it to move further from the wall and result in larger stress jumps.Fig. 13 Results of the mean stress (χ) and relative angle of principal stress (ψ) in case of ϕ = 30°, λ = 0, β = − 30° for various degrees negative wall friction.

The stress state in the vicinity of the wall in a negative sloping backfill is significantly influenced by the degree of negative wall friction. As shown in Fig. 14, an increase in negative wall friction results in an increase in the magnitude of the vertical stress at the wall while decreasing the horizontal stress at the wall. Notably, the stress state in the vicinity of the wall undergoes a transition from a passive stress state for smooth walls to an active stress state for perfectly rough walls. Furthermore, it is interesting to observe an abrupt change in the magnitude of the vertical stress as the value of negative wall friction equals the internal friction angle of soils. These findings shed light on the complex stress distribution and behavior of cohesionless materials subjected to negative wall friction, with significant improvement over current considerations in design practices.Fig. 14 Profiles of rectangular stress components at a height of interest in case of (ϕ = 30°, λ = 0, β = − 30°) for various degrees of negative wall friction (δ): (a) normalized vertical stress σz/γz; (b) normalized horizontal stress σx/γz; (c) normalized shear stress σx/γz.

The phenomenon of soil arching is evident in the negative sloping backfill, as shown in Fig. 15. The degree of wall roughness has a significant impact on the span of the arch. As the roughness of the wall increases, the arch expands. Moreover, the trajectories of the principal stresses highlight the change in the major principal stress direction from a horizontal orientation to a diagonally upward orientation. This observation suggests that the negative wall friction not only affects the magnitude of stresses in the vicinity of the wall, but also the direction of the stresses in the backfill. Such findings provide valuable insights for the design of retaining structures under negative wall friction.Fig. 15 The trajectories of principal stresses in negative sloping backfill in case of (ϕ = 30°, λ = 0, β = − 30°) with various degrees of negative wall friction (the stress discontinuity is shown in a red dashed line): (a) perfectly smooth wall (δ = 0); (b) rough wall (δ = − 15°); (c) perfectly rough wall (δ = − 30°).

Vertical wall with positive sloping backfills

In contrast to the vertical wall with horizontal and negative sloping backfills, the line of stress discontinuity disappears in the positive sloping backfill, but the line of stress characteristic divides the backfill into two regions: namely the Rankine and non-Rankine regions. See Fig. 16 for the numerical results of the mean stress (χ) and relative angle of principal stress (ψ). The magnitude of the negative wall friction causes the direction of the major principal stress to rotate counterclockwise while also reducing the magnitude of the mean stress at the wall’s rear face.Fig. 16 Results of mean stress (χ) and relative angle of the major principal stress (ψ) in case of (ϕ = 30°, λ = 0, β = 30°) for various degrees of negative wall friction (δ).

The stress state in the cohesionless backfill of the positively sloping wall undergoes a similar transformation to that of the previous negative sloping wall, i.e., a change from a passive stress state to an active stress state as the degree of negative wall friction increases. This transition to an active stress state was observed for a perfectly rough wall as shown in Fig. 17. Note that the Rankine region vanished from the backfill as the backfill is bounded by the line of stress characteristic on the sloping surface and the wall rear face (ref. Fig. 18). Also note that the major principal stress direction rotates in a counterclockwise direction as the magnitude of negative wall friction increases, resulting in a reduction of the mean stress at the wall rear face (Fig. 17).Fig. 17 Cartesian stress components at a height of interest in case of (ϕ = 30°, λ = 0, β = 30°) for various degrees of negative wall friction (δ): (a) δ = 0; (b) δ = − 15°; (c) δ = − 30°.

Fig. 18 The trajectories of principal stresses in negative sloping backfill in case of (ϕ = 30°, λ = 0, β = 30°) with various degrees of negative wall friction: (a) δ = 0; (b) δ = − 15°; (c) δ = − 30°.

Design charts and tables

Comprehensive design charts are presented in Fig. 19. They can be used to determine passive earth pressure coefficients with negative wall friction for various backfill geometries and internal soil friction angles. These charts also illustrate the effects of the angle of sloping backfills, the degree of negative wall friction, and the internal friction angle on the magnitude of the passive earth pressure coefficient. Notably, the presence of negative wall friction reduces the effect of the sloping angle on the earth pressure coefficient, and the value of the coefficient remains almost constant as the sloping angle increases. These findings have many practical implications for engineering applications and therefore they are presented in tabular form in Tables 3, 4 and 5.Fig. 19 Design charts for determining the passive earth pressure coefficient subjected to negative wall friction.

Table 3 Passive earth pressure coefficients for ϕ = 20° with various angles of sloping backfill (negative wall friction).

ϕ = 20°	β (degree)	
δ (degree)	− 20	− 15	− 10	− 5	0	5	10	15	20	
− 0	0.771	1.225	1.517	1.787	2.031	2.302	2.505	2.765	2.964	
− 5	0.727	1.114	1.364	1.567	1.765	1.983	2.166	2.336	2.463	
− 10	0.677	0.973	1.199	1.368	1.509	1.644	1.807	1.915	2.017	
− 15	0.620	0.869	1.019	1.135	1.258	1.330	1.410	1.488	1.566	
− 20	0.549	0.713	0.787	0.835	0.866	0.872	0.888	0.876	0.883	

Table 4 Passive earth pressure coefficients for ϕ = 30° with various angles of sloping backfill (negative wall friction).

ϕ = 30°	β (degree)	
δ (degree)	− 30	− 25	− 20	− 15	− 10	− 5	0	5	10	15	20	25	30	
− 0	0.607	1.072	1.414	1.782	2.173	2.523	3.008	3.507	3.963	4.503	5.052	5.605	6.192	
− 5	0.573	0.983	1.269	1.572	1.894	2.194	2.514	2.939	3.349	3.722	4.139	4.569	5.005	
− 10	0.540	0.892	1.140	1.385	1.640	1.883	2.150	2.397	2.738	3.063	3.348	3.639	3.961	
− 15	0.503	0.800	1.004	1.200	1.391	1.612	1.802	2.018	2.258	2.424	2.613	2.836	3.045	
− 20	0.467	0.707	0.876	1.015	1.174	1.302	1.452	1.614	1.730	1.859	2.009	2.189	2.344	
− 25	0.423	0.624	0.740	0.842	0.947	1.046	1.139	1.209	1.309	1.355	1.429	1.517	1.600	
− 30	0.376	0.508	0.578	0.629	0.671	0.704	0.722	0.738	0.737	0.751	0.743	0.756	0.744	

Table 5 Passive earth pressure coefficients for ϕ = 40° with various angles of sloping backfill (negative wall friction).

ϕ = 40°	β (degree)	
δ (degree)	− 40	− 35	− 30	− 25	− 20	− 15	− 10	− 5	0	5	10	15	20	25	30	35	40	
− 0	0.437	0.851	1.202	1.569	2.000	2.533	3.105	3.796	4.606	5.442	6.500	7.792	9.210	10.856	12.695	14.654	17.092	
− 5	0.415	0.783	1.079	1.396	1.755	2.172	2.649	3.177	3.809	4.492	5.256	6.212	7.262	8.478	9.816	11.234	12.920	
− 10	0.394	0.716	0.962	1.233	1.533	1.859	2.229	2.653	3.138	3.667	4.229	4.956	5.657	6.508	7.474	8.461	9.640	
− 15	0.370	0.653	0.860	1.081	1.324	1.581	1.874	2.183	2.553	2.918	3.372	3.849	4.326	4.953	5.601	6.264	7.055	
− 20	0.348	0.588	0.765	0.941	1.136	1.332	1.555	1.795	2.052	2.333	2.650	2.961	3.334	3.735	4.122	4.550	5.059	
− 25	0.325	0.529	0.674	0.812	0.961	1.110	1.272	1.448	1.628	1.821	2.011	2.226	2.476	2.718	2.963	3.222	3.527	
− 30	0.300	0.468	0.583	0.695	0.800	0.908	1.024	1.134	1.262	1.373	1.494	1.645	1.768	1.924	2.068	2.201	2.364	
− 35	0.274	0.407	0.495	0.573	0.636	0.721	0.795	0.846	0.930	1.001	1.069	1.149	1.204	1.271	1.333	1.412	1.490	
− 40	0.239	0.334	0.388	0.429	0.474	0.496	0.523	0.541	0.558	0.571	0.576	0.581	0.583	0.582	0.581	0.577	0.587	

Soft computing approaches

Sensitivity analysis using soft computing approach

Numerical solutions presented in the previous section have demonstrated rigor, accuracy, and efficacy, and a series of design charts and tables been also provided for practical uses. However, further insights into the load transfer mechanism of cohesionless backfill materials required more effort thus the statistical analysis is conducted by employing four well-known and robust machine learning models e.g., XGBoost32,33, decision tree (DT)34, long and short-term memory neural network (LSTM)35,36, and feedforward neural network (FNN)37,38.

The dataset used for conducting soft computing consisted of 529 samples, which were obtained by the analytical solution for both positive (found in Nguyen15) and negative (this paper) wall frictions. The statistical properties of the dataset, such as mean, standard deviation, minimum, and maximum values, are presented in Table 6. The distribution of normalized input and output variables is given in Appendix A1. These properties provide insights into the range of values for the input parameters and the resulting values for the passive earth pressure coefficient.Table 6 The statistical description of the input and output variables.

Variables	Unit	Notation	Min	Mean	Std	Skewness	Max	
Internal friction angle (ϕ)	Degree	X1	20	35	7.0271	− 0.5908	40	
Ratio sloping angle over internal friction angle (β/ϕ)	− 	X2	− 1	0	0.6366	0	1	
Ratio wall friction angle over internal friction angle (δ/ϕ)	− 	X3	− 1	0	0.6111	0.0649	1	
Passive earth pressure coefficient (Kp)	− 	Y	0.2393	2.2625	8.2270	3.9827	68.4312	

Among the machine learning models employed, XGBoost consistently demonstrated superior performance, as detailed in Appendix A3. The merits of XGBoost are evident through its exceptional performance metrics. For instance, the R2 values for the training data, testing data, and the entire dataset stand at 1, 0.998, and 1, respectively, as illustrated in Appendix A2.

XGBoost, introduced by Chen39, excels in handling large, complex datasets with high dimensionality. It outperforms other machine learning algorithms by effectively managing missing data and various data types (numeric, categorical, and ordinal variables). Notably, XGBoost is known for its speed and scalability, as confirmed by Chen and Guestrin32. In this study, therefore, XGBoost is employed to develop a soft computing model for determining passive earth pressure coefficients under varying loads and geometries. This model also facilitates sensitivity analyses, providing valuable insights into the problem’s physical behavior.

Figure 20 presents the average relative importance and the gain-based relative importance of each feature via the use of F-score. The results indicate that the internal friction angle (ϕ) has the greatest influence on the passive earth pressure coefficient, with a relative importance of over 40%, followed closely by the ratios of wall–soil friction angle over internal friction angle (δ/ϕ) and sloping angle over internal friction angle (β/ϕ) with a contribution of approximately 30%. Figure 20b shows that the gain of ϕ is approximately 12, while the gains of δ/ϕ and β/ϕ are approximately equal to 8, indicating that ϕ has the highest importance among the input parameters in predicting the target variable. This information is valuable in practice as it can quantify the influence of several design parameters on a specific design outcome such as Kp in this paper.Fig. 20 Feature importance analysis of the XGBoost model: (a) Average relative importance. (b) The gain based relative importance of each feature.

The use of partial dependence plots (PDPs) would further show the interaction between the target response and an input feature of interest. The PDPs in Fig. 21 reveal that all input parameters have a proportional impact on the output variable (passive earth pressure coefficient, Kp). That is, an increase in the value of any input parameter results in a nonlinear increase in the value of Kp, despite of the various levels of impact in each parameter. Specifically, for the internal friction angle (ϕ), the magnitude of Kp increases sharply after a value of 30 degrees (ref. Fig. 21a). The positive sloping angle (β/ϕ) also has a significant impact on increasing the value of Kp (ref. Fig. 21b).Fig. 21 PDPs of variables.

Finally, the impact of wall–soil friction is elucidated through the statistical analysis, as depicted in Fig. 22c. It becomes evident that while a positive wall friction angle (δ/ϕ) has a pronounced effect in significantly increasing the value of Kp, the presence of negative wall–soil friction diminishes the magnitude of the passive earth pressure exerted on the wall. Specifically, the value of Kp is approximately 4 for a smooth wall, but this value reduces to less than 2 when the wall–soil friction angle becomes negative with a magnitude half that of the friction angle. This observation underscores the critical need for investigating the influence of negative wall–soil friction on passive earth pressure. Such investigations are essential in our pursuit of conservative solutions that accurately account for these nuanced effects.Fig. 22 (a) Graphical User Interface (GUI) of RW: Automatic parameter setting. (b) Graphical User Interface (GUI) of RW: Model training, prediction, convergence history plotting, feature importance calculation, and partial dependence analysis. (c) Graphical User Interface (GUI) of RW: Predicted outcome.

In summary, the statistical analysis employing the soft computing model not only provides valuable insights into the behavior of passive earth pressure coefficients but also equips engineers and designers with quantifiable knowledge to optimize their designs and consider the full spectrum of factors influencing these critical structural parameters.

Graphical user interface (GUI)

As the final step of the proposed hybrid soft computing process, a graphical user interface (GUI) is implemented in this study. The developed graphical user interface (GUI), called RW, provides a user-friendly visual interface (e.g., graphical elements such as icons, buttons, menus, and windows) that allows one to interact with, rather than using traditional text-based commands. Several functions, such as automatic parameter setting, model training, prediction, convergence history plotting, feature importance calculation, and partial dependence analysis, are all created, as shown in Fig. 22a,b, and c, which provides a step-by-step guide to using the GUI.

Using the RW GUI, as shown in Fig. 22, numerical output from the user-friendly visual interface is validated with those of Kumar and Rao2 and Krabbenhoft40 as shown in Fig. 23. It is pleased to see that the outputs of the GUI are in excellent agreement with those of Krabbenhoft40 in Fig. 23a and Kumar and Rao2 in Fig. 23b with R2 = 0.997 and 0996 respectively. It can therefore be concluded that the RW GUI makes it easier for one to perform complex analyses and allows friendly interactive analyses in a more intuitive approach. It would also increase productivity in engineering design.Fig. 23 Validation with independent data (a) Krabbenhoft40 with ϕ = 20°, 25°, 30°, 35°, and 40°; δ/ϕ = 0, 1/3, ½, 2/3, and 1, β/ϕ = 0 (b) Kumar and Rao2 with ϕ = 40°, g− 1 ≤ δ/ϕ ≤ 1, β/ϕ = 0.

Conclusion

This paper has successfully studied the passive earth pressure exerted on rigid wall with negative wall friction. A subtle numerical procedure was developed to explore the statically admissible stress field in the backfill of the rigid retaining wall with various degrees of negative wall friction and backfill geometries using the fully plastic solution. In addition to the rigorous solution for the passive earth pressure on a rigid wall under negative wall–soil friction, the mechanics of load transfer in the cohesionless backfill material was also presented.

Negative wall friction has substantial effects on the load transfer mechanism, which is very distinct from that of positive wall friction. In particular, the magnitude of the negative wall friction is disproportional to the magnitude of the passive earth pressure. When the wall is fully rough, the value of Kp remains unchanged irrespective of the backfill sloping angle, and the values of Kp are all less than unity. The admissible stress fields have also revealed the mechanism of load transfer in the backfill, in that the negative friction angle would change the stress state in the non-Rankine region from passive to active stress state.

The obtained results of the passive earth pressure under negative wall friction and sloping backfill have been provided with design charts and tables to assist in engineering design practice. In addition, the incorporation of soft computing models in this study has the potential to greatly enhance our understanding of passive earth pressure and improve the design and analysis of retaining structures. The success of the study is a great contribution to the field of granular mechanics in studying the load transfer in a planar wedge, uplift anchor structures, and silo structure. It provides insight into the mechanism of load transfer of the upward relative movement between the wall and the backfill material.

In the future work, the proposed framework coupling statically admissible stress field and soft computing will be applied to investigate the seismic passive resistance in soil for negative wall–soil friction and compared with that of Choudhury and Rao11 using the limit equilibrium method.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-72094-8.

Acknowledgements

This research is funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number NCUD.02-2022.08. The authors extend heartfelt thanks to the anonymous reviewers for their invaluable time, insightful comments, and constructive feedback, which have significantly enhanced our work. Finally, gratitude is also extended to Prof. Thirapong Pipatpongsa (former Associate Professor at Kyoto University, Japan) for his valuable discussions with the corresponding author.

Author contributions

T.N.: Conceptualization, Methodology, Investigation, Validation, Visualization, Writing—Reviewing and Editing, Writing—Original draft preparation; J.S.: Methodology, Investigation, Validation, Writing—Reviewing and Editing; T.B.-N.: Methodology, Data curation, Writing- Original draft preparation.

Data availability

The data that support the findings of this study are available on request from the corresponding author.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Savage SB Yong RN McInnes D Stress discontinuities in cohesionless particulate materials Int. J. Mech. Sci. 1969 11 595 602 10.1016/0020-7403(69)90058-7
Savage, S. B., Yong, R. N. & McInnes, D. Stress discontinuities in cohesionless particulate materials. Int. J. Mech. Sci. 11, 595–602. 10.1016/0020-7403(69)90058-7 (1969).10.1016/0020-7403(69)90058-7
2. Kumar J Rao KSS Passive pressure coefficients, critical failure surface and its kinematic admissibility Géotechnique 1997 47 185 192 10.1680/geot.1997.47.1.185
Kumar, J. & Rao, K. S. S. Passive pressure coefficients, critical failure surface and its kinematic admissibility. Géotechnique 47, 185–192. 10.1680/geot.1997.47.1.185 (1997).10.1680/geot.1997.47.1.185
3. Nguyen T Shiau J Passive earth pressure in sand on inclined walls with negative wall friction based on a statically admissible stress field Acta Geotech. 2024 10.1007/s11440-024-02278-z
Nguyen, T. & Shiau, J. Passive earth pressure in sand on inclined walls with negative wall friction based on a statically admissible stress field. Acta Geotech.10.1007/s11440-024-02278-z (2024).10.1007/s11440-024-02278-z
4. Nguyen T Shiau J Ly D-K Enhanced earth pressure determination with negative wall-soil friction using soft computing Comput. Geotech. 2024 167 106086 10.1016/j.compgeo.2024.106086
Nguyen, T., Shiau, J. & Ly, D.-K. Enhanced earth pressure determination with negative wall-soil friction using soft computing. Comput. Geotech. 167, 106086. 10.1016/j.compgeo.2024.106086 (2024).10.1016/j.compgeo.2024.106086
5. Choudhury D Rao KSS Seismic uplift capacity of strip anchors in soil Geotech. Geol. Eng. 2004 22 59 72 10.1023/B:GEGE.0000014003.69378.6a
Choudhury, D. & Rao, K. S. S. Seismic uplift capacity of strip anchors in soil. Geotech. Geol. Eng. 22, 59–72. 10.1023/B:GEGE.0000014003.69378.6a (2004).10.1023/B:GEGE.0000014003.69378.6a
6. Meyerhof GG Adams JI The ultimate uplift capacity of foundations Can. Geotech. J. 1968 5 225 244 10.1139/t68-024
Meyerhof, G. G. & Adams, J. I. The ultimate uplift capacity of foundations. Can. Geotech. J. 5, 225–244. 10.1139/t68-024 (1968).10.1139/t68-024
7. Shiau JS Augarde CE Lyamin AV Sloan SW Finite element limit analysis of passive earth resistance in cohesionless soils Soils Found. 2008 48 843 850 10.3208/sandf.48.843
Shiau, J. S., Augarde, C. E., Lyamin, A. V. & Sloan, S. W. Finite element limit analysis of passive earth resistance in cohesionless soils. Soils Found. 48, 843–850. 10.3208/sandf.48.843 (2008).10.3208/sandf.48.843
8. Liu S Xia Y Liang L A modified logarithmic spiral method for determining passive earth pressure J. Rock Mech. Geotech. Eng. 2018 10 1171 1182 10.1016/j.jrmge.2018.03.011
Liu, S., Xia, Y. & Liang, L. A modified logarithmic spiral method for determining passive earth pressure. J. Rock Mech. Geotech. Eng. 10, 1171–1182. 10.1016/j.jrmge.2018.03.011 (2018).10.1016/j.jrmge.2018.03.011
9. Kérisel J Absi E Tables for the Calculation of Passive Pressure, Active Pressure and Bearing Capacity of Foundations 1990 Gauthier-Villard
Kérisel, J. & Absi, E. Tables for the Calculation of Passive Pressure, Active Pressure and Bearing Capacity of Foundations (Gauthier-Villard, 1990).
10. Terzaghi K Theoretical Soil Mechanics 1943 Wiley
Terzaghi, K. Theoretical Soil Mechanics (Wiley, 1943).
11. Choudhury D Rao KSS Seismic passive resistance in soils for negative wall friction Can. Geotech. J. 2002 39 971 981 10.1139/t02-023
Choudhury, D. & Rao, K. S. S. Seismic passive resistance in soils for negative wall friction. Can. Geotech. J. 39, 971–981. 10.1139/t02-023 (2002).10.1139/t02-023
12. Poncelet, J. (Memorial de l’officiere de genie, Paris, 1840).
13. Budhu M Soil Mechanics and Foundations 2010 Wiley
Budhu, M. Soil Mechanics and Foundations (Wiley, 2010).
14. James RG Bransby PL Experimental and theoretical investigations of a passive earth pressure problem Géotechnique 1970 20 17 37 10.1680/geot.1970.20.1.17
James, R. G. & Bransby, P. L. Experimental and theoretical investigations of a passive earth pressure problem. Géotechnique 20, 17–37. 10.1680/geot.1970.20.1.17 (1970).10.1680/geot.1970.20.1.17
15. Nguyen T Passive earth pressures with sloping backfill based on a statically admissible stress field Comput. Geotech. 2022 149 104857 10.1016/j.compgeo.2022.104857
Nguyen, T. Passive earth pressures with sloping backfill based on a statically admissible stress field. Comput. Geotech. 149, 104857. 10.1016/j.compgeo.2022.104857 (2022).10.1016/j.compgeo.2022.104857
16. Bouchaud JP Cates ME Claudin P Stress distribution in granular media and nonlinear wave equation J. Phys. 1995 I 5 639 656 10.1051/jp1:1995157
Bouchaud, J. P., Cates, M. E. & Claudin, P. Stress distribution in granular media and nonlinear wave equation. J. Phys. I(5), 639–656. 10.1051/jp1:1995157 (1995).10.1051/jp1:1995157
17. Wittmer JP Cates ME Claudin P Stress propagation and arching in static sandpiles J. Phys. 1997 I 7 39 80 10.1051/jp1:1997126
Wittmer, J. P., Cates, M. E. & Claudin, P. Stress propagation and arching in static sandpiles. J. Phys. I(7), 39–80. 10.1051/jp1:1997126 (1997).10.1051/jp1:1997126
18. Savage SB Herrmann HJ Hovi J-P Luding S Modeling and granular material boundary value problems Physics of dry granular media 1998 Springer 25 96
Savage, S. B. Modeling and granular material boundary value problems. In Physics of dry granular media (eds Herrmann, H. J. et al.) 25–96 (Springer, 1998).
19. Booker, J. R. Application of Theories of Plasticity to Cohesive Frictional Soils. (1969).
20. Rankine WJM On the stability of loose earth Philos. Trans. R. Soc. London 1857 147 9 27 10.1098/rstl.1857.0003
Rankine, W. J. M. On the stability of loose earth. Philos. Trans. R. Soc. London 147, 9–27 (1857).10.1098/rstl.1857.0003
21. Sokolovskii VV Statics of Granular Media 1965 Pergamon
Sokolovskii, V. V. Statics of Granular Media (Pergamon, 1965).
22. Lancellotta R Analytical solution of passive earth pressure Géotechnique 2002 52 617 619 10.1680/geot.2002.52.8.617
Lancellotta, R. Analytical solution of passive earth pressure. Géotechnique 52, 617–619. 10.1680/geot.2002.52.8.617 (2002).10.1680/geot.2002.52.8.617
23. Nguyen T An exact solution of active earth pressures based on a statically admissible stress field Comput. Geotech. 2023 153 105066 10.1016/j.compgeo.2022.105066
Nguyen, T. An exact solution of active earth pressures based on a statically admissible stress field. Comput. Geotech. 153, 105066. 10.1016/j.compgeo.2022.105066 (2023).10.1016/j.compgeo.2022.105066
24. Lee, E. in Elasticity: Proceedings of the Third Symposium in Applied Mathematics of the American Mathematical Society. 213 (McGraw-Hill).
25. Shield R Mixed boundary value problems in soil mechanics Q. Appl. Math. 1953 11 61 75 10.1090/qam/54512
Shield, R. Mixed boundary value problems in soil mechanics. Q. Appl. Math. 11, 61–75 (1953).10.1090/qam/54512
26. Fehlberg E Low-Order Classical Runge-Kutta Formulas with Stepsize Control and Their Application to Some Heat Transfer Problems NASA Technical Report 315 NASA Washington 1969 DC
Fehlberg, E. Low-Order Classical Runge-Kutta Formulas with Stepsize Control and Their Application to Some Heat Transfer Problems NASA Technical Report 315 NASA Washington (DC, 1969).
27. Cash JR Karp AH A variable order Runge-Kutta method for initial value problems with rapidly varying right-hand sides ACM Trans. Math. Softw. 1990 16 201 222 10.1145/79505.79507
Cash, J. R. & Karp, A. H. A variable order Runge-Kutta method for initial value problems with rapidly varying right-hand sides. ACM Trans. Math. Softw. 16, 201–222 (1990).10.1145/79505.79507
28. Press WH Teukolsky SA Adaptive stepsize Runge-Kutta integration Comput. Phys. 1992 6 188 191 10.1063/1.4823060
Press, W. H. & Teukolsky, S. A. Adaptive stepsize Runge-Kutta integration. Comput. Phys. 6, 188–191 (1992).10.1063/1.4823060
29. Nguyen T Pipatpongsa T Plastic behaviors of asymmetric prismatic sand heaps on the verge of failure Mech. Mater. 2020 151 103624 10.1016/j.mechmat.2020.103624
Nguyen, T. & Pipatpongsa, T. Plastic behaviors of asymmetric prismatic sand heaps on the verge of failure. Mech. Mater. 151, 103624. 10.1016/j.mechmat.2020.103624 (2020).10.1016/j.mechmat.2020.103624
30. Kiusalaas J Numerical Methods in Engineering with MATLAB® 2005 Cambridge University Press
Kiusalaas, J. Numerical Methods in Engineering with MATLAB® (Cambridge University Press, 2005).
31. Michalowski RL Park N Admissible stress fields and arching in piles of sand Géotechnique 2004 54 529 538 10.1680/geot.2004.54.8.529
Michalowski, R. L. & Park, N. Admissible stress fields and arching in piles of sand. Géotechnique 54, 529–538. 10.1680/geot.2004.54.8.529 (2004).10.1680/geot.2004.54.8.529
32. Chen, T. & Guestrin, C. in 22nd acm sigkdd international conference on knowledge discovery and data mining 785–794.
33. Nguyen T Ly D-K Huynh TQ Nguyen TT Soft computing for determining base resistance of super-long piles in soft soil: A coupled SPBO-XGBoost approach Comput. Geotech. 2023 162 105707 10.1016/j.compgeo.2023.105707
Nguyen, T., Ly, D.-K., Huynh, T. Q. & Nguyen, T. T. Soft computing for determining base resistance of super-long piles in soft soil: A coupled SPBO-XGBoost approach. Comput. Geotech. 162, 105707. 10.1016/j.compgeo.2023.105707 (2023).10.1016/j.compgeo.2023.105707
34. Myles AJ Feudale RN Liu Y Woody NA Brown SD An introduction to decision tree modeling J. Chemom. 2004 18 275 285 10.1002/cem.873
Myles, A. J., Feudale, R. N., Liu, Y., Woody, N. A. & Brown, S. D. An introduction to decision tree modeling. J. Chemom. 18, 275–285. 10.1002/cem.873 (2004).10.1002/cem.873
35. Ma X Tao Z Wang Y Yu H Wang Y Long short-term memory neural network for traffic speed prediction using remote microwave sensor data Transp. Res. Part C: Emerg. Technol. 2015 54 187 197 10.1016/j.trc.2015.03.014
Ma, X., Tao, Z., Wang, Y., Yu, H. & Wang, Y. Long short-term memory neural network for traffic speed prediction using remote microwave sensor data. Transp. Res. Part C: Emerg. Technol. 54, 187–197. 10.1016/j.trc.2015.03.014 (2015).10.1016/j.trc.2015.03.014
36. Nguyen T Truong TT Nguyen-Thoi T Van Hong Bui L Nguyen T-H Evaluation of residual flexural strength of corroded reinforced concrete beams using convolutional long short-term memory neural networks Structures 2022 46 899 912 10.1016/j.istruc.2022.10.103
Nguyen, T., Truong, T. T., Nguyen-Thoi, T., Van Hong Bui, L. & Nguyen, T.-H. Evaluation of residual flexural strength of corroded reinforced concrete beams using convolutional long short-term memory neural networks. Structures 46, 899–912. 10.1016/j.istruc.2022.10.103 (2022).10.1016/j.istruc.2022.10.103
37. Betkier I Oszczypała M A novel approach to traffic modelling based on road parameters, weather conditions and GPS data using feedforward neural networks Expert Syst. Appl. 2024 245 123067 10.1016/j.eswa.2023.123067
Betkier, I. & Oszczypała, M. A novel approach to traffic modelling based on road parameters, weather conditions and GPS data using feedforward neural networks. Expert Syst. Appl. 245, 123067. 10.1016/j.eswa.2023.123067 (2024).10.1016/j.eswa.2023.123067
38. Nguyen T Ly D-K Shiau J Nguyen-Dinh P Optimizing load-displacement prediction for bored piles with the 3mSOS algorithm and neural networks Ocean Eng. 2024 304 117758 10.1016/j.oceaneng.2024.117758
Nguyen, T., Ly, D.-K., Shiau, J. & Nguyen-Dinh, P. Optimizing load-displacement prediction for bored piles with the 3mSOS algorithm and neural networks. Ocean Eng. 304, 117758. 10.1016/j.oceaneng.2024.117758 (2024).10.1016/j.oceaneng.2024.117758
39. Chen T Introduction to boosted trees Univ. Washington Comput. Sci. 2014 22 14 40
Chen, T. Introduction to boosted trees. Univ. Washington Comput. Sci. 22, 14–40 (2014).
40. Krabbenhoft K Static and seismic earth pressure coefficients for vertical walls with horizontal backfill Soil Dyn. Earthq. Eng. 2018 104 403 407 10.1016/j.soildyn.2017.11.011
Krabbenhoft, K. Static and seismic earth pressure coefficients for vertical walls with horizontal backfill. Soil Dyn. Earthq. Eng. 104, 403–407. 10.1016/j.soildyn.2017.11.011 (2018).10.1016/j.soildyn.2017.11.011
