岩土力学 ›› 2024, Vol. 45 ›› Issue (8): 2492-2501.doi: 10.16285/j.rsm.2023.1461

• 数值分析 • 上一篇    下一篇

基于CFD-DEM方法的饱和砂土场地液化模拟研究

许文昊1, 2,王志华1, 2,申志福1, 2,高洪梅1, 2,刘殷强1, 2,张鑫磊1, 2   

  1. 1. 南京工业大学 岩土工程研究所,江苏 南京 211816;2. 河海大学 江苏省交通基础设施安全保障技术工程研究中心,江苏 南京 211816
  • 收稿日期:2023-09-26 接受日期:2023-12-21 出版日期:2024-08-10 发布日期:2024-08-12
  • 通讯作者: 高洪梅,女,1982年生,博士,教授,主要从事土动力学与地基处理方面的研究和教学工作。E-mail: hongmei54@163.com
  • 作者简介:许文昊,男,1997年生,博士研究生,主要从事土动力学方面的研究工作。E-mail: yzxwh1997@163.com
  • 基金资助:
    国家自然科学基金面上项目(No. 51678300,No. 52178336,No. 52108324);江苏省高等学校自然科学研究重大项目(No. 18KJA560002);江苏省高校“青蓝工程”中青年学术带头人;江苏省研究生科研与实践创新计划项目(No. KYCX22_1354)。

Simulation of saturated sand site liquefaction based on the CFD-DEM method

XU Wen-hao1, 2, WANG Zhi-hua1, 2, SHEN Zhi-fu1, 2, GAO Hong-mei1, 2, LIU Yin-qiang1, 2, ZHANG Xin-lei1, 2   

  1. 1. Institute of Geotechnical Engineering, Nanjing Tech University, Nanjing, Jiangsu 211816, China; 2. Jiangsu Province Engineering Research Center of Transportation Infrastructure Security Technology, Hohai University, Nanjing, Jiangsu 211816, China
  • Received:2023-09-26 Accepted:2023-12-21 Online:2024-08-10 Published:2024-08-12
  • Supported by:
    This work was supported by the General Program of National Natural Science Foundation of China (51678300, 52178336, 52108324), the Major Research Project in Natural Sciences for Higher Education Institutions of Jiangsu Province (18KJA560002), the Young Academic Leader in the “Qinglan Project” of Universities in Jiangsu Province and the Jiangsu Provincial Graduate Research and Practice Innovation Program (KYCX22_1354).

摘要: 砂土液化是常见的地震灾害,目前应用于研究砂土液化动力特性的室内试验以及模型试验还不能全面反映土体液化全过程。计算流体动力学(computational fluid dynamics,CFD)与离散元法(discrete element method,DEM)耦合模拟方法能够准确地模拟各类水土耦合问题。通过二次开发的CFD-DEM流固耦合模块实现离散元软件PFC3D与计算流体力学软件OpenFOAM之间的力学信息交互,利用颗粒水下自由沉降验证该方法的可行性。利用PFC3D软件模拟室内循环三轴试验标定出具有真实饱和砂土动力特性的数值砂样。根据已有的参数信息以及耦合模拟方法建立了饱和砂土的场地液化模型。模拟结果表明,离散元法能够复现室内砂土液化试验,标定参数可应用于场地液化模拟;单颗粒沉降速度与理论解一致验证了CFD-DEM耦合方法的准确性;峰值加速度0.25g下不同深度处土体均会发生液化,液化时超孔压比无法达到1,超孔压累计值由浅层往深层递增;液化后土体强度自下而上逐渐恢复,再固结的场地土体结构呈现均匀化发展趋势。

关键词: 砂土液化, CFD-DEM耦合方法, 颗粒自由沉降, 场地液化模拟

Abstract: Soil liquefaction is a prevalent seismic hazard. However, current indoor and model experiments studying the dynamic characteristics of sand liquefaction struggle to accurately represent the actual soil liquefaction process. The computational fluid dynamics (CFD) coupled with discrete element method (DEM) simulation method can accurately simulate various soil-water coupling problems. The CFD-DEM flow-solid coupling module facilitated the exchange of mechanical information between the commercial discrete element software PFC3D and the open-source computational fluid dynamics software OpenFOAM. The feasibility of this approach was confirmed through particle underwater free settling experiments. Calibration of numerical sand specimens with dynamic characteristics of real saturated sand was conducted using PFC3D software through simulated indoor cyclic triaxial tests. Based on the existing parameter information and the coupled simulation method, a site liquefaction model of saturated sand was established. The simulation results indicate that the discrete element method (DEM) can replicate indoor sand liquefaction experiments, and the calibrated parameters can be applied to site liquefaction simulations. The consistency between the settling velocity of individual particles and theoretical solutions validates the accuracy of the CFD-DEM coupling method. Under a peak acceleration of 0.25g, liquefaction occurs at various depths, and the ratio of excess pore pressure does not exceed 1 during liquefaction. The cumulative excess pore pressure increases from shallow layers to deep layers. After liquefaction, the soil strength gradually recovers from bottom to top, and the soil structure in the re-consolidated site shows a trend of homogenization.

Key words: soil liquefaction, CFD-DEM coupling method, particle free settling, site liquefaction simulation

中图分类号: TU 43
[1] 李金宇, 王伟, 王浩宇, 杨研科, 徐凯放, 张晓庆, 熊文. 地震序列作用下液化场地条件影响研究[J]. 岩土力学, 2024, 45(5): 1551-1559.
[2] 王晓磊, 刘理腾, 刘润, 刘历波, 董林, 任海. 地震历史对各深度土体抗液化性影响的振动台试验研究[J]. 岩土力学, 2023, 44(9): 2657-2666.
[3] 吴宏, 叶治, 张宇亭, 刘华北, . 穿越不同密实度饱和砂土地层的盾构隧道地震响应三维数值分析[J]. 岩土力学, 2023, 44(4): 1204-1216.
[4] 王维铭, 陈龙伟, 郭婷婷, 汪云龙, 凌贤长, . 基于中国砂土液化数据库的标准贯入试验液化 判别方法研究[J]. 岩土力学, 2023, 44(1): 279-288.
[5] 董青, 周正华, 苏杰, 李小军, 郝冰, . 基于对数动骨架的可考虑液化大变形本构[J]. 岩土力学, 2021, 42(7): 1903-1910.
[6] 庄海洋, 付继赛, 陈 苏, 陈国兴, 王雪剑, . 微倾斜场地中地铁地下结构周围地基液化与变形特性振动台模型试验研究[J]. 岩土力学, 2019, 40(4): 1263-1272.
[7] 魏 星, 张 昭, 王 刚, 张建民, . 饱和砂土液化后大变形机制的离散元细观分析[J]. 岩土力学, 2019, 40(4): 1596-1602.
[8] 付海清,袁晓铭,王 淼,. 基于现场液化试验的饱和砂土孔压增量计算模型[J]. , 2018, 39(5): 1611-1618.
[9] 宫凤强,李嘉维,. 基于PCA-DDA原理的砂土液化预测模型及应用[J]. , 2016, 37(S1): 448-454.
[10] 王良民,叶剑红,朱长歧. 近海欠密实砂质海床内波致渐进液化特征研究[J]. , 2015, 36(12): 3583-3588.
[11] 陈育民 ,高 星 ,刘汉龙,. 砂土液化流动变形的简化方法[J]. , 2013, 34(6): 1567-1573.
[12] 方火浪 ,张轶群 ,郭 婧 ,银 鸽 . 深厚砂质覆盖层土坝的弹塑性地震反应分析[J]. , 2013, 34(11): 3197-3204.
[13] 刘年平 ,王宏图 ,袁志刚 ,刘竟成. 砂土液化预测的Fisher判别模型及应用[J]. , 2012, 33(2): 554-557.
[14] 庄海洋,黄春霞,左玉峰. 某砂土液化大变形本构模型参数的敏感性分析[J]. , 2012, 33(1): 280-286.
[15] 侯龙清,徐红梅,曹振中,袁晓铭. 汶川地震液化土层类型验证及土性分析[J]. , 2011, 32(4): 1119-1124.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!