Experimental study on influences of wetting-drying cycles on microstructure and water-retention characteristics of clay
-
摘要: 长期干湿循环作用会削弱路基土的工程力学性质，产生工程灾害，有必要研究干湿循环对其微观结构及持水性能的影响。针对既有线有砟轨道路基土，开展了一系列考虑干湿循环影响的压汞试验、扫描电镜试验、压力板试验与饱和盐溶液蒸汽平衡试验。试验结果表明：在反复的干湿循环过程中，土体结构遭到破坏，内部胶结物质不断减少。土体微孔及小孔（< 5 μm）数量逐渐增多，中孔及部分大孔范围（5~30 μm）内的峰值孔径逐渐减小，分布密度逐渐增大。数次干湿循环后微裂隙（> 100 μm）逐渐增多，导致试样开裂。干湿循环会对土体的持水性能产生影响，初始干密度较小时，经历过1次干湿循环持水性能略有增强，而经历过3次及5次干湿循环持水性能持续减弱。Abstract: The Long-term wetting-drying cycles can weaken the engineering characteristics of subgrade fillers, resulting in engineering disasters. It is necessary to study the effects of the wetting-drying cycles on their microstructure and water-retention performance. A series of mercury intrusion tests, scanning electron microscopy tests, pressure plate tests and saturated salt solution vapor equilibrium tests considering the influences of the wetting-drying cycles are carried out on the subgrade filler of the existing ballasted track roads. The experimental results show that during repeated wetting-drying cycles, the soil structure is damaged, and the internal cementitious materials are continuously reduced. The number of micropores and small pores (< 5 μm) in the soil gradually increases, while the peak pore size in the range of mesopores and some large pores (5~30 μm) gradually decreases, and the distribution density gradually increases. After several wetting-drying cycles, the number of microcracks (> 100 μm) gradually increases, leading to cracking of the soil specimens. The wetting-drying cycles have effects on the water-retention characteristics of the soil. When the initial dry density is low, the air-entry value after experiencing one wetting-drying cycle is slightly higher than that without wetting-drying cycles, while the air-entry value of the specimen after three and five wetting-drying cycles continues to decrease.
-
0. 引言
中国大规模的基础设施建设工程每年产生大量的建筑垃圾，建筑垃圾可通过破碎、清洗、筛分、分级等二次处理后生产出再生骨料[1-2]。再生骨料可用于道路铺筑、工程回填、地基加固等岩土工程领域[3-4]。
强夯碎石桩因其具有提高地基承载力、减小地基沉降和加速地基固结排水等优点常应用于软土地基加固中，施工工序如图 1所示：首先将钢套筒打至设计深度并筒内清泥，套筒内分层填入碎石并使用夯锤夯实，上拔套筒一定高度并再次填料和夯实，直到形成串珠状的散体桩复合地基[5]。许多学者对碎石桩的承载性能、荷载传递、破坏模式等进行了研究[6-8]。受桩长、桩周土体、加载方式的影响，碎石桩可能出现的破坏有3种：鼓胀破坏、刺入破坏和剪切破坏。基于有限单元法和有限差分法等已广泛应用于碎石桩宏观特性的研究[9-10]，然而碎石桩由散体材料组成，连续方法不能准确地模拟其特性。
通过建立三维离散-连续耦合数值模型对再生骨料强夯桩的承载性能进行模拟，研究了再生骨料强夯桩的承载变形机理、桩体破坏模式及荷载传递规律。并分析了不同桩长和不同孔隙率的再生骨料强夯桩的承载性能。以期能为再生骨料强夯桩的推广应用提供参考。
1. 离散-连续耦合数值模型
1.1 数值模型建立
再生骨料强夯桩按图 1所示方法进行成桩，具体细节详见文献[11]。试验选用的再生骨料相对质量密度为2.62，内摩擦角为39度，粒径范围为0.6~9.5 mm，桩周土体为天津路基地表土，相对质量密度为2.65，厚度为60 cm，含水率为15%，土体黏聚力为20 kPa，内摩擦角为32度。成桩后桩长为60 cm，桩径为6 cm。成桩完成后进行加载试验，采用位移加载的方式施加荷载，加载速度为1 mm/min，加载至5 cm时停止加载。
采用FLAC-PFC耦合方法进行承载模拟，如图 2（a）所示，耦合数值模型包括3部分：有限差分模型（桩周土）、离散元模型（散体桩）和耦合墙。模型的建立步骤与室内模型试验一致，分为土基填筑、成桩和加载三步。首先建立桩周土的有限差分模型，采用的本构模型为莫尔-库仑，所用参数见表 1，接下来先把散体桩区域进行挖空，然后在桩底、桩顶和桩侧生成墙形成封闭区域使得颗粒在此区域内生成，桩底墙直径为25 cm，桩顶墙直径为6 cm，桩侧耦合墙直径为6 cm、高为60 cm，颗粒之间采用抗转动线性接触模型（rrlinear model），颗粒细观参数通过直剪试验进行标定后列于表 2。
表 1 土体参数Table 1. Parameters of soilsc/kPa ϕ/(°) E/MPa ν ρ/(g⋅cm-3) 20 32 10 0.4 1830 表 2 再生骨料标定参数Table 2. Calibration parameters of recycled aggregates接触类型 Ball-ball Ball-wall 有效模量E*/Pa 2.0×107 5.0×107 刚度比k* 1.5 1.5 抗转动摩擦系数μr 0.5 — 摩擦系数μ 0.5 0.5 阻尼比 0.7 0.7 孔隙比 0.86 — 最后在再生骨料强夯桩上方生成与桩径相同的加载板，以1 mm/min的恒定速率对碎石桩进行轴向加载，直至位移达到5 cm。如图 2（b）所示，在桩内布置19个测量圆，球体中心垂直位置位于-3~-57 cm处（土表z=0），获得测量圆内应力状态（轴向和径向应力）和微观特征（孔隙率和配位数）。在桩周土有限差分模型中布置21个监测节点，获取连续-离散界面节点处的侧向变形。
1.2 数值模型验证
图 3（a）为模型试验与数值模拟的荷载-沉降曲线对比图，两者曲线变化趋势相同，吻合程度较好。从图 3可以看出，在加载初期沉降达到5 mm之前曲线基本呈线性增加，桩顶应力在较小的沉降量下即迅速增加到约750 kPa；当沉降超过10 mm时，非线性变形明显发展；当桩顶应力超过1500 kPa后沉降迅速发展，沉降超过20 mm后p-s曲线几乎竖直，预示着再生骨料强夯桩丧失了承载能力。
图 3（b）为模型试验与数值模拟的桩体径向变形对比图，从图中可以看出，两者曲线整体变化趋势相同，桩体径向变形都主要发生在桩身20 cm（约3倍桩径）及以上，桩身20 cm以下变形较小，最大变形量均为9 mm左右，桩体呈现上部凸出中下部小的形态。
图 3（c）为模型试验与数值模拟的桩身应力传递系数对比图，两条曲线均为桩体最终破坏时测得，从图中可以看出随埋深增加应力传递系数逐渐减小，超过3倍桩径后传递系数均小于0.2，表明大部分荷载均由碎石桩上部承担，模拟和试验具有较好的一致性。数值模拟中上部第一个和第二个数据监测点异常减小是因为此时颗粒大部分被加载板下压，在测量球中的颗粒较少，所以其平均应力较小使得传递系数较小。
2. 结果分析
2.1 桩周土变形和应力
图 4给出了不同桩顶沉降下（18，34，50 mm）桩周土体的变形和应力云图。从图 4（a）可以看出桩体压缩主要产生在桩体上部，随着沉降增加该段区域内碎石产生剧烈鼓胀变形，从而引起桩周土体产生水平位移。土体水平应力在桩体鼓胀导致的水平挤压下增长剧烈。
![]()
图 4 不同沉降下桩周土变形和应力Figure 4. Deformations and stresses of soil around column under different settlements图 4（b）显示桩周上部土体附近产生了较大的竖直应力增量，增加的竖向应力主要由桩侧摩阻力及变形后的桩体与土体挤压产生，并通过土体剪应力以一定扩散角度传递到远离桩体的土体区域。图 4（c）显示桩周土体在靠近桩体压缩段也出现了较大的竖向位移，但位移量明显小于碎石沉降量。因此可以判断在桩体压缩强烈的范围内，桩体和桩周土体产生了较大相对位移。加载结束后土地表面产生了明显的隆起。
图 4（d）为不同沉降下桩周土体剪切应变分布，土体剪切应变大小直接反映桩体变形对土体的扰动强弱。随着桩体沉降量的增加，土体剪切扰动区域在桩体水平鼓胀推动下向四周扩展，最终观察到贯通地表的连续剪切滑动面。
2.2 桩体变形和应力
图 5（a）为不同沉降时桩体鼓胀变形沿桩身分布图，桩体鼓胀变形整体呈现上部突出中下部小的形态。随着沉降增加鼓胀变形逐渐增加，主要鼓胀变形发生在桩体3倍桩径的埋深范围内，径向应变最大为14%，最大鼓胀变形位置大约在1.5倍桩径处。图 5（b）为桩体水平应力沿桩身分布图，从图中可以看出水平应力随沉降增大而增大，随埋深增大发生较大衰减。与径向应变变化规律一致。
![]()
图 5 不同沉降下桩体变形和应力Figure 5. Column deformations and stresses under different settlements图 6显示了再生骨料颗粒在不同沉降下接触力的分布。随沉降增大，桩顶处直接承受荷载的强接触力链越来越明显，且其接触力值越来越大。接触力链主要集中在桩体的上部，大概3倍桩径深度内，这与径向应变一致。在加载过程中，强接触力链承担桩顶荷载后再通过若干弱接触力链向下部及桩周土体传递荷载。
![]()
图 6 不同沉降下颗粒间接触力链Figure 6. Contact force chains between particles under different settlements2.3 破坏模式
再生骨料强夯桩的破坏模式示意图如图 7所示。随着荷载的增加，再生骨料强夯桩最大鼓胀变形（1.5倍孔径）附近的土体首先进入屈服状态，接下来的荷载会使石柱产生较大的鼓胀变形并挤压石柱周围土体，造成土体塑性区以半球形式向周边扩展。随着荷载继续增加，桩周土体在桩体鼓胀变形的挤压下形成贯通地表的连续滑裂面，地表发生隆起变形，影响范围约为18 cm（3倍桩径）。桩体则最终因为鼓胀段（3倍桩径）径向应变无法限制而丧失承载力。
3. 参数分析
3.1 桩体长度的影响
其他条件不变，只改变桩体长度，得到不同长度下再生骨料强夯桩的荷载-沉降曲线如图 8所示。可以看出桩长过小时，再生骨料强夯桩的承载力很低，只有540 kPa，在桩顶沉降20 mm时便发生破坏，破坏形式为刺入破坏。随着桩长的增加，承载力逐渐提高。桩长7d时与全长桩（10d）时承载力接近，这说明全长桩底部3倍桩径并不承担荷载，有效长度为7d。从5d开始，再生骨料强夯桩发生鼓胀破坏。
3.2 孔隙率的影响
其他条件不变（桩长为10d），只改变孔隙率，得到不同孔隙率下再生骨料强夯桩的荷载-沉降曲线如图 9所示。
随着孔隙率减小，再生骨料强夯桩的承载性能逐渐提高。在孔隙率为0.3，0.35时，曲线随着沉降增加先为线性增加，然后快速发展，最后曲线趋于竖直，而孔隙率为0.4，0.45，0.5时，曲线近似线性变化，原因是桩体孔隙率较大时，颗粒间咬合力较弱，颗粒向桩周土位移时，桩周土能提供较强的侧限作用，当增加到一定沉降后，桩体被压密，此时颗粒不断挤向土体，因侧限作用失效产生破坏。
4. 结论
通过建立三维离散-连续耦合数值模型对再生骨料强夯桩的承载性能进行模拟，研究了再生骨料强夯桩的的承载变形机理、桩体破坏模式及荷载传递规律。并分析了不同桩长和不同孔隙率的再生骨料强夯桩的
承载性能。主要得到以下3点结论。
（1）桩周土体受颗粒挤压变形导致桩体上部位置位移和应力发展迅速。在水平方向表现为鼓胀变形，在竖直方向上地表会有隆起和抬升。鼓胀变形主要集中在3倍桩径范围内。
（2）竖向荷载使桩体产生较大的鼓胀变形并挤压石柱周围土体，造成土体塑性区以半球形式向周边扩展。随着荷载继续增加，桩周土体在桩体鼓胀变形的挤压下形成贯通地表的连续剪切滑裂面。
（3）桩长小于3d，桩体发生刺入破坏，大于3d时发生鼓胀破坏。全长桩底部3倍桩径并不承担荷载，有效长度为7d。随着孔隙率减小，再生骨料强夯桩的承载性能逐渐提高。
-
![]()
图 2 压力板法和饱和盐溶液蒸汽平衡法
Figure 2. Pressure plate method and saturated salt solution vapor equilibrium method
![]()
图 3 1~5次干化终点试样外形
Figure 3. Appearances of end point specimens after drying cycles of 1~5
![]()
图 4 饱和状态及干化终点试样孔隙比变化
Figure 4. Changes in porosity of specimens in saturated and drying states
![]()
图 5 干湿循环对土体孔隙结构的影响（压汞试验）
Figure 5. Influences of wetting-drying cycles on pore structure of soil mass (mercury intrusion tests)
![]()
图 6 土体微观结构随干湿循环的变化
Figure 6. Evolution of soil microstructure with wetting-drying cycles
![]()
图 7 初始干密度对SWCC的影响
Figure 7. Influences of initial dry density on soil-water characteristic curve (pressure plate method)
![]()
图 8 低吸力范围内干湿循环对SWCC的影响
Figure 8. Influences of wetting-drying cycles on soil-water characteristic curve in low suction range
表 1 试验黏土的基本物理参数
Table 1 Basic physical parameters of clay used in tests
GS 液限wL/% 塑限wP/% 塑性指数IP 最优含水率wop/% 最大干密度γd，max/(g·cm-3) 2.73 46.1 22.4 23.7 20.4 1.635 
下载: 导出CSV
-
[1] ZHAO G T, ZOU W L, HAN Z, et al. Evolution of soil-water and shrinkage characteristics of an expansive clay during freeze-thaw and drying-wetting cycles[J]. Cold Regions Science and Technology, 2021, 186: 103275. doi: 10.1016/j.coldregions.2021.103275
[2] DELAGE P, LEFEBVRE G. Study of the structure of a sensitive Champlain clay and of its evolution during consolidation[J]. Canadian Geotechnical Journal, 1984, 21(1): 21-35. doi: 10.1139/t84-003
[3] CUI Y J, LOISEAU C, DELAGE P. Microstructure changes of a confined swelling soil due to suction controlled hydration[C]// Proc of the Third Inter Conf on Unsaturated Soils, UNSAT, 2002.
[4] 叶为民, 钱丽鑫, 陈宝, 等. 高压实高庙子膨润土的微观结构特征[J]. 同济大学学报(自然科学版), 2009, 37(1): 31-35. YE Weimin, QIAN Lixin, CHEN Bao, et al. Characteristics of micro-structure of densely compacted gaomiaozi bentonite[J]. Journal of Tongji University (Natural Science), 2009, 37(1): 31-35. (in Chinese)
[5] 赵天宇, 王锦芳. 考虑密度与干湿循环影响的黄土土水特征曲线[J]. 中南大学学报(自然科学版), 2012, 43(6): 2445-2453. ZHAO Tianyu, WANG Jinfang. Soil-water characteristic curve for unsaturated loess soil considering density and wetting-drying cycle effects[J]. Journal of Central South University(Science and Technology), 2012, 43(6): 2445-2453. (in Chinese)
[6] KONG L W, SAYEM H M, TIAN H H. Influence of drying– wetting cycles on soil-water characteristic curve of undisturbed granite residual soils and microstructure mechanism by nuclear magnetic resonance(NMR) spin-spin relaxation time (T2) relaxometry[J]. Canadian Geotechnical Journal, 2018, 55(2): 208-216. doi: 10.1139/cgj-2016-0614
[7] 张俊然, 许强, 孙德安. 多次干湿循环后土-水特征曲线的模拟[J]. 岩土力学, 2014, 35(3): 689-695. ZHANG Junran, XU Qiang, SUN De'an. Simulation of soil-water characteristic curves during drying and wetting cycles[J]. Rock and Soil Mechanics, 2014, 35(3): 689-695. (in Chinese)
[8] YI F C, WANG Z, ZHOU M J. Research on the soil water characteristic curve and prediction of compacted bentonite[J]. Advanced Materials Research, 2012, 518/519/520/521/522/523: 2785-2791.
[9] 蔡国庆, 刘祎, 徐润泽, 等. 全吸力范围红黏土干湿循环土–水特征曲线[J]. 岩土工程学报, 2019, 41(增刊2): 13-16. doi: 10.11779/CJGE2019S2004 CAI Guoqing, LIU Yi, XU Runze, et al. Experimental investigation for soil-water characteristic curve of red clay in full suction range[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(S2): 13-16. (in Chinese) doi: 10.11779/CJGE2019S2004
计量
- 文章访问数: 107
- HTML全文浏览量: 17
- PDF下载量: 31
