In an effort to study undrained post-liquefaction shear deformation of sand, the discrete element method (DEM) is adopted to conduct undrained cyclic biaxial compression simulations on granular assemblies consisting of 2D circular particles. The simulations are able to successfully reproduce the generation and eventual saturation of shear strain through the series of liquefaction states that the material experiences during cyclic loading after the initial liquefaction. DEM simulations with different deviatoric stress amplitudes and initial mean effective stresses on samples with different void ratios and loading histories are carried out to investigate the relationship between various mechanics- or fabric-related variables and post-liquefaction shear strain development. It is found that well-known metrics such as deviatoric stress amplitude, initial mean effective stress, void ratio, contact normal fabric anisotropy intensity, and coordination number, are not adequately correlated to the observed shear strain development and, therefore, could not possibly be used for its prediction. A new fabric entity, namely the Mean Neighboring Particle Distance (MNPD), is introduced to reflect the space arrangement of particles. It is found that the MNPD has an extremely strong and definitive relationship with the post-liquefaction shear strain development, showing MNPD’s potential role as a parameter governing post-liquefaction behavior of sand.

Granular materials deposited under gravity usually exhibit inherent fabric anisotropy, which leads to profoundly varying strength under different loading directions in regard to the bedding plane that has been observed through experiments and numerical simulations. Classical isotropic strength criteria that are only stress dependent cannot describe such anisotropy in strength. This study introduces an anisotropic strength variable Λ that measures the distance between the bedding plane and the maximum shear stress ratio plane, thus incorporating both the stress tensor and the material orientation, and guarantees objectivity. The strength variable Λ is adopted to construct a general formulation that allows for the extension of many existing isotropic strength criteria to become anisotropic. The formulation requires only the following two parameters: kf0, reflecting the minimum strength of the material, and a, reflecting the intensity of inherent anisotropy. The two parameters can be conveniently cal...

Suction caissons are widely used foundations in offshore engineering. The change in excess pore pressure and seepage field caused by penetration and suction significantly affects the soil resistance around the caisson wall and tip, and also affects the deformation of the soil within and adjacent to the caisson. This study uses Arbitrary Lagrangian–Eulerian (ALE) large deformation solid-fluid coupled FEM to investigate the changes in suction pressure and the seepage field during the process of the suction caisson installation in sand. A nonlinear Drucker-Prager model is used to model soil, while Coulomb friction is applied at the soil-caisson interface. The ALE solid-fluid coupled FEM is shown to be able to successfully simulate both jacked penetration and suction penetration caisson installation processes in sand observed in centrifuge tests. The difference in penetration resistance for jacked and suction installation is found to be caused by the seepage and excess pore pressure generated during the suction caisson installation, highlighting the importance of using solid-fluid coupled effective stress-based analysis to consider seepage in the evaluation of suction caisson penetration.