Numerical framework for calculating stiffness of liquefied sands during post-shaking re-solidification
File(s) 1-s2.0-S0266352X2600368X-main.pdf (12.27 MB)
Published version
Author(s)
MA, shengjie
Taborda, David
Kontoe, Stavroula
Type
Journal Article
Abstract
Accurate prediction of earthquake-induced settlement in saturated sands remains challenging, primarily due to difficulties in capturing the coupled evolution of excess pore pressure, hydraulic conductivity, and compressibility during and after seismic liquefaction. In particular, quantifying stiffness recovery and consolidation behaviour during the post-shaking re-solidification phase remains problematic, as this stage is characterised by high hydraulic gradients and partially drained conditions. This study presents a physics-constrained numerical and inverse analysis framework to characterise post-shaking re-solidification in liquefied sand deposits. A hybrid explicit-implicit finite difference method (HEI-FDM) is developed to solve the nonlinear one-dimensional consolidation equation governing excess pore pressure dissipation and settlement efficiently and robustly. Building on this forward solver, a joint two-parameter inversion framework is proposed to identify the evolution of constrained modulus and hydraulic conductivity by simultaneously matching excess pore pressure time histories at multiple depths and surface settlement records. The framework is applied to four level-ground free-field dynamic centrifuge tests, enabling systematic investigation of stiffness recovery, hydraulic conductivity correction, and consolidation front propagation. Error-surface mapping and multi-start optimisation demonstrate that the inverse problem is well-constrained, yielding a unique and physically interpretable parameter set. Results show that stiffness inferred directly from conventional compression tests overestimates the stiffness of liquefied sand under
partially drained conditions, while joint optimisation of stiffness and hydraulic conductivity accurately reproduces observed pore pressure dissipation and settlement. The proposed framework advances predictive modelling of post-liquefaction re-solidification by providing a mechanistically interpretable representation of hydro-mechanical coupling beyond empirical fitting approaches.
partially drained conditions, while joint optimisation of stiffness and hydraulic conductivity accurately reproduces observed pore pressure dissipation and settlement. The proposed framework advances predictive modelling of post-liquefaction re-solidification by providing a mechanistically interpretable representation of hydro-mechanical coupling beyond empirical fitting approaches.
Date Issued
2026-10-01
Date Acceptance
2026-05-13
Citation
Computers and Geotechnics, 2026, 198
ISSN
0266-352X
Publisher
Elsevier
Journal / Book Title
Computers and Geotechnics
Volume
198
Copyright Statement
© 2026 The Author(s). Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
License URL
Publication Status
Published
Article Number
108262
Date Publish Online
2026-05-29
