有限元强度折减法.docxVIP

实用文档 有限元强度折减法 背景 1974 年，Smith Hobbs[1] 使用有限元方法分析了 φu=0 条件下的边坡稳定性并与 Taylar[2] 的结果进行对比，得到了很好的一致性； 1975 年， Zienkiewicz 等 [3] 考虑 c’、φ’进行有限元边坡稳定性分析，其结果与圆弧滑面解有较好吻合； 1980 年 Griffiths [4] 验证了一系列具有不同材料特性和形状的边坡稳定性并通过与 Bishop Morgenstern[5] 的结果进行了对比确定了数据的可靠性；此后也有研究证实了利用有限元方法进行边坡稳定性分析的可靠性 [6,7,8,9]；在文献 [9] 中，引入一些 案例证明了有限元强度折减法的准确性， 并证明了有限元强度折减法在分析非均 质边坡时相对于传统方法的优越性。 2001 年，郑颖人等 [10] 把有限元强度折减法引入国内，并对此进行了后续研究 [11,12,13,14]。 相较于一些传统的边坡稳定型分析方法， 有限元强度折减法有以下几个优点 [9] ： 不必假设滑面的位置和形状，当土体自身强度不足以抵抗剪应力时土体失稳会自然发生。 由于有限元强度折减法中没有条分的概念，因此也不必假设条间力，在整体失稳之前土体都处于整体稳定状态。 使用有限元方法能够查看破坏过程。 有限元强度系数折减法 1.模型参数 边坡模型主要包括六个参数，分别是：膨胀角 ψ、内摩擦角 φ’、黏聚力 c’、弹性模量 E’、泊松比 υ’、重度 γ。 膨胀角影响土体屈服后的体积变形， 若 ψ0，则土体屈服后体积减小， 若 ψ0 则体积增大， ψ=0则体积不变。 ψ=φ的情况被称之为关联流动法则，但是此时 ψ 值通常高于实验观测值， 特别是在侧限条件下会提高土的承载力预测值。 边坡稳定型问题通常是处于无侧限条件下， 此时膨胀角的选取不再重要 [9] ，因此文献 [9] 选取 ψ=0条件下的非关联流动法则，并且通过案例分析可以得出此膨胀角的选 取可以得出准确的安全系数以及滑动面。 c’和 φ’指 Mohr-Coulomb 准则中边坡土体的有效黏聚力和内摩擦角； E’和 υ’ 是土体材料的弹性参数，这两个参数对土体稳定性分析的影响较小； γ是土体的 重度。应用有限元方法进行边坡稳定性分析中最重要的三个参数是 c’、 φ’、和 γ。 2.屈服条件 (1)Mohr-Coulomb 准则 Mohr-Coulomb 准则用大小主应力表示如式 (1)所示： (1) 其中， 、 分别指土中一点的大小主应力。在主应力空间中，如果不考 虑 、 、 之间的大小关系， 屈服面是一个不等角六棱锥， 在 π平面上是一 个等边不等角六边形。 (2)D-P 准则 标准文案 实用文档 D-P 准则可以写成式 (2)形式： (2) 其中 I1 为第一应力不变量、 J2 为第二偏应力不变量， β和 kf 为试验常数。在主应力空间中其屈服面为一个圆锥，在 π平面上是一个圆形。 3)D-P 准则转换为 Mohr-Coulomb 准则 首先引入参数 b，如式 (3)所示： (3) 则， 和 分别可转化为式 (4)： (4) 其中 将其带入 (2)，得式 (5)： (5) 与式 (1)对比可知两个准则之间的转换关系如式 (6)所示： (6) 因此，当 b=0 时，即外角点外接 DP 圆的两个试验常数分别如式 (7)所示，当 b=1 时，即内角点外接 DP 圆的两个试验常数分别如式 (8)所示。 ， (7) ， (8) =23 σ b = 0 σσ2 σ3 1 σσσ 1≥ ≥ 1≥ 3 ≥ 2 σ 外角点外接D P 圆 b = 1 σσσ σσσ 3≥1 ≥ 2 2≥1 ≥ 3 内角点外接D P圆 σ σ 2 3 ≥ 3 ≥ 1 σ 3 ≥ ≥ 1 2 2 3.安全系数的定义 (1)Mohr-Coulomb 准则中的安全系数 [15] 1955 年， Bishop 首先在边坡稳定性分析中提出了抗剪强度折减的概念， 在有限元强度折减法中通过将坡体的强度参数： 黏聚力 c 和内摩擦角 φ同时除一 标准文案 实用文档 当计算不收敛时，对应的 Ft 即为所求的安全系数，此时坡体达到极限状态，发生剪切破坏。 c’ =c/Ft φ’ =arctan(tant) φ /F (2)D-P(Drucker-Prager)准则中的安全系数 取 Ft 为 D-P 准则中的强度折减系数，则 D-P 准则可以表示为式 (9)， (9) [13] (3)不同屈服条件下安全系数转换 首先引入 Mohr-Coulomb 等面积圆屈服准则，在 π平面上，其屈服面是一个圆，并且面积与 Mohr-Coulomb 准则的不等角六边形相等， Mohr-Coulomb 等面积圆屈服准则中的试验参数