强度折减时边坡稳定性分析..doc

强度折减时边坡稳定性分析 E. M. DAWSON,W. H. ROTH and A. DRESCHER 关键字：堤坝;滑坡;极限状态设计分析数值模拟和分析;可塑性斜坡。 在这里Hc指斜坡的临界高度，r是土体重度，c是粘聚力。 数值Ns很直观的把斜坡倾角，摩擦角和孔隙水应力系数联系起来。Bishop (1954)提出来孔隙水应力系数，并把它定义为孔隙水压力与上覆总压力之比。距离地表z处的孔隙水压力u可以用下式表示： 稳定性分析的例子 假设有一个均质堤坝，高度为H=10m，斜坡倾角为，内摩擦角，土体单元 ，粘聚力，。具有这些土体特征的斜坡用Chen的极限分析法算出的安全系数为正好为1.0。 把堤坝用FLAC在平面应力和小应变的情况下进行模拟。这部分土体被认为是线弹性或是完全塑性材料，符合Mohr-Coulomb屈服条件和相应的流动规则。在表格1表示的数字网格，宽和高均为20个单位。水平位移固定在沿着横向和纵向的节点上，而水平位移和垂直位移沿着底部边界进行。为了用强度折减法来进行斜坡稳定性分析，进行模拟时要对一系列的实验安全系数按照以下有关c和进行修正： （3） （4） 图表1均质路堤的数值网格 图2表示由剪应力小幅度降低引起的标准化的不平衡力，其中从0.8开始。图中表明随着的增大，土壤中的剪应力降低。当实验安全系数达到1.02时，模拟集中体现不平衡力为数量级的一种平衡。但是，当到达10.3时，不再表现一种平衡了，不平衡力达到了数量级。这似乎有点荒谬，当增大到一定值时滑坡就会发生。但是，F的定义为一个系数，剪应力按照这个系数来分解从而导致斜坡的破坏。当应力进一步减小时，不平衡力继续近似线性增加。因此，强度折减法给定的安全系数的范围是1.02到1.03之间。为1.03时，斜坡破坏时沿着按有限分析法得到的临界对数螺线表面的速度场如表3所示。表1中不平衡力的剧烈转折表明斜坡破坏时的实验安全系数很明确。这是应用了线弹性和完全塑性模型的结果，这种模型存在直接从线性到塑性的突然转变。一些强度折减法的前期应用，例如Matsui San(1992)曾经用过双曲线型的基本模型，这种表现了一种从线性到塑性的平缓的转变。应用这种模型时来判断极限状态时比较困难。 标准化不平衡力 试验安全系数 图表2试验安全系数小规模增加时的不平衡力 图表3 临界对数螺线表面的速度场 滑坡发生时的值用交叉和平分法可以比较容易得到。首先，建立一个较高支点和一个较低支点，最初的较低支点可以是模拟表现的任意一个，而较高支点是模拟不能表现的任意。接着，测试较高支点和较低支点之间的中间点。加入试验可以表明，则用这个点代替较低支点。假如试验不能表明，则代替掉较高点。重复这个过程，直到较高点和较低点的差距在特定的可容许范围内。 基准研究 为了判断强度折减法的精度，模拟要参数很广的范围进行。模拟的堤坝的倾角的范围为到，给定的土体的值要在到的范围，的值要取0.0,0.25，和0.5。每一个以上指标的组合，Chen (1975) or Michalowski (1995a)给出的稳定性数字，来计算斜坡的高度，内聚力和单元重度，堤坝将会得到一个接近1.0的稳定系数。 实际上，Mohr-Coulomb屈服条件常常会与一个张力切断相结合。为了与极限分析法作对比，这里就不再使用张力截断。相反的，Mohr-Coulomb屈服准则只在拉伸的正应力区域有效。换句话说，破坏包络线与正应力的曲线相交。假如应用零拉力截断法，这样土体的单轴拉应力为零，拉裂破坏发生在堤坝的顶端，计算出的安全系数会小几个百分点。 模拟既要以细网格（60*60单元）又要以粗网格（20*20）进行。安全系数要以交叉和平分法进行直到较高支点与较低指点差距小于。报告中填写的安全系数是较高点和较低点的平均值。以为0.0,0.25和0.5算出的安全系数分别在表一、表二、表三中。强度折减下地安全系数通常会与极限分析法得到的结果相差几个百分点。带陡边坡和较大的内摩擦角时这种差别是最大的。就像我们预想的一样，细网格得到的结果比粗网格要好。 垂直边坡，分别为0.25和0.5时，用极限分析法得到的结果很相近（见于表二、三的最后一行）。下面是由于高孔隙水应力和零水平总应力而引起的坡脚处破坏的模拟。例如，=0.5时，坡脚处斜坡表面的孔隙水应力就会超过土体的抗拉强度。 通过图表对比，可以很方便的表现出稳定值的数字性结果。相应计算出的安全系数的安全值可以用方程（1)和（3）得到。在=0.0时，稳定值对应的内摩擦角