Landslides can cause major economic damage and a large number of casualities as it is possible to see from past events occurred all over the world. Being able to predict these kind of hazards would then suppose the achievement of great benefits. Here a model that combines a depth integrated description of the soil-pore fluid mixture together with a set of 1D models dealing with pore pressure evolution within the soil mass is presented. The mathematical model is based on the Biot-Zienkiewicz equations, from where a depth averaged model is derived. Concerning the material behaviour, the approach used is the one suggested by the Perzyna viscoplasticity, which has been extensively used in the past to model solid behaviour prior to failure. In this framework, a simple shear rheological model is derived, providing the basal friction needed in depth integrated models. The Smoothed Particle Hydrodynamics (SPH) has been the numerical technique chosen to spatially discretised the depth integrated equations of the mathematical model. The purpose of this work is to apply the SPH depth integrated numerical model, together with the sub-model that predicts the evolution of the pore water pressure inside the landslide, to simulate the propagation phase of the Aberfan flowslide occurred in 1966.
The choice of a pure cohesive or a pure frictional viscoplastic model to represent the rheological behaviour of a flowslide is of paramount importance in order to obtain accurate results for real cases. The principal Goal of the present work is to clarify the influence of the type of viscous model—pure cohesive versus pure frictional—with the numerical reproduction of two different real flowslides that occurred in 1966: the Aberfan flowslide and the Gypsum tailings impoundment flowslide. In the present work, a depth-integrated model based on the v-pw Biot–Zienkiewicz formulation, enhanced with a diffusion-like equation to account for the pore pressure Evolution within the soil mass, is applied to both 1966 cases. For the Aberfan flowslide, a frictional viscous model based on Perzyna viscoplasticity is considered, while a pure cohesive viscous model (Bingham model) is considered for the case of the Gypsum flowslide. The numerical approach followed is the SPH method, which has been enriched by adding a 1D finite difference grid to each SPH node in order to improve the description of the pore water evolution in the propagating mixture. The results obtained by the performed simulations are in agreement with the documentation obtained through the UK National Archive (Aberfan flowslide) and the International Commission of large Dams (Gypsum flowslide).