Numerical simulation of negative pressure penetration of suction anchor foundation based on CFD-DEM fluid solid coupling method
-
摘要: 吸力锚基础作为一种新型高效的海洋结构物基础形式在海洋工程中得到了广泛应用。对于砂土中的吸力锚沉贯研究,传统试验或有限元数值模拟存在一定局限性。采用基于离散单元法与计算流体动力学理论(CFD-DEM)流固耦合方法对吸力锚在砂土中的吸力贯入过程进行数值模拟分析;通过与室内物理模型试验、沉贯阻力理论解析计算结果进行对比分析,验证CFD-DEM流固耦合方法的有效性和准确性。进一步地,从土颗粒细观尺度上深入分析吸力锚在砂土中沉贯特性。数值模拟捕捉到贯入过程中锚内土塞和砂层变化现象,发现砂层呈现出中间向上凸起的弧状分布,说明贯入产生挤土效应,其造成的土体位移和膨胀也是土塞产生的原因。同时模拟得到贯入过程中的超孔隙水压力等势线分布变化情况和锚内砂层水力梯度的变化情况,证实了砂土中吸力贯入的“安全机制”。最后将数值所得的锚外负压比变化与Houlsby和Byrne理论模型计算结果进行对比,得到锚内土体渗透系数随沉贯的变化规律。通过本研究验证,采用离散元法结合计算流体动力学理论有效地模拟和捕捉吸力锚沉贯过程中的细观土水相互作用、宏观土体变形及渗流场分布等现象,该方法可作为吸力锚沉贯特性研究的一种有效手段。
-
关键词:
- CFD-DEM耦合方法 /
- 吸力锚 /
- 土塞 /
- 负压沉贯
Abstract: As a new-type and efficient marine structure foundation form, the suction anchor is widely used in ocean engineering. However, the traditional laboratory tests and the finite element numerical simulations have limitations in studying the penetration of suction anchors in granular materials such as sand. The fluid-solid coupling method based on the discrete element method and computational fluid dynamics theory (CFD-DEM) is used to numerically simulate the penetration of the suction anchors in sand. The effectiveness and accuracy of the CFD-DEM fluid-solid coupling method are verified by comparing with these of the laboratory suction penetration tests and the theoretical analysis of penetration resistance. The penetration mechanism of the suction anchors in sand is explained from the mesoscopic scale of soil particles. The changes of soil plug and sand layer inside the anchor during the penetration process are obtained, which shows that the sand layer presents a convex arc distribution in the middle, indicating the compaction effects. The soil displacement and expansion caused by the compaction effects are also part of the reasons for the soil heave. Simultaneously, the equipotential distribution of the excess pore water pressure changes and the changes in the hydraulic gradient of the sand layer during the penetration process are simulated, which confirms the "safety mechanism" of suction penetration in the sand. Finally, the change of negative pressure loss ratio is compared with the results of Houlsby and Byrne's theoretical model, and the variation law of permeability coefficient of soil inside the anchor with penetration is obtained. In this study, it is proved that the CFD-DEM is an efficient method to analyze the suction installation of the suction anchors in sand.-
Keywords:
- CFD-DEM coupling method /
- suction anchor /
- sand plug /
- negative pressure penetration
-
0. 引言
淤堵按其产生的原因可分为机械淤堵、化学淤堵和生物淤堵等3种类型。岩土工程中,机械淤堵多被称为“渗透淤堵”,它指的是细小土颗粒在渗流作用下进入岩土多孔介质中滞留下来进而降低其渗透性的现象[1]。渗透淤堵在某些工况下具有积极的作用,如坝基淤堵能减少水库的渗漏量、基坑止水帷幕淤堵能提高其截水效果等[2-3]。但另一方面,渗透淤堵也常有不利影响甚至造成严重灾害,如减压井等工程排渗系统因淤堵而失效、斜坡体内淤堵引发滑坡、堰塞坝淤堵诱发漫顶溃决等[4-6]。因此,研究岩土多孔介质的渗透淤堵特性对于工程实践中“趋其利、避其害”具有重要意义。
截至目前,国内外学者大多通过物理模拟试验研究岩土体或土工织物的渗透淤堵作用,如径向渗透模型试验、梯度比试验、砂槽模型试验等[3-4,7-8]。物理模拟试验可获得试样淤堵过程中的流量、水力梯度和淤堵后的渗透系数等宏观水力参数,但无法洞悉淤堵后滞留颗粒在多孔介质中的空间分布特征及其与孔-喉结构的关系,从而也难以明晰渗透淤堵的内在机制。
岩土多孔介质中的渗透淤堵过程本质上是渗流-颗粒的相互作用即流-固耦合,数值模拟技术是研究此类问题的重要手段。Dabir等[9]基于孔隙网络模拟了多孔介质中颗粒运移和滞留造成的淤堵。但孔隙网络模拟技术将多孔介质复杂的孔-喉结构简化为规则的球体-圆柱网络,与实际情况差别较大。李识博等[10]基于商业离散元软件PFC中的流体动力学-离散单元法耦合(CFD-DEM)模块模拟了坝基粗粒土中的渗透淤堵作用。CFD-DEM中的流体计算基于宏观连续介质假设并采用经验或半经验的渗流模型(如达西定律)进行求解,且采用的是粗糙网格和局部平均的方案,难以进行更为精细的孔隙尺度的模拟[11-12]。
将基于介观的格子Boltzmann方法(LBM)与离散单元法耦合形成的LBM-DEM能有效克服上述CFD-DEM的不足,更适合从孔隙尺度模拟复杂的流体-颗粒系统[13-15]。近年来,Wang等[16-17]采用LBM-DEM耦合算法成功地模拟了接触侵蚀和水力劈裂等复杂的岩土流-固耦合问题。为此,本文拟采用LBM模拟岩土多孔介质中的渗流,采用DEM模拟土颗粒在多孔介质中的相互作用与运动,并采用浸入运动边界法(IMB)处理移动的固体颗粒与流体之间的相互作用,基于LBM-IMB-DEM对岩土多孔介质的渗透淤堵过程开展孔隙尺度的模拟,探究多孔介质骨架颗粒分布、骨架颗粒粒径、骨架颗粒定向性、孔隙率及水压力等因素对其渗透淤堵特性的影响。
1. 数值算法
1.1 LBM基本原理
本文基于介观的LBM[18-19]模拟岩土多孔介质中的渗流。在LBM中,流体的计算域由边长为h(格子步长)的立方体格子所划分,在每个节点上定义一个粒子分布函数
fα(r,t) 表示在r 位置处t 时刻以eα 速度运动的流体粒子的数目。在离散的每个时间步Δt 内,流体粒子以离散的速度eα 运动到相邻的节点上,紧接着在该节点处与其它的粒子发生碰撞,碰撞后粒子的分布函数将更新,该过程称为迁移与碰撞,可用公式分别表示为fα(r+eαΔt, t+Δt)=fα(r,t), (1) fα(r,t)=Ω(fα(r,t)), (2) 式中,
Ω 为碰撞算子,它规定了流体粒子的碰撞行为。通过迁移与碰撞的不断循环,流体粒子逐步趋向平衡态分布feqα 。本文采用最为常用的碰撞算子即单松弛时间BGK碰撞模型[20]:
Ω(fα(r,t))=fα(r,t)−Δtτ(fα(r,t)−feqα(r,t)), (3) 式中,
feqα(r,t) 为局部平衡态分布函数,τ为无量纲弛豫时间。本文LBM采用的离散速度模型为常用的D2Q9模型[21],其离散速度、平衡态分布函数及相关计算参数为
eα={(0,0) (α=0)c{cos[(α−1)π2],sin[(α−1)π2]} (α=1,2,3,4)√2c{cos[(2α−1)π4],sin[(2α−1)π4]} (α=5,6,7,8), (4) feqα=ρωα[1+eαuc2s+(eαu)22c4s−u22c2s], (5) c=h/Δt, (6) c2s=c2/3, (7) ωα={4/9 (α=0)1/9 (α=1,2,3,4)1/36 (α=5,6,7,8) , (8) 式中,c为格子速度,
cs 为格子声速,ωα为权系数,ρ为流体密度,u为宏观速度。根据质量和动量守恒定律,模型的宏观密度ρ,速度u和压力p分别定义为
ρ=8∑α=0fα, (9) u=1ρ8∑α=0fαeα, (10) p=ρc2s。 (11) 此外,流体的运动黏滞系数υ与无量纲弛豫时间τ及格子参数间关系为
υ=c2s(τ−12)Δt。 (12) 在LBM中,还需要处理好流体的边界,如压力边界和速度边界等。另外,为计算方便一般将物理量做无量纲处理,即采用格子单位。关于边界处理及单位转换可参考文献[13,18,19]。
采用LBM计算低黏度的流体如空气和水时常出现数值不稳定,将工程中常用的大涡模拟(LES)引入LBM中可有效解决该问题,具体可参见文献[13]。
1.2 IMB方法
本文采用IMB方法[22]处理流体流动与移动的固体颗粒之间的相互作用。在IMB中,需要先将移动的固体颗粒用流体节点进行离散并分别标记出固体内部节点、固体边界节点和流体边界节点,具体标记方法可参考文献[13]。针对这些节点,在LBM的碰撞算子
Ω 中要引入1个附加碰撞项ΩSα 来表征固体颗粒对流体的作用,其公式为:Ω=−Δtτ(1−B)[fα(r,t)−feqα(r,t)]+BΩSα, (13) ΩSα=[f−α(r,t)−feq−α(ρ,u)]−[fα(r,t)−feqα(ρ,uP)], (14) 式中,B为附加碰撞项的权函数,
−α 为方向α 的反方向,u ,uP 分别为节点r 处流体和固体的速度。B值由节点的固体占比ε (节点处控制单元被固体覆盖的体积分数)计算得到,B=ε(τ/Δt−0.5)(1−ε)+(τ/Δt−0.5)。 (15) 另一方面,流体对移动的固体颗粒的作用由流体流动施加在颗粒上的力和力矩表征,
Ff=chn∑k=1(Bk8∑α=0ΩSαeα), (16) Tf=chn∑k=1[(rk−rP)×(Bk8∑α=0ΩSαeα)], (17) 式中,n为该固体颗粒所标记节点的数量,
(rk−rP) 为从颗粒中心指向第k个节点的位置向量。1.3 DEM基本理论
本文采用基于软球模型的DEM[23]模拟固体颗粒之间的相互作用,即允许相接触的两个颗粒间有少量重叠。采用线性接触刚度模型和摩擦滑移模型模拟相接触的两颗粒间的本构行为(力-位移法则),
Fc=Fnc+Ftc=Fnn+Ftt, (18) Fn=knδn, (19) Ft=ktδt≤μFn, (20) Tc=R×Ftc, (21) 式中,
Fc 为颗粒间接触力,n ,t 为接触的法向和切向单位向量,Fnc ,Fn 为法向接触力矢量及其大小,Ftc ,Ft 为切向接触力矢量及其大小,kn ,kt 分别为接触法向刚度和切向刚度,δn ,δt 分别为接触法向重叠量和切向相对位移,μ 为颗粒间摩擦系数,Tc 为接触力产生的力矩,R 为从颗粒中心指向接触点的位置向量。固体颗粒可能受到的力包括颗粒间的接触力、颗粒与固体边界墙的接触力、流体的作用力以及重力等,在这些力的作用下颗粒的运动遵循牛顿定律(运动法则),
ma=Ff+∑iFc i+Fb, (22) I¨θ=∑iTci+Tf, (23) 式中,
m ,I 分别为颗粒的质量和惯性矩,a ,¨θ 分别为颗粒的线加速度和角加速度,Fc i ,Tci 分别为第i个接触产生的接触力和力矩,Fb 表示其它作用力如重力等。颗粒的运动方程采用显式中心有限差分法求解,以获得每个时刻颗粒的速度和位置。1.4 LBM-IMB-DEM数值计算平台及其验证
基于1.1—1.3节所述算法的基本理论,运用Fortran语言开发了LBM-IMB-DEM耦合数值计算平台,其中LBM模拟流体流动、DEM模拟固体颗粒间的相互作用与运动、IMB模拟移动的固体颗粒与流体之间的相互作用,耦合程序的计算流程如图1所示。
需要说明的是,LBM与DEM均采用显式时步循环算法且有各自的时间步长,一般来说DEM的时间步长
ΔtD 比LBM的时间步长Δt 要小,因此可将DEM的时间步长适当缩短使得Δt/ΔtD 为整数,进而可将DEM循环作为LBM的子循环(如图1所示)以达到耦合计算时间的统一。本文通过一系列经典算例对LBM-IMB-DEM数值计算平台进行了多方面的验证,主要包括:模拟二维Poiseuille流和圆柱绕流验证了LBM的准确性;模拟砂堆和双轴试验并与PFC模拟的结果对比验证了DEM的准确性;模拟单颗粒沉降和两个颗粒沉降DKT的过程验证了IMB的准确性。限于篇幅,这里只给出单颗粒下沉模拟的参数和结果。
单颗粒下沉算例常被用来验证颗粒-流体耦合算法的准确性。本文建立的单颗粒下沉模型中流体计算域为宽2 cm、高6 cm的长方形区域,将1个直径为0.25 cm的圆形固体颗粒放置于坐标为(1 cm,4 cm)的位置(原点位于计算域的左下角点)。固体颗粒的密度和法向刚度分别为1.25 g/cm3和1.0×107 N/m,计算域四周为固体边界墙且其刚度为1.0×108 N/m,不考虑颗粒与边界墙接触时的切向作用。流体的密度与运动黏度分别为1.0 g/cm3和1.0×10-5 m2/s,流体计算域的格子步长取为0.01 cm,弛豫时间为0.65,区域四周为不透水边界并由标准反弹准则[18]进行处理。LBM-IMB-DEM计算过程中,固体颗粒在重力作用下开始向下运动并受到向上的流体力作用,随后颗粒一直下落直到底边界并逐渐稳定下来,将该过程中颗粒竖向位置随时间的变化与Wan等[24]基于CFD模拟的结果进行对比,如图2所示。可以看到,两者吻合较好,这说明LBM-IMB-DEM耦合算法能较准确地模拟颗粒-流体耦合问题。
2. 渗透淤堵的LBM-IMB-DEM模拟
实际的岩土体中颗粒粒径、颗粒形态、颗粒定向性与孔隙等内部结构异常复杂,为了便于系统地研究其淤堵特性,本文建立了简化的岩土多孔介质渗透淤堵二维模型,如图3所示。模型总长50 mm、宽20 mm,其中多孔介质区域长为40 mm(图3中2条虚线之间),它由一定尺寸的骨架颗粒按一定孔隙率随机分布所构成,这些骨架颗粒在后续的数值计算过程中保持固定。图3中左侧绿色线框内的灰色圆盘为添加的拟淤堵颗粒,其直径为0.8 mm。模型的上、下边界不透水且为固体颗粒的墙边界,水流从左侧进入、右侧流出,左、右边界为压力边界并由非平衡态反弹法[18]处理,上、下不透水边界及固定的骨架颗粒由标准反弹法[18]处理。基于LBM-IMB-DEM数值计算平台对模型进行计算以模拟水流作用下拟淤堵颗粒在多孔介质中的迁移和沉积过程,主要计算参数列入表1。
表 1 渗透淤堵LBM-IMB-DEM模拟参数Table 1. LBM-IMB-DEM simulation parameters for seepage clogging固体参数 参数值 流体参数 参数值 颗粒密度/(kg·m-3) 2650 流体密度/(kg·m-3) 1000 法向刚度/(N·m-1) 1×107 运动黏度/(m2·s-1) 1×10-6 切向刚度/(N·m-1) 5×106 弛豫时间 0.51 摩擦系数 0.3 格子步长/m 8×10-5 DEM子循环步数 4 LBM时间步长/s 2.13×10-5 本文定义“淤堵渗透系数比”来表征多孔介质渗透淤堵的严重程度,其定义为多孔介质淤堵后渗透系数与淤堵前渗透系数的比值,该值越小表示淤堵越严重。另外,还定义了“滞留体积百分比”来表征淤堵后多孔介质中滞留颗粒数量的多少,其定义为多孔介质中滞留颗粒的总体积与孔隙体积之比(百分数)。
待颗粒运动与流场都稳定后停止淤堵过程的模拟。随后,针对淤堵前、后的多孔介质模型(图3中2条红色线之间区域),在保持其中颗粒(包括淤堵颗粒)位置不变的条件下,将其入口与出口密度分别设为1000.1和1000.0(雷诺数较小,基本为层流)并进行迭代计算,当流场收敛到稳态后停止迭代并由达西定律计算其渗透系数,据此可进一步得到该模型的淤堵渗透系数比。需要指出的是下文中所涉及的流场图也均是在此条件下获得的。
3. 模拟结果分析
本文主要模拟了不同骨架颗粒分布、骨架颗粒粒径、多孔介质孔隙率、骨架颗粒定向性及水流入口压力(水流入口密度)等条件下多孔介质的渗透淤堵,下面就这些条件下的模拟结果进行分析。需要说明的是,格子单位与物理单位的比值为一个固定的参考量,也就是说格子单位的变化规律与物理单位的变化规律是完全一致的,一般为了计算方便在计算过程中统一取格子单位。为方便起见,在下文的分析中也直接采用了格子单位(l.u.)下得到的模拟结果。
3.1 骨架颗粒分布的影响
在骨架颗粒是粒径(直径)为37.5(3 mm)的圆形,孔隙率为0.7,水流入口与出口密度分别为1001.0和1000.0的条件下,通过改变随机数种子生成了5种不同骨架颗粒随机分布的多孔介质模型并进行渗透淤堵过程的LBM-IMB-DEM数值模拟。模拟结果如下:5个模型的淤堵渗透系数比分别为0.41,0.36,0.43,0.37和0.43,平均值为0.40;5个模型的滞留体积百分比分别为3.82%,2.53%,2.51%,2.82%和2.55%,平均值为2.85%。其中第1,2种骨架颗粒分布的模型淤堵前后的流场与滞留颗粒分布情况如图4所示。
由模拟结果及图4可知:渗透淤堵发生后,尽管多孔介质中滞留颗粒体积只占总孔隙体积很小一部分(小于5.0%),但这些颗粒多驻留在孔喉处,致使流动通道数显著减少,因而淤堵后的渗透系数仅为淤堵前的40%;骨架颗粒分布1的模型中滞留颗粒相对于分布2时的较多,但部分滞留颗粒位于离入口较近的骨架颗粒表面和非流动区域,而分布2的滞留颗粒基本都位于孔喉处,所以分布2的淤堵渗透系数比反而较小,即淤堵更为严重。
由上述分析可得,多孔介质中骨架颗粒的空间分布不同,其孔隙结构也不同,渗透淤堵特性也有所区别。为此,在下文分析不同骨架颗粒粒径等条件下(除不同水流入口压力的情形)多孔介质渗透淤堵时,淤堵前渗透系数、淤堵渗透系数比和滞留体积百分比均采用5个随机分布模型的平均值。
3.2 骨架颗粒粒径的影响
在孔隙率为0.7,骨架颗粒是圆形,水流入口与出口密度分别为1001.0和1000.0的条件下,分别生成骨架颗粒粒径为25,37.5和50的多孔介质模型并进行渗透淤堵过程的LBM-IMB-DEM数值模拟。各模型淤堵后的淤堵渗透系数比、滞留体积百分比及淤堵前的渗透系数如图5所示,淤堵前、后的流场和滞留颗粒分布情况如图6所示。
由图5,6可知,随着骨架颗粒粒径的增大,淤堵前多孔介质中的流动通道数虽然减少,但其宽度显著增加,透水性也随之增大,细颗粒在水流携带下也更容易通过,故淤堵后滞留体积百分比逐渐减小,淤堵渗透系数比不断增大。因此,骨架颗粒粒径越小,多孔介质越容易发生淤堵,淤堵的程度也更严重,如图6(a)中细颗粒很快就堵塞了入口附近的孔喉,后续颗粒继续堆积在入口处,最终所有流动通道都被完全堵塞。
3.3 多孔介质孔隙率的影响
在骨架颗粒是37.5的圆形、水流入口与出口密度分别为1001.0和1000.0的条件下,分别生成孔隙率为0.6,0.7和0.8的多孔介质模型并进行渗透淤堵过程的LBM-IMB-DEM数值模拟。各孔隙率的模型淤堵前渗透系数及淤堵后的淤堵渗透系数比和滞留体积百分比如图7所示,淤堵前、后的流场和滞留颗粒分布情况如图8所示。
由图7,8可见,随着孔隙率的增加,多孔介质中孔隙变大,连通性也更好,故其透水性增强,细颗粒在水流携带下也更容易通过,滞留体积百分比随之减小,淤堵渗透系数比逐渐增大。因此,孔隙率越小,多孔介质越容易发生淤堵。
3.4 骨架颗粒定向性的影响
为了探究多孔介质中骨架颗粒定向性对其渗透淤堵特性的影响,本文建立了如图9所示的条形颗粒,其由3个同面积的圆组合而成且颗粒面积等于直径为37.5的圆形颗粒的面积。在骨架颗粒为条形、多孔介质孔隙率为0.7、水流入口与出口密度分别为1001.0和1000.0的条件下,分别生成颗粒长轴方向平行于流动方向、垂直于流动方向和随机的3种多孔介质模型。
通过渗透淤堵过程的LBM-IMB-DEM数值模拟得到这3种模型淤堵后的淤堵渗透系数比和滞留体积百分比以及淤堵前渗透系数,如图10所示,淤堵前、后的流场和滞留颗粒分布情况如图11所示。
由图10,11可见:骨架颗粒长轴平行于流动方向时,多孔介质中流动通道较宽、连通性较好且迂曲度较小,故其透水性最强,细颗粒在水流携带下也最容易通过,淤堵渗透系数比最大;骨架颗粒长轴垂直于流动方向时,多孔介质中流动方向的孔喉直径较小且渗流通道的迂曲度较大,故其淤堵前透水性最弱,最后被细颗粒完全堵塞;骨架颗粒长轴方向随机时多孔介质的淤堵渗透系数比居于前两者之间。
从图5,7,10可以看出,表征多孔介质渗透淤堵程度的淤堵渗透系数比的变化规律与淤堵前渗透系数的变化是一致的,这也表明透水性强的多孔介质更不容易发生淤堵。然而,淤堵渗透系数比和滞留体积百分比的变化规律并不完全一致,如前所述的图4中两种情形,图10中也出现了局部的不一致情况,这表明多孔介质中滞留颗粒越多并不一定代表其淤堵越严重。
3.5 水压力的影响
为了探究水压力对多孔介质渗透淤堵特性的影响,将图4(b)中所示多孔介质模型分别在水流入口密度为1000.5,1000.75和1001.0(出口密度均为1000.0)的条件下进行渗透淤堵过程的LBM-IMB-DEM模拟。这3种条件下多孔介质淤堵后的淤堵渗透系数比和滞留体积百分比如图12所示。
由图12可见,当水流出口密度保持不变,随着水流入口密度的增加,多孔介质淤堵后的滞留体积百分比逐渐增加,淤堵渗透系数比逐渐减小。究其原因,水流入口的密度越大即水压力越大(式(11)),流速也就越大,在较大的流速下入口处的细颗粒运动速度较快并趋于在多孔介质中分散开来,这样其驻留在孔喉处的概率更大,最终滞留的颗粒更多,淤堵也更为严重。如图13所示,入口密度为1001.0时细颗粒堵塞了多孔介质上下边界附近的流动通道,但入口密度为1000.5时这些部位并未发生堵塞。
4. 结论
本文基于LBM-IMB-DEM耦合算法从孔隙尺度模拟了岩土多孔介质渗透淤堵的过程,分析了多孔介质骨架颗粒分布、骨架颗粒粒径、孔隙率、骨架颗粒定向性以及水压力等因素对多孔介质淤堵后滞留体积百分比和淤堵渗透系数比的影响,并通过滞留颗粒在多孔介质中的分布及相应流场的变化情况探究了各因素影响的内在机制。
(1)LBM-IMB-DEM耦合算法能较准确地模拟颗粒-流体动力耦合问题,在岩土流-固耦合精细模拟等研究领域具有广阔的应用前景。
(2)岩土多孔介质中骨架颗粒的空间分布不同,其孔隙结构特征也不同,渗透淤堵特性也有所区别。
(3)岩土多孔介质的骨架颗粒粒径或孔隙率越小,渗透淤堵后的滞留颗粒越多,淤堵也更为严重。
(4)岩土多孔介质中骨架颗粒的定向性对其渗透淤堵特性有一定的影响,当骨架颗粒为条形且其长轴与流动方向垂直时淤堵最为严重。
(5)随水压力增大,渗透淤堵后滞留颗粒在多孔介质中更趋于分散,淤堵更为严重。
(6)岩土多孔介质渗透淤堵的程度与其淤堵前的渗透性密切相关,表现为透水性强的多孔介质相对不易发生淤堵。
后续研究还应进一步考虑其他因素如拟淤堵颗粒的粒径与形状、骨架颗粒的粒径分布,骨架颗粒的形状、流体的黏度等对多孔介质渗透淤堵特性的影响,并应进一步开发三维的LBM-DEM耦合计算程序以便开展渗透淤堵过程的三维模拟。
-
表 1 模型试验中吸力锚模型尺寸参数
Table 1 Dimensions of model suction anchor
锚长度H/
mm内径Di/
mm外径Do/
mm壁厚t/
mm长径比 150 150 158 4 1 -
[1] 鲁晓兵, 郑哲敏, 张金来. 海洋平台吸力式基础的研究与进展[J]. 力学进展, 2003, 33(1): 27-40. doi: 10.3321/j.issn:1000-0992.2003.01.005 LU Xiaobing, ZHENG Zhemin, ZHANG Jinlai. Progress in the study on the bucket foundations of offshore platform[J]. Advances in Mechanics, 2003, 33(1): 27-40. (in Chinese) doi: 10.3321/j.issn:1000-0992.2003.01.005
[2] BYRNE B, HOULSBY G, MARTIN C, et al. Suction caisson foundations for offshore wind turbines[J]. Wind Engineering, 2002, 26(3): 145-155. doi: 10.1260/030952402762056063
[3] SHONBERG A, HARTE M, AGHAKOUCHAK A, et al. Suction bucket jackets for offshore wind turbines: applications from in situ observations[C]// Proceedings of TC 209 Workshop – 19th International Conference on Soil Mechanics and Geotechnical Engineering. Seoul, 2017.
[4] STURM H. Design aspects of suction caissons for offshore wind turbine foundations[C]// Proceedings of TC 209 Workshop – International Conference on Soil Mechanics and Geotechnical Engineering. Seoul, 2017.
[5] 陈飞, 练继建, 马煜祥, 等. 粉砂中筒型基础沉贯过程筒-土作用机理试验研究[J]. 岩土工程学报, 2015, 37(4): 683-691. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201504017.htm CHEN Fei, LIAN Jijian, MA Yuxiang, et al. Model tests on bucket-soil interaction during installation of bucket foundation in silt sand[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(4): 683-691. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201504017.htm
[6] 王胤, 朱兴运, 杨庆. 考虑砂土渗透性变化的吸力锚沉贯及土塞特性研究[J]. 岩土工程学报, 2019, 41(1): 184-190. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201901025.htm WANG Yin, ZHU Xingyun, YANG Qing. Installation of suction caissons and formation of soil plug considering variation of permeability of sand[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(1): 184-190. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201901025.htm
[7] 李大勇, 张雨坤, 高玉峰, 等. 中粗砂中吸力锚的负压沉贯模型试验研究[J]. 岩土工程学报, 2012, 34(12): 2277-2283. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201212022.htm LI Dayong, ZHANG Yukun, GAO Yufeng, et al. Model tests on penetration of suction anchors in medium-coarse sand[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(12): 2277-2283. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201212022.htm
[8] TRAN M N, RANDOLPH M F, AIREY D W. Study of seepage flow and sand plug loosening in installation of suction caissons in sand[C]// Proceedings of The Fifteenth International Offshore and Polar Engineering Conference. Korea, 2005.
[9] RAGNI R, BIENEN B, STANIER S, et al. Observations during suction bucket installation in sand[J]. International Journal of Physical Modelling in Geotechnics, 2020, 20(3): 132-149. doi: 10.1680/jphmg.18.00071
[10] ZEINODDINI M, MOUSAVI S A, ABDI M R. Simulation of suction caisson penetration in seabed using an adaptive mesh technique[J]. Procedia Engineering, 2011, 14: 1721-1728. doi: 10.1016/j.proeng.2011.07.216
[11] 吕阳, 王胤, 杨庆. 吸力式筒形基础沉贯过程的大变形有限元模拟[J]. 岩土力学, 2015, 36(12): 3615-3624. doi: 10.16285/j.rsm.2015.12.035 LÜ Yang, WANG Yin, YANG Qing. Large deformation finite element simulation of suction bucket foundation during penetration[J]. Rock and Soil Mechanics, 2015, 36(12): 3615-3624. (in Chinese) doi: 10.16285/j.rsm.2015.12.035
[12] 丁红岩, 韩雪松, 张浦阳, 等. 离散元法模拟吸力锚施工中生成的土塞研究[J]. 海洋技术, 2004, 23(1): 52-56. doi: 10.3969/j.issn.1003-2029.2004.01.012 DING Hongyan, HAN Xuesong, ZHANG Puyang, et al. The simulation of soil plug formed in constructing of suction anchor by discrete element method[J]. Ocean Technology, 2004, 23(1): 52-56. (in Chinese) doi: 10.3969/j.issn.1003-2029.2004.01.012
[13] CHEN B F, HUANG T T. On fluid-filled mixture model for suction pile foundation analysis[J]. Ocean Engineering, 2019, 188: 106306. doi: 10.1016/j.oceaneng.2019.106306
[14] 王胤, 艾军, 杨庆. 考虑粒间滚动阻力的CFD- DEM流-固耦合数值模拟方法[J]. 岩土力学, 2017, 38(6): 1771-1780. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201706027.htm WANG Yin, AI Jun, YANG Qing. A CFD-DEM coupled method incorporating soil inter-particle rolling resistance[J]. Rock and Soil Mechanics, 2017, 38(6): 1771-1780. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201706027.htm
[15] ANDERSON T B, JACKSON R. Fluid mechanical description of fluidized beds. equations of motion[J]. Industrial & Engineering Chemistry Fundamentals, 1967, 6(4): 527-539.
[16] AI J, CHEN J F, ROTTER J M, et al. Assessment of rolling resistance models in discrete element simulations[J]. Powder Technology, 2011, 206(3): 269-282. doi: 10.1016/j.powtec.2010.09.030
[17] CHENG K, WANG Y, YANG Q, et al. Determination of microscopic parameters of quartz sand through tri-axial test using the discrete element method[J]. Computers and Geotechnics, 2017, 92: 22-40. doi: 10.1016/j.compgeo.2017.07.017
[18] DI FELICE R. The voidage function for fluid-particle interaction systems[J]. International Journal of Multiphase Flow, 1994, 20(1): 153-159. doi: 10.1016/0301-9322(94)90011-6
[19] WANG Y, ZHOU L X, YANG Q. Hydro-mechanical analysis of calcareous sand with a new shape-dependent fluid-particle drag model integrated into CFD-DEM coupling program[J]. Powder Technology, 2019, 344: 108-120.
[20] CLIFT R, GAUVIN W H. Motion of entrained particles in gas streams[J]. The Canadian Journal of Chemical Engineering, 1971, 49(4): 439-448.
[21] WADELL H. The coefficient of resistance as a function of Reynolds number for solids of various shapes[J]. Journal of the Franklin Institute, 1934, 217(4): 459-490.
[22] 吴野. 天然岩土颗粒材料液体中拖曳力系数试验研究[D]. 大连: 大连理工大学, 2017. WU Ye. Experimental Study on the Drag Force Coefficient of Natural Geotechnical Particle Materials in the Liquid[D]. Dalian: Dalian University of Technology, 2017. (in Chinese)
[23] 王胤, 周令新, 杨庆. 基于不规则钙质砂颗粒发展的拖曳力系数模型及其在细观流固耦合数值模拟中应用[J]. 岩土力学, 2019, 40(5): 2009-2015. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201905046.htm WANG Yin, ZHOU Lingxin, YANG Qing. New drag coefficient model for irregular calcareous sand particles and its application into fluid-particle coupling simulation[J]. Rock and Soil Mechanics, 2019, 40(5): 2009-2015. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201905046.htm
[24] KLOSS C, GONIVA C L. A new open source discrete element simulation software[C]//Proceedings of 5th international conference on discrete element methods.
[25] CHRLSTOPHER J G. The Open Source CFD Toolbox[R]. Paris: OpenCFD Ltd, 2009.
[26] 朱兴运. 吸力锚基础负压沉贯阻力及土塞形成机理研究[D]. 大连: 大连理工大学, 2018. ZHU Xingyun. Research on Penetration Resistance and Formation Mechanism of Soil Plug during Suction Anchor Installation[D]. Dalian: Dalian University of Technology, 2018. (in Chinese)
[27] TJELTA T. The suction foundation technology[C]// Proceedings of the Third International Symposium Frontiers in Offshore Geotechnics, Oslo, 2015.
[28] WANG Y, CHENG K, YANG Y F, et al. Microscopic mechanical analysis of sand production using a new arbitrary resolved-unresolved CFD-DEM model[J]. International Journal of Multiphase Flow, 2022, 149: 103979.
[29] HOULSBY G T, BYRNE B W. Design procedures for installation of suction caissons in sand[J]. Proceedings of the Institution of Civil Engineers - Geotechnical Engineering, 2005, 158(3): 135-144.
[30] ANDERSEN K H, JOSTAD H P, DYVIK R. Penetration resistance of offshore skirted foundations and anchors in dense sand[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2008, 134(1): 106-116.
[31] GONZÁLEZ-MONTELLANO C, FUENTES J M, AYUGA-TÉLLEZ E, et al. Determination of the mechanical properties of maize grains and olives required for use in DEM simulations[J]. Journal of Food Engineering, 2012, 111(4): 553-562.
[32] 姚涛, 余建星, 李志刚, 等. 筒型基础贯入过程中的渗流分析[J]. 海岸工程, 2012, 31(1): 8-15. https://www.cnki.com.cn/Article/CJFDTOTAL-HAGC201201003.htm YAO Tao, YU Jianxing, LI Zhigang, et al. Analysis of seepage field in the penetration of bucket foundation[J]. Coastal Engineering, 2012, 31(1): 8-15. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HAGC201201003.htm
[33] ERBRICH C T, TJELTA T I. Installation of bucket foundations and suction caissons in sand-geotechnical performance[C]// Proceedings of the 31th Annual Offshore Technology Conference. Texas, 1999.
[34] HARIRECHE O, MEHRAVAR M, ALANI A M. Suction caisson installation in sand with isotropic permeability varying with depth[J]. Applied Ocean Research, 2013, 43: 256-263.
-
其他相关附件