考虑3种因素影响的硫氧镁水泥固化土修正邓肯-张模型

朱剑锋，徐日庆，罗战友，饶春义

(1.浙江科技学院土木与建筑工程学院，浙江杭州，310023；2.浙江科技学院隧道与地下空间技术开发研究院，浙江杭州，310023；3.浙江大学滨海和城市岩土工程研究中心，浙江杭州，310058；4.浙江加州国际纳米技术研究院台州分院，浙江台州，318000；5.宁波大学土木与环境工程学院，浙江宁波，315211)

固化土的无侧限抗压强度主要受固化剂掺量、初始含水量、龄期等因素影响，人们对此开展了很多研究，例如：朱伟等[1]提出了水泥固化土无侧限抗压强度与水泥添加量之间的经验公式，但未考虑龄期的影响；董邑宁等[2]研究了掺入比、龄期对固化土强度的影响规律；KOLIAS 等[3]分析了不同配比下固化土的无侧限抗压强度随龄期的演化规律。软土初始含水量也是影响软土加固效果的重要因素，HORPIBULSUK等[4]研究了固化淤泥质土无侧限抗压强度随含水量、水泥掺量和养护龄期的变化规律；徐日庆等[5]建立了考虑初始含水量、固化剂掺量等因素影响的固化土综合抗压强度预测模型；杨爱武等[6]基于无侧限抗压强度试验结果，分析了含水率、龄期、固化剂掺量对吹填泥浆固化土强度的影响规律。除了无侧限抗压度指标外，固化土的抗剪强度、变形特征以及它们之间的关系(本构关系)也逐渐引起了学者们的广泛关注。童小东等[7]结合水泥固化土的试验结果，提出了水泥土的损伤理论模型；王军等[8]根据不排水三轴试验建立了考虑水泥土刚度软化的固化土模型；张涛等[9]建立了考虑胶结作用的木质素固化粉土的边界面模型。邓肯-张模型参数少、物理意义明确，且能够反映土体应力-应变关系的非线性，因此，胡亚元等[10]提出了考虑纤维掺入量的水泥土邓肯-张模型；孙凯等[11]基于水泥土的结构性，建立了符合水泥土力学特性的弹塑性本构模型；然而，在传统硅酸盐水泥型固化剂的生产过程中会产生大量的废气和灰尘，对环境造成严重的污染[12]。于是，朱剑锋等[13]研制了节能环保型镁质水泥复合固化剂(TZ18)来替代传统硅酸盐水泥并取得了良好的固化效果；在此基础上，饶春义等[14]提出了镁质水泥固化土的一维压缩模型。但是，上述固化土本构模型大多仅考虑了单一因素(如水泥掺量等)对固化土力学性质的影响，而无法预测多重因素影响下(如初始含水量、龄期、固化剂掺量等)固化土的力学特性。鉴于初始含水量(w)、固化剂掺量(Wg)和龄期(T)对水泥固化土的力学性能的显著影响[1-6,15-16]，本文作者基于岩土工程中广泛应用的邓肯-张模型[10,17-18]，通过开展不同初始含水量、固化剂掺量和龄期下硫氧镁水泥固化土不排水三轴试验，研究各因素对邓肯-张模型参数的影响规律，并综合考虑这3种因素的影响建立硫氧镁水泥固化土修正邓肯-张模型。

1 试验材料与方案

1.1 试验材料

选用宁波②2-2层淤泥质黏土为试验用土，其物理力学指标见表1，其中，w 为土样的天然含水量，γ为重度，e为孔隙比，wp为塑限，wL为液限，Es1-2为压缩模量，c为黏聚力，φ为内摩擦角。将淤泥质土烘干、碾碎、过孔径2 mm筛。试验采用镁质水泥复合固化剂TZ18[13]。

表1 土样的物理力学指标Table 1 Physical and mechanical properties of soil

1.2 试验方案

考虑初始含水量w、固化剂掺量Wg和龄期T对硫氧镁水泥固化土力学性能的影响，设置试验方案，如表2 所示，M0为基准配比试样，每种配比试样制备3个平行试样。

表2 硫氧镁水泥固化土三轴试验方案Table 2 Triaxial test plan for soft soil solidified by magnesium oxysulfate cement

首先，将淤泥质黏土烘干，然后按照设计配比称取过筛后的干土、水、改性硫氧镁水泥和外加剂(水玻璃、熟料、硅灰)充分搅拌成均匀的复合固化剂(TZ18)混合浆液，将其缓缓注入到试验淤泥中，并充分混合搅拌均匀。最后分3层装入三瓣膜(直径×高度为39.1 mm×80.0 mm)中，每层振捣2～3 min 排出气泡。试样在自然条件下养护2 d 后拆模，然后移至恒温通风环境下继续养护至设计龄期。由于水化反应作用，硫氧镁水泥固化土试样养护过程会有一定的收缩，因此，在三轴试验前需要将试样上下平面打磨光滑，并测量其直径和高度，后续计算以实测体积为依据。

1.3 试验仪器与方法

镁质水泥固化土的不排水三轴试验在GDS 高压与非饱和土动三轴仪器上进行，将试样装在加载底座上，然后加上压力室，采用GDSLAB 软件设置试验参数和数据收集，其中，轴力传感器与顶帽的接触轴力为5.6 N，试验的围压σ3分别为100，200，300 kPa，围压的加载速率为0.02 kPa/s，剪切速率为0.05%/min，当轴向应变达到20%或轴力达到2 MPa 时试验自动停止[19]。

2 试验结果

2.1 硫氧镁水泥固化土的破坏形态

硫氧镁水泥固化土侧向鼓胀破坏形态如图1所示。可见，硫氧镁水泥固化土的破坏形式主要有：1)侧向鼓胀破坏，试验后固化土表面没有明显的裂纹，大约在中间部位对称鼓起。这是由于TZ18固化剂水化反应生成的针状晶体[20]和胶体(由TZ18固化剂中的外加剂熟料、硅灰等反应生成)较少，且没有硬化胶结在一起形成一个结构体，固化土的强度低，具有较大的塑性，在较大的应变下发生破坏；2)斜向剪切破坏，试验后固化土表面有明显的破裂面，破裂面与水平方向大致成45°～60°(见图2)。这是因为水化反应生成的晶体和胶体把分散的土颗粒黏结在一起，贯穿整个试样的内部，形成一定的结构，从而能抵抗较大的荷载作用。

2.2 硫氧镁水泥固化土的应力-应变曲线

2.2.1 基准配比下的固化土三轴试验结果

图1 硫氧镁水泥固化土侧向鼓胀破坏形态Fig.1 Lateral bulging failure morphology of soft soil solidified by magnesium oxysulfate cement

图2 硫氧镁水泥固化土斜向剪切破坏形态Fig.2 Oblique shear failure morphology of soft soil solidified by magnesium oxysulfate cement

图3所示为基准试样M0的应力-应变曲线，其中，σ1-σ3为主应力差；ε1为轴向应变。由图3 可知：在基准配比下，硫氧镁水泥固化土应力-应变曲线呈双曲线形式，且与淤泥质土三轴试验结果类似[19]：围压越大，曲线的初始切线斜率越大并且渐进线越高。

图3 基准试样M0的应力-应变曲线Fig.3 Stress-strain curves of solidified soils sample M0

2.2.2 不同初始含水量的固化土三轴试验结果

图4 不同初始含水量下的硫氧镁水泥固化土的应力-应变曲线Fig.4 Stress-strain curves of soft soils solidified by magnesium oxysulfate cement with different initial water contents

图4所示为不同初始含水量淤泥固化后的应力-应变曲线关系，其中，Wg=15%，T=7 d。由图4可知：当初始含水量为40%～45%，围压一定时，随着固化土轴向应变ε1的增大，主应力差(σ1-σ3)先增大到峰值，然后下降，直到保持在稳定值附近，固化土的应力-应变曲线属于加工软化型；随着含水量的增加，硫氧镁水泥固化土的强度和刚度(初始切线模量)逐渐降低，应力-应变曲线逐渐由加工软化型转化为加工硬化型。这是因为w的增加提高了软土的孔隙比，降低了水化反应的离子浓度，抑止了硫氧镁水泥的水化反应，从而降低了硫氧镁水泥固化土的强度和刚度。另外，当w一定时，随着围压的增加，固化土逐渐被压密，从而提高了其刚度及强度。

2.2.3 不同固化剂掺量下的固化土三轴试验结果

图5所示为不同固化剂掺量的固化土在不同围压下的应力-应变曲线关系，其中，w=50%，T=7 d。可见：Wg介于5%～15%且围压一定时，固化土的主应力差(σ1-σ3)随着轴向应变ε1的增大而增大，应力-应变关系呈加工硬化型。随着Wg的增加，硫氧镁水泥固化土的强度和刚度逐渐提高，应力-应变曲线逐渐由加工硬化型转化为加工软化型。这主要是因为Wg增加提高了水化反应的离子浓度，生成了更多的固化产物，减小了淤泥中的孔隙比，并形成一定空间网格结构，使得硫氧镁水泥固化土的强度和刚度较原状淤泥质土均得到了显著提高。

2.2.4 不同龄期下的硫氧镁水泥固化土三轴试验结果

图6所示为不同龄期的固化土在不同围压下的应力-应变关系，其中，w=50%，Wg=15%。从图6可以看出：龄期较短时，固化土的主应力差(σ1-σ3)较小，随着龄期的增长固化土的主应力差有较大的变化。龄期在3～7 d 时，固化土的应力-应变曲线属于加工硬化型，随着龄期的增长，固化土固化产物逐渐增多，对原状淤泥质土的微观结构产生了较好胶结和填充，从而使得固化土的强度和刚度逐渐增大，应力-应变曲线逐渐转化为加工软化型。

图5 不同固化剂掺量下的硫氧镁水泥固化土的应力-应变曲线Fig.5 Stress-strain curve of soft soils solidified by magnesium oxysulfate cement with different curing agents contents

3 硫氧镁水泥固化土修正邓肯-张本构模型

3.1 邓肯-张模型简介

由图3～6可知：随着含水量、固化剂掺量和龄期的变化，硫氧镁水泥固化土的应力-应变曲线主要表现为加工硬化型和软化型，对于硬化型曲线及软化型曲线的峰前阶段，可采用邓肯-张模型来进行模拟：式中：a和b为试验常数，且a=1/Ei，b=1/(σ1-σ3)ult，(σ1-σ3)ult为双曲线渐近线对应的极限偏差应力，Ei为应力-应变曲线的初始切线模量，与围压存在以下的函数关系：

式中：K 和n 为模型参数；pa为标准大气压，本文取pa=101.3 kPa。定义破坏比Rf为[17]

其中：(σ1-σ3)f为破坏应力，根据摩尔-库仑强度准则，有：

综上，原始邓肯-张模型共有K，n，c，φ，Rf5个参数，可通过三轴试验进行标定。

3.2 硫氧镁水泥固化土的变形参数及破坏比

硫氧镁水泥固化土的Ei可近似用下式计算[17]：

图6 不同龄期下硫氧镁水泥固化土的应力-应变曲线Fig.6 Stress-strain curves of soft soils solidified by magnesium oxysulfate cement under different ages

其中，下标95%和下标70%分别代表当(σ1-σ3)为(σ1-σ3)f的90%和70%时的试验数据；对于图3～6所示的加工硬化型应力-应变曲线，取ε1=15%时的主应力差为(σ1-σ3)f[19]，对加工软化型，则取峰值处的主应力差为(σ1-σ3)f；(σ1-σ3)ult可近似通过下式计算：

根据式(5)和(6)可得试样M0～M12的初始切线模量Ei(见图7)。对式(2)两边均除以pa，并取对数可得：

图7 不同围压下硫氧镁水泥固化土的初始切线模量EiFig.7 Ei of soft soils solidified by magnesium oxysulfate cement under different confining pressures

将初始切线模量Ei和σ3代入式(7)可得试样M0～M12的参数K 和n，见图8(a)和(b)。从图8(a)和(b)可以看出：模型参数n在0.803～0.828之间变化，可近似认为n不受含水量、固化剂掺量以及龄期的影响，取其平均值(n=0.812)作为硫氧镁水泥固化土本构模型参数n的经验值。

根据式(3)和(6)并结合试验结果可得各试样的破坏比Rf，如图8(c)所示。由图8(c)可知：硫氧镁水泥固化土的Rf波动范围较小，取其平均值(Rf=0.822)作为硫氧镁水泥固化土修正邓肯-张模型参数Rf经验值。

3.3 硫氧镁水泥固化土的抗剪强度参数

试样M0～M12的黏聚力c 和内摩擦角φ 见图9。由图8和图9可知：原始邓肯-张模型参数n和Rf受w，Wg和T 影响较小，可近似取n=0.812 和Rf=0.822 作为硫氧镁水泥固化土的模型参数经验值。而K，c 和φ 对它们比较敏感。因此，有必要对原始邓肯-张模型进行修正，使其关键参数K，c和φ能体现w，Wg和T的影响，从而实现不同条件下的硫氧镁水泥固化土的力学特性精确预测。

3.4 硫氧镁水泥固化土的修正邓肯-张模型

硫氧镁水泥固化土中的水来源于淤泥的初始含水量和固化剂中的水，可以通过灰水比(C/W)来综合考虑初始含水量(w)和固化剂掺量(Wg)的固化土模型参数间的关系[21-23]：

式中：ws/Cs为固化剂的水灰比，且ws/Cs=0.155[13]。通过计算得出基准试验方案的硫氧镁水泥固化土的灰水比C0/W0=0.421，其中，C0和W0分别为基准配比下固化剂中固态材料和水的质量分数，C0为12.68%，W0为52.33%。

3.4.1 参数K的经验公式

根据图8(a)和式(8)可得归一化处理后参数K与灰水比的关系曲线，如图10(a)所示。由图10(a)可知：当灰水比较小时，K变化不大，随着灰水比增大，K 缓慢增大；当灰水比较大时，K 增长非常快，经拟合K与灰水比之间近似呈下述函数关系：

式中：K0为基准配比的镁质中水泥固化土的原始邓肯-张模型参数K；d1和f1为综合考虑灰水比对参数K影响的参数，本文取d1=1.018，f1=4.360。

由图8可知：固化土的参数K不仅与灰水比有关，与养护龄期也有很大的关系。以基准配比的硫氧镁水泥固化土的模型参数K0和灰水比(C/W)0为基础，对归一化后的K与龄期进行拟合，结果见图10(b)。由图10(b)可知：随着龄期增长，K 呈幂函数增长，通过拟合得到参数K与T的经验关系为式中：T0为龄期归一化参考值，本文取T0=7 d；d2和f2为龄期对K 的影响参数，本文取d2=0.723，f2=4.071。

图8 硫氧镁水泥固化土的模型参数Fig.8 Model parameters of soft soils solidified by magnesium oxysulfate cement

图9 硫氧镁水泥固化土的强度参数Fig.9 Shear strength parameters of soft soils solidified by magnesium oxysulfate cement

假设龄期与硫氧镁水泥固化土参数K存在如式(10)所示的对数关系，即可建立综合考虑w，Wg和T影响的硫氧镁水泥固化土参数K的经验公式：

3.4.2 抗剪强度指标c和φ的经验公式

根据图9和式(8)可分别得到归一化处理后的参数c和φ与水灰比之间的关系曲线，见图11。由图11 可知：与K-C/W 关系类似，随着灰水比增加，固化土的c 呈幂函数形式增大，φ 呈线性增大。对图11所示结果进行拟合，可以分别得到c和φ与灰水比之间的经验公式：

式中：c0和φ0分别为基准配比的硫氧镁水泥固化土试样的黏聚力和内摩擦角；d3和f3为综合考虑灰水比对参数c影响的参数；d5为综合考虑灰水比对参数φ 影响的参数。本文取d3=1.024，f3=2.827，d5=1.082。

以基准配比的硫氧镁水泥固化土试样的黏聚力c0、内摩擦角φ0以及灰水比(C/W)0为基础，分别对归一化后的c 和φ 与龄期进行拟合，结果见图12。与图11类似，随着龄期增加，固化土的c和φ分别呈幂函数和线性形式增大，对图12 拟合，可得：

图10 K与灰水比和龄期的关系Fig.10 Relationship between K and cement-water ratio and age

图11 抗剪强度指标与灰水比的关系Fig.11 Relationship between shear strength parameters and cement-water ratio

图12 抗剪强度指标与龄期的关系Fig.12 Relationship between shear strength parameters and age

式中：d4和f4为龄期对c的影响参数；d6为龄期对φ的影响参数。本文取d4=0.937，f4=1.712，d6=0.961。

分别以式(12)和式(13)中的c 和φ 替换式(14)和式(15)中的c0和φ0，可建立综合考虑初始含水量、固化剂掺量和龄期影响的硫氧镁水泥固化土抗剪强度指标c和φ的经验公式：

3.4.3 修正邓肯-张模型

硫氧镁水泥固化土的力学性质受初始含水量w、固化剂掺量Wg和龄期T影响较大，其本构模型必须考虑上述因素影响才能准确描述固化土力学特性。根据以上试验结果，结合式(1)，(2)，(4)，(11)，(16)和(17)可建立如式(18)所示综合考虑w，Wg和T 影响的硫氧镁水泥固化土修正邓肯-张模型。该模型包含5 个模型参数c0，φ0，K0，Rf，n，3个基准参数c0，W0，T0以及10个拟合参数d1，d2，d3，d4，d5，d6，f1，f2，f3，f4，且所有参数通过三轴试验即可确定。

4 模型验证

为验证本文硫氧镁水泥固化土修正邓肯-张模型的合理性，选取宁波②2-2层淤泥质黏土制样，控制淤泥的初始含水量分别为40%，42%和50%，固化剂掺量分别为7%，10%和20%，养护龄期分别为5，7 和9 d，进行常规三轴不排水剪切试验，C0=12.68%，W0=52.33%，T0=7 d，其他模型参数见表3，实测数据和模型预算结果见图13。由图13(c)，(d)，(e)，(g)可知：对于加工硬化型曲线，本文修正邓肯-张模型的预测结果始终与试验结果一致。而对于加工软化型曲线，仅对峰值前的试验结果进行预测，结果见图13(a)，(b)，(f)，(h)。由图13(a)，(b)，(f)，(h)可知：修正邓肯-张模型可较好地预测峰值前的加工软化型应力-应变曲线。因此，本文硫氧镁水泥固化土修正邓肯-张模型较好地实现了任意初始含水量、固化剂掺量和龄期下固化土应力-应变关系的准确预测。

表3 硫氧镁水泥固化土模型参数Table 3 Model parameters of soft soils solidified by magnesium oxysulfate cement

图13 硫氧镁水泥固化土的应力-应变实测值与预测值Fig.13 Measured and predicted stress-strain curves of soft soil solidified by magnesium oxysulfate cement

5 结论

1)初始含水量、固化剂掺量和龄期对固化剂水化反应速度、水化产物的数量和凝胶材料硬化程度均产生不同程度影响，从而使得硫氧镁水泥固化土的强度和结构性有较大区别。硫氧镁水泥固化土的破坏形式主要分为侧向鼓胀破坏和斜向剪切破坏，硫氧镁水泥固化土的应力-应变关系分别呈加工硬化和软化2种类型。

2)初始含水量、固化剂掺量和龄期对邓肯-张模型参数c，φ，K影响显著，而对n和Rf基本没有影响。

3)通过引入灰水比的概念，并考虑龄期的影响，建立了硫氧镁水泥固化土的参数c，φ，K 与w，Wg和T 的经验公式，并提出了考虑3 种因素影响的硫氧镁水泥固化土的修正邓肯-张模型。算例验证表明，该模型能准确预测加工硬化型和峰值前的加工软化型硫氧镁水泥固化土的应力-应变曲线。

4)本文提出的修正邓肯-张模型无法预测硫氧镁水泥固化土的软化特性、剪胀性和塑性变形，在下一步研究中，将采用弹塑性或亚塑性模型来描述硫氧镁水泥固化土的力学特性。