Influences of soluble salt content on mechanical properties of Ili undisturbed loess
-
摘要: 为揭示非饱和伊犁原状黄土的湿、载屈服特性及水量变化特性,采用非饱和土三轴仪,开展不同易溶盐含量下的控制吸力的净平均应力增减试验、控制净平均应力的吸力增减试验和控制净围压和吸力的固结剪切试验,探究易溶盐含量对非饱和伊犁黄土变形、屈服和水量变化特性的影响规律。研究结果表明:含盐量对提高颗粒间的胶结作用存在临界值,在14~20 g/kg含盐量区间存在峰值屈服净平均应力、峰值屈服剪应力,且加载增湿屈服线和剪切屈服线最大;随含盐量增加,压缩性指标降低、固结排水量减小;含盐量对水量的影响主要在土水特征干燥曲线的过渡区,但对吸力增加屈服线影响不大;盐分过量和吸力降低均可使非饱和伊犁黄土产生屈服。Abstract: In order to reveal the loading collapse yield and moisture characteristics of the unsaturated Ili undisturbed loess, a series of net mean stress increase-decrease tests on control suction, suction increase-decrease tests on control net mean stress, consolidation shear tests on control net confining pressure and suction under different soluble salt contents are carried out by using the unsaturated soil triaxial apparatus. The effects of soluble salt content on the deformation, yield and moisture of the unsaturated Ili loess are studied. The results show that there is a critical value for salt content to improve the cementation between particles. There are the peak yield net mean stress and the peak yield shear stress between 14 and 20 g/kg of salt content, and the loading collapse yield curve and the shear yield curve are the maximum. The compressibility index and moisture change at the consolidation stage decrease with the increase of salt content. The influences of salt content on moisture are mainly in the transition region of soil water characteristic drying curve, but there are few effects on the suction increase yield curve. The unsaturated Ili loess can be yielded by being excessively salting and reducing suction.
-
Keywords:
- unsaturated soil /
- Ili loess /
- soluble salt /
- yield characteristic /
- moisture characteristic
-
0. 引言
标志着现代土力学理论开端的剑桥模型将固结和剪切有效的统一起来,成为广泛应用的土的弹塑性模型。然而剑桥模型基于饱和土建立,如何将其推广到非饱和土的弹塑性模型中去,已成为土力学理论的研究热点之一[1-9]。在非饱和土本构模型研究的发展过程中等向固结曲线,即孔隙比e或比体积(
v=1+e )与平均有效应力p′ 或平均净应力ˉp 的关系在建立非饱和土本构模型中有着非常重要的作用。其中,净应力ˉp=p−ua ,平均有效应力p′=ˉp−Srs ,基质吸力s=ua−uw ,ua ,uw 分别为孔隙气压力和孔隙水压力。大量的试验表明,即使是没有经过超固结的非饱和土,在常吸力条件下的正常固结线在孔隙比和对数平均压力空间上为一条曲线[4-5]。目前描述土体等向压缩线的方法中常用的有单对数线性模型[6-10]和双对数线性模型[4-5,9]。
Roscoe等[10]在饱和正常固结土和弱超固结土的试验基础上,提出了各向等压固结中孔隙比或比体积将会沿正常固结线(NCL)变化,并且近似满足:
e=N−λlnp′, (1) 式中,e,
p′ 分别为饱和土的孔隙比和平均有效应力,N 为NCL在p′=1 kPa 时对应的孔隙比,λ 为NCL在e−lnp′ 平面的斜率。卸载和回弹过程的孔隙比变化和平均有效应力的关系表示为
e=Γ−κlnp′, (2) 式中,
Γ 为回弹线(SL)在卸载到p′=1 kPa 时对应的孔隙比,κ 为SL在e−lnp′ 平面的斜率。为了建立非饱和土的弹塑性本构模型,许多学者[6-8]将式(1),(2)直接进行了扩展,即各向等压曲线和回弹曲线分别为
NCL:e=N(s)−λ(s)lnp′, (3) SL:e=Γ(s)−κ(s)lnp′, (4) 式中,
N(s) ,λ(s) ,Γ(s) 和κ(s) 的含义与式(1),(2)相同,但均为基质吸力s=ua−uw 的函数。式(3),(4)表示的NCL和SL在e−lnp′ 平面如图1所示。与式(1),(2)类似,相同基质吸力s条件下,NCL和SL上孔隙比与平均有效应力或平均净应力的对数成线性关系,然而这种近似有诸多不合理之处[11]。如式(3)表示在等向压缩条件下,孔隙比的变化量与前期固结压力无关,只与压力的变化范围有关;式(4)表示的回弹线SL上,当
p′ 增加到一定值后,孔隙比e 将成为负数,并且当p′→0 时,e→∞ ,在理论上和实际上均不合理。为了克服e−lnp′ 线性近似的不足,考虑饱和土lne−lnp′ 为双对数线性关系式:d(lne)=−λvpd(lnp′)=−λvpdpp−uw+λvpduwp−uw, (5) 式中,
p′=p−uw 为Terzaghi有效应力,λvp 为饱和土NCL的斜率。SHENG等[5,12]采用双应力变量将式(5)扩展到非饱和土,如图2所示给出了非饱和土的各向等压固结曲线方程:
d(lne)=−λvpdˉpˉp+s−λvs(s)dsˉp+s, (6) 式中,
λvp 是一个材料常数,与吸力无关,λvs 则是一个随吸力变化的参数,饱和时λvp=λvs 。该式把体积变化划分为两部分:①净应力的改变引起的,②吸力的改变引起的,两部分引起的变形会有所不同。类似于式(6),SL方程可写为d(lne)=−κvpdˉpˉp+s−κvs(s)dsˉp+s, (7) 式中,
κvs 为回弹曲线的斜率,与吸力有关,土体达到饱和时λvp=λvs 。虽然双对数模型可以解决单一对数模型中存的许多问题,但式(6),(7)中人为的将斜率
λ 或κ 分为λvp 或κvp 以及λvs(s) 或κvs(s) 两种情况,缺乏理论依据,并且由于采用了双应力变量模型,无法直接从非饱和土过度到饱和状态[13]。1. 各向等压固结的Logistic模型
土体的等向压缩线由于结构性和各向异性等原因基本呈现平缓变化段和陡峭段两个阶段(图3)。笔者通过对大量试验数据的整理发现在等向压缩过程中
e−lnp′ 曲线符合Logistic函数。需要说明的是研究已表明土体在弹性阶段,回弹曲线的斜率几乎不受吸力变化的影响,但塑性阶段的斜率明显受吸力的影响,但是在通过比较不同应力变量,净应力e−lnˉp 和平均有效应力e−lnp′ 平面中的等向固结线的变化情况,两者的变化趋势相同,因此本文均采用e−lnp′ 表示。1.1 模型建立
Logistic函数是一种常见的S形函数,在预测种群增长方面有着广泛的应用,其微分形式如下[14]
dNdt=kN(1−Nb)。 (8) 式中 N,t分别为所研究种群的数量和时间;
k ,b 为待定的模型参数。考虑到Logistic函数可以有效地描述种群数量的增长起初阶段为指数增长,然后随着变得饱和,增加变慢,最后达到成熟时增加停止的规律,这与土体的正常固结曲线相一致。因此,在半对数坐标系e−lnp′ 中引入Logistic函数来描述土体的正常固结特性,根据式(8)其微分形式可重写为d(lne)d(lnp′)=k(1−eb)。 (9) 可看出,当式(8),(9)中的
1/b=0 ,即为双对数线性模型。虽然双对数线性模型中考虑了孔隙变化与当前孔隙比的关系,但认为孔隙比的变化量仅仅与压力的对数变化范围有关,而未考虑当前压力的大小。对式(9)进行积分并整理可得
e=b1+aexp[−kln(p′)]=b1+ap′−k, (10) 式中,
a 为积分常数。对于卸载回弹线SL,仍采用半对数线性模型,但与SHEN[7]不同,不再人为的将SL斜率分为
κvp 和κvs(s) 两种状态。实际上由于吸力对SL斜率影响很小,因此考虑如下关系:ded(lnp′)=κ, (11) 式中,
κ 为卸载回弹线在e−lnp′ 平面的斜率。式(10),(11)给出了含有3参数的等压固结线Logistic模型及含有1参数的卸载回弹半对数线性模型。
1.2 模型参数确定
通过对式(10)进行分析可以得到,当等向压力
p′=0 时,参数b=e0 即为土体的初始孔隙比。当等向压力p′=1 kPa 时,参数a=b−e1e1=e0e1−1 ,其中e1 为土体在p′=1 kPa 时的孔隙比。由于在压力较小时孔隙比不容易准确测定,因此也可用拟合的方法确定。根据式(10)整理可得
ln(b−ee)=lna−klnp′。 (12) 式(12)表明,在
ln(b−ee) 和lnp′ 平面内为线性关系,且直线的斜率和截距分别为k 和lna ,如图4所示,这样就可以很方便地确定模型中的参数k 和a 了。式(11)中SL斜率κ 亦可采用类似的方法获得。1.3 先期固结压力
为了获得非饱和土各向等压固结线上的前期固结压力
p′c ,采用与传统相类似的三步骤方法[15],通过几何关系确定出坐标点B,C(图3)。式(10)的一阶导数和二阶导数分别为
ded(lnp′)=kabexp(−klnp′)[1+aexp(−klnp′)]2, (13) d2ed(lnp′)2=k2abexp(−klnp′)[aexp(−klnp′)−1][1+aexp(−klnp′)]3。 (14) 由式(14)可得到曲线的反弯点B点坐标为
(lna/k,b/2) ,将其代入式(13)中可得到曲线在B点处的斜率为kb/4 。则过B点的切线BC可表示为
e=kb4(lnp′−lnak)+b2。 (15) 考虑
e=e0=b ,且与切线方程式(15)相交,即可得到交点C处的前期固结压力为lnp′c=4+lnak−b2。 (16) 2. 分析与讨论
为了验证本文提出的非饱和土等向固结曲线模型,与Sharma等[16]、Bellia等[17]、Burland等[18]的试验结果进行比对分析。图5~7为在不同基质吸力s条件下的试验数据和本文模型预测值的对比,图中并给出了相应的模型参数。从图可看出,Logstic函数可以很好地描述土体在等向压缩固结条件下的孔隙比变化趋势。对于非饱和土,模型(10)中的参数a,k以及前期固结压力
p′c 和土中基质吸力s相关。图8(a)~8(c)依次给出了与Sharma[16]、Bellia[17]和Burland等[18]的试验数据相应的模型参数lna,k和p′c 随吸力的变化曲线。从图中可以看出,随着基质吸力s的变化,参数lna和k并没有一个统一的规律,这与其它模型的固结曲线斜率随吸力的变化结论相一致,而前期固结压力p′c 随着吸力的增大呈增加趋势。采用Burland等[18]的在基质吸力分别为0,400 kPa时的试验数据,图9,10分别对本文提出的Logistic模型和
e−lnp′ 半对数模型以及lne−lnp′ 双对数模型进行了比对。从图中可以看出Logistic模型较之可以更好地描述土体的等向固结行为。对于卸载回弹曲线SL,采用
e−lnp′ 半对数线性模型,即式(11)可以进行很好的描述。图11给出了Sharma等[16]的回弹数据与预测值的对比,可以看出非饱和土中吸力对回弹曲线的斜率影响不大。图 11 卸载回弹曲线对比Figure 11. Comparison between experimental data of Sharma et al[16] and Eq.(11)3. 结论
考虑到土体的体积变形量不仅与荷载的变化范围有关,还与当前的荷载大小以及孔隙比大小有关,通过引入Logistic函数,建立了土体在等向压缩条件下的Logistic曲线模型,对于卸载回弹曲线,采用了半对数线性模型。通过与已有试验数据的分析和对比,发现本文提出的等向压缩曲线模型和回弹曲线模型可以很好的描述非饱和土的体变特性,并且模型中的各参数方便获取。
-
表 1 试验方案及加荷顺序(以一个含盐量为例)
Table 1 Test schemes and loading sequence (taking a salt content as an example)
试验类型 应力状态 加荷过程/kPa 净平均应力循环 初始吸力 5/ 50/ 100/ 150 净平均应力 5→100→20→300→100→400 吸力循环 初始吸力 5 净平均应力 100 吸力 5→100→25→150 固结剪切 吸力 5/ 50/ 100/ 150 净围压 100/ 200/ 300 注: 净平均应力循环试验中净平均应力小于100 kPa时,荷载变化20 kPa/级;净平均应力大于100 kPa时,荷载变化50 kPa/级;吸力循环试验中荷载变化25 kPa/级。表 2 试样排水量的量测值与校正值比较
Table 2 Comparison between measured and corrected values of amount of water discharge from specimens
试验类型 含盐量θ/(g·kg-1) 量测值/cm3 校正值/cm3 差值/cm3 相对误差/% 吸力循环试验(p=100 kPa) 5 33.53 32.5 1.03 3.18 14 36.57 35.64 0.93 2.61 26 33.79 32.17 1.62 5.04 表 3 不同含盐量下的系数值
Table 3 Coefficients under different salt contents
系数 5 g/kg 8 g/kg 14 g/kg 20 g/kg 26 g/kg κ100均值 0.008 0.009 0.009 0.009 0.009 κ300均值 0.016 0.018 0.019 0.019 0.016 λs/(10-3) 1.3 0.8 0.4 0.7 0.8 表 4 应力交点处的
εv 和εw Table 4 Values of
εv andεw at intersection of stress类别 εv /%εw /%s/kPa 0 50 100 150 0 50 100 150 p/kPa 100 100 100 100 100 100 100 100 ① 9.7 3.5 2.5 2.4 7.8 1.6 3.0 0.4 ② 12.0 17.1 23.5 26.4 0.2 9.4 28.4 34.9 ③ 13.9 6.0 4.4 5 10.9 0.6 0.2 0.1 注: ①、②、③分别表示净平均应力循环、吸力循环和固结剪切试验。 -
[1] 张爱军, 邢义川, 胡新丽, 等. 伊犁黄土强烈自重湿陷性的影响因素分析[J]. 岩土工程学报, 2016, 38(增刊2): 117-122. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2016S2019.htm ZHANG Ai-jun, XING Yi-chuan, HU Xin-li, et al. Influence factors of strong self-weight collapsibility of Ili loess[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(S2): 117-122. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2016S2019.htm
[2] 冯忠居, 乌延玲, 成超, 等. 板块状盐渍土的盐溶和盐胀特性研究[J]. 岩土工程学报, 2010, 32(9): 1439-1442. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201009027.htm FENG Zhong-ju, WU Yan-ling, CHENG Chao, et al. Salt-dissolution and salt-heaving characteristics of plate-like saline soil[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(9): 1439-1442. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201009027.htm
[3] 闫亚景, 文宝萍, 黄志全. 可溶盐对兰州非饱和重塑黄土抗剪强度的影响[J]. 岩土力学, 2017, 38(10): 2881-2887. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201710015.htm YAN Ya-jing, WEN Bao-ping, HUANG Zhi-quan. Effect of soluble salts on shear strength of unsaturated remoulded loess in Lanzhou city[J]. Rock and Soil Mechanics, 2017, 38(10): 2881-2887. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201710015.htm
[4] WANG Y G, ZHANG A J, REN W Y, et al. Study on the soil water characteristic curve and its fitting model of Ili loess with high level of soluble salts[J]. Journal of Hydrology, 2019, 578: 1-10.
[5] 康安栋, 王春艳, 张方涛, 等. 硫酸盐渍土变形特性试验研究[J]. 科学技术与工程, 2015, 15(18): 211-214, 230. https://www.cnki.com.cn/Article/CJFDTOTAL-KXJS201518037.htm KANG An-Dong, WANG Chun-yan, ZHANG Fang-tao, et al. The experimental study of deformation behaviors of sulphate salt soil[J]. Science Technology and Engineering, 2015, 15(18): 211-214, 230. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KXJS201518037.htm
[6] 颜荣涛, 赵续月, 于明波, 等. 盐溶液饱和黏土的等向压缩特性[J]. 岩土力学, 2018, 39(1): 129-138. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801017.htm YAN Rong-tao, ZHAO Xu-yue, YU Ming-bo, et al. Isotropic compression characteristics of clayey soil saturated by salty solution[J]. Rock and Soil Mechanics, 2018, 39(1): 129-138. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801017.htm
[7] 徐安花, 房建宏. 盐渍土抗剪强度变化规律的研究[J]. 交通研究, 2005(11): 54-58. https://www.cnki.com.cn/Article/CJFDTOTAL-JTBH200511014.htm XU An-hua, FANG Jian-hong. Study on laws of anti-shearing intensity of salty soil[J]. Transport Research, 2005(11): 54-58. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JTBH200511014.htm
[8] 刘威. 罗布泊盐渍土静力特性试验研究[D]. 乌鲁木齐: 新疆农业大学, 2012. LIU Wei. Study on Static Characteristics Test of Saline Soil in Lop Nur[D]. Urumqi: Xinjiang Agricultural University, 2012. (in Chinese)
[9] 蔡正银, 吴志强, 黄英豪, 等. 含水率和含盐量对冻土无侧限抗压强度影响的试验研究[J]. 岩土工程学报, 2014, 36(9): 1580-1586. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201409005.htm CAI Zheng-yin, WU Zhi-qiang, HUANG Ying-hao, et al. Influence of water and salt contents on strength of frozen soils[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(9): 1580-1586. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201409005.htm
[10] 陈正汉, 孙树国, 方祥位, 等. 多功能土工三轴仪的研制及其应用[J]. 后勤工程学院学报, 2007, 23(4): 1-5. https://www.cnki.com.cn/Article/CJFDTOTAL-HQGC200704004.htm CHEN Zheng-han, SUN Shu-guo, FANG Xiang-wei, et al. Development and application of multi-function triaxial apparatus for soil[J]. Journal of Logistical Engineering University, 2007, 23(4): 1-5. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HQGC200704004.htm
[11] ALONSO E E, GENS A, JOSA A. A constitutive model for partially saturated soils[J]. Géotechnique, 1990, 40(3): 405-430.
[12] 陈正汉. 重塑非饱和黄土的变形、强度、屈服和水量变化特性[J]. 岩土工程学报, 1999, 21(1): 82-90. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC901.016.htm CHEN Zheng-han. Deformation, strength, yield and moisture change of a remolded unsaturated loess[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(1): 82-90. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC901.016.htm
[13] 方瑾瑾, 邵生俊, 李荣, 等. 真三轴条件下Q3黄土的屈服特性研究[J]. 岩石力学与工程学报, 2016, 35(9): 1936-1944. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201609023.htm FANG Jin-jin, SHAO Sheng-jun, LI Rong, et al. Yield characteristics of Q3 loess in true triaxial tests[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(9): 1936-1944. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201609023.htm
-
期刊类型引用(2)
1. 谢虎,杨占杰,张伟,何超林,谢型浪,马海鑫. 基于混合蝙蝠算法的高比例风电电力系统调度方法. 可再生能源. 2023(06): 804-809 . 百度学术
2. 蔡清池,谢汉康,余剑毅. 基于固结实验推求先期固结压力合理性研究. 江科学术研究. 2020(03): 31-35 . 百度学术
其他类型引用(7)