Abstract
The bearing capacity of unsaturated soft soil foundations is affected by natural or artificial hard crust, due to the stress diffusion and closure effect. Based on the triple-shear unified solution of shear strength for unsaturated soil, the analytical solutions of critical edge load and critical load of unsaturated soft soil with hard crust are derived with the consideration of the lateral pressure coefficient . The formulas are degraded and compared. Meanwhile, the influences of stress diffusion, increased load, matric suction, intermediate principal, and lateral pressure coefficient are conducted on the bearing capacity of unsaturated soft soil foundations with hard crust. The results show that the critical edge load and critical load increased, and the stress diffusion coefficient decreased with the increase of the ratio of elastic modulus between hard crust and soft soil. The bearing capacity of soft soil foundation increased with the consideration of the stress diffusion and closure effect of hard crust. The effect of matric suction is remarkable, and the critical edge load and critical load linearly increase with the increase of low matric suction and variate with the dual influence of high matric suction. The bearing capacity of unsaturated soft soil foundation is obviously higher than that of saturated soil. The critical edge load and critical load increased with the increase of intermediate principal stress and lateral pressure coefficient. The triple-shear unified solution for bearing capacity of unsaturated soft soil foundation with hard crust contains different solutions under various working conditions, which has important theoretical significance and widely engineering application.
1. Introduction
Soft soil is not saturated due to evaporation and transpiration of water molecules, and this phenomenon is particularly prevalent in arid and semiarid regions [1–3]. With the long-term effect of dehydration, drying, weathering, ion exchange, and other factors or the influence of human action, hard crust with different thickness is formed on the surface of soft soil [4]. The bearing capacity of soft soil foundation is effected by hard crust, due to the stress diffusion and closure effect [5]. Wen et al. [5, 6] deeply analyzed the influence of stress diffusion and closure effect of hard crust on soft soil foundation, the results show that the stress diffusion expands the load area on the top of soft soil foundation, and the maximum value decreases. The closure effect delays the load of soft soil into the plastic stage and improves the bearing capacity of the foundation. Chen et al. [7, 8] indicated that the in situ stress state of natural soil is usually inherently anisotropic, that is to say, the static lateral pressure coefficient in the actual stress state of soil. Fattah et al. [9–12] studied the shear strength properties for unsaturated soil with different matric suction values, and it is concluded that the matric suction is considered an important contributor to unsaturated shear strength properties. Wang [13] and Wen et al. [5] derived the critical edge load of soft soil foundation with hard crust based on Mohr-Coulomb criterion, the results take no account of intermediate principal stress, and the assumption is inconsistent with the actual stress state of soil. Zhu et al. [14] deduced the critical edge pressure and critical loading of soft soil foundation based on the unified strength theory from the static lateral pressure coefficient of soil under actual conditions, and the results consider the stress diffusion effect and the influence of intermediate principal stress but do not consider the actual unsaturated state of soft soil. Zhang et al. [15] derived the critical edge load and critical load of unsaturated soft soil foundation based on the unified solution of plane shear strength of unsaturated soil, taking account of the stress dispersion and self-weight berm of hard crust. However, the results is derived based on the twin shear unified strength theory, and the problem of double slip angle exists under certain stress [16]. Gao et al. [17] obtained the calculation formula of plastic load and critical load of soft soil foundation with hard crust based on the triple-shear unified strength criterion, considering the influence of static lateral pressure coefficient, stress diffusion, and closure effect of hard crust but do not consider the actual unsaturated state of soft soil.
The above results show that considering the intermediate principal stress effect can better develop the strength of soil and improve the bearing capacity of soft soil foundation. The static lateral pressure coefficient has significant influence on the critical edge load and critical load of soft soil foundation, and is more consistent with the actual situation of soil. Matric suction of unsaturated soil has different effects on bearing capacity of soft soil foundation under high and low suction. Therefore, the effects of stress diffusion, closure effect, intermediate principal stress, static lateral pressure coefficient, and matric suction of unsaturated soil should be comprehensively considered.
Based on the triple-shear unified solution of shear strength for unsaturated soil, with the consideration of stress diffusion and closure effect of hard crust, the analytical solution of the critical edge load and critical load of unsaturated soft soil with hard crust are deduced under the lateral pressure coefficient . The formulas are degraded and compared. Meanwhile, the influences of stress diffusion , increased load , matric suction , strength criterion parameter , and lateral pressure coefficient are conducted on the bearing capacity of unsaturated soft soil foundation with hard crust. The results can provide theoretical reference for the design and construction of soft soil foundation with hard crust and have important engineering application value.
2. Model Establishment and Verification
2.1. The Triple-Shear Unified Solution of Shear Strength for Unsaturated Soil
The triple-shear unified strength criterion [18] overcomes the problem of double slip angle under certain stress in twin-shear unified strength criterion, which considers the influence of intermediate principal stress effect. Based on the triple-shear unified strength criterion and shear strength of unsaturated soil in terms of two state stress variables, a triple-shear unified solution of shear strength for unsaturated soil is obtained with comprehensively considering the effects of intermediate principal stress, matric suction, strength disparity of material, and strength criterion [19]. The formula is as follows:
where
In the formula, is the shear strength of unsaturated soil, is the unified total cohesion of unsaturated soil, is the triple-shear unified effective cohesion, is the triple-shear unified effective angle of internal friction, is the triple-shear unified angle related to matric suction, is the total normal stress, is the pore gas pressure, is the pore water pressure, is the net normal stress, is the matric suction, is the effective cohesion intercept, is the effective angle of internal friction, is the angle of matric suction contribution to shear strength; is the lode parameter, is the parameter of different strength criteria, which reflect the influence of intermediate principal shear stress and normal stress on material failure, and the value of is related to the mechanical properties of the material itself, which is determined by the test, where is the tensile yield limit, is the compression yield limit, and is the shear yield strength.
Equation (1) take intermediate principal stress effect into account by strength criterion parameter and lode parameter to reflect the influence of strength disparity of material and strength criterion difference on its shear strength. Meanwhile, effect of matric suction on shear strength of unsaturated soil also takes into account by . When , equation (1) degenerates into shear strength of unsaturated soil based on Mohr-Coulomb strength criterion [20], without considering the influence of intermediate principal stress. When , it degenerates into shear strength of saturated soil based on the triple-shear strength criterion [21].
2.2. Stress Analysis
Ignoring the compression and shear deformation of hard crust, the stress of any point in the foundation is shown in Figure 1. In the figure, is the uniform strip load on foundation, is the load width, is the thickness of hard crust, , , and are gravity, elastic modulus, and Poisson ratio of the hard crust, respectively, , , and are gravity, elastic modulus, and Poisson ratio of the soft layer, respectively, and are effective cohesion intercept and effective angle of internal friction of the soft soil, and is the suction angle related to matric suction of soft soil.

The hard crust has stress diffusion and closure effect on the underlying soft soil layer. Stress diffusion is generalized by stress diffusion coefficient ; so, the uniform strip load transmits to the top of soft soil and becomes [22]. Closure effect can be simplified by calculating that the increments of horizontal stress and vertical stress are the same, and both are the load imposed on the top of hard crust in the first stage [5]. Therefore, the stress at any point in the foundation consists of additional stress due to stress diffusion, self-weight stress, and stress generated by considering closure effect.
2.2.1. Additional Stress considering Stress Diffusion
The stress diffusion of hard crust to the underlying soft soil is related to the thickness of hard crust, the width of the load, and the compression modulus of the hard crust and soft soil. Xu [23] obtained the stress diffusion coefficient by fitting the finite element calculation, and its expression is
where
In the formula, and are intermediate parameters of elastic modulus, and and are fitting parameters.
Considering the stress diffusion, the uniform strip load on the top of soft soil is . Then, the additional stress of any point in soft soil foundation is [24]
2.2.2. Self-Weight Stress
The self-weight stress of any point in soft soil is
In the formula, is the static lateral pressure coefficient of soft soil, and generally, is less than one.
2.2.3. Additional Stress considering Closure Effect
The closure effect of hard crust has the same increment of the horizontal stress and the vertical stress, and both are the load imposed on the top of hard crust in the first stage. Then, the additional stress at any point in the soft soil is
2.2.4. The Total Stress at Point of Foundation
According to formulas ((5)–(11)), the total stress at any point in the soft soil foundation is
2.2.5. The Principal Stress at Point of Foundation
The principal stress at point is obtained from formulas ((12)–(14)).
The stress in the plastic zone is established within the depth range of /4 under the foundation edge [25], which meets the application conditions of the mathematical approximation formula [26]. So, equation (15) can be simplified as
2.3. Triple-Shear Unified Solution of Foundation Bearing Capacity
According to the triple-shear unified strength criterion, when point is in the limit equilibrium state, the major and minor principal stresses satisfy the following equilibrium conditions:
Substituting equation ((16) and (17)) into equation (18) can be obtained
According to the approximate formula , the upper formula
Let , and substituting equation (20) into equation (19) can be obtained development depth of plastic zone
When the properties of hard crust and soft soil are certain, is just a function of ; so, the maximum development depth of plastic zone can be obtained by derivation of by equation (21), namely,
So, , and substitute into equation (21) to obtain
When the maximum development depth of the plastic zone , the critical edge load of unsaturated soft soil foundation with hard crust can be obtained by equation (23).
When the maximum development depth of the plastic zone , the critical load of unsaturated soft soil foundation with hard crust can be obtained by equation (23).
In the formula, , , , and are the coefficients of foundation bearing capacity, and their expressions are, respectively,
From formula (24) and formula (25), it can be seen that the critical edge load of unsaturated soft soil foundation with hard crust is composed of three components: one is related to the weight of hard crust , one is related to the unified total cohesion of soft soil , and the other is related to the increasing load with the consideration of closure effect of hard crust; the critical load increases one more related to the weight of hard crust . From equation (26), it can be seen that the weight coefficients of critical edge load and critical load are related to the stress diffusion coefficient , the static lateral pressure coefficient , and the triple-shear unified effective angle of internal friction , which are the functions of strength criterion parameters and lode parameter.
Therefore, the triple-shear unified solution formula (24) and formula (25) of the bearing capacity of unsaturated soft soil foundation with hard crust takes into account the combined effects of the stress diffusion (stress diffusion coefficient ), the closure effect (increasing load ), the characteristics of unsaturated soil (matric suction ), the intermediate principal stress effect (strength criterion parameters and lode parameter), and the actual stress state (static lateral pressure coefficient ). It has important theoretical significance and broad engineering application prospects.
3. Analysis
3.1. Formula Degradation and Comparative Analysis
3.1.1. Formula Degradation Analysis
When and , equation (24) and equation (25) degenerate into ignoring the stress diffusion and closure effect, and the expression is
Formula ((27), (28)) is basically the same as the expression in literature [27], and the main difference is that the load on the soft soil in literature [27] is . Therefore, is added to the expressions of critical edge load and critical load, and the linear distribution of matric suction is also added, while other terms are the same.
When , , , and , it degenerates into the Mohr-Coulomb solution of saturated soft soil foundation without considering the stress diffusion and closure effect, and the expression is
Formula ((29), (30)) is basically the same as the expression in literature [28], and the main difference is that the load on the soil in literature [28] is . So, is added to the expressions of critical edge load and critical load. Meanwhile, the literature is single-layer soil, while this paper is double-layer soil (hard crust and soft soil), and other items are the same.
When , , and , it degenerates into the Mohr-Coulomb solution of saturated soft soil foundation with hard crust in the hydrostatic self-weight stress field and is expressed as
Formula (31) is the same as the expression in literature [5], which considering the stress diffusion and closure effect, but not considering the matric suction of unsaturated soil, intermediate principal stress effect, and .
When , , , , , and , it degraded to
Equation ((33), (34)) is the same as the formula of critical edge load and critical load of homogeneous foundation with zero burial depth in soil mechanics; so, the critical edge load and critical load of homogeneous foundation in classical soil mechanics are special case of this paper.
From the above analysis, it can be seen that the analytical solutions equation (24) and equation (25) obtained in this paper contain different solutions under various working conditions, which can be reasonably selected according to the specific conditions in practical engineering and have wide applicability.
3.1.2. Comparative Analysis
Based on the unified solution of plane strain shear strength of unsaturated soil, Zhang et al. [15] established the critical edge load and critical load of unsaturated soft soil foundation with hard crust, taking account of the stress dispersion and self-weight berm of hard crust. The analysis idea and derivation process of this paper are similar, and the obtained expression is basically the same. Comparing literature [15] with equation (24) and equation (25) in this paper, it can be seen that there are three differences in the expressions of critical edge load and critical load: first, this paper believes that the increment value of horizontal stress and vertical stress by closure effect is the same, and both are the loads imposed on the top surface of hard crust in the first stage; so, a term related to the increase load of the closure effect is added to the expression. Second, when calculated the total stress of in literature [15], the normal stress and shear stress of self-weight stress are obtained in the direction of additional stresses and and then superimposed, but the direction of total stress is not the direction of principal stress. In this paper, the stress in and directions are calculated, respectively, and then superimposed. At this time, the stress direction is consistent, and the principal stress should be sought after superposition. Third, the meaning of and is different. In literature [15], unified effective cohesion and effective angle of internal friction are analogy obtained based on the twin shear unified strength theory, while in this paper, triple-shear unified effective cohesion and effective angle of internal friction are obtained based on the triple-shear unified strength criterion. The twin shear unified strength criterion exist double slip angle under some specific stress, the analytical solution of this paper adopts the triple-shear unified strength criterion, which overcomes the problem of double slip angle.
3.2. Analysis of Influencing Factors
In this paper, the critical edge load and critical load of unsaturated soft soil foundation with hard crust are established based on the triple-shear unified solution of shear strength for unsaturated soil, taking into account the effects of the stress diffusion, the closure effect, the characteristics of unsaturated soil, the intermediate principal stress effect, and the actual stress state. Therefore, the influence of the above five factors is analyzed in this section.
The basic parameters of unsaturated soft soil foundation with hard crust are [7, 14] as follows: load width , thickness of hard crust , gravity , elastic modulus , and Poisson ratio . Here are the following parameters: the gravity of soft soil , elastic modulus , Poisson ratio , effective cohesion intercept , effective angle of internal friction , and under plane strain state [19].
3.2.1. The Stress Diffusion (Stress Diffusion Coefficient )
It can be seen from equation (3) that the ratio of elastic modulus between hard crust and soft soil directly affects the stress diffusion, namely, the stress diffusion coefficient , and then has a certain influence on the critical edge load and critical load . Figures 2 and Figures 3, respectively, show the variation rules of , , and with and , with when , , , and .

(a) Stress diffusion coefficient

(b) The critical edge load and critical load

It can be seen from Figures 2 and 3 that with the increase of the ratio of elastic modulus between hard crust and soft soil , the stress diffusion coefficient gradually decreases, and the critical edge load and critical load gradually increase. When increases from 3.0 to 7.0, decreases by 9.9%, and and increase by 11.0%.
Recommend non-segmenting, when the elastic modulus of hard crust and soft soil are quite different, the stress diffusion effect of hard crust on soft soil increases and then leads to the increase of the bearing capacity of soft soil foundation. increased by 8.52 kPa and increased by 12.03 kPa when decreased from 0.90 to 0.81, both of which increased by 11%. In other words, the diffusion effect of hard crust on soft soil foundation can improve the bearing capacity of soft soil foundation, but the effect is not significant. Therefore, reasonable natural or artificial hard crust should be selected to treat soft soil foundation in practical engineering, and the stress diffusion effect of hard crust on soft soil should be fully utilized to increase the bearing capacity of soft soil foundation.
3.2.2. The Closure Effect (Increasing Load )
According to equations (24) and (25), it can be seen that the influence of closure effect on the critical edge load and critical load is reflected in . Figure 4 shows the variation rules of and with when , , , and .

It can be seen from Figure 4 that the critical edge load and critical load increase linearly with the increase of (load applied on the top of hard crust in the first stage). When increases by 1 kPa, the critical edge load and critical load increase 0.66 kPa. When increases from 0 kPa to 10 kPa, the critical edge load and critical load increase 7.91% and 5.60%, respectively. Therefore, considering the closure effect of hard crust, the bearing capacity of soft soil foundation increases, and the increase depends on the load applied on the top of hard crust in the first stage.
3.2.3. The Characteristics of Unsaturated Soil (Matric Suction )
From equation (24) and (25), it can be seen that the critical edge load and critical load are related to the unified total cohesion of soft soil, which is the function of effective cohesion and matric suction, and matric suction directly affects the bearing capacity of soft soil foundation. Figure 5 shows the variation rules of critical edge load and critical load with matric suction when , , , , and the suction angle of matric suction adopts the hyperbolic model in literature [15].

(a) The critical edge load

(b) The critical load
It can be seen from Figure 5 that the variation of critical edge load and critical load with matric suction is related to , which is the intercept of the hyperbolic model of suction angle under high suction, reflecting the change rate of suction angle. When is in the range of low suction, and are independent of the intercept and increase linearly with the increase of . When is the air-entry value of unsaturated soft soil, and reach the peak point () or the turning point (), and the larger the intercept , the less obvious the turning trend. When is in the range of high suction, and are closely related to the intercept . When , and increase linearly with the increase of . With the increase of , it gradually turns to nonlinear increase. This curve reflects the nonlinearity of unsaturated soil strength in the range of high matric suction, namely, the double influence of matric suction on unsaturated soil strength. When , the decrease of suction angle makes the effect of matric suction on strength reduction more obvious than that of matric suction increase, and the transition is obvious. With the increase of , the improvement effect of matric suction gradually increases, which exceeds the influence of the decrease of suction angle, occupies a dominant position, and the turning trend is less obvious. When , unsaturated soft soil becomes saturated soft soil, and are less than the bearing capacity of unsaturated soft soil foundation. Therefore, the double influence of matric suction should be considered in the design and construction of soft soil foundation and fully utilize the bearing capacity of unsaturated soft soil significantly higher than that of saturated soil.
3.2.4. The Intermediate Principal Stress Effect (Strength Criterion Parameters )
Under plane strain state , parameter reflects the effect of intermediate principal stress, which is the selection parameter of different strength criteria. Figure 6 shows the variation rules of critical edge load and critical load with strength criterion parameters when , , , and .

It can be seen from Figure 6 that with the increase of , and increase nonlinearly. When increases from 0 to 1, increase from 83.05 kPa to 120.69 kPa, with an increase of 50.77%, and increase from 117.36 kPa to 170.53 kPa, with an increase of 45.31%. In other words, considering the intermediate principal stress effect can give full play to the strength potential of soft soil and improve its bearing capacity. Therefore, the effect of intermediate principal stress should be reasonably considered in the design of soft soil foundation, which can give full play to the self-bearing capacity of materials, save materials, and reduce costs.
3.2.5. The Stress State of Foundation Soil (Static Lateral Pressure Coefficient )
The static lateral pressure coefficient reflects the initial stress state of foundation soil. The general soil mechanics textbook assumes that ; that is, the self-weight stress of the soil is a hydrostatic stress state, which artificially increases the lateral stress of the foundation soil and overestimates the bearing capacity of the foundation. But in practice, most of the soil is . Figure 7 shows the variation rules of critical edge load and critical load with static lateral pressure coefficient when , , , and .

It can be seen from Figure 7 that with the increase of , and increase significantly. When ( reaches its minimum when while ), the difference between critical edge load and critical load is small, but the difference gradually becomes obvious with the increase of . When increases from 0.76 to 1.0, the critical edge load increases 75.18%, and the critical load increases 144%. That is to say, the artificial setting of significantly improves the bearing capacity of the foundation, which is seriously inconsistent with the actual. Therefore, the real static lateral pressure coefficient of the foundation should be reasonably measured in the design of soft soil foundation.
4. Conclusions
(1)The triple-shear unified solution of the foundation bearing capacity of unsaturated soft soil with hard crust is established in this paper, considering the comprehensive effects of the stress diffusion (stress diffusion coefficient ), the closure effect (increasing load ), the characteristics of unsaturated soil (matric suction ), the intermediate principal stress effect (strength criterion parameters ), and the actual stress state (static lateral pressure coefficient ). It contains different solutions under various working conditions, which is widely applicable to more engineering situations and has important theoretical significance and broad engineering application prospects(2)The critical edge load and critical load increase with the increase of the ratio of elastic modulus between hard crust and soft soil and increase with the decrease of the stress diffusion coefficient, while linearly increase with the increase of (the load applied on the top of the first stage of the hard crust). Therefore, the foundation bearing capacity of unsaturated soft soil is increased by considering the stress diffusion and closure effect of hard crust(3)The influence of matric suction on critical edge load and critical load is related to the values of matric suction. and increase linearly with the increase of matric suction in the range of low suction and reach the peak point or the turning point when the matric suction equal to the air-entry value of unsaturated soft soil, while changes with the relative magnitude of the double influence due to the matric suction in range of the high suction(4)The critical edge load and critical load increase nonlinearly with the increase of strength theoretical parameter , and increase with the increase of static lateral pressure coefficient . The calculation result of the foundation bearing capacity of unsaturated soft soil by Mohr-Coulomb criterion is too conservative, and the calculation by is seriously inconsistent with the actual. Therefore, the intermediate principal stress effect should be reasonably considered, and the real static lateral pressure coefficient of foundation soil should be measured
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study is supported by the National Natural Science Foundation of China (No. 12102379), Scientific Research Project of Education Department of Shaanxi Province (No. 21JK0977, No. 21JC035), and Doctoral Research Start-Up Foundation of Yan’an University (No. YDBK2018-01).