Citation

BibTex format

@article{MA:2026:10.1016/j.compgeo.2026.108262,
author = {MA, S and Taborda, D and Kontoe, S},
doi = {10.1016/j.compgeo.2026.108262},
journal = {Computers and Geotechnics},
title = {Numerical framework for calculating stiffness of liquefied sands during post-shaking re-solidification},
url = {http://dx.doi.org/10.1016/j.compgeo.2026.108262},
volume = {198},
year = {2026}
}

RIS format (EndNote, RefMan)

TY  - JOUR
AB - 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 hy
AU - MA,S
AU - Taborda,D
AU - Kontoe,S
DO - 10.1016/j.compgeo.2026.108262
PY - 2026///
SN - 0266-352X
TI - Numerical framework for calculating stiffness of liquefied sands during post-shaking re-solidification
T2 - Computers and Geotechnics
UR - http://dx.doi.org/10.1016/j.compgeo.2026.108262
UR - https://www.sciencedirect.com/science/article/pii/S0266352X2600368X?via%3Dihub
VL - 198
ER -

Contact Geotechnics

Geotechnics
Civil and Environmental Engineering
Skempton Building
51³Ô¹ÏÍø
South Kensington Campus
London, SW7 2AZ

Telephone:
+44 (0)20 7594 6077
Email: geotechnics@imperial.ac.uk 
Alternatively, you can find a member of Geotechnics staff on the Department of Civil and Environmental Engineering website

Follow us on Twitter: 

We are located in the Skempton Building (building number 27 on the South Kensington Campus Map). How to find us