路堑边坡传递系数法与改进

王宗林，梁天意，梁 冰

(辽宁工程技术大学矿业学院，辽宁阜新123000)

1 传递系数法

传递系数法的全名为不平衡推力传递系数法，在公路路堑边坡及其他类型边坡的稳定性分析上应用比较广泛，是基于有限平衡原理的稳定系数定量求解公式。

传递系数法是针对折线型移动破坏面条件提出的，假定块条间力方向与上一块条底滑面平行，根据力的平衡，逐条向下推算出上一块条对下一块条的推力，同时不断设定系数大小，使最后一条没有向下的推力，由此算出边坡的稳定性系数［1］。一般常用的传递系数法有强度储存法和超载法。

1.1 强度储存法

传递系数法假定:①边坡稳定属平面应变问题，即可取其某一横剖面为代表进行分析计算;②滑动面为折线形，即在横剖面上滑动面为直线;③整个滑动面上的稳定安全系数是一样的。

取坡体的其中一个块条作为研究对象 (见图1)，根据竖向与横向力的平衡条件得［2］:

将公式 (3)代入公式 (1)、(2)中，利用三角函数，并消去Ti、Ni得:

图1 传递系数法求解示意Fig.1 Transfer coefficient method calculation diagram

式中:Wi——第i块坡条重度，kN/m3;Pi——第i+1块坡条传上来的力，kN;Pi-1——第i-1块坡条传上来的力，kN;αi——第i块坡条移滑面与水平面的角度，(°);αi-1——第i-1块坡条移滑面与水平面的角度，(°);ci——第i块坡条的粘聚力，kPa;φi——第i块坡条的内摩擦角，(°);li——第i块坡条移滑面长度，m;c'i——强度折减后第i块坡条的粘聚力，kPa;φ'i——强度折减后第 i块坡条的内摩擦角，(°);Fs——坡体稳定性系数;φi-1——传递系数。

利用公式 (5)，首先假定一个稳定系数Fs，从边坡最上面的块条开始计算它的不平衡推力，之后依次向下进行计算，直到最后一块，如果最下面的边坡块条的推力求解后为0，则之前假定的Fs即为所求的稳定性系数。

1.2 超载法

如果依托不平衡下滑力的定义，即不平衡下滑力等于下滑力乘以Fs减去抗滑力，可得:

求解Fs的条件是Pi=0，则可得出一个Fs的一次函数。整个块条计算推力的过程中，之间的力如果出现方向向上，即为拉力，而这是不符合条件的力，因此下一块条的上推力就取Pi=0。

1.3 超载法的累积式

通过对超载法的公式逐步递推，可以得到累积式:

2 改进的传递系数法

传递系数法考虑了块条界面上剪应力的作用，计算简单且应用性较好。但是传递系数法中Pi的方向被假定为与上邻块条的底移滑面方向相同，当滑面较陡与较缓或考虑地下水作用时，条间力方向不能保证总是与底移滑面方向相同，导致计算有所偏差。

为此，本文对传统的传递系数法进行修改，不硬性规定块条之间相互作用力的方向，从强度储存法及超载法两个方面推导出新的传递系数法公式，同时将坡体受到水作用后的水压力［3］考虑进去，得出水作用下的传递系数法。

2.1 改进的传递系数法块条力学模型

传递系数法是在一定前提下将移滑面定为折线型的滑坡推算出的，改进的传递系数法仍延用该假定前提。首先将滑面分成n个块条，对其中第i个块条进行分析，将块条间的相互作用力的方向由以前的平行于底滑面改为与底滑面的角度为βi;同时由于原有传递系数法的块条力学模型中没有充分考虑外界力，因此本文在原有的块条力学模型上增设2个力，分别为坡条上的外界力Qi及水平地震力KsWi，并建立新的块条力学模型 (见图2)。新建的块条力学模型中Pi的方向与水平方向的角度为βi，那么Pi-1与水平方向的角度为βi-1。

图2 新传递系数法的计算示意图Fig.2 New transfer coefficient method calculation diagram

2.2 改进的传递系数法求解公式推导

根据新建块条的力学模型，通过水平方向与竖直方向的力平衡可得到:

将公式 (3)代入公式 (9)，(10)中，利用三角函数，并消去Ti、Ni得:

式中:Qi——坡条上的外界力，kN;δi——坡条上的外界力与竖直方向的角度，(°);Ks——地震作用系数，KsWi——平行地震力，kN;βi——Pi与水平方向的夹角，(°);βi-1——Pi-1与水平方向的夹角，(°)。

对于βi的取值，一般认为:

①如果块条i底滑面的角度αi≤45°，且相邻块条的角度没有太大变化，则块条间的作用力方向与折减后块条的内摩擦角及块条i底移滑面的角度的关系为

②如果块条i底滑面的角度αi＞45°，或者相邻块条的角度有大的变化，则块条间的作用力方向与折减后块条的内摩擦角及块条i底滑面的角度的关系为。

3 水力作用下的改进传递系数法

水对边坡的破坏性作用很大，水不仅能够通过增加坡体的重度，加大下滑力，而且水渗入岩石之后，无论是滞留在孔缝、断层之中，还是沿着孔隙继续扩大运动，都将让原有的孔隙水压力增大或者从无到有。这种孔隙水压力将对坡体产生力学、化学等多方面的作用。

静止的水对坡体周围的作用力称为静水压力，而流动的水对坡体中的骨架有渗流力，把这部分水对单位体积骨架的作用力称为动水压力，动水压力也称为体积力［4］。

3.1 静水压力下的传递系数法

同样取假定滑面的任意一块条，对其进行分析。由于水的作用，块条被分为两部分，上面为正常的块条，下面为有水作用的块条 (见图3)。

利用前面推导出的公式，可以得出:

将公式 (15)带入公式 (13)，(14)中，利用三角函数，并消去Ti、Ni得:

图3 静水压力作用Fig.3 Hydrostatic pressure effect

式中:W1i——第i块坡条中水润线以上块条的重度，kN/m3;W2i——第i块坡条中水润线以下块条的重度，kN/m3;Pai——第i块坡条左边界静水压力，方向与块条面垂直，kN;Pbi——第i块坡条右边界静水压力，方向与块条面垂直，kN;Ui——第i块坡条底移动破坏面静水压力，kN。

静水压力的计算如图4所示。根据流线与等势线垂直的流网性质来确定条块周边的静水压力。过D点和A点分别做EC的垂直线GD和EA，再过G点和E点分别做垂线GH、EF垂直于CD和BA。这样得到A、D点的水头高度分别为AF和DH，所以:

则左边和右边静水压力的总力为:

式中:ηi——水润线与水平线之间的角度，(°);ha——第i块坡条左边侵润线到底滑面的距离，m;hb——第i块坡条右边侵润线到底滑面的距离，m;γw——水的容重，kN/m3。

3.2 动水压力下的传递系数法

同样取假定滑面的任意一块条，对其进行分析。由于水的原因，块条被分为两部分，上

图4 静水压力计算示意Fig.4 Hydrostatic pressure calculation diagram

面为正常的块条，下面为有水作用的块条 (见图5)。

图5 动水压力作用Fig.5 Hydrodynamic pressure effect

同理可得:

将公式 (15)带入公式 (23)，(24)中，利用三角函数，并消去Ti、Ni得:

式中:Di——第i块坡条的动水压力，kN/m3。

4 改进型传递系数法的可靠性

为了验证改进型传递系数法的实际可行性，本文利用改进的传递系数法对假定边坡进行了稳定性系数的计算。对假定边坡的滑面进行分条，总共为10个部分 (见图6)，坡体的岩体力学参数见表1。

图6 边坡滑动面划分示意图Fig.6 The slope slip plane division schematic diagram

表1 岩体物理力学参数Table 1 Physical and mechanical parameters of rocks

目前的边坡稳定性系数计算主要采用的是传递系数法，bishop法、FLAC数值模拟等方法。其中传统的传递系数法计算的稳定系数偏高，偏高值随φ的增大而增加，一般情况下，偏高10%～15%以内。而简化的bishop法则考虑了条间推力，使受力更合理，计算结果更接近实际。本文对传统传递系数法进行改进，把水的作用考虑进去，在块体受力分析上更加全面［5～11］。

本文运用传递系数法、改进的传递系数法、bishop法、FLAC数值模拟方法对此工程进行了稳定性系数的计算，计算结果见表2。在计算的数值上，改进的传递系数法与其他计算方法的结果比较贴近，并且比传统的方法更精确。因此，改进后的传递系数法有着一定的精度及可靠性。

表2 不同方法下的稳定系数值Table 2 Stability coefficient values of different methods

5 结论

简述了传统传递系数法两种解法的基本原理及计算公式，同时也说明了传统传递系数法假定滑坡块条间相互作用力方向与上一块条底滑面方向相同的不足。针对此不足提出了改进的传递系数法，并推导出改进的传递系数法计算公式，同时考虑了在水作用情况下的传递系数法，并推导出静水和动水压力下传递系数法的计算公式，最后通过改进的传递系数法对假定边坡稳定性系数进行计算，与其他稳定性系数计算法的计算结果进行对比，可以得出改进的传递系数法有一定的精度及可靠性。改进的传递系数法具有在实际生产过程中应用的条件。

［1］ 陈祖义.岩质边坡稳定性分析——原理方法程序［M］.北京:中国水利水电出版社，2005.CHEN Zu-yi.Stability analysis of rock slope:Principle，method and procedure［M］.Beijing:China Water＆ Power Press，2005.

［2］ 景锋.岩质高边坡稳定性研究［D］.武汉:武汉大学，2008.JING Feng.Stability study on rock high slope［D］.Wuhan:Wuhan University，2008.

［3］ 冯文娟，李建林，杨学堂.岩土体边坡稳定性分析计算方法的研究进展［J］.黑龙江水专学报，2005，32(2):10～13.FENG Wen-juan，LI Jian-lin，YANG Xue-tang.Study on the calculation methods of the slope stability analysis［J］.Journal of Heilongjiang Hydraulic Engineering College，2005，32(2):10～13.

［4］ 黄志斌.濒危岩质高陡边坡稳定性及处置技术研究［D］.长沙:中南大学，2007.HUANG Zhi-bin.Study on stability of dangerous rock high slope and processing technology［D］.Changsha:Central South University，2007.

［5］ 赵尚毅，郑颖人，时卫民，等.用有限元强度折减法求边坡稳定安全系数［J］.岩土工程学报，2002，24(3):343～346.ZHAO Shang-yi，ZHENG Ying-ren，SHI Wei-min，et al.Analysis on safety factor of slope by strength reduction FEM ［J］.Chinese Journal of Geotechnical Engineering，2002，24(3):343～346.

［6］ 王肇慧，肖盛燮，刘文方.边坡稳定性计算方法的对比分析 ［J］.重庆交通学院学报，2005，24(6):99～103.WANG Zhao-hui，XIAO Sheng-xie，LIU Wen-fang.Comparison of analytical methods for slope stability ［J］.Journal of Chongqing Jiaotong University，2005，24(6):99～103.

［7］ 朱大勇，李焯芬，黄茂松，等.对3种著名边坡稳定性计算方法的改进 ［J］.岩石力学与工程学报，2005，24(2):183～194.ZHU Da-yong，LI Chao-fen，HUANG Mao-song，et al.Modifications to three well-known methods of slope stability analysis［J］.Chinese Journal of Rock Mechanics and Engineering，2005，24(2):183 ～194.

［8］ 刘才华，陈从新，冯夏庭，等.地下水对库岸边坡稳定性的影响［J］.岩土力学，2005，26(3):419～422.LIU Cai-hua，CHEN Cong-xin，FENG Xia-ting，et al.Effect of groundwater on stability of slopes at reservoir bank ［J］.Rock and Soil Mechanics，2005，26(3):419～422.

［9］ 邓科.边坡稳定性计算方法的适用条件及分析比较［J］.灾害与防治工程，2005，(2):12～16.DENG Ke.Comparison of analytical methods for slope stability ［J］.Disaster and Control Engineering，2005，(2):12～16.

［10］ 李领先.基于极限平衡法的边坡稳定性计算 ［J］.科技信息，2009，(12):258～259.LI Ling-xian.Slope stability calculation based on limit equilibrium method ［J］.Science and Technology Information，2009，(12):258～259.

［11］ 邓海政.公路高边坡稳定性计算方法的对比分析［J］.黑龙江科技信息，2008，(20):222.DENG Hai-zheng.Comparison analysis of calculating methods for high slope stability along highway ［J］.Heilongjiang Science and Technology Information，2008，(20):222.