考虑中主应力效应强度折减法的边坡稳定性分析.docVIP

考虑中主应力效应强度折减法的边坡稳定性分析 基于考虑中主应力效应的屈服准则编制强度折减法 程序，通过有限差分软件FLAC3D对边坡稳定性问题进行分 析，通过分析验证表明：计算结果与Bishop理论计算结果 更为接近，优于基于Mohr-Coulomb理论的计算结果，更具 有工程应用价值。 关键词：中主应力效应；强度折减法；安全系数；边坡 稳定性 U213. 1+3 A 0引言 在土坡稳定分析中，土坡在外荷载和土重作用下失稳 破坏会形成一个贯通滑裂带，整个土坡将沿滑裂带滑动［1］。 由于，滑裂带以上土体为机动土体，数值计算解不会收敛，为 解决这个问题强度折减法应运而生。 而在边坡稳定性分析中，技术人员通常采用的强度折减 理论，一般的建立在Mohr-coulomb强度准则的基础之上， 它们主要缺点是不能考虑中主应力效应的影响。相比之下， SMP屈服准则具有更强的物理推导过程，各变量物理意义明 确，考虑了中主应力的影响［2］,对于土质边坡更具有优势。 此外，传统的边坡稳定性分析方法中，为了便于分析计算的 进行，做了许多近似假设[3,4],如不考虑土体内部应力-应 变关系、不考虑支挡结构的作用等。因此，传统分析方法不 能得到滑体内的应力、变形分布状况。而有限差分法克服了 传统分析法的不足，不仅满足力的平衡条件，而且还考虑了 土体应力、变形关系，能够得到边坡在荷载作用下的应力、 变形分布，模拟出边坡的实际的滑移面。基于这些优点，近 年来有限差分法广泛应用于边坡稳定性分析。基于这些原 因，本文基于考虑中主应力效应的屈服准则建立强度折减法 程序，通过有限差分软件FLAC3D对边坡稳定性问题进行分 析，以期能更好地解决工程计算问题。 1平面条件下的SMP准则 0本名古屋工业大学的Matsuoka和Nakai于1974年提 出的SMP准则[5],考虑了中主应力的影响，建立的三维条 件下的土的强度公式，该破坏准则认为当剪应力与正应力的 比值达到某一数值时材料破坏，表达式为 (1) 式中为内摩擦角，K为材料常数，，，分别是应力的一次, 二次，三次不变量，用下式表示： 图1 Mohr-Coulomb和SMP准则在平面上的破坏面 令A二则由式(1)得 (3) 其中，为角隅函数的平均值，即：， 则 以上SMP准则的公式只适用于摩擦性材料而不适用于C- 材料，Mastuoka[6]于1990年做了如下修改，引入了黏结应 力o 0,表达式为 修改的SMP准则表达式为 (4) 式中： (5) 将公式(4) (5)联立并根据相关联流动法则可得黏性 材料的SMP屈服准则的表达式为 (6) 通过常规三轴实验求出土体的内摩擦角?与黏聚力C, 求出土体的K和A值，进而可由式(6)求出土体考虑中主应 力效应的SMP屈服准则。 2强度折减中的安全系数 抗剪强度折减系数的概念是由Zienkiewicz[7]等1975 年在土工弹塑性数值分析一文中首次提出的，由此所确定的 强度储备安全系数与Bishop[8]在极限平衡法中所给的稳定 安全系数在概念上是一致的。工程中常用的安全系数是在材 料强度保持一定情况下，把荷载提高一定的倍数，相当于荷 载增大系数或加载系数；而强度折减系数是把材料强度降低 一定倍数，两者在实质上是一致的。 强度折减法通过调整岩土体的强度指标c和e,表达式 为 (7) 通过逐步调整系数k,得到不同的CF将?F,,反复分 析边坡，直至坡体达到临界状态，此时得到的折减系数即为 安全系数Fs。运用强度折减差分法分析边坡的稳定性，采用 解的不收敛作为破坏标准。当在指定的收敛准则下算法不能 收敛，表示应力分布不能满足土体的破坏准则和总体平衡要 求，意味着出现破坏。 3工程应用 某边坡工程为均质天然边坡，边坡高度H = 10m,坡度为 tana=l： 1,其相关材料参数见表1,坡顶均布矩形超载q= 120 kPa,超载的宽度为2m,超载边缘离坡顶边缘的水平距 离 b= 1. 5 mo 表1材料参数 P /kg. m-3K/MPa G/MPa o t/MPaC/kPa(I)/ordm； 2000 10030 100000. 3223 图2分析模型示意图 图3考虑中主应力下的安全系数、剪切应变增量云图 图4基于Mohr-coulomb准则的安全系数、剪切应变增 量云图 表2不同方法对应的安全系数 屈服准则与方法M-CSMPBishop 安全系数 Fsl. 121.041.031 由图3、4可知,基于Mohr-Coulomb准则与SMP准则的 结果云图中，当达到极限状态下二者均出现了贯通的塑性 区，而基于SMP准则的计算云图的塑性区发展更为显著，主 要是由于其考虑了中主应力影响的缘故，从这个角度可知 Mohr-Cou 1 omb准则在计算中会偏于保守。 另外，从表2的计