土石混合料高填方边坡稳定性分析_童艳光.pdf

广东土木与建筑GUANGDONG ARCHITECTURE CIVIL ENGINEERING2023年2月第30卷 第2期FEB 2023Vol.30 No.2DOI:10.19731/j.gdtmyjz.2023.02.005作者简介：童艳光（1968-），男，硕士，高级工程师，主要从事岩土及地质工程的勘测设计及监测测试工作。E-mail：0引言随着时代的高速发展，城市的不断扩张，交通道路设施的建设不断从平地向地形复杂的山林地区修建，形成大量的高填方边坡13。但高填方工程设计和施工要求相对严格，时有发生的高填方边坡滑坡灾害，使得高填方边坡研究显得很有必要4。关于高填方边坡稳定性问题，诸多学者做了大量的研究，如徐坤等人5采用多种极限平衡法对涉及填料的物理力学特性、填筑过程等复杂因素的高填方边坡进行稳定性分析；杨大志等人6用有限元分析土工格室加筋在高填方边坡的应用；侯俊伟等人7用有限元和大型直剪试验研究分析高填方边坡变形失稳机制；叶帅华等人8针对黄土高填方边坡运用有限元探究其影响因素及变形规律。然而如今关于高填方边坡稳定性的研究很少会侧重对于土石混合料加筋土、抗滑桩锚索作用以及土石混合料浸水饱和的研究。因为若由于地形复杂和地质条件的限制，会使得许多高填方边坡无法进行自然放坡9，此时土石混合料加筋土边坡可以发挥其因地制宜、减少放坡率、提高边坡整体稳定性优等点10，满足在复杂地形地质下的工程要求；若高边坡下方土体性质较好，在抗滑桩处设置锚索也能对高填方边坡稳定性起到一定作用11；若在降雨工况下，浸水使土石混合料饱和会使土体强度降低，从而使边坡稳定性下降12。因此，研究土石混合料间的土工格栅、抗滑桩的锚索和浸水饱和的土石混合料对高填方边坡稳定性的影响也是一个值得关注的方向。本文以广州市某热力电厂工程为项目依托，通过有限元数值模拟，研究分析在此工程中，土工格栅、抗滑桩的锚索与土石混合料浸水饱和状态下对土石混合料高填方边坡稳定性的影响。文中采用专业岩土工程有限元软件Midas GTS，对工程中的边坡进行数值模拟13，通过工程资料确定边坡有限元模型的计算参数，通过控制变量的方法去探究土工格栅、抗滑桩的锚索以及土石混合料在浸水饱和状态下的对边坡稳定性的影响，为以后相关地质条件开展的边坡工程提供参考。土石混合料高填方边坡稳定性分析童艳光1，蔡增蛟1，张堃1，梁松鸿1，刘文召2，黎森宇3（1、广州环投花城环保能源有限公司广州510830；2、中国建筑一局（集团）有限公司北京100161；3、广东工业大学土木与交通工程学院广州510006）摘要：随着时代的高速发展，城市的不断扩张，交通道路设施的建设不断从平地向地形复杂的山林地区修建，形成大量的高填方边坡。运用Midas GTS三维数值模拟计算软件分析土石混合料高填方边坡中土工格栅、抗滑桩的锚索和浸水饱和的土石混合料对其稳定性的影响。结果表明：土工格栅和抗滑桩的锚索的设置有利于提升高填方边坡的稳定性；浸水饱和的土石混合料不利于高填方边坡的稳定性，边坡的设计要做好防水措施，防止雨水入渗。以实际工程作为依托，研究成果可为以后相关地质条件开展的边坡工程提供的参考。关键词：高填方边坡；边坡稳定性；土工格栅；Midas GTS；数值模拟中图分类号：U213.1+58文献标志码：A文章编号：1671-4563（2023）02-019-04Stability Analysis of Soil-rock Mixture High-filling SlopeStability Analysis of Soil-rock Mixture High-filling SlopeTONG Yanguang1，CAI Zengjiao1，ZHANG Kun1，LIANG Songhong1，LIU Wenzhao2，LI Senyu3（1、GZEPI Huacheng Environmental Protection Energy Co.，Ltd.Guangzhou 510830，China；2、China Construction First Group Corporation LimitedBeijing 100161；3、Faculty of Civil and Transportation Engineering，Guangdong University of TechnologyGuangzhou 510006，China）AbstractAbstract：With the rapid development of the times and the continuous expansion of cities，the construction of traffic road facilities hasgradually been built from flat to the complex terrain，forming a large number of high-filled slope.Using the Midas GTS 3D value simulation calculation software to analyze the effects of geogrid，anti-slip pile anchor and soil-rock mixture after water saturation in the high-fillingslope of the soil-rock mixture on its stability.The results show that：The setting of geogrid and anti-slip pile anchor is conducive to improving the stability of high-filling slope；The soil-rock mixture after water saturation is not conducive to the stability of high-filling slope andthe design of the slope should be taken to prevent waterproof measures to prevent rainwater from penetrating.Relying on actual projects，research results can be carried out in the future related geological conditions Reference provided by slope engineering.Key wordsKey words：high-filling slope；slope stability；geogrid；Midas GTS；numerical simulation19童艳光，等：土石混合料高填方边坡稳定性分析FEB 2023 Vol.30 No.22023年2月 第30卷 第2期1工程概况广州市某热力电厂土石混合料高填方边坡工程位于广州市花都区，距S114省道直线距离约2 km，在一期工程的南侧和北侧进行扩建，具体位置如图1所示。高填方边坡场地内广泛分布人工填土、粉质黏土、强风化砂岩和泥岩，岩石完整程度大多为极破碎，总体而言，场地内岩土性能情况较为恶劣，本文模型的建立选取高填方边坡P17-P17剖面一处，具体边坡剖面及土层分布如图2所示，其中高架桥为未建部分。2三维有限元计算模型2.1模型建立为了模拟真实的边坡状态，分析其稳定性，根据工程实际，建立了如图3所示的整体三维有限元模型，模型大小取109 m60 m52 m（长宽高），边坡坡度为1 1，充足考虑了边界效应所带来的影响因此在边坡的右侧延伸了一定的距离。锚索如工程实际设置在抗滑桩后且嵌入强风化粉砂质泥岩层，如图4所示；土工格栅在土石混合料之间分层设置。模型的边界条件：底部为固定约束边界，模型四周施加法向约束，其余为自由边界。有限元数值模拟计算中，锚索的模拟采用1D单元植入式桁架去模拟；土工格栅的模拟采用2D单元土工格栅（2D）去模拟；岩土体采用3D单元实体去模拟。锚索、抗滑桩及土工格栅等均采用线弹性本构模型，岩土体采用弹塑性莫尔-库伦本构模型。土工格栅在工作状态下应力应变关系位于弹性范围内，根据土工格栅的拉伸试验结果14，可确定其计算参数。其中，土石混合料的抗剪强度指标由现场大型直剪试验确定15，其试验结果如图5所示，根据抗剪强度包线，可得土石混合料在天然状态下以及浸水饱和状态下的强度参数，其余岩土体强度参数取值均参照项目工程所提供的勘探报告。具体有限元模型中各种材料物理力学参数如表1所示。2.2有限元强度折减法边坡稳定性是指抗滑力与滑动力的两者之比16。Midas GTS的边坡稳定分析采用了基于有限单元法的强度折减法17。有限元强度折减法（SRM）用于计算边坡在事故点或失稳点的安全系数，其原理是：在外荷载保持不变的情况下，不断对边坡抗剪强度参数c，进行折减，使强度不断调整，直到系统达到不稳定状态时，有限元计算将不收敛，即边坡处于临界平衡图1土石混合料高填方边坡工程具体位置Fig.1Specific Location of Soil-rock MixtureHigh-filling Slope Engineering边坡位置图2边坡剖面及其土层分布Fig.2Slope Profile and Soil Layer Distribution（mm）高架桥（未建）油泵房及油罐区110008000200080003000100080008000土工袋反包42.0001 1.0034.0001 1.00结构填土（掺入6%水泥）50.0002500400080004000500 mm级配碎石排水层冠梁3.001.80m底部连梁地面线228加筋材料极限强度300kN/m长22 m，间距1 m加筋材料极限强度300kN/m长22m，间距0.5m桩（3.002.00m）间距2.80m30158根17717.81 860锚索长28 m，自由段8 m，锚固段20 m，标准值300 kN长36 m，自由段10 m，锚固段26 m，标准值300 kN8根17717.81 860锚索50高程/m48464442403836343230282624222018图3Midas GTS有限元模型Fig.3Finite Element Model of Midas GTS土石混合料（其间可设置土工格栅）抗滑桩（桩后部可设置锚索）素填土强风化砂岩中风化砂岩粉质黏土（可塑）粉质黏土（硬塑）强风化砂岩素填土强风化粉砂质泥岩全风化粉砂质泥岩图4锚索在模型中的布置Fig.4Layout of Anchor Cable in the Model图5现场大型直剪试验土石混合料抗剪强度包线15Fig.5Shear Strength Envelope of Soil-rock Mixture inLarge-scale Direct Shear Test on Site050100法向应力n/kPa抗剪强度f/kPa天然状态饱和状态天然状态线性拟合饱和状态线性拟合04015012080200160f=ntan34.7+44.0R2=0.9862f=ntan7.7+45.7R2=0.973250cm厚黏土素填土20广东土木与建筑FEB 2023 Vol.30 No.22023年2月 第30卷 第2期状态，此时折减系数Fs即边坡安全系数18。强度折减法仅用于安全系数大于1的边坡，但是其不需要作假定，分析环境更理想，误差较小19。安全系数计算具体公式如下：c=c/Fstan=tan/Fs=c+tan式中：Fs为强度折减系数；c为折减之后的粘聚力；为折减之后的内摩擦角；为折减之后的抗剪强度。2.3有限元模型的变量为研究分析在此边坡中，土工格栅、抗滑桩的锚索和土石混合料在浸水饱和状态下对土石混合料高填方边坡稳定性的影响。本次有限元模拟设置了3个变量：天然状态下还是浸水饱和状态下的土石混合料；有无设置土工格栅；有无设置抗滑桩的锚索。为探究三者的关系，分别模拟了5种工况，其具体工况设置如表2所示。3数值模拟结果与分析3.1不同工况条件下边坡的稳定系数对5个不同工况下的高填方边坡进行数值模拟计算，得到最终安全系数Fs分别为1.57、1.44、1.50、1.43、1.23。其中在天然土石混合料状态下，即工况 14中，这 4种工况下他们的变量是有无设置土工格栅有无设置抗滑桩的锚索，且他们的滑裂面发生的位置相近，形状类似，以工况1为例，如图6所示。在土石混合料浸水饱和状态下，即工况5，其滑裂面与其他4种工况截然不同，如图6所示。3.2结果分析对土石混合料的不同状态进行分别讨论。3.2.1天然状态下土石混合料从上述结果，可以看出天然状态下土石混合料的4种工况安全系数都处在稳定状态下（Fs1.30），且其滑裂面会发生在土石混合料的下方与右侧，除了土石混合料的右下部有一小部分会形成外，滑裂面其余大部分发生在土质较软的素填土和粉质黏土处，且其形状显得较长且较宽。同时由于有抗滑桩的阻挡作用，无论软土部分还是土石混合料部分，抗滑桩的两侧土体均发生了一定的塑性应变，不过这对于高填方边坡的整体稳定性不会有太大的影响。此时，对比工况2和工况4，其安全系数差别并不大，两种工况所控制的唯一变量是有无土工格栅，有土工格栅的工况相比无土工格栅其安全系数仅提升了0.01，即提升了0.7%，这是由于滑裂面并没有发生在土石混合料填方的边坡，即土工格栅分层设置的地方，因此此时土工格栅的作用并没有被太大地体现出来。在土石混合料天然状态条件下起到影响安全系数或边坡稳定性的，其实更多体现得到的是锚索和抗表1有限元模型各材料物理力学参数Tab.1Physical and Mechanical Parameters of Materials in the Finite Element Model岩土层/结构素填土粉质黏土（可塑）粉质黏土（硬塑）全风化粉砂质泥岩强风化粉砂质泥岩中风化砂岩强风化砂岩天然状态下土石混合料浸水饱和状态的土石混合料土工格栅抗滑桩锚索土层编号11232-重度/kNm-319.020.119.119.920.020.020.020.520.51.025.025.0泊松比0.370.400.430.330.310.250.270.300.300.200.200.28粘聚力/kPa15.020.315.025.526.680.036.644.045.7-内摩擦角/10.018.212.015.919.430.026.534.77.7-标贯击数/击10.714.57.130.565.3-57.3-弹性模量/MPa29.640.219.784.5180.919 000.0158.7800.0800.01 000.030 000.025 000.0工况编号12345模型分析对应的条件天然状态/浸水饱和状态天然状态下天然状态下天然状态下天然状态下浸水饱和状态下土工格栅有有无无有锚索有无有无有表2不同工况下对应设置的条件Tab.2Corresponding Setting Conditions under DifferentWorking Conditions工况1，Fs=1.57图6土石混合料天然状态下有效塑性应变Fig.6Effective Plastic Strain of Soil-rock Mixtureunder Natural Conditions工况5，Fs=1.2321童艳光，等：土石混合料高填方边坡稳定性分析FEB 2023 Vol.30 No.22023年2月 第30卷 第2期滑桩的作用，由图4可以清晰看出，锚索所设置的位置会贯穿滑裂面的下侧，对于抗滑桩的抗倾倒以及边坡滑裂的发生起到一定的作用，对比工况3和工况4，锚索的存在安全系数提升了0.07，即提升了4.9%，因此，锚索对此边坡的稳定性起到很大作用。最后，对比工况1和工况3，设置土工格栅可以使安全系数得到进一步提升，因此，在此高填方边坡中，同时设置锚索以及土工格栅是一个正确的做法。3.2.2浸水饱和状态下的土石混合料从上述结果，可以看出浸水饱和状态下的土石混合料工况下，即工况5，其安全系数仅为1.23，处于基本稳定状态（1.05Fs1.30），其滑裂面发生在抗滑桩顶端的右侧，跨越了整个土石混合料部分，且其形状显得较短且较狭窄，相比天然状态土石混合料工况下的滑裂面有着明显的区别。由于滑裂面大部分发生在抗滑桩的上侧，因此抗滑桩的左右两侧并没有发生明显的塑性应变。对比工况1和工况5，土石混合料浸水饱和状态下与天然状态下相比，安全系数降低了0.34，即降低了21.7%。导致其安全系数有明显的降低，是由于浸水饱和后，土石混合料的土体强度下降明显，在工况5中，由于滑裂面发生的位置横跨土石混合料部分，使得土工格栅的作用得以体现，再次印证了土工格栅的设置对于维持边坡的稳定有着重要的作用。同时，从结果中可以推断出，当降雨的发生，土石混合料不断浸水饱和，其强度随之不断减少，会使得边坡滑裂面不断上移，且长度变短，宽度也变窄，安全系数不断降低，影响高填方边坡的稳定性。在此过程中，根据滑裂面的移动，影响边坡稳定性的主导构件也将会从锚索逐渐转至土工格栅中。4结论广州市某热力电厂土石混合料高填方边坡工程，其工程地质条件复杂，填方高度较大，对合理的边坡支护体系选择提出了较高要求，通过采用压实的加筋土石混合料作为填筑料，联合抗滑桩+锚索的支护体系，确保了高填方工程的顺利实施和安全稳定，得到的主要结论如下：当土石混合料处于天然状态时，高填方边坡处于稳定状态，其稳定安全系数大于1.4，边坡滑裂面出现在原状土层内部且埋置深度和延展宽度较大。当土石混合料浸水饱和状态时，高填方边坡的稳定安全系数降低至1.23，边坡滑裂面出现在土石混合料与原状土层的界面附近且埋置深度和延展宽度较小。土工格栅和锚索的设置均有利于提高填方边坡的稳定性。同时，为确保土石混合料填方边坡的稳定性具有足够的安全储备，需要做好土石混合料的排水防渗工作，避免土石混合料或其他土体浸水软化，强度降低，诱发边坡的失稳破坏。参考文献1 姚仰平，祁生钧，车力文.高填方地基工后沉降计算 J.水力发电学报，2016，35（3）：1-10.2姚仰平，车力文，祁生钧，等.高填方地基蠕变沉降计算方法研究 J.工业建筑，2016，46（9）：25-31.3 吴志轩，张大峰，孔郁斐，等.基-填界面开挖台阶对顺坡填筑高边坡稳定性影响研究 J.工程力学，2019，36（12）：90-97.4杨校辉，朱彦鹏，周勇，等.山区机场高填方边坡滑移过程时空监测与稳定性分析 J.岩石力学与工程学报，2016，35（S2）：3977-3990.5徐坤，周正飞.承德民用机场高填方边坡稳定性分析 J.公路工程，2014，39（2）：251-253+257.6杨大志，崔浩，薛鑫，等.土工格室在机场高填方边坡上的应用研究 J .四川建筑科学研究，2015，41（2）：148-149+153.7侯俊伟，唐秋元，李杨秋，等.西南某山区机场高填方边坡稳定性研究 J.重庆交通大学学报（自然科学版），2016，35（3）：82-88.8 叶帅华，张玉巧，房光文.黄土高填方边坡的稳定性影响因素及其变形规律 J.兰州理工大学学报，2021，47（3）：120-126.9 廖鸿，徐超，杨阳.某机场飞行区土工格栅加筋高边坡优化设计 J.水文地质工程地质，2021，48（6）：113-121.10 何必伍，徐国元，黄文通，等.碎石土混合料在加筋高边坡中的应用 J.河南科技大学学报（自然科学版），2020，41（4）：52-60.11 任洋，李天斌，杨玲，等.基于离心模型试验与数值计算的超高陡加筋土填方边坡稳定性分析 J.岩土工程学报，2022，44（5）：836-844.12 孔郁斐，宋二祥，杨军，等.降雨入渗对非饱和土边坡稳定性的影响 J.土木建筑与环境工程，2013，35（6）：16-21.13 帅红岩，韩文喜，赵晋乾.Midas/GTS软件在边坡三维稳定分析中的应用 J.地质灾害与环境保护，2009，20（3）：104-107.14 郭铭倍.基于有限元数值分析的土工格栅加筋路堤优化设计研究 D.石河子：石河子大学，2020.15 郭明鑫.土石混合料强度特性的大型直剪试验研究 D.广州：广东工业大学，2021.16 孙玉科，李建国.岩质边坡稳定性的工程地质研究 J.地质科学，1965（4）：330-352.17 赵智强，昌思.基于Midas/GTS高陡岩质边坡稳定性分析J.广东土木与建筑，2018，25（3）：54-56.18 代雪，张家明.某场地边坡稳定分析方法的比较研究 J.中国安全生产科学技术，2021，17（11）：119-124.19 李春辉，刘天鹏，王洪洋，等.边坡稳定性分析方法研究进展及展望 J.东北水利水电，2020，38（2）：63-65.22