用EXCEL迭代求解Morgenstern-Price法方法及应用.docVIP

用EXCEL迭代求解Morgenstern-Price法方法及应用 　　摘要：首先推导Morgenstern-Price法差分格式，在Excel表格中编写Morgenstern-Price差分求解公式；在Autocad中对计算剖面分条，量测各土条侧边顶底点坐标，把所得坐标复制到Excel指定表格中，计算得到稳定系数、水平条间力等。算例分析表明，该方法简单易行，可操作性强，计算结果满足工程要求。 　　关键词：极限平衡法，EXCEL迭代，Morgenstern-Price法，稳定系数 　　 　　1．引言 　　 目前，建立在极限平衡原理基础之上的条分法一直是边坡稳定性分析计算的主要手段[1] 。在极限平衡法中，Morgenstern-Price法全面考虑了力与力矩平衡[2、3]，是一种严格的计算稳定性的方法，且适用于任意形状的滑动面，但该方法要用Newton-Raphson迭代法进行迭代求解，编写Newton-Raphson迭代程序较为复杂。在工程实际中还没有广泛应用。 　　 L.H.Chen[4] 在Excel中编写Newton-Raphson求解法，使用Autocad软??加载Autolisp程序采集图形几何信息，得出稳定系数等计算结果，但Autolisp及Excel程序编写仍然较复杂。 　　 鉴于此，本文使用Excel电子表格编制Morgenstern-Price差分法公式，得到稳定系数、水平条间力等稳定性指标，有利于Morgenstern-Price法进一步广泛应用。 　　2．Morgenstern-Price法的差分格式的推导 　　 在图1的坐标系下，Morgenstern-Price法的力和力矩平衡公式分别如式(1)和式(2)： 　　 　　 　　 　　图1Morgenstern-Price法条块分析示意图 　　 　　 　　 　　 (1) 　　 (2) 　　其中：E为条块间的水平推力，X为条块间的垂直剪力，c 为有效粘聚力，为有效摩擦角，y为滑动面函数，yE为推力线位置函数，u为滑面上的孔隙水压力。 　　 在求解过程中，该方法假定土条垂直剪力与法向作用力为： 　　 　　 　　底滑面方程为： 　　 　　土条重为 ： 　　 　　对侧向力之间关系函数假定为： 　　 　　 　　 　　将已知条件代入（1）式推导后可得： 　　 　　 　　该式仅在[xi,xi+1]内成立，属于边值问题，边界条件为： 　　 　　 　　 　　该微分方程的解为： 　　 　　其中C为待定常数。 　　 将上式中的x看成局部坐标，以条块左边界为原点，当x=0时，，则 　　 (3) 　　 将（2）式变换得： 　　 　　 通过积分推导后可得： 　　 (4) 　　 采用差分格式，若用x表示每一条的宽度，则(3)和(4)式中有关项中的x用x代替，得 　　(5) 　　 (6) 　　式 　　 　　(5)和(6)为Morgenstern-Price法的显式求解公式。式中包含两个未知量F和，初始边界条件是E0=0，M0=0；最后一个土条的边界为En=0，Mn=0。首先给定E0=0，M0=0，逐条依次求解，最终求得En和Mn，调节F和，使En=0，Mn=0，此时的F和就是要求的收敛解，该求解步骤可在Spreadsheet下很容易实现。 　　3．在EXCEL下的求解格式 　　 在Excel表格中，首先输入基本数据，即滑体中各土条侧边顶底点坐标、分层线坐标和水位的坐标，及滑坡体重度和强度指标。各土条的几何尺寸可通过Autocad直接量测，然后在文本窗口中，复制各点坐标，粘贴到电子表格中。 　　 在Excel中编写上述公式，见图2，计算滑体各土条的宽度及水位线上下各土条的面积，求出土条重，计算式中的p、q、k、m值，再求出上式的K、L、M、N，即Ei值递推得解。同样，力矩公式也可按此法编写。 　　图2Morgenstern-price法Excel求解表达式 　　4．实例计算 　　 沿江路小商品市场滑坡位于三峡库区巴东县西壤坡新城区境内，总体地势南高北低，呈上宽下窄的梯形展布，后缘宽140m，前缘宽78m，滑坡轴向长250m，相对高差74m，主滑方向307。斜坡为顺向坡，坡度25°~45°。 　　 滑体主要为第四系松散堆积层及三叠系中统巴东组(T2b2)地层。滑体内有二层滑带（图3），浅层滑带以上以碎石土和块石土为主，浅层滑带和深层滑带之间为碎裂岩。 　　 浅层滑带为青灰色泥化带，厚度较薄，该层滑体平均厚度13m，最厚20.5m，滑带后陡前缓，后部坡度为20，中部及前缘坡度为9。 　　 深层滑带青灰色泥化带，厚度不均，在破碎岩体中，顺层滑动，平均埋深30m，滑带后陡前