基于CFD-DEM的级配缺失砂土渗流破坏规律研究.pdfVIP

第43卷第4期佛山科学技术学院学报（自然科学版）Vol.43No.4

2025年7月JournalofFoshanUniversity(NaturalSciencesEdition)Jul.2025

文章编号：1008-0171(2025)04-0008-06

基于CFD-DEM的级配缺失砂土渗流破坏规律研究

郭家仪1,程马遥1*，张煜雄1,黄碧琦1,陈林生2

（1.佛山大学土木与交通学院，广东佛山528225；

2.广州肖宁道路工程技术研究事务所有限公司，广东广州511452）

摘要：为探讨级配缺失砂土渗透破坏的发生发展机制，基于CFD-DEM建立渗透数值模型，对不同细颗粒含量

和渗流速度组合下的缺级配砂土进行了数值试验，分析了渗流区域、渗流路径和接触力链的变化情况。研究结

果表明：在渗流初始阶段，细小颗粒开始发生转动和微移；随着渗透力的增大，细颗粒呈抱团式运移；渗透破坏

后期，细颗粒通过孔隙从上游移至下游，导致土体呈现镂空状态。适当增加细颗粒含量可降低渗透破坏风险，

增强接触力链的连贯性，但细颗粒含量过高会发生过度堆积而导致下游堵塞。渗流速度对细颗粒迁移路径也

会产生影响。在低流速条件下，细颗粒主要沿着中间区域的渗流通道进行迁移，而当流速升高时，细颗粒首先

向两侧扩散，随后再沿中间的渗流通道进行移动。

关键词：级配缺失砂土；渗透破坏；CFD-DEM；数值模拟

中图分类号：TU441文献标志码：A

地下工程在水力作用下的渗透破坏具有频在传统试验研究方面。通过砂土管涌试验

发、隐秘、复杂、难以监测和处理等特点。渗透破发现，细颗粒含量会影响渗透破坏及发展过程，

坏的形式多种多样，一般包括潜蚀、管涌、流砂。细颗粒含量与临界水力梯度和临界流速之间存

对于级配不均特别是级配空白的砂土，典型渗透在负相关关系。张冬梅等[2]基于侵蚀试验中的水

破坏现象为潜蚀及其发展而成的管涌，土体内较土流失变化，将间断级配砂土中侵蚀过程分为快

细的颗粒通过较大颗粒所形成的孔道移动，其导速侵蚀阶段、稳定侵蚀阶段和收敛阶段，并总结

致的结果是减弱土体结构性，使得土骨架镂空，每个侵蚀阶段细颗粒损失量及渗透率等参数的

从而形成渗流通道。目前，针对级配缺失砂土渗阈值规律。朱秦等[3]结合细颗粒迁移情况，建立了

透破坏机制和规律的研究主要采取传统试验方渗透系数变化与土体细颗粒迁移行为的对应关

法和数值模拟方法。传统试验方法主要包括室内系，并提出3种宽级配土的细颗粒迁移模式。田

试验和现场试验，揭示砂土的力学特性和破坏模大浪等[4]利用自主研发的刚性壁渗透仪,对不同

式；而数值模拟方法则通过细观分析土体内部的细颗粒含量的间断级配砂砾石土在恒定水头渗

变形和破坏过程，进而预测砂土的力学行为和破流作用下进行全过程渗透变形试验，分析了渗透

坏特征。试验结束后土体的颗粒级配空间分布变化。

收稿日期：2024-10-15

基金项目：广东省自然科学基金地区培育项目（2024A1515140052）；广东省住房和城乡建设厅科技创新计划项目（2023-

K33-071697)

作者简介：郭家仪（1998-），女，广东广州人，佛山大学硕士研究生。

通信作者：程马遥（1988-），女，江苏宿迁人，佛山大学讲师，博士。

第4期郭家仪等：基于CFD-DEM的级配缺失砂土渗流破坏规律研究9

在数值模拟研究方面。刘先珊[5]基于水压力

级配缺失砂土渗透破坏的数值模拟

2

梯度力和拖曳力的流固耦合共同作用，分析了单

向颗粒的运动模式，并与理论结果相比较。2.1级配缺失砂土模型

Nguyen[6采用离散元法建立不同级配特性和孔隙根据WAN[12]提出的土体内部稳定与不稳定

率的模型,并结合Navier-Stokes方程的计算