加筋土桩板墙联合支护的高填方边坡稳定性计算

0 引言

在实施山区工程的选址、选线、选站时，常遇高填方边坡。由于填方高度过大，不能进行自然放坡，因此必须采取支护措施进行加固。加筋土挡墙由于造价低廉，施工方便，在高填方边坡中得到了广泛应用。加筋土技术于20世纪60年代由法国工程师Henri.Vidal提出，90年代以后，国际工程界和学术界又继续在加筋机理研究、完善加筋土结构设计计算理论、拓宽应用范围、数值模拟分析、影响因素分析等方面做了大量的工作也取得了许多重要成果［1-8］。目前的加筋土设计理论方面，应用最多的是极限平衡法。用极限平衡法计算需假定土体为理想刚塑性体，并假定滑动面的位置，考虑筋材的拉力对土体稳定的贡献，根据静力平衡计算潜在滑动面的稳定性，在分析计算中并未考虑筋土相互作用，计算结果精度不高［9］。而有限元强度折减法不需任何简化和假定，可以直接得到破裂面的位置和安全系数，并能考虑筋土相互作用，因此更加适用于加筋土挡墙的计算［10-12］。

随着工程条件越来越复杂，单一的支护结构很难满足支护要求，因此经常采用多种支护方式联合支护，如加筋土挡墙和桩板墙联合支护。国内外学者分别围绕加筋土挡墙和桩板墙开展了很多工作，但是关于加筋土桩板墙联合支护的研究很少。唐晓松等［13］通过有限元强度折减法计算加筋土挡墙-抗滑桩组合支挡结构，得到治理工程的稳定安全系数，进行支挡结构的受力分析，并通过研究计算参数对工程稳定性的影响，对支挡方案进行优化。薛鹏鹏等［14］采用有限分软件FLAC3D建立加筋土挡墙-抗滑桩组合支差挡结构数值分析模型，着重分析了不同面板浇筑方式对加筋土挡墙墙面水平位移、墙背土压力、桩身水平位移、桩身弯矩和土工格栅应力分布的影响。但前人在分析加筋土和桩板墙联合支护结构的边坡稳定性时，主要关注稳定安全系数，没有综合分析各支挡结构的内力、位移等；在分析稳定安全系数的影响因素时，多考虑筋带的影响，较少考虑抗滑桩桩长对结构稳定性及结构内力、位移的影响。

本文结合实际工程，采用PLAXIS有限元软件对加筋土桩板墙联合支护的高填方边坡进行分析计算。主要目的是对高填方加筋土桩板式挡墙结构的破坏模式进行研究，分析边坡的破坏形式、稳定安全系数、土工格栅内力以及桩板墙的内力和位移等，探讨影响高填方边坡稳定性的各类影响因素。从而对高填方边坡的优化加固方案选择提出建议，为高填方加筋桩板式挡墙结构的整体稳定性分析和治理设计提供参考。

1 工程概况

山西某山区工程场地复杂，形成高35 m的人工填方边坡。由于填方高度过大，采用加筋土和桩板墙联合支护的方式进行加固。高填方边坡共分三级，每级高8m，均为1∶1放坡。边坡土体加筋，筋带采用有纺土工布，长度为15m，垂直间距为0.5m，筋带抗拉强度为200 kN／m。在边坡下部设置了桩板墙，总长度为22m，悬臂段长11m，抗滑桩截面尺寸为1.8m×2.4m，桩间距为4m。桩板墙侧土体加筋，筋带的长度为26m，垂直间距为0.5m，筋带抗拉强度为200 kN／m。

图1 工程现状示意图 Fig.1 Sketch photo of present situation of the engineering

2 计算方法和计算参数

2.1 计算方法

相对于极限平衡法，有限元强度折减法不需对土体和破裂面进行简化和假定，可以直接得到破裂面的位置和安全系数，并能够考虑筋土相互作用，因此本文利用PLAXIS有限元程序，采用有限元强度折减法进行计算。

桩板墙的位移、内力分布如图6所示。分析可知，桩板墙最大位移发生在悬臂段顶端，最大值为0.03 m，在允许范围之内。抗滑桩剪力最大值为4114 kN，而弯矩最大值为8442 kN·m，并且桩的内力最大值发生在抗滑桩11 m长度处，并有明显的应力集中。

PLAXIS具有自动进行有限元强度折减的功能，所谓强度折减［14］就是将坡体的真实抗剪强度指标c和φ除以一个系数，即折减系数Fs，然后用折减后的抗剪强度代替土体的实际强度，土体达到极限临界平衡状态时的折减系数即为边坡的安全系数。本文利用PLAXIS程序内嵌的Phi-c折减过程，对所有岩土体进行折减，通过不断降低土体抗剪强度参数，使有限元平衡方程反复计算，直到计算不收敛，得到边坡最危险滑裂面，从而得到边坡稳定系数。

在天然湿度状态下的土层中开挖沟槽，且地下水位低于槽底时可开直槽，不设支撑，但对槽深有限制要求：砂土和砂砾土土层的深度不大于1.0 m；亚砂土和亚黏土土层的深度不大于1.25 m；黏土土层的深度不大于1.5 m。

PLAXIS程序中系数∑Msf定义为强度的折减系数，其表达式为：

为研究土体性质的影响，改变填土c、φ值，分析边坡安全系数、结构内力及位移随填土c、φ值的变化规律。为研究土体c值的影响，取φ值为38°，c值从6 kPa增至12 kPa，计算结果如表7所示。由表可见，随着c值的增大，边坡安全系数明显提高，由1.35提高到1.47，桩板墙位移基本不变，剪力和弯矩略有增大；为研究土体φ值的影响，取 c值为 12 kPa，φ 值从30°增至 42°，计算结果如表8所示。由表可知，随着φ值的增大，安全系数显著提高，由1.06提高到1.45，桩板墙位移不变，剪力和弯矩略有增大。

2.2 材料参数

根据实际工程，选择弹塑性的摩尔-库仑材料模型进行数值模拟，土层分布如图2，岩土体物理力学参数的选取如表1所示，既有地基土和坡体回填土（分层碾压碎石土，土石比8∶2）的参数根据工程勘察报告选取。

表1 岩土物理力学参数 Tab.1 Physical and mechanical parameters of rock and soil

材料名称 重度g／kN·m-3弹性模量E／MPa泊松比μ粘聚力c／kPa内摩擦角φ／°回填土 18 30 0.3 5 38地基土 20 30 0.3 12 35中风化石灰岩 26 30000 0.3 220 35

图2 土层分布示意图 Fig.2 Sketch map of soil distribution

2.3 土工格栅的模拟

筋带用土工格栅单元模拟，土工格栅是具有轴向刚度的细长形结构，只能承受拉力，不能承受压力，因此土工格栅单元的本构关系近似为线弹性单元。筋土之间的相互作用通过界面单元来实现。土工格栅唯一的材料性质是弹性轴向刚度EA，本工程采用的土工布极限抗拉强度为F＝200 kN／m，当应变为10%时基本达到其极限强度，由EA＝F／10%可以算得，土工布的轴向刚度为2000 kN／m。

在脓毒症的免疫抑制阶段，由于胸腺萎缩及外周血中T淋巴细胞减少可导致机体免疫监视功能减弱，对病原体的易感性增强，增加二次感染的可能性[32]。在多细菌诱导的脓毒症动物模型中，炎症早期使用S1P类似物(FTY-720)抑制T淋巴细胞从淋巴器官迁出，对小鼠生存率无显著影响[33]。然而，在柠檬酸杆菌诱导的脓毒症小鼠中，FTY-720长期治疗的小鼠体内病原体清除延迟，且其脾内的细菌含量显著增加[33-34]。

2.4 桩板墙的模拟

由于桩侧土体已有桩板墙的支护作用，因此考虑桩侧土体不加筋，计算模型如图7所示。计算结果如表4，计算得到边坡的安全系数为1.38，同加筋时相同，这主要是因为，土体沿第一级边坡底部滑出，潜在破裂面如图4所示，桩侧土体并未破坏，因此加固桩侧土体对提高边坡安全系数并无作用。桩板墙位移不变，桩板墙剪力和弯矩均增大3%，但能节省筋带44%左右，桩板墙内力可通过提高配筋率等措施降低，因此综合考虑，桩侧可不加筋。

表2 抗滑桩参数 Tab.2 Parameters of anti-slide pile

材料名称 截面宽度b／m截面高度h／m弹性模量E／MPa抗弯刚度EI／kN·m抗滑桩 1.8 2.4 30000 60×e6

3 稳定性分析

3.1 稳定安全系数分析

建立PLAXIS2D模型，如图3所示。计算得边坡稳定安全系数为1.43，大于规范要求的1.35，满足稳定性要求。滑裂面的位置如图4所示。从图中可以看出，滑裂面位置出现在加筋土体的边缘，边坡破坏时，土体沿第一级加筋土挡墙的底部滑出。

图3 有限元模型图 Fig.3 Finite element model

图4 滑裂面位置 Fig.4 Location of sliding surface

3.2 土工格栅内力分析

通过建立有限元模型，利用有限元强度折减法，计算得到了边坡稳定安全系数和结构内力及位移，并探究了高填方边坡稳定性的各类影响因素。研究初步结论如下：

图5 筋带内力沿高度分布图 Fig.5 Distribution of internal force of geogrid along the height

3.3 桩板墙位移及内力分析

PLAXIS是荷兰开发的岩土工程有限元软件，该程序能够模拟复杂的工程地质条件。该程序能够计算两类工程问题，即平面应变问题和轴对称问题，本文的高填方边坡建立模型时即按照平面应变问题考虑。

一是引才要更加精准务实。引进人才要结合本地实际，明确引才目标和优惠条件，切忌盲目贪大贪洋。大连是老龄化严重和产业转型升级任务繁重地区，亟待引进年轻人才、技能型人才、创新创意人才、跨界复合型人才、新学科发展人才，不仅瞄准高端人才，更要储备大量中高级专业化技能人才，创造条件将大连高等院校培养的大量人才优先吸引在连就业创业，切实改善区域人才结构和补充数量短板。

图6 桩板墙位移及内力分布 Fig.6 Displacement and internal force of pile slab wall

4 影响因素分析

4.1 筋带强度

分析筋带强度对桩板墙稳定性的影响，从而对设计方案进行优化。目前工程中常用的几种土工布强度为 200 kN／m、150 kN／m、100 k／m，表 3是三种不同筋带强度下的计算结果。从表中可以看出，随着筋带强度的降低，边坡安全系数减小，桩板墙位移不变，剪力和弯矩略有增大。筋带强度为100kN／m时，与筋带强度为200 kN／m时相比，桩板墙剪力仅增大0.8%，弯矩仅增大0.4%，安全系数为1.38，满足规范要求的1.35，因此选用抗拉强度为100 kN／m的土工布即可满足要求。

4.2 筋带长度

桩板墙采用PLAXIS中的板单元模拟，板的材料性质中最重要的参数是抗弯刚度EI和轴向刚度EA。本工程中抗滑桩的参数如表2所示，需将桩的抗弯刚度除以桩间距转换成板单元的抗弯刚度。计算得到板单元的抗弯刚度EI为15×e6 kN·m，轴向刚度EA为32.4×e6kN／m。

表3 不同筋带强度下的计算结果 Tab.3 Calculation results under different geogrid strength

筋带强度／kN·m-1 安全系数 桩板墙最大位移／m 桩板墙最大剪力／kN 桩板墙最大弯矩／kN·m 200 1.43 0.03 4057 8449 150 1.40 0.03 4069 8454 100 1.38 0.03 4088 8488

表4 桩侧加筋和无筋计算结果对比 Tab.4 Comparison of calculation results of reinforcement and non-reinforcement on side of pile

桩侧加筋情况 安全系数 桩板墙最大位移／m 桩板墙最大剪力／kN 桩板墙最大弯矩／kN·m桩侧有筋 1.38 0.03 4088 8488桩侧无筋 1.38 0.03 4211 8744

图7 桩侧不加筋有限元模型图 Fig.7 Finite element model diagram of pile side without reinforcement

由图4可知，边坡破坏时土体沿第一级加筋土底部滑出，从图5可见第一级加筋土挡墙筋带内力较大，因此其长度不再优化，对第二级和第三级筋带长度进行优化。如表5所示，随着筋带长度的减小，边坡安全系数降低，桩板墙位移不变，桩板墙剪力和弯矩基本不变。优化后的筋带长度为：第一级挡墙筋带长度15 m，第二级挡墙筋带长度14 m，第三级挡墙筋带长度13 m，这时的稳定安全系数为1.37，小于规范要求的1.35，满足要求。

天门冬蛋白酶（Aspartic proteinase）在动植物中均有存在，动物体内主要集中于乳汁中，在土豆和黑麦的叶中均有提取价值。土豆叶天门冬蛋白酶的分子量约40 000，在pH为3时活性最好。

4.3 桩长

为考虑桩长的影响，在筋带长度优化的基础上减小桩长（优化后的筋带长度为：第一级挡墙筋带长度15 m，第二级挡墙筋带长度14 m，第三级挡墙筋带长度13 m），研究桩长的影响，从而对设计方案进行优化。计算结果如表6所示。由计算结果可知，随着桩长的减小，边坡稳定安全系数降低，桩板墙位移不变，桩板墙剪力和弯矩减小。当桩长减小到悬臂段和埋入段均为9 m时，安全系数降低到1.25，此时已不满足规范要求的1.35。因此桩长取悬臂段和嵌固段均为10 m，总桩长为20 m。因此，优化后边坡安全系数为1.35，能够满足规范要求。

表5 不同筋带长度下的计算结果 Tab.5 Calculation results under different geogrid length

筋带长度／m 安全系数 桩板墙最大位移／m 桩板墙最大剪力／kN 桩板墙最大弯矩／kN·m二级15三级15 1.38 0.03 4211 8744二级14三级14 1.38 0.03 4211 8744二级13三级13 1.34 0.03 4210 8743二级14三级13 1.37 0.03 4211 8745二级14三级12 1.32 0.03 4211 8745

表6 不同桩长下的计算结果 Tab.6 Calculation results under different pile lengths

桩长／m 安全系数 桩板墙最大位移／m 桩板墙最大剪力／kN 桩板墙最大弯矩／kN·m悬臂11 m，嵌固11 m 1.37 0.03 4211 8745悬臂10 m，嵌固10 m 1.35 0.03 3741 7799悬臂9 m，嵌固9 m 1.25 0.03 3634 7572

4.4 土体性质

其中，tanφinput和 cinput为实际输入的强度参数值；tanφreduced和creducedt为经过折减后的强度参数值。

表7 不同土体c值下的计算结果 Tab.7 Calculation results under different cof soil

粘聚力c／kPa 安全系数 桩板墙最大位移／m 桩板墙最大剪力／kN 桩板墙最大弯矩／kN·m 6 1.35 0.03 3791 7873 7 1.37 0.03 3817 7928 8 1.38 0.03 3849 7994 9 1.41 0.03 3872 8042 10 1.43 0.03 3893 8084 11 1.46 0.03 3911 8121 12 1.47 0.03 3927 8155

表8 不同土体φ值下的计算结果 Tab.8 Calculation results under different φ of soil

内摩擦角φ／° 安全系数 桩板墙最大位移／m 桩板墙最大剪力／kN 桩板墙最大弯矩／kN·m 30 1.06 0.03 3700 7685 32 1.12 0.03 3684 7651 34 1.17 0.03 3700 7685 36 1.28 0.03 3729 7746 38 1.35 0.03 3755 7799 40 1.40 0.03 3777 7845 42 1.45 0.03 3794 7881

因此，由上述分析可见，土体性质对边坡稳定安全系数影响显著，当土体c、φ值提高时，安全系数明显提高，因此工程中要优先选用性质较好的填土，当高填方边坡的填土选用受限时，可采用分层碾压等施工工艺改善土体性质。

5 结论

图5给出了加筋土施工完成后土工格栅单元内力沿高度分布（地面至35 m高度处），土工格栅受力最大的位置发生在11 m高处，即桩板墙顶部位置。

本研究涉及的研究对象为中专起点的高职护生，他们以前系统学习过护理学基础理论和专业课程，绝大多数（94.64%）已经通过了护士执业资格考试，在学习目的上更多是提高学历和学习新知识理论。护生多数有过临床经验，在教学方法方面接受案例教学的护生占100.00%，体现其对能力锻炼的要求；也有71.43%的护生接受以教师讲授为主的教学方法。可见护生在临床实践中，充分认识到基础知识的扎实程度直接影响工作成果，对系统学习新方法技术的要求很迫切。这一点也提示我们要在日常教学中分析不同教学内容，选用适合的教学方法[2]。

（1）本工程优化后的方案为桩侧无筋，第一级挡墙筋带长度为15 m，第二级挡墙筋带长度14 m，第三级挡墙筋带长度13 m，筋带选用抗拉强度为100 kN的土工布。抗滑桩长度20 m，悬臂段10 m。优化后的边坡稳定安全系数为1.35，满足规范要求。

河道周边的绿化，要做到层次化和空间化,种植密度不宜过大，树木种类要较为丰富些，不同的流域要有不同的景观。通常会以高达的乔木作为背景，将亚乔木、灌木、草花及地被作为组团。水生植物和湿地植物要注重色彩的搭配，要体现出水体的美感。设计人员在设计城市河道道路景观时，要坚持自然、淳朴的思想理念，做好将城区运输道路与休闲道路分隔开，将河道周围的道路与植物区域、小区公园等相融合在一起，设计出丰富多变的景观区域，让前来休闲放松的人们能够享受到更加舒适的自然美。本项目市政道路改造措施为，新建8条市政道路，总里程约5km；新建5座景观桥梁。

（2）边坡安全系数随筋带强度的提高而增大，随筋带长度的提高亦增大，但当筋带长度过大时，对提高安全系数作用不大，所以在实际工程中要合理选择筋带长度。

（3）边坡安全系数随抗滑桩长度增大而增大，但长度过大时，对安全系数提高作用很小，因此应合理设计桩长。

对比工况一和工况二可以发现，幅流风机可以有效地扰动满载地铁车厢内的气流，能将空调出风更加均匀地送至“气流死区”。而对于工况二和工况三，即有幅流风机的两个工况，在幅流风机的作用下车厢内的风速度分布几乎无区别，最大风速均出现在幅流风机出风口处，风速约为2.8 m/s，截面平均风速约为0.5 m/s。

（4）通过研究发现，土体性质对边坡稳定安全系数影响显著，因此工程中要优先选用性质较好的填土。当高填方边坡的填土选用受限时，可采用分层碾压等施工工艺改善土体性质。

参考文献

［1］ Andrzej Sawicki.Creep of geosynthetic reinforced soil retaining walls［J］.Geotextiles and Geomembranes，1999， 17：51-65.

［2］ Manoucher Nejad Ensan， Isam Shahrour. A simplified elastoplastic macroscopic model for the reinforced earth material［J］.Mechanics.Research Communications， 2000， 27（1）： 79-86.

［3］ M.Bennis，P.De Buhan.A multiphase constitutive model of reinforced soils accounting for soil-inclusion interaction behaviour ［J］. Mathematical and Computer Modeling， 2003， 37： 469-475.

［4］ 闫澍旺，Ben Barr.土工格栅与土相互作用的有限元分析［J］. 岩土工程学报，1997，19（6）：56-61.

［5］ 张孟喜，孙钧.土工合成材料加筋土应变软化特性及弹塑性分析［J］.土木工程学报，2000，33（3）： 104-107.

［6］ 刘华北，Ling H I.土工格栅加筋土挡土墙设计参数的弹塑性有限元研究［J］.岩土工程学报，2004，26（5）： 668-673.

［7］ 周世良，刘占芳，王多垠，等.格栅加筋土挡墙数值分析的复合材料方法［J］.岩石力学与工程学报， 2006（11）：2327-2334.

［8］ 介玉新.加筋土不同计算方法之间的关系［J］.岩土力学， 2007，32（1）： 43-48.

［9］ 陈祖煜，弥宏亮，汪小刚.边坡稳定三维分析的极限平衡方法［J］. 岩土工程学报，2001 （5）：525-529.

［10］ 宋雅坤，郑颖人，刘一通，等.土工格栅加筋土挡墙稳定性影响因素敏感性分析［J］.后勤工程学院学报，2010，26（3）： 1-7.

［11］ 唐晓松，王永甫，郑颖人.土工格栅加筋土挡墙在复杂工程地质条件下的应用［J］.防灾减灾工程学报，2013，33（S1）： 88-93.

［12］ 陈金锋，宋二祥，徐明.强度折减有限元法在昆明新机场高填方边坡稳定分析中的应用［J］.岩土力学，2011，32（S1）：636-641+703.

［13］ 唐晓松，郑颖人，王永甫.加筋土挡墙与抗滑桩联合支挡的分析计算［J］.土木建筑与环境工程，2011 （S2）：42-46.

［14］ 薛鹏鹏，郑俊杰，曹文昭，等.加筋土挡墙-抗滑桩组合支挡结构数值模拟［J］.长江科学院院报，2017 （2）：75-79.

［15］ 秦海燕.边坡稳定性影响因素敏感性的两种分析结果比较［J］. 防灾科技学院学报，2008，10（4）：32-36.