基于离散元理论的砂土管涌数值模拟.docxVIP

基于离散元理论的砂土管涌数值模拟 1 基料-滤层系统的数值模型 管理河流已成为防止水库基础、河流廊道、基础等工程失衡的重要因素。在渗流过程中水土的相互作用贯穿管涌发展的全过程。管涌的渗透破坏具有复杂性、随机性、不易觉察性、且破坏的严重性已引起人们的重视。 如何有效的防治管涌是工程界所关注的,自1922年太沙基提出著名的反滤层防止管涌的理论后,砂土管涌的基料-滤层防治措施很快被工程界接受并应用。鉴于基料-滤层防治的复杂性,而关于该系统机理特性的理论研究很少。其中具有代表性的有:Payatakes把滤层孔隙中的渗流简化为解决一系列管径大小随机分布的管中流动问题,通过模型可模拟多孔介质渗流过程中的颗粒运动规律;Aberg基于颗粒和孔隙几何特性,建立数学模型描述颗粒的冲蚀和土体自滤现象;Indraratna基于极限平衡状态的临界水力梯度理论建立数学模型描述基料-滤层系统中的反滤现象;Lakshmi根据水流在筒状管中流动的基本理论,建立基料-滤层系统数值模型,分析由于流失颗粒滞留于滤层中而引起渗透系数减小现象;Mark Locke根据颗粒在不同水力梯度下的运动特性,提出描述粒状材料反滤的数值模型。 然而,由渗流引起的管涌现象是个复杂的流固耦合过程,涉及固体大变形和颗粒介质的流动理论。而当前对滤层的研究大多停留在宏观角度,通过建立颗粒移动和滞留数学模型对滤层特性进行研究,未能从细观角度深入探讨渗流过程中颗粒运动引起基料-滤层系统中孔隙率和渗透系数的变化,未能考虑到颗粒的运动引起整个土体结构和水力特性的变化,也就无法从本质上揭示滤层的工作机理。 本文采用离散元方法对渗流作用下的砂土管涌的基料-滤层系统进行数值模拟,分析不同层间系数(D15/d85s)对管涌发展和滤层有效性的影响;从细观角度揭示管涌发展过程中颗粒的运动特性和滤层防治机理,为管涌防治机理的更深入研究奠定理论基础。 2 基于pfc3的精度改进 砂土由一系列不连续的土颗粒构成,在离散单元法模型中,采用球颗粒来模拟砂土颗粒介质,采用圆颗粒主要是基于现有计算机运算速度和模型理想化考虑,进一步工作可通过PFC3D内置FISH语言开发非圆颗粒(CLUMP)来模拟砂土。颗粒间的作用采用接触本构模型来模拟。 利用PFC3D内置FISH语言定义流固之间的压力梯度方程和作用力方程,求解不可压缩流体中两相介质的Navier-Stokes连续方程和运动方程。流体的作用力作为体力施加给砂土体颗粒,同时砂土颗粒的运动和孔隙率等也影响流体运动的参数及其状态,二者充分考虑流固之间的耦合作用。 2.1 砂体流体特性方程 式中:n为孔隙率;t为时间;?表示为梯度,?={?/?x,?/?y,?/?z}T;?p为流体压力梯度;u为流速矢量;ρf为流体密度;τ为流体黏性应力张量;g为重力加速度矢量;fint为单位体积内颗粒与流体的相互作用力,主要包括流固之间拖拽力等,通过流体压力梯度方程求得,其一般形式: 当n0.8时,根据Ergun方程推导可得: 当n≥0.8时,根据Wen and Yu方程推导可得: 式中:?pj为流体压力梯度分量;υf为黏滞系数;为平均粒径;uj为平均流速;为砂土颗粒的平均速度;CD为拖曳系数。 从式(4)、式(5)可见流固间相互作用力与颗粒间孔隙率、流体颗粒之间的相对速度、流体密度、粒径等因素密切相关,可见水土相互作用贯穿渗流的整个过程。 2.2 颗粒的接触和接触、接触力 式中:mp为颗粒的质量;vp为颗粒速度矢量;wp为颗粒转动速度矢量;为颗粒速度矢量对时间的导数;为颗粒转动速度矢量对时间的导数;Ip为颗粒的转动惯量;fg为重力加速度矢量;fc为接触c(c=1,2,…)处颗粒之间的接触力;rc为接触c处指向颗粒中心的方向矢量;fd为流体施加于颗粒p的作用力,包括浮力及固液之间的相互作用力,大小与压力梯度有关。 2.3 求解方法 3 预聚焦滤波层数值模型的构建 3.1 数值模型的建立 通过配置级配不良的管涌型砂样为基料,改变滤料的粒径来满足不同层间系数比的需要,并生成试样。因为制作的数值试样级配较宽,细颗粒数量多。考虑到现有的计算机运行速度的局限性,且控制细颗粒的数量在合理范围,这里运用土工离心机试验中的相似性原理,将数值模型的尺寸定为20mm×20mm×40mm(长×宽×高),模型的侧壁是由4道不透水的刚性墙组成,顶部和底部用网状墙代表透水边界条件,以便能在底部施加水压,并能够在表面自由排水以及流失颗粒。与实际砂样相同,颗粒粒径范围为0.1~5.0mm,具体试验与数值的砂样级配曲线见图2。流体单元的尺寸为5mm×5mm×5mm,总共4×4×8单元。数值模型参数见表1。 颗粒间的接触刚度模型采用线性接触模型(The Linear Contact Model),由于模拟砂土情况,无黏结力,颗粒间