The adjustment of nonlinear stress is not taken into account in the limit equilibrium analysis of rigid body, and it is difficult to get sliding safety factor with finite element method. To solve these problems, based on multi-grid method, structure grid for nonlinear finite element method and sliding surface grid for calculating safety factor are established separately. So the sliding safety factor of arbitrary sliding surface or sliding blocks can be obtained. Thus limit equilibrium analysis can be integrated with nonlinear finite element method. This algorithm has been implemented in TFINE, a code for 3-D nonlinear finite element method. The calculation efficiency and the influence on result precision by grid are analyzed. Then it is applied to the stability analysis of Jinping high slope. The result shows that multi-grid method is effective and more reasonable.