Explicit meshfree u - p w solution of the dynamic Biot formulation at large strain.

In this paper, an efficient and robust methodology to simulate saturated soils subjected to low-medium frequency dynamic loadings under large deformation regime is presented. The coupling between solid and fluid phases is solved through the dynamic reduced formulation u - p w (solid displacement - pore water pressure) of the Biot's equations. The additional novelty lies in the employment of an explicit two-steps Newmark predictor-corrector time integration scheme that enables accurate solutions of related geomechanical problems at large strain without the usually high computational cost associated with the implicit counterparts. Shape functions based on the elegant Local Maximum Entropy approach, through the Optimal Transportation Meshfree framework, are considered to solve numerically different dynamic problems in fluid saturated porous media.

Modelling of Fluidised Geomaterials: The Case of the Aberfan and the Gypsum Tailings Impoundment Flowslides.

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 - p w 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).