Next Article in Journal
Conjunctive Water Management under Multiple Uncertainties: A Case Study of the Amu Darya River Basin, Central Asia
Next Article in Special Issue
Influence of Fines Content on the Stability of Volcanic Embankments under Rainfall and Earthquake
Previous Article in Journal
Health-Aware Economic MPC for Operational Management of Flow-Based Networks Using Bayesian Networks
Previous Article in Special Issue
Function of a Deep-Buried Isolated Trench and Its Effect on Cracking Failure Characteristics of a Slope under Artificial Rainfall
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Algorithm Implementation of Equivalent Expansive Force in Strength Reduction FEM and Its Application in the Stability of Expansive Soil Slope

1
Institute of Geotechnical Engineering, Xi’an University of Technology, Xi’an 710048, China
2
Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100081, China
3
School of Management, Xi’an University of Architecture & Technology, Xi’an 710055, China
4
Xi’an Jianchuang Geotechnique Technology Co., Ltd., Xi’an 710075, China
5
Huangshi City Construction Market Management Station, Huangshi 435000, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(10), 1540; https://doi.org/10.3390/w14101540
Submission received: 8 April 2022 / Revised: 6 May 2022 / Accepted: 9 May 2022 / Published: 11 May 2022
(This article belongs to the Special Issue Rainfall-Induced Geological Disasters)

Abstract

:
Loss of matric suction during rain has an important effect on the instability of expansive soil slope. A new device was designed for testing the expansive force in order to propose and determine the equivalent expansive force corresponding to the matric suction. First, the internal relation between the matric suction and the corresponding equivalent expansive force of unsaturated expansive soil was analyzed. Then, numerical algorithm implementation of the equivalent expansive force was discussed, and the equivalent expansive force was introduced into the strength reduction finite element method (FEM). Finally, the equivalent effects of the matric suction and the equivalent expansive force were compared and analyzed by evaluating the stability of an unsaturated expansive soil slope. The results show that the new testing device significantly improves the accuracy of the expansive force test and shortens the testing time. The relation between the matric suction and the equivalent expansive force with the change in initial water content is obvious. The equivalent expansive force can reflect the macro contribution of matric suction to unsaturated expansive soil, and the developed strength reduction FEM based on the equivalent expansive force can be used to evaluate the rainfall-induced instability of an expansive soil slope caused by the decrease in matric suction resulting from the rainfall infiltration.

1. Introduction

Rainfall-induced landslides often occur in expansive soil regions and always lead to regional geological disasters [1]. Expansive soil is a special type of soil with an obvious swell–shrink characteristic, which is distributed extensively [2,3]. The rainfall-induced instability of unsaturated expansive soil slopes is a problem that urgently needs to be solved in geotechnical engineering [4,5]. For one thing, rainfall-induced landslide hazard and risk analyses should be evaluated from a macroscopic perspective at the regional scale, and some probabilistic methods are used to assess regional susceptibility to landslides. For another, the triggering mechanism and kinematic process of landslides caused by rainfall infiltration should be researched [6].
Plenty of natural expansive soil landslides are often accompanied by rainfall infiltration, excavation, or earthquake action. Under the condition of rainfall infiltration, the water content of an expansive soil slope increases and the matric suction decreases [7], which indicates that the safety factor of an expansive soil slope also decreases gradually [8,9]. It is noticed that the measurement of matric suction is time-consuming, and the technique is complex and intricate [10,11,12]. For this reason, it is necessary to propose a simple, reasonable, and relatively advanced test method that can not only reflect the characteristics of unsaturated expansive soil but can also avoid the difficulties in the measurement of matric suction. In the view of Khan et al. [13], matric suction has a relation with the swell–shrink characteristic of unsaturated expansive soil. The expansive force can be used to effectively avoid the trouble in the measurement of matric suction.
At the same time, LU Zhaojun’s unsaturated expansive soil strength [14,15] disclosed the correlation between the matric suction and the expansive force. Therefore, the computational method for research on the rainfall-induced instability of unsaturated expansive soil, based on the equivalent substitute of matric suction by the expansive force, is also worthy of attention.
Owing to landslides of unsaturated soil slopes often occur during rainfall infiltration, the mechanics of unsaturated soil is also worthy of research. It is important to note that a reasonable evaluation of landslides should be based on the reasonable strength theory of unsaturated soil. Bishop et al. proposed the effective stress strength formula for unsaturated soil [16]. Based on the effective stress strength formula for unsaturated soil proposed by Bishop et al., Fredlund et al. suggested the double stress strength formula for unsaturated soil [17].
Through Bishop’s effective stress strength formula and Fredlund’s double stress strength formula, along with the analysis aspect of matric suction, Bishop’s slices method was applied to study the influence of matric suction on unsaturated expansive soil slope stability [18,19]. A plastic upper bound limit analysis was also conducted to analyze the influence of swelling deformation [20,21]. Yao et al. [22] showed that rainfall duration could also strongly affect the transient flows and reduce the stability of unsaturated expansive soil slopes. In the studies by Noorany Iraj [23] and Diyaljee et al. [24], the lateral extension in expansive soil slopes was also discussed. Based on a model test, Xiao et al. [25] evaluated and compared the reinforcement effect of geogrid in expansive soil slopes, and some failure mechanisms of expansive soil slopes were discussed from the aspects of strength characteristics [26,27]. These studies were advantageous to research on the stability principles of unsaturated expansive soil slopes. However, these methods only considered matric suction on sliding surfaces, whereas the actual influence of matric suction out of the sliding surface in unsaturated soil slope was underestimated, that is, these studies have not been involved in the non-uniform distribution of matric suction. However, it is worth noting that the strength reduction finite element method (FEM) based on Bishop’s effective stress strength formula and Fredlund’s double stress strength formula was developed and applied in the stability analysis of unsaturated slopes according to a non-uniform distribution of matric suction in the slopes [28].
Matric suction has a relation with the swell–shrink characteristic of unsaturated expansive soil [14,29], so the amount of swelling resulting from expansibility transforms into an expansive force when its bulk is restricted or its volume remains constant. Thus, the suctional strength contributed by matric suction should have a relationship with the equivalent expansive force (EEF), which is defined as the maximum expansive force for the expansive soil sample from a soaking state to a strictly saturated state with a constant volume under a certain initial water content. The inner correlation between the matric suction and the equivalent expansive force could provide a new analysis perspective to study the strength characteristic of unsaturated expansive soil. Lu et al. proposed a comprehensive formula for the strength of unsaturated expansive soil [15].
Aiming at the slope engineering of expansive soil, the strength reduction FEM based on Lu’s comprehensive formula, which will act as a new stability analysis method of unsaturated expansive soil slopes, will effectively widen the scope and thought for the stability analysis of unsaturated slopes.
In this study, the theoretical equivalent relationship between matric suction and the expansive force in the shear strength of unsaturated expansive soil is discussed. Then, some tests of matric suction and the expansive force as well as direct shear tests on expansive soil from Hanyin County in China are conducted. Studies on the correlation between matric suction and the equivalent expansive force with varying initial water contents of unsaturated expansive soil are also performed. The numerical algorithm implementation of the equivalent expansive force is discussed, and the equivalent expansive force is introduced into the strength reduction FEM. Equivalent effect research on evaluating the stability of unsaturated expansive soil slopes should be verified further based on the relation between matric suction and the equivalent expansive force for the sake of actual applications in engineering.

2. Equivalent Relation of Matric Suction and Expansive Force

This study focused on analyzing the theoretical equivalent relationship between matric suction and the expansive force in the shear strength of unsaturated expansive soil. Based on the deductive analysis of the relation between matric suction and the expansion force, the internal equivalent relation between them in the strength of unsaturated expansive soil can be reflected, and it is the theoretical basis of the core algorithm analysis for the stability analysis of expansive soil slopes.
Bishop’s effective stress strength formula for unsaturated soil is as follows,
τ f = c + ( σ u a ) tan φ + χ ( u a u w ) tan φ
where the shear strength of unsaturated soil is τf, the effective cohesion is c′, the angle of frictional resistance is φ’, the normal stress is σ, the pore water pressure is uw, the pore air pressure is ua, and the parameter of χ is determined by the saturation degree, Sr.
Fredlund’s double stress strength formula for unsaturated soil is as follows,
τ f = c + ( σ u a ) tan φ + ( u a u w ) tan φ b
where φ b is the suctional internal friction angle.
In further discussion under the condition of equal effective cohesion and equal effective stress, it can be seen that, in view of Formulas (1) and (2), Fredlund et al. proposed the equivalent relation of matric suction to the unsaturated soil strength as follows [30]:
χ ( u a u w ) tan φ = ( u a u w ) tan φ b
The simplified expression is Formula (4):
χ = tan φ b tan φ
There is an inner correlation between the matric suction and the equivalent expansive force, as reported by Lu’s comprehensive formula for the strength of unsaturated expansive soil.
τ f = c + ( σ u a ) tan φ + m p s tan φ
where m was the effective coefficient of the expansive force. ps was the expansive force of the expansive soil under a constant volume during saturation.
For Formulas (1), (2) and (5), under the condition of equal effective cohesion and equal effective stress, the equivalent relation of matric suction and expansion force to unsaturated soil strength are proposed in this paper as follows:
χ ( u a u w ) tan φ = ( u a u w ) tan φ b = m p s tan φ
Therefore, the simplified expressions are Formulas (7) and (8),
χ = m p s ( u a u w )
m = ( tan φ b tan φ ) / ( p s ( u a u w ) )
Through the above deductive analysis of the relation between matric suction and the expansion force, it can be concluded that there is a certain internal equivalent relation between them in the strength of unsaturated expansive soil, and this internal equivalent relation is the theoretical basis of the core algorithm analysis for the expansive soil slope in this paper.

3. Methodology

The testing methods for the traditional expansive force are numerous. Common traditional methods include the direct expansion-back pressure, preloading expansion, and gradually balancing pressure methods.
The direct expansion-back pressure method applies loads on a soil sample to return to its initial volume after fully absorbing water and freely expanding. The soil sample is generally stressed while expanding to its maximum limit during the experiment until it returns to its initial volume. To a certain extent, the force measured by the expansion-back pressure method should be the consolidation pressure rather than the true expansion pressure. When applying a huge pressure on the soil sample in the preloading expansion method, expansive soil is compressed, and its microstructure is disturbed and changed. Thus, the test value of the expansion force and the actual value often exhibit some differences. In the gradually balancing pressure method, loads should be applied continuously on the soil sample to return the pointer of the dial indicator to its initial position, i.e., during soaking and expanding in the consolidation apparatus. The test processes are complicated, and the values of the applied loads cannot be predicted. Moreover, complete saturation of the soil sample cannot be controlled under this method.
In order to consider the equivalent expansive force according to the requirements of physical significance, the microstructure of soil mass should be maintained in its original state during the test, and the corresponding measurements should be simple, accurate, and highly efficient. Based on these requirements, a new device that can control the saturation and constant volume in the test for expansive force was developed in the present work. The device had received a patent for invention [31] and could be adopted in research on the equivalent expansive force from a soaking state to complete saturation under a constant volume.
The structure of the testing device is presented in Figure 1. The rigid cutting ring system mainly consisted of a screw, a beam, a plunger, a porous disc, a cutting ring, and a cutting pedestal. The porous disc, cutting ring, and cutting pedestal were rigid. The positions of the beam and the cutting pedestal were fixed. The top cutting ring secured the position of the porous disc by fastening a screw to form a rigid cutting ring system with a fixed position. The main function of this system was to maintain the constant volume of the soil sample. Simultaneously, the screw, beam, plunger, and porous disc constituted the preloading system. The cutting ring and its surrounding parts as well as the rubber mat constituted the sealing system of the device. The vacuum saturation system was composed of the vacuum pump and the connecting pipe. Its main function was to control the saturation state of the soil sample in the cutting ring through the exhaust process of the vacuum pump. The measurement acquisition system consisted of an earth pressure sensor and an externally connected digital display. The former could directly measure the pressure on the earth pressure cell, and the data could be directly read by an externally connected digital display. The actual expansive force of the expanding soil could be calculated.
This study simulated the instability of filling an embankment with an expansive soil in the G316 National Highway. The G316 National Highway in Shaanxi Province was 603 km long. The length of expansive soil along this highway was approximately 200 km, whereas that in Hanyin County of Ankang was 64 km. Most of soil in the Hanjiang River basin and Ankang basin was basically typical expansive soil with a slope. The expansive soil sample for the test was collected from Hanyin County. Hanyin County belongs to the northern subtropical monsoon humid climate zone. The annual average temperature is 15.1 °C, the frost-free period is 258 days, precipitation is 782 mm, and the climate is mild and humid. For the equivalent expansive force test experiments under different initial water contents by the developed testing device, the change curves of the expansive force with time from a soaking state to a strictly saturated state are illustrated in Figure 2.
Further analysis of Figure 2, based on the subsequent experiment analyses on the unsaturated expansive soil, the relation between the equivalent expansive force, ps, and the initial water content, w, was found to follow a power function,
p s = A w λ
where A and λ were the fitting parameters according to the experimental data of the expansive soil.
The testing device had several advantages, including its simple structure, convenient operation, and highly efficient measurement process. It could solve problems on large disturbances in soil samples and the uncontrollable saturation in the direct expansion-back pressure, preloading expansion, and gradually balancing pressure methods.

4. Algorithm Implementation of Equivalent Expansion Force in Strength Reduction FEM

Effective strength parameters of the cohesion (c′) and angle of friction (φ′) are set to be the variables in the Mohr–Coulomb criterion in the strength reduction FEM. The reduced strength parameters c F and φ F can be defined as
c F = c / F
φ F = tan 1 ( tan φ / F )
The reduced strength parameters, respectively, replaced the strength parameters of c′ and φ′ of Mohr–Coulomb criterion. The strength reduction FEM was used to calculate the stress and strain of soil. First, the initial reduction coefficient, F, was set to be small so that the soil was in an elastic state. Then, the value of F was gradually increased until the destruction of the slope occurred, which showed that the strength reduction FEM calculation diverged under a physically real convergence criterion. The safety factor at failure was the value of F of the previous step.
In the strength reduction FEM, the shape of the sliding surface was not necessary to be assumed in this computation, so the potential sliding surface and the corresponding safety factor could be acquired by the calculation and its post-treatment. The strength reduction FEM was gradually becoming an indispensable tool in the stability analysis of slope. The details of the application of this method in the stability analysis of unsaturated expansive soil slopes is elaborated in this study based on the equivalent expansive force and Formula (3).
In fact, because of the uneven distribution of the matric suction, the corresponding distribution of the equivalent expansive force was also non-uniform, and according to Formula (3), the unsaturated strength of the expansive soil should also be non-uniform in matching with the distribution of the equivalent expansive force. By this token, the new non-uniform characteristic of effective stress could be produced by the non-uniform distribution of the equivalent expansive force whether the soil itself was homogeneous or not.
According to Formula (3), it could be found that the strength parameters of c′ and φ′ were not changed, and the Mohr–Coulomb criterion was employed as the strength failure surface in effective stress space. Because the matric suction was the effective stress, the equivalent expansive force that was the equivalent effect of the matric suction also had its significant contribution to the effective stress. In strength reduction FEM calculation, the macro effect of the equivalent expansive force could be obtained from the product of m and ps in every element, and then a new effective stress state could be formed when the macro effect mps of the equivalent expansive force was added to the gravitational effective stress in every element, which included the effects of the equivalent expansive force, as shown in Figure 3 within the three-dimensional effective stress space and Figure 4 within the meridian plane. Seen from that in the effective stress space, the stress state of the point A represented the initial or gravitational stress state in a saturated state, and mps represented the contribution of the equivalent expansive force to the increment of the effective stress. The stress state of the point B was a new effective stress state that integrated the gravitational effective stress with the mps in the unsaturated soil. The algorithm principle of the computational method on the stability of unsaturated expansive soil may indirectly reflect the contribution of matric suction using the corresponding equivalent expansive force, in which the matric suction is equivalently substituted by the equivalent expansive force.
In an analysis of strength reduction FEM, each element in the slope was assigned a value of the equivalent expansive force, so as to simulate its non-uniform distribution within the slope engineering. For the sake of effectively performing the implementation of the analysis on an unsaturated expansive soil slope in the strength reduction FEM, two distinctively different cases should be considered and treated separately:
(1) For those soil elements without the equivalent expansive force, as with the traditional strength reduction FEM, the effective stress strength parameters (c′ and φ′) were selected as the variables in the Mohr–Coulomb criterion.
(2) For those soil elements in which the equivalent expansive force existed, a new effective stress state could be formed when the mps was added to the gravitational effective stress in each element, while the effective stress parameters (c′ and φ′) were still set to be the variables in the Mohr–Coulomb criterion in the computation of strength reduction FEM.
Since the displacement in the strength reduction FEM calculation was only of relative significance, the displacement increment field was obtained by the difference between the two displacement fields and the adjacent value, F. It was assumed that the position of the potential sliding surface was determined according to the position of the maximum gradient of the displacement increment contour. In this program, it was one of the prominent advantages of this method that the position of a potential sliding surface could be determined by the contour density of the displacement increment.
According to the above discussion, through the algorithm implementation, which the equivalent expansive force introduced into the strength reduction FEM, the corresponding FORTRAN codes were developed and implemented for the stability analysis of unsaturated expansive soil slopes.

5. Results and Discussion

The expansive soil sample for the test was collected from Hanyin County of Shaanxi Province. The soil was brown, and its colour was deeper when it was under humid conditions and lighter when it was dry. The mineral contents of the expansive soil were identified via X-ray diffraction. The contents were mainly montmorillonite, quartz, and illite, including 51% montmorillonite and 31% quartz (Figure 5). The test employed a liquid–plastic limit combined device to measure the liquid limit (wL = 38.2), plastic limit (wP = 14.7), and shrinkage limit (wS = 12.3), respectively, and to calculate the plasticity index (IP = 23.5). The expansive soil had a high liquid limit and plasticity index. Based on the laboratory routine test measurement, the natural water content (w) of the expansive soil was 17.90%, soil natural density (ρ) was 1.94 g/cm3, and dry density (ρd) was 1.65 g/cm3.
According to the measured expansive force under different initial water contents in Figure 2, the equivalent expansion force was obtained by finding and determining the maximum expansion force, and the relation curve between the equivalent expansive force (ps)and initial water content (w) is illustrated in Figure 6.
As seen in Figure 6, the equivalent expansion force of the expansive soil sample decreased with the increase in the initial water content at the same dry density, which is consistent with the conclusions of some studies [14,15]. According to a theoretical deduction, if the initial water content was sufficiently large, then its expansive force would be zero. When the equivalent expansive force and initial water content data were matched, it was found that a certain power function relationship existed between the equivalent expansive force and initial water content, as shown in Formula (12). The equivalent expansive force (ps) can be described as follows:
p s = A 1 w λ 1
where A1 and λ1 are the fitting parameters. In this study, these parameters of the expansive soil could be determined as: A1 = 108492 and λ1 = −2.8339.
The effect parameters of the equivalent expansive force were determined by the slow direct shear tests on unsaturated and saturated expansive soil (Table 1) according to Formula (13),
m = c T c p s tan φ
By contrast, the effective cohesion force, c′, and the effective internal friction angle, φ′, were measured by the direct shear test under a completely saturated state of the expansive soil sample. Meanwhile, the total cohesion force, cT, and the equivalent expansive force, ps, were measured by the direct shear and expansive force tests under an unsaturated state of the expansive soil sample.
An effective segment of the soil and water characteristic curve of the expansive soil under humid conditions is shown in Figure 7. Therefore, a section of the fitting curve between the water content and matric suction is shown as follows:
S = A 2 w λ 2
where A2 and λ2 are the fitting parameters. In this study, these parameters of expansive soil could be determined as: A2 = 8 × 1014 and λ2 = −10.644.
Comparing Figure 6 and Figure 7, the expansive force could comprehensively reflect the macro contributions of unsaturated matric suction. In the calculation of the unsaturated expansive soil slope, the range of matric suction in the soil and water characteristic curve is generally considered for parts below 120 kPa, which affects the macroscopic effect of expansive soil slopes. Therefore, the equivalent expansive force and matric suction can reflect the principal mechanical properties of unsaturated expansive soil. In previous research, the equivalent expansive force was introduced into the strength reduction FEM. Moreover, this research quantitatively demonstrated the relevance and effect of the expansive force and matric suction to maintain the stability of an actual unsaturated expansive soil slope.
The G316 National Highway in Shaanxi Province was 603 km long. The length of expansive soil along this highway was approximately 200 km, whereas that in Ankang Hanyin was 64 km. Most of soil in the Hanjiang River basin and Ankang basin was basically typical expansive soil with a slope. The following landslide example simulated the instability of filling an embankment with an expansive soil along the G316 National Highway. The slope consisted of two parts: the lower part of the slope was 10.0 m high and had a gradient of 1.0:1.5, whereas the upper part was 3.0 m high and had a gradient of 1:10 (Figure 8). The foundation depth was 13 m, and the underground water table was located at a depth of 7.0 m.
The corresponding finite element meshes are seen in Figure 9. The elements considered were 4-node quadrilateral elements. There were 1157 elements and 1236 nodes. The boundary conditions of the displacement were as follows. (1) The bottom was a fixed boundary. (2) The normal displacement of the left-to-right lateral boundary in the horizontal direction was set to zero. The constitution of the slope soil materials adopted the Mohr–Coulomb elastic–plastic strength criterion. Table 2 provides the mechanical parameters of the soil measured during the test. m was the effective coefficient of the equivalent expansive force, χ was the parameter of suction, and φb was the internal friction angle of suction.
The calculation of Case 1 applied the non-uniform distribution of matric suction and the equivalent expansive force. Meanwhile, the non-uniform distribution scheme of matric suction was used for comparative analysis. To compare the differences and characteristics of the calculation results, Case 1-1 and Case 1-2 adopted the computing method of the strength reduction FEM, respectively, based on Bishop’s unsaturated strength formula and Fredlund’s unsaturated strength formula [28], while Case 1-3 adopted the new method of the strength reduction FEM based on the equivalent expansive force in this paper, as shown in Table 3. The stability of the expansive soil slope after rainfall was simulated and analyzed using the matric suction loss in case 2.
The distribution of the actual initial water content was obtained from the actual expansive soil slope (Figure 10). According to corresponding relationship of the equivalent expansive force curve as well as the characteristic curve of soil and water in unsaturated expansive soil in Figure 6 and Figure 7, the distributions of matric suction and the equivalent expansive force were provided, respectively (Figure 11 and Figure 12).
For the non-uniform distribution of matric suction or the equivalent expansive force, Figure 13a shows the contours of the displacement increment for Case 1-1. The safety factor (Fs = 1.64) was calculated by the strength reduction FEM based on Bishop’s unsaturated strength formula. The displacement increment contours for Case 1-2 are shown in Figure 13b. The safety factor (Fs = 1.63) was calculated by the strength reduction FEM based on Fredlund’s unsaturated strength formula. The contours of the displacement increment for Case 1-3 are illustrated in Figure 13c. The safety factor (Fs = 1.60) was calculated by the strength reduction FEM based on the equivalent expansive force.
On one hand, the slope safety factors based on three kinds of unsaturated soil strength were almost identical when the calculation results and the contours of the displacement increment in Cases 1-1, 1-2, and 1-3 were compared. On the other hand, the smaller the water content of the slope is, the larger the matric suction or the equivalent expansive force is, which leads to a greater strength of the expansive soil, placing the slope in a stable state. The safety factors were 1.64, 1.63, and 1.60 in Cases 1-1, 1-2, and 1-3, respectively.
The calculation results indicated that the strength reduction FEM based on the equivalent expansive force was successful. Therefore, the calculation results completely demonstrated the equivalence between the equivalent expansive force and matric suction in evaluating the stability of the unsaturated expansive soil slope.
After matric suction or the equivalent expansive force was lost in rainfall, the contours of the displacement increment and the potential sliding surface in Case 2 are shown in Figure 14. The safety factor was 1.08, so the slope stability was almost in a critical sliding state, and the sliding surface passed through the toe of the slope. Compared with Cases 1-1, 1-2, and 1-3, the calculated safety factor was reduced significantly. The shallow sliding mass near the slope surface is usually unsaturated. With increasing water content during rainfall conditions, its volume will change significantly, and the matric suction will decrease. As reported in the literature [8], the variation in the safety factor of unsaturated expansive soil slopes over time is studied by taking account of the shear strength reduction related to the loss of matric suction.
In summary, by means of introducing the equivalent expansive force of expansive soil into the strength reduction FEM, a calculation method was yielded for the stability of unsaturated expansive soil slopes based on the equivalent expansive force and can be used to equivalently evaluate the rainfall-induced instability of expansive soil slopes caused by the rainfall infiltration instead of using matric suction.

6. Conclusions

For three kinds of unsaturated soil strength, through the deductive analysis of the relation between matric suction and expansion force, it can be seen that there is a certain internal equivalent relation between them in the strength of unsaturated expansive soil, and this internal equivalent relation is the theoretical basis of the core algorithm analysis for the stability analysis of expansive soil slopes.
The testing device for the equivalent expansive force, which can achieve a completely saturated state under a constant volume, was developed for unsaturated expansive soil. The equivalent expansive force and matric suction tests for expansive soil were investigated under different water contents. The correlations between the equivalent expansive force and matric suction with the change in initial water content were observed, and the equivalent expansive force was found to reflect the macroscopic effects of matric suction.
Through the algorithm implementation, the equivalent expansive force was introduced into the strength reduction FEM, and the algorithm principle of the computational method on the stability of unsaturated expansive soil may indirectly reflect the contribution of matric suction using the corresponding equivalent expansive force, in which the matric suction was equivalently substituted by the equivalent expansive force.
Thus, the stability of a practical unsaturated expansive soil slope was evaluated and the equivalence between the equivalent expansive force and matric suction was demonstrated in evaluating the stability of an unsaturated expansive soil slope under the condition of non-uniform distribution of matric suction or equivalent expansive force.
According to the measurement of expansive soil expansive force by the new test device, the described method of suction-equivalent expansive force was established. Aiming at expansive soil slope engineering, the strength reduction FEM based on the equivalent expansive force, which acts as a new stability analysis method of unsaturated expansive soil slopes, can be used to evaluate the rainfall-induced instability of expansive soil slopes caused by the decrease in matric suction resulting from the rainfall infiltration, and may widen the scope and thought for the stability analysis of unsaturated slopes. Based on the established new evaluation method of the disaster-causing mechanism of expansive soil slopes with suction-equivalent expansion force, it provides certain technical support for changing slope-to-terrace construction in southern Shaanxi, so as to improve and perfect the construction technology of slope terraces in this region.

Author Contributions

Conceptualization, Q.Y. and R.L.; methodology, Q.Y.; software, R.L.; validation, S.Z. and W.B.; formal analysis, Q.Y.; investigation, S.Z.; resources, R.L.; data curation, Q.Y. and R.L.; writing—original draft preparation, Q.Y.; writing—review and editing, Q.Y., H.X. and W.B.; supervision, R.L.; project administration, R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Key R & D program of Shaanxi Province (2020ZDLGY07-03) and the China Geological Survey (DD20190268).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are all available in the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hou, T.; Xu, G.; Shen, Y.; Wu, Z.; Zhang, N.; Wang, R. Formation mechanism and stability analysis of the Houba expansive soil landslide. Eng. Geol. 2013, 161, 34–43. [Google Scholar] [CrossRef]
  2. Yang, Z.; Lv, J.; Shi, W.; Zhang, Q.; Lu, Z.; Zhang, Y.; Ling, X. Model Test Study on Stability Factors of Expansive Soil Slopes with Different Initial Slope Ratios under Freeze-Thaw Conditions. Appl. Sci. 2021, 11, 8480. [Google Scholar] [CrossRef]
  3. Kyokawa, H. A double structure model for hydro-mechano-chemical behavior of expansive soils based on the surface phenomena of mineral crystals. Eng. Geol. 2021, 294, 106366. [Google Scholar] [CrossRef]
  4. Lei, W.; Dong, H.; Chen, P.; Lv, H.; Fan, L.; Mei, G. Study on Runoff and Infiltration for Expansive Soil Slopes in Simulated Rainfall. Water 2020, 12, 222. [Google Scholar] [CrossRef] [Green Version]
  5. Xie, C.; Ni, P.; Xu, M.; Mei, G.; Zhao, Y. Combined measure of geometry optimization and vegetation for expansive soil slopes. Comput. Geotech. 2020, 123, 103588. [Google Scholar] [CrossRef]
  6. Xu, W.; Hu, R.; Li, T. Formation mechanism and stability analysis on the Hekou Sugar-refinery landslide in Xinping county Yunnan Province. Disaster Adv. 2013, 6, 20–25. [Google Scholar]
  7. Dai, Z.; Zhang, C.; Wang, L.; Fu, Y.; Zhang, Y. Interpreting the influence of rainfall and reservoir water level on a large-scale expansive soil landslide in the Danjiangkou Reservoir region, China. Eng. Geol. 2021, 11, 106110. [Google Scholar] [CrossRef]
  8. Qi, S.; Vanapalli, S.; Yang, X.; Zhou, J.; Lu, G. Stability analysis of an unsaturated expansive soil slope subjected to rainfall infiltration. Geomech. Eng. 2019, 19, 1–9. [Google Scholar]
  9. Guo, Z.; Zhao, Z. Numerical analysis of an expansive subgrade slope subjected to rainfall infiltration. Bull. Eng. Geol. Environ. 2021, 80, 5481–5491. [Google Scholar] [CrossRef]
  10. Kim, H.; Ganju, E.; Tang, D.; Prezzi, M.; Salgado, R. Matric suction measurements of compacted subgrade soils. Road Mater. Pavement Des. 2015, 16, 358–378. [Google Scholar] [CrossRef]
  11. Kim, H.; Tang, D.; Ganju, E.; Prezzi, M.; Salgado, R. Experimental investigation of matric suction in compacted fine-grained soils. Int. J. Pavement Eng. 2019, 20, 53–60. [Google Scholar] [CrossRef]
  12. Vandoorne, R.; Grbe, P.; Heymann, G. Soil suction and temperature measurements in a heavy haul railway formation. Transp. Geotech. 2021, 31, 100675. [Google Scholar] [CrossRef]
  13. Khan, M.; Nobahar, M.; Stroud, M.; Ferguson, S.; Ivoke, J. Performance Evaluation of a Highway Slope on Expansive Soil in Mississippi. Int. J. Geomech. 2022, 22, 05021005. [Google Scholar] [CrossRef]
  14. Lu, Z.; Zhang, H.; Chen, J.; Feng, M. Shear strength and swelling pressure of unsaturated soil. Chin. J. Geotech. Eng. 1992, 14, 330. [Google Scholar]
  15. Lu, Z.; Wu, X.; Sun, Y.; Yin, W. Role of swelling pressure in the shear strength theory of unsaturated soils. Chin. J. Geotech. Eng. 1997, 19, 20–27. [Google Scholar]
  16. Bishop, A.; Blight, G. Some aspects of effective stress in saturated and partly saturated soils. Geotechnique 1963, 13, 177–197. [Google Scholar] [CrossRef]
  17. Fredlund, D.; Morgenstern, N.; Widger, R. The shear strength of unsaturated soils. Can. Geotech. J. 1978, 15, 313–321. [Google Scholar] [CrossRef]
  18. Huang, R.; Wu, L. Stability Analysis of unsaturated expansive soil slope. Earth Sci. Front. 2007, 14, 129–133. [Google Scholar] [CrossRef]
  19. Liu, H.; Yin, Z. Research on analytical method of stability of expansive soil slope. Rock Soil Mech. 2010, 31, 1545–1549. [Google Scholar]
  20. Yin, H.; Xu, Q.; Li, Z. Effect of Swelling Deformation on Stability of Expansive Soil Slope. Rock Soil Mech. 2009, 30, 2506–2510. [Google Scholar]
  21. She, J.; Lu, Z.; Yao, H.; Fang, R.; Xian, S. Experimental Study on the Swelling Behavior of Expansive Soil at Different Depths under Unidirectional Seepage. Appl. Sci. 2019, 9, 1233. [Google Scholar] [CrossRef] [Green Version]
  22. Yao, H.; Zheng, S.; Li, W.; Chen, S. Parametric study on the effect of rain infiltration on stability of unsaturated expansive soil slope. Chin. J. Rock Mech. Eng. 2002, 21, 1034–1039. [Google Scholar]
  23. Noorany, I. Lateral extension of slopes in expansive soils. Geotech. Spec. Publ. 2013, 231, 642–656. [Google Scholar]
  24. Diyaljee, V. Discussion of “Lateral Extension of Compacted-Fill Slopes in Expansive Soils” by Iraj Noorany and Curt Scheyhing. J. Geotech. Geoenviron. Eng. 2015, 141, 07015034. [Google Scholar] [CrossRef] [Green Version]
  25. Xiao, J.; Yang, H.; Fu, H.; Tang, X. Comparison model test on effect of geogrid reinforcement for expansive soil slope. China J. Highw. Transp. 2014, 27, 24–31. [Google Scholar]
  26. Lu, T.; Wu, J.; Yang, S. Study on mechanism of expansive soil slope failure and numerical simulation. In Proceedings of the First International Symposium on Pavement and Geotechnical Engineering for Transportation Infrastructure, Nanchang, China, 5–7 June 2011; pp. 162–174. [Google Scholar]
  27. Dai, Z.; Guo, J.; Luo, H.; Li, J.; Chen, S. Strength Characteristics and Slope Stability Analysis of Expansive Soil with Filled Fissures. Appl. Sci. 2020, 10, 4616. [Google Scholar] [CrossRef]
  28. Li, R.; Yu, Y.; Deng, L. Discussion on stability analysis of unsaturated soil slope. Rock Soil Mech. 2007, 28, 2060–2064. [Google Scholar]
  29. Fattah, M.; Salim, N.; Irshayyid, E. Swelling Behavior of Unsaturated Expansive Soil. Transp. Infrastruct. Geotechnol. 2021, 8, 37–58. [Google Scholar] [CrossRef]
  30. Fredlund, D.G.; Rahardjo, H. Soil Mechanics for Unsaturated Soils; Wiley: New York, NY, USA, 1993. [Google Scholar]
  31. Li, R.; Chen, Y.; Duan, H. Undisturbed Soil Sample Expansion Force Testing Device Capable of Controlling Saturation; China National Intellectual Propety Administration: Beijing, China, 2012.
Figure 1. Structure of the device developed for testing EEF.
Figure 1. Structure of the device developed for testing EEF.
Water 14 01540 g001
Figure 2. Process curves of expansive force during wetting.
Figure 2. Process curves of expansive force during wetting.
Water 14 01540 g002
Figure 3. Strength surface in 3-D effective stress.
Figure 3. Strength surface in 3-D effective stress.
Water 14 01540 g003
Figure 4. Strength surface in p′-q space.
Figure 4. Strength surface in p′-q space.
Water 14 01540 g004
Figure 5. Mineral composition of expansive soil sample.
Figure 5. Mineral composition of expansive soil sample.
Water 14 01540 g005
Figure 6. Relation between the equivalent expansive force and initial water content.
Figure 6. Relation between the equivalent expansive force and initial water content.
Water 14 01540 g006
Figure 7. Characteristic curve of soil and water content.
Figure 7. Characteristic curve of soil and water content.
Water 14 01540 g007
Figure 8. Typical expansive soil slope. (a) Slope. (b) Slope profile (unit: m).
Figure 8. Typical expansive soil slope. (a) Slope. (b) Slope profile (unit: m).
Water 14 01540 g008
Figure 9. The finite element mesh.
Figure 9. The finite element mesh.
Water 14 01540 g009
Figure 10. Contours of the water content (unit: %).
Figure 10. Contours of the water content (unit: %).
Water 14 01540 g010
Figure 11. Contours of matric suction (unit: kPa).
Figure 11. Contours of matric suction (unit: kPa).
Water 14 01540 g011
Figure 12. Contours of the equivalent expansive force (unit: kPa).
Figure 12. Contours of the equivalent expansive force (unit: kPa).
Water 14 01540 g012
Figure 13. Contours of the displacement increment in Case 1. (a) Case 1-1. (b) Case 1-2. (c) Case 1-3.
Figure 13. Contours of the displacement increment in Case 1. (a) Case 1-1. (b) Case 1-2. (c) Case 1-3.
Water 14 01540 g013
Figure 14. Contours of the displacement increment in Case 2.
Figure 14. Contours of the displacement increment in Case 2.
Water 14 01540 g014
Table 1. The parameters of m.
Table 1. The parameters of m.
ParameterWater Content of 17.9% Water Content of 23.48%Saturated Soil
c/kPa31.922.813.1
φ/(°)16.115.415.0
ps/kPa3116
m2.12.2
Mean value of m2.15
Table 2. Parameters of the soil.
Table 2. Parameters of the soil.
Soil ParameterSoil Above the Water LevelSoil Below the Water Level
γ/(kN/m3)21.3019.40
E/MPa20.3818.40
ν0.300.35
c′/kPa13.1013.10
φ′/(º)15.0015.00
φb/(º)10.6015.00
χ0.701.00
m2.15-
A1108492-
λ1−2.8339-
A28 × 1014-
λ2−10.644-
Table 3. Description of cases and safety factors.
Table 3. Description of cases and safety factors.
Case Matric Suction Equivalent
Expansive Force
Calculation MethodSafety Factor
Fs
Case 1-1Non-uniform distribution-Strength Reduction FEM based on Bishop’s unsaturated strength 1.64
Case 1-2Non-uniform distribution-Strength Reduction FEM based on Fredlund’s unsaturated strength 1.63
Case 1-3-Non-uniform distributionStrength Reduction FEM based on the equivalent expansive force1.60
Case 2Suction loss-Strength Reduction FEM1.08
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, Q.; Li, R.; Zhang, S.; Li, R.; Bai, W.; Xiao, H. Algorithm Implementation of Equivalent Expansive Force in Strength Reduction FEM and Its Application in the Stability of Expansive Soil Slope. Water 2022, 14, 1540. https://doi.org/10.3390/w14101540

AMA Style

Yang Q, Li R, Zhang S, Li R, Bai W, Xiao H. Algorithm Implementation of Equivalent Expansive Force in Strength Reduction FEM and Its Application in the Stability of Expansive Soil Slope. Water. 2022; 14(10):1540. https://doi.org/10.3390/w14101540

Chicago/Turabian Style

Yang, Qiang, Rongjian Li, Shibin Zhang, Rongjin Li, Weishi Bai, and Huiping Xiao. 2022. "Algorithm Implementation of Equivalent Expansive Force in Strength Reduction FEM and Its Application in the Stability of Expansive Soil Slope" Water 14, no. 10: 1540. https://doi.org/10.3390/w14101540

APA Style

Yang, Q., Li, R., Zhang, S., Li, R., Bai, W., & Xiao, H. (2022). Algorithm Implementation of Equivalent Expansive Force in Strength Reduction FEM and Its Application in the Stability of Expansive Soil Slope. Water, 14(10), 1540. https://doi.org/10.3390/w14101540

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