温度对自膨胀高聚物浆液在裂隙中扩散行为影响试验研究

贾赫扬; 李晓龙; 曹东东; 王珊珊; 桂云祥; 钟燕辉; 张蓓

doi:10.11779/CJGE20221530

温度对自膨胀高聚物浆液在裂隙中扩散行为影响试验研究 English Version

贾赫扬^{1, 2,},
李晓龙^{1, 2, ,},
曹东东³,
王珊珊³,
桂云祥^{1, 2},
钟燕辉^{1, 2},
张蓓^{1, 2}

1.
郑州大学水利与土木工程学院, 河南郑州 450001
2.
重大基础设施检测修复技术国家地方联合工程实验室, 河南郑州 450001
3.
河南水利勘测设计研究有限公司, 河南郑州 450016

基金项目:

国家自然科学基金项目 52178401

河南高校科技创新团队支持计划项目 23IRTSTHN014

河南省2021年水利科技攻关项目 72

中原英才计划项目 234200510014

详细信息

作者简介:
贾赫扬(1993—)，男，博士，主要从事高聚物裂隙注浆理论与试验等方面的研究工作。E-mail: JiahyChina@outlook.com

通讯作者:
李晓龙, E-mail: lxl1977@zzu.edu.cn

中图分类号: TU45
计量
- 文章访问数: 372
- HTML全文浏览量: 32
- PDF下载量: 79
出版历程
- 收稿日期: 2022-12-11
- 网络出版日期: 2023-05-23
- 刊出日期: 2024-03-31

Experimental study on influences of temperature on diffusion behaviors of self-expanding polymer slurry in fractures

JIA Heyang^{1, 2,},
LI Xiaolong^{1, 2, ,},
CAO Dongdong³,
WANG Shanshan³,
GUI Yunxiang^{1, 2},
ZHONG Yanhui^{1, 2},
ZHANG Bei^{1, 2}

1.
College of Hydraulic and Civil Engineering, Zhengzhou University, Zhengzhou 450001, China
2.
National Local Joint Engineering Laboratory of Major Infrastructure Testing and Rehabilitation Technology, Zhengzhou 450001, China
3.
Henan Water & Power Engineering Consulting CO., Ltd., Zhengzhou 450016, China

摘要

摘要: 为了研究温度对膨胀型聚氨酯高聚物注浆材料在裂隙中扩散特性的影响，设计了平板裂隙注浆模型试验装置，以一种工程中应用的聚氨酯高聚物浆液为对象，开展了高聚物注浆系列试验，测试了不同浆液预热温度、环境温度、注浆量条件下高聚物浆液在平板裂隙中的扩散特性，得到了不同工况下浆液扩散速率、扩散范围以及浆液和裂隙壁内温度分布随时间变化情况。分析结果表明：浆液预热温度对高聚物扩散过程具有重要影响，相同注浆量条件下，预热温度越高，浆液扩散速率越快，初凝时间越早，浆液扩散半径随着预热温度的升高逐渐增大，低温下浆液扩散半径显著低于高温下的试验结果；浆液扩散过程与其化学反应进程不完全同步，其化学反应终止时间晚于浆液初凝时间；浆液温度随时间呈现先升高然后逐渐下降、变化速率由快至慢的演化特征，同一时刻注浆孔中心处温度最高，沿径向逐渐递减；浆液升温速率和温度峰值随预热温度和环境温度的升高而增大；浆液扩散过程中高聚物与裂隙壁间存在明显热传导效应。试验结果对于高聚物裂隙注浆机理研究和注浆技术研发具有参考作用。
- 高聚物 /
- 裂隙注浆 /
- 扩散机理 /
- 温度 /
- 试验
Abstract: In order to study the influences of temperature on the fracture diffusion characteristics of expansive polyurethane polymer grouting materials, a plate fracture grouting model test device is designed. With a polyurethane polymer slurry applied in engineering as the research object, a series of polymer grouting tests are carried out. The diffusion characteristics of polymer slurry in plate fractures under different slurry preheating temperatures, ambient temperatures and grouting amounts are tested. The changes of slurry diffusion rate, slurry diffusion range, and temperature distribution of slurry and fracture wall under different working conditions are obtained. The analysis results show that the preheating temperature of slurry has an important influence on the diffusion process of polymer. Under the same grouting amount, the higher the preheating temperature, the faster the slurry diffusion rate and the earlier the initial gel time. The slurry diffusion range gradually increases with the increase of the preheating temperature. The slurry diffusion radius at low temperatures is significantly lower than the test results at high temperatures. The slurry diffusion process is not synchronized with the chemical reaction process, and the chemical reaction termination time is later than the initial gel time. The slurry temperature firstly rises and then gradually decreases with time, and the change rate is from fast to slow. At the same time, the temperature at the center of the grouting hole is the highest and gradually decreases along the radial reaction. The rising rate and peak value of slurry temperature increase with the increase of the preheating temperature and ambient temperature. There is an obvious heat conduction between the slurry and the fracture wall during the process of slurry diffusion. The test results can be used as a reference for the researches on the polymer fracture grouting mechanism and grouting technology.
- polymer /
- fracture grouting /
- diffusion characteristic /
- temperature /
- test

HTML全文

0. 引言

随着经济的飞速发展,全世界对地下空间的不断开发,地铁等地下工程常遇到砂土交错、土壁塌方等软弱地层施工难题^[1]。常规地层加固方法很难彻底防治水害、解决失控塌方等技术难题。由此,遇到一系列工程难题和生命财产损失事故,制约经济发展。地层冻结技术是冻结施工的一种新技术,在英、美、法、日、俄兴起较早^[2-3]。在中国,地层冻结技术实际于20世纪70年代初开始广泛应用于土层稳定工程中^[4-5],并逐渐在城市地下交通工程^[6]、煤炭矿山开采工程^[7]和紧急事故处理工程^[8]中发挥加固和稳定作用。与常规方法不同的是,地层冻结技术通过人工制冷使土壤均匀降温,将土壤内的水冻结成冰,从而形成冻土,可极大提高地层整体强度,保证施工安全。地层冻结制冷技术通常利用氨（-33.4℃）、干冰（-78.5℃）、液氮（-195.8℃）等进行冻结,液氮化学性质稳定,冻结速度比传统冻结方式可提高10倍,被广泛应用于地层冻结工程^[9]。

随着前沿科学和新技术不断发展,工程材料与结构在超常规温度、速度、场强和恶劣天气等极端服役环境中的力学响应规律成为研究热点^[10]。一般超低温是指-50℃以下^[11-12],在冻土工程学中,研究冻土超低温强度特性主要为了掌握冻土地基及应急工程力学强度指标、破坏及变化规律,以保证构筑物的安全稳定。采用冻结法对土层进行冻结后,由于温度极低,原有结构发生破坏,会对地层土体产生强烈的冻结作用,改变岩土的物理力学性质。因此,探究冻土在超低温状态下的物理力学性质,尤其是超低温状态下冻土抗压强度的变化对于保证工程施工安全具有重要意义。

目前,国内外针对温度对冻土抗压强度的影响研究,主要集中在0℃～-30℃的力学强度、弹性模量以及应力应变关系等方面^[13-21],但针对超低温条件下冻土的力学性质的研究尚且不足。同时,对于温度与冻土抗压强度的影响关系式尚未达成一致的认识,付伟等^[22]和姚兆明等^[23]认为冻土单轴抗压强度随温度降低呈线性增加；陈士威等^[24]和李怀鑫等^[25]认为冻土抗压强度与温度的关系可用指数函数来表示。综上所述,尽管就温度对土体力学性质的影响开展了一些研究,但超低温条件下土体的抗压强度研究认识不足,且由于温度区间较小,温度与土体抗压强度关系式未能达到共识,这些研究还不能满足超低温工况下,土层稳定工程建设的需要。经过多年的研究和发展,很多专家学者对于普通冻土做了大量的试验与理论研究^[26-29],冻结土体在低温冻结的条件下力学强度得到提高,土体内部由原本的三相复合体变为四相复合体,有别于常规土体的应力-应变关系方程,但传统冻土的应力-应变关系方程在超低温条件下是否适用有待研究。

本文在东北农业大学水利与土木工程学院低温环境模拟实验室对低液限黏土进行冻土单轴抗压强度试验,探求超低温条件下,土体抗压强度及应力-应变关系变化规律,对超低温冻土力学性质进行补充,完善冻土温度与抗压强度关系式,提出一个与冻土温度有关的复合型方程,为冻土力学的基本理论及工程实践提供理论依据与技术应用。

1. 材料与方法

1.1 试验土样

本文试样所用土壤为低液限黏土,取自黑龙江省南部引嫩工程,为保证土体的均匀性,取样深度为2 m,将取回的土样进行风干、碾压、过2 mm筛备用,按照《土工试验方法标准》（GTB—50123—2019）进行物理性质试验测试,测得基本物理性质指标如表1所示。

表 1 试样的基本性质

Table 1. Basic properties of soil samples

土壤类型	天然状态		相对质量密度G_s	液限/%	塑限/%	塑性指数I_P	颗粒级配/%
土壤类型	含水率/%	干密度 $ρ_{d}$	相对质量密度G_s	液限/%	塑限/%	塑性指数I_P	>0.05 mm	0.05~0.005 mm	<0.005mm
低液限黏土	17.67	1.71	2.73	36.5	17.5	19.0	35.50	35.00	29.50

下载: 导出CSV

| 显示表格

1.2 试样制备与试验方法

冻土单轴抗压强度试验的试样直径为61.8 mm,高度为125 mm的圆柱形试件,采用轻型击实仪控制含水率和击实次数,制成干密度为1.70 g/cm³,含水率为17 %,20 %,23 %的待测土样,共84个试样。

为防止液氮侵蚀,将试样用保鲜膜密封并冷冻,如图1（a）。为防止快速冷冻导致试样产生裂纹,将制备好的试样利用不同梯度降温策略控制试样温度（图1（b））,首先,采用常规冰柜将试样从常温降至-20℃；然后,采用低温冰柜将试样从-20℃降至-40℃；接着,采用超低温冰柜将试样从-40℃降至-80℃；最后,采用液氮将试样从-80℃至-180℃。使用WDW-100低温电子万能试验机对不同温度的试样进行冻土单轴压缩试验,试验仪器如图1（c）,加载速率为6 mm/min,试验时确保试样中心轴与荷载传感器轴线重合,设置试验方案和应力路径,进行压缩试验,直至试样破坏或者应变达到一定标准,停止试验,记录试验数据及试样破坏情况。试验过程设置两组平行样作为对照组,为减少试验误差,在不同温况下设置陪样,将温度传感器插入土样中心,连接巡检仪检测试样温度状况如图1（d）,同时在试验箱中放入液氮,用温度传感器监测箱体温度与试样温度一致,在快加载速率条件下加压破坏,确保试验数据有效可靠。

图 1 超低温冻土单轴压缩试验过程

Figure 1. Uniaxial compression test process of ultra-low temperature frozen soil

下载: 全尺寸图片幻灯片

2. 试验结果与分析

2.1 抗压强度

（1）冻土的破坏情况

图2为不同温度下冻土的破坏情况,-10℃～-20℃的土样破坏时中部呈现明显的鼓胀形态,并伴随竖向裂纹,此种破坏形态与试样中孔隙水相态有关,试样承受荷载能力较弱,土体向较高应力状态处滑移,表面土体颗粒间胶结力不足以抵抗挤压,因此产生裂缝,同时受试验装置两端限制,两端变形有限,呈腰鼓形态破坏^[30]；-30℃～-70℃时,随着土样温度降低,抵抗破坏的能力,主要来源于冰土胶结作用,土样变形不明显,随温度继续降低,一些胶结在一起的冰晶体被压碎,且端部受到土体的约束较小,两端开始产生细小裂纹；-80℃～-120℃,土样表面随温度降低裂纹逐渐增加,微小裂隙逐渐延伸,形成贯穿裂缝；-140℃～-160℃,在荷载作用下,平行于加载轴的裂隙发生贯通,将试样分割成了条块状,破坏形式发展为多条裂缝断面,土样呈劈裂破坏；温度达到-180℃时,试样在荷载作用下易爆裂为碎块,破坏后试样产生碎屑较多,呈压碎破坏。

图 2 不同温度冻土破坏形态

Figure 2. Failure modes of frozen soil under different temperatures

下载: 全尺寸图片幻灯片

（2）不同温度条件下冻土抗压强度关系

将含水率为17%,19%,21%的待测土样经不同降温策略进行降温处理后进行冻土单轴抗压强度试验,将不同温度下最大轴向应力作为抗压强度绘制如图3。试验结果表明：土温高于-80℃时,冻土的抗压强度随土温降低而增大,抗压强度与土温的绝对值成正比,而在土温低于-80℃时,抗压强度增长趋势变缓。一方面,未冻土转化为冻土时,液态水在负温环境中相变成固态冰^[31],导致了土体强度随温度而变。一般毛细水的冰点低于0℃,弱结合水在-20℃～-30℃才完全冻结,强结合水在-78℃仍不冻结^[32]。资料表明^[32]土体在过了剧烈相变区（大约为-10℃）后,存在于大孔隙和毛细管中所有的自由水已经冻结,部分弱结合水也已经冻结。随温度降低,土中冰晶增加,在冻土中形成冰胶结物,而冰胶结物内聚力在土冻结的情况下产生并依赖于温度,冰的强度随温度降低而增加,冰-壤胶结面强度增大；另一方面,受冰晶体结构的影响,即氢离子活动性降低,冰更加坚硬密实。随着温度的逐渐降低,土体中的未冻水含量趋于平稳,温度降低变化为1℃,冰内水的相变不会超过0.1%,此段的未冻水主要是土颗粒表面的薄膜水,薄膜水随着温度的降低逐渐变薄,越薄越不易冻结,所以未冻水含量逐渐趋于平稳,接近完全冻结,冰由六方晶体变为立方晶体^[33],冰强度达到极限,冰-壤的抗压强度变化不明显。

图 3 不同温度和含水率冻土抗压强度拟合曲线

Figure 3. Fitting curves of compressive strength at different temperatures and moisture contents

下载: 全尺寸图片幻灯片

土中的含水率会直接影响冻土强度的大小^[34]。如图4所示,冻结强度受土体含水率的影响,在一定范围内,含水率越高,含冰量越多,冰-壤胶结力增强,抵抗荷载的能力提升,土体的抗压强度越大。对不同含水率的抗压强度曲线进行拟合：

图 4 不同含水率降温过程弹性模量变化曲线

Figure 4. Variation curves of elastic modulus in cooling process under different water contents

下载: 全尺寸图片幻灯片

$y = a /(1+exp(b \cdot (x + c)))$ ,

(1)

式中,y为抗压强度,x为土体负温,将式中各参数如表2所示。结果可知,冻土抗压强度和负温的拟合程度较高,相关性较好,能够很好的描述温度变化下不同含水率抗压强度变化规律。

表 2 抗压强度与负温拟合曲线参数

Table 2. Parametric solutions of fitting curve between compressive strength and negative temperature

含水率/%	a	b	c	R²
17	25.25	0.061	42.18	0.9903
21	27.59	0.061	43.26	0.9924
23	30.04	0.056	44.05	0.9845

下载: 导出CSV

| 显示表格

（3）冻土的弹性模量

屈服应力极限与相应的应变之比定义为弹性模量,冻土的弹性模量通常取极限轴向应力的1/2与其对应的应变的比值^[25],即

$E = \frac{σ}{2 ε_{1 / 2}}$ ,

(2)

式中,E为冻土弹性模量（MPa）, $σ$ 为抗压强度（MPa）,ε_1/2为极限轴向应力的一半对应的应变值。图5为不同含水率冻土降温过程的弹性模量。试验结果表明,弹性模量虽然有些上下起伏,但总体趋势随温度降低,弹性模量呈上升趋势。由图3可知,以-80℃为分界点,低于-80℃时,土体的抗压强度变化不再明显增加,而弹性模量并没有明显受温度影响,还是持续上升,说明低于-80℃时随温度降低,虽然试样的抗压强度不再变化,但屈服点提前,弹性模量继续增加；且试样的含水率越大,弹性模量越大。

图 5 不同温度下冻土的应力-应变关系曲线

Figure 5. Relation curves between stress and strain of frozen soil under different temperatures

下载: 全尺寸图片幻灯片

2.2 超低温冻土的应力-应变关系曲线及方程

（1）冻土应力-应变关系曲线

将不同温度条件下,17%含水率的低液限黏土的应力-应变关系曲线绘制如图5,从图5中可以看出,温度对土体力学性质的影响十分显著^[34]。土温在-10℃时,随着土温逐渐降低,冻土中的含冰量迅速增大,存在于大孔隙和毛细管中的自由水已经冻结,应力-应变关系为黏弹塑类型,无明显的弹性屈服点,即应力-应变关系曲线呈应变硬化—软化^[33]。土温为-20℃～-80℃时,这时弱结合水已冻结,对冻结强度起主导作用的为冰中氢离子活动性这一因素,应力-应变关系为弹塑性类型,有明显的弹性屈服点和峰值强度,应力-应变关系曲线呈弹性—硬化—软化。土温低于-80℃时,土中强结合水基本冻结,未冻水含量极低,应力-应变曲线开始呈线性上升,荷载达到峰值后发生瞬间破坏,土样经历了微裂隙的萌生发展与贯通过程,应力-应变关系曲线呈明显的脆性破坏。

（2）冻土应力-应变关系方程

由以上应力-应变关系曲线可见,不同情况下,应力-应变关系曲线有很大差别,很难用同一类型的应力-应变方程来描述所有类型的应力-应变形式。

a）Vialov提出了幂函数形式应力-应变关系^[26]：

$σ = A ε^{m}$ 。

(3)

式中, $σ$ 为应力, $ε$ 为应变,A和m为参数。这种简单幂函数形式一直被广泛应用于描述冻土的应力-应变性状。对-30℃、-60℃、-80℃、-120℃、-180℃冻结温度下典型应力-应变关系曲线计算求解,求解参数如表3。

表 3 简单幂函数参数解

Table 3. Parametric solutions of simple power function

温度/ºC	A	m	R²
-30	30.54	0.510	0.86430
-60	231.79	0.840	0.91558
-80	1.9478	2.092	0.94367
-120	75889.56	2.300	0.96976
-180	21943.28	1.910	0.92273

下载: 导出CSV

| 显示表格

将计算结果和试验测得的应力-应变关系曲线进行对比,如图6所示。

图 6 简单幂函数试验值和计算值对比

Figure 6. Comparison of experimental and calculated values of simple power function

下载: 全尺寸图片幻灯片

由图6可见,通过简单幂函数计算出来的应力-应变关系曲线与试验值有明显差异,且温度越高相关性越差,不能够描述出应力-应变曲线的峰值及屈服阶段,超低温冻土应力-应变关系并非服从幂函数定律。

b）考特钠（Kondner）建议采用双曲线公式^[27]描述 $(σ_{1} - σ_{3})$ - $ε$ 试验曲线,即当 $σ_{3}$ 为常数时：

$σ_{1} - σ_{3} = \frac{ε}{a + b ε}$ 。

(4)

本试验为无侧限单轴抗压强度试验,没有围压作用,即 $σ_{3}$ 恒为0,将式（4）改写成：

$σ = \frac{ε}{a + b ε}$ ,

(5)

式中, $a$ , $b$ 为参数,所得直线在 $ε_{1} / (σ_{1} - σ_{3})$ 坐标轴上的截距就是 $a$ ,所得直线的斜率为 $b$ 。

对-30℃、-60℃、-80℃、-120℃、-180℃冻结温度下典型应力-应变关系曲线计算求解,求解参数如表4。

表 4 双曲线公式参数解

Table 4. Parametric solutions of hyperbolic equation

温度/ºC	a	b	m	n	k	R²
-30	0.074	0.002	-0.07	-1.66	0.002	0.9996
-60	-286	1343	1.809	1.749	0.184	0.9997
-80	0.0372	0.01796	3.3526	2.912	-2.048	0.9992
-120	3139	2369	4.305	5.669	-11.1	0.9991
-180	6222	57.88	-37.6	-0.26	0.007	0.9857

下载: 导出CSV

| 显示表格

将计算结果和试验测得的应力-应变关系曲线进行对比,如图7所示。

图 7 双曲线方程试验值和计算值对比

Figure 7. Comparison of experimental and calculated values of hyperbolic equation

下载: 全尺寸图片幻灯片

由图7可见,通过双曲线计算出来的应力-应变关系曲线相关性比幂函数要好,但在-180℃时R²值较低,同样不能够描述出应力-应变曲线的峰值及屈服阶段,幂函数和双曲线方程都是建立在增量广义胡克定律基础上变模量的弹性模型,无法反映土的剪胀性。

c）王丽琴等^[28]通过试验求解与验证,发现复合幂指数非线性模型可以很好的描述各种类型土体的应力-应变曲线,且计算准确度较传统非线性模型进一步提高,下面对超低温冻土应力-应变关系是否适用进行验证。王丽琴等^[28]提出函数表达式为

$q = σ_{1} - σ_{3} = [(a ε_{1}^{m} - k) e^{- b ε_{1}^{n}} + k] p_{a}$ ,

(6)

式中,q为偏应力（MPa）, $p_{a}$ 为标准大气压（ $p_{a}$ =101.3 kPa）,a,b,m,n,k均为试验参数。本试验为无侧向单轴抗压强度试验,没有围压的作用,即 $σ_{3}$ 恒为0。对-30℃、-60℃、-80℃、-120℃、-180℃冻结温度下典型应力-应变关系曲线计算求解,求解参数如表5。

表 5 复合幂指数模型参数解

Table 5. Parametric solutions of composite power exponential model

温度/ºC	a	b	R²
-30	0.0031	0.0807	0.94097
-60	0.0018	0.0166	0.94052
-80	0.3438	-0.0603	0.91577
-120	0.0037	-0.0811	0.94251
-180	0.0026	-0.0500	0.89904

下载: 导出CSV

| 显示表格

将计算结果和试验测得的应力-应变关系曲线进行对比,如图8所示。

图 8 复合幂指数模型试验值和计算值对比

Figure 8. Comparison of experimental and calculated values of composite power exponent model

下载: 全尺寸图片幻灯片

由拟合曲线与试验曲线对比图可以看出,在超低温条件下,复合幂指数模型不仅可以很好的描述单轴抗压应力-应变关系,还可以计算出单轴抗压强度极限值,但在-120℃及-180℃应力-应变关系的峰值拟合精度较低,图8所示,复合幂函数模型能够很好的描述冻土弹塑性破坏过程,对于超低温度的脆性破坏过程拟合精度偏低,复合幂指数对于破坏段的应力-应变曲线并不适应,只能反应试验破坏前冻结温度与应力-应变之间的非线性关系,模型仍有一定的局限性。

d）朱元林等^[29]根据大量的冻土单轴压缩试验资料,提出可将冻土的应力应变性状分为9类,利用冻土黏弹塑性I型（VEP-I）针对超低温条件下的冻结土体进行拟合可用下式描述：

$σ / σ_{\max} = {(ε / ε_{f})}^{n} \exp [n (1 - ε / ε_{f})]$ ,

(7)

式中, $σ$ 为应力, $σ_{\max}$ 为峰值强度, $ε$ 为应变, $ε_{f}$ 为破坏应变,n为参数。

为检验黏弹塑性方程对于超低温冻土应力-应变关系曲线的适用性,对-30℃、-60℃、-80℃、-120℃、-180℃冻结温度下典型应力-应变关系曲线计算求解,求解参数如表6。

表 6 黏弹塑性方程参数解

Table 6. Parametric solutions of viscoelastic plasticity equation

温度/ºC	$σ_{\max}$ /MPa	$ε_{f}$ /%	n	R²
-30	8.30	0.1059	0.5087	0.99866
-60	18.61	0.0550	0.8372	0.99981
-80	24.21	0.0276	2.0919	0.99991
-120	24.24	0.0301	2.3048	0.99995
-180	25.89	0.0265	1.9069	0.99995

下载: 导出CSV

| 显示表格

将计算结果和试验测得的应力-应变关系曲线进行对比,如图9所示。

图 9 不同温度黏弹塑性方程试验值和计算值对比

Figure 9. Comparison of experimental and calculated values of elastoplastic equation under different temperatures

下载: 全尺寸图片幻灯片

由图9可以看出黏弹塑性方程和试验测得的应力-应变关系曲线相关性最好,且温度越低拟合效果越好,不仅能够很好地描述冻土单轴应力-应变关系曲线,还可以计算出冻土的峰值强度。

黏弹塑性方程虽然可以描述应力-应变关系曲线,但不能反映应力-应变与冻结温度的关系。为了进一步研究冻土应力-应变和温度的关系,在原有方程的基础上引入温度函数,提出一个与超低温有关的复合方程描述冻土应力-应变关系曲线：

$ε_{f} = a \cdot exp(- b | T |)$ ,

(8)

$σ / σ_{\max} = {[\frac{ε}{ε_{f}}]}^{n} \exp [n (1 - \frac{ε}{ε_{f}})] ，$

(9)

式中,T为负温值（-10℃≥T≥-180℃）,a,b,n为参数, $σ_{max}$ 为试验实测值。再次检验冻结条件下,应力-应变和温度相关的复合方程,对方程求解,参数如表7所示。

表 7 温度复合参数解

Table 7. Parametric solutions of composite temperature model

温度/℃	$σ_{\max}$ /MPa	n	a	b	R²
-30	8.30	0.509	0.176	-0.105	0.99903
-60	18.61	0.837	1.682	-0.466	0.99981
-80	24.21	2.0919	0.8833	-0.0054	0.99992
-120	24.24	2.305	1.752	-0.009	0.99995
-180	25.89	1.907	0.031	0.0003	0.99997

下载: 导出CSV

| 显示表格

由表7可见,该方程随温度降低,其R²值越大,相关性越好。将计算结果和试验测得的应力-应变关系曲线进行对比,如图10所示。

图 10 不同温度复合方程试验值和计算值对比

Figure 10. Comparison of experimental and calculated values of composite equation under different temperatures

下载: 全尺寸图片幻灯片

由图10可见,引入温度的复合型应力-应变关系方程拟合精度较高,和传统的黏弹塑性方程拟合效果相比,拟合R²值更大,相关性更好,提出的带有温度参数的应力-应变关系方程具有较好的适用性,可以给工程技术的设计与实践提供理论支持。该复合方程可以很好的描述黏性土应力应变关系,而由于碎石土成分复杂、颗粒形状多样、级配组成多变,致使碎石土的力学特性异于土体与岩石,对于碎石土类是否可以推广还有待研究。

3. 结论

以黑龙江省南部引嫩工程低液限黏土为研究对象,对冻土力学特性和应力-应变关系进行了研究,开展了超低温不同含水率条件下的冻土单轴抗压强度试验,基于黏弹塑性方程对冻土单轴抗压应力-应变关系曲线进行计算拟合,并对黏弹塑性方程进一步引入温度函数,分析了对不同温度的复合方程应力-应变关系曲线的适用性。主要结论如下：

（1）在超低温条件下进行冻土单轴抗压强度试验,得到不同温度下土体的破坏情况,土温高于-80℃时,土体呈弹塑性破坏,低于-80℃时土体呈脆性破坏。

（2）对试验数据整理分析,土温高于-80℃时,冻土的抗压强度随土温降低而增大,低于-80℃时,抗压强度基本保持不变；对温度和抗压强度关系曲线进行拟合,相关系数R²达到0.98以上。

（3）冻土随温度降低,弹性模量均增大；且含水率越高,抗压强度越大,弹性模量越大。

（4）比较发现黏弹塑性方程可以统一描述超低温冻土应力-应变曲线的变化关系,通过引入温度函数,对超低温冻土应力-应变关系曲线拟合和参数的求解计算,提出一个与超低温冻土温度有关的复合型方程,表明提出的复合型方程可以很好地描述冻结温度与应力-应变曲线关系。

图 1 试验装置示意图

Figure 1. Schematic plot of test device

下载: 全尺寸图片幻灯片

图 2 可视化裂隙注浆平台

Figure 2. Visual fracture grouting platform

下载: 全尺寸图片幻灯片

图 3 温度采集系统

Figure 3. Temperature acquisition system

下载: 全尺寸图片幻灯片

图 4 热电偶布置示意图

Figure 4. Schematic plot of thermocouple layout

下载: 全尺寸图片幻灯片

图 5 高聚物浆液在裂隙中扩散过程

Figure 5. Fracture diffusion process of polymer slurry

下载: 全尺寸图片幻灯片

图 6 注浆量500 g，预热温度40℃，环境温度30℃不同测点处高聚物温度随时间变化曲线

Figure 6. Variation curves of polymer temperature with time at different measuring points under grouting amount 500 g, preheating temperature 40℃ and ambient temperature 30℃

下载: 全尺寸图片幻灯片

图 7 注浆量500 g，预热温度40℃，环境温度30℃不同时刻高聚物浆液温度沿径向变化曲线

Figure 7. Radial variation curves of polymer temperature with time at different time under grouting amount 500 g, preheating temperature 40℃ and ambient temperature 30℃

下载: 全尺寸图片幻灯片

图 8 不同预热温度条件下高聚物扩散半径随时间变化曲线

Figure 8. Variation curves of polymer diffusion radius with time under different preheating temperature conditions

下载: 全尺寸图片幻灯片

图 9 不同注浆量时高聚物扩散半径-预热温度变化柱状图

Figure 9. Histograms of polymer diffusion radius-preheating temperature change under different grouting amounts

下载: 全尺寸图片幻灯片

图 10 不同注浆量时高聚物扩散半径、温度随时间变化曲线

Figure 10. Variation curves of polymer diffusion radius and temperature with time under different grouting amounts

下载: 全尺寸图片幻灯片

图 11 注浆量750 g不同预热温度、不同位置处浆液温度随时间变化曲线

Figure 11. Variation curves of slurry temperature with time at different preheating temperatures and different positions under grouting amount 750 g

下载: 全尺寸图片幻灯片

图 12 注浆量750 g不同预热温度、不同时刻浆液温度沿径向分布曲线

Figure 12. Distribution curves of slurry temperature along radial direction at different preheating temperatures and different time under grouting amount 750 g

下载: 全尺寸图片幻灯片

图 13 注浆量750 g，预热温度50℃不同环境温度不同位置处浆液温度随时间变化曲线

Figure 13. Variation curves of slurry temperature with time at different ambient temperatures and different positions under grouting amount 750 g and preheating temperature 50℃

下载: 全尺寸图片幻灯片

图 14 注浆量750 g，预热温度50℃不同环境温度不同时刻浆液温度沿径向分布曲线

Figure 14. Distribution curves of temperature along radial direction at different ambient temperatures and different time under grouting amount 750 g and preheating temperature 50℃

下载: 全尺寸图片幻灯片

图 15 注浆量750 g，预热温度50℃，环境温度30℃时高聚物温度、接触面温度及裂隙壁内温度随时间变化曲线

Figure 15. Variation curves of polymer temperature, contact surface temperature and wall temperature with time under grouting amount 750 g, preheating temperature 50℃ and ambient temperature 30℃

下载: 全尺寸图片幻灯片

图 16 注浆量750 g，预热温度50℃，环境温度30℃时不同时刻高聚物温度、接触面温度及壁面内温度沿径向分布曲线

Figure 16. Distribution curves of polymer temperature, contact surface temperature and wall temperature along radial direction under grouting amount 750 g, preheating temperature 50℃ and ambient temperature 30℃

下载: 全尺寸图片幻灯片

表 1 高聚物浆液材料配比

Table 1 Proportioning of polymer slurry materials

组分	原料	质量分数/%
A组分	多亚甲基多苯基多异氰酸酯	100
B组分	硬泡聚醚多元醇体系	35
	聚酯多元醇（PS-3152）	15
	阻燃剂（磷酸三乙酯）	10
	催化剂体系	2
	表面活性剂（L6950）	0.5
	物理发泡剂	5
	化学发泡剂	0.5
	其它	32

下载: 导出CSV

表 2 试验各影响因素取值

Table 2 Values of influencing factors of tests

编号	注浆量/g	浆液预热温度/℃	环境温度/℃
^#1	500	30	20
^#2	625	40	30
^#3	750	50	—

下载: 导出CSV

参考文献(15)

[1]	边学成, 程翀, 王复明, 等. 高速铁路路基沉降高聚物注浆修复后动力性能及长期耐久性的试验研究[J]. 岩土工程学报, 2014, 36(3): 562-568. doi: 10.11779/CJGE201403020 BIAN Xuecheng, CHENG Chong, WANG Fuming, et al. Experimental study on dynamic performance and long-term durability of high-speed railway subgrade rehabilitated by polymer injection technology[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(3): 562-568. (in Chinese) doi: 10.11779/CJGE201403020
[2]	张连震, 张庆松, 刘人太, 等. 基于浆液-岩体耦合效应的微裂隙岩体注浆理论研究[J]. 岩土工程学报, 2018, 40(11): 2003-2011. doi: 10.11779/CJGE201811006 ZHANG Lianzhen, ZHANG Qingsong, LIU Rentai, et al. Grouting mechanism in fractured rock considering slurry-rock stress coupling effects[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(11): 2003-2011. (in Chinese) doi: 10.11779/CJGE201811006
[3]	FUNEHAG J, THRON J. Radial penetration of cementitious grout-Laboratory verification of grout spread in a fracture model[J]. Tunnelling and Underground Space Technology, 2018, 72(feb.): 228-232.
[4]	SUI W, LIU J, HU W, et al. Experimental investigation on sealing efficiency of chemical grouting in rock fracture with flowing water[J]. Tunnelling and Underground Space Technology incorporating Trenchless Technology Research, 2015, 50: 239-249.
[5]	李术才, 张霄, 张庆松, 等. 地下工程涌突水注浆止水浆液扩散机制和封堵方法研究[J]. 岩石力学与工程学报, 2011, 30(12): 2377-2396. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201112002.htm LI Shucai, ZHANG Xiao, ZHANG Qingsong, et al. Research on mechanism of grout diffusion of dynamic grouting and plugging method in water inrush of underground engineering[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(12): 2377-2396. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201112002.htm
[6]	刘人太, 张连震, 张庆松, 等. 速凝浆液裂隙动水注浆扩散数值模拟与试验验证[J]. 岩石力学与工程学报, 2017, 36(增刊1): 3297-3306. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2017S1020.htm LIU Rentai, ZHANG Lianzhen, ZHANG Qingsong, et al. Numerical simulation of crack grouting process of quick setting slurry with running water and its experimental verification[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(S1): 3297-3306. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2017S1020.htm
[7]	WANG K, WANG L, REN B, et al. Understanding the Effect of cementitious grouting pressure on micro-fracture permeability for rock reinforcement underground: a lab study[J]. Energies, 2020, 13: 4225. doi: 10.3390/en13164225
[8]	LIANG Y, SUI W, QI J. Experimental investigation on chemical grouting of inclined fracture to control sand and water flow[J]. Tunnelling and Underground Space Technology, 2019, 83: 82-90. doi: 10.1016/j.tust.2018.09.038
[9]	石明生, 夏威夷, 王复明, 等. 高聚物锚固体与粉土间黏结性能试验研究[J]. 岩土工程学报, 2014, 36(4): 724-730. doi: 10.11779/CJGE201404017 SHI Mingsheng, XIA Weiyi, WANG Fuming, et al. Experimental study on bond performance between polymer anchorage body and silt[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(4): 724-730. (in Chinese) doi: 10.11779/CJGE201404017
[10]	李晓龙, 王复明, 钟燕辉, 等. 自膨胀高聚物注浆材料在二维裂隙中流动扩散仿真方法研究[J]. 岩石力学与工程学报, 2015, 34(6): 1188-1197. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201506013.htm LI Xiaolong, WANG Fuming, ZHONG Yanhui, et al. Simulation on flowing and diffusing of expansible grouting material of polymer in two-dimensional fractures[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(6): 1188-1197. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201506013.htm
[11]	HAO M, WANG F, LI X, et al. Numerical and experimental studies of diffusion law of grouting with expansible polymer[J]. Journal of Materials in Civil Engineering, 2018, 30(2): 04017290. doi: 10.1061/(ASCE)MT.1943-5533.0002130
[12]	LI X, HAO M, ZHONG Y, et al. Experimental study on the diffusion characteristics of polyurethane grout in a fracture[J]. Construction and Building Materials, 2020, 273: 121711.
[13]	LIANG J, DU X, FANG H, et al. Numerical and experimental study of diffusion law of foamed polymer grout in fracture considering viscosity variation of slurry[J]. Tunnelling and Underground Space Technology, 2022, 128: 104674. doi: 10.1016/j.tust.2022.104674
[14]	石明生, 夏威夷. 聚氨酯高聚物注浆固化温度试验研究[J]. 化工新型材料, 2014, 42(7): 133-135. https://www.cnki.com.cn/Article/CJFDTOTAL-HGXC201407044.htm SHI Mingsheng, XIA Weiyi. Experimental research on curing temperature of polyurethane polymer[J]. New Chemical Materials, 2014, 42(7): 133-135. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HGXC201407044.htm
[15]	郝梅美. 自膨胀高聚物浆液在充水裂隙中的扩散机理研究[D]. 郑州: 郑州大学, 2019. HAO Mei-mei. Study on the Diffusion Law of Expansible Polymer Grout in the Water-Filled Rock Mass Fracture[D]. Zhengzhou: Zhengzhou University, 2019. (in Chinese)

施引文献(32)

期刊类型引用(13)

1.	王涛，凡红，王康任，周国庆，王亮亮. 高温冻土双屈服面统一本构模型及其验证. 岩土工程学报. 2025(01): 135-143 . 本站查看
2.	姚兆明，唐赛，昌语，李鹏辉. 冻结改良土抗压特性分数阶模型可靠性分析. 河南城建学院学报. 2025(01): 44-51 . 百度学术
3.	梁秀玲，王彬，张子浩，杨炳瑶，吴嘉骏. 考虑蠕变特性的多冷媒非均质人工冻结壁黏弹性分析. 冰川冻土. 2025(01): 163-178 . 百度学术
4.	谭智勇，王超林，龙安发. 外部水源作用下岩石液氮冻结试验研究. 岩土工程学报. 2024(02): 415-425 . 本站查看
5.	汪恩良，任志凤，王储，刘君巍，刘兴超，田野，邹猛，卢孜筱，张伟伟，姜生元. 基于灰色关联分析模拟月壤抗压强度性能试验. 吉林大学学报(工学版). 2024(07): 2015-2025 . 百度学术
6.	姚兆明，宋梓豪，陈军浩，左维亚. 人工冻土分数阶导数应力-应变指数模型参数确定及验证. 煤炭学报. 2024(S1): 285-294 . 百度学术
7.	田金博，张勇敢，鲁洋，马文鑫，刘斯宏，王柳江，刘瑾. 考虑初始饱和度影响的冻结渠坡膨胀土力学特性. 哈尔滨工业大学学报. 2024(11): 123-131 . 百度学术
8.	梁靖宇，沈万涛，路德春，齐吉琳. 考虑沉积角影响的冻结砂土单轴压缩试验研究. 岩土力学. 2023(04): 1065-1074 . 百度学术
9.	刘勤龙，李旭，姚兆明，吴永康，蔡德钩. 冻土强度特性及其主控因素综述. 冰川冻土. 2023(03): 1092-1104 . 百度学术
10.	宋梓豪，姚兆明. 人工冻土复合幂-指数非线性强度模型. 河南城建学院学报. 2023(05): 37-42 . 百度学术
11.	曾宇，白瑶，孙鹏，韩天宇. 人工冻结软黏土力学特性试验研究. 工业建筑. 2023(10): 105-111 . 百度学术
12.	刘铭，鲍硕超. 热-湿-应力耦合作用下框架锚杆支撑的边坡冻融特性分析. 科学技术与工程. 2023(32): 13953-13962 . 百度学术
13.	汪恩良，田雨，刘兴超，任志凤，胡胜博，于俊，刘承前，李宇昂. 基于WOA-BP神经网络的超低温冻土抗压强度预测模型研究. 力学学报. 2022(04): 1145-1153 . 百度学术