关于土渗透系数颗粒流细观参数研究.pdfVIP

低 温 建 筑 技 术 2012年第 1期(总第 163期) 关于土渗透系数颗粒流细观参数研究 白若虚 ，程雪松 ，郑 刚 (1．滨海土木工程结构与安全教育部重点实验室。 天津 300072；2．天津大学土木工程系。 天津 300072) 【摘 要】 利用颗粒流程序PFC2D对地下水的渗流作用进行了模拟，给出了颗粒流数值模拟确定渗透系数 的方法，并对与渗透系数有关的细观参数进行了参数分析 ，并与渗透系数建立了关系。文中结论为颗粒流理论模 拟地下水渗流相关的岩土工程问题提供了依据和理论基础。 【关键词】 颗粒流；渗流；达西定律；参数分析 【中图分类号】 TU441．33 文【献标识码】 B 【文章编号】 1001—6864(2012)01—0o64—04 目前对于渗流耦合作用的数值模拟主要采用有限单元 平均半径。 法和有限差分法 ，但都有其局限性，即只能模拟小变形下简 2 模型建立 单边界的连续变形问题。国内外学者 ¨“]利用颗粒流PFC 本文将砂土颗粒简化成圆形颗粒，故可直接采用PFC。 (ParticleFlowCode)建立数值模型模拟砂土的渗透作用，证 进行模拟。模型尺寸为20mmx20mm，为了模拟土颗粒的不均 明了达西定律的适用性。但是尚未有学者对颗粒流模拟流 匀性，颗粒单元半径 采用从R 到R～ 的均匀分布，根据计 固耦合中的相关细观参数如何确定做出研究，而这些细观 算效率和计算精度要求，取R =0．35mm，R～ =0．50mm。 参数在模拟研究土体渗流破坏、水力劈裂、劈裂注浆和基坑 先按照模型尺寸生成四面墙，待系统平衡后再将其删去，并固 突涌等问题时都非常关键。由于颗粒流不能直接给出模型 定上下边界颗粒。然后调用程序 内置 FISH语言 dom．fis和 介质的物理力学参数 J，土的宏观参数只能靠调整基本粒 dom1．fis生成 “域”和 “管道”所组成的网络。模型示意图见图 子的微观参数获得。本文即从细观层面上对土的渗流耦合 1，浅蓝色小球为土颗粒，黑色小球为 “域”，“域”之间的黑色 进行数值模拟，试图建立细观参数与渗透系数之间的关系。 通道为 “管道”，白色线代表两颗粒间的接触。 1 水土耦合的颗粒流模拟 在PFC。中采用圆形颗粒代表土颗粒，进行流固耦合模 拟时其实并没有实实在在的流体存在，因此引入了 “域”的 概念。 “域”就是一系列封闭的颗粒链，在链上的每个链接 都是一个连接接触。“域”相当于颗粒间空隙的 “水库”，“水 库”可以储存水压，各水库之间存在通道 ，这些通道用在颗 粒接触处相切于两个颗粒的 “管道”来模拟。 一圈 PFC。定义的用于流固耦合计算的流动方程和压力方 程以及求解条件如下 。 图1 模型示意图 图2 压力圆分布示意图 (1) 流动方程。相邻域间的每一个链接都看作是一 条潜在的水力通道。通道长度为 ￡、孔径为o，在垂直平面方 在模型左侧施加水头压力 (均取3X10Pa)，右侧水头 向上为单位厚度，管道内的流速(单位时间内的体积)为： 压力固定为零，则在模型左右两侧形成水头压力差。水头压 q=Ka(P2一P1)儿 (1) 力加压范 围左右两侧各取 1．5mm，故实际渗流路径为 式中，为传导系数；L为周围颗粒半径和；(P 一p。)为 17mm。选取模型中段 10mm范围监控其计算过程 中所有 “域”的总流量Q。循环稳定后压力分布如图2所示。图中白 两个相邻域的压力差。 色圆点大小代表 “域”内压力大小，圆的直径大小代表压力 (2) 压力方程。周围管道流