A u-p form finite element, material point and finite volume coupled method (FEM-MPM-FVM) is developed for stability and accuracy when simulating large deformations in saturated soil. In the coupled FEM-MPM-FVM, solid phase with small and large deformation is discretized by FEM and MPM respectively, while fluid phase is discretized by FVM. The numerical implementation of an explicit sequentially coupled algorithm is developed, taking into consideration the balance between computational cost and numerical stability. The computational efficiency and accuracy of simulating soil without large deformation can be improved by FEM-MPM hybrid contact method. The inflow/outflow boundary conditions of pore fluid can be accurately simulated using FVM within the Eulerian framework. The method is validated through the simulation of the Selborne landslide experiment. The simulation achieves excellent agreement with the experiment measurements and observations.