基于有限单元法和极限平衡理论的边坡稳定性分析研究.pdfVIP

基于有限单元法与极限平衡理论的边坡稳定性分析研究 郑霞霞 (泾阳县灌溉管理工作站 陕西 泾阳 713700) 摘 要 本文将边坡稳定性分析的极限平衡法与应力变形计算的有限元单元法相结合，对边坡的应力、 变形分析采用有限元单元法，求出边坡应力场与位移场，然后利用有限元对应力的分析计算结 果，通过插值积分处理，求出土条底部的正向与法向力，并根据极限平衡法原理得出边坡的稳定 安全系数。该方法既能反映边坡的稳定和变形之间的关系，又能用工程界所熟悉的安全系数来 评价边坡的稳定性，具有较广泛的应用前景。 关键词边坡稳定；有限元；极限平衡理论；应力 中图分类号：TU457 文献标识码：A 1引言 表示如下： 系数方程： dv+p J，…’【c】田…a)=b』，N1fA ∑【c’BR+(N—uB)Rtan‘P1 目前对边坡稳定性分析的理论方法主要 NTdA+{R) (1) Fm-——(4) 有两类：①极限平衡法，②有限单元法。极限 式中：【B】表示应变一位移矩阵；【C】表示本 ∑Wx一∑Nf±∑Dd 平衡法是建立在刚体极限平衡理论之上，假 构关系矩阵；(a)表示节点增量的x．y位移列 定滑坡体为刚体，不考虑岩土体本身的变形 一 向量；N表示插值函数的行向量；A表示沿 22 对边坡稳定性的影响，假定条件较多，使得分 单元边界的面积；v表示单元体积；b表示单位 聒——(5) 析结果准确性欠佳。但这类方法能给出的边 体力强度；p表示增加的表面压强；{Fn)表示 ∑Nsinot一∑Dcosto 坡稳定安全系数物理意义明确，同时能够给 集中节点增加荷载。 出可能的滑动面，因此这类方法成为工程界 式中：c7为有效粘聚力；夺’有效摩擦角； 对于二维平面应变问题，认为所有的单 对边坡稳定性分析中应用最广泛的方法。以 U孔隙水压力；N土条底部法向力；W土条重 元都是单位厚度，对常单元厚度t，上述方程 有限元法为代表的数值计算方法是一种新生 可写为： 量；D集中荷载；B、R、x、f．d、∞为几何参数； 的方法，这类方法能够考虑边坡岩土体本身 仪为土条底面倾角。 的变形对边坡稳定性的影响，能够给出边坡 t f。【B】7[Cl[BldA{a)=btf。N’dA+PtIL 本文对土条底部的法向力和切向力采用 岩土体的应力场与位移场，以及计算边坡的 NML (2)有限元计算，通过在土条底面进行积分求得。 破坏发生和发展过程等。但是这种方法也不 有限元方程的简化形式为： 1=‘F 是完美的，它很难给出一个具有明确意义的 【K】‘a 1=‘Fb1+‘F3’+‘Fn’(3) 2工程概况 稳定性安全系数以及可能的滑动面，限制了 r【c】【B】 式中：【k】为刚度矩阵，[Kl=t}．【B 其工程应用。虽然许多学者致力于解决此问