冻融循环下非饱和膨胀土一维土柱模型试验研究

崔郁雪，杨忠年，时 伟 ，凌贤长,2，涂志斌

(1.青岛理工大学 土木工程学院，山东 青岛 266033；2.哈尔滨工业大学 土木工程学院，黑龙江 哈尔滨 150090)

冻土在世界范围内分布十分广泛，全球陆地面积有23%的范围内分布着多年冻土，季节性冻土区的占比高达56%[1]．我国是世界上第三大冻土国，中国陆地面积中有75%的范围分布着多年冻土和季节性冻土，主要分布在西北、华北和东北地区[2]．其中季节性冻土冬季冻结夏季融化，是一种含有水分和冰晶体的特殊土水体系，季节性冻融循环对地表的能量和水分平衡有显著影响，我国季节性冻土的面积为5.137×106km2，占我国陆地面积的53.5%[3-4]．在这些季节性冻土区中大量分布着膨胀土等非饱和土，非饱和土是一种由气、液、固三相组成的体系，其力学性质和工程特性与饱和土相比更为复杂，极易造成工程冻害问题．在岩土工程实践中，大多数问题均属于非饱和土力学的问题，而太沙基经典土力学已经不能有效解决非饱和土的问题，因此研究非饱和膨胀土对于实际工程和理论革新具有重要意义[5]．膨胀土具有胀缩性、超固结性、裂隙性等不良特性，物理力学性质极为复杂[6-7]．在冻融的影响下，非饱和膨胀土中的水分发生相态的转变，导致土体产生冻胀融沉和水分迁移等现象，大大降低了土体的强度和稳定性，对处于季冻区的公路、铁路、水利、边坡等多种工程造成严重的破坏和经济损失[8-11]．

目前对冻融循环条件下非饱和膨胀土的物理力学性质的研究已有很多．许雷等[7,12-13]以南阳膨胀土为研究对象，探究了不同含水率、不同冻结温度和不同冻融循环次数下膨胀土力学性质的变化情况．Tang等[14]对经历了最多9次冻融循环的膨胀土进行固结不排水(CU)三轴试验，探究了其应力应变、弹性模量、孔隙水压力等力学参数的变化情况．Yang等[11]对重塑膨胀土进行动三轴试验，探究了单次冻融和多次冻融下试样竖向位移和动力参数的变化规律．尽管目前针对冻融循环下膨胀土力学特性的试验已有很多，但少有研究人员通过一维土柱模型对非饱和膨胀土进行冻融循环试验．土柱试验在岩土工程、环境科学、生态保护等领域有着广泛的应用，是研究土体中水分迁移、溶质运移、水分入渗及蒸发等现象变化规律的重要手段．覃小华等[15]利用自行开发的一维土柱装置，模拟并探讨了4种不同降雨速率下一维土柱的入渗规律，并提出了计算非饱和土体渗透系数的新方法．Huang等[16]制作了2个长度为12.5 m的土柱，通过水平放置并施以稳定饱和水流的试验条件，探究土柱内部溶质运移的情况．Muhammad等[17]进行了一系列饱和黏土土柱的室内试验，探究在低渗透介质中溶质运移的情况．

由此可见，土柱试验可以很好地探究土体中水分迁移、溶质运移等现象的变化规律，但冻融循环下非饱和膨胀土的一维土柱试验尚待研究．实际上，冻融循环作用下的土体不仅涉及温度场的问题，在冻融过程中发生的冻胀融沉还会导致水分场和力学场的变化，因此土体的冻融往往是一个多场耦合的综合性问题[18-19]．研究发现，实际环境中土体的冻结是一维冻结问题[20]，本文通过一维土柱模型，并在不同高程处安装相应传感器，通过一维冻融试验探究非饱和膨胀土在不同冻融循环次数和不同深度下土体温湿度、土压力、孔隙水压力和竖向位移的变化规律．同时，一维土柱模型使用透明的亚克力材料制成，该材料具有良好的透光性、化学稳定性和耐候性，可在试验过程中直观看到土壤结构的宏观变化情况．使用一维土柱模拟土体的冻融循环过程，可对土体在单向冻融条件下物理力学特性的变化规律和影响机制进行探究，更加贴合实际，具有重要的工程意义．

1 试验材料与方法

1.1 土样的制备

本试验所用膨胀土为弱膨胀土，其基本物理性质如表1所示．制样前首先将膨胀土进行烘干、碾碎、过1 mm的开孔筛，然后加入蒸馏水配制成含水率为30%的土样，放入密封袋中润湿一昼夜，保证样品中的水分均匀分布．试样尺寸采用直径140 mm、高400 mm、压实度为85%的圆柱形土样．

表1 土样的基本物理性质Tab.1 Basic physical properties of soil samples

1.2 模型的安装

本试验通过温湿度一体传感器、土压力传感器、应变式微型渗压计和直线位移传感器(图1)对试验过程中温湿度、土压力、孔隙水压力和竖向位移进行监测．试验所用一维土柱模型由自制的亚克力圆筒柱体、位移传感器安装板以及位移传感器导杆组成，如图2所示．柱体为单侧封口圆筒，柱身内径140 mm、高600 mm、壁厚10 mm，底部亚克力封板直径160 mm，厚10 mm．为方便传感器的埋设，在距离柱底50 mm和350 mm的位置处分别设置有两个大小不同的钻孔．其中较小的一个直径为10 mm，用于埋设土压力传感器和应变式微型渗压计；较大的一个直径为40 mm，用于安装土壤温湿度一体传感器．柱体上侧装有自制位移传感器安装板，其中每个板面上设置有4排5列共20个直径5 mm的钻孔，以便于直线位移传感器灵活安装．直线位移传感器与下部导杆通过攻丝连接，从而间接测得土柱不同深度处的竖向位移变形．

图1 传感器Fig.1 Various sensors

装配一维土柱模型时，首先在模型底面和筒壁内侧薄涂上一层凡士林，然后进行分层制样，层厚取为2 mm，层与层交界处进行刮毛处理，如图2(a)所示．当高度达到60 mm时，按照土压力传感器、应变式微型渗压计和直线位移传感器导向杆的尺寸挖出一组深度为10 mm的凹槽进行第一组传感器和导向杆的埋设(图2(b))，然后继续进行分层压实，当土体高度达到360 mm时，挖出第二组凹槽进行第二组传感器和导向杆的埋设，然后继续压实土体达到预定高度400 mm，填充完的土柱如图3(a)所示．土柱填充完成后，在柱体顶部安装上位移传感器安装板，导杆与直线位移传感器通过活动螺栓连接，然后再将二者共同固定到位移传感器安装板上，通过土体对导向杆作用产生的竖向位移间接监测冻融循环下土体竖向胀缩位移的变化情况．然后在两个直径40 mm的预留孔处进行温湿度一体传感器的安装，注意安装时需要将传感器探针部分全部插入土体之中，安装好的传感器通过485转换器与电脑连接进行温湿度数据的采集．最后，将柱顶开口处用保鲜膜进行密封防止柱体内水分的蒸发．此时，土柱各个部分已制作安装完成．其中距离柱底350 mm的温湿度一体传感器、土压力传感器、应变式微型渗压计和直线位移传感器分别记为TW1、E1、P1、D1，距离柱底50 mm的温湿度一体传感器、土压力传感器、应变式微型渗压计和直线位移传感器分别记为TW2、E2、P2、D2．

图2 模型的装配Fig.2 Installation of soil column model

图3 一维土柱模型Fig.3 One-dimensional soil column model

1.3 试验方案

本试验使用自制工程热物理试验系统对一维土柱进行冻融循环的系统试验．该设备为一个绝热、保温、保湿、高温与低温交变湿热试验室，可满足试验要求的温度值和湿度值．为了满足土柱一维冻融的目的，采用10 mm厚的橡塑保温棉对柱体底部和四周进行两层包裹．然后将包裹完成的土柱放入工程热物理试验系统中，静置24 h后进行3次冻融循环．考虑我国冬季最低气温的变化情况，本试验确定-10 ℃和-20 ℃两个典型温度为试验温度．同时在土柱中安置了温度传感器，试验前已通过预试验获得土柱中不同深度温度场在各阶段达到稳定所需的时间，并据此确定了试验过程中的冻结和融化时间，单次冻融循环方案如表2所示．

表2 单次冻融方案Tab.2 Single freeze-thaw program

2 结果与分析

2.1 温度变化及分析

温度是影响冻土强度的一个重要因素，土体的温度变化往往会导致土体中水分发生变化，从而导致土体内应力重分布及土体变形的发生[21]，是土体未冻水含量、土压力、孔隙水压力和竖向位移变化的根本原因．根据冻结温度可以判定土体是否处于冻结状态，同时也是确定寒区冻结深度和人工冻土冻结壁厚度的依据[2,22]．本试验通过埋设在土柱不同深度处的2个温湿度一体传感器对试验过程的温度数据进行实时监测，得到了整个冻融循环过程中非饱和膨胀土不同深度的温度-时间曲线，如图4所示．从图中可以看出，当温度降低到-10 ℃，土柱的温度首先迅速下降到0 ℃以下，然后出现小幅的跳跃回升，回升温度仍低于0℃．随着冻结时间的继续进行，土柱温度继续降低，降低速度较跳跃点之前变慢并逐渐趋于稳定．整个冻结阶段土柱不同位置处温度均高于此阶段的设定温度(-10 ℃)．当温度降低到-20 ℃时，温度首先快速降低然后趋于平缓，并无明显拐点出现，整个过程中土柱的温度始终高于-20 ℃．融化阶段，工程热物理试验系统温度上升到30 ℃，土柱的温度迅速上升，短时间内达到一个温度高点，但最高温度仍低于30 ℃．在人工控制冻融温度的情况下，3次冻融循环的温度变化情况大致相似．

图4 整个冻融过程中土柱温度-时间曲线Fig.4 The temperature-time curve of the soil column during the whole freeze-thaw process

一般情况下，土体的整个冻结过程可分为过冷、跳跃、恒定、递降4个阶段(图5)[23]．在冻结的第一个阶段(F1)，土体温度首先快速降低进入过冷阶段．该阶段中土体中的水分还未发生冻结，孔隙水处于一个亚稳状态，土柱中无冰晶体出现[24]．随着时间的增加，温度-时间曲线出现一个低值点，即过冷点．此时土体中的水分开始发生相变形成冰晶并释放一定热量，导致土柱内温度小幅跳跃升高，达到冻结点．从图中可以得到非饱和膨胀土的冻结温度为-1.33 ℃．冻结点过后，自由水发生冻结，冰晶稳定生长，温度保持小范围稳定，冻结过程进入恒定阶段．值得注意的是，TW1并未出现明显的稳定阶段，这是因为TW1处距离土层表面的距离很近，受温度影响大，因此在冻结的情况下土体降温速率较快，开始冻结后温度持续降低，并无明显恒定阶段．随后结合水开始冻结，进入递降阶段，土水体系逐步冻结．融化时，随着时间的增加，土体温度迅速增加，并无明显的过热或突变现象出现．

图5 湿土典型降温曲线 Fig.5 Typical cooling curve of wet soil

对比不同深度处土体的温度变化情况可以发现，TW2较TW1相比存在滞后现象．这是因为在降温时浅层土体受温度影响大、先发生冻结，随着温度的一维传导，深层土体也开始受负温的影响，而在传导时会发生能量的损耗，因此达到稳定时TW1温度比TW2略高一些．融化过程依旧存在能量损耗的现象，浅层土体升温较快，温度曲线整体发生了偏移．在实际工程中，土体表面因工程建设及其他人为因素作用下，地表对于太阳辐射热量的吸收和蒸发散热的不平衡也会导致深层土体的温度变化滞后于浅层土[25]．

2.2 含水率变化及分析

含水率是影响冻土强度的另一重要因素[26]．研究表明，由于土壤颗粒表面力和土颗粒之间空隙的存在，在0 ℃以下土体中仍有部分水未发生冻结，即土体在冻融过程中始终有未冻水存在[27-29]．未冻水的存在可显著影响冻土强度、泊松比、模量等物理力学性质，从而对工程稳定性和环境等造成显著影响[2,30,31]．本试验通过温湿度一体传感器对未冻水含量进行监测，将整个冻融过程未冻水含量的变化绘制成图．从图6可以看出，温度变化对未冻水含量影响显著．冻结开始时含水率小幅降低，随后出现一个明显的拐点，该拐点出现的时间与温度的跳跃点出现的时间大致相同，随后含水率继续降低并逐渐趋于稳定．随后，温度下降到-20 ℃，土柱含水率继续减少，但减少幅度明显降低，大约为第一冻结阶段的50%．融化时，温度升高50 ℃达到30 ℃，含水率迅速增加并逐渐趋于初始含水率．对比整个冻融过程含水率的变化情况可以看出，随着冻融循环次数的增加，含水率变化幅值增大．FT=1、2、3时，含水率变化范围分别为21%～34%、18%～33%、17%～32%，3次冻融结束后大约损耗了2%的水分．

图6 整个冻融过程中土柱含水率-时间曲线Fig.6 Soil column moisture content-time curve during the entire freeze-thaw process

(1)单次冻融对含水率的影响

受土颗粒表面能的影响，土壤中的毛细水、弱结合水、强结合水均受到土颗粒表面静电引力场的作用，且随距离越小，受到的分子引力越大，冻结温度就越低(图7)．毛细水的冻结点最高，但也略低于0 ℃；弱结合水在-0.5 ℃时外层开始发生冻结，-30～-20 ℃时才能全部冻结；强结合水冻结温度最低，在-78 ℃时仍不冻结[32]．以FT=1为例．刚开始冻结时，土柱表层有少量冰晶体产生，含水率缓慢降低2%左右．到达冻结点之后，土柱内的自由水开始大量冻结成冰晶体，含水率快速减少．随着冻结时间的增加，自由水逐步冻结完成，固态冰阻碍了土柱中水分的运移，未冻水含量降低速度变缓，从快速结晶阶段变为缓慢结晶阶段．当冻结进入到F2阶段时，温度与F1阶段相比更低，此时，在第一冻结阶段未冻结的部分自由水和距离土颗粒表面较远的弱结合水进一步向冰晶体转换，含水率进一步降低．对比两个不同冻结温度下未冻水含量的变化情况可以发现，温度越低，水分迁移的越迅速，未冻水含量降低速率就越快．融化时，土体温度迅速增加，未冻水含量也迅速增加，并无明显突变现象出现．

图7 土颗粒与水分子的相互作用Fig.7 The interaction between soil particles and water molecules

(2)多次冻融对含水率的影响

未经冻融时，土颗粒和水分分布较为均匀，土体结构处于一个相对稳定的状态．第一次冻结后，土颗粒首先冻结收缩，土体中水分形成冰晶体发生冻胀，孔隙增大，土颗粒受到挤压，土柱原有的稳定结构被破坏，如图8所示．新形成的冰晶体、未冻水和土颗粒通过空间排列组合形成一种新的冻土结构．融化时，冰晶体受热转化成水，对新产生的较大孔隙进行填充，融化完成后土体重新恢复到一个新的稳定状态．当进行第二次冻融时，由于结构发生了变化，因此在冻结时更多的水分转变成了冰晶体，导致第二次冻结后未冻水的含量比第一次低4%左右．随着冻融循环次数的增加，结构进一步遭到破坏，第三次冻结后未冻水的含量比第二次低1%左右．在冻结过程中始终存在一定含量的未冻水，这是因为膨胀土的土颗粒非常小，比表面积大，因此颗粒外层不易冻结的强结合水较多，在冻融多次后仍能保持一定的未冻水量．另外，在冻融循环过程中还存在一定量的蒸发损耗，这也是峰值随冻融循环次数增加小幅降低的原因．

图8 多次冻融对土水结构的影响Fig.8 The influence of repeated freezing and thawing on soil-water structure

(3)深度对含水率的影响

对比TW1和TW2曲线可以发现，在冻结过程中浅层土体(TW1)的未冻水含量总是比深层土体中(TW2)的未冻水含量低1%～2%左右．这是因为在一维冻融过程中，负温首先到达浅层土体并发生冻结，在毛细作用和冻结区吸引力的影响下，深层未冻结水分逐渐上升，使冻结区得到水分的补给，深层易冻结的自由水减少．随着冷端温度的垂直向下深入，冻结区逐步向下移动(图9)．同时深层土体受负温影响时间更短，在这些因素的共同作用下使得深层土体的未冻水含量高于浅层土体的未冻水含量．在融化阶段，温度迅速升高50℃，整个土柱迅速融化且不同深度的土体仍存在小幅的差距．一维冻融过程中土体可以分为主动区和被动区，其中被动区主要进行热传导，而主动区是冰透镜发展的热传导、流体流动和应力发展的耦合过程[33]，因此整个冻融过程中始终伴随着土压力和孔隙水压力的变化．

图9 土柱的冻结过程Fig.9 Freezing process of soil column

2.3 土压力变化及分析

温度和水分的变化是产生冻胀和冻胀力的根本原因．由于土体温度的变化，土体中的水分在冻结时发生迁移并产生冰晶体从而产生冻胀，融化时水分融化孔隙被压缩产生融沉，冻融循环产生的体积变化对土体的应力应变状态有显著影响[34-35]．整个冻融过程中土压力的变化情况如图10所示．

图10 整个冻融过程中土柱土压力-时间曲线Fig.10 The earth pressure-time curve of the soil column during the whole freeze-thaw process

(1)单次冻融对土压力的影响

以FT=1为例，在冻结第一阶段，温度降到-10 ℃，孔隙中的水分逐渐冻结成冰体积膨胀并对颗粒之间的孔隙进行填充，随着冰晶体的逐渐增多，土颗粒受到挤压(图8(b))．土柱原有的稳定结构发生改变，土中的应力状态发生改变，作用于一定面积水平面上的总应力会增大，其中E1增长了0.4 kPa，E2增长了1.4 kPa．随着冻结温度的继续降低和冻结时间的增加，冰晶逐渐将土颗粒包围，此时土颗粒变成了独立的“冰晶土”，颗粒与颗粒之间的作用力减小，其中E1降低了3.2 kPa，E2降低了3.4 kPa．融化时，温度迅速升高，冰晶体融化，土颗粒接触面积增加，同时土颗粒吸收自由水和受热也会发生膨胀，三者的共同作用使得土颗粒之间的作用力增大，E1和E2增加，并略大于初始压力．

(2)多次冻融对土压力的影响

FT=1、2、3时，土压力变化幅值分别为3.3%、1.8%、1.4%，呈递减趋势．这是因为土压力的大小与土颗粒骨架直接相关：未冻融时，土体结构较为密实，土颗粒之间的接触点较多，力的传导较均匀．开始冻融后，土颗粒发生重分布，土体结构先有的稳定状态发生破坏，土压力变化显著．随着冻融循环次数的增加，土体颗粒逐渐均一化、土颗粒之间的孔隙大小逐渐稳定，在反复冻融的影响下土压力变化的幅值逐渐变小．第二次冻融和第三次冻融在刚开始冻结时土压力出现了与第一次冻融不同的变化情况，这是因为经历一次冻融后的土柱中大孔隙体积增加和裂隙的出现大大降低了土颗粒之间的接触点，冻结过程中冰晶体更易发生冻结且优先对大孔隙和裂隙处进行填充，从而导致土压力降低．

(3)深度对土压力的影响

一维冻融时，土柱由浅到深依次受温度影响．冻结时，表层土体(E1)开始发生冻结并逐渐向深层深入，由于模型三向约束的条件限制，土体只能向上发生自由冻胀，而对下层产生挤压作用，使下层土颗粒接触更为紧密，深层土(E2)的土压力增加．同时深层土所承受上层土的自重更大，二者的共同作用下使得深层土(E2)的土压力大于浅层土的(E1)土压力．融化时，温度快速增加，冰晶体由浅入深融化成水，土颗粒之间的接触点增多，作用力增加，但在自重应力作用下深层土体(E2)的土压力仍大于浅层土体(E1)．多次冻融后，土体内产生了大孔隙和裂隙，冻结时冰晶体优先对大孔隙和裂隙进行填充，产生的冻胀力对其他土颗粒造成挤压，使得刚冻融时浅层土的土压力比深层的大．当冻结进入第二阶段(F2)，土柱深处也发生冻结，土颗粒逐渐被冰晶体包围，土颗粒之间的作用力减小．此时冻胀力依旧存在但变得很小，土体的自重应力占主导作用，因此冻结后期深层土(E2)的土压力大于浅层土(E1)．

2.4 孔隙水压力变化及分析

非饱和土在冻融过程中因相变会导致毛管势和吸附势发生显著变化，固体基质在冻融过程中会产生变形，二者的共同作用是冻融过程中非饱和土孔隙水压力变化的主要原因．孔隙水压力对研究冻土退化、冻胀融沉等现象的机理具有重要意义[36-37]．整个冻融过程中孔隙水压力的变化情况如图11所示．

图11 整个冻融过程中土柱孔隙水压力-时间曲线Fig.11 The pore water pressure-time curve of the soil column during the whole freeze-thaw process

(1)单次冻融对孔隙水压力的影响

以FT=1为例．在冻结的第一阶段(F1)，P1、P2随着温度的降低小幅降低了0.7 kPa左右，这是因为温度突然降低，土颗粒发生收缩开裂，有效孔隙增加，因此孔隙水压力小幅下降．由于温度一维传导的滞后性，P2比P1降低量略小．达到冻结温度之后，土柱内部形成大量冰晶体，孔隙水压力快速减小，P1、P2降低显著．当进入第二个冻结阶段时，温度降低到一个更低的温度，部分自由水和弱结合水继续向冰晶体转化，P1随温度的降低继续降低到-3 kPa，但速率变缓．而深层土体在上层土冻结后失去自由面，随着冰晶体的产生，下层土体受到冻胀压力的作用，因此P2在F2阶段小幅升高．融化阶段，冰晶体随温度的升高快速融化成水，因此P1、P2在融化阶段呈快速上升趋势．但P2在快速上升后又出现了一小段降低的情况，这是因为在一维融化的影响下上层土体先发生融化，对下层产生的冻胀力降低，孔隙水压力降低．

(2)多次冻融对孔隙水压力的影响

观察整个冻融过程可以发现，初次冻融时不同深度处的孔隙水压力随温度的降低均呈降低趋势，而多次冻融时浅层土体(P1)的孔隙水压力随温度的降低首先小幅降低，而深层土体(P2)则呈现相反的趋势．这是因为初次冻融时，土体随温度的降低发生收缩开裂，有效孔隙增加，因此不同深度土体的孔隙水压力均小幅降低．多次冻融时，前一次的冻融完成后仍存在部分未融化的冰晶体，当再次开始降温，土体温度低于冻结温度，冰晶体继续融化成水，因此P2小幅增大．观察不同次冻融对孔隙水压力的影响，可以发现第一次冻融时孔隙水压力的变化幅值最大，随着冻融循环次数的增加变化逐渐减小．以深层土体(P2)为例，FT=1时，孔隙水压力变化幅值约为-3.8～1.3 kPa；FT=2时，变化幅值为-2～1 kPa；FT=3时，变化幅值为-3～-0.5 kPa．这是因为孔隙水压力的大小与孔隙结构密切相关，第一次冻融之前土体内部孔隙处于一个稳定状态，尚未形成有效通道，具有良好的持水性能．当开始第一次冻融时，水分变化主要以原位相变为主，因此孔隙水压力变化幅度较大．随着冻融的进行，土体的孔隙结构遭到破坏，土体内部逐渐形成细小的水流通道，持水性能降低，孔隙水压力变化幅值减小．

(3)深度对孔隙水压力的影响

对比整个冻融过程P1、P2的变化情况可以发现，P1的变化波动与P2相比更为剧烈．这是因为一维冻融时，浅层土体受温度影响更为剧烈，因此温度梯度相对较大，温度变化速率更快，因此土体结构受到的扰动更为剧烈，孔隙水压力变化波动更大．另外一维冻结时，上层土体冻结过程产生的冻胀力及自重应力作用于下层土体，因此P2波动变化较缓．

2.5 位移变化及分析

当经历冻融循环的土体达到土中水的冻结点时，便逐渐形成冰晶体，导致土体体积增大，即冻胀的发生．此时从宏观层面上可以发现土体在纵向产生位移(图9)，对土体强度和稳定性造成影响，当位移超过一定的限值就会对工程造成安全隐患[38-39]．本试验通过直线位移传感器对位移数据进行监测，整个冻融循环过程中的位移-时间曲线如图12所示，其中位移为正表示隆起，位移为负表示沉降．

图12 整个冻融过程中土柱位移-时间曲线Fig.12 The displacement-time curve of the soil column during the whole freeze-thaw process

(1)单次冻融对位移的影响

在冻融作用下，非饱和膨胀土的胀缩变形主要由以下几个部分组成：土体中的水分发生相态的转变，不同相态的体积差造成土体的胀缩变形；冻结时，未冻水转化成冰晶体，土颗粒失水收缩，融化时冰晶体又转化为未冻水，土颗粒吸水膨胀；土颗粒本身及土中的矿物成分受温度的影响发生热胀冷缩．三者的共同作用造成竖向位移的变化，其中水分相态的转变对位移的影响最大，土颗粒及矿物成分的热胀冷缩带来的影响最小．

以FT=1为例，在冻结的第一个阶段(F1)，两个位置的土体中的水分依次发生冻结，未冻水转化为冰晶体造成土体发生冻胀，位移增加．由图8可知，F1阶段结束后未冻水含量并未降到最低，当进入冻结第二阶段时，部分自由水和距离土颗粒表面较远的弱结合水进一步向冰晶体转换，因此F2阶段位移进一步增加．融化时，冰晶体融化成水体积减小，位移减小．整个冻融过程土颗粒及矿物成分造成的胀缩变形可忽略不计．

不管是位移的增大还是减小阶段，竖向位移的变化情况始终以波浪形发生变化．这是因为一维土柱模型内侧和土体之间为直接接触状态，在位移变形时产生了摩擦阻力，限制了土体的竖向位移变形．

(2)多次冻融对位移的影响

观察整个冻融过程的位移-时间曲线可以发现，不同深度的土体在第一次冻融过程中位移的变化幅度最大．这是因为冻融之前的土体结构处于一个相对稳定的状态，在冻融的作用下，土体的冻胀融沉和水分迁移使得原有的孔隙结构发生破坏、土颗粒之间的连结作用发生改变，因此位移变化幅度较大．经历第一次冻融循环后土体的稳态结构已经发生破坏，土体内部形成一定数量的“水分迁移通道”，随着冻融循环次数的增加逐渐趋于稳定，位移变化幅值也逐步减小．

FT=1时，土柱因冻胀产生的位移约为6 mm，而融沉产生的位移为-4 mm，冻胀大于融沉．随着冻融循环次数的增加，冻胀位移逐渐小于融沉位移．刚开始冻融时，膨胀土土颗粒之间的粘聚力较大，具有较好的抵抗冻融变形的能力，因此第一次冻融循环时冻胀位移大于融沉位移．多次冻融后冰晶体融化造成土颗粒及孔隙重新分布，土体结构发生破坏，粘聚力减弱，土体抵抗变形能力降低，在上层土体自重的作用下进一步沉降变形，同时非饱和膨胀土中存在吸力促进土体产生固结，因此整个冻融过程中融沉位移逐渐大于冻结位移，这与王建平等的研究结果相同[20]．

(3)深度对位移的影响

当FT=1时，浅层土体(D1)的位移变化幅值明显大于深层土体(D2)．研究发现，土体表面的冻胀量是下层冻结区各土层产生的冻胀量总和，而土体表面的隆起量是冻结区冻胀量及未冻结区固结的总和．浅层土体下部土体较多，产生的冻胀量与固结量较大，因此位移变化较大，深层土体反之．随着冻融循环次数的增加，土体结构发生破坏并逐渐达到一个新的稳定结构，不同深度处位移差距减小．

3 结 论

本文通过自制一维土柱模型对非饱和膨胀土进行冻融循环试验，通过在2个不同深度处埋设温湿度一体传感器、土压力传感器、应变式微型渗压计和直线位移传感器对不同深度的温湿度、土压力、孔隙水压力和竖向位移的变化情况进行实时监测并对这些参数的变化情况进行的规律分析，得到以下几个结论：

(1)一般情况下，土体的冻结过程可分为过冷、跳跃、恒定、递降和稳定5个阶段．从过冷点开始，土体中的水开始形成冰晶，到冻结点之后冰晶稳定发展，本试验中所用非饱和膨胀土的冻结温度为-1.33℃．而土体的融化过程中并无明显的过热或突变现象．

(2)在冻结过程中，未冻水含量的变化存在一个明显的转折点，发生时间与冻结点相对应，整个冻结过程始终有未冻水存在．随着冻融循环次数的增加，未冻水含量有降低趋势．由于毛细作用和冻结缘对自由水的吸引，在冻结过程中浅层土体的未冻水含量总是比深层的低．

(3)温度场和水分场的变化使土体孔隙结构发生变化，进而使土体内部的应力状态发生改变，这是土压力和孔隙水压力变化的根本原因．

(4)整个冻融过程中，竖向位移以波浪形发生变化．FT=1时，土体的冻胀量大于融沉量．随着冻融的进行，土体的冻胀融沉和水分迁移使得原有的孔隙结构发生破坏、土颗粒之间的连结作用发生改变，粘聚力和摩擦力减弱，土体抵抗沉降的能力降低，沉降量逐渐大于冻胀量，多次冻融后土体相对于初始状态表现为整体下沉．