Dynamic properties of undisturbed sand in deep overburden of a dam foundation
-
摘要: 超深厚覆盖层地基是当前水利水电工程建设所面临的巨大挑战, 覆盖层土体的动力力学特性是进行地基土液化判别和坝体稳定分析的基础。本文以某拟建大型土石坝工程超深厚覆盖层地基中的砂土层为研究对象，对原状土样（通过深围井直接开挖取得的）和重塑土样开展了一系列的高围压条件下的不排水循环三轴试验。研究表明：原状砂土的动强度高于重塑砂土，且随着围压的增大二者之间的差别逐渐减小；原状土样和重塑土样的孔压、轴向应变发展规律类似；固结应力比对动强度指标及土样的破坏模式都有明显的影响，随着固结应力比的增大，动强度增大，土样破坏模式逐渐由液化破坏转变为塑性应变累积破坏。Abstract: The ultra-deep overburden foundation is a great challenge for the construction of water conservancy and hydropower projects. The dynamic mechanical properties of overburden soil are the basis for discriminating liquefaction of foundation soil and analyzing dam stability. The sand layer in the ultra-deep overburden foundation of a proposed large-scale earth rock dam project is taken as the research object. A series of undrained cyclic triaxial tests are conducted on the undisturbed soil samples (directly extracted from an ultra-deep in-situ test well) and remolded soil samples under high confining pressure conditions. The laboratory results demonstrate that the dynamic strength of the undisturbed sand is higher than that of the remolded sand, and the difference between them gradually decreases with the increase of the confining pressure. The evolution laws of pore pressure and axial strain of the undisturbed and remolded soil samples are similar. The consolidation stress ratio has obvious influences on both the dynamic strength index and the failure mode of the soil samples. With the increase of the consolidation stress ratio, the dynamic strength increases, and the failure mode of the soil samples gradually changes from liquefaction failure to plastic strain accumulation one.
-
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。随着孔隙率减小，再生骨料强夯桩的承载性能逐渐提高。
-
![]()
图 1 深围井及现场取样图
Figure 1. Formation of in-situ test well[13] and photos of field sampling of soil materials
![]()
图 4 原状土样在不同固结条件下的动强度曲线
Figure 4. Dynamic strength curves of undisturbed soil samples under different consolidation conditions
![]()
图 5 原状土样和重塑土样动强度曲线对比（10%标准）
Figure 5. Comparison of dynamic strength curves between undisturbed and remoulded soil samples (strain standard of 10%)
![]()
图 6 CRR（20）和Ks(20)随围压变化图
Figure 6. Variation of CRR (20) and Ks(20) with confining pressure
![]()
图 8 不同固结应力比条件下土样的u/uf-N/Nf关系曲线（以σ3=1600 kPa为例）
Figure 8. u/uf-N/Nf relationship curves of soil samples under different consolidation stress ratios (taking σ3=1600 kPa as an example)
![]()
图 9 典型的有效应力路径（Kc=1.0）
Figure 9. Typical effective stress paths under consolidation condition of Kc=1.0
![]()
图 10 典型的有效应力路径(Kc=1.5)
Figure 10. Typical effective stress paths under consolidation condition of Kc=1.5
![]()
图 11 典型的有效应力路径(Kc=2.0)
Figure 11. Typical effective stress paths under consolidation condition of Kc=2.0
表 1 试验方案表
Table 1 Summary of test schemes
土样类型 固结应力比Kc 围压σ3/kPa 原状样 1.0，1.5，2.0 200，400，800，1600 重塑样 1.0 200，400，800，1600 说明：每个围压条件下4个试样，振动破坏周次分别在10，20，50，100左右。 
下载: 导出CSV
-
[1] 余挺, 邵磊. 含软弱土层的深厚河床覆盖层坝基动力特性研究[J]. 岩土力学, 2020, 41(1): 267-277. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202001032.htm YU Ting, SHAO Lei. Study of dynamic characteristics of dam foundation on deep riverbed overburden with soft soil layer[J]. Rock and Soil Mechanics, 2020, 41(1): 267-277. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202001032.htm
[2] HAN B, ZDRAVKOVIC L, KONTOE S, et al. Numerical investigation of the response of the Yele rockfill dam during the 2008 Wenchuan earthquake[J]. Soil Dynamics and Earthquake Engineering, 2016, 88: 124-142. doi: 10.1016/j.soildyn.2016.06.002
[3] WANG W B, HÖEG K, ZHANG Y B. Design and performance of the Yele asphalt-core rockfill dam[J]. Canadian Geotechnical Journal, 2010, 47(12): 1365-1381. doi: 10.1139/T10-028
[4] 李鹏, 杨兴国, 薛新华, 等. 瀑布沟心墙堆石坝地震响应分析[J]. 世界地震工程, 2015, 31(1): 217-223. https://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201501032.htm LI Peng, YANG Xingguo, XUE Xinhua, et al. Earthquake response analysis of Pubugou corewall rockfill dam[J]. World Earthquake Engineering, 2015, 31(1): 217-223. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201501032.htm
[5] YU X, KONG X J, ZOU D G, et al. Linear elastic and plastic-damage analyses of a concrete cut-off wall constructed in deep overburden[J]. Computers and Geotechnics, 2015, 69: 462-473. doi: 10.1016/j.compgeo.2015.05.015
[6] 杨正权, 刘启旺, 刘小生, 等. 厄瓜多尔CCS水电站深厚覆盖层火山灰沉积土动强度特性试验研究[J]. 水利学报, 2014, 45(增刊2): 161-166. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB2014S2026.htm YANG Zhengquan, LIU Qiwang, LIU Xiaosheng, et al. Study on dynamic strength characteristics tests for deposited soils of volcanic ash[J]. Journal of Hydraulic Engineering, 2014, 45(S2): 161-166. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB2014S2026.htm
[7] 杨玉生, 刘小生, 李小泉, 等. 固结应力状态对超深厚覆盖层深埋砂土动强度参数的影响[J]. 水利学报, 2016, 47(4): 518-526. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201604007.htm YANG Yusheng, LIU Xiaosheng, LI Xiaoquan, et al. Effects of effective confining stresses on cyclic resistance ratio of deep buried sands in deep alluvial soils[J]. Journal of Hydraulic Engineering, 2016, 47(4): 518-526. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201604007.htm
[8] EL-SEKELLY W, DOBRY R, ABDOUN T. Assessment of state-of-practice use of field liquefaction charts at low and high overburden using centrifuge experiments[J]. Engineering Geology, 2023, 312: 106952. doi: 10.1016/j.enggeo.2022.106952
[9] 周燕国, 谭晓明, 陈捷, 等. 易液化深厚覆盖层地震动放大效应台阵观测与分析[J]. 岩土工程学报, 2017, 39(7): 1282-1291. doi: 10.11779/CJGE201707015 ZHOU Yanguo, TAN Xiaoming, CHEN Jie, et al. Observations and analyses of site amplification effects of deep liquefiable soil deposits by geotechnical downhole array[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(7): 1282-1291. (in Chinese) doi: 10.11779/CJGE201707015
[10] OKA L G, DEWOOLKAR M M, OLSON S M. Liquefaction assessment of cohesionless soils in the vicinity of large embankments[J]. Soil Dynamics and Earthquake Engineering, 2012, 43: 33-44. doi: 10.1016/j.soildyn.2012.06.012
[11] JANA A, STUEDLEIN A W. Dynamic in situ nonlinear inelastic response of a deep medium dense sand deposit[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2021, 147(6): 04021039. doi: 10.1061/(ASCE)GT.1943-5606.0002523
[12] WANG L, YANG Z Q, ZHAO J M, et al. Transmitting characteristics of seismic motion in super-deep overburden layer ground[J]. Shock and Vibration, 2021, 2021: 8898012.
[13] YU T, JIANG B N, HE S B, et al. Study on the seepage control of ultra-deep in situ test wells in deep overburden of riverbeds[J]. Arabian Journal of Geosciences, 2022, 15(9): 834. doi: 10.1007/s12517-022-10111-8
[14] WANG L, YANG Z, ZHAO J, et al. Seismic response analysis of earth-rock fill dam on deep overburden under viscoelastic boundary condition[J]. Revista Internacional De Métodos Numéricos Para Cálculo y Diseño En Ingeniería, 2021, 37(1): 1-8.
[15] 邵磊, 余挺, 徐强. 无限元在超深厚覆盖层坝基动力分析中的应用[J]. 地下空间与工程学报, 2017, 13(6): 1537-1543. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201706016.htm SHAO Lei, YU Ting, XU Qiang. Dynamic analysis on super-deep overburdened layer foundation based on infinite element method[J]. Chinese Journal of Underground Space and Engineering, 2017, 13(6): 1537-1543. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201706016.htm
[16] 罗启迅, 丁庆, 李建国, 等. 用于批量快速获取原状土样的施工方法: CN114739729A[P]. 2022-07-12. LUO Qixun, DING Qing, LI Jianguo, et al. Construction Method for Rapidly Obtaining Undisturbed Soil Samples in Batches: CN114739729A[P]. 2022-07-12. (in Chinese)
[17] 丁庆, 罗启迅, 田先忠, 等. 用于批量获取原状土样的砂浆材料及取样施工方法: CN114656229B[P]. 2023-04-28. DING Qing, LUO Qixun, TIAN Xianzhong, et al. Mortar Material for Batch Acquisition of Undisturbed Soil Samples and Sampling Construction Method: CN114656229B[P]. 2023-04-28. (in Chinese)
[18] 水电水利工程土工试验规程: DL/T 5355—2006[S]. 北京: 中国电力出版社, 2006. Code for Soil Tests for Hydropower and Water Conservancy Engineering: DL/T 5355—2006[S]. Beijing: China Electric Power Press, 2006. (in Chinese)
[19] SEED H B, IDRISS I M, LEE K L, et al. Dynamic analysis of the slide in the lower San Fernando Dam during the earthquake of February 9, 1971[J]. Journal of the Geotechnical Engineering Division, 1975, 101(9): 889-911. doi: 10.1061/AJGEB6.0000195
计量
- 文章访问数: 124
- HTML全文浏览量: 55
- PDF下载量: 24
