基于离散元分析的各向异性砂土循环荷载下的宏微观液化特性

展开
  • (1. 中海石油(中国)有限公司北京研究中心,北京 100028; 2. 中国海洋大学 海底建设与保护山东省工程研究中心,山东 青岛 450007; 3. 大连海事大学 交通运输工程学院,辽宁 大连 116026)
李书兆(1985 — ),女,博士,副研究员,E-mail:lishuzhao_ouc@163.com;崔春义*(1978 — ),男,博士,教授,博士生导师, E-mail:cuichunyi@dlmu.edu.cn。

网络出版日期: 2025-01-04

基金资助

国家自然科学基金资助项目(42302304;42025702);山东省自然科学基金资助项目(ZR2022QD017)

Macroscopic and microscopic liquefaction characteristics of anisotropic sand under cyclic loading based on discrete element analysis

Expand
  •  (1. Beijing Research Center, China National Offshore Oil Co., Ltd., Beijing 100028, China; 2. Shandong Engineering Research Center of Marine Exploration and Conservation, Ocean University of China, Qingdao 266100, China; 3. College of Transportation Engineering, Dalian Maritime University, Dalian 116026, China)

Online published: 2025-01-04

摘要

利用离散元法分析初始各向异性砂土在非对称循环荷载下的液化行为。通过生成椭球颗粒制备各向异性砂土试样,采用常体积法开展不排水循环三轴剪切模拟,分析循环加载过程中力学配位数、冗余指数和组构各向异性的演化,揭示砂土在微观层面的液化特征。结果表明,离散元模拟有效反映了非对称循环加载条件下密砂更高的抗液化能力。试样在完全液化时,力学配位数和冗余指数分别降至0.35和2.4,这与是否对称循环加载无关。不同加载条件下,颗粒长轴的组构随着循环加载由各向异性向各向同性转变。在应力反转条件下,液化发生后,颗粒长轴组构各向异性迅速减小,同时,接触法向组构发生显著变化并出现方向切换。初始液化均发生在接触法向组构各向异性变量从负值逐渐接近0的过程中。研究结果有助于深入理解液化的微观力学机制,有助于构建复杂加载条件下各向异性砂土循环液化的宏观本构模型。

本文引用格式

李书兆, 孙国栋, 张安, 王栋, 崔春义 . 基于离散元分析的各向异性砂土循环荷载下的宏微观液化特性[J]. 大连海事大学学报, 2025 , 51(2) : 115 -124 . DOI: 10.16411/j.cnki.issn1006-7736.2025.02.013

Abstract

The discrete element method (DEM) was used to analyze the liquefaction behavior of initially anisotropic sand under nonsymmetric cyclic loading. Ellipsoidal particles were generated to prepare anisotropic sand samples, and undrained cyclic triaxial shear simulations were conducted by using constant volume method. The evolution of mechanical coordination number, redundancy index and fabric anisotropy during cyclic loading was studied, and the liquefaction characteristics of sand at the microscopic level was revealed. The results show that the DEM simulations effectively reflect that dense sand exhibits higher liquefaction resistance under nonsymmetric cyclic loading conditions. When sand undergoes full liquefaction, the mechanical coordination number and redundancy index decrease to 0.35 and 2.4, respectively, regardless of whether the cyclic loading is symmetrical. Under different loading conditions, the particle orientation fabric transitions from anisotropic to isotropic with cyclic loading. Under stress reversal conditions, once liquefaction occurs, the anisotropy of the particle orientation fabric rapidly decreases, while the contact normal fabric changes significantly and shows a directional switch. Initial liquefaction occurs as the contactnormal fabric anisotropy variable transitions from a negative value to near zero. These findings provide an in-depth understanding of the micromechanical mechanisms of liquefaction, contributing to the development of macroscopic constitutive models for cyclic liquefaction of anisotropic sands under complex loading conditions. 


参考文献

[1]DE GROOT M B, KUDELLA M, MEIJERS P, et al. Liquefaction phenomena underneath marine gravity structures subjected to wave loads[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 2006, 132(4): 325-335.
[2]LIU T, WEI G L, KOU H L, et al. Pore pressure observation: pressure response of probe penetration and tides[J]. Acta Oceanologica Sinica, 2019, 38(7): 107-113.
[3]HYODO M, MURATA H, YASUFUKU N, et al. Undrained cyclic shear strength and residual shear strain of saturated sand by cyclic triaxial tests[J]. Soils and Foundations, 1991, 31(3): 60-76.
[4]YANG Z X, PAN K. Flow deformation and cyclic resistance of saturated loose sand considering initial static shear effect[J]. Soil Dynamics and Earthquake Engineering, 2017, 92: 68-78.
[5]SIVATHAYALAN S, HA D. Effect of static shear stress on the cyclic resistance of sands in simple shear loading[J]. Canadian Geotechnical Journal, 2011, 48(10): 1471-1484.
[6]YANG J, SZE H Y. Cyclic behaviour and resistance of saturated sand under non-symmetrical loading conditions[J]. Géotechnique, 2011, 61(1): 59-73.
[7]PAN K, ZHOU G Y, YANG Z X, et al. Comparison of cyclic liquefaction behavior of clean and silty sands considering static shear effect[J]. Soil Dynamics and Earthquake Engineering, 2020, 139: 106338.
[8]LIU X, QIN Z H, YANG J. Undrained cyclic behaviour of loess with initial shear stress: a focus on failure mode[J]. Soil Dynamics and Earthquake Engineering, 2023, 171: 107971.
[9]CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Géotechnique, 1979, 29(1): 47-65.
[10]YANG S Y, HUANG D R. Understanding fabric evolution and multiple liquefaction resistance of sands in the presence of initial static shear stress[J]. Soil Dynamics and Earthquake Engineering, 2023, 171: 107962.
[11]WEI J T, XU T J, HE J X. Effect of static shear stress and cyclic loading direction on cyclic behaviors of granular soils by DEM analysis[J]. Computers and Geotechnics, 2024, 167: 106112.
[12]WU Q X, PAN K, YANG Z X. Effects of initial static shear on undrained cyclic behavior of granular materials: energy evolution and micromechanical interpretation[J]. Granular Matter, 2023, 25(1): 7.
[13]ZHOU X H, ZHOU Y G, CHEN Y M. DEM analysis of undrained cyclic behavior and resistance of saturated dense sand without stress reversals[J]. Granular Matter, 2023, 25(2): 36.
[14]SZE H Y, YANG J. Failure modes of sand in undrained cyclic loading: impact of sample preparation[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2014, 140(1): 152-169.
[15]ZHANG A, JIANG M J, WANG D. Effect of fabric anisotropy on the cyclic liquefaction of sands: insight from DEM simulations[J]. Computers and Geotechnics, 2023, 155: 105188.
[16]ITASCA CONSULTING GROUP. Particle Flow Code in Three Dimensions: User's Manual, Version 5.0[M]. Minneapolis, USA: 2015.
[17]JIANG M J, SHEN Z F, WANG J F. A novel three-dimensional contact model for granulates incorporating rolling and twisting resistances[J]. Computers and Geotechnics, 2015, 65: 147-163.
[18]JIANG M J, KONRAD J M, LEROUEIL S. An efficient technique for generating homogeneous specimens for DEM studies[J]. Computers and Geotechnics, 2003, 30(7): 579-597.
[19]DA CRUZ F, EMAM S, PROCHNOW M, et al. Rheophysics of dense granular materials: discrete simulation of plane shear flows[J]. Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 2005, 72(2): 021309.
[20]YANG M, TAIEBAT M, RADJAÏ F. Liquefaction of granular materials in constant-volume cyclic shearing: transition between solid-like and fluid-like states[J]. Computers and Geotechnics, 2022, 148: 104800.
[21]JIANG M J, ZHANG A, FU C. 3-D DEM simulations of drained triaxial tests on inherently anisotropic granulates[J]. European Journal of Environmental and Civil Engineering, 2018, 22(sup.1): s37-s56.
[22]SEED R B, HARDER L F, Jr. SPT-based analysis of cyclic pore pressure generation and undrained residual strength[C]//Proceedings of the H. Bolton Seed Memorial Symposium. Vancouver, B. C., Canada: BiTech Publishing, 1990, 2: 351-376.
[23]HARDER JR, L. F., BOULANGER, R. W. Application of Kσ and Kα correction factors[C] // Proceedings of the NCEER Workshop on Evaluation of Liquefaction Resistance of Soils. Buffalo, New York: National Center for Earthquake Engineering Research, State University of New York at Buffalo, 1997: 167–190.
[24]THORNTON C. Numerical simulations of deviatoric shear deformation of granular media[J]. Géotechnique, 2000, 50(1): 43-53.
[25]SATAKE M. Fabric tensor in granular materials[C]// Proceedings of the IUTAM Symposium on Deformation and Failure of Granular Materials. Amsterdam: A.A. Balkema, 1982: 63-68.
[26]LI X S, DAFALIAS Y F. Anisotropic critical state theory: role of fabric[J]. Journal of Engineering Mechanics, 2012, 138(3): 263-275.
[27]HU N, YU H S, YANG D S, et al. Constitutive modelling of granular materials using a contact normal-based fabric tensor[J]. Acta Geotechnica, 2020, 15(5): 1125-1151.
[28]WANG R, CAO W, XUE L, et al. An anisotropic plasticity model incorporating fabric evolution for monotonic and cyclic behavior of sand[J]. Acta Geotechnica, 2021, 16(1): 43-65.

文章导航

/