Abstract: Large landslides often cause huge human casualties and property losses. It is of great significance to carry out numerical simulation of the entire process of landslide for quantitative landslide risk assessment. This paper establishes a two-dimensional material point method(MPM)numerical model for Shenzhen"12 ·20"landslide. This model incorporates the updated Lagrangian governing equation with the ideal elastoplastic model under the Drucker-Prager yield criterion. The equivalent internal friction angle is derived by considering the effect of pore water pressure. The K-means clustering algorithm is used to identify the potential slip surface in combination with the relative displacement of the material points. The computed results show that the landslide slip surface,instability and movement characteristics,impact strength,and deposition feature of the Shenzhen landslide obtained by the MPM simulation analysis are in good agreement with the actual observation results,which implies that MPM is a reliable numerical simulation tool to reproduce the entire life cycle of landslide stability,instability,movement,impact and deposition process.