高速公路深路堑边坡应力及位移的有限元分析

D0I:10.13374/i.issm1001053x.2005.02.036 第27卷第2期 北京科技大学学报 Vol.27 No.2 2005年4月 Journal of University of Science and Technology Beijing Apr.2005 高速公路深路堑边坡应力及位移的有限元分析 乔兰刘宝许李长洪 北京科技大学土木与环境工程学院，北京100083 摘要为了研究山区高速公路深路堑边坡的稳定性问题，基于三维非线性应力及位移有 限元计算，对边坡稳定性进行了数值分析，并与Sarma法的计算结果进行了对比.结果表明， 有限元法解决边坡问题可以起到其他方法不可替代的重要作用 关键词深路堑，有限元；稳定系数；弹塑性分析；应力；位移 分类号U416.14 山区高速公路岩质边坡稳定性的分析方法 行了分析研究，大致确定出边坡的薄弱环节，进 可概括为工程地质类比法、结构分析法和力学计 而获得边坡稳定状态的全面评价，所得结果不仅 算法三大类.工程地质类比法是一种经验方法， 为京秦高速公路工程，也为同类山区高速公路的 在中小型边坡工程的设计中是很通用的一种方 深路堑边坡工程的设计提供参考. 法.结构分析方法是在综合分析坡体中起控制作 用的结构面或软弱带的空间组合特征的基础上， 1基本原理 可以粗略地确定边坡可能的变形破坏形式，极限 平衡法是目前广泛应用的力学分析方法，但在极 有限元法分析边坡稳定性问题是基于应力 限平衡法别中，滑动面形状是根据具体情况假定 场分析四的一种方法，能够全面分析岩体边坡的 的，而不是按塑性理论计算的，岩体作为刚体处 受力状态，其控制方程为： 理，不能反映岩体内的真实应力一应变关系.稳 [K☒{o}={R. 定系数是滑动面上的平均值，带有一定的假定 式中，闪为整体刚度矩阵，{δ}为整体结点位移， 性，也无法考虑累积破坏对稳定的影响.近20年 {R}为整体等效荷载. 来，有限单元法开始应用于研究边坡稳定问题， 对于岩质边坡其力学行为是非线性的，因 并发展为一种强有力的分析计算工具.岩土工程 此，刚度矩阵[)与材料的变形状态和应力状态 有限元分析的日益成熟为定量研究边坡工程中 有关，方程[闪{)={R}属于非线性方程. 的应力场、位移场和稳定性评价提供了强有力的 计算使用的是为岩土工程应用而开发的三 手段 维连续介质有限单元法计算程序3DG,它主要 北京至秦皂岛高速公路青龙连接线所经区 模拟岩土工程结构在三维应力、应变状态条件下 域地处燕山山脉东段，属山岭重丘区，穿越的地 的力学行为及动态变化规律.在该边坡稳定性计 形、地貌及地质条件环境复杂，其中某路段属于 算分析中，建立三维有限元模型，采用非线性弹 大挖方路段，最终将形成近50m的深路堑岩质 塑性模型，经单元划分，形成计算网格，共计8372 高边坡.为了研究公路边坡在开挖完成后的长期 个等参单元.由有限元计算结果可求出单元中点 稳定性问题，本文在有限元分析计算的基础上， 的应力，5，m从而计算出单元上的法向应力和 对边坡中的重点部位的关键点在开挖过程中以 切向应力.则每个单元上的安全系数S等于岩体 及整个边坡形成后的应力和位移的变化规律进 抗剪强度t,和岩体所受的最大剪应力值xm的比 值，由莫尔一库仑准则可得： 收稿日期：20040429修回日期：200407-19 S:=Citoutand (1) 基金项目：国家自然科学基金资助课题①No.50074002) 作者简介：乔兰(1963-一)，教授，博士生导师 式中，S,C,,中，分别是该单元上的安全系数

第 2 7 卷 第 2 期 2 0 0 5 年 4 月 北 京 科 技 大 学 学 报 Jo u r n a l o f U n iv e r s ity o f S e i e n e e a n d Te e h n o l o gy B e ij in g V b L 2 7 N o. 2 A P r. 2 00 5 高速公路深路堑边坡应 力及位移的有限元分析 乔 兰 刘 宝 许 李长 洪 北 京科 技大 学土 木与 环境工程 学院 , 北 京 10 0 0 83 摘 要 为 了研 究 山 区 高速公路深 路 堑边 坡 的稳 定性 问题 , 基 于三 维非 线性 应 力及位 移 有 限 元计算 , 对边 坡 稳定 性进 行 了数值 分析 , 并 与 S ~ a 法 的计 算 结果 进行 了对 比 . 结果 表 明 , 有 限元法 解 决边坡 问题 可 以起 到其他 方 法不 可替 代 的重 要作 用 . 关键 词 深路 堑 ; 有 限元 ; 稳 定系数 ; 弹 塑性 分 析; 应 力 ; 位移 分类 号 U 41 6 . l 滩 山 区高速 公 路 岩 质 边 坡 稳 定 性 的 分 析 方 法 可 概括 为 工程 地质 类 比 法 、 结 构分 析法 和 力学 计 算法 三 大类 〔1 . 工 程地 质类 比法是 一种 经验 方法 , 在 中小 型 边 坡工 程 的设 计 中是很 通 用 的一 种 方 法 . 结构 分析 方法 是在 综合 分 析坡 体 中起控制 作 用 的 结构面 或 软弱 带 的空 间组 合特 征 的基础 上 , 可 以粗 略地 确 定边 坡可 能 的变 形破 坏形 式 . 极 限 平衡 法 是 目前广 泛 应用 的力 学 分析 方法 . 但 在 极 限平 衡 法阴 中 , 滑 动 面 形状 是根据具 体情 况 假 定 的 , 而 不是 按 塑性 理论 计 算的 . 岩 体作 为刚 体 处 理 , 不 能反 映岩 体 内的真 实 应力 一应变 关 系 . 稳 定 系 数 是滑 动 面 上 的 平 均值 , 带 有 一 定 的 假 定 性 , 也无 法考 虑 累积 破 坏对 稳 定 的影 响 . 近 20 年 来 , 有 限单 元法 开 始 应用 于 研 究边 坡稳 定 问题 , 并发 展 为一 种 强有 力 的分析 计算 工具 . 岩土 工程 有 限 元 分析 的 日益 成 熟 为 定量 研 究 边 坡 工程 中 的应 力 场 、 位 移 场 和稳 定性 评价 提供 了强有 力 的 手 段日 一 61 北 京 至 秦皇 岛 高速 公路 青 龙 连 接 线 所 经 区 域地 处 燕 山 山 脉 东段 , 属 山 岭重 丘 区 , 穿 越的地 形 、 地 貌 及地 质 条件 环 境 复杂 , 其 中某 路 段 属于 大挖 方 路 段 , 最 终将 形 成近 50 m 的深 路 堑岩 质 高边 坡 , 为 了研 究 公路边 坡 在开 挖完 成后 的长期 稳定 性 问题 , 本 文在有 限元 分析 计算 的基 础 上 , 对边 坡 中 的重 点 部位 的关 键 点在 开挖 过程 中 以 及 整 个边 坡 形 成 后 的应 力和 位 移 的 变 化规 律 进 行 了分析 研究 , 大致 确 定 出边 坡 的薄 弱环 节 , 进 而 获 得边坡稳 定状 态 的全 面 评价 , 所 得 结 果不 仅 为 京秦 高速 公路工 程 , 也为 同类 山区 高速 公路 的 深 路 堑边坡 工程 的设 计 提供 参 考 , 收稿 日期 : 20 0 4 一0 4 es 2 9 修 回 日期 : 2 0 04一7 一 19 基金 项 目 : 国家 自然 科学 基金 资助课 题困 。 . 5 0 74 0 2) 作者 简介 : 乔兰 ( 1 9 6 3es ) , 教 授 , 博 士 生 导师 1 基 本 原 理 有 限元 法 分 析 边 坡 稳 定 性 问题 是 基 于 应力 场 分析 `7] 的一 种 方法 , 能够 全 面 分 析岩 体 边 坡 的 受 力状 态 , 其控 制 方 程 为 : 因 {咨} 二 { R } . 式 中 , 因 为整 体 刚 度矩 阵 , {升 为 整体 结 点位 移 , {R } 为整 体等 效 荷 载 . 对 于 岩 质 边 坡 其 力 学 行 为是 非 线 性 的 , 因 此 , 刚度矩 阵因 与 材料 的变 形 状态 和 应 力状 态 有 关 , 方 程因 {升二扭 }属 于 非线 性 方程 . 计 算 使 用 的 是 为 岩 土 工 程 应 用 而 开 发 的三 维连 续 介 质 有 限单 元 法 计算 程 序 3 D 、 , 它主 要 模 拟 岩土 工程 结 构在 三 维应 力 、 应 变状 态 条件 下 的力 学行 为 及动 态变化 规律 . 在 该 边坡 稳 定性 计 算分 析 中 , 建立 三 维有 限元 模 型 , 采 用 非线 性 弹 塑性 模型 , 经 单 元划 分 , 形成 计 算 网 格 , 共 计 8 3 72 个 等 参单 元 . 由有 限元 计算 结果 可 求 出单 元 中点 的应 力 氏氏几 , 从 而计 算 出单元 上 的法 向应 力 和 切 向应 力 . 则 每个 单 元 上 的安全 系 数 SF 等 于岩 体 抗 剪 强度 几 和 岩 体所 受 的 最 大剪 应 力 值瑞ax 的 比 值 , 由莫 尔 一 库仑准 则【别 可 得 : 凡 = c 汁丙沐a n 功 ` ( l ) 式 中 , 凡 , 。 , 丙 。 味i , 价 , 分 别 是该 单 元上 的安全系 数 、 DOI: 10. 13374 /j . issn1001 -053x. 2005. 02. 036

Vol.27 No.2 乔兰等：高速公路深路堑边坡应力位移有限元分析 ·143 内聚力、法向应力、切向应力以及内摩擦角.将最 算，虽然在有水的状态，边坡的稳定性要受到一 大、最小主应力和代入上式并对整个滑动面 定影响，但边坡的安全系数均大于安全规程规定 进行加权处理，得到整个滑动面的安全系数为： 的13~1.5，因此所分析的边坡是稳定的. cco吟p S=五= 2.2非线性有限元法 (2) 4-9 (1)本构模型. 2 -2c:cos2-(a+☑)sinp 在计算分析中，采用非线性弹塑性模型，， 01G3 用常规方法无法考虑边坡介质的塑性和蠕变性 当S=1时，单元处于极限状态；当S>1时，单 质，必须采用有限元法进行分析，用该模型可得 元处于稳定状态：当S

V6 1.2 7 N O . 2 乔 兰等 : 高速 公路 深路 堑边 坡应 力位 移 有限 元分 析 内聚力 、 法 向应 力 、 切 向应 力 以及 内摩擦 角 . 将最 大 、 最 小主 应 力 , l和 仍 代 入上 式 并对 整 个 滑动 面 进 行加 权处 理 , 得 到整 个 滑动 面 的 安全 系 数为 : 凡 = 兰 几ax al + 仍 . 、 c ’ ` U 即一 丁 一 sl qn, 口 1一 氏 ( 2 ) 2 = Z c · c o呻一 ( 。 ,+ 仍 ) 5 1呻 吓 1一 氏 当吕二 l 时 , 单 元 处于 极 限状 态 ; 当凡>1 时 , 单 元处 于稳 定状 态 ; 当凡l< 时 , 单 元进 入 塑性 状态 . 工 程 范 围涉及 的岩 体为 大理 岩 和花 岗岩 , 属 弹塑 性材 料 , 适用 于 莫 尔一 库 仑破 坏 准 则 : 呼 一 粤 S ionp +c cop 、 一 :一糯 一 c2 丫恶黯 ( 3 ) ( 4 ) 式 中 , a , , 伪 分 别是 最 大和 最 小主 应 力 ; 。 , 沪分别 是 材 料 内聚 力和 内摩擦 角 不 是 破坏 判 断系 数 , 当厂 ) 0 时 , 材料 处 于 塑性 流 动状 态 ; 当不< O 时 , 材料 处 于 弹性 变 形 阶段 . 在拉 应 力状 态 下 , 如果 拉应 力 超过 材 料 的抗 拉 强度 , 材料 发 生 拉破 坏 . 算 , 虽 然 在有 水 的状 态 , 边坡 的稳定 性 要受 到 一 定影 响 , 但 边坡 的安 全系 数均 大 于安全 规 程规 定 的 1 . 3一 1 . 5 , 因此 所 分 析 的边 坡 是稳 定 的 . .2 2 非 线性 有 限 元法 ( l) 本 构模 型 . 在 计 算分 析 中 , 采 用 非 线性 弹 塑性 模 型`,’ , , Z J , 用 常规 方 法 无法 考 虑 边 坡 介质 的 塑性 和 蠕 变性 质 , 必 须采 用 有 限元 法 进行 分 析 , 用 该模 型 可得 到边坡 体 内的应 力场 、 位移 场 以及破 损 区 发育和 演 化 规律 . 塑 性 状 态 下 的应 力 和 应 变 只 能 建 立应 力增 量 与应 变增 量 之 间 的关系 . 当应 力产 生 一无 限 小 增 量 时 , 假 设 应变 的变 化可 分 成弹 性和 塑性 两个 部 分 : ds , = d 。笋d搏 ( 5) 弹 性 应 力 增 量 与 弹 性 应 变增 量 之 间可 由弹 性 矩 阵D 联系 , 塑 性应 变增 量 由塑 性势理 论 给 出 . 对 弹 塑性 介 质 , 存在 塑 性 势 函 数 Q , 它 是 应 力状 态 和塑 性 应变 的 函数 , 使 得 d` 一 、织 . 口氏 2 计 算 结 果 对 比 京 秦 高速 公路 青 龙 连 接 线某 段 为深 挖 高填 工程 , 该 段工 程 完 工后 形 成近 50 m 的深 路 堑边 坡 , 工 程 范 围涉 及 的主 要 岩 体 为 大 理岩 和 花 岗 岩 , 属 弹 塑性 材 料 . 下 面 用 非线 性 有 限元 和极 限 平衡 分 析方 法 — S ~ a 法对 边 坡 稳 定状 态分 别 进 行 分 析 . 2 . 1 S a r m a 法 S ~ a 法可 以用 于评 价各 种 类 型滑 坡 的稳 定 性 , 如 平 面滑 动 、 楔 体 滑动 、 圆弧滑 动 和非 圆弧 滑 动 等各 种 复杂剖面 岩土 斜坡 , 且它 无 需条 块边 界 垂 直 , 即无 须 垂直条 分滑 体 , 从而 可 以对各 种 特 殊 的边 坡 破 坏结 构 进行 稳 定 性分 析 ` 90J] . 本 文 用 E v ert H oe k 编 制 的非 垂直 分 条 稳定 性 分析 S ~ a 程 序进 行边坡稳 定 性分 析 . 计算 中采用 的岩 石 力 学 参数 为 : 密度 p二 2 8 90 k g n/ 13 , 纵 波速 度 v 二 6 52 . 09 m s/ , 抗 压 强度 =ac 2 . 82 M P a , 抗拉 强 度几 = .0 9 80 M P a , 弹性模 量于9 . 3 4 9 G P a , 动切 变模 量 G 二 1 1 . 8 8 8 GP a ,泊 松 比尸= 0 2 8 3 , 内聚 力 e 二 0 3 3 17 M p a , 内摩 擦 角价= 4 30 . 用 S a n l l a 法 计 算 时分别 考 虑 边坡 无水 和 边坡 充满 水两 种情 况 , 计 算结果 是无 水 时边坡 安全 系 数 为 3 . 39 7 , 有水 时安全 系数 为 3 . 2 86 . 经过 分析 计 式 中 从是 一 正 的有 限量 , 称 塑 性系 数 , 它 的具 体 数值 和 材 料 硬 化 有 关 . 对 于 稳 定 的应变 硬 化 材 料 , Q 取与 后 继 屈 服 函 数 F 相 同的形 式 . 当 Q 二 F 时 , 上 式可 表 示 为 : 刁F a E二二 几, 苏一一 . 口 丙 则 总应变 增 量 可表 示 为 : 螂 二 D 一 ’ d氏 . 。 DF 十说 不犷一 0 丙 由一 致 性条 件 可推 出塑 性系 数 : 1 a F , 又二 ~ 冬 蔫二 - d a ` 八 一 (6 、 A 口氏 ~ 歹 甲少 对 理想 弹 塑 性材 料 , 无 应变 硬 化 时A = 0 ; 对 于 加 式硬 化 材料 , 采用 塑 性 功加 式硬 化 定律 时 , 日F ~ 日O 己F 己O A = 一决 土刀 D丢兰 es 一 书舟 口 决兰 ` 几 一 a 峪口 D确 a u 口 犷 a 丙 式 中 , u 为 塑 性功 , 则应 力 应 变关 系表 达 式 为 : ` 一 ’ 1 初 即 、 · _ 、 ds 办 二 (D 一 , 十 吃一 若岑 ` -攀 一 )d氏 (7 ) 一 “ 甲 A 口丙 d 丙 j 一 “ (2 ) 计 算模 型 . 根据 的现场 实 际情 况 , 并 加 以适 当简 化后 建 立三 维 计算 模型 , 该边 坡水 平方 向沿 公 路走 向长 18 0 m , 开 挖部 分 长 80 m , 垂 直 公 路走 向方 向取 15 3 m , 铅 直 方 向取 98 m . 经 单 元划 分 , 形成 计 算 网格 , 共计 8 3 7 2 个 等参 单元 , 36 5 87 个节 点 . 边 坡 剖面 的网格 划 分 见 图 1 , 图 2 为 开挖 完 成后 剖 面 上 的 网格情 况

·144 北京科技大学学报 2005年第2期 (a) 图1边坡开挖前网格划分图 Fig.1 Meshing of the slope before excavation (b) (c) 图2边坡开挖结束的网格划分图 Fig.1 Meshing of the slope after excavation 计算采用非线性弹塑性模型，计算材料参数 为：弹性模量E=9349MPa,泊松比4=0.283，容重 图3计算结果.(a)最大主应力等值线（单位：100kPa:(b)最小 y=0.0289MNm,内聚力c=3.317MPa,内摩擦角 主应力等值线（单位：100kPa:(c)水平位移等值线（单位：mm) 中=43°，抗拉强度a。=0.98MPa. Fig.3 Calculated results.(a)Isoclines of the maximal stress (unit: 100 kPa);(b)Isoclines of the minimum stress (unit:100 kPa);(c) (3)计算结果分析· Isoclines of horizontal displacement (unit:mm) 计算过程中按以下步骤进行：首先形成初始 2.3方法比较 应力场，然后从边坡最高处开始每6m一个台阶 进行开挖至设计平面，计算结束.计算得到的边 比较两种计算方法的计算结果可以看出：有 限元法与Sarma法计算山区高速公路深路堑边坡 坡最大主应力、最小主应力和位移见图3. 的稳定性得出的结论大致相同，显然有限元法的 由图可见，最大主应力在坡角处出现应力集 计算结果精确度更高一些，原因就是有限元法更 中现象，随着高程的增加逐渐减小，而在上部的 加真实地反映了岩土体的应力一应变关系，计算 各个台阶上由于卸载作用，其应力值要比同一水 平的岩体内部的应力偏低，甚至在台阶的平台上 模型不仅满足力的平衡条件，而且满足边坡岩土 有拉应力产生，但其数值最大为0.39MPa,低于 体的应力一应变关系，并且可以对边坡的非线性 弹塑性进行分析，所以计算结果精确度更高，更 岩体的抗拉强度.在每个台阶的内部，最小主应 可靠 力都有一个拉应力区，范围在距离坡面2~3m内， 但拉应力的数值较小，最大为0.17MPa,远低于 3实际验证 岩体的抗拉强度，因而不会产生受拉破坏.岩体 在x方向的位移分布情况，从趋势上看，在最终开 该段高速最终将建成了一段47.6m的深路 挖的平面上，位移最大，因而是最危险的部分.但 堑岩质高边坡，从建成到现在三年多的运营情况 从数值上看，最大的位移只有0.78mm,是一个非 来看，该边坡稳定状态保持良好，未发生滑坡和 常小的数值，因而也不会产生破坏.在整个计算 坍塌现象，尤其是经过三个雨季山洪多发期和两 范围内，没有产生塑性区，表明边坡是稳定的. 个冬季冻融期的检验后，该边坡仍保持较好的稳

北 京 科 技 大 学 学 报 2 0 0 5 年 第 2 期 ùōóó一ùō ōù一ùōōóù 一熟一 F ig , 1 图 1 边坡 开 挖 前 网 格划 分图 M e s h i . g o f t h e s l o P e b e lb r e xe e av at i o n 撇皿咖óó姗ù哪踱洲叫卿 图 2 边 坡开 挖结 束 的 网 格划 分 图 iF g · 1 M e s h i n g o f ht e s l o p e a ft e r e x e a v a t i o n 计算 采用 非 线 性弹 塑 性模 型 , 计 算材 料 参数 为 : 弹 性模 量 E 二 9 3 4 9 MaP , 泊松 比声= .0 2 83 , 容重 少二 0 . 02 8 9 M N 八n 3 , 内聚 力 c = 3 . 3 17 M P a , 内摩擦 角 砂= 43 “ , 抗 拉 强度 几 = 0 . 98 M aP . (3 ) 计 算 结 果分 析 . 计算 过程 中按 以 下 步骤 进行 : 首 先形 成初 始 应力 场 , 然后 从 边坡 最 高 处开 始 每 6 m 一 个 台 阶 进行 开 挖 至设 计 平面 , 计算 结束 . 计算 得 到 的边 坡最 大 主应 力 、 最 小主 应 力和 位移 见 图 3 . 由图可见 , 最大 主 应力 在 坡角 处 出现 应 力集 中现 象 , 随着 高 程 的增 加逐 渐 减 小 . 而 在 上部 的 各个 台阶 上 由于 卸载 作用 , 其 应 力值要 比 同一 水 平 的岩体 内部 的应力 偏低 , 甚至 在 台阶 的平 台上 有拉 应 力产 生 , 但 其 数值 最 大 为 .0 39 M p a , 低 于 岩 体 的抗 拉强度 . 在 每 个 台 阶 的 内部 , 最 小主 应 力都 有一 个拉应 力 区 , 范 围在距 离坡 面 2一3 m 内 , 但 拉 应 力 的数 值较小 , 最 大 为 0 . 17 M P a , 远 低 于 岩 体 的抗 拉 强 度 , 因而 不会 产 生 受拉 破 坏 . 岩 体 在 x 方 向的位移 分 布情 况 , 从趋 势 上看 , 在 最 终开 挖 的平面 上 , 位移 最大 , 因而 是最 危 险的部 分 . 但 从 数 值上 看 , 最 大 的位 移 只有 0 . 78 11 1和 , 是 一个 非 常 小的数 值 , 因而也 不 会产 生破坏 . 在 整个 计 算 范 围 内 , 没有 产 生塑 性 区 , 表 明边 坡 是 稳定 的 . 图 3 计 算结 果 . a() 最大 主应 力等 值线 (单位 : 10 kP a) ; (b )最 小 主应 力等值 线 (单位 : 10 kP a) ; c() 水平 位移 等值 线 (单位 : m m ) F ig . 3 C a l e u l a t e d 代s u l t s . ( a ) I s o e l i n e s o f t h e 爪a x i m a l s t概 s (u n i t : 1 0 0 kP a ) ; ( b ) I s o c il n es o f t h e m i n i m u m s t er s s ( u n it : 1 0 0 kLP a ) ; c( ) I s o e il n e s o f h o山 o n t a l d坛 P l a e e m e n t ( u n i t : m m ) 2 . 3 方 法 比较 比 较 两种 计 算方法 的计算结 果 可 以 看 出 : 有 限元 法与 S a n 们 a 法 计算 山 区 高速 公路 深路 堑边 坡 的稳 定性 得 出的结论 大 致相 同 , 显 然有 限元法 的 计 算结果 精 确度 更 高一 些 , 原 因就 是 有 限元法 更 加真 实地 反 映 了岩土 体 的应 力一应 变 关系 , 计 算 模 型 不仅 满足 力 的平 衡条件 , 而且 满足边 坡岩 土 体 的应 力一 应变关系 , 并且 可 以对边坡 的非线 性 弹塑 性进 行 分 析 , 所 以计 算结 果精 确 度 更高 , 更 可靠 . 3 实际 验 证 该段 高速最 终 将 建成 了一 段 47 . 6 m 的深 路 堑岩质 高 边坡 . 从建 成 到现 在三 年 多 的运 营情 况 来 看 , 该边坡 稳 定状 态 保持 良好 , 未 发 生滑 坡 和 坍 塌现 象 , 尤 其 是经 过三 个 雨季 山洪 多发 期和 两 个 冬季 冻 融期 的检 验 后 , 该 边坡 仍保 持 较好 的稳

Vol.27N0.2 乔兰等：高速公路深路堑边坡应力位移有限元分析 ·145· 定状态，目前边坡的植被长势很好，也为该边坡 法具有传统极限平衡方法不能做到的优越性，值 保持长期稳定性发挥了重要作用. 得在工程中推广使用， 4结论 参考文献 [1】邵龙潭，堂洪祥，韩国城.有限元边坡稳定分析方法及其 (1)计算边坡岩体范围内没有出现塑性区，即 应用.计算力学方法，2001,21(1)：81 在岩体内部没有出现塑性破坏，边坡是稳定的. [2林士炎，李长洪，乔兰，等.爆破震动对高速公路边坡影响 (2)最大主应力在边坡的坡角处出现应力集 的数值模拟.北京科技大学学报，2003,25(6：507 中，位移也在最终开采平面为最大，所以在最终 [3]郑颖人，赵尚毅，时卫民，等.边坡稳定分析的一些进展 地下空间，21(4)：262 形成路基的平台的坡角处是整个边坡稳定性最 [4]汪益敏.有限元法在边坡岩体稳定分析中的应用.西安公 薄弱的部位.但这个位置的应力和位移均在允许 路学院学报，1994,14(2)：13 范围内，不会产生破坏. [5)曾开华，陆兆臻.边坡变形破坏预测的混沌与分形研究.河 海大学学报，1999,27(3)：10 (3)有限元法可以模拟计算边坡介质的塑性 句吴梦军.极限分析上限法在公路边坡稳定分析中的应用 和蠕变特性，通过非线性有限元计算，可以得到 重庆交通学院学报，2002,21(3：52 边坡内任一点的位移和应力状态，将传统的以某 [刀吴继敏，孙少锐，陆小敏.台阶状结构面岩体边坡非线性 一平面来计算确定边坡的稳定性系数深化到以 有限元分析.红水河，2002,21(2)：48 ⑧]时卫民，郑颗人，张鲁豫.岩石高边坡的有限元分析及其 一点来计算它的安全系数，使问题得到更加精确 简化分析方法.地下空间，2001,21(5)：455 的解答， [9)孙玉科，牟会庞，嬈宝魁边坡稳定性分析，北京：科学出版 (4)本文的研究实例表明，对于山区高速公路 社，1999 深路堑岩质边坡的稳定性问题，可以用有限元的 [10)周萃英，滑坡灾害系统的自组织.中国地质大学学报，1996， 21(6):34 方法进行分析研究.基于有限元法的计算结果进 [1]蔡美锋，何满朝，刘东燕。岩石力学与工程.北京：学出版 行深路堑岩质边坡稳定性分析是切实可行的，能 社，2002.180 够符合实际的判别边坡稳定状态，可用应力、位 [12]夏琼，杨永海.利用有限元法分析路是边坡稳定性.甘肃 移等指标确定边坡的极限破坏区，同时，有限元 科技，2002,18(12)：48 Stress and deformation analysis of freeway deep cutting slope by finite element method QIAO Lan,LIU Baoxu,LI Changhong Civil and Environmental Engineering School,University of Science and Technology Beijing,Beijing 100083,China ABSTRACT Based on calculations of nonlinear stress-deformation by the 3-directional finite element method, this paper numerically analyzed the stability of a mountainous freeway deep cutting slope.The safety factors obta- ined by this method and Sarma's method were compared with engineering practice.The results show that the finite element method is more effective for analyzing the stability of a slope than traditional methods. KEY WORDS deep cutting;finite element method;stability factor;elasto-plastic analysis;stress;deformation

V b l . 2 7 N o 一 2 乔兰 等 : 高 速公 路深 路堑 边坡 应 力位 移有 限元 分析 1 4 5 定状 态 . 目前边 坡 的植 被 长 势很 好 , 也 为 该边 坡 保持长 期稳 定性 发挥 了重 要 作用 . 法 具有 传 统极 限平 衡方 法不 能做 到 的优越 性 , 值 得 在 工程 中推 广 使 用 . 4 结论 ( l) 计 算边 坡岩 体 范 围 内没 有 出现塑 性 区 , 即 在 岩 体 内部 没有 出现 塑性 破 坏 , 边 坡是 稳 定 的 . (2 ) 最 大 主应 力在 边 坡 的坡 角 处 出现 应 力集 中 , 位移 也 在最 终 开采 平 面 为最 大 , 所 以在 最 终 形 成 路基 的平 台 的坡 角 处是 整 个 边 坡 稳 定性 最 薄 弱 的部位 . 但这 个位 置 的应 力和 位移 均 在 允许 范 围 内 , 不会 产 生破 坏 . (3 ) 有 限元 法 可 以模 拟计 算 边坡 介 质 的塑 性 和蠕 变特性 , 通过 非 线 性有 限元计 算 , 可 以得 到 边坡 内任 一 点 的位 移和 应 力状 态 , 将传 统 的 以某 一平 面 来 计 算确 定边 坡 的稳 定 性 系数 深 化 到 以 一 点来计 算它 的安 全系 数 , 使 问题 得到 更 加精 确 的解 答 . (4 )本 文 的研究 实例 表 明 , 对 于 山 区 高速 公 路 深 路堑 岩质 边 坡 的稳 定 性 问题 , 可 以用 有 限元 的 方 法进 行分 析研 究 . 基 于有 限元 法 的计 算结 果进 行 深路 堑岩 质边 坡稳 定 性分 析是 切 实可行 的 , 能 够 符合 实 际 的判 别边 坡 稳 定状 态 , 可用 应 力 、 位 移 等指 标确 定边 坡 的极 限破 坏 区 . 同 时 , 有 限元 参 考 文 献 【1] 召肠龙潭 , 堂洪 祥 , 韩 国城 . 有 限元边坡 稳定 分析 方法及 其 应用 . 计 算力 学方法 , 2 0 0 1 , 2 1 ( 1) : 8 1 2[ ] 林 士 炎 , 李长洪 , 乔 兰 , 等 . 爆破 震动 对高速 公 路边 坡影 响 的数值 模拟 . 北 京科技 大学 学报 , 2 0 0 3 , 2 5( 6 ) : 50 7 价] 郑颖 人 , 赵 尚毅 , 时卫 民 , 等 . 边坡 稳定 分析 的一些 进展 . 地 下 空 间 , 21 (4 ) : 2 62 4[ l 汪益 敏 . 有 限元法 在边 坡岩体 稳定 分析 中的应用 . 西安 公 路 学院 学报 , 19 9 4 , 14仅) : 13 「5] 曾开华 , 陆兆 臻 . 边坡 变形破 坏预测 的混 沌与分 形研究 . 河 海大 学 学报 , 1 99 9 , 2 7 ( 3) : 10 师] 吴梦 军 . 极限分 析上 限法在 公 路边坡 稳定 分析 中的应 用 . 重庆 交通 学院 学报 , 2 0 02 , 2 1( 3) : 5 2 「71 吴继 敏 , 孙少 锐 , 陆 小敏 . 台阶状结 构面岩 体 边坡 非线性 有 限元分 析 . 红水 河 , 2 0 2 , 21 ( 2 ) : 48 8[ ] 时卫 民 , 郑颖 人 , 张鲁 豫 . 岩石 高边 坡 的有 限元 分析及 其 简化 分析方 法 , 地下 空 间 , 2 0 01 , 2 1( 5) : 4 5 9[ ] 孙 玉科 , 牟 会庞 , 姚宝 魁 . 边坡 稳定性 分析 . 北京 : 科学 出版 社 , 19 9 9 【10 」周 萃英 . 滑 坡灾害 系统 的 自组织 . 中国 地质大 学学报 , 19 9 6, 2 1 ( 6 ) : 3 4 【1 ]蔡美 锋 , 何满 朝 , 刘东燕 . 岩 石力 学与工程 . 北 京 : 学 出 版 社 , 2 0 0 2 . 1 80 【12 ] 夏 琼 , 杨 永海 . 利 用有 限元法 分析 路堤 边坡稳 定性 . 甘 肃 科 技 , 2 0 0 2 , 18 ( 12 ) : 4 8 S tr e s s an d d e fo rm at i o n an a ly s i s o f fr e e w ay d e e P c ut i n g s l o P e b y if n it e e l e m e n t m e t h o d 口侧 口 L a n, LI U B ao x 玩 lL hC a喇 o gn C i v i l an d E vn i r o曲 e nt a 1 E l l g i n e e r in g S hc o o l , nU iV er s ity o f S e i e n e e an d eT e 恤 o l o gy B e ij in g , B e ij i n g 10 0 0 8 3 , Ch in a A B S T R A C T B a s e d o n e al c u l iat on s o f n o n li n e ar s tre s s 一 d e of rm at i o n b y ht e 3 一 d i r e ct i o n a l if in t e e l e m e nt m e t h o d, iht s Pap er imnt e ir c al ly an a ly z e d ht e s t ab i li yt o f a m o l l n t aln o u s fr e e w ay d e e P e ut i gn s l o P e . hT e s afe yt af ct o r s o b t a - in e d 勿 iht s m e ht o d an d S a r Fn a ’ 5 m e th o d w e r e e o m Par e d w iht e n g i n e ier gn P r a c t i c e . hT e r e s ult s s ho w ht at ht e if n i t e e l em en t m e ht o d 1 5 m o r e e fe e ti v e of r an aly z ign ht e s t ab il iyt o f a s l op e ht an tr ad iit o n a l m e ht o d s . K E Y W O R D S d e e P e ut in g: if in te e l e m e nt m e hot d: s t ab i li yt acf to r: e l a s t o 一 Pl a s it e an a ly s i s : s tr e s s : d e of mr at i o n