留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

桩承式加筋路堤有限元几何建模方法及边坡效应研究

姜彦彬 何宁 许滨华 张中流 耿之周 石北啸

引用本文:
Citation:

桩承式加筋路堤有限元几何建模方法及边坡效应研究

    作者简介: 姜彦彬(1989—),男,山东临沂人,博士研究生,主要从事地基处理方面的研究。E-mail:tumujyb@163.com.
  • 基金项目: 国家重点研发计划资助项目(2018YFC1508505,2017YFC0404801);国家自然科学基金面上项目(51579152,51679149);中央级公益性科研院所基本科研业务费专项基金项目(Y317008);国家留学基金资助项目(CSC No. 201808320413)
  • 中图分类号: TU472

Geometrical modelling with finite element method and slope effect of geosynthetic reinforced and pile-supported embankments

  • 摘要: 桩承式加筋路堤常用的有限元几何建模方法包括单桩模型(含轴对称模型和单桩3D模型)、平面应变模型和3D断面模型,这些模型各具特点且前提假设各不相同,但却缺乏对比研究。使用非线性有限元软件ABAQUS结合现场试验结果进行了几何建模,对比分析了基于大面积堆载假定的模型与带边坡3D断面模型应力分布及变形的差异与原因。计算表明:3D断面模型最能反映桩承式加筋路堤真实工况,而单桩模型及大面积3D断面模型均低估了沉降量。边坡效应引起的水平位移及坡外隆起使得桩土沉降差、平均桩土应力比、桩顶荷载分担系数及土工加筋横向拉应力均在路肩区域达到极大值,且在路面范围内均大于大面积堆载工况。
  • 图 1  轴对称模型

    Figure 1.  Schematic diagram of axisymmetric model

    图 2  正方形布桩单桩3D模型

    Figure 2.  Three-dimensional model for single pile with square distributed piles

    图 3  二维平面应变模型

    Figure 3.  Two-dimensional plane strain model

    图 4  3D断面模型

    Figure 4.  Full 3D piled embankments model

    图 5  大面积堆载3D断面模型

    Figure 5.  Large area surcharged full 3D piled embankments model

    图 6  地表桩、土沉降对比

    Figure 6.  Settlement comparison of soil and pile on subsurface level

    图 7  路中测点位置桩、土应力分担

    Figure 7.  Vertical stress distribution of soil and pile at measuring point near center line of embankment

    图 8  路中位置桩、土平均应力分担

    Figure 8.  Average vertical stress distribution of soil and pile at central location of embankment

    图 9  路堤堆载结束时刻3D断面模型位移矢量图和水平位移云图

    Figure 9.  Displacement vector diagram and horizontal displacement contour of full 3D piled embankments model after surcharge filling

    图 10  路堤堆载结束时路堤和地基土大主应力云图(单位:kPa)

    Figure 10.  Principal stress contour of soils and fills after surcharge filling (unit: kPa)

    图 11  路堤堆载结束时各桩位处沉降及荷载

    Figure 11.  Settlement and load distribution between each pile and its surrounding soil after surcharge filling

    图 12  路堤堆载结束时土工加筋横向拉应力云图(单位:kPa)

    Figure 12.  Lateral tensile stress contour of geogrid after surcharge filling (unit: kPa)

    表 1  数值模型所用本构模型及单元类型

    Table 1.  Constitutive models and element type applied in FEM modelling

    模型分部 本构模型 单元类型
    轴对称模型 3D模型
    地下水位之下的地基土 修正剑桥模型 CAX4P C3D8P
    路堤填土、碎石垫层、地表填土 摩尔库伦理想弹塑性模型 CAX4 C3D8
    PCC桩 线弹性模型 CAX4 C3D8
    土工加筋 线弹性模型 MAX1 M3D4
    下载: 导出CSV
  • [1] 娄炎, 何宁, 娄斌. 高速公路深厚软基工后沉降控制成套技术[M]. 北京: 人民交通出版社, 2011: 91-104.

    LOU Yan, HE Ning, LOU Bin. Complete settlement control technology for deep soft foundation of expressway[M]. Beijing: China Communications Press, 2011: 91-104. (in Chinese)
    [2] 姜彦彬, 何宁, 林志强, 等. 路堤深厚软基管桩复合地基数值模拟[J]. 水利水运工程学报,2018(2):43-51. (JIANG Yanbin, HE Ning, LIN Zhiqiang, et al. Numerical simulation of pipe pile composite foundation of deep soft foundation under embankment[J]. Hydro-Science and Engineering, 2018(2): 43-51. (in Chinese)

    [3]

    ZHENG J J, CHEN B G, LU Y E, et al. The performance of an embankment on soft ground reinforced with geosynthetics and pile walls[J]. Geosynthetics International, 2009, 16(3): 173-182. doi: 10.1680/gein.2009.16.3.173
    [4]

    JENCK O, DIAS D, KASTNER R. Three-dimensional numerical modeling of a piled embankment[J]. International Journal of Geomechanics, 2009, 9(3): 102-112. doi: 10.1061/(ASCE)1532-3641(2009)9:3(102)
    [5]

    YAPAGE N N S, LIYANAPATHIRANA D S. A parametric study of geosynthetic-reinforced column-supported embankments[J]. Geosynthetics International, 2014, 21(3): 213-232. doi: 10.1680/gein.14.00010
    [6]

    ARIYARATHNE P, LIYANAPATHIRANA D S, LEO C J. Comparison of different two-dimensional idealizations for a geosynthetic-reinforced pile-supported embankment[J]. International Journal of Geomechanics, 2013, 13(6): 754-768. doi: 10.1061/(ASCE)GM.1943-5622.0000266
    [7]

    ARIYARATHNE P, LIYANAPATHIRANA D S. Review of existing design methods for geosynthetic-reinforced pile-supported embankments[J]. Soils and Foundations, 2015, 55(1): 17-34. doi: 10.1016/j.sandf.2014.12.002
    [8]

    JAMSAWANG P, VOOTTIPRUEX P, BOATHONG P, et al. Three-dimensional numerical investigation on lateral movement and factor of safety of slopes stabilized with deep cement mixing column rows[J]. Engineering Geology, 2015, 188: 159-167. doi: 10.1016/j.enggeo.2015.01.017
    [9]

    ZHUANG Y, WANG K Y. Finite-element analysis on the effect of subsoil in reinforced piled embankments and comparison with theoretical method predictions[J]. International Journal of Geomechanics, 2016, 16(5): 04016011. doi: 10.1061/(ASCE)GM.1943-5622.0000628
    [10] 陈仁朋, 徐正中, 陈云敏. 桩承式加筋路堤关键问题研究[J]. 中国公路学报,2007,20(2):7-12. (CHEN Renpeng, XU Zhengzhong, CHEN Yunmin. Research on key problems of pile-supported reinforced embankment[J]. China Journal of Highway and Transport, 2007, 20(2): 7-12. (in Chinese) doi: 10.3321/j.issn:1001-7372.2007.02.002

    [11]

    ZHUANG Y, ELLIS E, YU H S. Three-dimensional finite-element analysis of arching in a piled embankment[J]. Géotechnique, 2012, 62(12): 1127-1131. doi: 10.1680/geot.9.P.113
    [12]

    ZHANG J, ZHENG J J, CHEN B G, et al. Coupled mechanical and hydraulic modeling of a geosynthetic-reinforced and pile-supported embankment[J]. Computers and Geotechnics, 2013, 52: 28-37. doi: 10.1016/j.compgeo.2013.03.003
    [13]

    GIROUT R, BLANC M, DIAS D, et al. Numerical analysis of a geosynthetic-reinforced piled load transfer platform-validation on centrifuge test[J]. Geotextiles and Geomembranes, 2014, 42(5): 525-539. doi: 10.1016/j.geotexmem.2014.07.012
    [14]

    LIU H L, NG C W W, FEI K. Performance of a geogrid-reinforced and pile-supported highway embankment over soft clay: case study[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2007, 133(12): 1483-1493. doi: 10.1061/(ASCE)1090-0241(2007)133:12(1483)
    [15]

    KHABBAZIAN M, KALIAKIN V N, MEEHAN C L. Column supported embankments with geosynthetic encased columns: validity of the unit cell concept[J]. Geotechnical and Geological Engineering, 2015, 33(3): 425-442. doi: 10.1007/s10706-014-9826-8
    [16]

    LIU K W, ROWE R K. Numerical study of the effects of geosynthetic reinforcement viscosity on behaviour of embankments supported by deep-mixing-method columns[J]. Geotextiles and Geomembranes, 2015, 43(6): 567-578. doi: 10.1016/j.geotexmem.2015.04.020
    [17]

    LIU K W, ROWE R K, SU Q, et al. Long-term reinforcement strains for column supported embankments with viscous reinforcement by FEM[J]. Geotextiles and Geomembranes, 2017, 45(4): 307-319. doi: 10.1016/j.geotexmem.2017.04.003
    [18]

    BHASI A, RAJAGOPAL K. Geosynthetic-Reinforced piled embankments: comparison of numerical and analytical methods[J]. International Journal of Geomechanics, 2015, 15(5): 04014074. doi: 10.1061/(ASCE)GM.1943-5622.0000414
    [19]

    BHASI A, RAJAGOPAL K. Numerical study of basal reinforced embankments supported on floating/end bearing piles considering pile-soil interaction[J]. Geotextiles and Geomembranes, 2015, 43(6): 524-536. doi: 10.1016/j.geotexmem.2015.05.003
  • [1] 杨剑黎冰鲍安琪马文昊 . 考虑土性参数空间变异性的单桩竖向承载力分析. 水利水运工程学报, 2019, (5): 85-90. doi: 10.16198/j.cnki.1009-640X.2019.05.011
    [2] 唐彤芝明经平关云飞赵维炳 . 混凝土芯砂石桩振动沉管成桩效应及其累积特性研究. 水利水运工程学报, 2019, (1): 18-27. doi: 10.16198/j.cnki.1009-640X.2019.01.003
    [3] 郭东王建华范怡飞 . 桩靴贯入砂土层时邻近桩挤土压力分析. 水利水运工程学报, 2019, (5): 91-100. doi: 10.16198/j.cnki.1009-640X.2019.05.012
    [4] 周成王一冰王涵谭昌明 . 含砾土坡坡面降雨和坡脚浸泡模型试验研究. 水利水运工程学报, 2019, (6): 96-106. doi: 10.16198/j.cnki.1009-640X.2019.06.011
    [5] 杨敏李宁李宏儒李国锋 . 锚索应力增量法评价边坡稳定性. 水利水运工程学报, 2019, (2): 8-15. doi: 10.16198/j.cnki.1009-640X.2019.02.002
    [6] 陈宇吕杨齐广政孙熙平 . 船舶撞击荷载作用下高桩墩的破坏分析. 水利水运工程学报, 2019, (2): 25-32. doi: 10.16198/j.cnki.1009-640X.2019.02.004
    [7] 冯云芬高树飞 . 基于位移的高桩码头地震易损性分析. 水利水运工程学报, 2019, (3): 76-84. doi: 10.16198/j.cnki.1009-640X.2019.03.010
    [8] 林国财谢兴华阮怀宁朱珍德卢斌徐晨城路晓刚 . 降雨入渗边坡非饱和渗流过程及稳定性变化研究. 水利水运工程学报, 2019, (3): 95-102. doi: 10.16198/j.cnki.1009-640X.2019.03.012
    [9] 徐磊颜士开施伟王刚金玉杰周歆宇孙世宏 . 虹吸式出水流道水力性能数值计算湍流模型适用性. 水利水运工程学报, 2019, (4): 42-49. doi: 10.16198/j.cnki.1009-640X.2019.04.006
    [10] 邱莉婷沈振中马福恒聂柏松 . MLR-Legendre多项式模型在混凝土坝裂缝开度预测中的应用. 水利水运工程学报, 2019, (2): 111-119. doi: 10.16198/j.cnki.1009-640X.2019.02.016
    [11] 李国英韩朝军魏匡民米占宽 . 考虑坝体-地基接触效应的特高心墙堆石坝结构安全性研究. 水利水运工程学报, 2019, (6): 107-115. doi: 10.16198/j.cnki.1009-640X.2019.06.012
  • 加载中
图(12) / 表(1)
计量
  • 文章访问数:  235
  • HTML全文浏览量:  163
  • PDF下载量:  8
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-12-06
  • 刊出日期:  2020-02-01

桩承式加筋路堤有限元几何建模方法及边坡效应研究

    作者简介: 姜彦彬(1989—),男,山东临沂人,博士研究生,主要从事地基处理方面的研究。E-mail:tumujyb@163.com
  • 南京水利科学研究院 水利部土石坝破坏机理与防控技术重点实验室,江苏 南京 210029

摘要: 桩承式加筋路堤常用的有限元几何建模方法包括单桩模型(含轴对称模型和单桩3D模型)、平面应变模型和3D断面模型,这些模型各具特点且前提假设各不相同,但却缺乏对比研究。使用非线性有限元软件ABAQUS结合现场试验结果进行了几何建模,对比分析了基于大面积堆载假定的模型与带边坡3D断面模型应力分布及变形的差异与原因。计算表明:3D断面模型最能反映桩承式加筋路堤真实工况,而单桩模型及大面积3D断面模型均低估了沉降量。边坡效应引起的水平位移及坡外隆起使得桩土沉降差、平均桩土应力比、桩顶荷载分担系数及土工加筋横向拉应力均在路肩区域达到极大值,且在路面范围内均大于大面积堆载工况。

English Abstract

  • 桩承式加筋路堤是深厚软基地区公路、铁路建设中常用的地基处理技术手段。该技术通过土拱效应和加筋效应将附加荷载向复合地基桩顶集中并传递至压缩性更小的深层地基,具有总沉降小、稳定性好、工期短等优点,在土性差、工后沉降要求高的软基处理工程中得到广泛应用[1-3]。桩承式加筋路堤涉及排水固结、边坡稳定和桩土相互作用等土力学问题,其理论研究仍落后于工程应用。有限元数值模拟是一种高效低耗的研究方法,二维及三维有限元模型已广泛应用于桩承式加筋路堤的研究中,其经济性和灵活性是模型试验所不能取代的。根据几何条件的不同,桩承式加筋路堤常用有限元模型包括轴对称模型、平面应变模型、单桩3D模型和3D横断面模型[4-9]。平面应变模型人为地使用“桩墙”简化了群桩的三维特性;轴对称模型及单桩3D模型为大面积堆载工况下的理想单元体模型,虽便于批量建模但无法反映边坡对桩承式加筋路堤系统应力及变形的影响。经过不同程度的简化,模型通常无法完全反映真实的系统响应,最接近现场原貌几何特征的三维模型因其计算成本高而极少被用于多参数对比分析。

    有限元几何建模方法的合理性研究是桩承式加筋路堤数值模拟的重要环节,目前不仅缺乏上述模型的对比研究,也缺乏对边坡效应及大面积假定合理性的探讨。本文使用大型非线性有限元软件ABAQUS,基于现场试验建立有限元数值计算模型,对比桩承式加筋路堤有限元几何建模方法,并探讨边坡效应。

    • 基于大面积堆载及边界对称假定,可使用简化的单桩模型研究群桩基础的受力与变形。轴对称模型和单桩3D模型均属于单桩模型。通常认为桩承式加筋路堤中心位置满足以上条件,取单桩及其加固范围内的地基土连同上部堆载为典型单元,通过等效轴对称简化或3D建模均可建立单桩模型。

      图1所示的轴对称模型[2, 10]图2所示的单桩3D模型[4,9,11-13],虽简易且耗时少,能反映大面积堆载的单桩三维工况,但都无法反映边坡效应和群桩效应。

      图  1  轴对称模型

      Figure 1.  Schematic diagram of axisymmetric model

      图  2  正方形布桩单桩3D模型

      Figure 2.  Three-dimensional model for single pile with square distributed piles

    • 在一定长度范围内,路堤和地基沿纵向形状不发生改变,将等间距规则布置的群桩理想化为垂直于路堤横断面的桩墙,即可将三维问题简化为二维平面应变问题[3, 5-6, 12],再以路中心线为对称轴可建立半幅路堤的平面应变模型(图3)。Ariyarathne等[6]分别从面积等效和强度等效角度出发对桩墙的等效简化方法进行了总结和对比。但平面应变模型忽略群桩基础纵向相互作用,无论采用模量等效还是面积等效都在本质上改变了群桩基础的相互作用维度,其计算结果存在明显不合理性,因此不推荐使用此种建模方法进行机理研究。

      图  3  二维平面应变模型

      Figure 3.  Two-dimensional plane strain model

    • 桩承式加筋路堤系统涉及群桩效应、土拱效应等三维属性明显的问题。单桩模型及平面应变模型均为理想化模型,即在一定程度上忽略了问题的复杂性从而便于批量建模进行参数分析。随着计算机运算速度的不断提升,完整反映现场几何特征的全断面三维有限元数值模型的运行成本逐渐降低,其报道已不鲜见[14-17]。如图4所示,桩承式加筋路堤的3D断面模型尺寸较大,建模时可在横向利用对称性取半幅路堤,在其纵向通常根据布桩方式取N倍(N≥0.5)的桩间距厚度。

      图  4  3D断面模型

      Figure 4.  Full 3D piled embankments model

    • Liu等[14]报道了位于上海北郊的桩承式加筋路堤现场试验及三维有限元研究,该报道数据详实,并被多次引用[6-7, 18-19]。基于该文献分别建立轴对称模型、单桩3D模型和3D断面模型,对比各模型应力、变形特点;同时,为分析大面积堆载假设的合理性,又在3D断面模型基础上建立了大面积3D断面模型(图5)。

      图  5  大面积堆载3D断面模型

      Figure 5.  Large area surcharged full 3D piled embankments model

      所述工程主要软弱土层为10.2 m厚的淤泥质粉质黏土,天然含水率接近液限,不排水抗剪强度为10~25 kPa,工程性质差。路段采用外径1.0 m的现浇混凝土薄壁管桩(PCC桩)复合地基加固处理,桩长16.0 m,桩端位于相对坚硬的粉砂层上,采用正方形布桩,间距3.0 m。路堤填高5.6 m(其中包括0.5 m厚的碎石加筋垫层),分55 d填筑完成。所有观测仪器在路堤施工前布设完毕,观测工作持续至路堤填筑完成后的125 d(现场观测时间共计180 d)。

      具体数值建模计算参数均参考文献[14],不同之处及相关建模要点概述如下:(1) 使用Embedded Region方法模拟加筋与碎石垫层相互作用;(2) 为避免3D断面模型中路堤边坡区域网格畸形而收敛困难,分别建立路堤和复合地基部件,并以摩尔-库伦破坏准则建立接触,摩擦角取地表填土的内摩擦角;(3) 根据Liu等[16]的研究,更改三维有限元模型单元类型为8节点实体单元以节约计算时间(表1);(4) 在ABAQUS中使用修正剑桥本构模型(搭配多孔介质模型),进行地应力平衡前必须设置正的初始地应力。

      表 1  数值模型所用本构模型及单元类型

      Table 1.  Constitutive models and element type applied in FEM modelling

      模型分部 本构模型 单元类型
      轴对称模型 3D模型
      地下水位之下的地基土 修正剑桥模型 CAX4P C3D8P
      路堤填土、碎石垫层、地表填土 摩尔库伦理想弹塑性模型 CAX4 C3D8
      PCC桩 线弹性模型 CAX4 C3D8
      土工加筋 线弹性模型 MAX1 M3D4
    • 现场实测结果是检验数值模型不可或缺的标准。对于地表沉降,由图6(a)可知,3D断面模型的路中及路肩位置桩土沉降计算结果与实测值均匹配良好。单桩模型及大面积3D断面模型的桩、土沉降计算结果相近且均小于实测值,而桩间土沉降表现明显(第180天测量结束时比实测值小26~32 mm);这表明采用基于大面积堆载假定的模型会低估桩承式加筋路堤的沉降量。

      图  6  地表桩、土沉降对比

      Figure 6.  Settlement comparison of soil and pile on subsurface level

      现场桩、土压力传感器均由土压力盒测得,图7为与现场传感器测点位置相对应的应力,图8为桩、土平均应力及平均应力比(轴对称模型须按照单元对应的半径进行积分)。

      图  7  路中测点位置桩、土应力分担

      Figure 7.  Vertical stress distribution of soil and pile at measuring point near center line of embankment

      图  8  路中位置桩、土平均应力分担

      Figure 8.  Average vertical stress distribution of soil and pile at central location of embankment

      对比图78可见:(1) 各模型桩间土测点位置及整体平均竖向应力计算结果与实际测值均吻合较好(图7(b)图8(b)),这说明桩位布置形心位置的桩间土测点具有良好代表性。(2) 多数模型在桩顶测点位置的竖向应力与实测值对应良好,而大面积3D断面模型计算结果明显高于实测值(图7(a));各模型桩顶平均竖向应力及桩土平均应力比均小幅高于实测值(图8(a))。根据姜彦彬等[2]的数值研究,桩顶竖向应力并非均匀分布,无论是数值模型还是现场测量,桩顶有限范围(或数量)的测点或单元均不一定能真实反映桩顶平均应力分担。因此,探讨桩、土平均应力分担更有实际意义。(3) 如图8所示,各模型桩顶及桩间土的平均应力、桩土平均应力比均随填土荷载增加而快速增加,至第55天路堤堆载结束后桩间土平均应力达到峰值,随后小幅缓慢降低,而桩顶平均应力则小幅缓慢上升。实际情况是,当填土完成后,桩、土沉降及其沉降差均随排水固结仍有所增长,土拱效应随之继续发展,一部分由桩间土承担的路堤荷载会被转移至桩顶。因此,第55天填土结束后桩土应力比同样表现为小幅稳定增长,至第180天观测结束时达到峰值(图8(c));此时所有模型中单桩3D模型的平均桩土比最大(为22),会高估桩所承担的荷载;大面积3D断面模型最小(为16),轴对称模型和3D断面模型居中(分别为19和18)。

    • 由于坡外没有附加荷载,在实际桩承式加筋路堤堆载过程中复合地基会发生向坡外的水平位移,同时坡外发生不同程度的隆起;隆起和水平位移大小与荷载水平、地基强度、堆载速度及布桩方式等相关。绘制3D断面模型位移矢量图及水平位移云图(图9),在路堤堆载结束时,地基最大隆起量为17 mm,位于边桩外10.5 m处;地基最大侧向位移56 mm,位于下边坡第二排边桩桩顶(对应图10的8#桩)。模型水平位移云图呈现以边坡为中心的发散状,路中心传感器埋设区域的水平位移量为1~10 mm。

      图  9  路堤堆载结束时刻3D断面模型位移矢量图和水平位移云图

      Figure 9.  Displacement vector diagram and horizontal displacement contour of full 3D piled embankments model after surcharge filling

      图  10  路堤堆载结束时路堤和地基土大主应力云图(单位:kPa)

      Figure 10.  Principal stress contour of soils and fills after surcharge filling (unit: kPa)

      基于大面积堆载假定的单桩模型和3D断面模型均忽略了边坡的存在,无法考虑坡外隆起和向坡外的水平位移,因此低估了桩承式加筋路堤沉降量(图6(a)),其计算结果偏于不安全。

    • 通过对比路堤堆载结束时有边坡及大面积堆载两种3D断面模型的大主应力云图(ABAQUS中以拉为正,土力学中的最大主应力对应其最小主应力)可见:桩顶位置应力明显集中,土拱效应直观(图10),本文形象称之为“蘑菇头”。大面积堆载工况每个桩顶都具有大小和形状相同且左右对称的“蘑菇头”,而带边坡的3D断面模型中“蘑菇头”的对称性自路堤中心向坡脚方向依次减弱,直至变得不规则。边坡的存在使路堤下的地基产生显著的水平位移,从而影响了桩土应力分布,并产生局部应力集中(带边坡3D断面模型大主应力极值超过大面积工况)。

      图11为堆载结束时各桩位处的桩土沉降差与桩土平均应力比、桩顶荷载分担系数E(即各桩加固范围内桩顶承担荷载占路堤总荷载的比值,边坡处除外)关系图。可见,(1) 大面积堆载工况下,桩土沉降差及桩土平均应力比均为常数(分别为45 mm和13);(2) 带边坡3D断面模型自路中心向坡外桩土沉降差先小幅增大,在接近路肩位置的4#桩、5#桩处达到最大值(58 mm)后快速减小;桩土应力比亦呈现先升后降的规律,并在路肩下的6#桩达到最大值29;(3) 除位于下边坡的8#和9#桩外,在路面及上边坡范围内,带边坡的3D断面模型的桩土沉降差及应力比均大于大面积3D断面模型。

      图  11  路堤堆载结束时各桩位处沉降及荷载

      Figure 11.  Settlement and load distribution between each pile and its surrounding soil after surcharge filling

      与之类似,路堤堆载结束时刻,大面积3D横断面模型的土工加筋横向拉应力呈现重复对称分布,而带边坡3D断面模型1#~7#桩位处拉应力均大于大面积堆载工况,且峰值出现在靠近路肩位置的5#桩(图12)。

      图  12  路堤堆载结束时土工加筋横向拉应力云图(单位:kPa)

      Figure 12.  Lateral tensile stress contour of geogrid after surcharge filling (unit: kPa)

      综上所述,路堤荷载作用下的复合地基具有显著的边坡效应,靠近路肩位置的桩间土发生绕流,坡外隆起,使桩土沉降差、平均应力比、桩顶荷载分担系数及土工加筋横向拉应力在路面范围内的桩位处均大于大面积堆载工况,且在路肩附近的桩位处达到极值。

      模型对比分析表明,3D断面数值模型在模拟现场工况受力变形的真实性上具有其他简化模型无可比拟的优势。然而,桩承式加筋路堤是一个复杂的三维系统,影响其工作性状的因素众多,包括几何属性(桩间距、桩径、桩土置换率、填土高度、软基深度和桩长等)、强度属性(桩身模量、填料摩擦角、加筋强度、软土地基模量等)和时间属性(软基渗透系数、路堤填筑速度)等多类因素。仍需进一步深入开展相关研究以提高3D断面数值模型参数分析的工作效率。

    • 结合现场试验总结并对比了桩承式加筋路堤有限元几何建模方法,基于带边坡及大面积堆载的3D断面模型探讨了边坡对桩土变形及应力分布的影响,主要结论如下:

      (1) 相比其他有限元几何模型,三维断面模型计算结果与现场试验桩土变形及应力分担更为契合,能反映桩承式加筋路堤真实工况。

      (2) 基于大面积堆载假定的单桩模型及3D断面模型均无法考虑桩承式加筋路堤的坡外隆起和水平位移,使其沉降计算结果偏小,从而低估了地基沉降量。

      (3) 桩承式加筋路堤边坡效应使得系统存在不可忽略的水平位移和坡外隆起,导致路面范围内桩位处桩土沉降差均大于大面积堆载工况,从而使得平均桩土应力比、桩顶荷载分担系数及土工加筋横向拉应力亦大于大面积堆载工况,且上述参数均在路肩附近的桩位处达到极大值。

      (4) 平均桩土应力较土压力盒所得散点实测值更有价值,应改进土压力盒(尤其是桩顶)布置方式或荷载测量方法以获取更真实的桩土应力分担。

      (5) 对桩承式加筋路堤几何建模方法及边坡效应的探讨均基于所述现场试验工况,更普适的结论要求考虑除几何模型外的本构关系、桩土接触等因素,并使用较为完备的无量纲参数进行三维断面建模分析。

参考文献 (19)

目录

    /

    返回文章
    返回