A comparison between the Lagrangian and Eulerian kernels used in meshfree approximation is illustrated. Then with the Lagrangian approximation an accelerated meshfree formulation for dynamic as well as static analysis of slope failure at finite deformation range is presented. The present formulation employs the stabilized conforming nodal integration (SCNI) to ensure the spatial stability and improve computational efficiency. The degradation of slope materials is characterized by the coupled elasto-plastic damage model via the rate form

where the Druker-Prager yield function and isotropic strain-based damage criterion are used. The non-local nature of meshfree approximation and SCNI scheme enables the current approach to effectively simulate the shear band arising from strain localization during slope failure process. Finally a numerical example is given to demonstrate the proposed method.