毕肖普法土坡稳定的程序计算法

1、毕肖普法土坡稳定的程序计算法姓名：翟慧君学号：63085217007毕肖普法是由毕肖普(A. W. Bish op, 19 55 )提出的进行土坡稳定分析的一种方法。我们知 道瑞典条分法在进行土坡稳定分析的时候，不考虑相邻土条间作用力的相互影响。一般说这 样得到的稳定安全系数可能偏低10%20% ,而且这种误差随着滑弧圆心角和孔隙水应力 的增大而增大，严重时可使算出的安全系数比其它较严格的方法的结果小一半。而毕肖普法 考虑了土条侧面的作用力，并且假定各土条底部滑动面上的安全系数均相同，即等于整个滑 动面的平均安全系数。因此毕肖普法是比较合理的土坡稳定分析方法。在土坡稳定安全系数的计算中，由于滑

2、动圆弧的圆心和半径都是任意假定的，计算出的安 全系数不一定是最小的安全系数，所以需要多次试算，假定多个滑裂面才能找到计算土坡的 最小安全系数。这就使得求解过程虽然不复杂，但计算量很大的土坡稳定安全系数的计算需 要花费大量的时间，因此人们的视线自然而然的转向了利用计算机来缩短计算时间这个方 向。这样的环境之下，考虑到窗口界面已成为程序设计的基本要求，优选visual- basic语言 计算土坡的稳定安全系数。本程序不需要输入公式，只要输入土体容重、内摩擦角、凝聚力、 土坡高度、土条宽度、坡脚等一些参数，即可计算出土坡的最小安全系数。1计算原理如图1所示，假定滑动面为一圆心为0,半径为R的圆弧。任

3、一土条中，其上的作用力 有土条自重Wi，土条底部的总法向力Ni和总切向力T i、条块间的法向力Ei、Ei+ 1和 切向力Xi、X i + 1。共有7个未知力。为使问题求解，毕肖普假定可忽略土条间的切向 力的作用。1. 1滑动面圆心位置的确定滑动面圆心位置的确定采用费伦纽斯确定最危险滑动面圆心的方法。如图所示:D点的位置距坡脚A点的水平距离为4.5 H，竖直距离为H。0点的位置为 从坡顶和坡底引出的与坡边坡和坡顶分别成B 1，B 2的两条直线的交点。最危险滑弧的圆 心在D 0直线的延长线附近。圆心01，0 2,，对应的圆弧分别求得稳定安全系数K 1，K 2,绘出K值曲线可得到最小安全系数值K m

4、 i n,其相应的圆心0 min即为最危险滑 动面的圆心。考虑到土坡的最危险滑动面圆心位置有时可能在D0直线的之外，因此通过 0 min点作D0线的垂线，在垂线上取几个点试算滑动面的圆心0c 1、0 c 2,并计算稳 定安全系数Kc 1、Kc 2,绘得Kc值曲线，最小的安全系数(Kc m i n)的对应的圆心，才 是最危险滑动面的圆心。土坡城度坡侑aPiX 1核3了1： L 5特：23亨Ik 311:4ModII; 5i邛371. 2土坡稳定安全系数的推求确定了圆心位置之后，如图1所示，根据任一土条i上的竖向力的平衡条件可得:图I最危险滑弧的确定N i cos A i = Wi - T i s

5、in A i ()假定土坡的稳定安全系数为K ,则i 土条滑动面上的实际产生的抗剪强度为 TOC o 1-5 h z T.二 皿由+(2)将公式(2)代入公式(1 )得：Ni ( Wi c7zsin t/K )/(cos 仕卜tail 砰sin q/K j( 3J按滑动体对圆心的力矩的平衡有：少二 IFMinM(4)将公式(2) , ( 3)代入(4)整理后可得：犬=志r晶心e +Witail 氽 Wisin %(5J其中M折二 ms(J- + tan 晚sin(h/K (6)公式(5)为简化的毕肖普法计算土坡稳定安全系数公式。由于式中MH i也含有K值, 因此公式(5)须用叠代法求解，即先假

6、定一个K值，按公式(6)求得M H i的值，代入公 式(5)求出K值，若此K值与假定不符，则用此K值重新计算M H i求得新的K值, 如此反复叠代，直至假定的K值与求得的K值相近为止。2计算机程序的编制及流程程序采用的编程环境是微软公司的VisualBasic 6. 0。程序具有通用性和结构化的特点。注: 在程序中未计入渗流作用。2. 1 程序编制计算程序分为四大部分：确定DO直线方程DO的直线方程为：X - ( 7()+ H) (Lan BVtsm QJ + 1 - 4_ 5 x tan (Hi + + tan Ss/tan ( Pi -k dj x (tan By tan aj + tan f 3i + d) + Lan P1+ 4. 5 x /点坐标为Xa - (tan 免 x H/Lan 1+ H )/Lan 3i + tan的匕二 m Lan f 3| + a) x tan fVtan a+ H j/f tan f PI + + Mfl A,求公式(5)中的各个未知量。公式(5)中的c