考虑土体结构劣化的高原山区桩-土体系数值分析

寇海磊 侯王相 荆皓 陈琦 李恒

摘要：高原山區土体常年受冻融循环作用，结构不断劣化，导致桩-土界面作用削弱，对构筑物稳定性产生不利影响.本文开展室内试验结合数值模拟，构建了土体在冻融循环作用下的结构劣化模型，建立了考虑土体结构劣化的水-热-力耦合数学模型并验证了其可靠性. 基于此，对服役期内高原山区桩-土体系稳定性进行了数值分析.试验结果表明，冻融循环条件下，土体渗透系数及孔隙比随冻融次数增加呈对数型增大，土体黏聚力随冻融次数增加呈指数型减小，内摩擦角随冻融次数增加呈对数型增大.服役期内，桩周土体最大冻胀位移随服役年限增加逐渐减小，最大融沉位移逐年增大，桩周土体整体呈融沉趋势；桩基冻胀位移随服役年限增加逐渐减小，融沉位移逐年增大，冻胀融沉增长速率不断减小，但桩基融沉位移大于其冻胀位移.试验结果可为高原山区桩基础设计提供技术支持与理论指导.

关键词：高原山区；桩基工程；冻融循环；土体结构劣化；水-热-力耦合模型

中图分类号：TU445；TU473文献标志码：A

Numerical Analysis of Pile-soil System in Plateau Mountainous Areas Considering Soil Structure Deterioration

KOU Hailei1，HOU Wangxiang1，JING Hao2，CHEN Qi1，LI Heng3

（1. College of Engineering，Ocean University of China，Qingdao 266100，China；2. District Power Supply Company of Lingcheng District，State Grid Shandong Electric Power Company，Dezhou 253500，China；3. Qingdao Construction Engineering Management and Service Center，Qingdao 266071，China）

Abstract：Due to the perennial freezing and thawing cycles，the soil structure in plateau mountainous area deteriorates continuously，which adversely affects the pile-soil interface and the stability of structures. In this paper，the indoor laboratory tests and numerical simulation were carried out，a freezing-thawing cycle degradation model of soil structure was established，a coupled mathematical model of water-thermal-mechanics considering soil structure degradation was then established，and its reliability was also verified. Based on this，the stability of pile-soil system in plateau mountainous area during service period was numerically analyzed. The experimental results show that the soilpermeability coefficient and pore ratio increase logarithmically with the number of freezing-thawing cycles，the soil cohesion decreases exponentially with the number of freezing-thawing cycles，and the internal friction angle increases logarithmically with the number of freezing-thawing cycles. During service period，the maximum frost heave displacement around soil gradually decreases with the increase of service life，and the maximum thaw settlement displacement also decreases. As a whole，the around soil shows a thaw settlement trend. The frost heave displacement of pile decreases gradually with time，the thaw settlement displacement increases as well，and the growth rate of frost heave and thaw settlement continuously decreases. However，the thaw settlement displacement of pile is larger than its frost heave displacement. Test results can supply some technical support and theoretical guidance for the design of pile foundation in plateau mountainous areas.

Key words：plateau mountain；pile foundation engineering；freeze-thaw cycles；deterioration of soil structure；water-thermal-mechanical coupling model

随着川藏铁路建设的持续推进，西部高原山区基础设施建设持续开展.高原山区桩基工程在服役过程中常年反复的冻融循环导致桩-土界面作用弱化，进而造成上部构筑物产生不均匀沉降、结构断裂等工程灾害，对上部结构物的稳定性带来了极大的挑战[1-4].

冻融循环作用下寒区土体的物理力学性质变化始终备受国内外学者的关注[5-8].Chamberlain等[9]通过室内试验对冻融条件下寒区土体的渗透性以及微观结构进行了研究，发现土体在冻融过程中孔隙增大，且土体的原生结构遭到破坏；Hotineanu等[10]通过室内模型试验探究了冻融条件下不同种类寒区土体力学性质变化规律，试验结果表明，冻融作用对土体的黏聚力影响较大，是影响土体稳定性的主要因素；苏永奇等[11]探究了冻融循环作用对青藏粉质黏土的力学参数影响，认为土体在冻融过程中发生了结构损伤，最终导致寒区土体力学性能产生不可逆的变化；赵茜等[12]通过室内三轴固结渗透试验探究了冻融循环条件下土体渗透性能的变化规律，结果表明土体的渗透系数随冻融循环次数的增加不断增大.虽然国内外学者对冻融条件下不同土体的物理力学性质变化进行了大量研究，但关于土体结构损伤对桩-土体系的影响及土体结构损伤模型的建立研究甚少.

国内外学者对高原山区桩-土体系相互作用进行了大量的物理模型试验[13-15].因数值手段具有便捷性，数值模拟一直受各国学者的青睐[16-19].Lai等[20]通过数值耦合模型对服役期内桩基冻胀过程进行了非线性分析；Lu等[21]提出了一种预测冻土中单桩冻胀效应的数学模型，阐明了土体冻结深度对桩身轴力以及竖向位移的影响，但其并未考虑土体结构劣化对桩土界面相互作用的影响；Aldaeef和Rayhani[22-23]研究了桩基础在冻结、非冻结黏土中的荷载传递规律，并提出了一种基于冻结温度以及黏聚力因子预测桩基承载力的方法，但该方法未考虑桩-土界面相互作用随冻融循环不断弱化的因素；邢爽等[24]基于非线性计算模型对地震作用下季冻区冻土-结构动力相互作用进行了研究，但未考虑土体结构随冻融作用的劣化特性；唐丽云等[25]探究了不同温度条件下地下水位对多年冻土区桩基承载性能的影响，但未考虑常年冻融作用对桩基承载特性的影响；陈坤等[17]对多年冻土地区灌注桩桩身温度分布特性以及桩-土体系导热过程进行了研究，但未考虑土体因温度变化发生的冻胀融沉特性.综上，已有研究虽对考虑多场耦合条件下桩-土体系承载特性变化规律进行了研究，但对冻融作用下土体结构劣化所导致的桩-土体系相互作用弱化研究较少.

本文通过开展变水头渗透试验以及温控三轴试验对冻融条件下高原山区土体力学参数劣化规律进行研究，构建了基于土体结构劣化的水-热-力耦合模型并对其可靠性进行了验证.在此基础上，以实际工程为依托，对常年反复冻融循环条件下桩周土体以及桩基冻胀融沉特性进行了数值模拟与分析.

1室内试验

1.1试验土样

试验所用土样取自西藏日喀则萨迦县至秋洛村公路工程曲洛3号中桥桩基钻孔，如图1所示.依据《土工试验方法标准》（GB/T 50123—2019）[26]对土样进行室内土工试验，确定试样土样为粉质黏土，具有一定可塑性，其基本物理性质参数指标见表1.

1.2试验过程

1.2.1变水头渗透试验

为探究冻融作用对土体渗透特性的影响，开展不同冻融循环条件下土样的变水头渗透试验.根据《土工试验方法标准》（GB/T 50123—2019）[26]制备渗透试验标准环刀土样.首先对试样进行冻融处理，依据日喀则当地工程地质资料5，设定冻结温度为-10 ℃，融化温度为20 ℃，冻结与融化时间均为12 h，以保证土样完全冻结与充分融化.整个冻融过程在TMS9018-250型冻融循环试验箱（图2）中进行，最终完成制備0次、1次、2次、4次、6次、8次、10次冻融循环后的试样.

冻融试样制备完毕后，将其置于真空饱和器中进行抽真空饱和2 h，待无气泡产生后表明土样已达到饱和状态.将饱和后的试样放入变水头渗透试验装置中，试验装置如图3所示.首先记录起始水头高度，每间隔5 min观测记录水头高度变化，每组试样连续计测3次.每组试样重复试验5次后取平均值作为其渗透系数，冻融循环条件下的渗透试验方案如表2所示.

1.2.2温控三轴试验

为研究冻融作用对土样力学特性的影响，开展不同冻融循环条件下土样的温控三轴试验.首先对制备好的试样进行冻融循环处理，在冻融循环箱中完成1次、2次、4次、6次、8次、10次冻融循环后将试样取出，采用保鲜膜密封以防止冻融过程中水分的蒸发.

三轴试验采用GDS液压动三轴测试系统（ELCTTS-2017）进行，如图4所示，该设备主要由温控装置、剪切装置和控制系统三部分组成.其可施加最大轴力为10 kN，最小轴力为0.001 kN，所产生的最大法向位移为90 mm，最小位移为0.001 mm，围压最大可施加至70 MPa.经过不同冻融循环后的试样分别在不同有效围压下进行固结不排水（CU）三轴剪切试验，设定固结不排水剪切速率为0.07 mm/min，试验方案如表3所示.

1.3结果与分析

1.3.1冻融循环对土体渗透特性的影响

冻融循环作用对土体渗透特性具有重要影响[12，28].试验土体渗透系数K_T可通过式（1）获得.

式中：A为试样的横截面积，cm²；a为变水头管的横截面积，cm²；L为试样高度，cm；t₁、t₂分别为测读水头高度的起始、终止时间，s；h₁、h₂分别为起始水头和终止水头的高度，cm.

图5表示土体渗透系数随冻融循环次数变化关系.由图5可知，土体渗透系数随冻融循环次数的增加呈对数型上升趋势.未经过冻融循环（N=0）土样，初始渗透系数K_T0=4.46×10^-6cm/s；1次凍融循环（N=1）后，土体渗透系数增大，达到6.36×10^-6cm/s，较未冻融处理土样（N=0）增大了42.6%；6次冻融循环（N=6）后，土体渗透系数达到1.05×10^-5cm/s，较未冻融处理土样增大了135%.此后，当N>6时，土体渗透系数变化趋于稳定，最终经历10次冻融循环后的土体（N=10），渗透系数达到1.2×10^-5cm/s，较未冻融处理土样（N=0）增大了169%.这主要是由于冻融循环作用下土颗粒作为土体骨架发生了不可逆的损伤破坏，导致土颗粒间的孔隙增大.

不同冻融循环条件下土体孔隙比可通过表征软黏土孔隙比与渗透系数的Taylor公式[29]求得：

式中：e和e₀分别为土体的当前孔隙比和初始孔隙比；k和k₀分别为土体的当前渗透系数和初始渗透系数.其中土体转换系数C_k=0.5e₀，土样的初始孔隙比e₀=0.25.

图6表示土体孔隙比与不同冻融循环次数之间的关系曲线.分析可知，土体孔隙比随冻融循环次数的增加呈对数型上升趋势.N=0时，土体的初始孔隙比e₀为0.249；1次冻融循环（N=1）后，土体的孔隙比增大至0.254，较未冻融处理土样增大了2%；经历6 次冻融循环（N=6）后，土体的孔隙比增大至0.261，较未冻融处理土样（N=0）增大了5%；当试样最终经历10次冻融循环过程后（N=10），土体的孔隙比增大至0.267，较初始土样（N=0）增大了7%.产生该现象的主要原因为，土体冻结过程中，孔隙水由液态相变成冰，其体积增大，这一过程中土体内部由于相变产生的应力作用使土体内的孔隙增大；而在土体融化的过程中，土体内的冰晶重新相变成液态，但此时土体骨架来不及收缩或已发生不可逆转的损伤，因而总体来看土体内的孔隙是在不断增大的.

1.3.2冻融循环对土体剪切特性的影响

图7表示不同冻融循环条件下土体固结不排水（CU）三轴剪切试验过程中应力-应变关系曲线.由图7可知，土体的抗剪强度随轴向应变的增大逐渐增大，均呈应变硬化特征.相同冻融条件下，围压对土体抗剪强度特征影响突出.以6次冻融循环后的土体应力-应变曲线为例（图7（c）），取轴向应变ε=15%时所对应的剪应力为土体的破坏强度.当σ₃= 100 kPa时，试样破坏强度为152 kPa；当σ₃=200 kPa 时，试样破坏强度为221 kPa，增大了45.4%；当σ₃= 300 kPa时，破坏强度则为335 kPa，较σ₃=100 kPa 时，其破坏强度增大了120.4%.值得注意的是，试样经历6次冻融循环后，其轴向应变发展出现了明显的超前特征.在剪应力还未发挥剪切作用时，试样就已经产生了变形，且随着冻融循环次数的增加，该特征更显著.

为更直观分析冻融循环次数对土体破坏强度的影响，绘制试样破坏强度与冻融循环次数之间的关系曲线，如图8所示.土体破坏强度随冻融循环次数增加呈现不断减小的趋势.当^3=100 kPa时，10次冻融循环后的试样（N=10）相较于未经过冻融循环试样（N=0）的破坏强度降低了40.2%；当σ₃=200 kPa 时，10次冻融循环后的试样（N=10）相比于未经过冻融循环试样（N=0）的破坏强度降低了32.1%；而当^3=300 kPa时，这一数值则降低了25.7%.寒区土体在冻融循环的过程中其强度不断降低，土体破坏强度受冻融循环影响不断减小，这表明土体结构在冻融循环过程中产生了明显的劣化.

根据摩尔-库仑理论，可计算得到不同冻融循环次数下土体破坏强度对应的黏聚力c、内摩擦角φ，汇总结果如图9、图10所示.

图9表示土体黏聚力随冻融循环次数变化关系曲线.可见，土体黏聚力随冻融循环次数的增加呈指数型下降.当N=0时，土体的初始黏聚力c₀=28.96 kPa；经历了1次冻融循环（N =1）土体的黏聚力降至25.79 kPa，较初始试样（N=0）降低了11%；经历了6 次冻融循环后（N=6）土体黏聚力降至14.45 kPa，较初始试样（N=0）降低了50.1%；而当经历了10次冻融循环后（N =10），土体的黏聚力最终降至11.05 kPa，较未冻融处理试样降低了61.8%.结果表明，第1次冻融循环后试样黏聚力下降幅度最为明显，随着冻融循环次数的增加，前6次冻融循环造成的累积效应对试样黏聚力的影响较大；后续随冻融次数的继续增加，对土体黏聚力的影响逐渐减弱.其原因主要是孔隙水冻结成冰导致土体体积膨胀，密度减小，土颗粒间的联结方式发生变化，导致土体黏聚力不断降低.随着冻融循环过程的继续发生，有部分土体持续发生着结构变化，但一部分原生土体结构已发生不可逆转的破坏，在此过程中表现为土体黏聚力的变化速率不断减小.

图10表示内摩擦角与冻融循环次数变化关系曲线.分析可知，土体的内摩擦角随着冻融循环次数的增加呈对数型缓慢增加.N=0时，土体的初始内摩擦角φ₀=16.37°；1次冻融循环过程后（N =1），土体的内摩擦角为17.31°，较未冻融处理试样（N=0）增大了5.7%；经历了6次冻融循环后（N=6），土体的内摩擦角为19.07°，较未冻融处理试样增大了16.5%；10 次冻融循环后（N =10），土体的内摩擦角为20.67°，较未冻融处理试样（N=0）增大了26.3%. 土体在冻融循环过程中，内摩擦角逐渐增大但其增长速率却在逐渐降低.这是因为孔隙水凝结成冰，孔隙体积增大，土颗粒间的有效接触面积减小，内摩擦角增大.随着冻融循环过程的继续，前一次融化过程中土颗粒深层孔隙中的冰晶体未完全融化，再一次冻结过程中孔隙水结成冰的相变过程较前一次有所削弱，最终表现为内摩擦角的增大速率逐渐降低.

通过对比分析冻融循环过程中黏聚力和内摩擦角的变化幅度，10次冻融循环后土体黏聚力降低了61.8%，而土体内摩擦角仅增大了26.3%.表明冻融循环对内摩擦角的影响远远小于其對土体黏聚力的影响.

2基于土体劣化桩-土体系数值模拟

2.1基于土体结构劣化水-热-力耦合数学模型

2.1.1土体渗透特性劣化特征

对不同冻融循环次数条件下渗透系数和孔隙比的劣化系数运用广义最小二乘法求得劣化系数的拟合曲线方程，拟合结果如图11所示.

2.1.2土体剪切特性劣化特征

由1.3.2小节的试验结果分析可得，土体的黏聚力c随冻融循环呈指数型降低，土体的内摩擦角φ随冻融循环呈对数型增加.以冻融循环次数N为变量，构建黏聚力c、内摩擦角φ的劣化模型，定义c、φ的劣化系数X_c、X_φ分别为：

式中：X_c、X_φ分别表示土体黏聚力、内摩擦角的劣化系数；c₀、φ₀分别表示未经过冻融循环处理试样的黏聚力以及内摩擦角参数值；c_n、φ_n分别表示经过n次冻融循环后试样黏聚力以及内摩擦角参数值.

将不同冻融循环次数条件下土体黏聚力和内摩擦角的劣化系数进行计算汇总，并运用广义最小二乘法进行劣化系数曲线方程的拟合，拟合结果如图12所示.

由此可推导出N次冻融循环后土体的黏聚力c_n以及内摩擦角φ_n表达式分别为：

c_n=c₀e^-0.106N（9）

φ_n=φ₀×（0.091 ln N+1.030 2）（10）

2.1.3土体水-热-力耦合控制方程修正

土体水-热-力耦合控制方程的修正，是在归纳总结国内外学者们的研究方法与理论基础上[30]，充分考虑土体劣化在桩服役期间的影响，对现有的控制方程进行的改进.

1）水分场控制方程

根据质量守恒定律建立的水分场控制方程：

式中：θ为未冻水体积分数；θ_i为冰体积分数；ρ₁为水的密度，g/cm³；ρ_i为冰的密度，g/cm³；v为水流通量；ψ为基质势；z为重力势；K_T（θ₁）为土体的渗透系数，cm/s.

2）温度场控制方程

根据能量守恒定律建立温度场控制方程：

式中：C为土体热容，kJ/（m³·K）；L为相变潜热，kJ/kg；T为温度，℃；A为导热系数，W/（m·K）.

在冻融循环过程中，温度场与水分场之间相互影响，具有强耦合关系.温度场控制方程与水分场控制方程耦合过程中存在θ₁、θ_i以及T三个未知函数，因此需要一个联系方程实现水热两场耦合并进行修正.通过引入白青波等[31]定义的固液比B（T）作为两场之间的联系方程：

式中：T_f为土体的冻结温度，℃；b为固液比系数，粉土取值0.47，黏土取值0.56.

3）应力场控制方程

由水-热两场耦合产生的体积应变作为土体内部应力场建立的控制方程如式（15）：

ε_vf=0.09（θ₀+Δ_θ-θ_u）+Δθ-e（15）

式中：θ₀为初始含水量；Δθ为水分迁移量；θ_u为未冻水含量；e为初始孔隙比.

在冻融循环过程中，孔隙比e随着冻融循环次数N的增加不断变化，将式（6）代入式/15）后可得修正后的应力场控制方程：

ε_vf=0.09（θ₀+Δθ-θ_u）+Δθ-e₀·（0.014 8ln N+1.0266）（16）

2.1.4模型验证

基于时伟等[32]研究冻融循环作用下膨胀土动力学特性的土柱试验进行数学模型验证.土柱直径为70 mm，高140 mm，如图13所示.由于实际模型试验中对土柱的底部和四周采取了保温措施，因此在构建模型时，假定土柱与外界之间的热量交换仅从顶部发生.水分场边界条件设置为零通量，温度场边界条件设置6次单向冻融循环试验，数值模型计算所需参数如表4所示.

图14为土柱顶部位移随时间的变化曲线.由图14可知，随着冻融循环次数的增加，冻融作用对土柱的冻胀融沉特性的影响开始减弱，这与前文中冻融作用对土体的渗透特性与剪切特性参数值的影响规律一致，证明本文所提出的劣化模型的合理性. 同时，考虑了土体结构劣化的模拟值更接近试验值，且模拟值与实测数据误差在5%以内.而未考虑结构损伤模型计算出的最终冻胀位移要比实测值小0.41 mm，误差约为28%.进一步表明了基于土体结构损伤的水-热-力耦合数学模型的可靠性.

2.2桩-土体系数值模拟

2.2.1工程概况

以日喀则萨迦县至秋洛村公路工程曲洛3号中桥桩基为实际工程背景，利用COMSOL有限元数值软件中系数型偏微分方程模块进行二次开发，构建二维桩-土体系数值模型，对桩基服役50年内桩- 土体系位移场进行预测分析.曲洛3号中桥桩长为20 m，桩直径为1.2 m，选用钻孔钢筋灌注桩，混凝土标号为C30，桩周土体以粉质黏土为主.

2.2.2模型建立

基于2.1.3节中建立的考虑土体结构劣化的水- 热-力三场耦合数学模型，利用有限元软件COMSOL系数性偏微分方程进行二次开发，将温度场、水分场以及应力场控制方程转换为系数性偏微分方程形式，如式（17）所示.设置桩周土体材料，将土体剪切特性劣化参数输入.

图15为二维桩-土体系数值模型，设置桩长为20 m，桩径1.2 m，通过试算设定土体的厚度为30 m，土体宽度为12 m（10倍桩径），使用自由三角形划分网格，桩-土界面使用自定义接触对，以固定节点定义接触.桩体及桩周土的相关热力学参数如表5、表6所示.

设置桩-土体系数值模型左右边界为绝热条件，无热传导.上边界设置为当地大气温度，预测50年后平均气温升高1～2.6 ℃[33]，定义温度解析函数；同时，上边界采用狄式边界条件，设置水分场边界为无通量，上边界为自由边界；左边界为对称边界，右边界为辊支撑，限制左右位移，底部为固定约束.

2.2.3预测结果分析

图16为桩周土最大位移随时间变化曲线.分析可知，随着桩基服役年限增加，土体的最大冻胀位移不断减小，冻胀速率变慢.桩基服役前期，随着冻融循环的进行，桩周土体原生结构破坏严重，土体上部含水量增长较快.冻结条件下，液态水的相变导致土体产生冻胀；桩基服役后期，土体结构受冻融循环影响变小，上层含水量增量减小，冻胀位移变化趋于平缓.另外，从图16中可看出，桩周土体的最大融沉量随桩基服役年限的增加逐渐增大，其主要原因为土体结构劣化导致冻融层含水量增加.与桩周土体冻胀位移变化规律类似，冻融循环后期冻融作用对土体融沉位移影响减小.整体来看，融沉位移要远大于冻胀位移，桩周土体呈下沉趋势.该下沉趋势对桩基产生一定的负摩阻力，同时桩-土界面冰膜的融化会进一步降低界面剪切强度.因此，土冻融循环会削弱桩基承载力，对长期服役桩基承载性能带来不利影响.

图17为桩基位移随服役年限的变化曲线.由图17可知，随着桩基服役年限增加，桩基冻胀位移变化趋于平缓，冻胀速率变慢.一方面，主要是因为土体结构劣化受冻融作用影响逐年降低；另一方面，桩基融沉量随服役年限的增加逐渐增大.该现象同样是由于土体结构劣化受冻融作用影响随时间减小，位移变化速率降低.总体而言，桩基服役期间，桩基融沉位移大于其冻胀位移，桩身整体呈融沉特征.对于高原山区实际工程而言，要正确看待桩基融沉现象并采取一定措施以减小其融沉位移的产生.

3结论

本文通过室内试验对高原山区土体渗透特性以及剪切力学特性随冻融循环作用的劣化规律进行了研究，构建了考虑土体劣化的水-热-力耦合数学模型并验证其可靠性.以实际工程为依托，通过数值分析对桩基服役期间位移场进行了预测分析，可得到如下结论：

1）高原山区土体的渗透系数以及孔隙比在冻融作用下呈对数型增大.

2）土体的黏聚力随冻融次数的增加呈指数型减小，内摩擦角随冻融循环次数呈对数型增大.

3）桩基服役过程中，桩周土体最大冻胀位移随服役年限逐渐减小，最大融沉位移随服役年限逐渐增大，土体整体呈融沉趋势，这会削弱桩基的承载能力.

4）桩基服役过程中，桩基础冻胀位移随服役年限的增加逐渐减小，桩基础融沉位移随服役年限逐渐增大，且桩基础融沉位移远远大于其冻胀位移，桩基础整体表现出融沉特征.

参考文献

[1]王万平，张熙胤，陈兴冲，等.考虑冻土效应的桥梁桩-土动力相互作用研究现状与展望[J].冰川冻土，2020，42（4）：12131219.

WANG W P，ZHANG X Y，CHEN X C，et al. Study on dynamic interaction between bridge pile and soil with permafrost effect：status and review[J]. Journal of Glaciology and Geocryology，2020，42（4）：1213-1219. （In Chinese）

[2]ALDAEEF A A，RAYHANI M T. Pile-soil interface characteristics in ice-poor frozen ground under varying exposure temperature [J]. Cold Regions Science and Technology，2021：103377.

[3]熊辉，杨佳.桩-土（软支撑）-结构非线性动力相互作用分析[J].湖南大学学报（自然科学版），2016，43（1）：97-102.

XIONG H，YANG J. Nonlinear dynamic interaction analysis of pile soil （soft support）—structure [J]. Journal of Hunan Univer- sity（Natural Sciences），2016，43（1）：97-102.（In Chinese）

[4]汪優，任加琳，李赛，等.土-结构接触面剪切全过程本构关系研究[J].湖南大学学报（自然科学版），2021，48（3）：144-152.

WANG Y，REN J L，LI S，et al. Study on the constitutive relationship of the whole shear process of soil structure interface[J]. Journal of Hunan University （Natural Sciences），2021，48（3）：144-152.（In Chinese）

[5]QI J L，MA W，SONG C X. Influence of freeze-thaw on engineering properties of a silty soil[J]. Cold Regions Science and Technol- ogy，2008，53（3）：397-404.

[6]REN J，VANAPALLI S K. Effect of freeze-thaw cycling on the soil-freezing characteristic curve of five Canadian soils [J]. Vadose Zone Journal，2020，19（1）：1-8.

[7]陈永，黄英豪，朱洵，等.冻融循环对膨胀土变形和力学特性的影响研究[J].水利水运工程学报，2021（5）：112-119.

CHEN Y，HUANG Y H，ZHU X，et al. Study on the influence of freeze-thaw cycle on the deformation and mechanical properties of expansive soil[J]. Journal of Water Resources and Transport En- gineering，2021（5）：112-119.（In Chinese）

[8]张玲玲，龙建辉，邢鲜丽，等.冻融循环作用下吕梁地区马兰黄土性质研究[J].太原理工大学学报，2021，52（4）：557-563.

ZHANG L L，LONG J H，XING X L，et al.Study on the properties of Malan loess in Liiliang area under freeze-thaw cycles [J]. Journal of Taiyuan University of Technology，2021，52（4）：557- 563.（In Chinese）

[9]CHAMBERLAIN E J，ISKANDAR I，HUNSICKER S E. Effect of freeze-thaw cycles on the permeability and macrostructure of soils [J]. Cold Regions Science and Technology，1990：145-155.

[10]HOTINEANU A，BOUASKER M，ALDAOOD A，et al. Effect of freeze-thaw cycling on the mechanical properties of lime- stabilized expansive clays[J]. Cold Regions Science and Technology，2015，119：151-157.

[11]苏永奇，马巍，钟秀梅，等.冻融循环对青藏粉质黏土动力非线性参数影响的试验研究[J].岩石力学与工程学报，2020，39（S1）：2973-2985.

SU Y Q，MA W，ZHONG X M，et al.Experimental study of influence of freeze-thaw cycles on dynamic nonlinear parameters of Qinghai-Tibet silty clay [J]. Chinese Journal of Rock Mechanics and Engineering，2020，39（S1）：2973-2985.（In Chinese）

[12]赵茜，苏立君，刘华，等.冻融循环对黄土渗透系数各向异性影响的试验研究[J].冰川冻土，2020，42（3）：843-853.

ZHAO Q，SU L J，LIU H，et al. Investigation on the influence of freezing-thawing cycle on the permeability coefficient anisotropy of loess[J].Journal of Glaciology and Geocryology，2020，42（3）：843-853.（In Chinese）

[13]黄旭斌，盛煜，黄龙，等.季节冻土区扩底单桩受力性能研究进展与展望[J].冰川冻土，2020，42（4）：1220-1228.

HUANG X B，SHENG Y，HUANG L，et al.Study of mechanical behaviors of pile foundation with enlarged end in seasonally frozen ground regions：progress and review[J]. Journal of Glaciology and Geocryology，2020，42（4）：1220-1228.（In Chinese）

[14]陳坤，俞祁浩，郭磊，等.多年冻土区灌注桩的人工冷却试验研究[J].中国公路学报，2020，33（9）：104-114.

CHEN K，YU Q H，GUO L，et al. Artificial cooling of cast-in- place piles in permafrost regions[J]. China Journal of Highway and Transport，2020，33（9）：104-114.（In Chinese）

[15]刘庆贺，王永涛，徐湘田，等.冻结粉质黏土-桩基接触面剪切特性试验研究[J].冰川冻土，2020，42（2）：491-498.

LIU Q H，WANG Y T，XU X T，et al. Experimental study on shear characteristics of the interface between the frozen silty clay and the pipe foundation[J]. Journal of Glaciology and Geocryology，2020，42（2）：491-498.（In Chinese）

[16]ZHUSUPBEKOV A，MONTAYEVA A S，MONTAYEV S A，et al. Ensuring chemical resistance of pile foundations when they are installed in permanently and seasonally frozen soils with aggressive environments[J]. IOP Conference Series：Materials Science and Engineering，2020，775（1）：012130.

[17]陈坤，俞祁浩，郭磊，等.基于灌注桩试验的多年冻土区桩-土传热过程分析[J].岩石力学与工程学报，2020，39（7）：1483- 1492.

CHEN K，YU Q H，GUO L，et al.Analysis of pile-soil heat transfer process based on field test in permafrost regions[J]. Chinese Journal of Rock Mechanics and Engineering，2020，39（7）：1483- 1492.（In Chinese）

[18]李仁杰，何菲，王旭，等.不同成樁工艺条件下冻结粉土中基桩承载性状试验研究[J].冰川冻土，2021，43（6）：1809-1817.

LI R J，HE F，WANG X，et al.Experimental research on bearing behavior of piles by using different pile formation technologies in the frozen silt[J]. Journal of Glaciology and Geocryology，2021，43（6）：1809-1817.（In Chinese）

[19]ZATAR W，XIAO F，CHEN G S，et al.Identification of viscoelastic property of pile-soil interactions with fractional derivative model [J]. Journal of Low Frequency Noise，Vibration and Active Control，2021，40（3）：1392-1400.

[20]LAI Y M，ZHANG X F，XIAO J Z，et al. Nonlinear analysis for frost-heaving force of land bridges on Qinghai-Tibet railway incold regions [J]. Journal of Thermal Stresses，2005，28（3）：317-331.

[21]LU J F，YIN J，SHUAIJ. A model for predicting the frost-heave effect of a pile embedded in the frozen soil[J]. Cold Regions Science and Technology，2018，146：214-222.

[22]ALDAEEF A A，RAYHANI M T. Load transfer of pile foundations in frozen and unfrozen soft clay[J].International Journal of Geotechnical Engineering，2020，14（6）：653-664.

[23]ALDAEEF A A，RAYHANI M T. A quick approach for estimating load transfer of conventional and helical piles in ice-rich frozen soils [J]. Geotechnical and Geological Engineering，2021，39（4）：2927-2944.

[24]邢爽，吴桐，李日兵，等.冻土-结构相互作用体系振动台试验及数值分析[J].岩土工程学报，2021，43（11）：2003-2012.

XING S，WU T，LI Y B，et al. Shaking table tests and numerical analysis of frozen soil-structure interaction system[J]. Chinese Journal of Geotechnical Engineering，2021，43（11）：2003-2012. （In Chinese）

[25]唐丽云，杨柳君，贾海梁，等.地下水对多年冻土区桩基承载性能的影响研究[J].铁道标准设计，2020，64（11）：91-98.

TANG L Y，YANG L J，JIA H L，et al. Effects of groundwater on the bearing behavior of pile foundation in permafrost regions[J]. Railway Standard Design，2020，64（11）：91-98.（In Chinese）

[26]土工试验方法标准：GB/T50123—2019[S].北京：中国计划出版社，2019.

Standard for geotechnical testing method：GB/T 50123—2019[S]. Beijing：China Planning Press，2019. （In Chinese）

[27]李栋梁，钟海玲，吴青柏，等.青藏高原地表温度的变化分析[J].高原气象，2005，24（3）：291-298.

LI D L，ZHONG H L，WU Q B，et al. Analyses on changes of surface temperature over Qinghai-Xizang plateau [J]. Plateau Me- teorology，2005，24（3）：291-298 （In Chinese）

[28]钟佳思，童富果，刘畅，等.非饱和黏土受冻融影响的渗透性演化试验研究[J].人民黄河，2018，40（7）：134-137.

ZHONG J S，TONG F G，LIU C，et al. Experimental study on the permeability evolution of unsaturated clay affected by freeze-thaw[J]. Yellow River，2018，40（7）：134-137 （In Chinese）

[29]TAYLOR D W. Fundamentals of soil mechanics [M]. New York：John Wiley and Sons Inc，1948：700.

[30]朱志武，宁建国，马巍.土体冻融过程中水、热、力三场耦合本构问题及数值分析[J].工程力学，2007，24（5）：138-144.

ZHU Z W，NING J G，MA W Constitutive model and numerical analysis for the coupled problem of water，temperature and stress fields in the process of soil freeze-thaw[J]. Engineering Mechanics，2007，24（5）：138-144 （In Chinese）

[31]白青波，李旭，田亚护，等.冻土水热耦合方程及数值模拟研究[J].岩土工程学报，2015，37（S2）：131-136.

BAI Q B，LI X，TIAN Y H，et al Equations and numerical simulation for coupled water and heat transfer in frozen soil[J]. Chinese Journal of Geotechnical Engineering，2015，37（S2）：131-136 （In Chinese）

[32]時伟，张亮，杨忠年，等.冻融循环条件下膨胀土力学特性试验研究[J].西安建筑科技大学学报（自然科学版），2019，51（4）：480-485

SHI W，ZHANG L，YANG Z N，et al Experimental study on mechanical properties of expansive soil of artificial preparation under freeze-thaw cycle conditions [J]. Journal of Xi'an University of Architecture & Technology （Natural Science Edition），2019，51（4）：480-485 （In Chinese）

[33]葛建军.气候变暖对青藏铁路多年冻土区路基影响分析[J].路基工程，2008（3）：6-8.

GE J J Impact of climate warming on subgrade in permafrost area of Qinghai-Tibet railway[J]. Subgrade Engineering，2008（3）：6- 8 （In Chinese）