Research PaperNon-associated plasticity analysis of slope stability under steady unsaturated flow conditions
Introduction
Many soil slopes are unsaturated in nature, and soil strength can be enhanced by the matric suction, which consequently contributes to slope stability (Fredlund and Rahardjo, 1993). To evaluate the characteristics of unsaturated soil, great efforts have been made to find a representation of the shear strength. Fredlund et al. (1978) proposed an extended Mohr-Coulomb (MC) criterion with the linear shear strength to simulate unsaturated soil. Then, Vanapalli et al., 1996, Fredlund et al., 1996 proposed the nonlinear shear strength of unsaturated soil based on the soil–water characteristic curves. Furthermore, Lu et al. (2010) predicted an effective MC model that can directly introduce the suction stress into the shear strength without new assumptions and parameters. Based on the MC criterions mentioned above, the effect of the matric suction on slope stability has been investigated widely (Lu and Godt, 2008, Vahedifard et al., 2016, Li and Yang, 2018, Yuan and Du, 2018).
Limit analysis is an attractive choice to estimate the stability of unsaturated soil slopes. It usually generates the upper bound (UB) and lower bound (LB) solutions by the statically admissible field and the kinematically admissible field (Michalowski, 1995, Zhou and Wang, 2017, Qin and Chian, 2017, Li and Jiang, 2020), respectively. Compared with the limit equilibrium method, the bound theorems of classical plasticity produce more meaningful solutions considering the constitutive relationship of soil. Moreover, unlike the numerical simulation, limit analysis is easier to operate and has a lower computational cost.
Chen (1975) made such comments: “The slip line solution is not necessarily the true solution. If the associated flow rule is employed and the resulting stress–strain rate equations can be integrated to yield a kinematically admissible velocity field, the slip line solution is an UB solution. If, in addition, the slip line stress field can be extended over the entire half-space of the soil domain such that the equilibrium equations, the stress boundary conditions and the yield condition are satisfied, the slip line solution is also a LB solution and is hence the actual solution.” Therefore, if the slip line field of unsaturated soil slopes that satisfies the stress boundary and the velocity boundary simultaneously can be found out, the accurate solution will be directly predicted, not the narrow range determined by the conventional limit analysis.
In terms of limit analysis, the plastic flow rule plays an important role (Mizuno and Chen, 1983). The conventional limit analysis usually applies the associated flow rule to estimate slope stability. In actual, soil, when modelled as a rigid-perfectly plastic solid, does not obey the associated flow rule, and the ultimate load is no longer unique (Drescher and Detournay, 1993). Meanwhile, slope stability involving the non-associated plasticity is difficult to resolve, especially when the friction angle is much larger than the dilatancy angle (Chen et al., 2019). To overcome the difficulty, Davis and Lee (1968) introduced the modified shear strength of MC model to convert the non-associated plasticity into the equivalently associated plasticity, and this approach has been widely adopted to evaluate slope stability with the non-associated flow rule (Yang and Huang, 2009, Oberhollenzer et al., 2018, Wang et al., 2019).
Although many studies mentioned above have made great contributions to the stability assessment of unsaturated soil slopes and the non-associated plasticity problem, the ultimate load of unsaturated soil slopes with the non-associated plasticity has not been addressed. The aim of this study is to search for a slip line field of unsaturated soil slopes with the non-associated plasticity using the method of characteristics combined with boundary value problems. The slip line field should satisfy the statically admissible stress field and the kinematically admissible velocity field simultaneously. Thus, the accurate solution could be directly predicted, which is more effective than the UB and LB solutions. Although the kinematically admissible velocity field is not strictly fulfilled here, the presented solution is between the UB and LB solutions. This solution can be considered to be near the accurate one. Meanwhile, the effective MC criterion proposed by Lu et al. (2010) is modified to analyze the suction stress of unsaturated soil, and the Davis approach (Davis and Lee, 1968) is applied to account for the non-associated flow rule. Then, the effect of the matric suction and the dilation on slope stability can be investigated, and the results are compared with the existing methods.
Section snippets
Matric suction of unsaturated soil
Based on the characteristic curve of the suction stress, Lu and Likos (2004) proposed the effective stress of both saturated and unsaturated soil, of which a unified representation can be written aswhere is the effective stress; is the total stress; is the pore air pressure; and is the suction stress. A closed-form function for has been established as
where is the inverse of the air entry pressure that varies
Non-associated plasticity analysis
The conventional limit analysis assumes that soil obeys the associated flow rule. Thus, the stress characteristics and the velocity characteristics in a slip line field are consistent (Li et al., 2019). However, the deformation and the failure property of soil can be better simulated under the non-associated flow rule (Drescher and Detournay, 1993). Davis approach (Davis and Lee, 1968) adopted the modified strength parameters ( and ) of the associated MC model to approximate the
The method of characteristics combined with boundary value problems
The non-associated MC model that accounts for the effect of the matric suction and the dilation is adopted in the proposed approach. Then, the following basic assumptions are made: (i) soil is an ideal rigid, perfectly plastic material; (ii) the plane strain condition is considered; and (iii) soil slopes obey the toe failure mechanism.
Comparison with existing methods
To validate the proposed approach, comparisons are conducted with the existing methods and are divided into three parts: (i) the comparison of the failure mechanism between the proposed approach and the finite element limit analysis (FELA), where soil is under the associated flow rule () and the no-suction condition; (ii) the comparison with the results proposed by Vahedifard et al. (2016), where a vertical steady unsaturated flow is considered and the effect of non-associated plasticity
Results and discussion
In this section, a parametric analysis is conducted to investigate the influence of the vertical steady flow and the non-associated plasticity on the ultimate load of unsaturated soil slopes. The non-dimensional ratio is introduced to present a variation in the ultimate load along with kN/m3, kN/m3 and m. In addition, the effective cohesion and the effective friction angle in the proposed approach should be more than zero. Thus, for sand slopes, the effective
Conclusions
This study adopts the method of characteristics combined with boundary value problems to present the ultimate load of unsaturated soil slopes considering the non-associated plasticity. To do this, the non-associated MC failure criterion is introduced to describe the shear strength of unsaturated soil with the non-associated plasticity. Moreover, three types of value boundary problems are resolved using the finite difference method to build the slip line filed that satisfies the static boundary
CRediT authorship contribution statement
Chengchao Li: Conceptualization, Methodology, Investigation, Data curation, Writing - original draft, Writing - review & editing. Pengming Jiang: Formal analysis, Writing - original draft, Writing - review & editing. Aizhao Zhou: Writing - original draft, Writing - review & editing.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgement
The authors would like to thank the editor and the anonymous reviewers for their valuable suggestions in improving this study. In addition, this study is financially supported by the National Natural Science Foundation of China (51979128) and the Excellent Scientific and Technological Innovation Team of Jiangsu Higher Education Institutions (Jiangsu Techer, 2017 No. 51).
References (34)
- et al.
A modified Davis approach for geotechnical stability analysis involving non-associated soil plasticity
Géotechnique
(2019) - et al.
Undrained seismic bearing capacity of strip footing placed on homogeneous and heterogeneous soil slopes by finite element limit analysis
Comput. Geotech.
(2019) - et al.
Rigorous solution of slope stability under seismic action
Comput. Geotech.
(2019) - et al.
Ultimate load of nonhomogeneous slopes determined by using the method of characteristics
Eng. Geol.
(2019) - et al.
Stability of 3D slope under steady unsaturated flow condition
Eng. Geol.
(2018) - et al.
Eccentric bearing capacity of embedded strip footings placed on slopes
Comput. Geotech.
(2020) - et al.
Plasticity analysis of slope with different flow rule
Comput. Struct.
(1983) - et al.
Finite element analyses of slope stability problems using non-associated plasticity
J. Rock Mech. Geotech. Eng.
(2018) - et al.
Three-dimensional stability analysis of slope in unsaturated soils considering strength nonlinearity under water drawdown
Eng. Geol.
(2018) - et al.
Effective stress-based upper bound limit analysis of unsaturated soils using the weak form quadrature element method
Comput. Geotech.
(2018)
Lower bound limit analysis of wedge stability using block element method
Comput. Geotech.
Limit Analysis and Soil Plasticity
Theories of plasticity and failure of soil masses
Limit load in translational failure mechanisms for associative and non-associative materials
Géotechnique
Soil Mechanics for Unsaturated Soils
The shear strength of unsaturated soils
Can. Geotech. J.
The relationship of the unsaturated soil shear strength to the soil-water characteristic curve
Can. Geotech. J.
Cited by (6)
Nonlinear Analysis of Slope Stability Considering Geomaterial Dilatancy
2024, International Journal of GeomechanicsMethod of Rigorous Characteristics for Seismic Bearing Capacity of Strip Footings Placed Adjacent to Homogeneous Soil Slopes
2022, International Journal of GeomechanicsOblique and eccentric bearing capacity of embedded strip footings using the method of rigorous characteristics
2022, International Journal for Numerical and Analytical Methods in GeomechanicsApproaches to Slope Stability Analysis Considering the Effects of Dilatancy and Strength Non-linearity: A Review
2022, Lecture Notes in Civil EngineeringA comprehensive investigation on bearing capacity of shallow foundations on unsaturated fly ash slopes adopting finite element limit analysis
2022, European Journal of Environmental and Civil Engineering