深埋大跨度地下硐室动态开挖支护的数值模拟与稳定性分析论文.pdfVIP

深埋大跨度地下硐室动态开挖支护的 数值模拟与稳定性分析 张巨伟1 柳郁2 (1．中国有色金属矿产地质调查中心2．北京科技大学土木与环境工程学院) 摘要本文采用了大变形有限差分数值分析软件FLAC，对某地下深埋大跨度硐室的分部开挖施工 过程和喷锚支护工艺进行了数值模拟，在此基础上，揭示了大跨度硐室动态开挖过程中围岩 的稳定性，为开挖和支护设计提供了理论依据。 关键词 地下工程深埋大跨度硐室开挖数值模拟稳定性分析 1 前言 理论分析和工程实践证明，地下硐室开挖的破坏方式主要有两种类型：一种为围岩整体 性失稳，这种失稳破坏主要是应力重新调整形成的次生应力超过围岩强度，使得硐室周围出 现压剪破坏区、塑性区或拉力区，如不及时控制变形破坏的发展，围岩将发生整体性失稳破 坏。另一种是围岩局部失稳破坏，该破坏方式是由不利的软弱结构面组合形成不稳定块体的 滑移或冒落[1]。无论哪种形式，硐室的稳定性与硐室开挖方式与施工工艺密切相关。因此， 借助于数值分析方法，对硐室的施工工艺和支护参数进行动态模拟，由此优化施工工艺和支 护参数，对地下工程的安全、经济开发与建设是十分重要的。 地下硐室开挖过程的稳定性评价，通常涉及围岩的失稳判据、失稳范围和稳定计算方 法。一般情况下，对于整体失稳破坏的硐室，硐室围岩的强度准则和变形量值常用作为围岩 稳定的判别条件。为此，首先确定硐室在开挖过程中的围岩和结构的应力分布以及基于强度 准则的围岩屈服破坏范围和变形特征；然后，基于围岩的塑性区范围和变形量值进行硐室的 稳定性评价；最后，根据不同开挖施工工艺的稳定状态的对比分析，进行硐室开挖工艺优化 和支护设计。 本文针对某水电站地下厂房硐室的开挖与支护的施工工艺，采用数值分析方法，使用 FLAC软件模拟硐室分部开挖过程，并分析围岩的应力状态和受损程度以及支护效果。 2工程背景及工程评价 道方向N60。E，与硐室轴线在电厂以南相交，进出隧道的入口位于泻水口地方。两上勘探隧 道从进出隧道分叉，第一条勘探隧道在硐室w／sw约35m地方；第二条通过硐室的拱顶。 这两条的方向大致平行于硐室轴线。硐室所处围岩主要为某花岗岩岩盘的一部分。该花岗岩 盘的组成是均匀中粒花岗闪长岩，含有许多石英和细晶岩脉。根据现场调查和初步分析，预 测其围岩属于IlI类围岩。 13】 3动态开挖模拟数值模型与参数 3．1计算模型 考虑到该硐室跨度为20．0m，高度为37．5m， 故其计算范围取150m×150m，共分为35×35= 1225个单元。计算模型网格如图1所示。 3．2边界条件 为了模拟水平构造应力场，采取赋初值对 原岩进行赋相应的应力场。采用应力边界条件， 模型的上表面施加均匀的垂直压应力，模型两 侧面施加随深度变化的水平压应力，模型下表 面垂直和水平位移固定。模型所有岩层均采用 摩尔一库仑模型。 原岩应力的垂直应力取自重应力，而水平 图1计算模型网格 应力取侧压系数A=1．1。原岩应力的边界条件为： 5．929MPa。 3．3模拟支护方案 初期支护采用5m短密锚杆，二次支护采用25m长锚索，在模拟时考虑支护时机对硐室 稳定性的影响。 3．4计算参数 (1)原岩岩体计算参数 比弘20．28。 (2)围岩计算参数 0．6MPa，凝聚力c=0．6MPa，内摩擦角∞=350． (3)锚杆(锚索)参数 为了进行锚杆、锚索和喷层的数值模拟，需采用相应的结构单元加以模拟。在FLAC 程序中，有专门开发的锚杆单元，可以模拟全长粘结锚杆或预应力锚索。加钢筋网的喷层可 采用梁单元。 根据锚杆的直径、锚固方式(端锚)、间排距等确定相应的计算参数。全长砂浆锚杆的 计算参数列于表1中，预应力锚索计算参数如表2所示，梁单元的计算参数由表3给出。 锚杆单元计算参数 表1 直径 截面积 排距 灰浆抗压强度 弹性模量 屈服