Accurately modeling soil-fluid coupling under large deformations is critical for understanding and predicting phenomena such as slope failures, embankment collapses, and other geotechnical hazards. This topic has been studied for decades and remains challenging due to the nonlinear responses of geotechnical structures, which typically result from plastic yielding and finite deformation of the soil skeleton. In this work, we comprehensively summarize the theory involved in the soil-fluid coupling problem. Within a finite strain framework, we employ an elasto-plastic constitutive model with linear hardening to represent the solid skeleton and a nearly incompressible model for water. The water content influences the behavior of the solid skeleton by affecting its cohesion. The governing equations are discretized by material point method and two sets of material points are employed to independently represent solid skeleton and fluid, respectively. The proposed method is validated by comparing simulation results with experimental results for the impact of water on dry soil and wet soil. The capability of the method is further demonstrated through two cases: (1) the impact of a rigid body on saturated soil, causing water seepage, and (2) the filling of a ditch, which considers the erosion of the foundation. This work may provide a versatile tool for analyzing the dynamic responses of fluid and solid interactions, considering both mixing and separation phenomena.

Marine structures are commonly situated near the mildly sloping sandy seabed characterized by the slope angles (alpha) not exceeding 10 degrees. The seabed liquefaction can be triggered due to the generation of the excess pore water pressure (EPWP), posing a threat to the stability of marine structures. This study focuses on the analysis of waveinduced liquefaction in the mildly sloping (MS) sandy seabed. A dynamic poro-elasto-plastic seabed model is developed to simulate the behavior of the MS sandy seabed under wave loading. The results indicates that the loading cycle required to trigger the initial liquefaction decreased as the position moved from the toe towards the crest of the MS sandy seabed. The amplitude of shear stress increases with the loading cycle and tends to increase with growing alpha before liquefaction, resulting in a slower accumulation of EPWP with larger alpha. Both the horizontal and vertical displacements induced by wave action reach the maximum at the crest of the sloping seabed. Notably, the horizontal displacement is much greater than the vertical displacement in the seabed under wave action. The displacement of the MS sandy seabed depends on not only the shear stress amplitude developed in the soils but also the accumulation of EPWP required to trigger the liquefaction in the seabed.

The 2011 off the Pacific Coast of Tohoku earthquake caused extensive liquefaction damage to reclaimed land along the Tokyo Bay coast, even though it was approximately 400 km from the epicenter. The characteristics of the liquefaction damage include the fact that liquefaction occurred in soils with a high percentage of fine particles and that the distribution of liquefied and nonliquefied areas was nonuniform. The factors contributing to such nonuniform liquefaction damage included the heterogeneity of the ground materials and their depositional conditions, and the effects of the long earthquake duration. Although these points are certainly valid as reasons for the occurrence of severe liquefaction damage, they do not fully explain the mechanisms of the liquefaction of the fine-grained soils, or the localized extent of the liquefaction. To elucidate the severe and nonuniform damage, seismic response analyses of a multi-layered ground were conducted focusing on the stratigraphic irregularities in the ground beneath Urayasu city. The results showed that the thicker and softer sedimentary layers amplify the slightly long-period component of the seismic motion and increase the shaking at the ground surface. Moreover, the wave propagation in the ground became very complicated owing to the focal effect caused by the refractions and reflections of body waves at the stratum boundary, surface wave excitation at the base of the slope, and amplified interference between body and surface waves. This complex wave propagation contributed to nonuniform surface ground shaking and severe liquefaction damage. In addition, surface waves, which consist primarily of slightly long-period components, can propagate far and wide; as such, they triggered extensive damage owing to delayed shaking phenomena that continue even after the earthquake. The analysis results suggested that multidimensional elasto-plastic seismic response analyses considering stratigraphic irregularities are important for detailed seismic evaluation.