考虑三维应力旋转隧洞衬砌支护时机的研究.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 考虑三维应力旋转的隧洞衬砌支护时机研究 张社荣 梁礼绘 天津大学建筑工程学院 天津 摘要 基于岩体裂隙的扩展机理对隧洞开挖过程进行数值模拟 通过对模型中计算点岩体的应力路径分析表明 随 着隧洞挖掘的深入 其顶部岩体的第三主应力方向由竖直向下旋转至沿洞壁的切线方向 隧洞边墙岩体的第三主 应力方向一直维持不变 当掌子面推进至超过计算点 左右时 顶部计算点岩体的应力旋转趋于稳定 对隧洞 顶部计算点岩体的位移路径分析表明 沿隧洞纵剖面观察 其顶部岩体的位移方向与竖直方向 向下为正 的夹角 随掌子面地推进逆时针从 变化至 同样也是在掌子面推进至超过计算点 左右时位移方向变化趋于稳定 因此应该在此时对围岩进行衬砌支护 关键词 隧洞开挖 数值模拟 应力路径 衬砌支护 应力旋转 中图分类号 文献标识码 研究背景 隧洞的开挖破坏了岩体应力的原平衡状态 洞周处初始地应力被释放 围岩应力重新分布 隧洞周 围岩石发生形变 随着开挖的不断进行 应力重新分布区域逐渐扩大 围岩的应力状态发生变化 对 围岩稳定造成不利影响 因此必须对隧洞进行合理的支护设计 充分发挥围岩的自承能力是现代支护 理论的一个基本观点 隧洞开挖后应尽早进行一次支护以使围岩与支护协同作用 形成统一的承载体 在选择二次支护时机时 必须考虑两个方面的问题 一是不能让围岩的变形过大而进入松动状态 二是 确保二次支护上的力在合理范围之内 以满足支护的强度要求 要满足这两个方面的要求 必须对隧洞 开挖后围岩的应力状态进行分析 围岩的应力状态包括应力的方向和大小两个方面 应力方向的旋转 常常被人们忽略 多年来 各国学者对最佳支护时机做了大量研究工作 辛全山等 曾利用二次耦合理 论确定二次支护的最佳时间 林勇等 从应变能的角度对合理的衬砌支护时间进行过研究 但都没有考 虑隧洞开挖过程中围岩的应力旋转的影响 主应力的旋转是一个常常被人们忽略的问题 实际上 大量 土工实验表明 纯主应力旋转 三主应力值不变 会导致土体出现明显的塑性变形 因而主应力旋转是 必须考虑的岩土工程的力学问题 本文在已有研究的基础上 对某引水隧洞开挖过程进行三维数值模拟分析 主要对掌子面推进过程 中围岩的应力状态进行分析和研究 在考虑三维应力旋转的基础上进行衬砌支护时机的选择 其结论可 以为相似工程提供参考 基本原理 计算原理 本文采用显式有限差分方法进行数值模拟 根据高斯定理 对于函数 有 收稿日期 作者简介 张社荣 男 天津人 博士 教授 主要从事水工结构静动力分析研究 式中 为函数求解域 或单元 的体积 为 的边界 为 的单位外法线矢量 定义梯度 的平均值为 式中 表示平均值 对于一个具有 条边的多边形 式 可写成对 条边求和的形式 式中 是多边形的边长 是 在 上的平均值 若以速度 代替式 中 且 取边两端节点 即差分网格的角点 和 的速度平均值 则 式中 分别为差分网格的对角线两点 为求和的基数 由几何方程可求得单元的平均应变增量 式中 为差分网格单元的平均尺寸 由广义虎克定律 各向同性材料的本构方程为 式中 为应变 为拉梅常数 为体积应变 因此 单元的平均应力增量可表达为 同时 若以应力表示应变 则其本构关系为 式中 为弹性模量 为应力第一不变量 为泊松比 为应变 这样 通过上述迭代求解 便可求出每一迭代时步相应各单元的应力和应变值 根据不同的塑性屈 服准则 通过相关计算可求出各域 单元 的应力 从而可以判定哪些单元屈服 由节点的运动方程 ? 式中 为密度 为总加速度 为重力加速度 对单元应力表达式沿积分路径积分得 式中 为某节点周围单元作用在该节点的集中力 式中 为质量 指作用在节点的合力 静力 利用中心差分 得某