有限元强度折减法.doc

有限元强度折减法1背景1974年，Smith & Hobbs［1 ］使用有限元方法分析了 血=0条件下的边坡稳定性并与Taylar⑵的结果进行对比，得到了很好的一致性；1975年，Zienkiewicz等［3］ 考虑c' ©进行有限元边坡稳定性分析，其结果与圆弧滑面解有较好吻合； 1980年Griffiths⑷验证了一系列具有不同材料特性和形状的边坡稳定性并通过与 Bishop & Morge nsterN5］的结果进行了对比确定了数据的可靠性；此后也有研究证 实了利用有限元方法进行边坡稳定性分析的可靠性 Q"9】；在文献［9］中,引入一些案例证明了有限元强度折减法的准确性，并证明了有限元强度折减法在分析非均 质边坡时相对于传统方法的优越性2001年，郑颖人等［10］把有限元强度折减法 引入国内，并对此进行了后续研究［11,12,13,14］相较于一些传统的边坡稳定型分析方法，有限元强度折减法有以下几个优点［9］:(1) 不必假设滑面的位置和形状，当土体自身强度不足以抵抗剪应力时土体 失稳会自然发生2) 由于有限元强度折减法中没有条分的概念，因此也不必假设条间力，在 整体失稳之前土体都处于整体稳定状态。

3) 使用有限元方法能够查看破坏过程2有限元强度系数折减法1. 模型参数边坡模型主要包括六个参数，分别是：膨胀角 职内摩擦角黏聚力c'弹性模量E'、泊松比u'重度丫膨胀角影响土体屈服后的体积变形，若书＜0则土体屈服后体积减小，若书＞0 则体积增大，书=0则体积不变书=©的情况被称之为关联流动法则，但是此时 书值通常高于实验观测值，特别是在侧限条件下会提高土的承载力预测值 边坡稳定型问题通常是处于无侧限条件下， 此时膨胀角的选取不再重要［9］，因此文献［9］ 选取书=0条件下的非关联流动法则，并且通过案例分析可以得出此膨胀角的选 取可以得出准确的安全系数以及滑动面c和©指Mohr-Coulomb准则中边坡土体的有效黏聚力和内摩擦角； E'和u、是土体材料的弹性参数，这两个参数对土体稳定性分析的影响较小； 丫是土体的重度应用有限元方法进行边坡稳定性分析中最重要的三个参数是 c' ©、和Y2. 屈服条件(1) Mohr-Coulomb 准则Mohr-Coulomb准则用大小主应力表示如式(1)所示:??-??32??+??'2?????????? ????????其中，？?、？?分别指土中一点的大小主应力。

在主应力空间中，如果不考虑??、??、??之间的大小关系，屈服面是一个不等角六棱锥，在n平面上是一个 等边不等角六边形2) D-P准则??-??3??-??3》(?? - ??)「2 + (??-??v??= (??-V3??)其中？?= V1- ?? ??将其带入(2),得式⑸：??- ??=1.5 ????? (?? + ?3?) + ?? ??2-????+巧 _-????+^3与式(1)对比可知两个准则之间的转换关系如式 ⑹所示：1.5sin ?? = ?? ??-? ■- ????二？ 2 ■■…心 -???2c?qos (??= ?? ??-2 ■…V3因此，当b=0时，即外角点外接DP圆的两个试验常数分别如式 ⑺所示，当 b=1时，即内角点外接DP圆的两个试验常数分别如式(8)所示2????????卩=v3(3-????????) 2????????卩=v3(3+???????，???=???=6??????????v3(3-????????)6??????????v3(3+????????)(8)D-P准则可以写成式(2)形式：-B1 + Vv2 = ??? (2)其中Ii为第一应力不变量、J2为第二偏应力不变量，B和kf为试验常数。

在 主应力空间中其屈服面为一个圆锥，在 n平面上是一个圆形3)D-P准则转换为Mohr-Coulomb准则 首先引入参数b,如式(3)所示：b =贝?1和二?分别可转化为式(4):_ 3(?? + ????=0》02> 001外角点外接D芦 T> 0203=0 b = 1=3b=0(,)2> 3 1 0》2 > 10》0 T 03 内角点外接D P圆02*3. 安全系数的定义(1)Mohr-Coulomb准则中的安全系数1955年, Bishop[15]首先在边坡稳定性分析中提出了抗剪强度折减的概念， 在 有限元强度折减法中通过将坡体的强度参数：黏聚力 c和内摩擦角©同时除一 个折减系数Ft,得到一组新的c'和©值，作为一个新的强度参数输入进行试算，当计算不收敛时，对应的Ft即为所求的安全系数，此时坡体达到极限状态，发生 剪切破坏c' =c/F ©' =arctan(tant) © /F(2) D-P(Drucker-Prager)准则中的安全系数取Ft为D-P准则中的强度折减系数，则 D-P准则可以表示为式(9),?? ???-?????+ "?=鬲(3) 不同屈服条件下安全系数转换［13］首先引入Mohr-Coulomb等面积圆屈服准则，在n平面上，其屈服面是• 圆，并且面积与 Mohr-Coulomb准则的不等角六边形相等，Mohr-Coulomb 积圆屈服准则中的试验参数如式(10)所示：sin © (9)等面・3 =V3(乜？？？？？?？？？？？?????？?？？)v3ccos ©(10)???=-v3????????????????????12 4 2 2 ?????? ?? 2.一？？？？？？？？？？^??????^ +1)(~1)]式中???= ???????????? 9 ?????? 3 —玄宁-1)简称外接圆屈服准则为DP1准则，其试验常数分别为 3,kf1 ； Mohr-Coulomb 等面积圆屈服准则为DP2准则，其试验常数分别为 血，kf2。

把DP1准则表示为 f1 = V??= ??!?+ ???1， DP2 准则可表示为 f2 = V??= ????+ ???2 o 令 f1 ?????+?????2n =3 2=kf1\kf2=f( ©,) f1 = ????+ ???1= ?????+ ????，所以 W =选??+???2 = ??= ??(??)由此可知，n是©的函数，当©取不同值时可以得到不同的 n值如表1所 列：表1不同内摩擦角时的 n值

郑颖人［14］提出，有限元的计算迭代过程就是寻找外力和内力达到平衡状态 的过程，整个迭代过程直到一个合适的收敛标准得到满足才停止 可见，如果边坡失稳破坏，滑面上将产生没有限制的塑性变形，有限元程序无法从有限元方程 组中找到一个既能满足静力平衡又能满足应力-应变关系和强度准则的解，此时 不管是从力的收敛标准，还是从位移的收敛标准来判断有限元计算都不收敛3案例分析例一，不含地基的均质边坡［9］该边坡如图1所示，有限元程序采用Mohr-Coulomb失效准则，建立平面应变条件下八节点四边形单元减缩积分计算模型， 其强度参数为©' =2氏’/ 丫 h=d05边坡坡度为26.57 ° (2:1)坡底水平，其边界条件为坡底约束竖直方向位移与水平 方向位移，左侧约束水平方向位移，其余面为自由面施加重力荷载后使安全系 数从0.8到1.4逐步变化直至计算不收敛1—— 仲-一一一+ 冷 巴詈常图1不含地基的均质边坡每一个安全系数对应的迭代次数如表 2所列，当真正的安全系数接近时需要 更多的迭代次数表2例一计算结果FOSE r^HUK/y h'Iterationso-8o n037921 000 3KIJO1-20042220b300453411-350 5447921-401 476lOOt)当安全系数为1.4时，无量纲位移E' nSax/ yH突变，并且此时计算无法收敛， 在此情况下有限元计算结果与 Bishop & MorgensterH5］给出的结果吻合良好，女口图2所示。

Bishop & Morgertstern (I960)FOS = 1 380oa 1 V2 ?4 ~V6FOS图2安全系数与无量纲位移边坡失稳时(F0S=1.4)节点位移矢量和网格变形如图 3(a)和图3(b)所示，由此可得到边坡的潜在滑动面9)㈣图3安全系数为1.4计算不收敛时边坡变形 (a)节点位移矢量(b)网格变形例二，有软弱层的不排水黏性土边坡在本案例中，使用Tresca准则((U=0)进行总应力分析边坡几何形状如图 4所示，地基厚度与边坡高度相同，该边坡有一个软弱层，在有限元计算中，令其 抗剪强度(Cu2)在一定范围内变化但其周围土体抗剪强度保持 Cui/ 丫 H=0.2不变利用有限元方法计算该边坡的安全系数结果如图 5所示，对于均质边坡情况，Cu2/Cu1=1，有限元计算结果与Taylor⑵的结论很接近，随着软弱层的强度逐 渐减小，在Cu2/Cu1- 0.6寸，结果发生了明显的变化分别假定圆弧滑面和穿过软 弱面的三段线滑面并利用Janbu法计算安全系数，可见在 Cu2/Cu1 - 0.6处也发生 了滑动机制的转换，当Cu2/Cu1>0.6时，潜在滑面形状为圆弧，当Cu2/Cu1<0.6时， 潜在滑面为结构软弱面。

图 6更加清晰的展示了这一现象，图 6(a)为均质边坡(Cu2/Cu1=1)时的潜在滑面，可见此时的滑面形状为圆弧滑面，与 Taylor[2]的预测相同；图6(c)为软弱层强度只有其周围土体 20%( Cu2/Cu1=0.2)时的潜在滑面，此 时潜在滑面沿软弱层发展；图6(b)为软弱层强度只有其周围土体60%( Cu2/Cu1=0.6) 时的潜在滑面，此时圆弧滑面和沿软弱层的三段线式滑面都有可能发展， 至少存在两种明显的滑动机制1 p3 -1614121 -0 80硏04 -02Taylor (1937)FOS = 1 -47唧YH = 02500图5不同软弱层强度时的安全系数图例三，不同坡度边坡安全系数计算[13]，验证Mohr-Coulomb等面积圆屈服准均质边坡，坡高H=20m,土容重丫 =25kN/m,黏聚力c=42kPa，内摩擦角 © =17求坡角B分别为30° ,35 ° ,40 ° ,45时边坡的安全系数。