APP下载

温度效应下的冻结砂土经验蠕变模型分析

2021-12-15马祖飞高亚飞

水力发电 2021年9期
关键词:砂土单轴斜率

杨 明,马祖飞,高亚飞,张 雯

(1.陕西省一八五煤田地质有限公司,陕西 榆林 719000;2.安徽理工大学土木建筑学院,安徽 淮南 232001)

0 引 言

随着我国隧道、深基坑等地下工程的不断发展,在地下工程建设过程中通常利用人工冻结法建立临时冻结壁,而温度是影响冻结壁强度与稳定性的主要因素,因此需要探究在低温条件下周围土体的强度特性、蠕变变形规律等物理力学行为[1- 4],以此建立合适的本构模型并对其进行预测,以保证工程土体在开挖过程中的稳定性和安全性。

近年来,许多学者针对冻土进行了大量的试验探究,通过建立相应的本构模型来描述其力学性质,并取得了一定的的研究成果。赵晓东等[5]进行不同温度梯度的三轴蠕变试验,探究了温度对蠕变变形规律的影响。郭梦圆等[6]利用粒子群算法对S-M模型进行参数识别,将模型计算曲线与试验曲线进行对比,曲线吻合度较好,为预测土体蠕变变形提供了一种新的思路。罗飞等[7]为研究冻结砂土的蠕变特性,在元件模型的基础上,引入损伤变量,建立新的蠕变模型,较好地描述蠕变的全过程,尤其是加速度蠕变阶段。陈军浩等[8]考虑传统整数阶模型的不足,在广义开尔文中引入分数阶理论,模拟出人工冻土的稳态蠕变变形规律。刘萌心等[9]通过不同温度下进行的加载试验,并在一维蠕变模型中考虑应力历史的作用,采用回归分析的方式确定模型参数,与试验结果进行对比发现具有较高的精度。马冬冬等[10]对冻结砂土进行主动围压试验,并结合应力-应变曲线各个阶段变化特点,分析砂土破坏特征形态。杜海民等[11]通过分析三轴压缩试验数据,探究含水量对屈服面的影响,建立考虑含水量影响的偏应力-剪应变增量型函数关系式,通过该函数关系计算的模拟值与试验结果吻合度较高。

本文主要在冻结砂土单轴抗压强度试验基础上,对其进行不同温度以及加载应力下的单轴压缩蠕变试验,并进行影响因素分析。引入S-M经验蠕变模型,考虑温度和加载应力的影响,对模型参数进行识别,得到改进的冻结砂土S-M经验蠕变模型,为人工冻结法更好地应用于实际工程提供一种理论支撑。

1 单轴压缩蠕变试验

1.1 试验方案

本文采用西北地区的砂土作为研究对象,将含水率为16%的重塑砂土制作成尺寸为直径50 mm、高100 mm的圆柱试样,分别在-5、-10、-15、-20 ℃这4种温度水平下,采用分级加载的方式进行单轴压缩蠕变试验。基于冻结砂土单轴抗压强度试验,每个温度下试验的单轴抗压强度峰值为σ,分级加载应力分别为0.1σ、0.2σ,具体加载应力强度见表1。从表1可知,单轴抗压强度峰值随温度的降低不断增大。试验仪器采用WDT-100型冻土压力试验机,当试验满足试验力下降20.0%、应变超过8.0%、试验时间超过10 h这3种情形之一时,试验自动停止。

表1 冻结砂土加载应力强度σ

1.2 试验结果分析

冻结砂土单轴抗压强度峰值与温度拟合关系见图1。从图1可知,温度T与加载应力强度σ之间近似呈线性变化,函数表达式为

图1 温度与加载应力强度关系

(1)

不同温度及加载应力下的冻结砂土单轴蠕变试验曲线见图2。从图2可知,当加载荷载为低应力0.1σ和0.2σ时,试验曲线均表现为衰减蠕变阶段和稳定蠕变阶段,未呈现出加速度蠕变趋势。温度对蠕变变形量影响较大,在0.2σ加载等级下时较为明显。随着温度逐渐越低,试样中的未冻水逐渐减少,抗压强度增大,导致蠕变变形量逐渐减小。同一温度下,随着加载应力的增加,蠕变变形量也相应增加。不同温度及加载应力下的蠕变曲线均在较短时间内达到稳定状态。

图2 不同温度及加载应力下的应变与时间的关系

2 冻结砂土经验蠕变模型

2.1 S-M模型

Singh等[12]总结提出了关于软黏土的S-M经验模型,以此预测土的蠕变变形特性。该模型还未应用于冻结砂土的蠕变变形中,本文尝试将其用于描述冻结砂土的单轴蠕变特性,其方程式如下

(2)

应力σ。对式(2)进行积分,不考虑初始应变,得

(3)

式中,ε为应变;σ为加载应力,令C=At1/(1-m)、B=1-m,得出S-M经验蠕变模型

(4)

2.2 模型参数识别

对不同温度下的蠕变试验曲线,取单位时间t1=1 h,对式(4)两边同时取对数,得

lnε=lnC+ασ+Blnt

(5)

图3为不同温度下lnε-lnt关系曲线。从图3可知,曲线近似呈线性关系,且同一温度下各加载应力作用下的拟合曲线斜率接近平行。经过拟合,得出不同温度以及加载应力作用下直线的平均斜率B,见表2。从表2可知,温度对斜率值B影响不大。因此,对不同温度下的斜率再取平均值,得出最终斜率为0.197 6。

图3 不同温度下lnε-lnt关系

表2 不同温度下的斜率值B

取t=t1=1 h时作为参考时间,对式(4)两边同时取对数,得

lnε=lnC+ασ

(6)

表3为不同温度下的lnε的试验结果,当温度为-5 ℃时,将表1和表3中各加载应力和对数应变数据代入式(6),得

表3 不同温度下的lnε值

-0.448 8=lnC+0.497α

(7)

0.164 7=lnC+0.994α

(8)

联立式(7)、(8),求解出C和α值。同理,可求得其他温度下对应的参数C和α值,得出不同温度下的参数值,见表4。

表4 不同温度下的C和 α值

通过拟合发现温度T和参数C值之间呈线性关系,见图4,函数关系式为

图4 温度T和C之间的关系

C=0.016 8T+0.433 7 (R2=0.966 0)

(9)

根据表4拟合发现,温度T和α值之间呈二次函数关系,见图5,函数关系式为

图5 温度T和α之间的关系

α=0.004 6T2+0.136 9T+1.787 3 (R2=0.951 7)

(10)

由此得出,考虑温度和加载应力为自变量的改进冻结砂土S-M经验蠕变模型

ε=C·eασ·t0.196 7

(11)

根据式(11)建立的改进冻结砂土S-M经验蠕变模型,计算出不同温度和加载应力作用下的冻结砂土的蠕变变形值,将计算结果与试验结果绘制出相应的曲线进行对比(见图6)。通过对比发现,该经验模型能较好地预测试验曲线整体趋势,表明该经验蠕变模型在冻结砂土应用方面具有一定的实用性。

图6 不同温度下应变-时间关系对比

3 结 语

本文以西北地区的砂土作为研究对象,通过冻结砂土单轴压缩蠕变试验,建立温度和加载应力影响下改进的冻结砂土S-M经验蠕变模型,得出以下结论:

(1)在低应力不同温度作用下,蠕变曲线表现为衰减型蠕变,曲线过程中并未涉及加速度蠕变阶段,且温度越高,加载应力越大,砂土蠕变变形量越大。

(2)模型参数通过考虑温度的影响并利用计算的方式进行识别,相对于直接拟合得出参数更具有实际物理意义。对比模型计算曲线与试验曲线发现,该模型可以较好地描述冻结砂土的衰减蠕变阶段和稳定蠕变阶段,为人工冻结法更好地应用于实际工程提供一种理论支撑。

猜你喜欢

砂土单轴斜率
水泥土换填法在粉质砂土路基施工中的应用研究
饱和砂土地层输水管道施工降水方案设计
单轴压缩条件下岩石峰后第Ⅱ种类型应力——应变曲线的新解释
废旧轮胎橡胶颗粒——黏土的单轴抗压特性
椭圆中关联斜率的一个优美性质
物理图像斜率的变化探讨
不同低温温度下砂土物理力学特性试验研究
龙之中华 龙之砂土——《蟠龙壶》创作谈
单轴旋转惯性导航系统的水平姿态角误差修正方法
求斜率型分式的取值范围