CFD-DEM simulation of microbially treated sands under undrained consolidated cyclic triaxial tests
-
摘要: 微生物诱导碳酸钙沉积是一种新型的地基处理技术,处理后的土体可以看成一种结构性土。首先,在已有三维含颗粒抗转动和抗扭转模型及三维胶结破坏准则的基础上,通过考虑颗粒碰撞接触过程中颗粒本身的塑性变形及率相关性的接触黏滞阻尼,建立考虑循环荷载作用下的三维胶结模型；然后,参考已有研究,建立了反硝化反应在加固砂土中的时效性关系。并引入CFD-DEM耦合程序,用以模拟分析不同胶结含量以及不同气泡含量下,微生物处理砂土在固结不排水循环剪切试验中的力学特性；最后,从宏微观角度分析生物胶结与生物气泡对砂土抗液化性能的影响及其作用机理。研究表明,胶结和气泡共同存在对抗液化能力的提升并没有起到“1+1=2”的效果；胶结的存在提高了非饱和砂土的抗液化能力,明显抑制孔压比和轴向应变的发展,力学配位数得到了提升；而气泡的存在却降低了胶结砂土的抗液化能力,使得胶结砂土达到初始液化的振次减少,轴向应变向受拉方向大幅增长,力学配位数下降明显。
-
关键词:
- 微生物诱导碳酸钙沉积 /
- CFD-DEM耦合 /
- 不排水循环三轴剪切 /
- 地基液化 /
- 离散单元法
Abstract: The microbially induced calcite precipitation is a promising technology to improve ground, and the treated soil can be regarded as the structural one. In this study, firstly, based on the three-dimensional (3D) contact model for granulates incorporating rolling and twisting resistances and 3D bonds failure criteria, and considering both the slight plastic deformation of particles during collisions and the rate-dependency, a cyclic bonded contact model is established. A time-dependent relationship is then proposed to describe the denitrification reaction in reinforced sand. Next, the mechanical responses of microbially treated sands at different cementation and bubble contents are investigated by the coupled CFD-DEM in undrained-consolidated cyclic triaxial tests. The effects of biological bond and biological bubbles on the liquefaction resistance of sands are analyzed in link with the mechanism from macroscopic and microscopic scales. The results show that the coexistence of cementation and bubble does not increase the liquefaction resistance as expected in the form of "1+1=2". The presence of cementation enhances the liquefaction resistance of unsaturated sands evidenced by the decrease of excess pore water pressure ratio and axial strain, and the increase of coordination number. However, the presence of bubbles reduces the liquefaction resistance of cemented sands, where the number of cycles to the initial liquefaction decreases, the axial strain increases significantly in the tensile direction, and the coordination number decreases significantly. -
0. 引言
黄土是指在第四纪地质时代中以风力作用形成的黄色粉土沉积物,是一种在特定环境中形成的具有特殊性质的土。由于黄土地区处于季节冻土区,黄土受季节性冻融作用的影响显著,每年春季发生的冻融灾害非常频繁。此外,对于盐蚀作用诱发的黄土边坡剥落等病害也不容忽视。笔者对陕北黄土高原区边坡盐蚀剥落病害进行了实地调研,每处调研点由表及里（0～30 cm）每隔5 cm取一份土样,用于室内易溶盐分析试验。试验结果表明,边坡表层主要富集Na2SO4等易溶盐,在冻融循环等条件下极易发生反复溶解和结晶的盐蚀作用,使被侵蚀黄土结构损伤扩展,劣化破坏,从而诱发边坡盐蚀剥落病害[1]。因此,开展冻融环境下含盐黄土的盐蚀劣化研究具有重要意义。
目前关于冻融循环作用对盐渍土物理力学性质的影响研究已有大量文献资料。高江平[2]系统开展了盐渍土冻胀、盐胀机理的试验与数值仿真研究工作。李振等[3]以某干渠渠基盐渍土为研究对象,研究了冻胀量与各影响因素之间的变化规律。牛玺荣等[4]推求了考虑冻胀和盐胀耦合作用下硫酸盐渍土体积变化模型并采用自行研制的试验系统对模型进行了验证。黄雪峰等[5]基于室内单次降温试验,对宁夏扶贫扬黄灌溉工程沿线硫酸盐渍土盐胀特性及影响因素进行了试验研究。张莎莎等[6]开展了多因素交互作用下砾类亚硫酸盐渍土的盐胀试验并建立了其盐胀率简化预测模型。此外,部分学者对冻融作用下黄土强度变形特性进行了研究。胡再强等[7]探讨了冻融循环条件下杨凌重塑黄土强度变化规律并构建了其强度劣化系数拟合方程。谷琪等[8]探究了冻融循环过程中湿陷性黄土的变形及湿陷情况。然而关于冻融作用对含盐原状黄土影响的研究较少,含盐量、冻融循环次数及其耦合效应导致土体强度衰减程度的定量化关系尚不明确。
研究学者对黄土及盐渍土微细观结构演化规律亦开展了大量研究工作。陈正汉等[9]结合CT技术和三轴剪切设备,系统研究了黄土在多种应力路径及浸水湿陷过程中的细观结构演化特性并进一步基于CT数据构建了其细观结构演化方程。张伟等[10]基于CT–三轴剪切试验,探究了剪切过程中硫酸盐渍土的应力应变关系和CT图像及CT数的变化规律。叶万军等[11]研究了冻融环境下黄土体结构在不同观测尺度下的损伤演变规律。但关于冻融作用下盐渍原状黄土微结构损伤扩展演化模式的研究尚未见有专门研究报道。
盐渍原状黄土是一种具有典型天然结构强度的特殊土,冻融循环作用对其结构强度影响较大。基于此,本文选取西安Q3原状黄土,人工制备不同Na2SO4含量的盐渍原状黄土试样,通过冻融循环条件下的三轴剪切及CT扫描试验,研究Na2SO4盐渍原状黄土冻融过程盐蚀劣化规律及微细观结构损伤演化机制。研究成果对探究黄土盐蚀作用诱发边坡剥落等病害的成灾机理具有重要的参考意义。
1. 试验材料与试样制备
1.1 试验黄土
试验黄土为取自陕西西安某基坑的Q3原状黄土,基坑深度为10 m,取样深度为8～10 m。原状土样用削土刀在基坑侧壁上刻取。按照GB/T50123—2019《土工试验方法标准》对原状土样进行基本物理性质分析,结果如表1所示。试样粒组含量：>0.05 mm(5%),0.01～0.05 mm(52%),0.005～0.01 mm(24%),<0.005 mm(19%)。
表 1 试样物理性质指标Table 1. Physical properties of undisturbed loess土粒相对密度Gs 干密度ρd/(g·cm-3) 孔隙比e 液限wL/% 塑限w P/% 塑性指数Ip 2.70 1.42 0.92 34.3 19.6 14.7 1.2 试样制备
三轴剪切和CT扫描试样尺寸相同,均削制成直径39.1 mm、高度80 mm的标准圆柱试样。然后通过向原状土样中浸入不同浓度的Na2SO4盐水来人工制备含水率w为20%,Na2SO4含量η分别为0.0%,0.5%,1.0%,1.5%的盐渍原状黄土试样。关于人工制备盐渍原状黄土试样的方法,其具体操作流程参见笔者已发表相关文献[12],在此不再赘述。为进一步验证上述方法的制样效果,特选取3个制备好的标准圆柱试样,其编号分别为T-1（w=20%,η=0.5%）、T-2（w=20%,η=1.0%）、T-3（w=20%,η=1.5%）。将每一个标准试样均等分成4块20 mm高的小圆柱体,以代表试样的4个不同位置,其编号依次为1,2,3,4。然后分别测定每个小圆柱体试块的含水率和含盐量,试验结果如图1所示。由图1可见,试样不同部位含水率和含盐量变化曲线平缓,趋于一条直线,即试样内部含水率和含盐量的分布较为均匀。此外,试样不同部位的含水率和含盐量与目标值基本一致,因而人工制备盐渍原状黄土试样的方法可以满足试验要求。
2. 试验方案
2.1 冻融试验
利用保鲜膜将制备好的含盐原状试样包裹,构成一个不补（散）水的密封环境,随后放入高低温试验箱进行冻融循环试验。由于三轴试样尺寸较小,土样端部和侧面换热条件虽有所差异,但影响不大。因此本次冻融循环试验为封闭系统下的多向快速冻融循环试验,以保证冻融时试样水分迁移较少。
图2所示为陕北榆林和延安地区近十年月平均最低气温变化曲线。从图中可以看出,环境温度变化区间在-20℃～20℃之间,因此以当地最低气温（约-20℃）为冻结温度,最高气温（约20℃）为融化温度,进行冻融循环试验。试验方案如下：为保证试样充分冻融,低温-20℃条件下冻结12 h,高温20℃条件下融化12 h,即一个冻融循环为24 h；冻融循环次数N分别为0,1,2,5,10。
2.2 三轴剪切试验
试样经受不同冻融循环次数后,为避免排水固结对冻融损伤后试样结构状态扰动,进行不固结不排水（UU）三轴剪切试验。试验过程中剪切速率设定为0.4 mm/min；考虑到盐渍原状黄土试样浅层分布特点,其侧向压力相对较小,因而三轴剪切试验采用较低围压,分别为50,100,150,200 kPa。剪切过程中若应力–应变关系表现为有峰值的软化型曲线,当峰值后轴向应变达到3%～5%时,结束试验；反之,若应力应变关系表现为应变硬化型曲线,则以轴向应变达到15%作为剪切终止条件。
2.3 CT扫描试验
CT试验采用Brilliance 16螺旋CT机,其空间分辨率为0.208 mm,CT值范围为-1024～+3071。试验时将经受不同冻融循环次数的含盐原状黄土试样每次扫描3个断面,每个断面之间均间隔20 mm,3个断面从土样底部到顶部依次标记为a,b,c截面。通过CT扫描得到试样的细观结构图像,然后根据CT成像原理,基于物质吸收系数运算可以得到材料内部任意一点的CT数,进一步统计分析可以获得图像任意一个区域CT数均值ME。ME值体现扫描断面内所有物质点的平均密度,密度越小则试样ME值越小,因而其数值变化可以反映试样内部损伤扩展过程[13]。
3. 试验结果与分析
3.1 三轴试验结果分析
（1）应力–应变关系曲线
图3所示为试样应力–应变曲线变化规律。从图中可以看出,冻融循环过程对黄土试样应力–应变曲线的类型及特征无明显影响,冻融循环前后均表现为应变硬化型。此外,相同冻融循环次数条件下,含盐量变化对应力–应变曲线的变化特征亦无显著影响,亦均表现为“硬化型”的变化趋势。
为深入揭示冻融条件下Na2SO4盐渍原状黄土强度劣化特性,取轴向应变
ε1 为15%时的偏应力作为破坏偏应力(σ1−σ3)f ,亦即试样的强度值,如图4所示。从图中可以看出,冻融循环对强度产生明显的劣化效应,破坏偏应力随着冻融循环次数增加逐渐减小,但降低速率逐渐减小,表现出显著的减速劣化特征；冻融条件下Na2SO4含量亦对强度产生显著的盐蚀劣化效应,破坏偏应力随着含盐量增加逐渐减小,且呈现显著的线性或加速劣化特征。（2）黏聚力
图5（a）所示为黏聚力随冻融循环次数变化规律曲线。由图可见,随着冻融循环次数增加,黏聚力逐渐减小,但衰减主要集中在前5次冻融循环,之后逐渐趋于稳定,表现出减速劣化特征,这与前述三轴破坏偏应力随冻融循环次数的劣化规律是一致的（图4（a））。分析其原因,冻融循环作用下孔隙水相态变化对土颗粒产生挤压力,破坏了土颗粒之间的原生结构强度,导致黏聚力降低。多次冻融循环条件下土颗粒排列趋于平衡状态,颗粒联结强度达到稳定残余强度,黏聚力亦趋于稳定。
黏聚力随含盐量变化规律曲线如图5（b）所示。从图中可以看出,除了未经受冻融循环（N=0）试样的黏聚力随含盐量增加无明显变化外,其它试样的黏聚力均随含盐量的增加逐渐减小,且衰减速率逐渐增大,表现出加速劣化特征,这亦与前述三轴破坏偏应力随含盐量的劣化特性表现出了很好的一致性（图4（b））。分析其原因,冻融循环条件下试样内部Na2SO4可溶盐的结晶—溶解—重结晶过程的反复盐蚀作用使得黄土体结构受到破坏并变得较为松散,强度显著降低。未经受冻融循环黄土试样内部的硫酸钠无相态变化亦即无法表现出盐蚀作用,因而其黏聚力无明显变化。
（3）内摩擦角
图6所示为内摩擦角随冻融循环次数及Na2SO4含量变化规律。由图6可见,部分试样的内摩擦角随冻融循环次数增加逐渐减小,部分试样的内摩擦角表现出波动变化特征,但整体变化幅值均较小；随含盐量增加大体上表现出先减小而后增大的特征,但变化幅值亦很小,可以认为内摩擦角随含盐量增大无显著变化。总体而言,笔者分析认为冻融作用下Na2SO4盐渍原状黄土试样内摩擦角变化幅值较小且无显著变化规律。这是由于内摩擦角大小主要取决于土颗粒间的摩阻力和咬合作用,而冻融循环和盐蚀作用对上述因素并无显著影响。
（4）冻融与盐蚀劣化作用解耦分析
冻融循环条件下Na2SO4盐渍原状黄土不仅受到孔隙水相变引起的冻融劣化作用,且要受到Na2SO4易溶盐相态改变引起的盐蚀劣化作用。为定量分析冻融与盐蚀劣化规律及相互关系,以下基于冻融循环条件下黏聚力变化规律,对冻融与盐蚀劣化作用进行解耦分析。
图7所示为冻融与盐蚀劣化作用解耦路径示意图。图中路径a-b表示在5次冻融循环后仅冻融作用引起的劣化值；路径b-c表示5次冻融循环后,含盐量0.5%时仅盐蚀作用引起的劣化值；路径b-d表示5次冻融循环后,含盐量1.0%时仅盐蚀作用引起的劣化值；路径b-e表示在5次冻融循环后,含盐量1.5%时仅盐蚀作用引起的劣化值。
为定量化揭示任意损伤过程中冻融与盐蚀劣化作用的贡献,采用归一化处理方法分别定义冻融与盐蚀劣化因子,可通过下式计算：
D1=c0−ci−0c0−ci−j, (1) D2=ci−0−ci−jc0−ci−j。 (2) 式中 D1为冻融劣化因子,亦即冻融损伤占比；D2为盐蚀劣化因子,亦即盐蚀损伤占比；
c0 表示初始状态试样黏聚力；ci−0 表示冻融循环次数为i、含盐量为0.0 %时黏聚力；ci−j 表示冻融循环次数为i、含盐量为j %时黏聚力。（其中i取1,2,5,10；j取0.5,1,1.5）。进一步由式（1）,（2）,可推得以下关系式：
D1+D2=1。 (3) 式（3）表明,任意损伤过程归一化处理后,冻融和盐蚀劣化因子D1,D2即分别定量化揭示了冻融和盐蚀损伤的贡献,其物理意义直观明了。
图8（a）所示为冻融与盐蚀劣化因子的比值随冻融循环次数变化规律。从图中可以看出,冻融与盐蚀劣化因子的比值随冻融次数增加逐渐增大,亦即冻融作用对试样强度劣化的影响逐渐增强；但随着冻融次数持续增大,其增速逐渐减小,趋于一个稳定数值。值得注意的是,冻融作用下含盐量较高试样（1.5%）冻融与盐蚀劣化因子的比值始终小于1,亦即冻融劣化效应弱于盐蚀劣化效应。分析其原因,由前述盐蚀劣化因子D2的计算公式可知,对于含盐量较高试样,其计算值显著增大,从而导致冻融与盐蚀劣化因子的比值急剧减小。
![]()
图 8 冻融与盐蚀劣化因子比值的变化规律曲线Figure 8. Variation of ratio of freeze-thaw deterioration factor to that by salt erosion冻融与盐蚀劣化因子的比值随含盐量变化规律如图8（b）所示。由图可见,冻融与盐蚀劣化因子的比值随含盐量增大逐渐减小,亦即盐蚀作用对试样强度劣化的影响逐渐增强；随着含盐量持续增大,其衰减速率逐渐减小,趋于一个较小数值。此外,仅冻融循环1次条件下,冻融劣化因子与盐蚀劣化因子的比值始终小于1,亦即冻融劣化效应弱于盐蚀劣化效应。分析其原因,由前述冻融劣化因子D1的计算公式可知,对于仅冻融循环1次试样,其计算值显著减小,从而导致冻融与盐蚀劣化因子的比值亦显著减小。
3.2 CT扫描试验结果分析
（1）冻融循环次数对细观结构的影响规律
图9所示为Na2SO4含量1.0%试样冻融过程不同截面的CT细观结构图像。图中白色区域代表试样密度较大处,黑色区域代表试样中孔洞、裂隙发育的位置。从图中可以看出,冻融循环条件下不同扫描断面的细观结构损伤表现出相似的变化规律：冻融前黄土试样内部存在一定的孔洞及微裂纹且无序分布,从而导致试样产生一定的初始损伤效应。随着冻融循环次数增加,受冰水相变的影响,微小孔洞逐渐扩大,且有互相贯通之势；初始微裂纹逐渐发展成裂隙,其长度及宽度均有不同程度扩展。
图10所示为CT数ME值随冻融循环次数变化规律曲线。由图可见,冻融初期阶段试样的CT数ME值衰减速率较大,随着冻融次数持续增大,CT数ME值逐渐趋于稳定,表现出减速衰减特征。上述变化规律反映出冻融初始阶段试样细观结构损伤演化速率较大,亦即试样内部微裂隙增长较快；多次冻融循环后,试样细观结构损伤劣化作用减弱。这很好的解释了前述三轴破坏偏应力和黏聚力随冻融循环次数的变化关系（图4（a）,5（a））。
（2）Na2SO4含量对细观结构的影响规律
图11所示为冻融循环5次条件下不同Na2SO4含量盐渍原状黄土试样c截面CT细观结构扫描图像。由图可见,冻融作用下不含盐试样内部大孔隙或裂隙发育相对较少。随着含盐量增加,CT图像灰度有一定程度的加深,特别是图像下部裂隙有一定的扩展和发育,出现新的细小孔隙及裂隙,盐蚀劣化程度较不含盐试样显著增大；对含盐量为1.5%试样,裂缝扩展程度较大,且中部产生较大孔隙。
图12所示为CT数ME值随含盐量变化曲线。由图可见,冻融条件下CT数ME值随含盐量增大呈现出近似线性或加速衰减特征。上述变化规律反映出冻融条件下试样细观结构损伤演化速率随含盐量增大有增大的趋势,亦即试样内部微裂缝扩展速率逐渐增大,盐蚀劣化作用增强。这亦很好的解释了前述三轴破坏偏应力和黏聚力随含盐量的变化规律（图4（b）,5（b））。值得注意的是,未经受冻融循环（N=0）试样的CT数ME值随含盐量增加无显著变化,这是由于未冻融条件下试样内部Na2SO4易溶盐无相态变化,盐蚀劣化作用无法产生,因而CT数ME值无明显变化。
（3）CT细观损伤变量DME
根据CT原理[14]可以得到
ρ=ρ0(1000+H)/(1000+H0)。 (4) 式中 ρ为试样的密度；ρ0为初始状态下试样的密度；H,H0分别为与ρ和ρ0所对应的试样CT数ME值。
令
Δρ=ρ−ρ0, (5) 将式（4）代入式（5）,则有
Δρ=(1000+H1000+H0−1)×ρ0。 (6) 根据密度损伤变量公式[15],可以得到CT细观结构损伤变量DME的表达式如下：
DME=−1m20Δρρ0, (7) 式中,m0为CT机空间分辨率。
将式（6）代入式（7）,则可以最终得到基于CT数ME值的损伤变量DME表达式：
DME=1m20×(H0−H1000+H0), (8) 式中,DME数值越大表示试样损伤幅值越大,DME=0代表试样初始状态。
图13所示为CT细观损伤变量DME变化规律。从图中可以看出,细观损伤变量随着冻融循环次数增加逐渐增大,但增速逐渐减缓,表现出减速劣化特性；细观损伤变量随着含盐量增加近似线性或加速增大,表现出等速或加速劣化特性。
进一步根据细观损伤变量DME随冻融循环次数变化曲线,发现损伤变量DME与冻融循环次数N具有如下关系：
DME(N)=Na+bN, (9) 式中,N为冻融循环次数,a,b均为拟合参数。
将损伤变量DME随冻融循环次数变化结果按N/DME–N的关系进行拟合分析,结果如图14（a）所示。从图中可以看出,二者近似呈线性关系,其中拟合参数a为直线的截距,b为直线的斜率。
对式（9）进一步求导可得
d(DME)dN=a(a+bN)2。 (10) 在曲线的起始点,N=0,则式（10）可表示为
IDME=1a, (11) 式中,
IDME 即为损伤变量DME的初始斜率。当N→∞时,从式（9）中可得出
(DME)ult=1b, (12) 式中,
(DME)ult 为损伤变量DME的极限值,表示试样达到冻融破坏极限时损伤变量的峰值。由此可以看出,a代表损伤变量
DME 的初始斜率IDME 的倒数；b为损伤变量DME 极限值(DME)ult 的倒数。图14（b）所示为拟合参数与含盐量变化关系。从图中可以看出,拟合参数a,b随含盐量可近似表示为线性变化关系。因此,进一步考虑含盐量的影响,建立拟合参数a,b与含盐量关系式：
a=αη+a0, (13a) b=βη+b0, (13b) 式中,η为试样含盐量（%）；a0,b0和α,β分别为图中直线的截距与斜率。
将式（13）代入式（9）中,可得到细观损伤变量DME在不同冻融循环次数及含盐量下的多变量演化方程：
DME(N,η)=N(αη+a0)+(βη+b0)N, (14) 式中,所有参数均通过试验结果拟合得到,a0=4.941,b0=3.715,α=-0.915,β=-1.087。
图14（c）是含盐量分别为0%,0.5%,1.0%,1.5%下的损伤变量试验值与计算值对比图。从图中可以看出,损伤变量试验值与计算值均匀分布于直线y=x两侧,拟合相关性较好,表明该模型可较好预测Na2SO4盐渍原状黄土冻融过程细观结构损伤演化规律。
4. 讨论
Na2SO4盐渍原状黄土宏观损伤特征表现为强度劣化,劣化过程中必然伴随着试样内部细观裂隙的产生和发展,CT扫描试验可以定量化研究试样细观结构损伤演化规律。基于此,分别定义基于黏聚力的宏观损伤变量λC和CT数ME值的细观损伤变量DME,以对比分析宏细观损伤演化规律的相互关系。CT细观损伤变量
DME 已通过前述式（8）得到,基于黏聚力宏观损伤变量λC可通过下式计算：λC=c0−ci−jc0, (15) 式中,λC数值越大表示试样损伤幅值越大,λC=0代表试样初始状态,λC=1代表完全损伤状态；
c0 ,ci−j 的含义如前述公式（1）和（2）所示。宏观损伤变量λC和细观损伤变量DME变化规律如图15所示。由图可见,宏细观损伤变量表现出相似的变化规律：冻融初始阶段宏细观损伤变量均增速较大,随着冻融次数持续增加,逐渐趋于稳定,表明冻融作用会导致试样的减速劣化效应。此外,宏细观损伤变量随含盐量增加其增速均逐渐增大,表明盐蚀作用会导致试样的加速劣化效应。综上,宏细观损伤变量表现出的一致变化规律表明CT细观结构损伤变量准确揭示了宏观三轴剪切力学强度指标的劣化机理。
![]()
图 15 宏细观损伤变量变化规律曲线Figure 15. Variation of macroscopic and mesoscopic damage variables5. 结论
（1）冻融作用对应力–应变曲线的类型及特征无明显影响,均表现为应变硬化型。破坏偏应力随着冻融次数增加逐渐减小,但降低速率逐渐减小,表现出减速劣化特征；冻融条件下破坏偏应力随着含盐量增加逐渐减小,且表现出线性或加速劣化特征。
（2）黏聚力呈现与破坏偏应力相似的劣化特征；内摩擦角变化幅值较小,无显著变化规律。
（3）冻融与盐蚀劣化因子的比值随冻融次数增加逐渐增大但增速逐渐减小,随含盐量增大逐渐减小且衰减速率逐渐减小。
（4）CT数ME值表现出与破坏偏应力和黏聚力相似的劣化规律；构建了CT细观损伤变量演化方程,可较好预测冻融过程细观结构损伤演化规律。
（5）宏细观损伤变量表现出一致的变化规律,表明CT细观结构损伤变量准确揭示了宏观三轴剪切强度指标的劣化机理。
-
![]()
图 6 等应力幅循环剪切加载过程
Figure 6. Loading program for uniform-amplitude cyclic triaxial compression tests
![]()
图 7 微生物处理砂土试样循环应力比-液化振次曲线
Figure 7. Cyclic stress ratio-cycles to initial liquefaction curves of sand samples treated by microorganism
![]()
图 8 胶结与气泡共同作用与单一作用砂土试样液化线对比（箭头代表试样在该振次下仍未破坏）
Figure 8. Comparison of liquefaction lines of sand samples under single and combined actions of cementation and bubbles (arrows represents that sample is still intact under this vibration)
![]()
图 9 胶结与气泡共同作用和单一作用砂土试样孔压比对比
Figure 9. Pore pressure ratios of sand samples under single and combined actions of cementation and bubbles
![]()
图 10 胶结与气泡共同作用与单一作用试样轴向应变对比
Figure 10. Axial strains of samples under single and combined actions of cementation and bubbles
![]()
图 11 胶结与气泡共同作用和单一作用试样力学配位数对比
Figure 11. Mechanical coordination number of samples under single and combined actions of cementation and bubbles
表 1 胶结颗粒模型微观参数
Table 1 Microscopic parameters of model for bond grains
组成 参数名称 数值 颗粒部分 颗粒模量Ep/Pa 3×108 颗粒法向与切向刚度比 ξ 1.5 颗粒摩擦系数 μ 0.5 颗粒抗转动系数 β 0.2 颗粒局部破损系数 ξc 2.1 颗粒法向接触加卸载刚度比 en 0.95 胶结部分 胶结半径比例系数 η (0.5)* 胶结模量Eb/Pa 2×108 胶结泊松比 νb 0.2 胶结抗拉强度 σt /Pa1.6×106 胶结抗压强度 σc /Pa8×107 黏滞阻尼 法向临界阻尼比 βn 0.6 切向临界阻尼比 βs 0.6 注： 胶结半径比例系数η 在标定过程中取0.5,在此后的模拟中为控制胶结含量的研究变量。
下载: 导出CSV
表 2 数值模拟试样微生物处理过程时效性
Table 2 Time-dependent characteristics of microbial processing of numerical simulation samples
反应时间/d 微生物处理生成胶结 微生物处理生成气泡 胶结体积/mm3 胶结半径比例系数η 气泡体积/mm3 饱和度Sr/% 孔隙流体压缩系数C 0 0 0 0 100 4.5×10-7 2 0.718 0.231 2.46 98.3 3.87×10-5
下载: 导出CSV
表 3 微生物处理砂土循环三轴数值模拟工况与结果列表
Table 3 Simulation program and main results of cyclic triaxial compression tests on different microbial treatment time on sand
试样编号 模拟工况 模拟结果 土体种类 微生物处理时间/d 胶结半径比例系数η 孔隙流体压缩系数C 循环应力比CSR 破坏振次NL 破坏模式 L0 净砂 0 0 0.15 1 I.L. MU2-1 非饱和MICP胶结砂土 2 0.273 3.87×10-5 0.30 45+ — MU2-6 非饱和MICP胶结砂土 2 0.273 3.87×10-5 0.60 3 I.L. MP2-2 MICP胶结砂土 2 0.273 0.60 5 I.L. US2-4 非饱和砂土 2 3.87×10-5 0.30 2 I.L.
下载: 导出CSV
-
[1] 沈珠江. 理论土力学[M]. 北京: 中国水利水电出版社, 2000. SHEN Zhu-jiang. Theoretical Soil Mechanics[M]. Beijing: China Water and Power Press, 2000. (in Chinese)
[2] GAO Y F, HANG L, HE J, et al. Mechanical behaviour of biocemented sands at various treatment levels and relative densities[J]. Acta Geotechnica, 2019, 14(3): 697-707. doi: 10.1007/s11440-018-0729-3
[3] WANG Y, SOGA K, DEJONG J T, et al. A microfluidic chip and its use in characterising the particle-scale behaviour of microbial-induced calcium carbonate precipitation (MICP)[J]. Géotechnique, 2019, 69(12): 1-9.
[4] INAGAKI Y, TSUKAMOTO M, MORI H, et al. A centrifugal model test of microbial carbonate precipitation as liquefaction countermeasure[J]. Japanese Geotechnical Journal, 2011, 6(2): 157-167. doi: 10.3208/jgs.6.157
[5] HAN Z G, CHENG X H, MA Q. An experimental study on dynamic response for MICP strengthening liquefiable sands[J]. Earthquake Engineering and Engineering Vibration, 2016, 15(4): 673-679. doi: 10.1007/s11803-016-0357-6
[6] 刘汉龙, 肖鹏, 肖扬, 等. MICP胶结钙质砂动力特性试验研究[J]. 岩土工程学报, 2018, 40(1): 38-45. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201801003.htm LIU Han-long, XIAO Peng, XIAO Yang, et al. Dynamic behaviors of MICP-treated calcareous sand in cyclic tests[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(1): 38-45. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201801003.htm
[7] REBATA-LANDA V, SANTAMARINA J C. Mechanical effects of biogenic nitrogen gas bubbles in soils[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2012, 138(2): 128-137. doi: 10.1061/(ASCE)GT.1943-5606.0000571
[8] HE J. Mitigation of liquefaction of sand using microbial methods[D]. Singapore: Nanyang Technological University, 2013.
[9] CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Géotechnique, 1979, 29(1): 47-65. doi: 10.1680/geot.1979.29.1.47
[10] KHOUBANI A, EVANS T M, MONTOYA B M. Particulate simulations of triaxial tests on bio-cemented sand using a new cementation model[C]//Proceedings of GeoChicago: Sustainability, Energy, and the Geoenvironment. Chicago, 2016: 1-10.
[11] FENG K, MONTOYA B M, EVANS T M. Discrete element method simulations of bio-cemented sands[J]. Computers and Geotechnics, 2017, 85: 139-150. doi: 10.1016/j.compgeo.2016.12.028
[12] ANDERSON J D, WENDT J. Computational Fluid Dynamics[M]. New York: McGraw-Hill, 1995.
[13] CHANG C Y, SCHMIDT J, DÖRENKÄMPER M, et al. A consistent steady state CFD simulation method for stratified atmospheric boundary layer flows[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2018, 172: 55-67. doi: 10.1016/j.jweia.2017.10.003
[14] PAREKH J, RZEHAK R. Euler-Euler multiphase CFD-simulation with full Reynolds stress model and anisotropic bubble-induced turbulence[J]. International Journal of Multiphase Flow, 2018, 99: 231-245. doi: 10.1016/j.ijmultiphaseflow.2017.10.012
[15] EL SHAMY U, ZEGHAL M. A micro-mechanical investigation of the dynamic response and liquefaction of saturated granular soils[J]. Soil Dynamics and Earthquake Engineering, 2007, 27(8): 712-729. doi: 10.1016/j.soildyn.2006.12.010
[16] 蒋明镜, 张望城. 一种考虑流体状态方程的土体CFD-DEM耦合数值方法[J]. 岩土工程学报, 2014, 36(5): 793-801. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201405002.htm JIANG Ming-jing, ZHANG Wang-cheng. Coupled CFD-DEM method for soils incorporating equation of state for liquid[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(5): 793-801. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201405002.htm
[17] 蒋明镜. 现代土力学研究的新视野—宏微观土力学[J]. 岩土工程学报, 2019, 41(2): 195-254. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201902002.htm JIANG Ming-jing. New paradigm for modern soil mechanics: Geomechanics from micro to macro[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(2): 195-254. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201902002.htm
[18] JIANG M J, SHEN Z F, WANG J F. A novel three- dimensional contact model for granulates incorporating rolling and twisting resistances[J]. Computers and Geotechnics, 2015, 65: 147-163. doi: 10.1016/j.compgeo.2014.12.011
[19] 金树楼. 结构性砂土三维微观接触力学试验及离散元数值模拟[D]. 上海: 同济大学, 2016. JIN Shu-lou. Three Dimensional Experimental and Numerical Study on Micro- and Macro-mechanical Behaviors of Structural Sands[D]. Shanghai: Tongji University, 2016. (in Chinese)
[20] THORNTON C, CUMMINS S J, CLEARY P W. An investigation of the comparative behavior of alternative contact force models during inelastic collisions[J]. Powder Technology, 2013, 233(3): 30-46.
[21] SHEN Z F, JIANG M J, THORNTON C. DEM simulation of bonded granular material: Part I contact model and application to cemented sand[J]. Computers & Geotechnics 2016, 75: 192-209.
[22] 谭亚飞鸥. 考虑循环荷载的三维微观胶结模型及微生物处理砂土循环三轴CFD-DEM耦合模拟[D]. 上海: 上海理工大学, 2018. TAN Ya-fei-ou. A Novel Three-dimensional Bonded Contact Model Incorporating the Effect of Cyclic Loads and CFD-DEM Simulation of Microbially Treated Sands Under Undrained Consolidated Cyclic Triaxial Tests[D]. Shanghai: University of Shanghai for Science and Technology, 2018. (in Chinese)
[23] ZHANG A, JIANG M J. Numerical simulation of undrained triaxial tests for granular soil using a coupled CFD-DEM method with moving mesh[J]. Acta Geotechnica, 2019. (to be submitted)
[24] ZHAO T. Investigation of Landslide-induced Debris Flows by the DEM and CFD[D]. Oxford: University of Oxford, 2014.
[25] VAN PAASSEN L A, DAZA C M, STAAL M, et al. Potential soil reinforcement by biological denitrification[J]. Ecological Engineering, 2010, 36(2): 168-175.
[26] AL QABANY A, SOGA K, SANTAMARINA C, et al. Factors affecting efficiency of microbially induced calcite precipitation[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2012, 138(8): 992-1001.
[27] SCHUURMAN I E. The compressibility of an air/water mixture and a theoretical relation between the air and water pressures[J]. Géotechnique, 1966, 16(4): 269-281.
[28] 刘侃, 朱小军, 张帆舸, 等. 含气泡土的孔隙流体压缩系数计算分析[J]. 岩土工程学报, 2017, 39(增刊2): 120-123. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2017S2031.htm LIU Kan, ZHU Xiao-jun, ZHANG Fan-ge, et al. Calculation of coefficient of compressibility for air-water mixture in gassy soils[J]. Chinese Journal of Geotechnical Engineering,2017, 39(S2): 120-123. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2017S2031.htm
[29] DORSEY N E. Properties of Ordinary Water-Substance[M]. New York: Reinhold Publishing Corporation, 1940.
[30] FENG K, MONTOYA B M. Influence of confinement and cementation level on the behavior of microbial-induced calcite precipitated sands under monotonic drained loading[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016, 142(1): 04015057.
[31] JIANG M J, KONRAD J M, LEROUEIL S. An efficient technique for generating homogeneous specimens for DEM studies[J]. Computers and Geotechnics, 2003, 30(5): 579-597.
[32] KUHN M R, RENKEN H E, MIXSELL A D, et al. Investigation of cyclic liquefaction with discrete element simulations[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2014, 140(12): 04014075.
-
期刊类型引用(16)
1. 张志超，游志浪，许健，彭洋，GAGARIN L. 莫斯科—喀山高速铁路沿线大气冻融指数时空变化特征. 长沙理工大学学报(自然科学版). 2025(01): 15-27 . 
百度学术
2. 孔元元，梁晨涛，王智慧，谢柏涵，成忠焱. 考虑水盐运移的分散性盐渍土冻融破坏机制. 长安大学学报(自然科学版). 2024(02): 1-11 . 百度学术
3. 包卫星，吴倩，吴谦，秦川，侯天琪. 冻融循环作用下伊犁盐渍化黄土力学特性. 岩石力学与工程学报. 2024(07): 1775-1787 . 百度学术
4. 武立波，杨嘉伟，刘惠阳，陈宏信. 煤气化粗渣改良季节冻土区黄土填料的冻融特性及其微观机理研究. 冰川冻土. 2024(03): 930-943 . 百度学术
5. 樊培培，张凌凯，丁旭升. 干湿-冻融循环条件下湿陷性黄土剪切及压缩特性的劣化规律. 岩土力学. 2024(07): 2050-2060 . 百度学术
6. 包卫星，秦川，李伟，刘亚伦，陈锐. 冻融循环下低液限盐渍化粉土力学特性研究. 地下空间与工程学报. 2024(04): 1223-1235 . 百度学术
7. 张卫兵，李晓，雷过，刘臻祥，周鑫磊. 冻融—干湿循环下硫酸盐渍土的微观孔隙研究. 地下空间与工程学报. 2023(02): 465-473 . 百度学术
8. 王晓强，李明玉，高富强. 硫酸钠盐渍黄土物理力学特性演变规律及机理研究. 公路. 2023(06): 392-396 . 百度学术
9. 孙杰龙，王弘起，李盛斌，李大卫，邱明明. 冻融作用下高填方黄土抗剪强度劣化特性分析. 岩土工程技术. 2023(05): 609-613 . 百度学术
10. 王昊，韦金城，宋晓辉，李雪和，吴文娟，刘海峰，任瑞波. 黄泛区粉土工程特性及其改性固化研究进展. 粉煤灰综合利用. 2023(05): 107-119 . 百度学术
11. 李品良，许强，刘佳良，何攀，纪续，陈婉琳，彭大雷. 盐分影响重塑黄土渗透性的微观机制试验研究. 岩土力学. 2023(S1): 504-512 . 百度学术
12. 包卫星，李伟，毛雪松，陈锐，秦川，刘亚伦. 冻融循环条件下盐渍化风积沙力学特性. 交通运输工程学报. 2023(06): 114-124 . 百度学术
13. 王亮亮，王照腾，方薇，田建胜. 气候边缘地带膨胀土强度特性随冻融循环劣化规律. 中南大学学报(自然科学版). 2022(01): 288-295 . 百度学术
14. 张莎莎，刘亚超，杨晓华，陈伟志，金好乾. 粗粒盐渍土区既有高铁路基变形特性试验研究. 建筑科学与工程学报. 2022(02): 135-142 . 百度学术
15. 刘华，胡鹏飞，王梦南，刘乃飞，胡文乐. 冻融循环对酸污染黄土抗拉特性劣化试验研究. 西安建筑科技大学学报(自然科学版). 2021(04): 493-501 . 百度学术
16. 程卓，崔高航，高原昊，刚浩航，高泽宁，杨政，张鑫. 季冻区粉煤灰加固路基土力学性能试验研究. 硅酸盐通报. 2021(11): 3854-3864+3875 . 百度学术
其他类型引用(20)
计量
- 文章访问数: 1168
- HTML全文浏览量: 68
- PDF下载量: 757
- 被引次数: 36
