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
and apparent cohesion
distributions of fine sand, silt, and clay at different times
t and with different water table depths
l. It is assumed that
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
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
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
compared to the middle soil. As a result, the upper and lower parts of the soil layer exhibit a higher saturation degree
, while the middle part remains relatively low in saturation degree
.
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
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
generally decreases over time. Moreover, the trend in
along the depth is similar to that of
. When the water table drops from 3 m to 6 m, the clay’s
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
initially increase and then decrease with increasing depth. However, when the water table drops to 6 m, the profile of
becomes disordered. These differences in cohesion variations among different soils are primarily attributed to variations in the values of
and
, which result in different saturation degrees
.
The apparent cohesion
can be obtained by substituting the saturation into Equation (20). Differentiating Equation (20) with respect to
K yields:
When the saturation exceeds 1/e, the apparent cohesion
decreases as the saturation degree
. increases. Conversely, when
is less than 1/e, the opposite trend occurs, as shown in
Figure 13. Therefore, the maximum value of
occurs when
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
caused by transient seepage is an independent item, unrelated to the effective cohesion
, surcharge load
, and unit weight of soil
. The dimensionless bearing capacity factors
,
, and
, 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
. By analyzing Equation (39), it can be determined that the parameters influencing
are the permeability coefficient (
), infiltration intensity (
), water storage capacity (
), 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.,
.
Figure 14,
Figure 15 and
Figure 16 depict the results of the additional bearing capacity
for the three soil types, involving different infiltration times
t, water table depths
l, and infiltration ratios
. From the graphs, it can be observed that at
t = 0, the value of
due to matric suction in unsaturated soils does not vary with changes in
. Among these three soil types, clay exhibits the highest value of
, 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
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
becomes disordered for fine sand and silt. However, the value of
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
on the additional bearing capacity
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.,
. Previous studies [
2,
30,
31] have typically assumed that
to simplify the analysis. In practical engineering, due to uncertainties in rainfall and surface water, the infiltration ratio
can take any value between 0 and 1. Therefore, this study considers six scenarios by extending the investigation to different infiltration ratios of
= 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
, the
for fine sand and silt show a similar trend over time, while the
for clay monotonically decreases over time. Specifically, for fine sand and silt, when
, their
value initially increases and then decreases as time goes by. Moreover, smaller values of
result in larger maximum values of
and longer corresponding infiltration times
t. The maximum value of
at
is the highest among all scenarios, with magnitudes of 91.15 kPa and 70.44 kPa, respectively. The maximum values of
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
generates a slower rate of decrease. When
, the
of fine sand and silt monotonically increase with time, but the rate of increase gradually slows down. Furthermore, a larger
leads to a greater magnitude of increase in
. The response of clay’s
to the development of transient seepage is completely different. As time goes by, the
corresponding to all values of
decreases, and the range of reduction increases with an increase in
.
The reason for the different trends corresponding to different infiltration ratios 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 . On the other hand, clay has a saturation degree greater than 1/e initially, so its apparent cohesion will decrease regardless of the infiltration flux.
Assuming
and
,
Figure 17,
Figure 18 and
Figure 19 depict the variation in additional bearing capacity
with respect to the effective internal friction angle
for different values of the water table depth
l, water storage capacity
, and desaturation coefficient
. As the water table decreases,
slightly increases. For smaller values of
, an increase in
also leads to an increase in
, 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
continuously increases, and the speed of the increase accelerates. The influence of
on
becomes more evident for smaller values of
. For example, when
,
, and
l = 8 m, the value of
increases from 27.05 kPa to 1161.52 kPa with increasing
from 10° to 40°. From the equation
, it can be determined that an increase in
can directly increase the
, thereby enhancing
. 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
to the desaturation coefficient
.