By simulating the sliding face with point contact

a numerical method for slope stability analysis with a prescribed sliding plane is proposed combining with the FEM and strength reduction theory. For a specific slope

the shear and normal forces of the sliding plane are represented by the forces of the shear and normal springs. Combined with the theory of strength reduction and the Mohr-Coulomb strength criterion

numerical results can be obtained by iteration of the spring stiffness under equilibrium condition. Compared with limit equilibrium method

the present method avoids the limits due to the slice-force assumption. Compared with shear strength reduction method in a finite element model

the present method overcomes the sliding piles and the convergence difficulty. The present method can provide body stress

boundary force and stability coefficient under limit state. The amplitude

direction and action point of the landslide thrust on an arbitrary section of the landslide

can also be obtained under a given safety factor. The present method

with strict theoretical basis

provides some guidance to landslide harness in practical engineering applications.