Poro-elastic and poro-elasto-plastic modeling of sandy seabed under wave action

Abstract

The seabed could be liquefied due to the wave-induced accumulation of residual excess pore water pressure (EPWP), resulting in the instability of offshore and marine structures. Thus, it is of importance to consider both the buildup of EPWP and the associated deformation behavior of the seabed. In this study, a dynamic poro-elasto-plastic seabed model, integrated with a state-dependent plasticity model, is developed in COMSOL Multiphysics. The numerical model of seabed is validated against the analytical solutions and model tests. The elliptical trajectories of soil particle at different depths of the poro-elastic seabed under the progressive wave are illustrated. The ratio of the amplitude of vertical displacement to that of horizontal displacement generally decreases as the depth increases for the poro-elastic seabed. In the poro-elasto-plastic seabed, however, the soil particles tend to move following the direction of the wave propagation. Moreover, the soil particles at the shallow part of the seabed go downwards before the initial liquefaction, while the soil particles at the lower position initially go upwards and then turn to go downwards when the initial liquefaction occurs. Particularly, the surface seabed with lower relative density is more prone to lose the stability, resulting in a dramatic upward vertical displacement.