边坡稳定分析Sarma法中条块间作用力方向的两种可能性.docVIP

论边坡稳定分析Sarma法中条块间作用力方向的两种可能性 提 要在各种边坡稳定分析方法中，Sarma法使用了具有倾斜边界的土条，假定侧面上也达到极限平衡。研究发现，在这些倾斜边界上，存在着两种条块间相对滑动方向，即左条块相对右条块向上或向下滑动，如果对这两种可能性不予以区分，则可能导致Drucker准则被破坏。算例表明，只考虑一个剪切方向的计算结果有时不合理，会导致侧面的凝聚力越大，安全系数越小的反常情况。 关键词 边坡稳定分析；极限平衡；Drucker准则。 本文于1996年11月7日收到，系国家自然科学基金资助重大项目。 1 前言 传统的边坡稳定极限平衡分析法采用垂直条分法，这个方法不能很好地考虑条块侧面的力的特性，特别是岩质边坡的断层节理特征。Sarma(1979)提出对滑坡体进行斜分条的极限平衡分析法［1］，该法假定沿条块侧面也达到了极限平衡，这样，通过静力平衡条件即可唯一地确定边坡的安全系数或加载系数。这个方法受到Hoek教授的极力推崇(Hoek,1983)［2］。在开发这一方法过程中，发现沿斜条块的摩擦方向有两种可能性，即左条块相对右条块可能向上或向下运动。Sarma法只考虑了常见的一种，即左条块相对右条块向上运动。在某些情况下，仅考虑这一种情况可能导致错误的计算成果，本文对这一问题作一简要的讨论。 图1 Sarma法和能量法的计算力 Fig.1 Sarmas method and the energy method 2 求解安全系数的两种途径 对于一个多块体组成的破坏模式，可以通过求解静力平衡方程计算安全系数，也可以通过建立虚功原理的方程求解，现以图1所示两个块体的例子说明此问题。 分析作用于左、右两个处于极限平衡的块体的力，根据摩尔库仑破坏准则，作用在滑面上的力可分成两部分，一部分是由法向N和由它决定的摩擦阻力Ntgφe构成的一个“摩擦反力”P，它与滑面的法线夹角为φe，另有凝聚力ce提供的“凝聚阻力”C，其值为ceL,L为滑面长度。其中安全系数隐含于对强度指标修正中，即tanφe=tanφ/F,ce=c/F。同样，在条块的界面AB上右条块作用在左条块的力为Pj和Cj，左、右条块的静力平衡方程分别为（式（1）至式（7）的推倒过程，均指向量运算）： Wl+Pl+Cl+Pj+Cj=0 (1) Wr+Pr+Cr-Pj-Cj=0 (2) 其中W为条块重量，上述两个矢量方程如在x，y方向投影，可建立四个方程式，其中的未知数分别为Pl,Pr,Pj的大小(其方向已知，与法线方向夹φe角)，另外一个未知数是包含在ce和φe中的F值，这样，就可以求得F值。 采用虚功原理求解，则我们分别在左、右条块赋一个虚位移Vl和Vr，有 Wl·Vl+Pl·Vl+Cl·Vl+Pj·Vl+Cj·Vl=0 （3） Wr·Vr+Pr·Vr+Cr·Vr-Pj·Vr-Cj·Vr=0 （4） 两式相加得 Wl·Vl+Wr·Vr+Pl·Vl+Pr·Vr+Cl·Vl+Cr·Vr+Pj·Vj+Cj·Vj=0 （5） 其中：??? Vj=Vl-Vr （6） 为左条块相对右条块移动的位移。虚功原理对任意一组协调的位移场都适用。式(5)说明，对滑动体加一组虚位移时，重力所作的功等于沿左、右滑面和界面的内能耗散。 现在施加一组特殊的位移Vl、Vr和Vj，它们分别与滑面的夹角为φe，见图1，由于Pl、Pr和Pj与相应滑面的法线夹角为φe，故Pl和Vl，Pr和Vr，和Pj和Vj正交，它们的积应该分别为零，式(5)简化为 Wl·Vl+Wr·Vr+Cl·Vl+Cr·Vr+Cj·Vj=0 （7） 使用位移协调条件(6)式，则一旦Vl确定，Vr和Vj即可表达为Vl的线性函数。式(6)代入式(7)后Vl被消去，式(7)中仅剩下隐含在ce,φe中的F值。因此，根据虚功原理建立的式(7)使求解边坡稳定安全系数成为通过一个方程解一个未知数的十分简捷的方法。它的解题步骤是先计算一个与各条块底滑面夹角为φe的协调的位移场，然后根据外力功等于内能耗散的虚功方程式(7)求解F值。使用这一方法带来的另一个好处是求解过程中存在着坡体的位移图式，这正好为论述本文提出的命题提供了理论基础。 3 边坡稳定分析的能量法 近年来Donald和Chen提出了的求解边坡稳定分析的能量法［3］［4］就是建立在上节介绍的虚功原理方法上的。对于相邻条块，作者推导出由左侧条块速度确定右侧条块速度的方程(图2a)。 Vr=Vlsin(θl-θj)/sin(θr-θj) （8） Vj=Vlsin(θr-θl)/sin(θr-θj) （9） （10） 其中D为内能耗散，上标j和s分别代表界面和滑面，其数值就是式(7)左边的第3，4，5项的负值。 从现在开始，将“位移”改称为“速度”，即单位外力增量作用下的塑性位移。 在Do