岩溶地层隧道爆破开挖的数值模拟研究.pdfVIP

·8· 矿山建设工程新进展一2006全国矿山建设学术会议文集(下册) 岩溶地层隧道爆破开挖的数值模拟研究 梁为民 杨小林 战 军 (河南理工大学土木工程学院 河南焦作454000) 摘要：本文以圆梁山深埋特长隧道复合岩层为对象，通过数值计算模拟了不同岩层爆破应力波传播 规律，以及溶洞与隧道空间关系对隧道爆破开挖的影响，指出邻近溶洞一侧隧道的拱脚与直墙连接处为危 险区域，应加强支护。 关键词：岩溶地层；隧道爆破；数值计算 I 引言 068 圆梁山深埋特长隧道全长11 m，是渝怀线的关键性控制工程。隧道进口位于细沙 车站怀化端，车站伸于隧道949m，出口紧邻炭厂河～号大桥。隧道地处渝、鄂、黔三省市毗 连地区，为川东褶皱山地与鄂西山地、贵州高原的接触带，属中、低山地形。主要发育毛坝向 斜、桐麻岭背斜及伴生断裂，向斜区内发育较多横张断裂。地貌形态明显受构造和岩性控制， 具带状展布特征，以褶皱构造为骨架，形成北北东向山脉和纵向河谷相间。地形和地质条件 异常复杂，主要工程地质问题有：岩溶涌(突)水、高水压、高应力及煤层瓦斯。本文以圆梁山 深埋特长隧道为工程背景，对岩溶地层隧道爆破开挖进行了数值模拟，以期能为现场爆破施 工提供指导建议。 2圆梁山岩溶区隧道工程地质概况 圆梁山深埋特长隧道地质构造单元属扬子准地台东南部的上扬子台坳，为川东陷褶束 ～DK357+820为泥盆系(D)、志留系(S)非可溶岩地层构成。 m。 埋段，此段隧道埋深约150～180 3岩溶地层爆破开挖对岩体影响的数值模拟 目前在岩石力学中应用的计算方法和软件众多，目前可用于层状岩体中岩体破裂与岩 层移动的数值计算方法主要包括两类：一类是建立连续介质力学基础上的数值计算方法，诸 如有限单元法、边界元法等；另一类则是通常所说的非连续介质力学分析方法，其中最有代 表性的属离散元法。本文爆破数值模拟程序采用DYNA一2D程序，该程序能处理结构在高 速碰撞和高能炸药爆炸下的动态响应，使得该程序非常适合于进行爆破数值模拟计算。 DYNA2D程序中共包含有11种材料模型，其中对于岩石材料常用的两种材料模型为 弹性和弹塑性材料模型，炸药的材料模型为高能炸药燃烧模型。 (1)各向同性弹性材料 爆破工程 ·9· 输人参数为：弹性模量E，泊松比。 考虑到结构有大变形和大转动，采用Jaumann应变率。 J一畦 d暑=aij—dIka-2ki—d 式中，皿J=丢(薏一差)为旋转张量。对于各向同性材料，其本构关系为 略=C,j盯ek，=[蛾+产(以懿+欧)]％ 2)随动硬化、各向同性硬化及其组合硬化弹塑性 输入参数为：弹性模量E，泊松比￡，，初始屈服极限cr0，切线模量E丁，硬化参数p，(卢一0 时为随动硬化；卢=1时为等向硬化；O卢z时为混合硬化)。 各向同性硬化材料塑性加载时，屈服面中心固定不变而屈服面半径增大，随动硬化材料 塑性加载时屈服面半径不变而屈服面中心沿塑性应变方向平移。对于大多数实际材料，屈服 的硬化规律介于各向同性硬化和随动硬化之间，即为混合硬化材料。 (3)高能炸药燃烧 输入参数：爆速D，C—J压力Pc-J和燃烧系数F。 以上二种材料只有炸药燃烧需要状态方程，常用的为JWL高能炸药，下列状态方程通 常用于描述高能炸药及爆轰产物，其形式为： 户刊(1_南卜y+B(卜南)e_R∥+等 式中，A、B、R，、R：、1,O为输人参数。 3．1计算模型和参数 由于岩溶地区多为层状岩体，因此在数值模拟的过程中采用软硬相间的复合岩体构筑 模型进行模拟。计算模型简化为平面应变问题，装药简化为一个集中装药点，装药量根据隧 道的开挖面积，按照单位长度进行计算。岩石介质条件及爆源条件列于表1、表2。 表1 围岩力学参数