Fracture failure mechanism and fracture criterion of compacted clay under compression and shear action
-
摘要: 土体既有裂缝的开裂问题通常被认为是拉剪应力导致,而实际工程中压实黏土多处于压剪应力状态,为探究压剪作用下压实黏土中既有裂缝的起裂机理,对压剪闭合裂缝应力场和断裂准则进行了理论分析,开展了含中心裂缝压实黏土单轴压缩试验,研究了裂缝倾角、裂缝长度及裂缝形态对压剪断裂性状的影响,揭示了压实黏土中闭合裂缝压剪张拉断裂机理,充分验证了闭合裂缝压剪张拉断裂准则在压实黏土中的适用性。根据试验结果与理论预测曲线对比分析,本文压实黏土的临界尺寸rc为2 mm左右,明显小于常用经验公式估算值12 mm,讨论了rc经验估算公式存在的问题。此外,尝试性地引入T应力,建立了非闭合裂缝的压剪-张拉断裂准则,解释了闭合裂缝和非闭合裂缝的差异化起裂行为,并说明了其局限性及其产生原因。Abstract: Soil cracks are often considered to be caused by tensile-shear stress. In actual engineering, the compacted clay is mostly under compressive-shear stress state. In order to explore the crack initiation mechanism of the existing cracks in the compacted clay under compressive-shear stress state, the stress field and fracture criterion of the closed crack are analyzed theoretically. The uniaxial compression tests on the compacted clay with central cracks are carried out. The influence of dip angle, length and shape of cracks on the tensile fracture behavior of the compacted clay are investigated. The applicability of the compression-shear closed-crack tension fracture criterion considering T-stress in the compacted clay is verified, and the compression-shear tension fracture mechanism of closed cracks in the compacted clay is revealed. According to the comparative analysis of the test results and the theoretical prediction curve, the critical size rc of the compacted clay is about 2 mm, which is obviously smaller than the estimated value of 12 mm by the empirical formula used in the rock field. Through detailed discussion, it is shown that the formula applied in the rock field is not suitable for the compacted clay. In addition, the T-stress is introduced tentatively to establish a compression-shear-tension fracture criterion for open cracks. The differential initiation behavior of closed and open cracks is explained, and its limitations are shown.
-
Keywords:
- compressive shear stress /
- closed crack /
- fracture criterion /
- compacted clay /
- T-stress
-
0. 引言
随着中国对环境和生态的日益重视,全国主要的河流湖泊基本都已开展了底泥清淤和疏浚,例如太湖第一轮清淤已经将3600万m3污染底泥全部清除,而正在开展的太湖第二轮清淤,疏浚工程量仍然高达1800万m3。根据中国的底泥疏浚国情,固化方法仍然是处理大量疏浚淤泥的主要方法[1-3]。目前中国的淤泥固化处理主要以低含水率水平下固化再进行适当碾压为主,控制含水率在较低水平主要是为了减少固化材料的用量,但这种碾压方式的固化淤泥施工,经常会受限于固化设备的处理能力,并且碾压后固化淤泥的强度会有很大损失,不利于大规模固化淤泥施工[4-6]。
国外例如日本、美国,除了上述碾压方式的固化淤泥常用作筑堤以外[7-8],对于造岛、机场跑道等大规模的淤泥固化处理工程,更多的是采用流动固化的方式[9],该方法采用高浓度泥浆输送,经过调理以后加入固化剂,固化剂在管道输泥的过程中经过气流混合形成流动性的固化淤泥,直接浇筑在场地内形成人工岛[10]。此种方法从清淤到固化浇筑主要都在管道内完成,效率高,固化淤泥一次成型,因此强度大、变形小,是未来大规模淤泥固化处理技术的发展方向[11]。
国内对固化淤泥的流动性研究已有不少学者开始关注,丁建文等[12]研究了不同初始含水率和不同固化材料添加量下固化淤泥的流动性,建立了流动值的预测方法。黄英豪[11]、朱伟等[13]研究了淤泥固化之后流动性的变化模式,并提出了固化淤泥的流动浇筑和碾压施工的选择方法。顾欢达等[9]通过调整淤泥含水率,使得固化淤泥保持一定的流动性,从而实现淤泥的回填或者回灌处理的技术要求。陆强[14]研究了海水环境下固化淤泥的流动性。纪文栋等[15]采用高吸水材料来改善淤泥的流动性,但其目的是为了减少流动性以便于淤泥的运输和周转。
本文采用两种疏浚淤泥,研究了淤泥加入水泥固化之后流动性和黏滞性的变化,并建立了流动性和黏滞性之间的内在关系,以期对淤泥的流动固化施工技术发展有所助益。
1. 试验材料和方法
1.1 试验材料
研究所用淤泥分别取自淮安白马湖和太湖梅梁湖,淤泥的主要物理性质指标如表 1所示,根据土的分类方法,均属于高液限黏土(CH)。考虑到目前大规模淤泥固化处理采用的固化材料主要为水泥,为了减少影响因素,本文只采用#32.5普通硅酸盐水泥,没有和其他固化材料进行复合。
表 1 试验用淤泥物理性质指标Table 1. Physical properties of test soils淤泥名称 天然含水率wn/% 密度ρ/(g·cm-3) 相对质量密度Gs 黏粒含量/% 液限wL/% 塑限wP/% 有机质含量/% 白马湖 98 1.48 2.65 47.4 76.5 30 0.4 梅梁湖 80 1.53 2.68 26.3 56 27 3.5 1.2 试验方法
流动性测试采用的是新鲜搅拌的固化淤泥,即刚搅拌完成,水泥尚未完成初凝的固化淤泥,本文称之为新拌固化淤泥。试验时先按照设计配比将固化材料和淤泥在容积为6 L的砂浆搅拌机中强制搅拌10 min,使其混合均匀,刚搅拌完成的固化淤泥立即进行流动性和黏滞性的测试。
固化淤泥流动性测试目前在中国尚没有标准方法,学者多借鉴日本的圆柱筒容器测量水平扩展度。本文为了使得淤泥的流出更加顺畅,借鉴采用水泥砂浆流动性测试所采用的的截锥圆模装置,截锥圆模上底直径7 cm,下底直径10 cm,高6 cm。试验时把截锥模放置在水平玻璃板上,截锥模内壁和玻璃板上涂薄层凡士林,将新拌的固化淤泥装入到截锥圆模中,装入的过程不进行振捣,完全靠自身重力作用下流动充满至圆模,用刮泥刀削平截锥模上表面溢出的固化淤泥,将截锥圆模垂直向上轻轻提起,测定固化淤泥在互相垂直两个方向的流动直径,取其平均值作为新拌固化淤泥的流动度。
在流体力学中,流体的流变特性一般是采用旋转黏度计来测试,由于流体力学中的流变特性和土力学中常讲的土体的流变特性含义不尽相同,为了避免本专业人员的误解,本文采用黏滞特性一词来代替流变特性。黏滞性测试采用NXS-11A型旋转黏度计,它是一种通用的同轴圆筒上旋式黏度计,可用来测量各种牛顿流体和非牛顿流体的黏滞特性,为了测量不同黏度的的物质,仪器带有5个转子,选择不同的上旋转子和外筒及转速,可测量黏度不同黏度和不同剪切应力范围的的流体参数,如图 1所示。
为了考察不同含水率下新拌固化淤泥的流动性和黏滞性,将白马湖淤泥的初始含水率分别调至液限的不同倍数,具体为76%,98%,118%,153%,190%,228%,296%,456%,912%;将梅梁湖淤泥初始含水率分别调至56%,65%,87%,99%,130%,163%,343%,674%,856%。为了考察水泥添加量对流动性和黏滞性的影响,设定50,100,200 kg/m3三个水泥添加量。
2. 新拌固化淤泥的流动性
2.1 淤泥的流动性
这里含水率用w表示,并定义流动度D为淤泥或者新拌固化淤泥坍落稳定后在水平面内投影的直径(单位为cm,精确到mm)。图 2给出了梅梁湖淤泥典型的流动性的试验照片,(白马湖淤泥的照片类似,限于篇幅不再展示),对两种淤泥的流动性试验结果进行分析,发现有以下共同规律:①当含水率非常高,w≥6.0wL左右时,虽然还未提起截锥模,但已经有少量淤泥从截锥模的下底面流出,这主要是因为淤泥的稠度过小;②当含水率很高,2.5wL≤w<6.0wL时,流动稳定后淤泥的形态呈现非常规则的薄圆饼状,并且竖直方向的高度要远小于水平方向的扩展度;③当含水率较高,2.0wL≤w<2.5wL时,流动稳定后淤泥的形态虽然仍呈现圆饼状,但竖直方向的高度相比于前一阶段要显著增大;④当含水率较低,wL≤w<2.0wL时,淤泥的流动性非常小,流动稳定后的淤泥形态呈现圆台状;⑤当含水率很低,w<wL时,淤泥已经不具有流动性,并且从定义看,含水率小于液限的土体已经不属于淤泥的范畴。
图 3所示为白马湖淤泥和梅梁湖淤泥的流动度随含水率的变化,可见,两种泥的流动度随含水率降低都急剧减小,并且在相同含水率下两种泥的流动度相差不大,表明含水率是决定淤泥流动度的重要因子。另一方面,流动度随含水率变化的过程可以比较明显地分为两段,对白马湖淤泥当含水率w < 150%时,淤泥的流动度随着含水率的增加有很少量的增大,而当含水率w > 150%时,淤泥的流动度随着含水率的增加迅速增大。对梅梁湖淤泥,也同样存在这样的临界点,只不过临界含水率w′在110%左右。
如果将淤泥的含水率转化为其相应的液限倍数,从流动度随液限倍数的变化图来看,则流动度的变化规律更为明显,如图 3(b)所示。可以看出,无论对于白马湖淤泥还是梅梁湖淤泥,在w=2.0wL左右时,淤泥的流动度变化速率开始呈现差别。当w<2.0wL时,两种泥的流动度相差很小,并且随着含水率的变化改变不明显。当w>2.0wL时,两种泥的流动度都随着含水率增大而显著增加,但两种泥的流动度曲线并不重合,显示出液限对淤泥的流动度有一定影响,同样的液限倍数时,梅梁湖淤泥流动度要小于白马湖淤泥的流动度。这是因为白马湖淤泥的液限要高于梅梁湖淤泥的液限,使得液限倍数相同时白马湖淤泥的含水率要高于梅梁湖淤泥,后文将进一步从淤泥的黏滞性进行分析。
2.2 新拌固化淤泥的流动性
(1)初始含水率对新拌固化淤泥流动性影响
将两种淤泥调至不同含水率,分别加入50,100,200 kg/m3的固化材料,考察流动性的变化。如图 4所示,可以看出在水泥添加量相同时,初始含水率越高,新拌固化淤泥的流动度越大,和未加入水泥处理的原泥的流动度相比,加入水泥后的新拌固化淤泥的流动度要有很大程度的减小。这说明水泥的加入能显著增大淤泥的稠度,减小淤泥的流动性能。从流动度的变化趋势上来看,新拌固化淤泥的流动度变化也可以分为两段,只不过流动度变化的临界含水率点相比于原泥有明显增大,对于新拌白马湖固化淤泥,临界含水率w′s约为220%,对于新拌梅梁湖固化淤泥,临界含水率w′s约为150%,小于临界含水率时流动度变化很小。
如果将初始含水率用液限的倍数来表示,则新拌固化淤泥的流动性变化规律则更加明显,见图 4(c)所示。可以看出两种新拌固化淤泥的流动度都大致在同一临界点发生变化,该临界点的值约为w′s=2.8wL。小于该临界点时,两种固化淤泥的流动度相差很小,超过临界值之后,流动度迅速增大,并且同样的液限倍数下新拌白马湖固化淤泥的流动度要大于梅梁湖固化淤泥的流动度。
由前文分析可知,对未处理的淤泥,当w<2.0wL时,两种淤泥的流动度变化缓慢且相差不大,而对于新拌的固化淤泥,当w<2.8wL时,两种新拌固化淤泥的流动度变化缓慢且也相差不大。如果把两种原泥的流动度和水泥量为50 kg/m3的新拌固化淤泥的流动度随初始含水率的变化对比,如图 5所示。由图 5可见,所有的原泥的流动度都近似在一条直线上,而所有的新拌固化淤泥的流动度都在淤泥流动度的下方,并且也大致成一直线。对淤泥和新拌固化淤泥的流动度结果进行回归分析,可得如下关系:
D=11.34w(ac=0), (1) D=7.61w(ac=50)。 (2) 由式(1),(2)可知,对于流动度比较大的淤泥和新拌固化淤泥,流动度和初始含水率间都有很简单的线性正相关关系,比较而言,淤泥的流动度随初始含水率的变化速率要高于新拌固化淤泥的流动度变化速率。即淤泥经过固化之后,它的流动度会显著降低,在大流动度的情况下,新拌固化淤泥的流动度约为淤泥流动度的0.67倍,即加入50 kg/m3的水泥后,淤泥的流动度会减小1/3。
(2)水泥添加量对新拌固化淤泥流动性影响
由前文分析可知,初始含水率对新拌固化淤泥的流动度有显著影响。为了进一步考察不同水泥添加量对新拌固化淤泥流动度的影响,分别对白马湖淤泥和梅梁湖淤泥添加不同的水泥量,试验结果如图 6所示。由图 6可见,水泥量对流动度的影响主要发生在添加量小于50 kg/m3的阶段,而超过此水泥量后,继续添加水泥对流动度的影响很小。表明只要有少量水泥的添加,就会显著降低淤泥的流动性能。
3. 新拌固化淤泥的黏滞性
淤泥固化后的流动度有很大程度的降低,为了进一步探究为何淤泥固化后的流动性会变化如此之大,本节从旋转黏度计测得的淤泥的黏滞性方面进一步研究。
3.1 流体黏滞特性的划分
流体流动时有速度梯度存在,因此运动较慢的液层阻滞着较快层的运动,因此产生流动阻力,切应力的作用就是克服这流动阻力,以维持一定的速度梯度而流动。表达剪切力和剪切速率之间关系的数学方程称为黏滞方程,它们之间的图解关系称为黏滞曲线。流体的黏滞曲线主要可以分为4个类型:牛顿流体、宾汉姆流体、幂律流体、带屈服值的幂律流体,具体见图 7。图 7中,τ为切应力(Pa);μ为液体的黏度(Pa·s);˙γ为剪切速率(s-1);τ0为剪切屈服应力,通常也可称为为动切力;n为流动指数,n值小于1,是非牛顿性的量度,与1相差越多,非牛顿性越强。
3.2 淤泥的黏滞特性
和前述测试淤泥流动度的试验方案相对应,分别对白马湖淤泥和梅梁湖淤泥的黏滞性测试,典型的试验结果如图 8所示(限于篇幅,仅只给出3幅白马湖淤泥的典型结果图)。综合分析两种淤泥的试验结果可以发现无论对于白马湖淤泥还是梅梁湖淤泥,在淤泥的含水率由高到底的变化过程中,淤泥的黏滞特性逐渐发生变化,可分为3种不同的类型。
(1)当淤泥的含水率w>12.0wL,淤泥的黏滞曲线基本为通过原点的直线,表明此时可以认为淤泥是牛顿流体,黏性特征不显著。这主要是因为淤泥的含水率非常高,土颗粒在水中分布形成悬液,土颗粒之间尚没有形成连接作用,淤泥的剪切应力的大小只与剪切速率相关。
(2)当淤泥的含水率2.0wL≤w≤12wL,淤泥的黏滞曲线为不通过原点的直线,表明在此范围内淤泥为宾汉姆流体。在这个阶段,溶液中的土颗粒数量增多,使得土颗粒之间能够形成连接作用,淤泥要实现流动,首先必须克服这种连接作用,也就是要克服初始动切力τ0的作用。
(3)当淤泥的含水率w<2.0wL,淤泥的黏滞曲线呈现出带屈服值的伪塑性流体的特征,剪切应力和剪切速率之间基本为幂函数关系。在这个阶段,淤泥的稠度已经很大,土颗粒间的作用非常强烈。
一般工程实际中,需要固化处理的淤泥含水率不会大于12倍液限含水率。疏浚工程中采用绞吸式挖泥船清淤出的泥浆的含水率多为液限的几倍。因此,工程实际中的刚刚吹填的淤泥一般属于宾汉姆流体。能表征宾汉姆流体特征的两个黏滞参数是动切力τ0和黏滞系数μ,对白马湖淤泥和梅梁湖淤泥,当淤泥属于宾汉姆流体时的动切力和黏滞系数进行分析,图 9为梅梁湖淤泥的动切力τ0和黏滞系数μ随淤泥初始含水率的变化。
由图 9可见,当淤泥属于宾汉姆流体时,它的黏滞系数都比较小,并且两种淤泥的黏滞参数随含水率的变化都呈现指数关系变化:
μ=ae−0.02w (3) τ0=be−0.02w (4) 式中,e为自然常数,w为淤泥的含水率,a和b分别为常数。由关系式可知,指数都为-0.02,表明两种淤泥的黏滞参数随着含水率的升高而降低的速率相同,但参数a和b则不同,说明了不同淤泥种类对黏滞系数的影响。白马湖淤泥的黏滞系数参数a为0.70,而梅梁湖淤泥的黏滞系数参数a为0.91;白马湖淤泥的动切力参数b为768,梅梁湖淤泥的动切力参数b为560。另一方面,白马湖淤泥的黏粒含量为47.4%,而梅梁湖淤泥的黏粒含量为26.3%,可见参数a和黏粒含量负相关,而参数b和黏粒含量呈现正相关。
从流动度的试验中知道,当淤泥的含水率w<2.0wL时,淤泥的流动度随着含水率的增加变化不大,而当淤泥的含水率w>2.0wL时,淤泥的流动度随着含水率的增加迅速增大。从淤泥的黏滞性试验中知道,当淤泥的含水率w<2.0wL时,淤泥属于带屈服值的幂律流体,而当淤泥的含水率w>2.0wL时,淤泥属于宾汉姆流体。这表明,当淤泥的黏滞特性属于宾汉姆流体时,它的流动性会很大,而当淤泥属于屈服值的幂律流体时,它的流动性很小。
3.3 新拌固化淤泥的黏滞特性
(1)初始含水率对新拌固化淤泥黏滞性影响
对新拌固化淤泥黏滞特性的研究分为两种情况:①固定水泥添加量为50 kg/m3,研究不同初始含水率下固化淤泥的黏滞特性;②在确定新拌固化淤泥的黏滞曲线类型之后,选择有代表性的初始含水率,研究不同水泥添加量对黏滞特性的影响。典型的白马湖新拌固化淤泥黏滞性试验结果见图 10。可见在水泥添加量为50 kg/m3的情况下,新拌固化淤泥的黏滞特性显著受到初始含水率的影响并且对于白马湖淤泥和梅梁湖淤泥有一样的规律。在含水率由高到底变化过程中,新拌固化淤泥的黏滞曲线可以显著的分为两个阶段、两种类型。
a)对白马湖淤泥,当初始含水率w≥228%;对梅梁湖淤泥,当淤泥的初始含水率w≥163%,固化淤泥的黏滞特性为宾汉姆流体。可以近似认为当淤泥的初始含水率w≥2.8wL时,新拌固化淤泥的黏滞特性表现为宾汉姆流体。
b)当淤泥的初始含水率w<2.8wL时,无论对于白马湖淤泥,还是对于梅梁湖淤泥,新拌固化淤泥的黏滞性都已经不符合宾汉姆流体,近似为带屈服值的幂律流体。
由新拌固化淤泥的流动度试验知道,新拌固化淤泥的流动度显著受到淤泥初始含水率的影响,并可以分为两个阶段。当w<2.8wL时,固化淤泥的流动度随含水率的增加变化不大;当w>2.8wL时,固化淤泥的流动度随着含水率的增加迅速增大。对比固化淤泥的黏滞性试验结果可以得出,当w>2.8wL时,新拌固化淤泥属于宾汉姆流体,这个阶段固化淤泥的流动性很大,并且初始含水率对固化淤泥的流动性影响很显著;当w<2.8wL时,新拌固化淤泥属于带屈服值的伪塑性流体,这个阶段固化淤泥的流动性很小。
当新拌固化淤泥属于宾汉姆流体时,黏滞系数μ和动切力τ0是表征其流动特性的两个参数,新拌白马湖固化淤泥的黏滞系数μ和动切力τ0随初始含水率的变化见图 11,原泥的宾汉姆流体参数也放在了一起以便于比较分析。由图 11可以看出,固化淤泥的动切力和黏滞系数都要大于相应原泥的参数,说明淤泥经过固化之后,黏性增强,动切力增大。并且固化淤泥的动切力和黏滞系数随初始含水率增大也在减小,但减小的速率要明显大于未处理的原泥,说明初始含水率对新拌固化淤泥的黏滞特性有非常显著的影响。
(2)水泥添加量对新拌固化淤泥黏滞性影响
为了研究不同水泥添加量对固化淤泥黏滞特性的影响,对白马湖淤泥选择296%,190%,118%三个不同初始含水率,对梅梁湖淤泥选择130%的初始含水率,进一步研究分别添加100,200 kg/m3的水泥后的固化淤泥的黏滞性,上述含水率范围可以涵盖淤泥为宾汉姆流体和幂律流体两种情况。通过分析发现在初始含水率相同的情况下,不同水泥添加量会对淤泥的黏滞特性产生影响,但在水泥添加量从50 kg/m3提高到200 kg/m3时,固化淤泥黏滞曲线的形态并没有明显变化,仅列出初始含水率为296%,水泥添加量为100 kg/m3时的新拌白马湖固化淤泥的黏滞特性试验结果,可见呈现出典型宾汉姆流体的特征,见图 12。
虽然不同的水泥添加量下固化淤泥的黏滞曲线形态相同,但对黏滞参数的影响需要进一步明确,以新拌固化淤泥属于宾汉姆流体时来考察不同水泥添加量对宾汉姆流体参数的影响。图 13给出了初始含水率为296%的白马湖淤泥分别添加不同水泥量下新拌固化淤泥的黏滞系数和动切力的变化曲线。由图 13可见,水泥添加量对新拌固化淤泥黏滞系数和动切力的影响主要发生在水泥量小于50 kg/m3的范围内,水泥添加量继续增加,只能使固化淤泥的黏滞系数和动切力有少量的增大。因此,水泥添加量对淤泥黏滞性的影响,主要发生在水泥量小于50 kg/m3的范围。前面的研究已经表明,水泥对新拌固化淤泥流动度的影响也主要发生在水泥量50 kg/m3的范围,原因也就在于水泥对淤泥黏滞性的影响主要发生在这个范围内。
4. 基于滨汉姆流体时的流动性和黏滞性的关系
从淤泥和新拌固化淤泥的流动度试验结果知道,对于淤泥当含水率w>2.0wL时,淤泥的流动性会随着含水率的增加而快速增大,并且坍落稳定后的形态呈现圆饼状;对于新拌固化淤泥当淤泥的初始含水率w>2.8wL时,固化淤泥的流动特性表现出和上述淤泥相类似的特征。从黏滞性试验结果知道,淤泥和固化淤泥之所以会表现出比较大的流动性,原因在于此种状态下淤泥或者新拌固化淤泥属于宾汉姆流体。
表征宾汉姆流体特征的两个参数是黏滞系数μ和动切力τ0,黏滞系数可以认为是流体开始流动后流动快慢的表征,而动切力的大小则决定了流体在外力的作用下能否进行流动。本文中的流动度试验是淤泥或者固化淤泥完全在自身重力作用下的流动,如果淤泥(固化淤泥)的动切力小于自身重力所产生的应力时,那么提起坍落度筒(截锥圆模)后,淤泥在自身重力作用下会向四周流动。因此,黏滞参数是反映流体流动特性的内在原因,而流动度只是这种流动特性的宏观表现。对于流动性比较大的属于宾汉姆流体时的淤泥或者新拌固化淤泥,有理由相信黏滞参数和流动度之间有内在联系,换句话讲,如果知道了宾汉姆流体参数也就应该可以知道流动度的结果。
当淤泥或新拌固化淤泥的流动性很大时,稳定状态后的水平流动度要远大于竖直方向的坍落度。淤泥在坍落度筒中时有一个自由表面,对于这种有自由表面的流体的流动,可采用浅水长波方程求解。坍落度筒中的淤泥在自重作用下的流动,可以看作恒定体积的黏性液体在水平面上的释放,初始淤泥的形状是截锥状,可以看作轴对称问题来求解,采用如图 14所示柱坐标系统(O,r,θ,z),并作以下假设:①流体是不可压缩的,密度是恒定不变的;②水平方向的扩展距离远大于竖直方向的坍落距离;③坍落流动过程中水固分离过程忽略不计;④在自重作用下流动稳定后,在水平面上的投影是规则的圆形。
按照浅水长波方程,可得平衡方程为
∂p∂r=∂τ∂z, (5) ∂p∂z=−ρg。 (6) 淤泥或者新拌固化淤泥在流动时,其黏滞特性满足Bingham方程,即满足
τ=τ0+μ˙γ, (7) 式中,ρ为淤泥的密度,g为重力加速度。τ为淤泥的剪切应力,τ0为淤泥的动切力,˙γ为剪切速率。则有
˙γ=∂Vr∂z。 (8) 假设淤泥压力p的分布符合流体静力学,则可得
p=ρg(h(r)−z), (9) 式中,h(r)为淤泥的深度,z=h(r)处的压力为大气压力,假设气液交界面的表面张力相对于动切力可以忽略不计。联立式(5),(6),并考虑边界条件z=0和z=h得
ρgh(r)dh(r)dr=−τ0−f(˙γ)|z=0。 (10) 在停止的时刻,˙γ=0,式(10)可以简化为
ρgh(r)dh(r)dr=τ0。 (11) 进一步考虑边界条件h(r=R)=0,可得坍落稳定后的淤泥高度为
h(r)=[2τ0(R−r)ρg]12。 (12) 淤泥的体积V可以通过积分得到:
V=∫2π 0∫R0h(r)rdθ=8√2τ0π R5215√ρg。 (13) 则可得坍落稳定后半径R的表达式为
R=(225ρgV2128π 2τ0)15。 (14) 流动度D为流动稳定后的直径,所以
D=2R=2(225ρgV2128π 2τ0)15。 (15) 由式(15)可知,流动度可以根据试样体积、屈服应力、重度来进行预测。同时也可以看出,当淤泥或者新拌固化淤泥满足宾汉姆流体方程时,屈服应力和流动度之间有唯一对应的关系,和淤泥的种类无关,并且对于新拌固化淤泥也有这样的对应关系。对于本文中所有的含水率在2wL以上的淤泥和初始含水率在2.8wL以上的新拌固化淤泥,如果把黏滞性试验中得到的动切力τ0作为已知条件,则可以根据式(15)对淤泥或者新拌固化淤泥的流动度进行预测,其中截锥的体积V=344 cm3。按照此种方法得到的流动度和动切力的关系见图 15。可见根据动切力预测的流动度和实测的流动度之间吻合得很好,说明此种方法可以有效地把两者之间联系起来。
5. 结论
采用白马湖和梅梁湖的淤泥,选择水泥作为固化材料,通过截锥圆模流动度试验和旋转黏度计试验,研究了淤泥和新拌固化淤泥的流动性和黏滞特性,得到以下6点结论。
(1)初始含水率对淤泥的流动性影响显著,并存在临界含水率w′≈2.0wL,当淤泥初始含水率高于2.0wL时,淤泥的流动性迅速增大。
(2)淤泥固化后,流动性降低,初始含水率相同时,水泥添加量越高,新拌固化淤泥流动性越小,且主要发生在水泥添加量ac≤50 kg/m3。初始含水率对新拌固化淤泥流动性的影响也存在临界含水率w′s,w′s和水泥添加量有关,当ac=50 kg/m3,w′s≈2.8wL,此时新拌固化淤泥的流动性约为未固化淤泥流动性的2/3。
(3)黏滞性试验表明,含水率w>12.0wL,淤泥为牛顿流体;2.0wL≤w≤12wL,为宾汉姆流体;w<2.0wL,为带屈服值的伪塑性流体。当淤泥的初始含水率w≥2.8wL时,新拌固化淤泥为宾汉姆流体;当w<2.8wL时,新拌固化淤泥近似为带屈服值的幂律流体。
(4)新拌固化淤泥的黏滞系数μ和动切力τ0都随初始含水率增大而减小,并且都大于对应原泥的黏滞参数,水泥添加量对对黏滞参数数的影响主要发生在水泥添加量ac≤50 kg/m3。
(5)对于淤泥,当含水率w≥2.0wL时,或者对于新拌固化淤泥,当初始淤泥含水率w≥2.8wL时,即当淤泥或者新拌固化淤泥满足宾汉姆流体方程时,根据浅水方程推导得到了屈服应力和流动度之间的唯一关系,即流动度可以由试样体积、重度和屈服应力获得,也说明黏滞特性参数是决定淤泥宏观流动性的内在原因。
(6)流动度和黏滞性试验联合表明,新拌固化淤泥流动性小于淤泥,根本原因在于固化材料的加入使得水泥和水发生化学反应从而成为振凝性流体,随着时间的延长,固化淤泥逐渐从流体变成为固体状态,从而实现了淤泥的固化作用。本文的规律性结论是基于两种淤泥,水泥作为单一固化材料通过大量试验得出的,后续需要对更多的淤泥进行研究,以检验结论是否有普适性,但无疑,本文对于大规模淤泥流动固化浇筑施工具有重要参考价值。
-
表 1 试验土料强度特性
Table 1 Strength properties of test clay
KIC/(kPa·m0.5) σt /kPa c/kPa φ/(°) 17 62 68 22 表 2 试验方案
Table 2 Test schemes
类型 裂缝倾角β/(°) 无量纲裂缝长度2a/W 闭合 0,15,30,45,60,75 0.5 闭合 45 0.2,0.35,0.65 非闭合 15,30,45,60 0.5 注: W为试样的边长,100 mm。 -
[1] 张丙印, 温彦锋, 朱本珍. 土工构筑物和边坡工程发展综述——作用机理与数值模拟方法[J]. 土木工程学报, 2016, 49(8): 1-15. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201608001.htm ZHANG Bing-yin, WEN Yan-feng, ZHU Ben-zhen, et al. Development review of geotechnical structures and slope engineering-action mechanism and numerical simulation method[J]. Chinese Journal of Civil Engineering, 2016, 49(8): 1-15. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201608001.htm
[2] 吉恩跃, 陈生水, 傅中志. 土心墙堆石坝坝顶裂缝扩展有限元模拟[J]. 岩土工程学报, 2018, 40(增刊2): 17-21. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2018S2006.htm JI En-yue, CHEN Sheng-shui, FU Zhong-zhi. Numerical simulation of crest cracks in an earth core rockfill dam using extended finite element method[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(S2): 17-21. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2018S2006.htm
[3] 凌华, 王伟, 王芳. 砾石土心墙料水力劈裂试验研究[J]. 岩土工程学报, 2018, 40(8): 1444-1448. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201808011.htm LING Hua, WANG Wei, WANG Fang. Experimental study on hydraulic fracture of gravelly soil core[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(8): 1444-1448. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201808011.htm
[4] CAMP S, GOURC J P, PLE O. Landfill clay barrier subjected to cracking: multi-scale analysis of bending tests[J]. Applied Clay Science, 2010, 48(3): 384-392. doi: 10.1016/j.clay.2010.01.011
[5] SKEMPTON A W, SCHUSTER R L, PETLEY D J, et al. Joints and fissures in the london clay at wraysbury and edgware[J]. Géotechnique, 1970, 20: 208-209. doi: 10.1680/geot.1970.20.2.208
[6] THUSYANTHAN N I, TAKE W A, MADABHUSHI S P G, et al. Crack initiation in clay observed in beam bending[J]. Géotechnique, 2007, 57(7): 581-594. doi: 10.1680/geot.2007.57.7.581
[7] WANG J J, ZHU J G, CHIU C F, et al. Experimental study on fracture toughness and tensile strength of a clay[J]. Engineering Geology, 2007, 94(1/2): 65-75.
[8] RAO Q, SUN Z, STEPHANSSON O, et al. Shear fracture (Mode II) of brittle rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2003, 40(3): 355-375. doi: 10.1016/S1365-1609(03)00003-0
[9] JI W, PAN P, LIN Q, et al. Do disk-type specimens generate a mode II fracture without confinement?[J]. International Journal of Rock Mechanics and Mining Sciences, 2016, 87: 48-54. doi: 10.1016/j.ijrmms.2016.05.010
[10] LIN H, YANG H T, WANG Y X, et al. Determination of the stress field and crack initiation angle of an open flaw tip under uniaxial compression[J]. Theoretical and Applied Fracture Mechanics, 2019, 104: 102358. doi: 10.1016/j.tafmec.2019.102358
[11] 郑安兴, 罗先启. 压剪应力状态下岩石复合型断裂判据的研究[J]. 岩土力学, 2015(7): 1892-1898. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201507014.htm ZHENG An-xing, LUO Xian-qi. Research on combined fracture criterion of rock under compression-shear stress[J]. Rock and Soil Mechanics, 2015(7): 1892-1898. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201507014.htm
[12] 高赛红, 曹平, 汪胜莲. 水压力作用下岩石中Ⅰ和Ⅱ型裂纹断裂准则[J]. 中南大学学报(自然科学版), 2012, 43(3): 1087-1091. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201203046.htm GAO Sai-hong, CAO Ping, WANG Sheng-lian. Analysis of fracture criterion of Ⅰ and Ⅱ type of crack in rock masses under hydraulic pressure[J]. Journal of Central South University (Science and Technology), 2012, 43(3): 1087-1091. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201203046.htm
[13] ZHAO Y, WANG Y, TANG L. The compressive-shear fracture strength of rock containing water based on Druker-Prager failure criterion[J]. Arabian Journal of Geosciences, 2019, 12(15): 452. doi: 10.1007/s12517-019-4628-1
[14] 周家文, 徐卫亚, 石崇. 基于破坏准则的岩石压剪断裂判据研究[J]. 岩石力学与工程学报, 2007(6): 1194-1201. doi: 10.3321/j.issn:1000-6915.2007.06.014 ZHOU Jia-wen, XU Wei-ya, SHI Chong. Investigation on compression-shear fracture criterion of rock based on failure criteria[J]. Chinese Journal of Geotechnical Engineering, 2007(6): 1194-1201. (in Chinese) doi: 10.3321/j.issn:1000-6915.2007.06.014
[15] 李世愚, 和泰名, 尹祥础. 岩石断裂力学导论[M]. 合肥: 中国科学技术大学出版社, 2010. LI Shi-yu, HE Tai-ming, YIN Xiang-chu. Introduction of Rock Fracture Mechanics[M]. Hefei: University of Science and Technology of China Press, 2010.
[16] ZHANG X, WONG L N Y. Cracking processes in rock-like material containing a single flaw under uniaxial compression: a numerical study based on parallel bonded-particle model approach[J]. Rock Mechanics and Rock Engineering, 2012, 45(5): 711-737.
[17] SMITH D J, AYATOLLAHI M R, PAVIER M J. The role of T-stress in brittle fracture for linear elastic materials under mixed-mode loading[J]. Fatigue Fracture of Engineering Materials and Structures, 2001, 24(2): 137-150. doi: 10.1046/j.1460-2695.2001.00377.x
[18] LIU H. Wing-crack initiation angle: a new maximum tangential stress criterion by considering T-stress[J]. Engineering Fracture Mechanics, 2018, 199: 380-391. doi: 10.1016/j.engfracmech.2018.06.010
[19] 朱俊高, 王俊杰, 张辉. 土石坝心墙水力劈裂机制研究[J]. 岩土力学, 2007(3): 487-492. doi: 10.3969/j.issn.1000-7598.2007.03.011 ZHU Jun-gao, WANG Jun-jie, ZHANG Hui. Study on mechanism of hydraulic fracturing in core of earth-rockfill dam[J]. Rock and Soil Mechanics, 2007(3): 487-492. (in Chinese) doi: 10.3969/j.issn.1000-7598.2007.03.011
[20] 王俊杰. 基于断裂力学的土石坝心墙水力劈裂研究[D]. 南京: 河海大学, 2005. WANG Jun-jie. Study on Hydraulic Fracturing in Core of Earth-Rock Fill Dam Based on Fracture Mechanics[D]. Nanjing: Hohai University, 2005. (in Chinese)
[21] WILLIAMS J G, EWING P D. Fracture under complex stress —the angled crack problem[J]. International Journal of Fracture Mechanics, 1972, 8(4): 441-446. doi: 10.1007/BF00191106
[22] ALIHA M R M, AYATOLLAHI M R. Two-parameter fracture analysis of SCB rock specimen under mixed mode loading[J]. Engineering Fracture Mechanics, 2013, 103: 115-123. doi: 10.1016/j.engfracmech.2012.09.021
[23] TANG S B, BAO C Y, LIU H Y. Brittle fracture of rock under combined tensile and compressive loading conditions[J]. Canadian Geotechnical Journal, 2017, 54(1): 88-101. doi: 10.1139/cgj-2016-0214
[24] SCHMIDT R. A microcrack model and its significance to hydraulic fracturing and fracture toughness testing[J]. US Symposium on Rock Mechanics, 1980, 18(5): 581-590.
[25] AYATOLLAHI M R, ALIHA M R M. On the use of Brazilian disc specimen for calculating mixed mode I-II fracture toughness of rock materials[J]. Engineering Fracture Mechanics, 2008, 75(16): 4631-4641. doi: 10.1016/j.engfracmech.2008.06.018
[26] 赵彦琳, 范勇, 朱哲明, 等. T应力对闭合裂纹断裂行为的理论和实验研究[J]. 岩石力学与工程学报, 2018, 37(6): 1340-1349. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201806003.htm ZHAO Yan-lin, FAN Yong, ZHU Ze-ming, et al. Analytical and experimental study on the effect of T-stress on behavior of closed cracks[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(6): 1340-1349. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201806003.htm
[27] WANG J J, HUANG S Y, GUO W L, et al. Experimental study on fracture toughness of a compacted clay using semi-circular bend specimen[J]. Engineering Fracture Mechanics, 2020, 224: 106814. doi: 10.1016/j.engfracmech.2019.106814
[28] 蒲成志, 杨仕教, 张春阳. 张开度影响的裂隙体破断机制探讨[J]. 岩土工程学报, 2019(10): 1836-1844. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201910009.htm PU Cheng-zhi, YANG Shi-jiao, ZHANG Chun-yang. Fracture mechanism of fracture body affected by opening degree[J]. Chinese Journal of Geotechnical Engineering, 2019(10): 1836-1844. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201910009.htm
[29] 郭少华. 岩石类材料压缩断裂的实验与理论研究[D]. 长沙: 中南大学, 2003. GUO Shao-hua. The Experimental and Theoretical Investigation on Fracture of Rock-Type Materials under Compressive Loading[D]. Changsha: Central South University, 2003. (in Chinese)
[30] MUSKHELISHVILI N. Some Basic Problems of the Mathematical Theory of Elasticity[M]. Leyden: Noordhoff, 1953.
[31] HUA W, DONG S, FAN Y, et al. Investigation on the correlation of mode II fracture toughness of sandstone with tensile strength[J]. Engineering Fracture Mechanics, 2017, 184: 249-258. doi: 10.1016/j.engfracmech.2017.09.009
[32] BOBET A. The initiation of secondary cracks in compression[J]. Engineering Fracture Mechanics, 2000, 66(2): 187-219. doi: 10.1016/S0013-7944(00)00009-6
-
期刊类型引用(1)
1. 谭明伦,仝令帅,周鸣亮,张乐,黄宏伟. 数物双驱动的水下隧道围岩稳定性分析. 现代隧道技术. 2024(S1): 183-193 . 百度学术
其他类型引用(1)