Soils with strain-softening behavior — manifesting as a reduction of strength with increasing plastic strain — are commonly found in the natural environment. For slopes in these soils,a progressive failure mechanism can occur due to a reduction of strength with increasing strain. Finite element method based numerical approaches have been widely performed for simulating such failure mechanism,owning to their ability for tracing the formation and development of the localized shear strain. However,the reliability of the currently used approaches are often affected by poor convergence or significant mesh-dependency,and their applicability is limited by the use of complicated soil models. This paper aims to overcome these limitations by developing a finite element approach using a local arc-length controlled iterative algorithm as the solution strategy. In the proposed finite element approach,the soils are simulated with an elastoplastic constitutive model in conjunction with the Mohr-Coulomb yield function. The strain-softening behavior is represented by a piece-wise linearrelationship between the Mohr-Coulomb strength parameters and the deviatoric plastic strain. To assess the reliability of the proposed finite element approach,comparisons of the numerical solutions obtained by different finite element methods and meshes with various qualities are presented. Moreover,a landslide triggered by excavation in a real expressway construction project is analyzed by the presented finite element approach to demonstrate its applicability for practical engineering problems. Soils with strain-softening behavior - manifesting as a reduction of strength with increasing plastic strain - are commonly found in the natural environment. For slopes in these soils, a progressive failure mechanism can occur due a reduction of strength with increasing strain. method based numerical approaches have been widely done for simulating such failure mechanism, owning to their ability for tracing the formation and development of the localized shear strain. , and their applicability is limited by the use of complicated soil models. This paper aims to overcome these limitations by developing a finite element approach using a local arc-length controlled iterative algorithm as the solution strategy. In the proposed finite element approach, the soils are simulated with an elastoplastic constitutive model in conjunction with the Mohr-C oulomb yield function. The strain-softening behavior is represented by a piece-wise linear relationship between the Mohr-Coulomb strength parameters and the deviatoric plastic strain. To assess the reliability of the proposed finite element approach, comparisons of the numerical solutions obtained by different finite element methods and meshes with various qualities are presented. Moreover, landslide triggered by excavation in a real expressway construction project is analyzed by the presented finite element approach to demonstrate its applicability for practical engineering problems.