1、 1 论文题目 : P hase 2在边坡分析中的应用 团队成员： 祝凯 徐振峰 陈凤 张杨 蒋宏伟 专业班级： 地 质 工 程 实 验 班 指导教师： 胡 斌 学 院 ： 工 程 学 院 提交日期： 2012年 10月 8日 2 Phase2在边坡分析中的应用 祝凯 ， 徐振峰 ， 陈凤 , 张 杨 , 蒋宏伟 摘要 ： 本文主要是介绍有限元法的原理在边坡稳定性分析上面的应用，通过建立简单的模型，利用数值模拟软件 phsae2 对模 型进行稳定性分析，模拟出边坡模型的各种变形和滑坡的情况。这种方法在分析边坡稳定性上比传统方法更为简单、精确，其分析结果对处理实际的边坡工程问题具有一定的指导意义
2、。 关键词 ： 有限元；稳定性；边坡；数值模拟； phase2 The application of phase2 in the slope analysis ZHU Kai , XU Zhenfeng , CHEN Feng , ZHANG Yang , JIANG Hongwei Abstract: This paper is mainly to introduce the use of the finite element method in the slope stability analysis . Through the establishment of a simple mode
3、l, we use phsae2 which is a numerical analysis software to analyze the stability of the model, so we can simulate the possible deformation and landslide of high slope . This method is more simple and precise than the traditional method in the analysis of slope stability , and the results of the anal
4、ysis can be very useful when actually dealing with the problem of slope stability . Key words: Finite Element; Stability; Slope; Numerical Simulation; phase2 0 引言 山区铁路与公路的建设、深基坑的开挖都会在一定程度上破坏和扰动原来较为稳定的岩土体而形成人工边坡 ,由此引起的边坡稳定性问题也不断增多。因此正确合理地选用边坡稳定性研究方法对边坡工程实际有着至关重要的意义。 有限 元 法 是边坡稳定性分析常用的 一 种方法 ，有限元法 把岩体
5、当作变形体，按照岩体的变形特性，计算岩体内的应 力分布，然后再根据 摩尔 -库伦强度理论，验算边坡的整体抗 滑稳定性 。 有限元的数值分析方法主要利用数值计算软件对边坡开挖、支护过程进行模拟，计算精度较高，应用较为广泛。本文主要是介绍利用 phase2软件对边坡进行稳定性分析。 Phase2是一款适用于地面和地表开挖、回填、设计和计算的二维弹塑性有限元分析软件 。它可以 广泛的用于工程项目中，包括脆弱岩体的复杂隧道问题、地下发电站 、露天采矿挖掘以及岩体或土壤上的斜坡。采用强度折减法， Phase2可以对边坡稳定进行有限元分析，它对于开挖产生的应力和变形均可以进行详细的计算及结果输出。 1.
6、计算模型设计 1.1 模型概要设计 本文采用的计算模型是根据工程常见的实例进行简化而后得到的，其模型用 AutoCAD 软件绘制如图 1 所示。 3 图 1.模拟模型 整个模型下部长 50m,高为 35m，边坡坡角为 40，由于五级岩体完整性较差，其与二级岩体接触的结构面呈折线型。边坡采用模拟材料是根据工程岩体分级标准 （ Standard for engineering classification of rock masses） 选取，中部选用的是二级岩体，左上角选用的是五级岩体。 1.2 材料确定 边坡采用的模拟材料是根据国家标准确定，参 考资料为工程岩体分级标准，岩体基本质量参数是根据
7、岩体基本质量的特征进行的定性分类，其标准如表 1 所示。 表 1.岩体基本质量分级 基本质 量级别 岩体基本质量定性特征 岩体基本质量指标（ BQ） I 坚硬岩，岩体完整 550 II 坚硬岩，岩体较完整； 较坚硬岩，岩体完整 550451 III 坚硬岩，岩体较破碎； 较坚硬岩或软硬互层，岩体较完整； 较软岩，岩体完整 450351 IV 坚硬岩，岩体破碎； 较坚硬岩，岩体较破碎 破碎； 较软岩或软硬互层，且软岩为主，岩体较完整 较破碎； 软岩，岩体完整 较完整 350251 V 较软岩，岩体破碎； 软岩，岩体较破碎 破碎； 全部极软岩以及全部极破碎岩 26.5 60 2.1 33 0.35
8、 本文设计的模型采用的二级岩体和五级岩体的物理力学参数选取的是表 2 参数中的高值，具体参数如表 3 所示。 表 3.二级岩体和五级岩体物理力学参数 岩体基本 质量等级 重力密度 （ KN/m） 抗剪强度 变形模量 E（ GPa） 泊松比 内摩擦角 （ ） 粘聚力 C（ MPa） 二级岩体 26.5 60 2.1 33 0.2 五级岩体 22.5 27 0.2 1.3 0.35 2.计算及分析 2.1 计算方法简介 本文采用的边坡稳定性分析的方法是基于 phase2 软件的有限元法，利用 phase2 软件对模拟的工程模型建立二维有限元模型，然后分别添加二级和五级岩 体，定义不同的材料参数，然
9、后通过软件计算得出应力、应变和 是塑性区 等结果，再对得到的图形和数据进行分析，然后得5 出该边坡是否稳定的结论。 Phase2 分析边坡稳定性有如下特点： 1.三角形或四边型等有限元网格的自动生成； 2.多样化的材料种类： Hoek-Brown、 Mohr-Coulomb、 Drucker-Prager等多种材料模型；各向同性、横观各向同性等弹性材料模型；本文采用的是应用比较广泛的 Mohr-Coulomb 材料模型。 3.支护形式锚杆（部分锚固，全部锚固，分离装置，用户自定义等） 4.地下水有效孔隙水压力分析； 该软件还具有很多其他方面的 优 点，不一一列举 。 图 2. phase2 软
10、件的操作界面 图 2即为 phase2软件的基本操作界面，根据工程实际测得参数对该软件进行操作，即可得到 分析 边坡稳定性 的数据和图像等。 2.2 计算结果简要分析 2.2.1 二维有限元模型建立 图 4 为利用 phase2 软件对本文的模型建立的二维有限元模型，灰色部分为材料性质为五级岩体，黑色部分材料为二级岩体。模型采用有限元单元划分方式为三角形划分（其具体采用的划分方式如图 3 所示），共划分 750 个三角形单元格， 410 个节点。 约束条件：模型上部为自然水平面，不考虑地面上的外加荷载，模型上部无约束条件，左右边界约束水平位移，底部边界固定。 图 3.有限单元划分方式 6 图
11、4.二维有限元模型 2.2.2 计算结果 2.2.2.1 基本情况分析 本文中的结果分析均采用应用较为广泛的摩尔 库伦准则（ Mohr Coulomb） ,计算得到应力（ 1 和 3 ）、位移（水平位移、竖向位移及变形后的整体图象）、塑 性 区等相应的图像和数据。 图 5. 主应力 1 分布图 图 5 为模型中主应力 1 分布图，由图中应力分布数据可以看出， 1 在坡面上是随着坡面的下降而不断增大的，在 二级和五级岩体中随着深度的增加也不断增大， 1 沿着五级岩体和二级岩体接触的结构面随着其下降先增大后减小， 1 在坡面拐角处应力突然增大，出现应力集中现象，这也表明该模拟结果比较准确。 7 图
12、 6. 主应力 3 分布图 图 6 为模型中主应力 3 分布图，由图中应力分布数据可以看出， 3 在坡面上是保持不变的，符合实际情况，在二级和五级岩体中随着深度的增加也不断增大， 3 沿着五级岩体和二级岩体接触的结构面随着其下降先增大后减小。 图 7. 水平位移分布图 图 7 表示的为水平位移分布图（ Horizontal Displacement）， 由图中水平位移分布数据可以看出， 水平位移 在坡面上是随着坡面的下降先增大后减小的，坡面大概中部位置水平位移量最大，五级岩体和二级岩体接触的结构面 水平位移 基本为零，五级岩体内部随着深度增加 水平位移量也是先 增大后减小的，二级岩体部分水平位
13、移量基本为零。因为二级 岩体较坚硬，完整性好，水平位移基本为零，而五级岩体由于较软，塑 性 较强而发生较大的水平位移，符合工程实际情况。 8 图 8. 竖向位移分布图 图 8 表示的为竖向位移分布图（ Vertical Displacement）， 由图中 竖向位移 分布数据可以看出， 竖向位移 在坡面上是随着坡面减小的，坡面顶部位置竖向位移量最大，五级岩体和二级岩体接触的结构面 水平位移 基本为零，五级岩体内部随着深度增加 竖向位移量也是减小的，二级岩体部分竖向位移量基本为零。五级岩体沿着结构面向下滑动 ，所以顶部竖向位移量最大，也表明该计算结果符合工程实际情况。 图 9.总位移分布图 图
14、9 表示的是模型的整个变形量及变性后的形状，由图中可以看出，五级岩体的边坡发生了很明显的变形，其左上角边界部分位移总量最大，整个五级岩体发生了较大滑动，与水平位移和竖向位移图分析的结果吻合，从图形上也表明了该边坡稳定性较差。 9 图 10. 塑 性 区分布图 图 10 表示的是整个模型的塑 性 区分布图，由图中可以看出，二级岩体 基本上为塑 性 区，其变形量很大，五级岩体基本不发生变形，该图也表明该边坡稳定性较差 。 2.3.2.2 添加 等水位线情况分析 图 11.添加等水位线图 10 图 12.添加等水位线前后位移总量对比 Phase2软件还能够模拟岩体在不同等水位线的情况下边坡的稳定性情
15、况，上图是模拟添加等水位线后边坡的稳定性改变情况， 图 11 中蓝色线条为等水位线。 两个小图表示添加等水位线前后边坡的整体位移情况，左图为下雨前模型的整体位移情况，右图为下雨后模型的整体位移情况。我们知道，下雨过后，岩体之间由于水的润滑作用，相互之间的摩擦力减小，岩体之间较下雨前更容易产生滑动，导致岩体产生滑坡的概率增大，我们比较上图也很容易发现，由于等 水位线的增加右图的边坡整体位移较左图增大，比较符合实际 。 3. 分析总结 （ 1） 本论文模型选用模拟材料为二级岩体和五级岩体， 二级岩体材料部分岩体较坚硬，整体性较好，粘聚力 、 内摩擦角 和变形模量 等参数较大， 理论上该 模型二级岩
16、体稳定性 应 较好 ；五级岩体材料部分岩体较软，岩体较破碎，粘聚力和内摩擦角等参数较小， 该 模型中五级岩体稳定性 应较 差 。计算结果表明二级岩体位移变化量基本为零，五级岩体位移变化量很大，稳定性较差， 该计算结果也比较符合 理论 实际情况。 （ 2）由该软件模拟计算得到的图型我们可以直观地看到应力、位移和塑性 区等分布图。在主 应力（ 1 和 3 ） 的应力分布图中，由于岩体自重的影响，主应力整体上是随着深度的增加不断增大，图 5 中还可以看到应力集中现象，与我们在弹塑性力学所学到的知识相对应；在 位移（水平位移、竖向位移及变形后的整体图象） 分布图中，我们可以看到模型不同位置的水平位移、
17、竖直位移和总位移情况，模型中二级岩体整体上位移变化量为零，表明二级岩体比较稳定，我们还可以观察到五级岩体边坡产生明显滑动，表明该边坡稳定性较 差；在 塑 性 区 分布图中可以看到五级岩体基本上为塑性区，这也表明该边坡稳定性差 。 （ 3） 该软件还能 模拟 添加等水位线 对边坡稳定性的影响 。下雨后 岩体之间由于水的润滑作用，相互之间的摩擦力减小，岩体之间较下雨前更容易产生滑动，导致岩体产生滑坡的概率增大， 添加等水位线后 可以观察到 边坡稳定性变差，模拟结果也符合工程实际情况。这说明数值分析软件 phase2 是可以比较精准的反映出边坡工程的实际情况 。 （ 4）通过数值模拟软件，我们可以比较清楚直接的看到边坡的整体应力和 应变 等 图，直观的反映出 了边坡的整体破坏情况，为我 们 合理 地 开挖和有效防护提供了有效依据， 体现出数值分析的优越性，是一种经济适用的好方法。
