On the Computational Precision of Finite Element Algorithms in Slope Stability Problems

Although the finite element method (FEM) has been used extensively to analyse the slope stability problems, the computational precision and definition of failure are still two main key concepts of finite element algorithms that attract the attention of researchers. In this paper, the modified Euler algorithm and the explicit modified Euler algorithm with stress corrections are used to analyse two dimensional (2D) slope stability problems with the associated flow rule, based on the shear strength reduction method. The rounded hyperbolic Mohr-Coulomb (M-C) yield surface is applied. Effects of the element type and various definitions of failure on the computational precision of 2D slope stability problems are evaluated. Conclusions can be drawn that the modified Euler scheme is applicable when the factor of safety (FOS) is small; however, the explicit modified Euler algorithm with stress corrections is more precise if the factor of safety is relatively large. The fully integrated quadrilateral isoparametric element is better than the triangular element in terms of the precision. With respect to the definition of failure, the displacement mutation of the characteristic point combining with the continuums of the plastic zone can be regarded as a reliable definition of failure and can be widely used to perform and analyse numerical simulations of slope stability problems.