Soil liquefaction describes a loss of strength of saturated sand upon sudden or cyclic loading. A slight disturbance of such a soil’s fabric might lead to severe damage, e.g. the collapse of sea dikes. Accurate modeling of the state transition between saturated soil and a liquefied soil-water mixture, as well as post-liquefaction phenomena, is crucial for the prediction of such damage. However, developing an appropriate numerical model remains a challenging problem, especially when the simulation involves dynamic large deformation processes. In order to make a first step towards an accurate simulation of soil liquefaction, a two-phase formulation of the finite element method (FEM) in conjunction with the elastoplastic UBC3D-PLM model is investigated. The performance of this approach is analyzed based on a shaking table benchmark.@en

Dynamic two-phase interaction of soil can be modelled by a displacement-based, two-phase formulation. The finite element method together with a semi-implicit Euler–Cromer time-stepping scheme renders a discrete equation that can be solved by recursion. By experience, it is found that the CFL stability condition for undrained wave propagation is not sufficient for the considered two-phase formulation to be numerically stable at low values of permeability. Because the stability analysis of the two-phase formulation is onerous, an analysis is performed on a simplified two-phase formulation that is derived by assuming an incompressible
pore fluid. The deformation of saturated porous media is now captured in a single, second-order partial differential equation, where the energy dissipation associated with the flow of the fluid relative to the soil skeleton is represented by a damping term. The paper focuses on the different options to discretize the damping term and its effect on the stability criterion. Based on the eigenvalue analyses of a single element, it is observed that in addition to the CFL stability condition, the influence of the permeability must be included. This paper introduces a permeability-dependent stability criterion. The findings are illustrated and validated with an example for the dynamic response of a sand deposit.@en