1. 引言
边坡稳定分析和加固处理一直是岩土工程领域内的热点,边坡稳定性的分析方法有很多,目前边坡稳定分析中常用的方法是极限平衡法和抗剪强度折减法。自从瑞典学者彼德森提出极限平衡法之后,很多学者在极限平衡法分析边坡稳定上展开了更深入的研究 [1]。周资斌 [2] 和苏利军 [3] 用极限平衡法对实际工程做了扩展性分析。卢坤林 [4] 等将离散后的条柱间作用力等效为滑面上的正应力,由平衡条件提出了一种适应空间形态的三维极限平衡法。邓东平 [5] 等在三维模型中,根据采用极限平衡方程的数目不同,提出了三维边坡极限平衡严格法、准严格法和非严格法。白桃 [6] 等比较了四种不同算法下的安全系数,提出只有同时采用安全系数和破坏概率两个指标对边坡安全评价才会更加有效。还有一些学者通过有限元软件分析了抗滑桩的桩长、桩径和抗滑桩的加固位置对边坡的临界滑动面和安全系数的影响 [7] [8] [9]。
以上学者仅通过理论分析或者数值模拟的方法对边坡的稳定性分析进行了研究,本文采用极限平衡理论和有限元强度折减理论研究了柬埔寨55号公路边坡采取不同方案下的加固效果,对分析了锚索加固和抗滑桩加固下的锚固深度加固位置的影响。
2. 工程地质概况
本文分析的边坡位于柬埔寨55号公路K145 + 660~K145 + 880路段路基左侧,该段山体的自然坡度约43˚,斜坡表层有大量的白垩纪强风化砂质泥岩少部分有块状砂岩。根据地质勘察报告,边坡顶部约1.7 m为坡积粉质黏土和全风化砂质泥岩,其下为强风化砂质泥岩,再向下为强风化砂岩,最下面为中风化砂质泥岩,以泥质结构为主,岩层结构面逆坡体方向倾斜,倾角约10˚,无外倾结构面,该边坡为类土质边坡。边坡开挖深度60.5 m,坡体中下部为中风化砂质泥岩,坡体节理裂隙发育,容易产生由于坡脚应力不足的坡脚压碎变形破坏。边坡原设计为6级,最大挖深为60 m,坡率为1:0.5、1:0.75、1:1、1:1.25、1:1.25、1:1.5,每级边坡高度为10 m,边坡间设置2 m宽平台。
3. 边坡稳定分析
对于一复杂的边坡分析中,一般情况下会取该边坡的剖面对其简化分析,常用的分析方法有极限平衡法和有限元强度折减法,本文将会采用这两种方法计算边坡的安全系数,分析边坡的稳定性。
3.1. 边坡稳定分析原理
3.1.1. 极限平衡理论
极限平衡理论首先在边坡土体中假设一破坏面,破坏面内的土体作为一脱离体计算,当脱离体达到静力平衡状态时脱离体所需要的抗力或者抗剪强度,与破坏面实际所能提供的抗力或者抗剪强度的比值,为边坡或者岩体的稳定安全系数。
极限平衡法可分为严格条分法和非严格条分法,非严格条分法满足力平衡或者力矩平衡,而且在计算时主要采用标准圆弧计算,但是在实际的滑坡中很少有圆弧形,所以非严格条分法在应用时有一定的不足。严格条分法不仅满足力平衡而且满足力矩平衡,而且在应用时严格条分法可以对任意的滑动形状进行求解。进过半个世纪的发展,不少学者提出了自己的方法,比如Morgenstern和Price (1965)提出了Morgenstern-Price法 [10],Chough (1986)提出了一种边坡稳定分析的通用的条分法(GLE) [11] [12],这些方法为边坡稳定安全系数的求解提供了丰富的理论。
3.1.2. 有限元强度折减理论
强度折减法的基本原理是将土体抗剪强度指标c、φ用强度折减系数Fr进行折减,分别得到cm、φr,cm、φr由公式(1)和(2)求得,当边坡土体达到临界破坏状态时所对应的强度折减系数Fr即为对应条件下的护坡最小稳定安全系数。
(1)
(1)
强度折减系数定义为:在外荷载保持不变的情况下,边坡坡体所能发挥的最大抗剪切强度与外荷载在边坡内所产生的实际剪应力之比。当假定边坡内所有坡体抗剪强度的发挥程度相同时,这种抗剪强度折减系数定义为边坡的整体稳定系数。
3.2. 安全系数和临界滑动面计算分析
本文使用通用条分法和有限元强度折减法进行计算安全系数和临界滑动面,根据地质勘察报告提供的数据,计算参数选用表1。计算结果参看表2,计算边坡滑动面参看图1。
![](Images/Table_Tmp.jpg)
Table 1. Calculation parameter table
表1. 计算参数表
![](Images/Table_Tmp.jpg)
Table 2. Comparison table of safety factor calculation
表2. 安全系数计算对比表
![](//html.hanspub.org/file/4-2751194x12_hanspub.png)
Figure 1. Comparison diagram of sliding surface between GLEmethod and strength reduction method for 2 m platform slope
图1. 2 m平台边坡GLE法和强度折减法的滑动面对比示意图
通过表2安全系数计算结果对比可知,有限元强度折减法计算出的结果要比GLE法中圆弧滑动面和折线滑动面都要小,探究其原因可能是在模型计算中网格划分时不够合理。而在GLE法中不同滑动面的计算结果对比可以看出,折线滑动面计算出的结果要比圆弧滑动面计算的要小,其主要原因在于折线滑动面能够搜索到更过的滑动面,而且搜索到的滑动面更接近滑坡真实滑动面。在图1中可以看到GLE强度折减法的滑动面要比另外两种方法的滑动面都要大。
![](//html.hanspub.org/file/4-2751194x13_hanspub.png)
Figure 2. Comparison diagram of sliding surface between GLE method and strength reduction method for 5 m platform slope
图2. 5 m平台边坡GLE法和强度折减法的滑动面对比示意图
三种方法计算出的安全系数相差不大,但是使用强度折减法计算出的安全系数小于1.25,这不满足《公路路基设计规范》中相关要求,而且用极限平衡法算出的安全系数也仅仅大于1.25,因此这样的边坡既没有足够的安全储备,也不满足规范要求。所以需要对边坡重新优化设计。
结合项目的实际特点,公路中存在大量填方路段,若增加开挖边坡上部土方,这样不仅能够减少上部荷载,而且能够减少工程造价。因此将第三级边坡与第四级边坡之间的平台调整为5 m,再次进行安全系数计算和稳定分析。计算结果见表3,分析滑动面见图2。
![](Images/Table_Tmp.jpg)
Table 3. Comparison table of safety factor calculation
表3. 安全系数计算对比表
由表3可知,当选用5 m宽平台时,极限平衡法与强度折减法算得的边坡稳定安全系数,均大于1.25,满足《公路路基设计规范》中的要求。
4. 边坡防护设计分析
边坡防护方案多种多样,在设计防护方案类型时,首先要保证边坡的稳定性,其次还要考虑经济、美观、适用等因素,按照边坡的水文、地质、岩质等情况,选择适当的防护方案。通常主要的防护方案是锚杆(锚索)加固、桩板墙及抗滑桩加固、锚喷加固等,再配合骨架防护、挡土墙、植被防护等,达到加固边坡,避免失稳破坏。
本文将选择锚杆加固和抗滑桩加固方案进行分析。
4.1. 锚杆加固边坡稳定性分析
锚杆加固边坡稳定性分析使用ABAQUS软件进行有限元强度折减法分析计算,锚杆采用非预应力锚杆,采用杆单元为T2D2,外径25 mm,抗拉强度为80 KN,锚杆逆坡体方向倾斜,倾角约20˚,锚固于二、三、四级边坡,锚杆与土的接触方式选择锚杆嵌入土中。
分别将锚杆嵌入素土边坡临界滑动面以下不同深度,对比分析锚固段不同深度对于边坡加固效果的影响。通过ABAQUS软件分析计算得到锚固深度为2 m时的边坡位移云图(图3)和其他不同锚固深度时的不利滑动面(图4),以及不同锚固深度下的安全系数(图5)。
![](//html.hanspub.org/file/4-2751194x14_hanspub.png)
Figure 3. Displacement nephogram of finite element isoline for 2 m anchored slope stability
图3. 2 m锚固边坡稳定性有限元等值线位移云图
![](//html.hanspub.org/file/4-2751194x15_hanspub.png)
Figure 4. Schematic diagram of the most unfavorable sliding surface of slope reinforced with different anchorage depths
图4. 不同锚固深度加固边坡最不利滑动面示意图
![](//html.hanspub.org/file/4-2751194x16_hanspub.png)
Figure 5. Diagram of variation of safety factor with anchorage depth
图5. 安全系数随锚固深度变化示意图
由图3可以看出,随着锚杆嵌入深度的增加,边坡最不利滑动面也随着向坡面逆向偏移。安全系数也随之增大,边坡的稳定性逐渐增大。在图4中开看出,在锚杆嵌入稳定岩层0~4 m的深度内,安全系数增加的速率是最快的。在4~16 m内的增速放缓。考虑整体节约经济的角度,选取锚杆长度时,可选取嵌入稳定岩层内4~8 m内的锚杆比较合适。
4.2. 抗滑桩加固边坡稳定性分析
使用ABAQUS软件分析不同深度下的抗滑桩的稳定性,抗滑桩视为理想弹性体,桩土间采用接触相互作用,接触方式选取硬接触,桩径为0.8 m,摩擦系数为0.3,采用CPE4单元。设置位置在坡高二分之一处,分不同的锚固深度,计算并提取各锚固深度模型的最不利滑动面与安全系数,计算结果如下。
![](//html.hanspub.org/file/4-2751194x17_hanspub.png)
Figure 6. Schematic diagram of slope sliding surface strengthened by anti-slide pile with different anchorage depth
图6. 抗滑桩不同锚固深度加固边坡滑动面示
![](//html.hanspub.org/file/4-2751194x18_hanspub.png)
Figure 7. Comparison of safety factors of reinforced slopes with different anchorage depths
图7. 抗滑桩不同锚固深度加固边坡滑动面示
在图6~7中对比可以看出,随着锚固深度的增加,抗滑桩边坡的最不利滑动面会向边坡逆方向移动,但滑动面比较接近。在0~3 m内,安全系数会随着抗滑桩的锚固深度而增加巨大,在3~6 m内,安全系数增加比较缓慢。由此可以看出,抗滑桩的锚固深度在2~4 m内比较合适。
抗滑桩的深度对边坡的加固效果并不是唯一影响,需要考虑抗滑桩在不同坡高处对加固效果的影响。在锚固深度选择4 m的情况下,研究在不同坡高设置抗滑桩对边坡加固效果的影响。在坡高
、
、
位置处设置抗滑桩,计算安全系数和边坡临界滑动面如图8和图9所示。
![](//html.hanspub.org/file/4-2751194x22_hanspub.png)
Figure 8. Schematic diagram of critical sliding surface of slope at different reinforcement positions of anti-slide pile
图8. 抗滑桩不同加固位置的边坡临界滑动面示意图
![](//html.hanspub.org/file/4-2751194x23_hanspub.png)
Figure 9. Comparison of slope safety factors at different reinforcement positions of anti-slide piles
图9. 抗滑桩不同加固位置的边坡安全系数对比图
观测抗滑桩在不同位置处的边坡稳定分析计算可知,随着抗滑桩加固位置的提高,安全系数随之降低,但三种加固位置处的边坡临界滑动面几乎是重合的。这表明在
处用抗滑桩加固的效果是最优的。
5. 结论
本文以柬埔寨55号公路边坡为依托工程,使用有限元强度折减法和极限平衡法对边坡进行了稳定性分析和计算。计算结果表明:
1) 安全系数会随着锚杆嵌入深度的增大而增加,而且在嵌入深度4 m前增加速率比较快,4 m之后增加缓慢。
2) 加固后的临界滑动面会随着嵌入深度的增加而逆向坡面扩展。
3) 若采用抗滑桩,在抗滑桩锚固深度为0~3 m内,安全系数会随着抗滑桩的锚固深度急剧增加,在锚固深度达到3 m后安全系数增速放缓。
4) 抗滑桩的加固位置不同,同样也会引起加固效果的变化,抗滑桩加固的位置越低,承受上部荷载越多,安全系数就大加固效果越好。
基金项目
陕西省自然科学基础研究计划,2019JM-100。
参考文献