Traditional geotechnical analyses for landslides involve failure prediction (i.e. onset of failure) and the design of structures that can safely withstand the applied loads. The analyses provide limited information on the post-failure behaviour. Modern numerical methods are able to simulate large mass movements and there is an opportunity to utilise such methods to evaluate the risks of catastrophic damage if a landslide occurs. In this paper, various large-deformation analysis methods are introduced and their applicability for solving landslide problems is discussed. Since catastrophic landslides often involve seepage forces, consideration of the coupled behaviour of soil and pore fluid is essential. Two approaches to model soil–pore fluid coupling in large-deformation analysis using the material point method (MPM) are introduced. An example simulation is presented for each approach; one on a model levee failure and the other on a natural cut slope failure (the Selborne experiment conducted by Cooper and co-workers in 1998). In the levee failure case, MPM simulation was able to capture a complex failure mechanism including the development of successive shear bands. The simulation was also able to predict excess pore pressure generation during the failure propagation and the subsequent consolidation stage. The simulations demonstrated the importance of the dilation characteristics of soil as well as changes in geometry for the post-failure behaviour. In the Selborne case, MPM was able to simulate the progressive failure of brittle, overconsolidated clay. The evolution of shear stresses along the failure surface was also captured by the MPM. The changes in the pore pressure and the actual shape of the failure surface were simulated by the MPM. The importance of accurately modelling the shear band within the MPM framework is highlighted.