渗水作用下膨胀土中单桩承载变形特性研究 渗水作用下膨胀土中单桩承载变形特性研究
刘建新1 邓洪权1 郭志广2 刘运思2
(1.中铁二十局集团有限公司 陕西西安 710016；2.湖南科技大学岩土工程稳定控制与健康监测省重点实验室 湖南湘潭 411201)
摘 要：目前膨胀土中桩基设计仍基于非膨胀土力学的设计原则。基于FLAC3D采用热-力耦合方法实现吸热膨胀来模拟膨胀土的吸湿膨胀，并重点阐述温度场模拟湿度场的关键参数取值方法。建立膨胀土中单桩三维实体模型，得出渗水作用下膨胀土中单桩的荷载变形特性、桩侧摩阻发挥特征，并分析了桩长、桩径、膨胀系数等对单桩承载特性的影响规律。研究结论：非扩底桩最小桩长宜大于2倍膨胀影响深度；细长桩比短粗桩更能有效降低膨胀土中桩顶的位移。该方法可获得渗水作用下任一时刻的桩身受力变形特性，为研究膨胀土中桩基受荷性能提供一种可行的变通手段。
关键词：桩土相互作用 膨胀土地基 承载变形特性 热力耦合
1 引言
膨胀土中黏粒主要为蒙脱石和伊利石，具有吸水膨胀、失水收缩且反复胀缩的特性。渗水作用下，极易导致上部或邻近建筑物损坏或事故。据统计，位于膨胀土区的常规铁路，路基完好率仅为25%[1]。为保证膨胀土区构筑物安全，常用桩进行加固。随着高速铁路的快速发展，通过膨胀土区已无法避免。但高铁无砟轨道路基“毫米级”的工后沉降要求，使得CFG桩[2]、钻孔灌注桩[3]、微型桩[4]等大量使用。
4）病虫防治。清除梨园枯枝落叶，剪除病虫枝，刮除老翘皮，并将其集中烧毁或深埋；入冬前，用谷草做1～2个小草把，诱引害虫在草把上越冬。开春时将草把集中烧毁。落叶后至发芽前，梨园喷2～3波美度石硫合剂1～2次，药渣可用于刷白树干。
现代青花艺术虽与传统青花瓷有明显的差异，但是在创作表达上依旧存在关联性。虽然说现代陶艺家创作时会因为生活背景、民族信仰等一系列因素受到影响，但同时也会对传统青花产生评判比较，其实这正是对传统青花艺术的继承创新。其实不仅是现代艺术家，哪怕是在过去，在技术、题材上后代工匠也是不停的对前人的创造进行改善。所以现代陶艺家在进行创作时，对传统青花的了解是必不可少的，这是推动现代青花发展的必要步骤之一。
桩与膨胀土相互作用机理迥异于常规桩土作用。目前，膨胀土中桩基受力机理研究还不够全面，研究方法主要有模型试验[5]、现场测试[6]和理论计算[7]等。而数值模拟方面，膨胀土因增湿导致膨胀的变化规律过于复杂，且膨胀土属于典型的非饱和土，开发膨胀土本构模型涉及参数多且工程应用难度较大。基于温度场与湿度场的相似性，缪协兴[8]等提出了渗水作用下膨胀岩体的湿度应力场理论。此后该方法被广泛用于计算分析膨胀区的隧道[9]、边坡[10]和基坑[11]受力特性。
本文首先介绍FLAC3D中采用热-力耦合模拟湿度应力场的方法；其次通过对比分析已有室内模型试验结果，对所提方法进行验证；最后进行参数敏感性分析，进一步探讨膨胀土中单桩的受力性状和桩土作用机理。
2 基于热力耦合的吸湿膨胀特性模拟方法
热传导与渗流问题在数学描述和定解条件以及各相关变量的物理含义等方面具有相似性[12]。采用热力耦合来模拟膨胀土的吸湿膨胀，需要对等吸热与吸湿产生的膨胀变形，以及对等升温与吸湿的速率。
式中，α为热膨胀系数(1/℃)；β为膨胀土的线膨胀系数；Δw为含水率变化量；ΔT为温度变化量；δ为Kronecker记号。
教师通过让学生熟读课文，之后让学生用自己的话来理清故事的发展脉络，有利于学生更好地把握文章的整体方向和描写的人物。《包身工》这篇文章主要描写的就是包身工的居住环境、饮食环境、工作环境等，通过对这些环境的描写来揭示包身工的命运。同时，文章的描写也体现了作者对于包身工满满的同情，以及作者对于当时制度的不满。
（二）默读能保持阅读情绪，有利于静心思考，能体现阅读的个性化。默读由于不出声，不会互相干扰，保证了环境的安静，更利于集中思考、理解所读的内容，而且不易疲劳，易于持久。平时上课时，我经常发现学生刚上课时均能精神饱满地大声朗读，可过一会儿后就会发现一部分同学变朗读为默看了，究其原因是疲劳造成的，可见朗读很难保持积极有效的阅读情绪。朗读时，人的视觉、听觉、发音系统都要高度集中，时间一久，自然就会疲劳；而默读则很好地弥补了这一缺陷，它能让学生以饱满的热情和精力，把通过眼睛获得的文字信息迅速处理并静心思考，能以个体为主，让学生真正做到独立感悟，且不会受到他人的影响。
令两者应变增量相等，则:
式(3)便为热胀模拟湿胀的对应关系。膨胀土浸润速度决定吸湿膨胀的速度，模拟时可由导热系数λ控制。参数转化后，可借助FLAC3D热-力耦合分析来模拟膨胀土含水率变化时的力学性能。
3 膨胀土中单桩模型案例分析
基于FLAC3D的热力耦合方法模拟膨胀土中单桩浸水模型试验，并与试验结果对比以验证其有效性。
3.1 室内模型试验概况
试验装置[5]为一直径50 cm、高90 cm的油桶，自底向上依次填筑0.1 m卵石、0.16 m细砂及0.58 m膨胀土(按最优含水率15%预配)。模型桩为φ50 PVC管，桩长65 cm(土中58 cm)。管壁上粘贴应变片量测桩身变形，桩端受力通过压力盒测量，桩顶位移和地面隆起通过位移计测量。
课程体系设置是高职院校人才培养目标的重要体现，高职院校通常以专业能力和就业导向为主，力图使旅游专业课程设置涵盖旅游相关的所有岗位类型，这就相对忽视了学生综合素质与创新能力的培养。所以，在旅游新常态背景下，高职院校旅游专业课程体系设置应致力于提升学生的综合素质和创新能力。
人工渗水，膨胀土浸水历时230 h直到桩顶位移和地面隆起趋于稳定。此时地面隆起41.2 mm，桩顶上升3.59 mm。
3.2 膨胀土单桩模型建立与结果分析
(1)有限差分网格及本构
依据模型试验，建立长0.5 m、宽0.25 m、高0.75 m的有限差分模型，桩和土体均采用三维实体单元，节点数4 890个、单元数4 056个。
桩和底部细砂采用弹性模型。膨胀土破坏符合摩尔库伦准则，采用各向同性热传导模型，其余设为非传热模型。膨胀土与桩体设置接触面单元，其剪切特性采用折线型，表达式见式(4)。
取消药品加成前后住院病人性别、年龄构成、医保身份构成，平均住院天数无差异。 26 550人次 ，男性10 832人次，女性15 718人次；平均年龄为56.25岁；医保身份城镇职工20 433人次，城镇居民6 117人次；平均住院天数10.10天。 27 304人次，男性10 796人次，女性16 508人次；平均年龄为56.45；医保身份城镇职工21 110人次，城镇居民 6 194人次；平均住院天数10.07天。
式中，τs为峰值剪应力；ws为相对τs时的极限位移(可取5～25 mm)；τ1为残余应力；ks为剪切刚度系数。峰值剪应力τs和残余剪应力τ1可由式(5)求得，而剪切刚度ks依据FLAC3D手册确定。
式中，cs和φs分别为黏聚力和摩擦角；c1和φ1分别为残余黏聚力和残余摩擦角。
模拟计算采用的相关参数见表1～表2。
表1 浸水试验材料物理力学参数
类别 密度/(g·cm-3)弹模/MPa泊松比黏聚力/kPa内摩擦角/(°)膨胀土 1.63 6.0 0.3 50 40细砂 1.80 50.0 0.3 - -桩体 1.50 3 000 0.2 - -
表2 桩与膨胀土接触参数
剪切刚度/(GPa·m-1)极限黏聚力/kPa极限内摩擦角/(°)残余黏聚力/kPa残余内摩擦角/(°)1.0 46.7 35.4 7.5 27.1
(2)边界条件及模拟过程
模型底部及侧面约束法向位移，顶面为自由面。热传导模拟渗水过程中，膨胀土初始温度设为0℃(对应初始含水率15%)，然后固定膨胀土顶面温度为100℃(对应饱和含水率26%)。经多次试算，热膨胀系数取值为1.68×10-4，导热系数取值为1.0×105。
在前处理中为压电陶瓷两极加载瞬时电压激励，每片压电陶瓷激励电压幅值为100 V。采用直接法在18 kHz～21 kHz范围内取30个点对整个模型进行积分运算，求解换能器稳态动力学响应。直接法不需要提取系统的特征模态，结果更为精确。图4为换能器频率响应曲线，从图中可以看到，在频率为18 761 Hz时换能器末端输出振幅达到了43 μm，满足超声加工的需要。
(3)模拟结果及对比分析
膨胀土和桩身位移云图见图1，竖向应力云图见图2。
图1 渗水作用下位移云图
图2 渗水作用下竖向应力云图
由图1可知，膨胀土顶面最终隆起量模拟值为4.05 cm(试验值为3.95 cm)，桩顶最终隆起量模拟值为4.01 mm(试验值为5.71 mm)，与实测结果相近。
由图2可知，由于桩顶无荷载，整个桩身均处于受拉状态，桩身轴向应力最大值(530 kPa)出现在桩身偏上位置。归一化桩身轴力模拟值与试验值(见图3)沿深度变化趋势一致，但实测的峰值点偏向桩的中下部，原因在于膨胀土样经击实后土层参数沿埋深不同，而模拟时视为均匀。
图4为计算结束时桩侧摩阻力分布曲线，可以看出膨胀将引起桩身向上隆起，而下部土体则起锚固作用，即呈现上拔下锚的状态，且桩侧摩阻为零的点处(中性点)，桩身轴力最大。
图3 桩身轴力分布对比
图4 桩侧摩阻分布曲线
4 膨胀土单桩荷载传递特性影响因素分析
从桩长、桩径、膨胀系数以及桩顶荷载方面探讨增湿膨胀对桩荷载传递的影响，主要侧重于膨胀对桩身位移和轴力的影响。膨胀土和桩等物理力学参数见表1。参数基准值为:桩长10 m，桩径0.2 m；桩顶无荷载作用；膨胀土影响深度均为4 m。
和通话者确定地点与患者症状表现，叮嘱先不要搬动患者，以平卧的方法保持其头部朝向一侧，同时立即出车，做好车上急救准备。赶赴到现场后，对患者病情做出初步评估与诊断，检查其脉搏、血压、心电图，建立静脉通道、吸氧，做好运送至医院的各种准备，取得家属的同意与支持。运送过程中除了要确保患者呼吸道顺畅外，还应和院内急诊科取得联系，告知详情，使院方能够做好接诊准备。
4.1 桩长影响分析
保持其它参数不变，桩长分别取4 m、6 m、8 m和10 m，即长径比分别为20∶1、30∶1、40∶1和50∶1。 计算得出的桩身位移和轴力沿埋深分布见图5～图6。
协同工作模式能充分减少不必要的时间消耗，使得工程师可以将更多的时间用于审核，从而提高工作的总体质量；通过高效的沟通和对接，还可促进高质量优化意见的产生。
图5 桩身位移随桩长变化曲线
图6 桩身轴力随桩长变化曲线
由图5可知，对于4 m桩，桩端因膨胀而上移(39.6 mm)，此时桩身位移最大(桩顶43.9 mm)，工程中应避免该工况。随桩长增加，桩身位移明显减小，当桩基嵌入稳定土层长度等于影响深度后再增加桩长并未能减少桩身位移。表明通过桩基嵌入非膨胀土层来抑制膨胀土的膨胀危害是可行的，且桩基嵌入非膨胀土层要足够深，但并不是越深越好，设计中存在临界值，本实例中临界值约为影响深度(8 m和10 m桩长桩身位移基本一致)，而Polous给出的建议值为2倍的影响深度值。
图6表明，影响深度中间部位(2 m)桩身轴力最大，沿上下方向轴力逐渐减小；随着桩长增加，桩端嵌入稳定土层，轴力最大值逐渐增加，至桩端嵌入深度等于影响深度，最大轴力保持不变。
4.2 桩径影响分析
保持其它参数不变，桩径分别为0.2 m、0.3 m、0.4 m、0.5 m 和 1 m，即桩径分别为 0.02L、0.03L、0.04L、0.05L和0.1L。计算得出桩身位移和轴力沿深度分布如图7～图8所示。
一日，县城“百里香肠铺”掌柜百里香来找坤二少爷。百里香本名田桂贤，也是油铺人，家在村东头，跟坤二少爷既是同年，又是同学，二人自幼交好，情同手足。
图7 桩身位移与桩径关系曲线
图8 桩身轴力与桩径关系曲线
由图7可知，不同桩径下的桩身位移均沿深度递减。随桩径增加，存在一个转折点(约2.0 m处)，使得不同桩径下的桩身位移在该点上下变化规律不同。桩径越大，该点以上桩身位移先减小后逐渐增大，而该点以下桩身位移逐渐增大。即桩径越大，膨胀变形影响深度越大。相反，大直径桩对减小桩顶抬升量效果不明显。如1 m桩径时只比0.2 m桩径时桩顶抬升量减少12%，而桩端位移反而呈百倍增加。图8表明:随桩径增加，桩身轴力增大，且桩身中下部轴力增加显着。可见，嵌入稳定土层较深的小直径桩能有效地降低膨胀土中桩顶位移。
4.3 膨胀系数影响分析
保持其它参数不变，影响深度范围内膨胀系数分别为1α(α ＝1.68 ×10-4)、2α、3α、4α 和 5α。 计算得出桩身位移和轴力沿深度分布见图9～图10。
图9 桩身位移与膨胀量关系曲线
图10 桩身轴力与膨胀量关系曲线
由图9可知，随着影响深度范围内膨胀量的增加，桩顶隆起逐渐增大，同时膨胀影响深度也逐渐加大。如膨胀系数为4α和5α时，桩端位置上升0.6 mm和2.0 mm。
由图10可知，随着影响深度范围内膨胀量的增加，桩身最大轴力(约2.5 m处)先增大，至膨胀系数达到3α后，最大轴力开始小幅减小且位置下移，而且影响深度逐渐加深。桩身最大轴力处至桩顶部位，由于桩顶部位土体和桩身位移有较大隆起量，达到接触面峰值强度后，桩身轴力具有减小的趋势。
5 结论
基于FlAC3D的热力耦合方法探讨了渗水作用下膨胀土中单桩的荷载传递特性，并进行影响因素分析。主要得到以下结论:
(1)基于温度场和湿度场相似的理论和数学描述，采用吸热膨胀来模拟膨胀土的吸湿膨胀，经模型试验对比分析，证明该方法有效可行。其中热膨胀系数α和导热系数λ是关键参数。
(2)桩长增加，桩身抬升位移减小。增加桩的嵌入深度可降低隆起变形，提高上部结构物的安全性。但桩长增加，会引起桩身拉力的增加，严重时会将桩拉断，因此，必须进行抗拉验算，必要时增加桩身配筋。建议一般土层中的非扩底桩，最小桩长不宜小于2倍的膨胀影响深度。
(3)埋入稳定土层的小直径桩能有效降低膨胀土中桩顶位移，而大直径桩对减小桩顶位移效果不大。
(4)随影响深度内膨胀量的增加，桩身位移和轴力不断增大，必须相应地加大桩长和增加桩身配筋来增强桩的抗拉性能，或者采取措施(比如隔离或对膨胀土进行改良)减弱膨胀的影响。
参考文献
[1] 周海.膨胀土(岩)地区常见铁路路基病害及整治研究[J].路基工程，(1):44-48.
[2] 冯研.高速铁路路基下膨胀土荷载传递规律与变形特性研究[D].成都:西南交通大学，.
[3] 王诗.基于正交试验方法优化膨胀土地区桩板结构路基设计参数的研究[D].兰州:兰州交通大学，.
[4] 陈伟志，李安洪，蒋关鲁，等.铁路路基下膨胀土长短微型桩抗隆起研究[J].岩石力学与工程学报，，38(2):409-423.
[5] 范臻辉.膨胀土地基胀缩特性及桩-土相互作用研究[D].长沙:中南大学，.
[6] 吴志伟，宋汉周，王宏宇，等.膨胀土中基桩胀拔力原型试验研究[J].岩土力学，，33(S2):173-177.
[7] Liu Yun-long，Vanapalli S.K.Load displacement analysis of a single pile in an unsaturated expansive soil[J].Computers and Geotechnics， ，106:83-98.
[8] 缪协兴，杨成永，陈至达.膨胀岩体中的湿度应力场理论[J].岩土力学，1993(4):49-55.
[9] 曾仲毅，徐帮树，胡世权，等.增湿条件下膨胀土隧道衬砌破坏数值分析[J].岩土力学，(3):871-880.
[10]郑俊杰，吕思祺，曹文昭，等.高填方膨胀土作用下刚柔复合桩基挡墙结构数值模拟[J].岩土力学，(1):395-402.
[11]李朝阳，谢强，康景文，等.基于摩尔-库仑准则的膨胀土弹塑性本构模型及其数值实现[J].土木建筑与环境工程，(2):92-99.
[12]范从友.合肥膨胀土吸湿膨胀特性数值模拟研究[J].土工基础，，31(6):753-757.
Research on Bearing and Deformation Characteristics of a Single Pile in Expansive Soil Under the Action of Seepage
LIU Jianxin1， DENG Hongquan1， GUO Zhiguang2， LIU Yunsi2
(1.China Railway 20thBureau Group Co.Ltd.， Xi’an Shaanxi 710016， China； 2.Hunan Provincial Key Laboratory of Geotechnical Engineering for Stability Control and Health Monitoring， Hunan University of Science and Technology， Xiangtan Hunan 411201， China)
Abstract:The existing problem is that the design principle of pile foundation in expansive soil foundation still follows the design principle of non-expansive soil mechanics.Based on FLAC3D，the thermal-mechanical coupling method is used to simulate the hygroscopic expansion deformation of expansive soils.The method of key parameters of temperature field simulation humidity field is given.A three-dimensional solid model of single pile in expansive soil is established.The load-deformation characteristics and the frictional resistance of single pile in expansive soil under seepage are obtained.The influence of pile length， pile diameter， expansion coefficient on the bearing characteristics of single pile is analyzed.The research conclusions are as follows:it is suggested that the minimum length of non-enlarged bottom pile should not be less than 2 times the depth of expansion effect；the displacement of pile top in expansive soil can be reduced more effectively by the long pile with small diameter than by the short pile with large diameter.The method can simulate the displacement and load transfer characteristics of piles at any time under the influence of seepage，and provides a feasible means for studying the load-bearing behavior of piles in expansive soil.
Key words:pile-soil interaction； expansive soil foundation； load-bearing and deformation characteristics； thermalmechanical coupling
中图分类号：TU475＋.5；TU473.1
文献标识码：A
DOI:10.3969/j.issn.1009-4539..04.003
收稿日期：-02-23
基金项目：湖南科技大学岩土工程稳定控制与健康监测省重点实验室开放基金资助项目(E21708)
作者简介：刘建新(1981-)，男，内蒙古商都人，高级工程师，主要从事隧道与地下工程施工及技术管理；E-mail:413907782@
