高边坡治理实体有限元数值模拟及方案比选.pdf

0引言对于土质边坡高度大于20 m、小于100 m或岩质边坡高度大于30 m、小于100 m的边坡，其边坡高度因素会对边坡稳定性产生重要影响，边坡防护加固工程需进行专项设计，这些边坡工程被称为高边坡工程。高边坡治理在公路、市政及房建工程中广泛存在，边坡稳定性一旦遭到破坏，将造成巨大的生命和财产损失。因此，高边坡治理成为影响工程顺利开展的重要环节。由于项目存在高风险性，设计时应对高边坡治理的多个方案进行综合论证。国内学者针对高边坡治理进行了大量研究。王恭先1对滑坡性质及治理方案进行深入研究，提出确定防治方案应考虑的 4 个因素和不同类型滑坡的治理方案，为预防及处理高边坡和路堤填方引起的“工程滑坡”提供了参考。石广斌等2根据边坡崩塌后揭露的地质构造特性，分析高边坡的崩塌机理，并在此基础上拟定加固方案，用极限平衡法对锚索锚固力进行优化，用二维接触非线性有限元法分析边坡岩体与混凝土框架之间的相互作用，不仅提高了边坡岩体的安全度，而且获得了良好的经济效益。目前，针对边坡治理大多采用二维有限元分析法，对于场地及边坡形式简单的施工场景，尚能较好地满足工程安全性的要求。但是，当场地、边坡形式及荷载的空间关系相对复杂时，二维有限元分析法则存在局限性，难以真实地模拟现场实际情况。本文介绍的工程项目的岩质边坡高度为30 m，在其上修建12.8 m的填土边坡及支挡结构，由于周边有道路及学生宿舍，空间关系复杂，安全性要求较高。为更准确地模拟该高边坡工程的治理方案，本文建立实体有限元模型，分析各方案的安全性及经济性，从而确定最合适的方案，并通过监测结果验证有限元模拟的正确性。1工程概况及治理方案设计本工程位于云南昭通某中学南侧场地（如图 1所示），场地东南角分别为2级高边坡及深沟谷，宿舍（5层）紧邻第二级边坡，最短平面距离约 8 m，因此工程对边坡稳定性及位移的要求非常高。第一级边坡坡顶标高为 1 742.00 m，第二级边坡坡顶标高为1 754.80 m。在第一级边坡上进行2次填土，形成2级平台，平台1的标高为1 746.00 m，平台2的标高为1 754.00 m，在平台2上布置8.0 m宽的道路。抗滑桩布置如图2所示（下文中的方案2）：第一排桩桩号为110#，桩径为1.4 m，桩长10.35 m，间距为2.0 m；第二排桩桩号为1119#，桩径为1.4 m，桩长5.15 m，间距为 2.0 m；第三排桩桩号为 2029#，桩径为 1.0 m，桩长12.2 m，间距为2.0 m。第一排桩与第二排桩的间距【作者简介】尹亮洲，男，湖南衡阳人，硕士，中国电建集团中南勘测设计研究院有限公司，注册一级结构工程师，工程师，研究方向：岩土工程；汤祖平，男，湖南株洲人，硕士，任职于中交中南工程局有限公司，工程师，研究方向：岩土工程。【引用本文】尹亮洲，汤祖平.高边坡治理实体有限元数值模拟及方案比选 J.企业科技与发展，2023（7）：91-93，108.高边坡治理实体有限元数值模拟及方案比选尹亮洲1，汤祖平2（1.中国电建集团 中南勘测设计研究院有限公司，湖南 长沙 410014；2.中交中南工程局有限公司，湖南 长沙 410002）摘要：鉴于高边坡治理工程地质、地形的复杂性及施工的高风险性，文章对某高边坡工程提出3种治理方案，并利用有限元分析软件，建立3个方案的三维实体有限元模型，从位移、内力、边坡稳定性及经济性等角度综合论证3个方案的优点和缺点，从中选取最合适的治理方案。分析表明：墙顶位移主要由墙体自身刚度决定，桩基数量对水平位移影响不大；排桩方案刚度较小，容易产生过大变形；扶壁式挡墙整体刚度大，自重轻，既能减少墙顶位移，又不会使桩基产生过大的内力。将最终方案数值模拟结果与后期施工监测结果进行对比，理论结果与监测结果趋势基本一致。关键词：高边坡治理；三维实体有限元模型；治理方案；施工监测中图分类号：TU418文献标识码：A文章编号：1674-0688（2023）07-0091-04尹亮洲，汤祖平.高边坡治理实体有限元数值模拟及方案比选企业技术实践91为3.2 m，第二排桩与第三排桩的间距为2.0 m，各排桩之间呈梅花形布置。第一级边坡1 742.00平台11 746.00平台21 754.00第二排桩第一排桩第三排桩填土2第二边坡1 754.80车道荷载压力P1填土2强风化玄武岩宿舍基底压力P0图1高边坡治理施工方案示意图（单位：m）20#图2宿舍、道路及桩位平面关系由于支挡结构位于填土1上，天然地基承载力较低，下面土层为强风化玄武岩及中风化玄武岩，填土地基上存在车道荷载P1，考虑宿舍基底压力P0对支挡结构的影响，经综合分析，拟采用自重较轻的扶壁式挡墙及桩板墙进行支护，初步拟定以下3个高边坡治理方案。方案1：填土到平台1标高，施工第一、第二排桩基础，然后浇筑上部扶壁式挡墙，挡墙扶臂厚0.4 m，悬臂厚0.4 m，底板厚0.5 m，最终填土到平台2标高。桩及扶壁式挡墙均采用C30混凝土。方案2：填土到平台1标高，施工第一、第二、第三排桩基础，然后浇筑上部扶壁式挡墙，挡墙扶臂厚0.4 m，悬臂厚0.4 m，底板厚0.5 m，最终填土到平台2标高。桩及扶壁式挡墙均采用C30混凝土。方案3：填土到平台1标高后继续填土到平台2标高，施工第一、第二、第三排桩基础，然后施工顶板及侧板，顶板厚1 m，侧板厚0.2 m。桩及板均采用C30混凝土。下文通过实体有限元模型对3个方案进行分析，从内力、位移、稳定性及经济性等方面对3个方案进行比选。2实体有限元数值模拟根据上文各支护方案，采用MIDAS/GTS有限元软件对高边坡治理进行有限元模拟。其中，填土、强分化玄武岩及中风化玄武岩均采用实体单元。方案1、方案2中的扶臂挡土墙的扶臂、侧板、底板及方案3中的顶板、侧板均采用板单元，桩基础采用梁单元。实体单元、板单元及梁单元之间彼此共节点耦合。5层宿舍的基底压力P0为60 kPa3，车道为8 m宽的双车道，车道荷载P1为2.625 kPa4，分别以面荷载形式作用在实际实体单元表面上。实体有限元模型如图3所示；各方案结构部分如图4所示。图3三维实体有限元模型（a）方案1（b）方案2（c）方案3图43个高边坡治理方案的结构图土体本构模型采用Mohr-Coulomb模型。根据岩土工程地质勘察资料，结合周边桩基施工监测到的桩身变形结果，对土体参数进行多次试算与反演，确定符合本工程的合理的主要岩土参数（见表1）。表1岩土力学参数表岩土类别填土碎石垫层强分化玄武岩中分化玄武岩弹性模量E（MPa）3201880重度（kN.m-3）17192428泊松比0.30.20.250.2粘聚力C（kPa）8060240内摩擦角（）2.5354050企业科技与发展，2023年，第7期，总第501期923方案计算结果对比3.1墙体位移经过计算，方案1的墙体最大位移为墙底部水平向位移，大小为1.47 cm。方案2的墙体最大位移为墙顶部水平向位移，大小为1.53 cm。方案3的墙体最大位移为侧板顶部水平向位移，大小为4.00 cm。分析表明：3个方案在平台1以下入土部分桩体的水平位移均较小，位移主要由上部结构的刚度决定。方案1、方案2的上部扶壁式挡土墙结构整体刚度比方案3的抗滑桩大，前二者位移相对较小。3.2墙体内力各方案墙体弯矩如图5所示。方案1、方案2的最大弯矩发生在墙趾、墙踵与扶臂交汇处，大小分别为-726358 kNm、-294123 kNm。方案3的最大弯矩发生在顶板中部和角部，大小为-892275 kNm。分析表明，方案1的最大墙身弯矩约为方案2的2.5倍。方案3虽然桩数较多，但是整体刚度小，顶板弯矩较大，最大弯矩约为方案2的3倍。各方案墙体剪力结果表明：方案1、方案2的最大剪力发生在墙踵与扶臂交汇处，大小分别为-1 1781 201 kNm、-838870 kNm。方案3中除个别边角处有应力集中外，剪力比较平均，大小为-781851 kN。分析表明：方案1的最大墙身剪力约为方案2的1.4倍，方案3与方案2的剪力大小相当。3.3桩体内力对3个方案的第一、第二排桩各方案的桩身内力及桩顶位移进行分析，得出以下结论：方案1承受较大的剪力和弯矩；方案2承受弯矩明显小于方案1；方案3承受的内力较小，但产生的位移较大。对3个方案第三排桩各方案的桩身内力及桩顶位移进行分析，得出以下结论：方案2承受较大的内力，但位移较小；方案 3 承受的内力较小，但产生的位移较大（见表2）。4方案选择表2是以方案2为参照的对比数据与结果，表2中的S、N、V、M分别表示桩顶位移、桩轴力、桩剪力、桩弯矩，字母的数字编号表示桩排号，例如S1表示第一排桩的桩顶位移。由上文方案计算结果及表2对比结果可知，方案3的墙顶位移为4 cm，变形过大，不满足工程要求，不予考虑。方案1的桩身内力偏大，对桩身承载力要求高，但桩数量少，可节约造价及工期。方案2的位移较小，同时桩身内力更均匀，桩数有所增加，但结构更安全，经济性满足要求，经反复论证，作为最终方案。表23个高边坡治理方案的对比结果（以方案2为参照）对比内容桩数（根）挡土结构墙体最大位移倍数墙身最大弯矩倍数桩顶位移（S）倍数桩轴力（N）倍数桩剪力（V）倍数桩弯矩（M）倍数抗倾覆稳定安全系数优点缺点方案119扶壁式挡墙0.962.51.19（S1）、1.19（S2）1.05（N1）、3.33（N2）1.57（V1）、2.33（V2）1.99（M1）、9.51（M2）2.56桩数量少，节约造价及工期；上部结构刚度大，水平位移小桩身剪力和弯矩较大方案229扶壁式挡墙0000002.98上部结构刚度大，水平位移小；桩体以受轴力为主，弯矩和剪力较小，受力性能好桩数量多，费用稍高方案329顶板、侧板体系2.6133.25（S1）、3.25（S2）、3.25（S3）0.71（N1）、1.50（N2）、0.65（N3）0.77（V1）、0.62（V2）、0.30（V3）1.15（M1）、4.37（M2）、0.34（M3）2.11一桩到顶，施工工艺简单，易操作变形过大，桩数量多，费用稍高，工期较长（a）方案1（b）方案2（c）方案3图53个高边坡治理方案的墙体弯矩图（下转第108页）尹亮洲，汤祖平.高边坡治理实体有限元数值模拟及方案比选93（见表4），可以看出整体偏差较小，其中项目总投资的误差控制在0.02%以内，在项目的允许偏差范围内。表4项目概算与估算偏差序号一二三四五工程或费用名称工程费用工程建设其他费用基本预备费运营期费用总投资估算（万元）120 096.769 131.656 461.422 484.02138 173.85概算（万元）120 094.719 111.816 460.322 484.02138 150.86偏差（%）-0.001 7-0.22-0.020-0.02可见，通过优化建设项目内容和结构，可提升项目估算的准确性，减小项目后期变更所带来的影响；建立人力、材料、机械设备等所需费用数据库，有利于资源的利用和分配，提高项目资金的利用率；通过划分项目阶段和任务，可提高项目计划和管理效率，提升项目风险管控的能力。3结语在项目估算过程中，通过对项目建设内容进行结构优化，一方面使能够防止估算过程中发生缺项、漏项，另一方面能使项目内容得到深层次的优化，使项目的整体框架更加清晰。此外，通过建立当地人力、材料、机械设备等所需费用机数据库，不仅能够及时调整项目的单价，还为后续开展评审和评估工作提供借鉴指导。根据项目的整体情况，划分项目阶段和任务，有效地避免各单位间因信息差造成的工程进度缓慢、资金浪费等问题，通过项目后期概算误差控制在0.02%以内，表明文章提出的管控策略可有效提升项目估算精确度。4参考文献1 刘少浪.设计院主导的EPC模式下工期与成本关系的研究与优化 D.广州：华南理工大学，2020.2 范伟.“EPC+O”总承包模式下政府监管的研究以江门市蓬江区为例 J.经济研究导刊，2021（30）：120-122.3 刘运鹄.基于IABC-LS-SVM的水利工程项目建设成本估算与控制方法 D.大连：大连理工大学，2021.4 唐光，洪朝霞.生态水利工程概（估）算编制难点及对策分析以某生态水利项目为例 J.工程造价管理，2023，（2）：36-43.5 郭创.工程“EPC+O”实施模式要点探讨基于流域水环境综合整治工程实践的分析 J.重庆建筑，2021，20（1）：18-20.6 何凯红，张伟，祁玉婷.商业建筑采用“EPC+O”模式的招标管控要点 J.建筑，2022，968（24）：32-34.5模拟数值与监测数据对比施工过程中对挡土墙顶水平及竖向位移进行监测，监测点为 A、B、C 3 个点，分别位于平面上10#、5#、1#桩对应的墙顶位置。将监测值与方案2的数值计算理论值进行对比，结果如图6所示。由图6可知，墙顶竖向位移中间大、两端小。由于荷载不对称，A点水平位移大于其余2点。监测值与理论值发展趋势一致，但水平位移和竖向位移均比理论值小。ABC-0.020-0.015-0.010-0.0050.0000.0050.0100.015位移（m）挡土墙顶位移监测点 水平位移理论值 水平位移监测值 竖向位移理论值 竖向位移监测值0.0150.0100.0050-0.005-0.010-0.015-0.020挡土墙顶位移监测点BAC水平位移理论值水平位移监测值竖向位移理论值坚向位移监测值位移（m）图6位移理论值与实测值对比6结论本文以某边坡治理为研究案例，结合项目实际情况，提出3种边坡治理方案。经过对比分析，得到如下结论：墙顶位移主要由墙身刚度决定，不同数量桩基对水平位移影响不大。排桩方案由于自身刚度小，在高边坡情况下容易产生过大变形，本工程高边坡临近学生宿舍，对变形控制要求高，故此方案不适用于本工程。扶壁式挡墙整体刚度大，自重轻，既能减少墙顶位移，又不会使桩基产生过大的内力。方案2上部结构采用扶臂式挡墙，基础采用3排桩基础，内力和位移都比较小，抗倾覆稳定性较好，经济安全，因此作为最终方案。7参考文献1 王恭先.滑坡防治方案的选择与优化 J.岩石力学与工程学报，2006，25（S2）：3867-3873.2 石广斌，杨经会，安盛勋，等.石泉水电站陡高边坡崩塌机理分析及加固设计 J.岩石力学与工程学报，2004（7）：1186-1192.3 肖世国.边（滑）坡治理中h型组合抗滑桩的分析方法及工程应用 J.岩土力学，2010（7）：2146-2152.4 许东俊，李小春，余毓良，等.大型露天矿山不稳定边坡治理初探 J.岩土力学，1993（1）：33-39.（上接第93页）企业科技与发展，2023年，第7期，总第501期108