Implementation of A Liquefaction Subroutine in ABAQUS Based on The Modified Hardin-drnevich Soil Backbone Curve Model
DOI:
https://doi.org/10.54691/nqfagb50Keywords:
Numerical simulation; Byrne incremental model; Modified Hardin-Drnevich; Liquefaction.Abstract
Liquefaction caused by earthquakes can affect the stability of underground structures. Numerical simulation of liquefaction is therefore necessary. By combining the modified Hardin-Drnevich soil constitutive model with Byrne's modified Martin pore pressure increment model, this paper develops an ABAQUS subroutine incorporating the modified Hardin-Drnevich soil constitutive model and the Byrne pore pressure increment model based on the implicit interface of ABAQUS. Numerical models in ABAQUS and FLAC3D are established respectively based on the undrained dynamic triaxial test results of Nanjing saturated fine sand, and the effective analytical capability of the proposed subroutine is preliminarily verified. Through comparing a total of 6 groups of numerical calculation results under different working conditions, this study verifies the influences of pore pressure growth trend, initial effective pore pressure and dynamic shear ratio on the liquefaction process. It is proved that the calculation results of the proposed subroutine in dynamic triaxial tests are superior to those of FLAC3D. The causes of numerical calculation errors are discussed, and the feasibility of developing the liquefaction numerical simulation subroutine for ABAQUS is verified.
Downloads
References
[1] Shen, Q., et al. (2015). Study on scheme for underwater shield tunnel through sand liquefaction zone in strong earthquake load. Highway Engineering, 40(2), 54–58.
[2] Gao, L. L., et al. (2019). Study on the impact of geological hazard sources and measures on survey of utility tunnel in Tianjin urban area. Mineral Exploration, 10(12), 3023–3032.
[3] Li, X. X., et al. (2024). Research on liquefaction-induced lateral spreading in the valley site of the utility tunnel crossing the Yangzi River. Journal of Disaster Prevention and Mitigation Engineering, 44(3), 670–677.
[4] Zhang, L. F., et al. (2015). Prospects and study on development of underground space in Hangzhou, China. Tunnel Construction, 35(4), 285–291.
[5] Wang, L., et al. (2016). The implementation of dynamic effective stress method based on Matasovic constitutive model in ABAQUS software. Journal of Disaster Prevention and Mitigation Engineering, 36(2), 231–238.
[6] Wang, G., & Zhang, J. M. (2007). Numerical modeling of liquefaction-induced deformation in sand. Chinese Journal of Geotechnical Engineering, 29(3), 403–409.
[7] Mao, H. Y., & Le, T. (2012). Numerical simulation method for seismic load action on tailings dams. Modern Mining, 27(5), 65–66.
[8] He, J. P., & Chen, W. Z. (2011). Numerical experiments and analysis for liquefaction in a non-free-field with gravel stone drainage layer around underground structures. Journal of Shandong University (Engineering Science), 41(3), 93–100.
[9] Xia, Z. F., et al. (2010). Numerical analysis of effect of waterproof measures on earthquake induced liquefaction of embankment. Rock and Soil Mechanics, 31(8), 2682–2688. https://doi.org/10.16285/j.rsm.2010.08.032
[10] Martin, G. B., Finn, W. D. L., & Seed, H. B. (1975). Fundamentals of liquefaction under cyclic loading. Journal of the Geotechnical Engineering Division, 101(GT5), 423–438. https://doi.org/10.1061/AJGEB6.0000164
[11] Hashash, Y. M. A., & Park, D. (2001). Non-linear one-dimensional seismic ground motion propagation in the Mississippi embayment. Engineering Geology, 62(1–3), 185–206. https://doi.org/10.1016/S0013-7952(01)00061-3
[12] Tokimatsu, K., & Seed, H. B. (1987). Evaluation of settlements in sands due to earthquake shaking. Journal of Geotechnical Engineering, 113(8), 861–878. https://doi.org/10.1061/(ASCE)0733-9410(1987)113:8(861)
[13] Byrne, P. M. (1991). A cyclic shear-volume coupling and pore pressure model for sand. In Proceedings of the 2nd International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics (pp. 47–55). Missouri, USA.
[14] Liao, H. J., et al. (2009). Comparative analysis of backbone curve models for loess soils. Rock and Soil Mechanics, 30(Suppl. 2), 17–21.
[15] Liyanathirana, D. S., & Poulos, H. G. (2002). Numerical simulation of soil liquefaction due to earthquake loading. Soil Dynamics and Earthquake Engineering, 22(7), 511–523. https://doi.org/10.1016/S0267-7261(02)00022-1
[16] Zhang, X. L., et al. (2023). Property evolution and liquefaction stage characteristics of saturated sand under cyclic loading. Engineering Mechanics, 40(2), 157–167. https://doi.org/10.6052/j.issn.1000-4750.2021.1632
[17] Ruan, B., et al. (2017). A new effective stress method using the integration algorithm of the modified Davidenkov constitutive model and the Byrne pore pressure increment model, and validation and calibration. Journal of Basic Science and Engineering, 25(5), 956–966. https://doi.org/10.16058/j.issn.1005-0930.2017.05.013
Downloads
Published
Issue
Section
License
Copyright (c) 2026 Frontiers in Science and Engineering

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.






