Numerical Simulation of Liquefaction Behaviour in Coastal Reclaimed Sediments

This study presents a validated numerical investigation into the seismic liquefaction potential of fine-grained reclaimed sediments commonly encountered in coastal, containment, and reclamation projects. Fine-grained reclaimed sediments pose a particular challenge for seismic liquefaction assessment due to their low permeability, high fines content, and complex cyclic response under earthquake loading. A fully coupled, nonlinear finite element model was developed using the Pressure-Dependent Multi-Yield (PDMY) constitutive framework, calibrated against laboratory Cyclic Direct Simple Shear (CDSS) tests and verified using in situ Cone Penetration Tests with pore pressure measurement (CPTu). The model effectively captured the dynamic response of saturated sediments, including excess pore pressure generation, cyclic mobility, and post-liquefaction behavior, under three earthquake ground motions: Livermore, Chi-Chi, and Loma Prieta. Results showed that near-surface layers (0–2.3 m) experienced full liquefaction within two to three cycles, with excess pore pressure ratios (Ru) approaching 1.0 and peak pressures closely matching laboratory data with less than 10% deviation. The numerical approach revealed that traditional CPT-based cyclic resistance methods underestimated liquefaction susceptibility in intermediate layers due to limitations in accounting for pore pressure redistribution, evolving permeability, and seismic amplification effects. In contrast, the finite element model captured progressive strength degradation, revealing strength gain in deeper layers due to consolidation, while upper zones remained vulnerable due to low confinement and resonance effects. A critical threshold of Ru ≈ 0.8 was identified as the onset of rapid shear strength loss. The findings confirm the advantage of advanced numerical modeling over empirical methods in capturing the complex cyclic behavior of reclaimed sediments and support the adoption of performance-based seismic design for such geotechnically sensitive environments.