基金项目:国家自然科学基金项目(51608436); 陕西省自然科学基础研究计划(2019JQ756); 陕西省教育厅专项科研项目(19JK0452); 中国博士后科学基金项目(2019M663648)
第一作者:刘华(1983-),男,博士,副教授,主要从事环境岩土工程方面研究.E-mail:liuhua029@xauat.edu.cn
(1.西安建筑科技大学 土木工程学院,陕西 西安 710055; 2.陕西省岩土与地下空间工程重点实验室,陕西 西安 710055)
(1.School of Civil Engineering,Xi'an Univ. of Arch. & Tech.,Xi'an 710055,China; 2.Shaanxi Key Laboratory of Geotechnical and Underground Space Engineering,Xi'an 710055,China)
undisturbed loess; freeze-thaw cycle; acid pollution; tensile strength; deterioration mechanism
DOI: 10.15986/j.1006-7930.2021.04.005
我国黄土地区因季节的更替而引发的周期性冻融循环,一定程度上改变了黄土的内部组构,从而对黄土的力学性能产生损伤,进一步影响黄土上建(构)筑物的安全性[1-3].由于气候的周期性变化,黄土体的特殊组构在降雨、蒸发和冻结、融化的交替变化下发生改变,黄土的抗拉强度也处于不断变化中,改变着黄土地表的形态特征.与此同时,工农业的发展以及城市现代化的不断建设,环境污染的问题愈来愈显著,孔隙水污染对黄土场地及边坡的稳定状态造成了一定程度的影响.因此对于冻融循环条件下污染黄土的抗拉强度的研究,对黄土场地稳定性的影响具有实际的工程意义.
近年来,研究者对于冻融循环效应下的黄土体特性研究越来越多.叶万军等[4-5]通过扫描电镜、CT扫描探究了不同冻融循环次数下黄土结构损伤,结果表明土体内部骨架连接方式发生改变,影响着冻融循环作用下黄土边坡剥落模式和结果.李宝平等[6]通过平面应变试验研究了冻融循环原状黄土的力学性质,得出了冻融的循环周期越大,土体损伤越严重.Zhou[7]基于巴西劈裂试验方法探究冻土内的冰含量对土的抗拉强度的影响,结果表明冻土的抗拉强度随着温度的降低而增加,与冻土中的冰含量紧密相关.赵茜等[8]研究了冻融循环对黄土渗透系数异性的影响,结果表明随着冻融次数的增加,黄土的渗透系数增大,各向异性受到弱化.陈涛等[9]提出冻融循环中大孔隙增多、细颗粒聚集对单轴抗压性能具有弱化影响.许健等[10-11]通过CT扫描、剪切试验,研究了原状黄土和Na2SO4盐渍黄土冻融过程劣化规律和微观结构的损伤演化机制,提出原状黄土粘聚强度劣化模型.郑勋等[12-13]探讨了冻融循环对黄土结构性影响的机理,定义了“冻融结构势”作为土体结构性的定量参数,并分析了冻融循环中土的三相组成在水分相变和迁移过程的影响,得出了冻融过程中水分变化对土颗粒及孔隙的反作用力,是对土结构性影响的根本原因之一.因此冻融循环作用对黄土结构性的损伤不仅对土体压-剪性能有影响,而且对土体的抗拉强度也有一定的影响.
在污染黄土的力学性质方面,王绪民等[14]分析了原状黄土受到酸性溶液侵蚀下黄土力学性质的演变规律,得出浸泡液酸性越强,钙质胶结物溶蚀速度越快,黏聚力越小.杨秀娟等[15]通过室内模拟的方法探究了酸性溶液对黄土工程性质的影响,结果表明随着酸溶液浓度的增大,抗剪强度降低,渗透系数增大,土体孔隙孔径增大.刘华等[16]通过渗透试验及电镜扫描探究了酸污染原状土的渗透微观特性演变规律,结果表明随着酸溶液浓度增大,污染土的渗透系数降低.朱春鹏等[17]通过三轴剪切试验和直剪试验探究了酸碱污染土的强度特性,得出了浓度和总应力强度之间的关系.基于冻融劣化和孔隙水变化的影响结果,周凤玺等[18]探究了稳定条件下非饱和土中温度-水-盐(THC)多场耦合效应,得出了温度边界、含水率、孔隙率在饱和土中各分场的演变规律.
目前多数学者都是对污染土的抗剪、抗压、渗透系数等基本力学特性进行研究,对冻融循环和孔隙水酸污染耦合条件下黄土的抗拉强度的变化及强度劣化机制的研究,而对于冻融循环作用下酸污染土的抗拉强度特性研究较少.基于此,本文选取陕西省铜川市原状黄土为研究对象,采用卧式单轴抗拉强度测试仪,探明冻融循环和酸污染耦合条件下原状黄土的抗拉强度演变结果,探讨了酸污染和冻融循环次数对黄土抗拉强度的耦合影响,进一步分析酸污染黄土在冻融循环作用下的损伤机制.对边坡剥落、地表地貌塑形以及稳定性预测,在工程建设和运营中具有重要的指导意义.
试验所用黄土试样取自陕西省铜川市耀州区某在建基坑,采用人工取土的方法取得原状土样,取土深度为2.0~3.0m.该土样为Q3黄土,土质较均匀,土体中有微量颗粒状钙质胶结物,少量植物根系和虫孔.取样位置及现场见图1,原状黄土试样的基本物理性质指标见表1,颗分曲线如图2所示.
根据《土工试验方法标准》GB/T50123-2019规定[19],原状土样制备:制备尺寸为直径39.1 mm、高80 mm的圆柱体试样,制备好的土样用保鲜膜包裹封装,备用.
试验所选酸液包括盐酸、硫酸和硝酸三种,以蒸馏水为基液配置浓度为0.1 mol/L、0.5 mol/L、1.0 mol/L和2.0 mol/L四个浓度梯度,经过调研大量文献资料,大量化工工厂以及造纸业产生的废水的排放、泄漏的未稀释的污水对土体造成的污染[20],生活垃圾常年堆积产生的酸液的排放以及部分地区酸雨的影响,在长期发展过程中,由空气挥发、地下水和雨水的稀释,致使酸污染浓度有所降低,因此所选取酸溶液低浓度为0.1 mol/L和0.5 mol/L,高浓度为1.0 mol/L和2.0 mol/L.通过抽气污染的方法,制备污染试样,抽气污染试验装置如图3所示.将4种不同浓度的污染液注入试样上端,观察气压表,调节饱和缸内气压,相同温度、气压条件下,制备完成所需试样,并用保鲜膜包裹密封静置48~72 h,进行养护,最终污染试样含水率为20%.
测试黄土抗拉强度的试验仪器和方法有多种,分为直接拉伸试验和间接拉伸试验[21-22],但对于抗拉强度测试仪器尚未形成统一标准,相关测试仪器仍在讨论和规范化过程中,并无明确标准可依.本文以铜川黄土为研究对象,采用自行装配黄土卧式单轴抗拉强度测试装置,如图4和图5所示.
试验仪器最大测量荷载为200 N,精度为0.01 N,经冻融循环污染后的试样放入圆形卡箍内后,圆形卡箍和试样之间用快干胶黏结,进行固定,保证在拉伸过程中试样不会出现拔脱,试样下方涂抹凡士林,消减试样和玻璃承台间的摩擦; 抗拉试验进行时,将百分表读数归零,施加砝码,试样受到端部拉应力,进入拉伸状态,S型拉力传感器将数据传输到数据采集器中,通过换算公式得到抗拉强度值.
为研究冻融循环次数与污染黄土抗拉强度之间的关系,选定0.1 mol/L、0.5 mol/L、1.0 mol/L和2.0 mol/L四个浓度梯度为控制点,进行0、1、3、5、10和15次冻融循环的污染土的抗拉强度试验,为了避免实验过程中的误差,保证数据的准确性和平衡性,同一个冻融循环下的每个摩尔浓度梯度制备3个试样,另外设置不同冻融循环下含水率相同的未污染试样,进行对比,因此,原状黄土抗拉强度试样总共为252个.
为真实还原土体所处实际环境,经调研相关资料,铜川属于暖温带大陆性气候,不同季节日均气温相差较大,最低气温约为-15℃,融化季节日均气温约为20℃,因此在冻融循环试验过程中,冻结融化温度分别设置为-15℃和20℃.将养护好污染均匀后的试样取出,放入冻融循环箱进行封闭保湿冻融循环试验(保证黄土体在冻结过程中含水率基本不变).每次冻融循环时间设置为24 h,以冻结12 h和融化12 h作为一次冻融循环,冻结融化温度分别为-15℃和20℃,不同浓度的污染试样分别冻融0、1、3、5、10、15次,冻融箱循环过程如图6所示.
冻融循环和酸污染耦合条件下黄土的拉应力-应变关系如图7、8、9所示,图中APFTN为冻融循环和酸污染耦合次数.从图中可以看出经过不同浓度的三种酸污染之后,黄土体的拉应力随拉应变的增加不断增大,直至试样拉伸破坏,其拉应力-拉应变曲线类型为应变硬化型,属脆性破坏; 在同一应变下,盐酸、硝酸和硫酸污染后的试样,随着冻融循环次数的增加,土体抗拉强度不断减小; 在经历3次冻融循环之后,第5、10和15次冻融循环得出的应力应变曲线较为接近,经过3次冻融循环之后,黄土抗拉强度降低,再进行冻融循环黄土时内部结构已破坏,抗拉强度变化量减小.
图7 冻融循环与硫酸污染耦合下黄土拉应力-应变曲线
Fig.7 Tension stress-strain curve of loess under the coupling of freeze-thaw cycles and sulfuric acid pollution
图8 冻融循环与硝酸污染耦合下黄土拉应力-应变曲线
Fig.8 Tension stress-strain curve of loess under the coupling of freeze-thaw cycles and nitric acid pollution
图9 冻融循环与盐酸污染耦合下黄土拉应力-应变曲线
Fig.9 Tension stress-strain curve of loess under the coupling of freeze-thaw cycles and hydrochloric acid pollution
同一冻融循环下,原状黄土污染试样的拉应力随着拉应变的增加呈非线性上升趋势,曲线斜率逐渐变缓,呈现出幂函数关系.在污染黄土的拉应力-拉应变试验关系曲线中,可以发现酸污染和未污染黄土拉伸破坏应力的对比.以盐酸溶液污染为例,经过0.1 mol/L盐酸污染后黄土抗拉强度较未污染下降了17.4%,而经过1次冻融循环后黄土抗拉强度较未污染下降了13.7%; 经过2.0 mol/L盐酸污染后黄土抗拉强度较未污染下降了33.4%,而经过15次冻融循环后黄土抗拉强度较未污染下降了25.8%.可以得出,经过酸污染后黄土抗拉强度较未污染土冻融循环后黄土抗拉强度下降较大即酸污染对于黄土强度的影响程度大于冻融循环对黄土强度的影响.分析其原因:黄土经酸液侵蚀后严重破坏了土体本身的结构性,酸溶液与土中物质发生反应,生成可溶性和难溶性盐,土体结构的破坏程度较大; 冻融的反复进行即温度的变化,致使孔隙水形态不断改变以及孔隙水在土体中不断移动,孔隙水的流动带动土体内细颗粒在孔隙中游离,土体原生结构变化较小,黄土强度损失相对较小.
通过黄土单轴抗拉强度试验得出不同冻融循环次数下酸污染黄土的抗拉强度和酸液浓度关系曲线如图 10所示,图中N为冻融循环次数.
图 10 冻融循环下土抗拉强度与酸液浓度的关系
Fig.10 The relationship between soil tensile strength and acid concentration under freeze-thaw cycles
由图 10可以看出,经过三种酸污染之后的土的抗拉强度均有所减小,但最大值都不超过12 kPa; 当盐酸浓度从0.1 mol/L升高至2.0 mol/L时,0次冻融循环的原状黄土的抗拉强度下降了19.26%,但经过1次冻融循环之后原状黄土的抗拉强度降低了23.98%,经过3次冻融循环之后原状黄土的抗拉强度降低了21.99%,经过第5次冻融循环之后原状黄土的抗拉强度降低了23.7%,经过第10次冻融循环之后原状黄土的抗拉强度降低了20.54%.
同一冻融循环次数下盐酸、硝酸和硫酸污染的黄土抗拉强度和酸液浓度呈非线性关系; 经过盐酸和硝酸污染之后的抗拉强度随浓度的增加呈非线性的减小,经过硫酸污染之后的抗拉强度随浓度是呈非线性的增大; 根据曲线拟合可得,相同冻融循环下原状黄土抗拉强度随酸液浓度的增大呈指数函数关系; 当浓度逐渐较大时,曲线斜率逐渐变缓,经过3次冻融循环之后,抗拉强度逐渐趋近.
原状黄土的抗拉强度源于土颗粒间吸附强度、原始黏聚力和加固黏聚力[23],土颗粒间吸附强度是由毛细压力和基质吸力形成,土体内多种碳酸盐及盐类薄膜经过长时间的化学胶结作用形成黄土中加固黏聚力,胶结物质的数量、形态和成分以及含水率对土颗粒间的黏聚力有不同影响.黄土中主要包括易溶盐(Na+、Mg+、Ca2+的氯盐、碳酸盐和硫酸盐)、中溶盐(CaSO4)和难溶盐(CaCO3)组成的化学胶结物以及游离氧化物(SiO2、Fe2O3、Al2O3)[24-26].
不同酸污染黄土的抗拉强度的变化主要是因为黄土受到酸溶液的侵蚀后,土体中含有的碳酸盐及其它化学胶结物质发生溶解后,胶结力降低,产生SiO2-3、Fe2+、Fe3+、AlO2-等游离离子,一部分和盐酸、硝酸反应生成了可溶性盐和其它新的物质,另一部分被其它化学胶结物质再吸附或游离到土颗粒孔隙中生成具有吸附强度的胶结物质.在土体未被污染时,土体中的化学胶结物质数量是不变的,所含化学物质成分是相同的,当污染液的浓度增大时,溶入土中的Cl-、NO3-数量增多,土中原有胶结物质的浓度减小,生成新的可溶性盐和其他新物质,破坏了土体原有的内部结构,土颗粒间的吸附强度降低,土体的胶结作用减小,黏聚力减弱,其抗拉强度下降.其中硫酸和土体中化学胶结氧化物反应之后,低浓度的硫酸溶液中SO2-4较少,生成难溶性盐较少,土颗粒间的孔隙变大,黏结结构发生破坏,土体抗拉强度降低.随着硫酸溶液的浓度增加,生成的难溶性硫酸盐逐渐增多,生成的硫酸盐填补了土颗粒间部分孔隙,其中硫酸铁是红褐色的黏稠液体,加强了土颗粒间的黏结作用,土体抗拉强度相对有所提高.
冻融循环作用会改变土颗粒大小,破坏土体颗粒间的连结,胶结力弱化,冻结生成的冰晶,导致土颗粒挤压移动,土体结构体系发生变化[27-29],结构损伤,对黄土抗拉强度产生一定影响.通过试验数据归纳整理,得出冻融循环次数与污染原状黄土抗拉强度关系曲线如图 11所示.由图 11可知,同一种酸液浓度污染的原状黄土的抗拉强度随着冻融循环次数的增加逐渐降低,但是随着冻融循环次数的逐渐增多,污染黄土抗拉强度值降低幅度逐渐减弱,直至趋于稳定.在经历1、3次冻融循环之后,污染黄土抗拉强度变化较大,其抗拉强度值衰减了约20%,经过5次冻融循环之后,其抗拉强度衰减幅度降低.可看出污染原状黄土抗拉强度加速衰减集中在3~5次冻融循环,之后劣化特征减缓,抗拉强度逐渐趋于稳定的,可得出经开挖裸露后的黄土地层的抗拉性能在3~5个冻融周期内的变化幅值是最大的推断.
图 11 冻融循环次数与污染黄土抗拉强度关系
Fig.11 The relationship between the number of freeze-thaw cycles and the tensile strength of contaminated loess
原状黄土经过酸液污染之后,再进行反复冻融作用,土体中水分进行不断的冻结和融化,孔隙水发生迁移,迁移力不断改变土体内部孔隙大小及分布; 冻融循环过程中孔隙水形态发生变化,液-固-液态交替进行,由于水冻结成冰,在颗粒孔隙之间约束作用下产生冻胀力作用在土颗粒上,破坏了土体原生内部结构.当冻融循环达到一定次数时,由于土体内水形态的变化,在土颗粒间产生的迁移力和冻胀力导致土体的胶结力不断受损,土体中较大粒组结构遭到破坏,改变了原始土颗粒的形态,进行重组排列,初始结构破坏,土体的黏聚力降低,黄土的抗拉强度降低.同时由于酸液的侵蚀,在土体内发生了化学反应,生成了可溶性盐,离子反应方程式为
冻融循环过程中产生的可溶性盐,发生溶解-结晶-溶解交替变化,形态反复的改变和不断重复侵蚀,黄土内部结构遭到破坏,土颗粒排列形态较为松散; 经过多次冻融,土颗粒排列趋于平衡状态,土颗粒间联结强度达到平稳状态即残余抗拉强度.
基于前述冻融循环次数对污染黄土抗拉强度的影响分析,建立酸侵蚀与冻融循环耦合抗拉强度劣化模型.
酸污染黄土抗拉强度与冻融次数呈指数函数关系,可用以下指数函数公式进行拟合分析.
σt=A·eBN+C (2)
式中:σt为抗拉强度; N为冻融循环次数:A、B、C为拟合参数(见表2).
基于表2拟合参数,考虑酸污染浓度对黄土抗拉强度产生的影响,对表2数据进行拟合,采用(3)式进行拟合,拟合效果较好.
A=A1n+A2
B=B1n+B2
C=C1n+C2 (3)
式中:n为酸溶液浓度; A1、A2、B1、B2、C1、C2、为拟合参数,拟合结果见表3.
将(3)式代入(2)式,可得到酸污染原状黄土抗拉强度与冻融循环次数的关系式.
式中:n为酸溶液浓度; N为冻融循环次数.
通过已有独立试验数据和模型计算值进行对比验证,此酸侵蚀与冻融循环耦合抗拉强度劣化模型得出计算值和试验值相差较小,表明式(4)能较好反映冻融循环条件下酸污染原状黄土抗拉强度劣化特性.
本文通过配置孔隙水酸污染的原状黄土,探讨了在冻融循环效应下抗拉强度退化的演变规律,得出了以下结论:
(1)酸污染原状黄土的抗拉强度较小,随着酸液浓度的增大抗拉强度呈非线性变化态势发展.
(2)在相同污染浓度下,抗拉强度随着冻融循环次数增加而减小,经过3~5次冻融循环后,抗拉强度劣化程度较大,5次冻融循环后,抗拉强度逐渐趋于稳定.
(3)通过室内试验得出在冻融循环效应下酸污染土的拉应力-拉应变演变规律,并分析了酸污染黄土抗拉强度损伤结果,建立了酸侵蚀与冻融循环耦合抗拉强度劣化模型.