Finite element analysis of slope stability by expanding the mobilized principal stress Mohr's circles - Development, encoding and validation

In recent years, finite element analysis is increasingly being proposed in slope stability problems as a competitive method to traditional limit equilibrium methods (LEMs) which are known for their inherent deficiencies. However, the application of finite element method (FEM) to slope stability as a strength reduction method (SRM) or as finite element limit analysis (FELA) is not always a success for the drawbacks that characterize both methods. To increase the performance of finite element analysis in this problem, a new approach is proposed in this paper. It consists in gradually expanding the mobilized stress Mohr's circles until the soil failure occurs according to a prescribed non-convergence criterion. The present approach called stress deviator increasing method (SDIM) is considered rigorous for three main reasons. Firstly, it preserves the definition of the factor of safety (FOS) as the ratio of soil shear strength to the mobilized shear stress. Secondly, it maintains the progressive development of shear stress resulting from the increase in the principal stress deviator on the same plane, on which the shear strength takes place. Thirdly, by introducing the concept of equivalent stress loading, the resulting trial stresses are checked against the violation of the actual yield criterion formed with the real strength parameters rather than those reduced by a trial factor. The new numerical procedure was encoded in a Fortran computer code called S(4)DINA and verified by several examples. Comparisons with other numerical methods such as the SRM, gravity increasing method (GIM) or even FELA by assessing both the FOS and contours of equivalent plastic strains showed promising results. (C) 2022 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V.