不同膨胀潜势等级的膨胀土残余强度环剪试验研究

周葆春; 王江伟; 单丽霞; 李颖; 郎梦婷; 孔令伟

doi:10.11779/CJGE20230225

不同膨胀潜势等级的膨胀土残余强度环剪试验研究 English Version

1.
信阳师范大学河南省非饱和土与特殊土工程技术研究中心/建筑与土木工程学院, 河南信阳 464000
2.
中国科学院武汉岩土力学研究所, 湖北武汉 430071

基金项目:

国家自然科学基金面上项目 11772290

详细信息

作者简介:
周葆春(1978—)，男，教授，主要从事非饱和土与特殊土力学及工程方面的研究工作。E-mail: zhoubc@xynu.edu.cn

中图分类号: TU411
计量
- 文章访问数: 292
- HTML全文浏览量: 38
- PDF下载量: 92
出版历程
- 收稿日期: 2023-03-14
- 网络出版日期: 2024-06-04
- 刊出日期: 2024-05-31

Torsional ring shear tests on residual strength of expansive soils with different swelling potentials

1.
Henan Unsaturated Soil and Special Soil Engineering Technology Research Center, College of Architecture and Civil Engineering, Xinyang Normal University, Xinyang 464000, China
2.
Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China

摘要

摘要: 为探讨不同膨胀潜势等级的膨胀土残余强度与其成分特征的相关性，采用经自主改装的环剪仪开展Denver强膨胀土、荆门中、弱膨胀土环剪试验研究，并与排水反复直接剪切试验结果进行系统对比分析，在此基础上，用实测结果探讨残余强度与土成分特征相关关系的适用性。结果表明：①环剪试验测得Denver强膨胀土、荆门中、弱膨胀土有效残余强度内摩擦角分别为11.3°，12.3°，23.7°，比排水反复直接剪切试验测得的低1.7°~3.3°，并无应力相关性且具很好精度，更适于在工程中应用。②相关经验关系和公式存在高估或者低估膨胀土有效残余强度内摩擦角的情况。③残余强度与土成分特征间定性关系非常明确，上述经验关系对一般土性状况下残余强度预测较为准确；但在土成分复杂状况下，定量准确预测较为困难，这也反映出实测的重要性与必要性。
- 膨胀土 /
- 残余强度 /
- 膨胀潜势 /
- 环剪试验 /
- 反复直接剪切试验 /
- 成分特征
Abstract: The correlation between residual strength and compositional characteristics of expansive soils with different swelling potentials is investigated. The torsional ring shear tests of Denver expansive soil with high swelling potential and Jingmen expansive soils with medium and low swelling potentials are carried out by using the modified torsional ring shear system, and the results are compared with those of the reversal direct shear tests proposed by the Chinese standard for geotechnical testing method (GB/T 50123−2019). The empirical correlations for the drained residual shear strength suggested by Mitchell and Soga (2005), Stark et al. (2013) and Lehmann et al. (2021) are assessed by the measured results. The conclusions are drawn as follows: (1) The residual friction angles of Denver expansive soil with high swelling potential and Jingmen expansive soils with medium and low swelling potentials measured by the torsional ring shear tests are 11.3°, 12.3° and 23.7°, 1.7°~3.3° lower than those measured by the reversal direct shear tests and have good accuracy, and are not stress-dependent, which means more suitable for application. (2) Mitchell and Soga (2005) and Stark et al. (2013) overestimate the residual friction angles of expansive soils, while Lehmann et al. (2021) underestimate them. (3) Although the qualitative relationship between the residual strength and the compositional characteristics of expansive soils is perfectly clear, the above empirical correlations are capable to predict the residual friction angles of expansive soils with normal compositions (Jingmen expansive soil with low swelling potentials) correctly. However, it is difficult to predict the residual friction angles of expansive soils with complicated compositions (Denver expansive soil with high swelling potential) quantitatively and accurately. In this case, the torsional ring shear tests are necessary.
- expansive soil /
- residual strength /
- swelling potential /
- torsional ring shear test /
- reversal direct shear test /
- compositional characteristic

HTML全文

0. 引言

黄土是指在第四纪地质时代中以风力作用形成的黄色粉土沉积物,是一种在特定环境中形成的具有特殊性质的土。由于黄土地区处于季节冻土区,黄土受季节性冻融作用的影响显著,每年春季发生的冻融灾害非常频繁。此外,对于盐蚀作用诱发的黄土边坡剥落等病害也不容忽视。笔者对陕北黄土高原区边坡盐蚀剥落病害进行了实地调研,每处调研点由表及里（0～30 cm）每隔5 cm取一份土样,用于室内易溶盐分析试验。试验结果表明,边坡表层主要富集Na₂SO₄等易溶盐,在冻融循环等条件下极易发生反复溶解和结晶的盐蚀作用,使被侵蚀黄土结构损伤扩展,劣化破坏,从而诱发边坡盐蚀剥落病害^[1]。因此,开展冻融环境下含盐黄土的盐蚀劣化研究具有重要意义。

目前关于冻融循环作用对盐渍土物理力学性质的影响研究已有大量文献资料。高江平^[2]系统开展了盐渍土冻胀、盐胀机理的试验与数值仿真研究工作。李振等^[3]以某干渠渠基盐渍土为研究对象,研究了冻胀量与各影响因素之间的变化规律。牛玺荣等^[4]推求了考虑冻胀和盐胀耦合作用下硫酸盐渍土体积变化模型并采用自行研制的试验系统对模型进行了验证。黄雪峰等^[5]基于室内单次降温试验,对宁夏扶贫扬黄灌溉工程沿线硫酸盐渍土盐胀特性及影响因素进行了试验研究。张莎莎等^[6]开展了多因素交互作用下砾类亚硫酸盐渍土的盐胀试验并建立了其盐胀率简化预测模型。此外,部分学者对冻融作用下黄土强度变形特性进行了研究。胡再强等^[7]探讨了冻融循环条件下杨凌重塑黄土强度变化规律并构建了其强度劣化系数拟合方程。谷琪等^[8]探究了冻融循环过程中湿陷性黄土的变形及湿陷情况。然而关于冻融作用对含盐原状黄土影响的研究较少,含盐量、冻融循环次数及其耦合效应导致土体强度衰减程度的定量化关系尚不明确。

研究学者对黄土及盐渍土微细观结构演化规律亦开展了大量研究工作。陈正汉等^[9]结合CT技术和三轴剪切设备,系统研究了黄土在多种应力路径及浸水湿陷过程中的细观结构演化特性并进一步基于CT数据构建了其细观结构演化方程。张伟等^[10]基于CT–三轴剪切试验,探究了剪切过程中硫酸盐渍土的应力应变关系和CT图像及CT数的变化规律。叶万军等^[11]研究了冻融环境下黄土体结构在不同观测尺度下的损伤演变规律。但关于冻融作用下盐渍原状黄土微结构损伤扩展演化模式的研究尚未见有专门研究报道。

盐渍原状黄土是一种具有典型天然结构强度的特殊土,冻融循环作用对其结构强度影响较大。基于此,本文选取西安Q₃原状黄土,人工制备不同Na₂SO₄含量的盐渍原状黄土试样,通过冻融循环条件下的三轴剪切及CT扫描试验,研究Na₂SO₄盐渍原状黄土冻融过程盐蚀劣化规律及微细观结构损伤演化机制。研究成果对探究黄土盐蚀作用诱发边坡剥落等病害的成灾机理具有重要的参考意义。

1. 试验材料与试样制备

1.1 试验黄土

试验黄土为取自陕西西安某基坑的Q₃原状黄土,基坑深度为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

土粒相对密度G_s	干密度ρ_d/(g·cm^-3)	孔隙比e	液限w_L/%	塑限w_P/%	塑性指数I_p
2.70	1.42	0.92	34.3	19.6	14.7

下载: 导出CSV

| 显示表格

1.2 试样制备

三轴剪切和CT扫描试样尺寸相同,均削制成直径39.1 mm、高度80 mm的标准圆柱试样。然后通过向原状土样中浸入不同浓度的Na₂SO₄盐水来人工制备含水率w为20%,Na₂SO₄含量η分别为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可见,试样不同部位含水率和含盐量变化曲线平缓,趋于一条直线,即试样内部含水率和含盐量的分布较为均匀。此外,试样不同部位的含水率和含盐量与目标值基本一致,因而人工制备盐渍原状黄土试样的方法可以满足试验要求。

图 1 试样含水率及含盐量分布

Figure 1. Water and salt profiles for standard specimens

下载: 全尺寸图片幻灯片

2. 试验方案

2.1 冻融试验

利用保鲜膜将制备好的含盐原状试样包裹,构成一个不补（散）水的密封环境,随后放入高低温试验箱进行冻融循环试验。由于三轴试样尺寸较小,土样端部和侧面换热条件虽有所差异,但影响不大。因此本次冻融循环试验为封闭系统下的多向快速冻融循环试验,以保证冻融时试样水分迁移较少。

图2所示为陕北榆林和延安地区近十年月平均最低气温变化曲线。从图中可以看出,环境温度变化区间在-20℃～20℃之间,因此以当地最低气温（约-20℃）为冻结温度,最高气温（约20℃）为融化温度,进行冻融循环试验。试验方案如下：为保证试样充分冻融,低温-20℃条件下冻结12 h,高温20℃条件下融化12 h,即一个冻融循环为24 h；冻融循环次数N分别为0,1,2,5,10。

图 2 榆林和延安地区月平均最低气温变化曲线

Figure 2. Curves of meteorological data for Yulin and Yan’an

下载: 全尺寸图片幻灯片

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所示为试样应力–应变曲线变化规律。从图中可以看出,冻融循环过程对黄土试样应力–应变曲线的类型及特征无明显影响,冻融循环前后均表现为应变硬化型。此外,相同冻融循环次数条件下,含盐量变化对应力–应变曲线的变化特征亦无显著影响,亦均表现为“硬化型”的变化趋势。

图 3 冻融作用下试样应力–应变曲线

Figure 3. Stress-strain curves under freeze-thaw action

下载: 全尺寸图片幻灯片

为深入揭示冻融条件下Na₂SO₄盐渍原状黄土强度劣化特性,取轴向应变 $ε_{1}$ 为15%时的偏应力作为破坏偏应力 ${(σ_{1} - σ_{3})}_{f}$ ,亦即试样的强度值,如图4所示。从图中可以看出,冻融循环对强度产生明显的劣化效应,破坏偏应力随着冻融循环次数增加逐渐减小,但降低速率逐渐减小,表现出显著的减速劣化特征；冻融条件下Na₂SO₄含量亦对强度产生显著的盐蚀劣化效应,破坏偏应力随着含盐量增加逐渐减小,且呈现显著的线性或加速劣化特征。

图 4 试样强度变化规律曲线

Figure 4. Variation of strength of specimens

下载: 全尺寸图片幻灯片

（2）黏聚力

图5（a）所示为黏聚力随冻融循环次数变化规律曲线。由图可见,随着冻融循环次数增加,黏聚力逐渐减小,但衰减主要集中在前5次冻融循环,之后逐渐趋于稳定,表现出减速劣化特征,这与前述三轴破坏偏应力随冻融循环次数的劣化规律是一致的（图4（a））。分析其原因,冻融循环作用下孔隙水相态变化对土颗粒产生挤压力,破坏了土颗粒之间的原生结构强度,导致黏聚力降低。多次冻融循环条件下土颗粒排列趋于平衡状态,颗粒联结强度达到稳定残余强度,黏聚力亦趋于稳定。

图 5 黏聚力变化规律曲线

Figure 5. Variation of cohesion

下载: 全尺寸图片幻灯片

黏聚力随含盐量变化规律曲线如图5（b）所示。从图中可以看出,除了未经受冻融循环（N=0）试样的黏聚力随含盐量增加无明显变化外,其它试样的黏聚力均随含盐量的增加逐渐减小,且衰减速率逐渐增大,表现出加速劣化特征,这亦与前述三轴破坏偏应力随含盐量的劣化特性表现出了很好的一致性（图4（b））。分析其原因,冻融循环条件下试样内部Na₂SO₄可溶盐的结晶—溶解—重结晶过程的反复盐蚀作用使得黄土体结构受到破坏并变得较为松散,强度显著降低。未经受冻融循环黄土试样内部的硫酸钠无相态变化亦即无法表现出盐蚀作用,因而其黏聚力无明显变化。

（3）内摩擦角

图6所示为内摩擦角随冻融循环次数及Na₂SO₄含量变化规律。由图6可见,部分试样的内摩擦角随冻融循环次数增加逐渐减小,部分试样的内摩擦角表现出波动变化特征,但整体变化幅值均较小；随含盐量增加大体上表现出先减小而后增大的特征,但变化幅值亦很小,可以认为内摩擦角随含盐量增大无显著变化。总体而言,笔者分析认为冻融作用下Na₂SO₄盐渍原状黄土试样内摩擦角变化幅值较小且无显著变化规律。这是由于内摩擦角大小主要取决于土颗粒间的摩阻力和咬合作用,而冻融循环和盐蚀作用对上述因素并无显著影响。

图 6 内摩擦角变化规律曲线

Figure 6. Variation of internal friction angle

下载: 全尺寸图片幻灯片

（4）冻融与盐蚀劣化作用解耦分析

冻融循环条件下Na₂SO₄盐渍原状黄土不仅受到孔隙水相变引起的冻融劣化作用,且要受到Na₂SO₄易溶盐相态改变引起的盐蚀劣化作用。为定量分析冻融与盐蚀劣化规律及相互关系,以下基于冻融循环条件下黏聚力变化规律,对冻融与盐蚀劣化作用进行解耦分析。

图7所示为冻融与盐蚀劣化作用解耦路径示意图。图中路径a-b表示在5次冻融循环后仅冻融作用引起的劣化值；路径b-c表示5次冻融循环后,含盐量0.5%时仅盐蚀作用引起的劣化值；路径b-d表示5次冻融循环后,含盐量1.0%时仅盐蚀作用引起的劣化值；路径b-e表示在5次冻融循环后,含盐量1.5%时仅盐蚀作用引起的劣化值。

图 7 解耦路径示意图

Figure 7. Schematic diagram of decoupling path

下载: 全尺寸图片幻灯片

为定量化揭示任意损伤过程中冻融与盐蚀劣化作用的贡献,采用归一化处理方法分别定义冻融与盐蚀劣化因子,可通过下式计算：

$D_{1} = \frac{c_{0} - c_{i - 0}}{c_{0} - c_{i - j}}$ ,

(1)

$D_{2} = \frac{c_{i - 0} - c_{i - j}}{c_{0} - c_{i - j}}$ 。

(2)

式中 D₁为冻融劣化因子,亦即冻融损伤占比；D₂为盐蚀劣化因子,亦即盐蚀损伤占比； $c_{0}$ 表示初始状态试样黏聚力； $c_{i - 0}$ 表示冻融循环次数为i、含盐量为0.0 %时黏聚力； $c_{i - j}$ 表示冻融循环次数为i、含盐量为j %时黏聚力。（其中i取1,2,5,10；j取0.5,1,1.5）。

进一步由式（1）,（2）,可推得以下关系式：

$D_{1} + D_{2} = 1$ 。

(3)

式（3）表明,任意损伤过程归一化处理后,冻融和盐蚀劣化因子D₁,D₂即分别定量化揭示了冻融和盐蚀损伤的贡献,其物理意义直观明了。

图8（a）所示为冻融与盐蚀劣化因子的比值随冻融循环次数变化规律。从图中可以看出,冻融与盐蚀劣化因子的比值随冻融次数增加逐渐增大,亦即冻融作用对试样强度劣化的影响逐渐增强；但随着冻融次数持续增大,其增速逐渐减小,趋于一个稳定数值。值得注意的是,冻融作用下含盐量较高试样（1.5%）冻融与盐蚀劣化因子的比值始终小于1,亦即冻融劣化效应弱于盐蚀劣化效应。分析其原因,由前述盐蚀劣化因子D₂的计算公式可知,对于含盐量较高试样,其计算值显著增大,从而导致冻融与盐蚀劣化因子的比值急剧减小。

图 8 冻融与盐蚀劣化因子比值的变化规律曲线

Figure 8. Variation of ratio of freeze-thaw deterioration factor to that by salt erosion

下载: 全尺寸图片幻灯片

冻融与盐蚀劣化因子的比值随含盐量变化规律如图8（b）所示。由图可见,冻融与盐蚀劣化因子的比值随含盐量增大逐渐减小,亦即盐蚀作用对试样强度劣化的影响逐渐增强；随着含盐量持续增大,其衰减速率逐渐减小,趋于一个较小数值。此外,仅冻融循环1次条件下,冻融劣化因子与盐蚀劣化因子的比值始终小于1,亦即冻融劣化效应弱于盐蚀劣化效应。分析其原因,由前述冻融劣化因子D₁的计算公式可知,对于仅冻融循环1次试样,其计算值显著减小,从而导致冻融与盐蚀劣化因子的比值亦显著减小。

3.2 CT扫描试验结果分析

（1）冻融循环次数对细观结构的影响规律

图9所示为Na₂SO₄含量1.0%试样冻融过程不同截面的CT细观结构图像。图中白色区域代表试样密度较大处,黑色区域代表试样中孔洞、裂隙发育的位置。从图中可以看出,冻融循环条件下不同扫描断面的细观结构损伤表现出相似的变化规律：冻融前黄土试样内部存在一定的孔洞及微裂纹且无序分布,从而导致试样产生一定的初始损伤效应。随着冻融循环次数增加,受冰水相变的影响,微小孔洞逐渐扩大,且有互相贯通之势；初始微裂纹逐渐发展成裂隙,其长度及宽度均有不同程度扩展。

图 9 不同冻融循环次数下CT图像

Figure 9. CT images under different freeze-thaw cycles

下载: 全尺寸图片幻灯片

图10所示为CT数ME值随冻融循环次数变化规律曲线。由图可见,冻融初期阶段试样的CT数ME值衰减速率较大,随着冻融次数持续增大,CT数ME值逐渐趋于稳定,表现出减速衰减特征。上述变化规律反映出冻融初始阶段试样细观结构损伤演化速率较大,亦即试样内部微裂隙增长较快；多次冻融循环后,试样细观结构损伤劣化作用减弱。这很好的解释了前述三轴破坏偏应力和黏聚力随冻融循环次数的变化关系（图4（a）,5（a））。

图 10 CT数ME值与冻融循环次数关系

Figure 10. Variation of ME with freeze-thaw cycles

下载: 全尺寸图片幻灯片

（2）Na₂SO₄含量对细观结构的影响规律

图11所示为冻融循环5次条件下不同Na₂SO₄含量盐渍原状黄土试样c截面CT细观结构扫描图像。由图可见,冻融作用下不含盐试样内部大孔隙或裂隙发育相对较少。随着含盐量增加,CT图像灰度有一定程度的加深,特别是图像下部裂隙有一定的扩展和发育,出现新的细小孔隙及裂隙,盐蚀劣化程度较不含盐试样显著增大；对含盐量为1.5%试样,裂缝扩展程度较大,且中部产生较大孔隙。

图 11 不同含盐量下CT图像

Figure 11. CT images of specimens with different salt contents

下载: 全尺寸图片幻灯片

图12所示为CT数ME值随含盐量变化曲线。由图可见,冻融条件下CT数ME值随含盐量增大呈现出近似线性或加速衰减特征。上述变化规律反映出冻融条件下试样细观结构损伤演化速率随含盐量增大有增大的趋势,亦即试样内部微裂缝扩展速率逐渐增大,盐蚀劣化作用增强。这亦很好的解释了前述三轴破坏偏应力和黏聚力随含盐量的变化规律（图4（b）,5（b））。值得注意的是,未经受冻融循环（N=0）试样的CT数ME值随含盐量增加无显著变化,这是由于未冻融条件下试样内部Na₂SO₄易溶盐无相态变化,盐蚀劣化作用无法产生,因而CT数ME值无明显变化。

图 12 CT数ME值与含盐量关系

Figure 12. Variation of ME with salt content

下载: 全尺寸图片幻灯片

（3）CT细观损伤变量D_ME

根据CT原理^[14]可以得到

$ρ = ρ_{0} (1000 + H) / (1000 + H_{0})$ 。

(4)

式中 ρ为试样的密度；ρ₀为初始状态下试样的密度；H,H₀分别为与ρ和ρ₀所对应的试样CT数ME值。

令

$Δ ρ = ρ - ρ_{0}$ ,

(5)

将式（4）代入式（5）,则有

$Δ ρ = (\frac{1000 + H}{1000 + H_{0}} - 1) \times ρ_{0}$ 。

(6)

根据密度损伤变量公式^[15],可以得到CT细观结构损伤变量D_ME的表达式如下：

$D_{ME} = - \frac{1}{m_{0}^{2}} \frac{Δ ρ}{ρ_{0}}$ ,

(7)

式中,m₀为CT机空间分辨率。

将式（6）代入式（7）,则可以最终得到基于CT数ME值的损伤变量D_ME表达式：

$D_{ME} = \frac{1}{m_{0}^{2}} \times (\frac{H_{0} - H}{1000 + H_{0}})$ ,

(8)

式中,D_ME数值越大表示试样损伤幅值越大,D_ME=0代表试样初始状态。

图13所示为CT细观损伤变量D_ME变化规律。从图中可以看出,细观损伤变量随着冻融循环次数增加逐渐增大,但增速逐渐减缓,表现出减速劣化特性；细观损伤变量随着含盐量增加近似线性或加速增大,表现出等速或加速劣化特性。

图 13 细观损伤变量D_ME变化规律曲线

Figure 13. Variation of mesoscopic damage variable

下载: 全尺寸图片幻灯片

进一步根据细观损伤变量D_ME随冻融循环次数变化曲线,发现损伤变量D_ME与冻融循环次数N具有如下关系：

$D_{ME} (N) = \frac{N}{a + b N}$ ,

(9)

式中,N为冻融循环次数,a,b均为拟合参数。

将损伤变量D_ME随冻融循环次数变化结果按N/D_ME–N的关系进行拟合分析,结果如图14（a）所示。从图中可以看出,二者近似呈线性关系,其中拟合参数a为直线的截距,b为直线的斜率。

图 14 细观损伤变量D_ME拟合分析

Figure 14. Fitting analysis of mesoscopic damage variable

下载: 全尺寸图片幻灯片

对式（9）进一步求导可得

$\frac{d(D_{ME})}{d N} = \frac{a}{{(a + b N)}^{2}}$ 。

(10)

在曲线的起始点,N=0,则式（10）可表示为

$I_{D_{ME}} = \frac{1}{a}$ ,

(11)

式中, $I_{D_{ME}}$ 即为损伤变量D_ME的初始斜率。

当N→∞时,从式（9）中可得出

${(D_{ME})}_{ult} = \frac{1}{b}$ ,

(12)

式中, ${(D_{ME})}_{ult}$ 为损伤变量D_ME的极限值,表示试样达到冻融破坏极限时损伤变量的峰值。

由此可以看出,a代表损伤变量 $D_{ME}$ 的初始斜率 $I_{D_{ME}}$ 的倒数；b为损伤变量 $D_{ME}$ 极限值 ${(D_{ME})}_{ult}$ 的倒数。

图14（b）所示为拟合参数与含盐量变化关系。从图中可以看出,拟合参数a,b随含盐量可近似表示为线性变化关系。因此,进一步考虑含盐量的影响,建立拟合参数a,b与含盐量关系式：

$a = α η + a_{0}$ ,

(13a)

$b = β η + b_{0}$ ,

(13b)

式中,η为试样含盐量（%）；a₀,b₀和α,β分别为图中直线的截距与斜率。

将式（13）代入式（9）中,可得到细观损伤变量D_ME在不同冻融循环次数及含盐量下的多变量演化方程：

$D_{ME} (N, η) = \frac{N}{(α η + a_{0}) + (β η + b_{0}) N}$ ,

(14)

式中,所有参数均通过试验结果拟合得到,a₀=4.941,b₀=3.715,α=-0.915,β=-1.087。

图14（c）是含盐量分别为0%,0.5%,1.0%,1.5%下的损伤变量试验值与计算值对比图。从图中可以看出,损伤变量试验值与计算值均匀分布于直线y=x两侧,拟合相关性较好,表明该模型可较好预测Na₂SO₄盐渍原状黄土冻融过程细观结构损伤演化规律。

4. 讨论

Na₂SO₄盐渍原状黄土宏观损伤特征表现为强度劣化,劣化过程中必然伴随着试样内部细观裂隙的产生和发展,CT扫描试验可以定量化研究试样细观结构损伤演化规律。基于此,分别定义基于黏聚力的宏观损伤变量λ_C和CT数ME值的细观损伤变量D_ME,以对比分析宏细观损伤演化规律的相互关系。CT细观损伤变量 $D_{ME}$ 已通过前述式（8）得到,基于黏聚力宏观损伤变量λ_C可通过下式计算：

$λ_{C} = \frac{c_{0} - c_{i - j}}{c_{0}}$ ,

(15)

式中,λ_C数值越大表示试样损伤幅值越大,λ_C=0代表试样初始状态,λ_C=1代表完全损伤状态； $c_{0}$ , $c_{i - j}$ 的含义如前述公式（1）和（2）所示。

宏观损伤变量λ_C和细观损伤变量D_ME变化规律如图15所示。由图可见,宏细观损伤变量表现出相似的变化规律：冻融初始阶段宏细观损伤变量均增速较大,随着冻融次数持续增加,逐渐趋于稳定,表明冻融作用会导致试样的减速劣化效应。此外,宏细观损伤变量随含盐量增加其增速均逐渐增大,表明盐蚀作用会导致试样的加速劣化效应。综上,宏细观损伤变量表现出的一致变化规律表明CT细观结构损伤变量准确揭示了宏观三轴剪切力学强度指标的劣化机理。

图 15 宏细观损伤变量变化规律曲线

Figure 15. Variation of macroscopic and mesoscopic damage variables

下载: 全尺寸图片幻灯片

5. 结论

（1）冻融作用对应力–应变曲线的类型及特征无明显影响,均表现为应变硬化型。破坏偏应力随着冻融次数增加逐渐减小,但降低速率逐渐减小,表现出减速劣化特征；冻融条件下破坏偏应力随着含盐量增加逐渐减小,且表现出线性或加速劣化特征。

（2）黏聚力呈现与破坏偏应力相似的劣化特征；内摩擦角变化幅值较小,无显著变化规律。

（3）冻融与盐蚀劣化因子的比值随冻融次数增加逐渐增大但增速逐渐减小,随含盐量增大逐渐减小且衰减速率逐渐减小。

（4）CT数ME值表现出与破坏偏应力和黏聚力相似的劣化规律；构建了CT细观损伤变量演化方程,可较好预测冻融过程细观结构损伤演化规律。

（5）宏细观损伤变量表现出一致的变化规律,表明CT细观结构损伤变量准确揭示了宏观三轴剪切强度指标的劣化机理。

致谢: 感谢美国国家地质调查局—科罗拉多矿业大学岩土工程联合实验室提供的Denver膨胀土土样。

图 1 SRS-150型环剪仪

Figure 1. SRS-150 residual ring shear system

下载: 全尺寸图片幻灯片

图 2 改装后的环剪仪

Figure 2. Modified residual ring shear system

下载: 全尺寸图片幻灯片

图 3 剪应力-剪切位移关系

Figure 3. Relationship between shear stress and shear displacement

下载: 全尺寸图片幻灯片

图 4 残余强度-有效竖向压应力关系

Figure 4. Relationship between residual shear strength and effective normal stress

下载: 全尺寸图片幻灯片

图 5 不同试验方法获得的荆门中膨胀土残余强度

Figure 5. Residual shear strengths of Jingmen expansive soil with medium swelling potential determined by different test types

下载: 全尺寸图片幻灯片

图 6 不同试验方法获得的荆门弱膨胀土残余强度

Figure 6. Residual shear strengths of Jingmen expansive soil with low swelling potential determined by different test types

下载: 全尺寸图片幻灯片

图 7 残余强度内摩擦角与活动度A及黏粒含量CF的关系^[4]

Figure 7. Composite relationships showing dependence of residual friction angle on soil compositions as represented by activity and clay size fraction^[4]

下载: 全尺寸图片幻灯片

图 8 荆门中膨胀土残余强度的Stark公式预测结果

Figure 8. Predicted results of residual strength of Jingmen expansive soil with medium swelling potential by Stark's equations

下载: 全尺寸图片幻灯片

图 9 荆门弱膨胀土残余强度的Stark公式预测结果

Figure 9. Predicted results of residual strength of Jingmen expansive soil with low swelling potential by Stark's equations

下载: 全尺寸图片幻灯片

图 10 比表面积SSA与黏粒含量CF的关系^[12]

Figure 10. SSA as a function of clay fraction for smectite compared to moderately active illite and inactive kaolinite clay^[12]

下载: 全尺寸图片幻灯片

图 11 残余强度内摩擦角与比表面积SSA的关系^[12]

Figure 11. Residual friction angles as a function of SSA^[12]

下载: 全尺寸图片幻灯片

表 1 试验用土的物理性质指标

Table 1 Physical property indexes of test soils

土样名称	土粒相对质量密度^[8]	颗粒含量/%^[8]				液限/%^[8]	塑限/%^[8]	塑性指数^[8]	活动度A	USCS定名^[13]	比表面积/(m²·g^-1) ^[14]	CEC(NH₄⁺)/(mmol·kg^-1) ^[8]	标准吸湿含水率/%^[15]	自由膨胀率/%^[1]	膨胀潜势^[1]
土样名称	土粒相对质量密度^[8]	> 0.075 mm	0.005~0.075 mm	0.002~0.005 mm	< 0.002 mm	液限/%^[8]	塑限/%^[8]	塑性指数^[8]	活动度A	USCS定名^[13]	比表面积/(m²·g^-1) ^[14]	CEC(NH₄⁺)/(mmol·kg^-1) ^[8]	标准吸湿含水率/%^[15]	自由膨胀率/%^[1]	膨胀潜势^[1]
Denver膨胀土	2.73	10.0	84.0	0.7	5.3	104	48	56	10.6	MH	506.8	976	18.6	111	强
荆门黄褐色膨胀土	2.75	2.1	47.4	21.5	29.0	63	26	37	1.3	CH	236.5	309	7.3	75	中
荆门棕褐色膨胀土	2.72	11.6	42.4	26.7	19.3	42	21	21	1.1	CL	110.4	254	5.1	40	弱

下载: 导出CSV

表 2 矿物成分^[8]

Table 2 Mineral compositions 单位: %

土样名称	原生矿物	石英	钾长石	斜长石	方解石	铁白云石	重晶石	黏土矿物	蒙脱石	高岭石	伊利石	绿泥石	伊/蒙混层
Denver膨胀土	22.8	4.2	—	—	18.6	—	—	77.2	73.4	2.9	—	0.9	—
荆门黄褐色膨胀土	57.5	45.6	3.2	5.6	—	0.5	2.5	42.5	—	0.5	7.0	0.4	34.6
荆门棕褐色膨胀土	74.1	62.9	—	11.2	—	—	—	25.9	—	1.0	7.8	0.8	16.2

下载: 导出CSV

表 3 试验用土的阳离子质量摩尔浓度

Table 3 Mass molarities of cation of test soils 单位: mmol/kg

土样名称	Na⁺	K⁺	Ca²⁺	Mg²⁺
Denver膨胀土	48.4	3.3	159.0	216.0
荆门黄褐色膨胀土	6.8	6.4	85.2	25.4
荆门棕褐色膨胀土	4.5	4.2	50.4	19.3

下载: 导出CSV

表 4 SRS-150型环剪仪技术参数

Table 4 Technical parameters of SRS-150 residual ring shear system

技术参数名称	参数值
剪切速率/(°·min^-1)	0.001 ~ 360
最大连续剪切扭矩/(N·m)	300
电机峰值剪切扭矩/(N·m)	820
最大竖向荷载/kN	10
竖向位移行程/mm	50
环状试样外圆半径R₂/mm	76.05
环状试样内圆半径R₁/mm	48.40
环状试样面积/ mm²	10810.4
最大竖向压应力/ kPa	925
最大连续剪应力/ kPa	439
峰值剪应力/ kPa	1199

下载: 导出CSV

表 5 3种膨胀土残余强度

Table 5 Residual shear strength of three typical expansive soils

${\sigma '_{\text{v}}}$ /kPa	Denver强膨胀土	荆门中膨胀土			荆门弱膨胀土
${\sigma '_{\text{v}}}$ /kPa	d_h /cm	τ_f /kPa	d_h /cm	τ_f /kPa	d_h /cm	τ_f /kPa
5	17.37	1.8	19.08	1.5	6.06	2.5
10	13.50	3.1	14.27	3.1	4.44	4.8
25	21.77	6.7	9.66	6.8	14.61	11.4
50	7.43	11.0	19.14	13.1	24.71	22.8
100	2.73	22.1	22.03	24.7	4.94	44.6
200	22.81	45.1	21.68	45.0	19.23	89.2
400	11.81	77.0	13.67	85.6	16.12	174.7

下载: 导出CSV

参考文献(18)

[1]	膨胀土地区建筑技术规范: GB 50112—2013[S]. 北京: 中国建筑工业出版社, 2013. Technical Code for Buildings in Expansive Soil Regions: GB 50112—2013[S]. Beijing: China Architecture & Building Press, 2013. (in Chinese)
[2]	徐永福, 程岩, 肖杰, 等. 膨胀土滑坡和工程边坡新型防治技术研究[J]. 岩土工程学报, 2022, 44(7): 1281-1294. doi: 10.11779/CJGE202207008 XU Yongfu, CHENG Yan, XIAO Jie, et al. New prevention and control technology for expansive soil slopes[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(7): 1281-1294. (in Chinese) doi: 10.11779/CJGE202207008
[3]	叶为民, 孔令伟, 胡瑞林, 等. 膨胀土滑坡与工程边坡新型防治技术与工程示范研究[J]. 岩土工程学报, 2022, 44(7): 1295-1309. doi: 10.11779/CJGE202207009 YE Weimin, KONG Lingwei, HU Ruilin, et al. New prevention and treatment techniques and their applications to landslides and engineering slopes of expansive soils[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(7): 1295-1309. (in Chinese) doi: 10.11779/CJGE202207009
[4]	MITCHELL J K, SOGA K. Fundamentals of Soil Behavior[M]. 3rd ed. Hoboken: John Wiley & Sons, 2005.
[5]	许成顺, 王馨, 杜修力, 等. 不同黏性土的残余强度及其抗剪强度指标特性研究[J]. 岩土工程学报, 2017, 39(3): 436-443. doi: 10.11779/CJGE201703006 XU Chengshun, WANG Xin, DU Xiuli, et al. Experimental study on residual strength and index of shear strength characteristics of different clay soils[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(3): 436-443. (in Chinese) doi: 10.11779/CJGE201703006
[6]	HAYDEN C P, PURCHASE-SANBORN K, DEWOOLKAR M. Comparison of site-specific and empirical correlations for drained residual shear strength[J]. Géotechnique, 2018, 68(12): 1099-1108. doi: 10.1680/jgeot.17.P.200
[7]	范志强, 唐辉明, 谭钦文, 等. 滑带土环剪试验及其对水库滑坡临滑强度的启示[J]. 岩土工程学报, 2019, 41(9): 1698-1706. doi: 10.11779/CJGE201909014 FAN Zhiqiang, TANG Huiming, TAN Qinwen, et al. Ring shear tests on slip soils and their enlightenment to critical strength of reservoir landslides[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(9): 1698-1706. (in Chinese) doi: 10.11779/CJGE201909014
[8]	土工试验方法标准: GB/T 50123—2019[S]. 北京: 中国计划出版社, 2019. Standard for Geotechnical Testing Method: GB/T 50123—2019[S]. Beijing: China Planning Press, 2019. (in Chinese)
[9]	STARK T D, EID H T. Drained residual strength of cohesive soils[J]. Journal of Geotechnical Engineering, 1994, 120(5): 856-871. doi: 10.1061/(ASCE)0733-9410(1994)120:5(856)
[10]	Standard Test Method for Torsional Ring Shear Test to Determine Drained Residual Shear Strength of Fine-Grained Soils: ASTM D6467−21[S]. 2021.
[11]	STARK T D, HUSSAIN M. Empirical correlations: drained shear strength for slope stability analyses[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2013, 139(6): 853-862. doi: 10.1061/(ASCE)GT.1943-5606.0000824
[12]	LEHMANN P, LESHCHINSKY B, GUPTA S, et al. Clays are not created equal: How clay mineral type affects soil parameterization[J]. Geophysical Research Letters, 2021, 48(20): e2021GL095311. doi: 10.1029/2021GL095311
[13]	Standard Practice for Classification of Soils for Engineering Purposes (Unified Soil Classification System): ASTM D2487-17[S]. 2017.
[14]	鲁如坤. 土壤农业化学分析方法[M]. 北京: 中国农业科技出版社, 2000. LU Rukun. Agriculture Chemical Analysis Methods of Soil[M]. Beijing: China Agricultural Science and Technology Press, 2000. (in Chinese)
[15]	公路土工试验规程: JTG E40—2007[S]. 北京: 人民交通出版社, 2007. Test Methods of Soils for Highway Engineering: JTG E40—2007[S]. Beijing: China Communications Press, 2007. (in Chinese)
[16]	马全国, 赵鑫鑫, 韩灵灵, 等. 荆门压实弱膨胀土的排水抗剪强度[J]. 信阳师范学院学报(自然科学版), 2019, 32(2): 333-339. doi: 10.3969/j.issn.1003-0972.2019.02.028 MA Quanguo, ZHAO Xinxin, HAN Lingling, et al. Drained shear strengths of Jingmen compacted expansive soil[J]. Journal of Xinyang Normal University (Natural Science Edition), 2019, 32(2): 333-339. (in Chinese) doi: 10.3969/j.issn.1003-0972.2019.02.028
[17]	周葆春, 李颖, 马全国, 等. 低应力下膨胀土的压缩与剪切行为[J]. 信阳师范学院学报(自然科学版), 2022, 35(1): 157-167. https://www.cnki.com.cn/Article/CJFDTOTAL-XYSK202201027.htm ZHOU Baochun, LI Ying, MA Quanguo, et al. Compressibility and shear strength of expansive soil at low stresses[J]. Journal of Xinyang Normal University (Natural Science Edition), 2022, 35(1): 157-167. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XYSK202201027.htm
[18]	LADE P. The mechanics of surficial failure in soil slopes[J]. Engineering Geology, 2010, 114(1/2): 57-64.

施引文献(36)

期刊类型引用(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 . 百度学术