岩土力学 ›› 2024, Vol. 45 ›› Issue (1): 302-312.doi: 10.16285/j.rsm.2023.0163

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

各向异性砂土渗流潜蚀流体动力学-离散元流固耦合分析

周闯1, 2, 3,钱建固1, 2,尹振宇3   

  1. 1. 同济大学 地下建筑与工程系,上海 200092;2. 同济大学 岩土及地下工程教育部重点实验室,上海 200092; 3. 香港理工大学 土木与环境系,香港
  • 收稿日期:2023-02-16 接受日期:2023-05-04 出版日期:2024-01-10 发布日期:2024-01-17
  • 通讯作者: 钱建固,男,1972年生,教授,博士生导师,主要从事岩土力学与本构理论、软土基础工程和交通岩土工程方面的研究。E-mail:qianjiangu@tongji.edu E-mail:1910309@tongji.edu.cn
  • 作者简介:周闯,男,1996年生,博士研究生,主要从事砂土渗流侵蚀方面的研究。
  • 基金资助:
    国家自然科学基金项目(No. 52178345)

Computational fluid dynamics-discrete element fluidsolid coupling analysis on suffusion in anisotropic sandy soils

ZHOU Chuang 1, 2, 3, QIAN Jian-gu 1, 2, YIN Zhen-yu3   

  1. 1. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China; 2. Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education, Tongji University, Shanghai 200092, China; 3. Department of Civil and Environmental Engineering, The Hong Kong Polytechnic University, Hong Kong, China
  • Received:2023-02-16 Accepted:2023-05-04 Online:2024-01-10 Published:2024-01-17
  • Supported by:
    This work was supported by the National the Natural Science Foundation of China (52178345).

摘要: 针对具有初始各向异性的间断级配砂土细粒潜蚀问题,引入可以考虑颗粒投影面积影响的Ganser拖曳力计算模型,实现非球形颗粒的计算流体动力学(computational fluid dynamics,简称CFD)和离散元(discrete element method,简称DEM)的双相耦合。通过与单颗粒下沉试验的对比,验证了该数值方法在解决异形颗粒与流体相互作用时的适用性。在此基础上,生成具有不同沉积方向和不同细粒含量的初始各向异性试样,模拟向上渗流潜蚀试验,并在试验中监测细粒流失量、强弱力链组成以及颗粒组构变化等宏微观特性,研究不同充填状态下(欠填充和过填充)不同组构各向异性土体渗流潜蚀特征。之后,对受潜蚀前后的试样进行了排水三轴试验,探究渗流对土体强度弱化的影响。结果表明,过填充试样质量损失随着颗粒沉积角度的增大而增大,而欠填充试样质量损失随沉积角度先增大后减小;欠填充试样细粒损失主要来源于低连通性细颗粒,而对于过填充试样,潜蚀则会导致低连通性和高连通性细颗粒数量同时减小。此外,三轴试验表明,潜蚀致土体峰值强度发生显著弱化,且峰值强度随沉积角度的变化也会受到土体充填状态的影响。

关键词: CFD-DEM, 渗流潜蚀, 异形颗粒 , 拖曳力模型, 各向异性, 微观组构

Abstract: For the problem of suffosion in gap-graded sand with initial anisotropy, the Ganser drag force model, which can take into account the effect of the projected area of particles, is introduced to achieve a two-phase coupling of computational fluid dynamics (CFD) and discrete elements method (DEM) for non-spherical particles. The applicability of the numerical method in solving the interaction between the non-spherical particles and fluid is verified by comparing with single particle settlement tests. On this basis, specimens with different bedding orientations and fine contents are further generated to simulate upward seepage suffosion tests, during which both macroscopic and microscopic properties, such as the fine loss, composition of strong and weak force chains, and changes in grain fabric, are monitored to explore the seepage suffosion characteristics of anisotropic soils with various fabrics under different filling states (underfilled and overfilled). Drained triaxial tests are carried out on specimens before and after erosion to investigate the effect of seepage on the weakening of soil strength. The results show that the mass loss of the overfilled specimens increases with increasing bedding angle, while the mass loss of the underfilled specimens firstly increases and then decreases with the bedding angle. The loss of fines in the underfilled specimens is mainly due to the low connectivity fines, whereas for the overfilled specimens, suffosion leads to a simultaneous reduction in the number of both low and high connectivity fines. In addition, the triaxial tests show that suffusion causes a significant weakening of the peak strength of the soil, and the change in peak strength with the bedding angle is also influenced by the soil filling state.

Key words: CFD-DEM, suffosion, non-spherical particles, drag force model, anisotropy, microscopic fabric

中图分类号: 

  • TU 43
[1] 张沛云, 木林隆, 黄茂松, . 基于计算流体力学-离散元耦合的相对密实度对断级配土体渗蚀效应的影响研究[J]. 岩土力学, 2024, 45(1): 267-283.
[2] 罗国立, 张科, 齐飞飞, 朱辉, 张凯, 刘享华, . 基于3D打印的裂隙岩体力学特性尺寸效应及各向异性初探[J]. 岩土力学, 2023, 44(S1): 107-116.
[3] 张革, 曹玲, 王成汤, . 考虑各向异性影响的冻土修正线性黏结接触模型开发及应用[J]. 岩土力学, 2023, 44(S1): 645-654.
[4] 余俊, 李东凯, 和振, 张志中. 带有两端防渗墙坝基的各向异性渗流解析解[J]. 岩土力学, 2023, 44(8): 2381-2388.
[5] 王智超, 彭乙芹, 秦云, 田英辉, 罗光财, . 基于统一屈服准则超固结土的应力诱导各向异性下负荷面模型[J]. 岩土力学, 2023, 44(7): 1891-1900.
[6] 王伟, 张宽, 曹亚军, 陈超, 朱其志, . 层状千枚岩各向异性力学特性与脆性评价研究[J]. 岩土力学, 2023, 44(4): 975-989.
[7] 黄娟, 和振, 余俊, 杨鑫歆. 渗透各向异性圆形围堰稳态渗流场解析解[J]. 岩土力学, 2023, 44(4): 1035-1043.
[8] 陈曦. 基于方向性粗糙度参数的节理峰值抗剪 强度理论模型[J]. 岩土力学, 2023, 44(4): 1075-1088.
[9] 梁晓敏, 杨朔成, 顾晓强, . 砂土应力诱发弹性波速各向异性的试验研究[J]. 岩土力学, 2023, 44(11): 3235-3240.
[10] 张涛, 徐卫亚, 孟庆祥, 王环玲, 闫龙, 钱坤, . 基于3D打印技术的柱状节理岩体试样 力学特性试验研究[J]. 岩土力学, 2022, 43(S2): 245-254.
[11] 包含, 陈志洋, 兰恒星, 裴润生, 伍法权, 晏长根, 陶悦, . 矿物定向排列致各向异性岩石渐进破坏强度特征 ——以黑云母石英片岩为例[J]. 岩土力学, 2022, 43(8): 2060-2070.
[12] 蒋中明, 肖喆臻, 唐栋, 何国富, 许卫, . 基于裂隙渗流效应的水封油库涌水量预测分析[J]. 岩土力学, 2022, 43(4): 1041-1047.
[13] 张虎元, 王赵明, 朱江鸿, 周光平, . 混合型缓冲材料砌块渗透性及其各向异性研究[J]. 岩土力学, 2022, 43(3): 573-581.
[14] 刘洋, 于鹏强, 徐硕. 基于微形态理论的各向异性散体波动分析[J]. 岩土力学, 2022, 43(3): 635-648.
[15] 潘文韬, 杨文波, 吴枋胤, 何川, 赵亮亮, 姚人杰, 傅舰锋, . 基于单三轴试验与直剪试验的层状软岩模拟[J]. 岩土力学, 2022, 43(12): 3437-3452.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 姚仰平,侯 伟. 土的基本力学特性及其弹塑性描述[J]. , 2009, 30(10): 2881 -2902 .
[2] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[3] 向天兵,冯夏庭,陈炳瑞,江 权,张传庆. 三向应力状态下单结构面岩石试样破坏机制与真三轴试验研究[J]. , 2009, 30(10): 2908 -2916 .
[4] 夏栋舟,何益斌,刘建华. 土-结构动力相互作用体系阻尼及地震反应分析[J]. , 2009, 30(10): 2923 -2928 .
[5] 张其一. 复合加载模式下地基失效机制研究[J]. , 2009, 30(10): 2940 -2944 .
[6] 陶干强,杨仕教,任凤玉. 崩落矿岩散粒体流动性能试验研究[J]. , 2009, 30(10): 2950 -2954 .
[7] 黄润秋,徐德敏. 岩石(体)渗透性测试的体变量法研究[J]. , 2009, 30(10): 2961 -2964 .
[8] 卢 正,姚海林,骆行文,胡梦玲. 公路交通荷载作用下分层地基的三维动响应分析[J]. , 2009, 30(10): 2965 -2970 .
[9] 陈 阵,陶龙光,李 涛,李海斌,王综勇. 支护结构作用的箱基沉降计算新方法[J]. , 2009, 30(10): 2978 -2984 .
[10] 王淑云,鲁晓兵,赵 京,王爱兰. 粉质黏土周期荷载后的不排水强度衰化特性[J]. , 2009, 30(10): 2991 -2995 .