Next Article in Journal
Modeling the Within-Host Dynamics of SARS-CoV-2 Infection Based on Antiviral Treatment
Next Article in Special Issue
Analytical Solution for Bearing Capacity of Reinforced Strip Footings on Unsaturated Soils under Steady Flow
Previous Article in Journal
A Note on Some Generalized Hypergeometric Reduction Formulas
Previous Article in Special Issue
Global Sensitivity Analysis Method for Embankment Dam Slope Stability Considering Seepage–Stress Coupling under Changing Reservoir Water Levels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analytical Framework for Assessing the Unsaturated Bearing Capacity of Strip Footings under Transient Infiltration

School of Civil Engineering, Central South University, Changsha 410075, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(16), 3480; https://doi.org/10.3390/math11163480
Submission received: 18 July 2023 / Revised: 5 August 2023 / Accepted: 9 August 2023 / Published: 11 August 2023

Abstract

:
The evaluation of the bearing capacity of strip footings generally assumes that the soil is either dry or fully saturated, which contradicts the actual condition in nature where the soil is often in a partially saturated state. Furthermore, infiltration has a significant impact on the shear strength of the soil. Following the upper bound theory of the limit analysis, this article provides a theoretical framework for assessing the bearing capacity under transient flow with linear variation in infiltration intensity for the first time. Firstly, the closed form of suction stress under linear transient infiltration is derived using Laplace transform and introduced into the Mohr–Coulomb criterion. A discrete failure mechanism with fewer variables and higher accuracy is provided to ensure kinematic admissibility. The upper bound solution for bearing capacity is obtained by solving the power balance equation. The present results are compared with results from the published literature and the finite element, confirming the validity and superiority of the theoretical framework provided. A parametric analysis is also conducted on three hypothetical soil types (fine sand, silt, and clay), and the results show that unsaturated transient infiltration has a positive influence on increasing the foundation bearing capacity. The magnitude of the influence is comprehensively controlled by factors such as soil type, saturated hydraulic conductivity, infiltration intensity, infiltration time, and water table depth. The increase in bearing capacity due to unsaturated transient infiltration can be incorporated into Terzaghi’s equation as a separate component presented in tabular form for engineering design purposes.

1. Introduction

For decades, the classical problem of assessing the bearing capacity of strip footings has been extensively studied in conjunction with various practical engineering factors such as soil types, characteristics of loads, seismic forces, and seepage, as well as reinforcement with steel bars. A variety of methods have been used in this research, including experimental testing, numerical simulation, limit equilibrium analysis, limit analysis, and probabilistic analysis [1,2,3,4,5]. Among these methods, limit analysis has gained widespread application due to its computational simplicity and ability to provide reliable theoretical solutions while considering various influencing factors.
Calculating the bearing capacity of a foundation using analytical methods requires first determining its failure mechanism. Initially, failure mechanisms were classified as Prandtl [6] and Hill [7] mechanisms based on whether the foundation’s roughness was considered. Both mechanisms have a plastic zone bounded by a logarithmic spiral, but the Prandtl mechanism, which takes into account foundation friction, yields a higher bearing capacity. Building upon these concepts, Michalowski [8] proposed a method to model the failure mechanism with a large number of triangular rigid blocks, forming a multi-block failure mechanism. His approach involved minimizing the failure load with respect to all possible kinematically admissible directions at the base of each selected block, thus introducing multiple variables. Subsequently, Soubra [9] and Zhu [10] used Michalowski’s method. Soubra [9] obtained the upper-bound solution for seismic bearing capacity using the multi-block mechanism, while Zhu [10] combined the mechanism with the limit equilibrium method. To simplify the optimization process for the multi-block mechanism, Kumar [11] assumed that all the triangular blocks’ outer edges still followed a logarithmic spiral. This idea was used by Kang et al. [12] in their analysis of the bearing capacity of rock foundations. However, in Kang et al.’s mechanism, the triangular blocks no longer shared the same logarithmic spiral, but each triangle’s base corresponds to a separate logarithmic spiral.
Research on foundation bearing capacity often assumes strict conditions of extreme saturation, that is, either fully saturated or completely dry. However, soils in the natural environment are often in a partially saturated state due to factors such as rainfall, evaporation, groundwater fluctuations, and topography. Previous studies have shown that the presence of gas or gas–water interfaces and capillary forces in unsaturated soils greatly enhance their shear strength [13,14]. In order to evaluate the mechanical properties of soils under realistic conditions, Bishop [15], Fredlund et al. [16], and Vanapalli et al. [17] proposed different effective stress equations based on matric suction. However, Professor Lu and his colleagues [18,19] pointed out the limitations of using matric suction in geotechnical problems and introduced the concept of suction stress as a fundamental component for assessing the shear strength of unsaturated soils. Suction stress can provide a reasonable explanation for the properties of unsaturated soils, climate variations (infiltration/evaporation), and fluid flow behavior. Moreover, it can unify the solutions for unsaturated and saturated soils. By using the effective stress expression based on suction stress, some classical geotechnical problems such as slope stability, tunnel face stability, and earth pressure have obtained more accurate solutions [20,21,22]. Considering the influence of suction stress, the bearing capacity of the foundation can be assessed more optimistically. It is worth noting that most of the relevant literature is still based on Terzaghi’s equation. Vahedifard and Robinson [23] and Tang et al. [24] made different extensions to Terzaghi’s bearing capacity equation, considering the non-linear relationship between bearing capacity and suction stress. Du et al. [25] and Xu and Zhou [26] obtained upper-bound solutions for the bearing capacity under seepage, which still applied the original Terzaghi’s equation. Roy and Chakraborty [27], using the Hill mechanism, added the contribution of suction stress to bearing capacity as a separate item to Terzaghi’s equation.
Research on the bearing capacity of shallow foundations under unsaturated conditions, including the literature mentioned above, typically considers one-dimensional steady-state flow. In the case of steady-state flow, it is assumed that the soil moisture movement reaches an equilibrium state, where the input and output rates of water are equal. However, in practical engineering, water flow rates and hydraulic head vary over time, and the strength of the soil exhibits significant spatiotemporal variations. Events such as rainfall, snowmelt, and changes in groundwater levels are typical examples of transient infiltration. During the gradual development of transient infiltration, the unsaturated components of shear strength change with time and depth. Recently, climate change has led to more extreme rainfall events [28]. Some scholars incorporated the one-dimensional transient flow analytical solution developed by Sriastava and Yeh [29] into their research on geotechnical-related problems [30,31]. The slip line method and lower bound finite element method were used by Tan and Vanapalli [32,33] and Fathipour et al. [2] to extend the calculation for foundation bearing capacity to unsaturated soil under transient flow conditions, respectively. However, these studies under transient flow assume that the infiltration intensity is constant, which is clearly not realistic. Taking rainfall as an example, its intensity tends to fluctuate significantly over time. Among them, rainfall with linear intensity changes is the most common, including advanced-peak type, delayed-peak type, and uniform type. No researchers have yet studied the bearing capacity response under such conditions.
This paper aims to obtain more realistic solutions for the bearing capacity of strip footings under unsaturated transient infiltration. Assuming that the infiltration intensity varies linearly with time, analytical solutions for the suction stress and apparent cohesion under three types of transient infiltration are derived. The modified M-C criterion with apparent cohesion is introduced into the theoretical framework of upper bound limit analysis [34,35,36]. A discrete multi-block mechanism with fewer variables and higher accuracy is proposed for the numerical integration of apparent cohesion. After verifying the correctness of the present results using the geotechnical analysis software Geo-studio, three hypothetical soils (fine sand, silt, and clay) are selected for parametric analysis and collapse range analysis. The increase in bearing capacity due to the transient infiltration is incorporated as a separate item in Terzaghi’s equation, which is presented in tabular form for ease of use.

2. Apparent Cohesion Subject to Transient Infiltration

2.1. Suction Stress and Apparent Cohesion

The stress distribution of unsaturated soil is a critical theoretical basis for the design and analysis of strip footings. The effective stress expression based on suction stress proposed by Lu and Likos [18] has been widely applied in the study of stress distribution in unsaturated soils. The unified form is:
σ = σ u a σ s
where σ = the effective stress, σ = the total stress, u a = the pore air pressure; and σ s = the suction stress.
The introduction of suction stress avoids the uncertainty in the effective stress coefficient χ in Bishop’s theory. Moreover, the magnitude of the suction stress completely depends on the matrix suction, and the expression is as follows:
σ s = { ( u a u w ) if ( u a u w ) < 0 S e ( u a u w )   if ( u a u w ) 0
where ( u a u w ) = matric suction and S e = effective degree of saturation, which can be obtained with:
S e = θ θ r θ s θ r
where θ , θ s , and θ r are the volumetric water content, saturated volumetric content, and residual volumetric water content, respectively. To describe the shear strength characteristics of unsaturated soil, an extended M-C criterion is used, which is defined as:
τ f = c + σ tan φ         = c + ( σ u a ) tan φ σ s tan φ
where τ f = the strength of the soil and c and φ represent the effective cohesion and friction angle, respectively. In general, the additional cohesion caused by suction stress σ s tan φ in Equation (4) is called apparent cohesion c a p p , which means c a p p = σ s tan φ . The introduction of apparent cohesion c a p p unified the approach for treating dry soil and variably saturated soil, which divides the cohesion into two parts in the analysis: apparent cohesion and effective cohesion. Of course, this approach needs to be based on the assumption that the effective shear strength, including c and φ , is independent of changes in apparent cohesion c a p p and remains constant during infiltration.

2.2. Analytical Solution for Unsaturated Transient Flow

In homogeneous and isotropic soil layers, the general assumption of seepage in the soil is subject to Darcy’s law [37,38,39,40]. In that case, the one-dimensional expression of Richard’s equation reflecting transient flow is:
θ ( h m ) t = z [ k ( h m ) ( h m z + 1 ) ]
where t = the time, z = vertical distance from groundwater level, with z = l representing the earth’s surface, h m = suction head, derived from h m = ( u a u w ) γ w , and γ w = the unit weight of water, generally taking a value of 10 kN / m 3 . In order to conveniently obtain the solution for the nonlinear partial differential Equation (5), the volume of water content θ and hydraulic conductivity k are calculated according to Gardner’s model [41]:
k = k s e α h m
θ = θ r + ( θ s θ r ) e α h m
where k s = hydraulic conductivity of saturated soils and α is a parameter for unsaturated soils that represents the decrease in hydraulic conductivity and water content as the suction head increases. By substituting Equations (6) and (7) to Equation (5), a linearized partial differential equation system can be obtained, namely:
α ( θ s θ r ) k s k ( h m ) t = 2 k ( h m ) z 2 + α k ( h m ) z
To solve this equation, two boundary conditions and an initial condition are required. This article uses a steady seepage field formed by a previous infiltration intensity of q A the initial state. The two boundary conditions refer to (1) the suction stress h m at the water table equals h 0 , with h 0 = 0 . In this case, the initial suction head h m is linearly distributed in the soil. (2) When t > 0 , transient infiltration with an intensity of q B acts on the surface of the foundation ( z = l ).
q B = q 0 + m t
where q 0 = initial infiltration intensity and m = infiltration intensity variation coefficient, where m = 0 represents uniform infiltration, m > 0 represents delayed-peak infiltration, and m < 0 represents advanced-peak infiltration, as shown in Figure 1. Herein, several parameters are introduced:
K = k ( ψ ) k s , Q A = q A k s , Q 0 = q 0 k s , Z = α z , L = α l , T = α k s t θ s θ r , M = m ( θ s θ r ) α k s 2
In this way, the control Equation (8) can be simplified as:
K T = 2 K Z 2 + K Z
The boundary conditions and initial conditions can be recorded as:
{ K ( Z ,   0 ) = Q A ( Q A 1 ) e Z [ K Z + K ] Z = L = Q B K ( 0 ,   T ) = 1
An analytical solution for the normalized hydraulic conductivity K is obtained using Laplace forward and inverse transformations, with a detailed derivation process referred to Sriastava and Yeh [29] and Qin et al. [42].
K = Q A ( Q A 1 ) e Z + e L Z 2 [ M G ( t ) + ( Q 0 Q A ) F ( t ) ]
where
G ( t ) = 2 Z cosh ( Z 2 ) e L 2 + sinh ( Z 2 ) [ 2 T e L 2 4 e L cosh ( L 2 ) 2 L e L 2 ]                       + n = 1 16 cos 2 ( λ n L ) sin ( λ n Z ) sin ( λ n L ) e ( λ n 2 1 4 ) T 1 + L 2 + 2 λ n 2 L
F ( t ) = e ( L Z ) 2 e ( L + Z ) 2 + n = 1 4 sin ( λ n Z ) sin ( λ n L ) e ( λ n 2 1 4 ) T 1 + L 2 + 2 λ n 2 L
and λ n refers to the nth positive root of the following pseudo-periodic characteristic equation:
tan ( λ n L ) + 2 λ n = 0
For the uniform type of transient infiltration, substituting M = 0 into Equation (12) yields an expression for K in this case, i.e.,:
K = Q 0 ( Q 0 1 ) e Z 4 ( Q 0 Q A ) e L Z 2 e T 4 n = 1 sin ( λ n Z ) sin ( λ n L ) e λ n 2 T 1 + L / 2 + 2 λ n 2 L
The saturation degree S e , water head h m , suction stress σ s , and apparent cohesion c a p p can all be expressed as functions of K:
S e = K
u a u w = γ w ln K α
σ s = γ w K ln K α
c a p p = γ w K ln K α tan φ

3. Bearing Capacity of Strip Footings on Unsaturated Soils

3.1. Problem Description and Solution Method

As shown in Figure 2, a rigid strip footing with a width of B is resting over homogeneous, isotropic, saturated soils, and the burial depth of the water table is l. In order to simplify the bearing capacity analysis, some basic assumptions are required: (1) the load on the foundation is vertical and concentric; (2) the soil in the damaged area follows the modified Mohr–Coulomb (M-C) failure criterion; (3) the interface between the foundation material and the soil is smooth to achieve conservative estimation; and (4) the water table exceeds the damaged area and is parallel to the ground.
The upper-bound method of limit analysis has been widely used in classic soil mechanics problems such as slope stability analysis, tunnel stability analysis, and foundation bearing capacity calculation [43,44,45,46,47,48,49]. No matter how complex the problem may be, the upper-bound method can be used to obtain an actual value for the failure load, which is very practical for engineers. This method is based on the principle of virtual work, which can be expressed as follows:
A T i v i d A + V F i v i d V = V σ i j ε i j d V
where v i = the velocity of failure block, σ i j and ε i j = the internal stress and the corresponding plastic strain rate, respectively, T i and F i = the surface force and body force acting on foundations, respectively, A and V = the area and volume of the integration region, respectively.

3.2. Discrete Failure Mechanism

The failure mechanism under a static state is generally symmetrical, and the most common types are the Prandtl mechanism and the multi-block mechanism. The Prandtl mechanism only divides the failure area into three regions with fewer variables and faster calculation speed. Nevertheless, the multi-block mechanism divides the failure area into many blocks, and a smaller bearing capacity value can be obtained by optimizing multiple variables. This article proposes a new discrete failure mechanism based on these two mechanisms, simultaneously meeting the advantages of fast computation speed and small upper solutions. As shown in Figure 3, the transition zone BB1Bn+1 in the Prandtl mechanism is discretized into n triangles, each with a top angle δ θ and a bottom angle π / 2 + φ L . The active zone ABB1 is an isosceles triangle with a base angle of α , and the passive zone B1 Bn+1C is assumed to have a base angle of β , so there are only three unknown variables ( α , β , φ L ) in this new failure mechanism. The sine theorem can be applied to calculate the side length of each triangle:
l i + 1 = l i cos φ L cos ( δ θ + φ L )                 1 i n
d i = l i   sin δ θ cos ( δ θ + φ L )                 1 i n
where l i and d i are the length of line BBi and line BiBi+1, respectively, with l 1 = B sec θ / 2 . Furthermore, the length of line BC and line Bn+1C can also be calculated as l n + 2 = | B C | = l n + 1   cos φ L cos ( β + φ L ) and d n + 1 = | B n + 1 C | = l n + 1   sin β cos ( β + φ L ) , respectively.
In order to form an admissible velocity field in kinematics, each triangular block slides to one side of the surface, and its velocity v i forms an angle of φ with the line d i . Moreover, the relative velocity [ v ] i direction of adjacent triangles also has an angle of φ with the shared line l i . The velocity field is shown in Figure 4. Because of the symmetry in the mechanism, only half of its velocity field needs to be studied. Assuming the isosceles triangle ABB1 moves vertically downwards at the speed of v 0 , the velocity field can be calculated as:
v i = {   cos ( θ φ ) cos ( φ L 2 φ )   v 0                 i = 1 cos ( φ L + δ θ 2 φ ) cos ( φ L 2 φ ) v i 1           2 i n + 1
[ v ] i = {   sin ( θ + φ φ L ) cos ( φ L 2 φ )   v 0                 i = 1 sin δ θ cos ( φ L 2 φ ) v i 1         2 i n + 1

3.3. Computation of the Working Rate

Due to the fact that the cohesion of unsaturated soils is composed of effective cohesion c and apparent cohesion c app , the total internal energy dissipation D int can also be divided into two parts, namely:
D int = D e + D a
where D e and D a are the dissipation caused by effective cohesion c and apparent cohesion c app , respectively. The effective cohesion c remains a constant as time and space change, thus:
D c = i = 1 n + 1 D OB i + i = 1 n D B i B i + 1 + D B n + 1 C = c cos φ i = 1 n + 1 ( l i [ v ] i + d i v i )
However, the apparent cohesion c app caused by one-dimensional transient infiltration varies nonlinearly along the depth, so numerical integration is needed to calculate D a .
D a = i = 1 n + 1 D OB i + i = 1 n D B i B i + 1 + D B n + 1 C           = i = 1 n + 1 0 l i sin θ i c a p p [ v ] i cos φ d y sin θ i 1 2 i = 1 n [ σ B i s + σ B i + 1 s ] d i v i sin φ             + 0 d n + 1 cos ( β + φ L ) c a p p v n + 1 cos φ d y cos ( β + φ L )
where σ B i s and σ B i + 1 s represent the suction stress of points B i and B i + 1 , respectively. The coordinate system of the failure mechanism takes point O as the origin, with the vertical downward direction being the y-axis positive direction and the horizontal rightward direction being the x-axis positive direction. It is worth noting that the y value for a certain point in the underground soil is its burial depth, which is different from the z value in Equation (5). The relationship between them is y = l z . For a large value of n, the length of B i B i + 1 is short enough that its apparent cohesion c app can be seen as linearly changing.
The total power caused by the external force W ext can be obtained by accumulating the power of each part of the external force, i.e.,:
W ext = W q u + W q s + W γ
where W q u , W q 0 , and W γ = external forces caused by the vertical load of the superstructure q u , the surcharge load q s , and the soil’s self-weight within the collapse mechanism, respectively. Their detailed calculation formulas are as follows:
W q u = 1 2 q u B v 0
W q 0 = q s l n + 2 v n + 1 cos ( θ + n δ θ )
W γ = 1 2 γ i = 1 n + 1 l i d i v i cos φ L cos [ θ + ( i 1 ) δ θ ] + 1 8 γ B 2 tan θ v 0

3.4. Bearing Capacity and Optimization

By equaling external work power to internal energy power, a functional balance equation can be expressed as:
W ext = D int
Substituting Equations (26)–(32) into Equation (33) the ultimate bearing capacity q u of shallow strip footings under transient infiltration can be derived as:
q u = 2 ( D e + D a W q s W γ ) B v 0
Indeed, this solution of q u is an upper-bound solution. To measure the impact of unsaturated infiltration on the bearing capacity of shallow strip footings, this paper modified the classic Terzaghi’s expression:
q u = 1 2 γ B N γ + q s N q + c N c + q u s u c t i o n
where N γ , N q , and N c = the bearing capacity parameters related to γ , q s , and c , respectively, and q u s u c t i o n = additional bearing capacity due to transient infiltration. It should be noted that the first three terms are not influenced by transient infiltration in unsaturated soils, and the fourth term is also independent of the first three terms’ parameters and only related to the matric suction profile. When calculating one of these four items separately, the relevant parameters related to the other items are assumed to be 0. According to Equation (34), the expressions for N γ , N q , N c , and q u s u c t i o n can be derived as:
N γ = 4 W γ γ B 2 v 0
N q = 2 W q s q s B v 0
N c = 2 D e c B v 0
q u s u c t i o n = 2 D a B v 0
From the derivation process for the upper bound solution above, it can be seen that the bearing capacity is only related to three variables, namely:
q u = f ( θ , δ θ , β )
The smaller the upper bound solution, the more engineering practicality there is. In this paper, the multi-objective sequential quadratic programming (SQP) algorithm is used for optimization, with the detailed process shown in Figure 5 and the constraints listed in Table 1. The current study considers three hypothetical soils: fine sand, silt, and clay, whose typical parameters are recorded in Table 2 [30]. In all data analyses below, it is assumed that the initial infiltration intensity q A = 0 , the infiltration intensity variation coefficient m = 0 , the width of the foundation B = 1   m , the surcharge load q 0 = 10   kPa , and the unit soil gravity γ = 20   kN / m 3 , if not stated otherwise.
A convergence study is needed to determine the optimal number n of discrete triangles in the passive zone. Apparently, a larger value of n means a more accurate solution, which means a greater computational workload at the same time. So, this convergence study needs to balance these two contradictory requirements. Figure 6 takes sand as an example and obtains its bearing capacity q u corresponding to different n values according to Equation (34). From Figure 6, it can be seen that as n increases from 5 to 40, q u continues to increase, whereas the rate of increase gradually slows down. When n increases from 25 to 30, the bearing capacity only increases by 0.19%. Therefore, throughout the entire analysis process, n can be reasonably taken as 30.
This paper provides two methods to calculate the ultimate bearing capacity q u of strip footings under the action of transient infiltration. Method 1 is an individual method in which the ultimate bearing capacity q u is calculated by superimposing the contributions of the effective cohesion, the surcharge load, the soil’s self-weight, and the apparent cohesion to the bearing capacity. Its detailed procedure is to calculate N γ , N q , N c , and q u s u c t i o n according to Equations (36)–(39) and then substitute them into Equation (35). Method 2 is a joint method, which applies the Equation (34) to computer q u directly. Assuming the soil is saturated or dry, Table 3 compares the bearing capacity calculated using these two methods with the results from previous studies in the literature [9,50,51]. As can be seen, the results obtained using method 2 are greater than those obtained used method 1 for all comparisons. Furthermore, the present results obtained using both methods are smaller than the previous results, thanks to the innovative discrete mechanism in this paper. Considering the unsaturated soil subjected to transient infiltration, Table 4 compares the foundation bearing capacity q u obtained using these two methods for different types of soils. It can be more intuitively seen that method 1 tends to give more conservative results, but the difference between these two methods does not exceed 6%. Compared to the complex optimization algorithm used in method 2, method 1 can be used to calculate q u easily by directly querying N γ , N q , N c , and q u s u c t i o n for different operating modes. Therefore, method 1 is recommended for the practical design of foundations.

4. Verification

Considering that no scholars have yet studied the impact of transient infiltration on the upper-bound solution of the foundation bearing capacity, the validation in this study mainly consists of two steps. Firstly, the feasibility and accuracy of the proposed framework in obtaining the foundation bearing capacity in saturated or dry soil are tested using published research. Then, the correctness of the analytical solution for unsaturated transient infiltration is verified. The apparent cohesion generated by transient infiltration is added to the functional equation as a separate internal energy dissipation term, so as long as these two verification steps are successful, the present method for computing bearing capacity is effective.

4.1. Verification of the Bearing Capacity under the No-Suction Condition

Assuming the strip footing is located on a dry or saturated soil layer, its N γ , N c , and N q values obtained using method 1 are compared with previous results [8,9,10,51,52,53], as listed in Table 5, Table 6 and Table 7. It is worth noting that the theory used by Vesic [52] and Zhu [10] is the limit equilibrium theory, while all other solutions are upper-bound solutions. For the values of N γ recorded in Table 5, the present results gradually become the least upper-bound solutions as the φ value increases, while Chen’s results are the largest from beginning to end. The difference in the N γ value is mainly caused by the variety of collapse mechanisms. Chen [53] applied the symmetric Prandtl mechanism, whereas Michalowski [8], Soubra [9], and Zhu [10] applied the multi-block mechanism, whose periphery does not assume any shape. Apparently, the bearing capacity solutions for the multi-block mechanism are all smaller than those for the Prandtl mechanism. Moreover, both Soubra [9] and Zhu [10] observed that the average reduction in upper-bound solutions calculated using the symmetric mechanism was 23.12% compared with those obtained using the unilateral mechanism. Similarly, the limit equilibrium solutions achieved using the symmetric mechanism were found to be, on average, 2.85% smaller than those attained using the unilateral mechanism. Therefore, the symmetrical mechanism is more suitable for the calculation of bearing capacity under a static state.
From Table 6 and Table 7, it is easy to see that for all φ values, the present values of both N c and N q are the smallest. The results of Wang et al. [51] are the largest, followed by Soubra [9]. This is because Wang et al. [51] used the multi-wedge discrete mechanism, whose wedge includes quadrilaterals. In contrast, the discrete blocks in this paper and those used by Soubra [9] contain only triangles, which are more advantageous. Overall, assuming that discrete blocks have the same top and bottom angles, the multi-block mechanism proposed in this study avoids the computational complexity caused by too many variables. At the same time, it also provides a good basis for numerical integration of apparent cohesion.

4.2. Verification of the Matric Suction Profile with Infiltration Time

The finite element analysis software Geo-Studio is applied to verify the correctness of the analytical solution for transient infiltration in unsaturated soils. As shown in Figure 7, a rectangular soil layer with a length of 5 m and a height of 3 m is first established. It is assumed that the buried depth of the water table is 3 m, and the soil is fine sand. A total of 172.8 mm of rain falls on the foundation for 24 h. Prior to this rainfall, the area where the foundation is located had not experienced rainfall for a long time, so q A = 0 . The seepage analysis module SEEP/W is used to apply a boundary condition of late infiltration intensity q B at the surface. To ensure that rainfall is the same, three scenarios for linear changing rainfall intensity are set: (1) scenario 1 (advanced-peak infiltration): q B decreases from 3 × 10 6   m / s to 1 × 10 6   m / s , which means q 0 = 3 × 10 6   m / s and m = 2.315 × 10 11   m / s 2 ; (2) scenario 2 (uniform infiltration): q B = 2 × 10 6   m / s remains unchanged, which means q 0 = 2 × 10 6   m / s and m = 0 ; and (3) scenario 3 (delayed-peak infiltration): q B increases from 1 × 10 6   m / s to 3 × 10 6   m / s , which means q 0 = 1 × 10 6   m / s and m = 2.315 × 10 11   m / s 2 .
For the matric suction profile of these three scenarios, Figure 8 compares the numerical solutions using the software Geo-Studio with the analytical solutions calculated according to Equations (12) and (18). It can be easily observed that no matter the type of rainfall, the magnitude and variation in matric suction obtained using the two methods are almost the same, indicating that the analytical expression for transient seepage derived in this paper is correct. In addition, the soil matric suction along the depth changes from linear to nonlinear during rainfall. As time goes on, the magnitude of matric suction in the middle and upper part of the foundation is decreasing, and the decreasing amplitude and the affected area’s depth are increasing. Comparing these three different rainfall types, when t = 0, the distribution of matric suction is identical. When t = 6 h, the magnitude of matrix suction at all depths is ranked as advanced-peak infiltration < uniform infiltration < delayed-pack infiltration. The reason for the largest reduction in matric suction in advanced-peak infiltration lies in its maximum rainfall amount at t = 6 h. However, over time, the differences in matric suction between the three types gradually decrease, as their rainfall amounts converge toward the same. Finally, a stress analysis module SIGMA/W is superimposed based on the seepage field of t = 24 h. The load of the foundation is applied as a function of the displacement boundary condition. Note that the left side of the model is the symmetry axis of the middle line, so the actual foundation width is 1 m. Figure 9 presents a velocity vector diagram for each point in the soil layer during the foundation failure. It is seen that the failure envelope in Figure 9 is highly similar to that of the present discrete failure mechanism in Figure 3.

5. Results and Discussion

5.1. Saturation Degree and Apparent Cohesion Profiles

To explore the response of unsaturated soil under transient seepage, Figure 10, Figure 11 and Figure 12 describe the saturation degree S e and apparent cohesion c a p p distributions of fine sand, silt, and clay at different times t and with different water table depths l. It is assumed that q 0 / k s = 1 and m = 0. Obviously, the development of transient flow is affected by soil type, water table depth, infiltration time, and distance from the study site to the ground. From Figure 10a, Figure 11a and Figure 12a, as the infiltration time t increases, the saturation degree S e of all soils increases, eventually converging to 1.0 at all depths. These three soils reach complete saturation at different times. Specifically, infiltration of fine sand is a rapid process, often reaching full saturation in less than 10 days, while infiltration of silt is much slower, usually taking tens of days, and clay is the slowest, taking hundreds of days. The water table depth l also affects the process of soil saturation under the action of transient flow. It can be seen that when the water table depth is increased from 3 m to 6 m, the time required to reach the full saturation state is almost doubled. And the deeper the water table, the greater the saturation degree S e of the soil at the same depth at the same time. Furthermore, as time progresses, the soil near the surface experiences a higher increase in saturation degree S e compared to the middle soil. As a result, the upper and lower parts of the soil layer exhibit a higher saturation degree S e , while the middle part remains relatively low in saturation degree S e .
The complexity of apparent cohesion variation is undeniable, yet there are still fundamental rules that govern it. From Figure 10b, Figure 11b and Figure 12b, the apparent cohesion c a p p of the three soils will eventually converge to 0 as time goes by. Notably, the convergence time to zero cohesion aligns closely with the time required for complete soil saturation. This is actually because there is no apparent cohesion in fully saturated soil, which reaffirms the validity of the findings presented in this study. For clay, regardless of whether the water table depth is 3 m or 6 m, its apparent cohesion c a p p generally decreases over time. Moreover, the trend in c a p p along the depth is similar to that of S e . When the water table drops from 3 m to 6 m, the clay’s c a p p at the same depth and time double. However, for fine sand and silt, when the water table is shallow (l = 3 m), their apparent cohesions c a p p initially increase and then decrease with increasing depth. However, when the water table drops to 6 m, the profile of c a p p becomes disordered. These differences in cohesion variations among different soils are primarily attributed to variations in the values of α and k s , which result in different saturation degrees S e .
The apparent cohesion c a p p can be obtained by substituting the saturation into Equation (20). Differentiating Equation (20) with respect to K yields:
d c a p p d K = γ w α tan φ ( ln K + 1 )
When the saturation exceeds 1/e, the apparent cohesion c a p p decreases as the saturation degree S e . increases. Conversely, when c a p p is less than 1/e, the opposite trend occurs, as shown in Figure 13. Therefore, the maximum value of c a p p occurs when S e equals 1/e, which is consistent with the trend depicted in Figure 10, Figure 11 and Figure 12.

5.2. Parametric Analysis

The previous discussion demonstrated that the additional bearing capacity q u s u c t i o n caused by transient seepage is an independent item, unrelated to the effective cohesion c , surcharge load q s , and unit weight of soil γ . The dimensionless bearing capacity factors N γ , N q , and N c , which are recorded in Table 5, Table 6 and Table 7, are only related to the effective internal friction angle φ . Therefore, the parametric analysis in this section is specifically focused on the additional bearing capacity q u s u c t i o n . By analyzing Equation (39), it can be determined that the parameters influencing q u s u c t i o n are the permeability coefficient ( k s ), infiltration intensity ( q B ), water storage capacity ( θ s θ r ), effective internal friction angle ( φ ), desaturation coefficient ( α ), rainfall duration (t), and water table depth (l). For the purpose of simplification, the analysis assumes uniform rainfall for transient infiltration, i.e., q B = q 0 .
Figure 14, Figure 15 and Figure 16 depict the results of the additional bearing capacity q u s u c t i o n for the three soil types, involving different infiltration times t, water table depths l, and infiltration ratios q 0 / k s . From the graphs, it can be observed that at t = 0, the value of q u s u c t i o n due to matric suction in unsaturated soils does not vary with changes in q 0 / k s . Among these three soil types, clay exhibits the highest value of q u s u c t i o n , followed by silt, and then fine sand. For example, when t = 0 days and l = 6 m, the bearing capacity of the foundation increases by 23.60 kPa for fine sand, 36.07 kPa for silt, and 143.80 kPa for clay. Furthermore, as the water table depth increases from 4 m to 8 m, the initial response of q u s u c t i o n differs for different soil types. Fine sand and silt present an increasing trend, with the values decreasing from 59.04 kPa to 7.96 kPa and from 61.01 kPa to 18.26 kPa, respectively. However, clay consistently increases from 115.73 kPa to 152.52 kPa. With the development of transient seepage, the influence of water table variations in q u s u c t i o n becomes disordered for fine sand and silt. However, the value of q u s u c t i o n for clay continues to increase as the water table drops, aligning with the trend in the apparent cohesion variation observed in Figure 10, Figure 11 and Figure 12.
It is important to pay attention to the influence of the infiltration ratio q 0 / k s on the additional bearing capacity q u s u c t i o n of the three soils. Theoretically, the rainfall flux rate is always less than or equal to the saturated hydraulic conductivity of the soil, i.e., q 0 / k s 1 . Previous studies [2,30,31] have typically assumed that q 0 / k s = 1 to simplify the analysis. In practical engineering, due to uncertainties in rainfall and surface water, the infiltration ratio q 0 / k s can take any value between 0 and 1. Therefore, this study considers six scenarios by extending the investigation to different infiltration ratios of q 0 / k s = 0.05, 0.25, 0.10, 0.50, 0.75, and 1. From Figure 14, Figure 15 and Figure 16, it can be observed that for a given water table depth and infiltration ratio q 0 / k s , the q u s u c t i o n for fine sand and silt show a similar trend over time, while the q u s u c t i o n for clay monotonically decreases over time. Specifically, for fine sand and silt, when q 0 / k s 0.50 , their q u s u c t i o n value initially increases and then decreases as time goes by. Moreover, smaller values of q 0 / k s result in larger maximum values of q u s u c t i o n and longer corresponding infiltration times t. The maximum value of q u s u c t i o n at q 0 / k s = 0.50 is the highest among all scenarios, with magnitudes of 91.15 kPa and 70.44 kPa, respectively. The maximum values of q u s u c t i o n for fine sand typically occur between 0.5 days and 1 day, while for silt, they occur between 2 days and 4 days. Additionally, in the descending part of the additional bearing capacity over time, a smaller q 0 / k s generates a slower rate of decrease. When q 0 / k s = 0.05 ,   0.25 ,   0.10 , the q u s u c t i o n of fine sand and silt monotonically increase with time, but the rate of increase gradually slows down. Furthermore, a larger q 0 / k s leads to a greater magnitude of increase in q u s u c t i o n . The response of clay’s q u s u c t i o n to the development of transient seepage is completely different. As time goes by, the q u s u c t i o n corresponding to all values of q 0 / k s decreases, and the range of reduction increases with an increase in q 0 / k s .
The reason for the different trends corresponding to different infiltration ratios q 0 / k s is that smaller rainfall flux rates continuously increase the saturation degree of fine sand and silt, but it takes a longer time to reach the inflection point of 1/e, resulting in an increasing trend. However, when the rainfall intensity is larger, the saturation degree quickly reaches 1/e and continues to increase, leading to an initial increase and subsequent decrease in apparent cohesion c a p p . On the other hand, clay has a saturation degree greater than 1/e initially, so its apparent cohesion c a p p will decrease regardless of the infiltration flux.
Assuming q 0 / k s = 0.5 and k s = 5 × 10 6 , Figure 17, Figure 18 and Figure 19 depict the variation in additional bearing capacity q u s u c t i o n with respect to the effective internal friction angle φ for different values of the water table depth l, water storage capacity θ s θ r , and desaturation coefficient α . As the water table decreases, q u s u c t i o n slightly increases. For smaller values of α , an increase in θ s θ r also leads to an increase in q u s u c t i o n , although the effect is minor. Clearly, the effective internal friction angle φ and desaturation coefficient α have more significant impacts on the additional bearing capacity. With an increase in φ , the q u s u c t i o n continuously increases, and the speed of the increase accelerates. The influence of φ on q u s u c t i o n becomes more evident for smaller values of α . For example, when θ s θ r = 0.5 , α = 0.2 , and l = 8 m, the value of q u s u c t i o n increases from 27.05 kPa to 1161.52 kPa with increasing φ from 10° to 40°. From the equation c a p p = σ s tan φ , it can be determined that an increase in φ can directly increase the c a p p , thereby enhancing q u s u c t i o n . Additionally, as one of the most important parameters of soil, the value of φ plays a critical role not only in bearing capacity neglecting suction stress but also in seismic bearing capacity [54]. Equations (6) and (7) indicate that the soil water characteristic curve and hydraulic conductivity are both influenced by α . This explains the high sensitivity of additional bearing capacity q u s u c t i o n to the desaturation coefficient α .

5.3. Collapse Mechanism and Design Tables

The outer contours of the proposed discrete failure mechanisms generated for three kinds of soils with different infiltration times at three different water table depths (l = 4 m, 6 m, and 8 m) are presented in Figure 20. Without considering suction stress, there is no doubt that the bearing capacity of these three kinds of soil is ordered as clay > silt > fine sand. The same is true for the collapse range ordering shown in Figure 20. In fact, considering the effect of transient seepage, the three kinds of soils remain unchanged when l = 4 m and 6 m in order of bearing capacity, as shown in Figure 20. When the water level drops to 8 m, the bearing capacity of fine sand surpasses that of silt, but it is still smaller than that of clay. Furthermore, it can be observed that the water table depth significantly affects the bearing capacity while exerting minimal influence on the collapse range.
Since the contribution of transient infiltration is a separate part of the ultimate bearing capacity q u , Table 8, Table 9 and Table 10 record the addition bearing capacity q u s u c t i o n for various combinations of φ , t, l, and ks with the infiltration ratio q 0 / k s being equal to 1.00, 0.75, and 0.50. Applying method 1, these q u s u c t i o n values can be substituted with N γ , N q , and N c in Table 5, Table 6 and Table 7 into Equation (35) to calculate the ultimate bearing capacity q u .

6. Conclusions

Transient infiltration, such as heavy rainfall, can lead to dramatic changes in the degree of saturation and matrix suction within the soil, thus affecting the foundation bearing capacity. In this study, a theoretical framework for considering the contribution of transient infiltration to foundation bearing capacity is provided for the first time. An analytical solution for the transient flow of linear infiltration is derived to express the soil saturation. Then, suction stress and apparent cohesion are introduced to modify the M-C criterion. Within the framework of the kinematic approach in limit analysis, a new discrete failure mechanism suitable for considering unsaturated effects is proposed, which possesses the characteristics of fewer variables and higher accuracy. After that, the power balance equation is established, and the SQP algorithm is applied to obtain the rigorous upper-bound solution for the bearing capacity. The effectiveness and rationality of the proposed theoretical framework are fully demonstrated using comparisons with the results of previous studies and finite element results. Due to the high uncertainty in unsaturated soil properties, three hypothetical soil materials are selected for analysis. According to the comparison and parametric study, several significant conclusions are drawn as follows:
(1)
Two methods for calculating the bearing capacity under transient infiltration are provided. Method 1 is an individual method, which adds an additional bearing capacity item q u s u c t i o n to the three bearing capacity items proposed by Tarzaghi and optimizes each item individually before superimposing it. Method 2 is a joint method, which directly optimizes the objective function in Equation (34) to obtain bearing capacity q u directly. A comparison of the results shows that the upper-bound solution for method 1 is smaller than that for method 2; therefore, method 1 is more conservative. The q u s u c t i o n values under different combinations of parameters are recorded in Table 8, Table 9 and Table 10.
(2)
The additional bearing capacity q u s u c t i o n due to the transient infiltration has significant temporal variability. In the initial state (t = 0), the unsaturated effect is greatest for clay, followed by silt and fine sand. For example, at l = 6 m, taking into account the effect of transient infiltration, the bearing capacity of the fine sand, silt, and clay increases by 23.60 kPa, 36.07 kPa, and 143.80 kPa, respectively. However, after a certain period of transient infiltration, the q u s u c t i o n will converge to zero for all three soils. This process takes only a few days for fine sand, tens of days for silt, and hundreds of days for clay. Notably, the speed of infiltration depends on the saturated hydraulic conductivity k s . Increasing the value of k s will accelerate the increase in the saturation degree and the gradual loss of apparent cohesion.
(3)
The infiltration ratio q 0 / k s also has a crucial effect on the additional bearing capacity q u s u c t i o n . For fine sand and silt subjected to uniform rainfall, the q u s u c t i o n increases and then decreases over time, thus presenting a clear local maximum. The maximum value occurs between 0.5 d and 1 d for fine sand and between 2 d and 4 d for silt. For the same infiltration time, a larger value of q 0 / k s means a smaller value of q u s u c t i o n in clay. The core effect of q 0 / k s on q u s u c t i o n lies its influence on the soil saturation degree S e . The apparent cohesion increases as S e increases from 0 to 1/e and decreases as S e increases from 1/e to 1. Therefore, the additional bearing capacity of soil achieves its peak value at S e = 1 / e .

Author Contributions

Methodology, S.X.; Resources, D.Z.; Data curation, S.X.; Writing—original draft, S.X.; Writing—review & editing, D.Z.; Supervision, D.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data available on request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bildik, S.; Laman, M. Experimental investigation of the effects of pipe location on the bearing capacity. Geomech. Eng. 2015, 8, 221–235. [Google Scholar] [CrossRef]
  2. Fathipour, H.; Tajani, S.B.; Payan, M.; Chenari, R.J.; Senetakis, K. Impact of transient infiltration on the ultimate bearing capacity of obliquely and eccentrically loaded strip footings on partially saturated soils. Int. J. Geomech. 2023, 23, 04022290. [Google Scholar] [CrossRef]
  3. Soufi, G.R.; Chenari, R.J.; Javankhoshdel, S. Conventional vs. modified pseudo-dynamic seismic analyses in the shallow strip footing bearing capacity problem. Earthq. Eng. Eng. Vib. 2021, 20, 993–1006. [Google Scholar] [CrossRef]
  4. Liu, J.; Xu, S.; Yang, X.L. Modified pseudo–dynamic bearing capacity of strip footing on rock masses. Comput. Geotech. 2022, 150, 104897. [Google Scholar] [CrossRef]
  5. Griffiths, D.V.; Fenton, G.A.; Manoharan, N. Bearing capacity of rough rigid strip footing on cohesive soil: Probabilistic study. J. Geotech. Geoenvironmental Eng. 2002, 128, 743–755. [Google Scholar] [CrossRef]
  6. Prandtl, L. Eindringungsfestigkeit und festigkeit von schneiden. Angew Math UMech. 1921, 1, 15–20. [Google Scholar] [CrossRef]
  7. Hill, R. Mathematical Theory of Plasticity; Clarendon Press: Oxford, UK, 1950. [Google Scholar]
  8. Michalowski, R. An estimate of the influence of soil weight on bearing capacity using limit analysis. Soils Found. 1997, 37, 57–64. [Google Scholar] [CrossRef] [Green Version]
  9. Soubra, A.-H. Upper-bound solutions for bearing capacity of foundations. J. Geotech. Geoenvironmental Eng. 1999, 125, 59–68. [Google Scholar] [CrossRef] [Green Version]
  10. Zhu, D. The least upper-bound solutions for bearing capacity factor Nγ. Soils Found. 2000, 40, 123–129. [Google Scholar] [CrossRef] [Green Version]
  11. Kumar, J. Effect of footing—Soil interface friction on bearing capacity factor Nγ. Géotechnique 2004, 54, 677–680. [Google Scholar] [CrossRef]
  12. Kang, X.D.; Zhu, J.Q.; Yang, X.L. Seismic bearing capacity of rock foundations subjected to seepage by a unilateral piece-wise log-spiral failure mechanism. Comput. Geotech. 2023, 158, 105363. [Google Scholar] [CrossRef]
  13. Oh, W.T.; Vanapalli, S.K. Interpretation of the bearing capacity of unsaturated fine-grained soil using the modified effective and the modified total stress approaches. Int. J. Geomech. 2013, 13, 769–778. [Google Scholar] [CrossRef]
  14. Oh, W.T.; Vanapalli, S.K. Modelling the applied vertical stress and settlement relationship of shallow foundations in saturated and unsaturated sands. Can. Geotech. J. 2011, 48, 425–438. [Google Scholar] [CrossRef] [Green Version]
  15. Bishop, A.W. The principal of effective stress. Tek. Ukebl. 1959, 39, 859–863. [Google Scholar]
  16. 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]
  17. 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]
  18. Lu, N.; Likos, W.J. Unsaturated Soil Mechanics; John Wiley and Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
  19. Lu, N.; Godt, J.W.; Wu, D.T. A closed-form equation for effective stress in unsaturated soil. A closed-form equation for effective stress in unsaturated soil. Water Resour. Res. 2010, 46, W05515. [Google Scholar] [CrossRef] [Green Version]
  20. Griffiths, D.V.; Lu, N. Unsaturated slope stability analysis with steady infiltration or evaporation using elasto-plastic finite elements. Int. J. Numer. Anal. Methods Geomech. 2005, 29, 249–267. [Google Scholar] [CrossRef]
  21. Gong, W.; Zhao, C.; Juang, C.H.; Zhang, Y.; Lu, Y. Coupled characterization of stratigraphic and geo-properties uncertainties—A conditional random field approach. Eng. Geol. 2021, 294, 106348. [Google Scholar] [CrossRef]
  22. Vahedifard, F.; Leshchinsky, B.A.; Mortezaei, K.; Lu, N. Active earth pressures for unsaturated retaining structures. J. Geotech. Geoenvironmental Eng. 2015, 141, 04015048. [Google Scholar] [CrossRef]
  23. Vahedifard, F.; Robinson, J.D. Unified Method for Estimating the Ultimate Bearing Capacity of Shallow Foundations in Variably Saturated Soils under Steady Flow. J. Geotech. Geoenviron. Eng. 2016, 142, 04015095. [Google Scholar] [CrossRef]
  24. Tang, Y.; Taiebat, H.A.; Senetakis, K. Effective stress based bearing capacity equations for shallow foundations on unsaturated soils. J. Geoeng. 2017, 12, 59–64. [Google Scholar]
  25. Du, D.; Zhuang, Y.; Sun, Q.; Yang, X.; Dias, D. Bearing capacity evaluation for shallow foundations on unsaturated soils using discretization technique. Comput. Geotech. 2021, 137, 104309. [Google Scholar] [CrossRef]
  26. Xu, S.; Zhou, D. Seismic Bearing capacity solution for strip footings in unsaturated soils with modified pseudo-dynamic approach. Mathematics 2023, 11, 2692. [Google Scholar] [CrossRef]
  27. Roy, S.; Chakraborty, M. Unsaturated bearing capacity of strip foundations by using the upper bound rigid block method. Comput. Geotech. 2023, 156, 105260. [Google Scholar] [CrossRef]
  28. Ragno, E.; Agha Kouchak, A.; Love, C.A.; Cheng, L.; Vahedifard, F.; Lima, C.H. Quantifying changes in future intensity-duration-frequency curves using multimodel ensemble simulations. Water Resour. Res. 2018, 54, 1751–1764. [Google Scholar] [CrossRef]
  29. Srivastava, R.; Yeh, T.C.J. Analytical solutions for one-dimensional, transient infiltration toward the water table in homogeneous and layered soils. Water Resour. Res. 1991, 27, 753–762. [Google Scholar] [CrossRef]
  30. Li, T.Z.; Yang, X.L. New approach for face stability assessment of tunnels driven in nonuniform soils. Comput. Geotech. 2020, 121, 103412. [Google Scholar] [CrossRef]
  31. Li, Z.W.; Yang, X.L. Three-dimensional active earth pressure under transient unsaturated flow conditions. Comput. Geotech. 2020, 123, 103559. [Google Scholar] [CrossRef]
  32. Tan, M.; Vanapalli, S.K. Foundation bearing capacity estimation on unsaturated soil slope under transient flow condition using slip line method. Comput. Geotech. 2022, 148, 104804. [Google Scholar] [CrossRef]
  33. Tan, M.; Vanapalli, S.K. Failure envelops for foundation subjected to inclined and eccentric loading considering steady state and transient flow conditions in unsaturated soils. Comput. Geotech. 2023, 157, 105315. [Google Scholar] [CrossRef]
  34. Yang, X.L.; Yin, J.H. Slope stability analysis with nonlinear failure criterion. J. Eng. Mech. 2004, 130, 267–273. [Google Scholar] [CrossRef]
  35. Yang, X.L.; Huang, F. Collapse mechanism of shallow tunnel based on nonlinear Hoek–Brown failure criterion. Tunn. Undergr. Space Technol. 2011, 26, 686–691. [Google Scholar] [CrossRef]
  36. Yang, X.L.; Li, L.; Yin, J.H. Seismic and static stability analysis for rock slopes by a kinematical approach. Geotechnique 2004, 54, 543–549. [Google Scholar] [CrossRef]
  37. Smirnov, N.N.; Nikitin, V.F.; Skryleva, E.I. Microgravity investigation of seepage flows in porous media. Microgravity Sci. Technol. 2019, 31, 629–639. [Google Scholar] [CrossRef]
  38. Dushin, V.R.; Smirnov, N.N.; Nikitin, V.F.; Skryleva, E.I.; Weisman, Y.G. Multiple capillary-driven imbibition of a porous medium under microgravity conditions: Experimental investigation and mathematical modeling. Acta Astronaut. 2022, 193, 572–578. [Google Scholar] [CrossRef]
  39. Smirnova, M.N.; Nikitin, V.F.; Skryleva, E.I.; Weisman, Y.G. Capillary driven fluid flows in microgravity. Acta Astronaut. 2023, 204, 892–899. [Google Scholar] [CrossRef]
  40. Smirnov, N.N.; Legros, J.C.; Nikitin, V.F.; Istasse, E.; Schramm, L.; Wassmuth, F.; Hart, D.A. Filtration in artificial porous media and natural sands under microgravity conditions. Microgravity-Sci. Technol. 2003, 14, 3–28. [Google Scholar] [CrossRef]
  41. Gardner, W.R. Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table. Soil Sci. 1958, 85, 228–232. [Google Scholar] [CrossRef]
  42. Qin, W.X.; Deng, C.X.; Hu, H.; Xiong, X.Y.; Liu, Z.C. Analytical solution of transient ultimate bearing capacity for foundation under delayed-peak pattern rainfall infiltration. Adv. Sci. Technol. Water Resour. 2023, 43, 33–38. (In Chinese) [Google Scholar]
  43. Zhang, Z.L.; Yang, X.L. Unified solution of safety factors for three-dimensional compound slopes considering local and global instability. Comput. Geotech. 2023, 155, 105227. [Google Scholar] [CrossRef]
  44. Zhang, Z.L.; Zhu, J.Q.; Yang, X.L. Three-dimensional active earth pressures for unsaturated backfills with cracks considering steady seepage. Int. J. Geomech. 2023, 23, 04022270. [Google Scholar] [CrossRef]
  45. Xu, S.; Liu, J.; Yang, X.L. Pseudo-dynamic analysis of a 3D tunnel face in inclined weak strata. Undergr. Space 2023, 12, 156–166. [Google Scholar] [CrossRef]
  46. Chen, B.H.; Luo, W.J.; Xu, X.Y.; Hu, R.Q.; Yang, X.L. Seismic bearing capacity of strip footing with nonlinear Mohr–Coulomb failure criterion. Int. J. Geomech. 2022, 22, 06022029. [Google Scholar] [CrossRef]
  47. Zhong, J.H.; Yang, X.L. Kinematic analysis of the three-dimensional stability for tunnel faces by pseudodynamic approach. Comput. Geotech. 2021, 128, 103802. [Google Scholar] [CrossRef]
  48. Yang, X.L.; Wang, J.M. Ground movement prediction for tunnels using simplified procedure. Tunn. Undergr. Space Technol. 2011, 26, 462–471. [Google Scholar] [CrossRef]
  49. Yang, X.L.; Yin, J.H. Upper bound solution for ultimate bearing capacity with a modified Hoek-Brown failure criterion. Int. J. Rock Mech. Min. Sci. 2005, 42, 550–560. [Google Scholar] [CrossRef]
  50. Jin, L.; Zhang, H.; Feng, Q. Ultimate bearing capacity of strip footing on sands under inclined loading based on improved radial movement optimization. Eng. Optim. 2021, 53, 277–299. [Google Scholar] [CrossRef]
  51. Wang, Y.J.; Yin, J.H.; Chen, Z.Y. Calculation of bearing capacity of a strip footing using an upper bound method. Int. J. Numer. Anal. Methods Geomech. 2001, 25, 841–851. [Google Scholar] [CrossRef]
  52. Vesic, A.S. Bearing Capacity of Shallow Foundations, Foundation Engineering Handbook; Winterkorn, F.S., Fand, H.Y., Eds.; Springer: Berlin/Heidelberg, Germany, 1975. [Google Scholar]
  53. Chen, W.F. Limit Analysis and Soil Plasticity; Elsevier: Amsterdam, The Netherlands, 1975. [Google Scholar]
  54. Zhong, J.; Hou, C.; Yang, X. Bearing capacity of foundations resting on rock masses subjected to Rayleigh waves. Soil Dyn. Earthq. Eng. 2023, 167, 107791. [Google Scholar] [CrossRef]
Figure 1. Finite element model for a foundation under rainfall infiltration.
Figure 1. Finite element model for a foundation under rainfall infiltration.
Mathematics 11 03480 g001
Figure 2. Transient seepage in the homogeneous foundation.
Figure 2. Transient seepage in the homogeneous foundation.
Mathematics 11 03480 g002
Figure 3. Discrete collapse mechanism for strip footings under unsaturated transient infiltration.
Figure 3. Discrete collapse mechanism for strip footings under unsaturated transient infiltration.
Mathematics 11 03480 g003
Figure 4. (a) Velocity field for the discrete collapse mechanism. (b) Velocity hodographs.
Figure 4. (a) Velocity field for the discrete collapse mechanism. (b) Velocity hodographs.
Mathematics 11 03480 g004
Figure 5. Flowchart for obtaining the optimal upper solution of bearing capacity.
Figure 5. Flowchart for obtaining the optimal upper solution of bearing capacity.
Mathematics 11 03480 g005
Figure 6. The magnitude of q c e versus different numbers of rigid triangular blocks n.
Figure 6. The magnitude of q c e versus different numbers of rigid triangular blocks n.
Mathematics 11 03480 g006
Figure 7. Finite element model for a foundation under rainfall infiltration.
Figure 7. Finite element model for a foundation under rainfall infiltration.
Mathematics 11 03480 g007
Figure 8. Distribution of matric suction along depth under (a) advanced-peak infiltration, (b) uniform infiltration, and (c) delayed-peak infiltration.
Figure 8. Distribution of matric suction along depth under (a) advanced-peak infiltration, (b) uniform infiltration, and (c) delayed-peak infiltration.
Mathematics 11 03480 g008
Figure 9. Soil velocity vector diagram during foundation failure.
Figure 9. Soil velocity vector diagram during foundation failure.
Mathematics 11 03480 g009
Figure 10. Profiles of (a) effective saturation Se and (b) apparent cohesion c a p p for fine sand.
Figure 10. Profiles of (a) effective saturation Se and (b) apparent cohesion c a p p for fine sand.
Mathematics 11 03480 g010
Figure 11. Profiles of (a) effective saturation Se and (b) apparent cohesion c a p p for silt.
Figure 11. Profiles of (a) effective saturation Se and (b) apparent cohesion c a p p for silt.
Mathematics 11 03480 g011
Figure 12. Profiles of (a) effective saturation Se and (b) apparent cohesion c a p p for clay.
Figure 12. Profiles of (a) effective saturation Se and (b) apparent cohesion c a p p for clay.
Mathematics 11 03480 g012
Figure 13. Apparent cohesion versus S e for different types of soils.
Figure 13. Apparent cohesion versus S e for different types of soils.
Mathematics 11 03480 g013
Figure 14. Additional bearing capacity q u s u c t i o n versus infiltration time t for l = 5 m in (a) fine sand, (b) silt, and (c) clay.
Figure 14. Additional bearing capacity q u s u c t i o n versus infiltration time t for l = 5 m in (a) fine sand, (b) silt, and (c) clay.
Mathematics 11 03480 g014
Figure 15. Additional bearing capacity q u s u c t i o n versus infiltration time t for l = 6 m in (a) fine sand, (b) silt, and (c) clay.
Figure 15. Additional bearing capacity q u s u c t i o n versus infiltration time t for l = 6 m in (a) fine sand, (b) silt, and (c) clay.
Mathematics 11 03480 g015
Figure 16. Additional bearing capacity q u s u c t i o n versus infiltration time t for l = 8 m in (a) fine sand, (b) silt, and (c) clay.
Figure 16. Additional bearing capacity q u s u c t i o n versus infiltration time t for l = 8 m in (a) fine sand, (b) silt, and (c) clay.
Mathematics 11 03480 g016
Figure 17. Additional bearing capacity q u s u c t i o n versus φ for different values of α for l = 4 m: (a) θ s θ r = 0.3 , (b) θ s θ r = 0.4 , and (c) θ s θ r = 0.5 .
Figure 17. Additional bearing capacity q u s u c t i o n versus φ for different values of α for l = 4 m: (a) θ s θ r = 0.3 , (b) θ s θ r = 0.4 , and (c) θ s θ r = 0.5 .
Mathematics 11 03480 g017
Figure 18. Additional bearing capacity q u s u c t i o n versus φ for different values of α for l = 6 m: (a) θ s θ r = 0.3 , (b) θ s θ r = 0.4 , and (c) θ s θ r = 0.5 .
Figure 18. Additional bearing capacity q u s u c t i o n versus φ for different values of α for l = 6 m: (a) θ s θ r = 0.3 , (b) θ s θ r = 0.4 , and (c) θ s θ r = 0.5 .
Mathematics 11 03480 g018
Figure 19. Additional bearing capacity q u s u c t i o n versus φ for different values of α for l = 8 m: (a) θ s θ r = 0.3 , (b) θ s θ r = 0.4 , and (c) θ s θ r = 0.5 .
Figure 19. Additional bearing capacity q u s u c t i o n versus φ for different values of α for l = 8 m: (a) θ s θ r = 0.3 , (b) θ s θ r = 0.4 , and (c) θ s θ r = 0.5 .
Mathematics 11 03480 g019
Figure 20. The critical slip surface of the three soils for (a) l = 4 m, (b) l = 6 m, and (c) l = 8 m.
Figure 20. The critical slip surface of the three soils for (a) l = 4 m, (b) l = 6 m, and (c) l = 8 m.
Mathematics 11 03480 g020
Table 1. Limitations to be embedded in the algorithm.
Table 1. Limitations to be embedded in the algorithm.
Constraints/Conditions
Geometric compatibility θ + n δ θ + β = π and φ L + β < π / 2
Kinematic admissibility v i + 1 > v i and v i > 0
Range of values 0 < θ ,   δ θ ,   β < π / 2
Table 2. Input parameters for the four types of soils used in present study.
Table 2. Input parameters for the four types of soils used in present study.
Soil Type α   ( m 1 ) k s   ( m s 1 ) θ s θ r c (kPa) φ
Fine sand0.7 5.0 × 10 6 0.410.05030°
Silt0.5 9.0 × 10 7 0.450.10525°
Clay0.13 5.0 × 10 8 0.580.051020°
Table 3. Comparison of bearing capacity with previous studies.
Table 3. Comparison of bearing capacity with previous studies.
MethodThis StudySoubra [9]Jin et al. [50]Wang et al. [51]
Method 1: q sup er   ( kN / m 2 ) 675.18680.58686.76704.89
Method 2: q u   ( kN / m 2 ) 718.30726.13728.03736.19
Note: γ = 18   kN / m 3 , φ = 30 , and c = 10   kPa , q s = 10   kPa .
Table 4. Comparison of bearing capacity using the two methods proposed in this study.
Table 4. Comparison of bearing capacity using the two methods proposed in this study.
Fine sand (t = 1 d)Silt (t = 4 d)Clay (t = 60 d)
φ Method 1Method 2DifferenceMethod 1Method 2DifferenceMethod 1Method 2Difference
q super ( kN / m 2 ) q u ( kN / m 2 ) (%) q super ( kN / m 2 ) q u ( kN / m 2 ) (%) q super ( kN / m 2 ) q u ( kN / m 2 ) (%)
1576.8179.112.91138.21142.162.78219.95225.102.29
20136.37143.124.71222.35233.164.64345.48359.073.79
25248.64265.405.92373.14396.615.92563.62592.794.92
30477.84508.495.86665.66707.915.97970.851025.205.30
35968.381027.255.731266.861347.005.951777.811886.005.74
402118.272241.965.522625.712791.375.933301.023367.011.96
Note: l = 5   m .
Table 5. Comparison of N γ with previous studies.
Table 5. Comparison of N γ with previous studies.
φ ( ) This StudyWang et al. [51] Soubra [9]Michalowski [8] Chen [53]Zhu [10]
SymmetricalOne-SidedSymmetricalOne-SidedSymmetricalSymmetricalSymmetricalOne-Sided
152.493.531.952.101.942.941.942.10
204.896.564.494.674.476.204.474.66
259.8912.269.8110.069.7712.969.7610.03
3021.4224.2121.5121.8821.3927.6721.3821.81
3548.3650.9449.0049.6248.6861.4748.6549.38
40117.17122.95119.84120.96118.83145.19118.76120.15
45315.98331.22326.59328.88322.84374.02322.62325.77
Table 6. Comparison of N c with previous studies.
Table 6. Comparison of N c with previous studies.
φ ( ) This StudyWang et al. [51]Soubra [9]Vesic [52]
SymmetricalOne-SidedSymmetricalOne-SidedSymmetrical
1510.9411.3810.9911.0010.98
2014.7515.0014.8614.8714.83
2520.5620.8020.7720.7820.72
3029.8330.2030.2430.2530.14
3545.5246.5046.3346.3546.12
4074.0775.9075.7775.8075.12
45131.12141.38134.99135.09133.88
Table 7. Comparison of N q with previous studies.
Table 7. Comparison of N q with previous studies.
φ ( ) This StudyWang et al. [51]Soubra [9]Vesic [52]
SymmetricalOne-SidedSymmetricalOne-SidedSymmetrical
153.943.943.953.953.94
206.406.476.416.416.40
2510.6710.7210.6910.6910.66
3018.4118.5018.4618.4618.40
3533.3033.3133.4433.4333.30
4064.2664.9064.5864.5564.20
45134.53141.38135.99135.91134.87
Table 8. Additional bearing capacity q u s u c t i o n due to transient infiltration ( q 0 / k s = 1 ).
Table 8. Additional bearing capacity q u s u c t i o n due to transient infiltration ( q 0 / k s = 1 ).
φ ( ) l ( m ) k s = 5 × 10 6 k s = 5 × 10 7 k s = 5 × 10 8
t = 0 dt = 0.5 dt = 1 dt = 2 dt = 4 dt = 4 dt = 8 dt = 12 dt = 16 dt = 60 dt = 120 dt = 180 dt = 240 d
10412.67 11.81 8.89 4.59 1.07 12.40 10.01 7.84 6.04 11.21 7.84 5.27 3.46
68.95 13.27 11.65 8.19 3.72 13.42 12.38 10.91 9.48 13.03 10.91 8.82 7.04
85.48 13.40 12.52 9.65 5.44 13.28 13.02 11.95 10.77 13.36 11.95 10.20 8.62
20447.02 43.88 33.26 17.26 4.04 45.94 37.40 29.40 22.67 41.72 29.40 19.81 13.03
633.66 48.87 43.53 30.87 14.09 49.12 46.06 40.88 35.67 48.17 40.88 33.21 26.58
820.72 48.74 46.58 36.35 20.64 47.77 48.18 44.63 40.44 48.96 44.63 38.37 32.55
304153.96 143.56 110.22 57.62 13.54 149.68 123.50 97.68 75.55 137.01 97.68 66.07 43.53
6113.51 157.69 144.24 104.29 48.07 157.23 151.40 136.26 119.88 156.55 136.26 111.93 90.08
870.70 153.20 152.92 122.61 70.65 147.87 156.26 147.79 135.45 155.89 147.79 129.00 110.34
404489.22 454.83 353.44 186.17 43.89 471.93 394.40 314.17 243.96 435.49 314.17 213.55 140.48
6443.68 571.38 540.94 404.51 189.74 566.48 560.66 516.22 460.70 571.33 516.22 432.36 351.58
8282.82 529.11 561.79 475.13 281.79 506.15 561.53 552.06 517.46 545.35 552.06 496.70 431.74
Note: α = 0.4   m 1 , θ s θ r = 0.4 , and m = 0 .
Table 9. Additional bearing capacity q u s u c t i o n due to transient infiltration ( q 0 / k s = 0.75 ).
Table 9. Additional bearing capacity q u s u c t i o n due to transient infiltration ( q 0 / k s = 0.75 ).
φ ( ) l ( m ) k s = 5 × 10 6 k s = 5 × 10 7 k s = 5 × 10 8
t = 0 dt = 0.5 dt = 1 dt = 2 dt = 4 dt = 4 dt = 8 dt = 12 dt = 16 dt = 60 dt = 120 dt = 180 dt = 240 d
10412.67 12.89 11.36 8.93 6.85 13.17 11.97 10.78 9.76 12.59 10.78 9.32 8.27
68.95 13.48 13.13 11.59 9.26 13.36 13.37 12.84 12.21 13.50 12.84 11.90 11.01
85.48 13.09 13.46 12.49 10.51 12.72 13.48 13.33 12.94 13.31 13.33 12.72 12.04
20447.02 47.43 41.83 32.68 24.78 48.41 44.09 39.67 35.83 46.35 39.67 34.17 30.18
633.66 49.23 48.42 42.88 34.11 48.62 49.18 47.42 45.16 49.44 47.42 44.01 40.73
820.72 47.14 49.37 46.23 38.88 45.42 49.22 49.06 47.80 48.21 49.06 47.03 44.59
304153.96 153.12 134.96 104.24 77.32 156.16 142.39 127.80 114.90 149.71 127.80 109.29 95.76
6113.51 157.52 157.11 140.06 110.64 154.98 158.91 154.31 147.40 158.79 154.31 143.73 132.98
870.70 146.81 158.54 151.03 127.22 140.14 156.68 158.52 155.60 151.47 158.52 153.43 145.97
404489.22 473.14 409.64 305.36 217.45 482.97 436.21 384.38 340.81 461.71 384.38 322.08 277.35
6443.68 568.39 574.45 517.48 403.89 558.79 577.96 566.52 543.77 574.27 566.52 530.80 491.00
8282.82 507.11 567.02 558.46 473.38 482.65 552.76 573.08 570.69 526.40 573.08 565.34 542.16
Note: α = 0.4   m 1 , θ s θ r = 0.4 , and m = 0 .
Table 10. Additional bearing capacity q u s u c t i o n due to transient infiltration ( q 0 / k s = 0.5 ).
Table 10. Additional bearing capacity q u s u c t i o n due to transient infiltration ( q 0 / k s = 0.5 ).
φ ( ) l ( m ) k s = 5 × 10 6 k s = 5 × 10 7 k s = 5 × 10 8
t = 0 dt = 0.5 dt = 1 dt = 2 dt = 4 dt = 4 dt = 8 dt = 12 dt = 16 dt = 60 dt = 120 dt = 180 dt = 240 d
10412.67 13.47 12.97 12.00 11.10 13.53 13.19 12.75 12.34 13.39 12.75 12.16 11.72
68.95 13.00 13.51 13.39 12.74 12.74 13.40 13.54 13.50 13.19 13.54 13.45 13.25
85.48 11.97 13.11 13.54 13.28 11.48 12.81 13.30 13.48 12.33 13.30 13.52 13.52
20447.02 49.36 47.44 43.66 40.15 49.59 48.29 46.60 45.01 49.05 46.60 44.30 42.56
633.66 47.39 49.47 49.14 46.67 46.36 49.03 49.63 49.53 48.13 49.63 49.35 48.64
820.72 42.93 47.69 49.61 48.74 40.98 46.43 48.49 49.33 44.40 48.49 49.52 49.62
304153.96 158.57 151.63 138.10 125.56 159.45 154.65 148.60 142.92 157.41 148.60 140.39 134.19
6113.51 151.69 159.17 158.54 150.08 148.25 157.48 159.90 159.76 154.23 159.90 159.23 156.88
870.70 133.80 151.46 159.64 157.40 127.15 146.57 154.69 158.23 139.05 154.69 159.12 159.94
404489.22 485.76 453.65 402.66 359.19 490.38 467.36 441.49 420.22 480.78 441.49 410.98 389.00
6443.68 550.68 577.74 575.96 540.72 539.38 571.20 580.80 580.64 559.41 580.80 578.67 569.24
8282.82 467.50 536.78 577.48 573.26 444.77 515.95 551.54 569.23 486.53 551.54 574.21 580.52
Note: α = 0.4   m 1 , θ s θ r = 0.4 , and m = 0 .
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, S.; Zhou, D. An Analytical Framework for Assessing the Unsaturated Bearing Capacity of Strip Footings under Transient Infiltration. Mathematics 2023, 11, 3480. https://doi.org/10.3390/math11163480

AMA Style

Xu S, Zhou D. An Analytical Framework for Assessing the Unsaturated Bearing Capacity of Strip Footings under Transient Infiltration. Mathematics. 2023; 11(16):3480. https://doi.org/10.3390/math11163480

Chicago/Turabian Style

Xu, Sheng, and De Zhou. 2023. "An Analytical Framework for Assessing the Unsaturated Bearing Capacity of Strip Footings under Transient Infiltration" Mathematics 11, no. 16: 3480. https://doi.org/10.3390/math11163480

APA Style

Xu, S., & Zhou, D. (2023). An Analytical Framework for Assessing the Unsaturated Bearing Capacity of Strip Footings under Transient Infiltration. Mathematics, 11(16), 3480. https://doi.org/10.3390/math11163480

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