计算机在边坡稳定分析中的应用.pdfVIP

第十四届全国工程设计计算机应用学术会议论文集 杭卅．2008 计算机在边坡稳定分析中的应用 李方范幸义 (重庆大学土木工程学院) 【摘要】边坡失稳产生的滑坡现象同地震和火山并称为全球三大地质灾害，不仅会造成重大的经济损失， 而且造成人员伤亡．由于边坡内部结构的复杂性和组成边坡岩石物质不同，造成边坡破坏具有不同模式， 对于不同破坏模式就存在不同形态的滑动面．因此也应采用不同的分析方法和计算公式来分析其稳定状 态，并采用不同的支护结构防止发生失稳。本文主要介绍边坡稳定分析的分析计算方法。以及各分析方法 在计算机方面的应用与发展。 I关键词】边坡失稳；分折计算方法；计算机应用 1边坡稳定分析的方法 边坡稳定分析计算方法是边坡工程设计以及治理的主要内容，也是边坡应用的理论基 础。现有的分析计算方法，归纳起来主要有极限平衡法、极限分析法(上、下限分析法)、 变形分析(DDA)以及数值流形法(NMM)等方法。 1．1极限平衡法 极限平衡法是基于边坡处于平面应变状态的假设，对边坡进行条分，假设各条块为剐 体，求解条块的静力平衡方程，最终求得安全系数的方法。其计算过程一般先假定边坡破 坏是沿某一确定的滑裂面滑动，再根据滑裂土体的静力平衡条件和莫尔一库仑破坏准则计 算沿该滑裂面滑动的可能性，即安全系数的大小或破坏概率的高低，然后系统地选取多个 可能的滑动面，用同样的办法计算稳定安全系数或是破坏概率。最后安全系数最小或者是 破坏概率最高的滑动面就是最可能的滑动面。 1．2极限分析法 极限分析方法是将土体看作服从流动法则的理想塑性材料，分析考虑塑性的应力．应变 关系，运用塑性力学中的上下限定理求解边坡稳定问题的方法。分为上限法和下限法两种。 上限法又称能量法，计算时先假设一个滑裂面，并将土体分成若干块，然后构筑一个协调 位移场，根据虚功原理求解滑体处于极限状态时的极限荷载活或稳定安全系数。极限分析 下限法在计算过程中需要构造一个合适的静力许可的应力分布，一般可用应力柱法或应力 不连续等方法来求得问题的下限解。 1．3有限元法(FEM) 目前根据有限单元法结果来评价边坡稳定性分析方法大致可以分为两类： 第十四届全国工程设计计算机应用学术会议论文集 杭州·2008 第一类是单纯的弹塑性分析方法，它是基于弹塑性的计算原理，首先通过有限元计算 出边坡内的应力场，然后依据数学规划方法和极限平衡原理确定临界滑动面及相应安全系 数。另一类方法是强度折减弹塑性有限元分析方法，它是将强度折减技术与弹塑性有限元 方法结合，通过不断降低边坡岩土体强度使边坡岩土体达到临界状态，得到边坡的安全系 数。 1．4离散元(DEM) 又称个别单元法，它是以软弱结构面切割而成的离散块体为基本单元，对块体逐个进 行受力分析，在块体间的接触约束下运用牛顿第二运动定律描述各块体的运动过程。分为 动态松弛法和静态松弛法两种。 1．5边界元(BEM) 边界元法是在有限元的基础上发展起来的一种分析连续问题的通用数值法，通过离散 问题的边界，先求得边界值的数值解，然后通过解析公式求得计算区域内任一点的解。 1．6非连续变形分析(DDA) DDA是与有限元的连续变形分析方法相平行的一种非连续变形分析方法。该方法以位 移为未知变量，并按有限元格式求解平衡方程。有完各的块体运动学理论及其数值实现、 全一阶位移近似、严格的平衡假定、正确的能量耗散过程和较高的计算效率。DDA主要用 于非连续块体系统，并严格遵守经典的力学规则。另外该方法不仅允许块体自身有位移和 变形，而且允许块体问有滑动、转动、张开等运动形式。该方法只要给定了各点和各方向 的所测位移值就可以计算出每个块体的所有位移、变形以及应变的最d,--乘最佳拟合结果， 以及块体截面的张开、闭合和滑动。 1．7小结 边坡稳定分析方法的研究是一个不断发展进步和完善的过程。传统方法理论不断完善， 新方法不断出现和更新，如数值流行方法，它统一解决解析法、有限元法和DDA方法的 计算问题，是具有广泛发展前景的最新的数值方法。但是由于边坡稳定分析是复杂的，确 定性计算方法有时无法概括其复杂性，一个经过抽象的数学模型所求得的确定性解