爆破作用下节理岩体应变能变化规律及影响研究.pdf

第 3 5卷第 3期 2 0 1 8年 9月 爆破 BLASTI NG Vo 1 ． 3 5 No ． 3 S e p． 2 01 8 d o i 1 0 ． 3 9 6 3 ／ j ． i s s n ． 1 0 0 1 4 8 7 X． 2 0 1 8 ． 0 3 ． 0 0 1 爆破作 用下节理岩体应变 能变化规律及 影响研 究 术 马小 满 ， 陈 a ,b 严 a ,b 刘 义佳 武汉大学 a ． 水资源与水电工程科学国家重点实验室； b ． 水工岩石力学教育部重点实验室， 武汉 4 3 0 0 7 2 摘要 节理岩体 的爆破破 裂机 制及特征是岩土工程领域 的研 究热 点之 一， 从岩体应 变能 角度进行研 究更 接近岩体破坏本质。建立节理岩体的爆破开挖模型， 采用波动理论推导出爆破引起的节理岩体的应变能密 度和节理张开位移的理论解。理论解表明 节理岩体的应变能与爆破荷载强度及初始地应力水平密切相关。 以白鹤滩水电站工程为实例， 采用理论解计算了爆破荷栽从 1 0～ 2 0 M P a ， 初始地应力从 0～1 0 MP a时的岩 体应变能密度与节理张开位移， 结果表明 爆破荷载作用下岩体应变能密度表现出明显的先积蓄和后释放特 征 ， 爆破荷载越大， 节理岩体积蓄的应变能越多， 释放后导致的岩体节理面的张开位移越大。初始地应力的 大小对岩体应变能和岩体节理面的张开位移的影响呈现相似规律。 关键词 岩体 ； 节理 面； 爆破 ； 应 变能；位移 中图分类号 0 3 8 2 文献标识码 A 文章编号 1 0 0 1 4 8 7 X 2 0 1 8 0 3 0 0 0 1 0 6 S t u d y o n S t r a i n En e r g y Pr o c e s s e s a n d I n flu e n c e o n J o i n t e d Ro c k M a s s u n d e r Bl a s t i n g Ex c a v a t i o n M A X i a o m a n 一 ， C H E N Mi n g 一， Y A N P e n g ， L I U Y i - j i a ， a ． S t a t e K e y L a b o f Wa t e r R e s o u r c e s a n d H y d r o p o w e r E n g i n e e r i n g S c i e n c e ； b． Ke y La b o f Ro c k Me c h i n Hy d r a u l i c S t r u e t Eng Mi n i s t r y o f Ed uc a t i o n， Wu h a n U n i v e r s i t y ， Wu h a n 4 3 0 0 7 2 ， C h i n a Abs t r a c tT h e r o c k f r a c t u r i n g me c h a n i s m a n d c h a r a c t e ri s t i c s u n d e r b l a s t i n g e x c a v a t i o n a r e h o t s p o t s i n g e o t e c h n i e a l e n gi n e e ring fie l d ．I t S c l o s e r t o t h e e s s e nc e t o s t ud y s t r a i n e n e r g y d en s i t y ．I n t h i s s t u dy， t h e t h e o r e t i c a l s t r a i n e l l e r g y o f j o i n t e d r o c k m a s s a n d j o i n t o p e n i n g d i s p l a c e m e n t a r e o b t a i n e d b y s o l v i n g t h e m a t h e m a t i c a l p h y s i c s p r o b l e m 0 f a s i m p l i fi e d me c h a n i c a l m o d e l b a s e d o n t h e e x c a v a t i o n o n j o i n t e d m a s s ． I t i n d i c a t e s t h a t s t r a i n e n e r g y p r o c e s s a n d t h e i n flue n c e ha v e a c l o s e r e l a t i o ns h i p wi t h t he i n t e n s i t y o f b l a s t i n g l o a d a nd i n i t i a l s t r e s s ． W i t h t h e t h e o r e t i c a l s o l ut i o n， r o c k e x c a v a t i o n i n B a i h e t a n h y d r o p o w e r s t a t i o n i s d i s c u s s e d a s a n e x a mp l e w h e r e t h e b l a s t i n g l o a d r a n g e s f r o m 1 0 【o 20 MPa wh i l e t h e i n s i t u s t r e s s fro m 0 t o 1 0 MPa． Th e r e s u l t s s h o w t h a t s t r a i n e n e r g y o f r o c k i s f i r s t a c c u mu l a t e d， t h e n r e l e a s e d u n d e r b l a s t i n g l o a d ． B i g g e r b l a s t i n g l o a d me a n s m o r e i n p u t e n e r gy， t h u s c r e a t e s b i g g e r j o i n t o p e n i n g d i s p l a c e m e n t ． T h e i n s i t u s t r e s s l e v e l h a s s i m i l a r e f f e c t o n s t r a i n e n e r gy a n d j o i n t o p e n i n g d i s p l a c e m e n t w i t h b l a s t i n g l o a d ． Ke y wo r d s r o c k m a s s ； j o i n t s ； b l a s t i n g ； s t r a i n e n e r gy；d i s p l a c e m e n t 收稿 日期 2 0 1 8 0 4 0 9 作者简 介 马小满 1 9 9 5一 ， 女 ， 甘肃西 和人 ， 硕士研究 生， 从事岩石 动力学方 向的研 究工 作 ， E m a i l m a x i a o ma n w h u ． e d u ． cno 通讯作者 陈明 1 9 7 7一 ， 男 ， 湖南末 阳人 ， 博士 、 教 授 ， 主要从事 水利水 电工程施工技术及岩石动力学方面的教 学与研究 ， E - m a i l w h u c h m w h u ． e d u ． c n 。 基金项 目 国家 自然科学项 目面上项 目 5 1 4 7 9 1 4 7、 5 1 2 7 9 1 4 6 我 国8 0 ％ 以上 的水能资源分布在地质 条件复 杂、 岩体节理发育的西南高山峡谷地 区， 新建的大型 水 电工程如二滩 、 锦屏二级水电站等在爆破开挖时 均揭露 出岩体卸荷松动 和大变形 现象 ， 因此 有 必要研究爆破开挖时节理岩体 的破坏机制和特征 。 2 爆破 2 0 1 8年 9月 节理 岩体在爆破荷 载和高地应力的耦合作用 下 ， 响应更 为显著 ， 具 有典型卸荷松动特 征 I4 J ， 损 伤范 围更 大 J 。在对其 爆破破 坏 机制 的研 究 中， A b u o v提出， 岩体爆破过程中， 岩体应力的快速释放 是造成附近岩体损伤的原 因E 6 ] ； T a o M提 出岩体开 挖卸荷过程是由应变能密度 S E D 所控制 ； 黎立 云认为岩体的破坏与岩体内可释放的应变能有 关 ； 从应变能角度入手建立 的破坏机制更接近其 本质 J 。现有研究认为 岩体 的爆破开挖是一种动 态破坏过程 ， 岩体在爆破过程 中存在能量 的积蓄和 释放 ， 会引起岩体变形 的挤压 和回弹 ； 岩体应变 能在爆破开挖时经历动态变化 ， 其释放过程与围岩 损伤有关⋯’ j ， 会造成节理 张开 ， 可 以用来 预测 爆破块度、 损伤和松动范围 。研究节理岩体开挖 过程中的应变能变化 ， 对 明确节理岩体爆破破碎机 理和特征有重要意义。 由于岩体位移和应变能变化过程在实际工程很 难直接测量 ， 本文采用理论分析方法研究爆破荷载 作用下节理岩体的应变能变化过程 ， 分析节理岩体 的松动特征和破坏机理 ， 为优化爆破效果提供理论 依据 。 1 节理岩体 应变能的分析方法 1 ． 1 节理岩体的应变能密度 材料的破坏是由储存在化学键中能量的释放所 驱动， 将承受动荷载的岩体作为研究体系， 从热力学 角度进行分析 ， 外荷载会导致岩体产生变形 ， 转换为 可释放 的岩体弹性应变能和岩体耗散能 W 1 式 中 为外界输入能量 ； 为岩体弹性应变 能； 为耗散能。 材料 的破坏一般从局部失效开始 ， 因此在探究 应变能时， 通常不研究应变能总量 ， 而是将材料看成 许多微小单元组成的连续块体 ， 研究每个单元储存 的能量 ， 即应变能密度 S E D 1 U 2 二 式 中 为应变能密度 ； 为单元应力 ； 为单 元应变 。 假定材料为理想弹性体， 其应变能密度 即弹性 应变能密度 E S E D U e 3 二 式 中 为弹性应变 能密度 ； E为材料 弹性模 量 1 ． 2 节理岩体爆破开挖的理论模型 岩体的S E D与单元的应力和应变有关， 在实际 工程难以直接测量 ， 因此对简化后 的工程模 型进行 理论求解 。将实际工程中常见的岩体 台阶爆破进行 简化 ， 如图 1 所示。 施加 面 图 1 节理岩体台阶爆破不 意图 F i g ．1 B e n c h b l a s t i n g o n j o i n t e d r o c k ma s s 图 1中， 开挖岩体被垂直节理和水平节理正交 切割 ， 分为基岩和开挖岩块两部分 ， 右端为爆破开挖 形成的直立面。假定基岩不可压缩 ， 岩块完全弹性 ， 节理面无粘结强度 ， 仅存在摩擦， 将外荷载等效施加 在开挖轮廓面上 ， 则该开挖模型简化为 图 2 ， 岩体被 垂直节理面 A B和水平节理面 B C所切割 ， 以垂直节 理外法线方向为正 向， 模型长为 Z ， 底部受地面摩擦 力． 厂 ， t 的作用 ， 自由端受均布荷载作用 。 A I 一 一I I I I P t I V I I f x ， f I ⋯ ．．．．一．．． ．一 图 2 简化后的节理岩体爆破开挖模型 F i g ．2 S i mp l i fi e d mo d e l o f b l a s t i n g e x c a v a t i o n o n j o i n t e d r o c k ma s s 1 ． 3 爆破作用下节理岩体的应变能求解 在爆破抛掷之前 ， 岩体 的变形满足位移连续性 4 dx 岩体为理想弹性体时 ， 则其应变能密度 S E D 为弹性应变能密度 E S E D ， 将式 4 代入式 3 即 求得应变能密度理论解。 分析简化的爆破开挖模型， 在不考虑泊松 比影 响的情况下 ， 近似按一维波动进行求解。根据岩块 是否 出现正向位移 ， 将岩块的运动划分为两个阶段 ①在外部压荷载的作用下岩块被压缩， 因基岩不可 压缩 ， 节理面闭合， 此时边界条件为一端 固定一端 自 第 3 5卷第 3期 马小满， 陈明， 严鹏， 等爆破作用下节理岩体应变能变化规律及影响研究 由； ②随外部压荷载减小， 岩块开始出现正向位移， 假定此时时刻为 t ， 此时岩块边界条件变为两端 自 由。两个阶段对应的解答如下 f a 鲁 l 十 八 ， l 0 x ‘P 一 ① { o ， ￡ 0 ； E J P z 5 ② ； 。 0 2 U z 磐 ，-t a t 十 1t ， a t L 8 P 0； E i ，。 i p1 o ， -t 一面 1 d 2 p 一a 2 P 广 - t T o O “讪 一 ，0 0 一 l面d P 1 1 从理论解 中可 以看 出， 通过岩体弹性模量 E、 密 度P、 岩块长度 Z 、 摩擦力_ 厂 ， 和外荷载 P ， 即可 计算岩块 的位移 ， 得到岩体的应变能密度。 6 2 节理岩体应变能变化的算例分析 式 中 为位移 ； E为岩块弹性模量 ； p为岩块密 度； 0为弹性波波速， “ 为①阶段结束时刻的位 移； 为①阶段结束时刻的速度函数； t t 。 根据一维波动理论， 可求得方程①的解为 毒．．[ ’ 。 sin i n d 小n ㈩ 式中 t 言 n d 一 1 ， 3， ⋯ 8 方程②的解为 ～∑ 1 2 ⋯ c o s s i n n ，， ⋯ L ”‘ 丁 s in d 丁 l c 0 s 9 式 中 i s 睾 。s 【 。 s 丁n l“fx 2 1 睾{io 丁c o sn 2 ． 1 计算参数 岩体应变 能 的影 响 因素有 E、 P 、 f 、 L ， ， t 和 P ， 以白鹤滩水电站坝址区柱状节理岩体参数及 爆破开挖设计参数为例进行研究 ， 岩体参数取 E 4 0 G P a ， P 2 7 0 0 k g ／ m ， 节理距开挖面的距离近 似取为柱体长度 ， 即 f _ 4 m。 摩擦力按滑动摩擦力计 ， 有 ， g 1 2 式 中 k为岩石动摩擦 系数 ， 根据 白鹤滩工程地 质报告 1 ， 取 0 ． 6 ； g为重力加速度 ， 取 9 ． 8 1 m ／ s 。 节理岩体爆破开挖时 ， 当爆破荷载大于炮孔连 线处的地应力时， 爆破荷载作用下炮孔 问的岩体才 可能爆破成缝 ， 甚至破碎 、 抛掷 。随着爆破荷载衰减 到和地应力水平相 当， 地应力 开始卸载 引， 假定地 应力与爆破荷载 同时卸除， 因此外荷载 P t 是爆破 荷载和地应力卸荷荷载的耦合荷载 ， 如图 3 。 ． ＼ - ．／一 T ／ 釜 P P ’ 荷载 力 图 3 模型外荷载示意图 F i g ．3 T h e e x t e r n a l l o a d o n t h e s i mp l i f i e d mo d e l 爆破荷载用三角形型式 ， 按 J WL方程计算。 P 。 A 一 B ， 一 南 1 3 式中 为爆炸产生的压力； V为比容； E 为炸 药 比内能 ； A、 B、 R 、 R 、 ∞为 J WL状态 方程参数 ， 对 于乳化炸药 ， 取值如表 1 。 根据 白鹤 滩 现 场试 验 ， 地 应 力 水 平 取 o r 5 MP a ， 预裂 孔 深 8 ． 5 m， 孔 距 0 ． 6 m， 炮 孔 孔 径 7 6 m m， 线装药密度 2 8 0 g ／ m， 单孔装药 3 ． 5 k g ， 将荷 4 爆破 2 0 1 8年 9月 载等效为面荷载， 计算得爆破荷载峰值为 P 1 3 ． 2 MP a 。 据杨建华对爆破历程 的研究成果 ， 计 算中爆破荷载上升时间 t ， 取为 1 m s ， 持续时间 t 取 为 6 ms 。 表 1 乳化炸药基本参数及 J WL状态方程参数 Ta b l e 1 Ba s i c P a r a me t e r s o f e mu l s i o n e x p l o s i v e a n d t h e J W L s t a t e e q u a t i o n 地应力按线性卸荷进行计算 ， 根据裂纹扩展速 度及爆轰波的传播速度进行估算 引， 近似认为卸荷 速率与爆破荷载相同。 2 ． 2 应变能变化过程 利用上述参数及式 3 、 4 、 7 、 9 可求解出 岩块上任意一点 的应 变能密度。以位于 自由面处 ， 、 3 1 ／ 4处、 中点 1 ／ 2 、 f ／ 4处和节理 处 0 的点为例， 绘出其应变能密度变化过程， 如 图 4所示。图4是模型不同点处 的应变能密度变化 过程。从应变能密度曲线的变化过程来看 ， 岩块整体 有 3 1 5 J ／ m 。 的初始应变能密度， 为初始地应力储存 在岩体中的能量。随后从爆破临空面 z 到节理 面 0 处依次开始上升， 达到峰值 ， 约 1 2 0 0 0 J ／ m 以上， 远超初始能量密度， 随后迅速下降， 逐渐衰减， 在约 1 5 ms 时刻应变能密度全部衰减到 0 。 0 2 4 6 8 1 0 1 2 1 4 1 6 1 8 20 时间 t ／ ms 图 4 模型不同点处应变能密度变化过程 F i g ．4 S E D C H I V E S o f d i f f e r e n t p o i n t s 对 比不同点的情况， 从右 f 到左 0 处 依次开始积蓄 ， 说明岩块从开挖边界到节理面逐渐 被压缩 ； 从右到左 S E D曲线峰值 逐渐增加 ， 说 明距 离开挖面距离越远 ， 积蓄的应变能越多。 z 处 的 点在 t 1 m s 时 S E D曲线达到峰值， t 6 In s 时释放 完毕 ， 与爆破荷载变化过程完全一致 ， 符合开挖边界 的特点。 0处点在约 1 m s时刻开始积蓄 ， 符合 应力波传递规律 ； 应变能密度率先衰减为 0， 符合边 界点特点 。应变能密度峰值在 t 2 ． 5 m s时达到 ， 此时整个岩块压缩达到了最大。 从机理讲 ， 随爆破荷载施加 ， 能量密度先增加后 减小 ， 表示岩体 中的能量先积蓄， 后释放。积蓄的能 量 随岩块的运动慢慢消耗 ， 最后基本达到零 ， 说明在 爆破破岩过程中， 初始应变能得到完全释放 。 不 同点处 的应变能密度变化规律基本是一致 的， 因此后文分析 中用节理面处的点 0 进行分 析 。 2 ． 3 不同爆炸荷载作用时的应变能变化 不考虑岩石参数的变化 ， 则影响节理岩体应变 能 的主要是爆破荷载和初始地应力。为探究爆破荷 载的影响， 峰值 P 取 1 0 M P a 、 1 2 ． 5 M P a 、 1 5 MP a 、 1 7 ． 5 MP a 和 2 0 MP a ， 计算在地应力 5 M P a下的 应变能变化情况如图 5所示 。 2 40 0 0 2 00 0 0 1 60 0 0 1 20 0 0 80 0 0 40 0 0 O 2 4 6 8 时间 t ／ ms 图 5 不同爆破荷载的应变能密度变化过程 Fi g ．5 SED c u r v e s u n d e r d i f f e r e n t b l a s t i n g l o a ds 从 图 5可以看 出， 不 同爆破荷载下应变能密度 的变化过程均有积蓄和释放规律 。一方面 ， 不 同爆 破荷载时 S E D曲线 到达峰值和降低 为 0的时刻不 变， 曲线形态与荷载强度无关 ， 说 明能量变化路径与 外荷载强度关系不大。另一方面 ， 应变能密度 的数 值与爆破荷载正相关 ， 爆破荷载越大 ， 应变能密度越 大， 岩体中积蓄的能量越多。 2 ． 4不同地应力水平时的应变能变化 岩体中赋存的地应力在爆破荷载的作用下会释 放， 地应力大小 、 释放速率 、 释放路径均会对岩体 的 ∞ ∞ ∞ ∞ ∞ ∞ ∞ ∞ 加 ∞ ∞ ∞ ∞ 加 糙涩 第 3 5卷第 3期 马小满 ， 陈明， 严鹏 ， 等爆破作用下节理岩体应变能变化规律及影响研究 5 破坏产生影 响。为研究不 同地应力水平 的影响 ， 保 持爆破荷载为 P 1 3 ． 2 MP a 不变 ， 计算地应力水平 为 0～1 0 MP a的应变能密度 图 6对 比不 同地应力水平下的 S E D曲线 ， 可以 看 出曲线形态与地应力水平无关 ， 说 明应变能的积 蓄与释放路径不受地应力水平影响。地应力水平则 对应变能密度数值影响较大 ， 地应力水平越高 ， 岩体 的初始应变能密度越大、 S E D曲线峰值也越大 ， 说明 地应力的存在会增加岩体中积蓄的应变能。在实际 爆破工程中， 这表示地应力水平越高， 岩体越容易达 到能量极限而发生破坏 。 1 2 00 0 舞 8 0 0 0 4 00 0 O 2 4 6 时 间t ／ ms 图 6 不同地应力水平下的应变能密度变化过程 Fi g ．6 S ED c u r v e s u nd e r di f f e r en t i n s i t u s t r e s s 3 节理岩体应变能变化 的影响 3 ． 1 节理岩体的位移变化过程 岩体应变能是岩体 中存储 的能量 ， 其对岩体 造 成 的直观影 响是岩体宏观损伤和松动位移 。根据式 7 及式 9 ， 可以计算岩体位移的变化过程 。 根据 3 ． 1节的计算参数 ， 可以得到 图 7所示 岩 体不同部位 的位移过程。图 7表示节理岩体 的不 同 部位位移变化过程 ， 简称 “ 位移图” 。从 中可 以看 出 岩块上不同点处的位移变化规律基 本一致 随爆破 荷载施加 ， 岩体开始有负向位移并不断增大， 达到某 一 值后 回弹向正向运动 ， 正 向位移不断增大并逐渐 达到稳定 ， 岩块稳定位移值为 2 5 mm左右。由图 2 的计算模型可知 ， 该正向位移表示节理的张开位移 ， 这是爆破作用对节理岩体造成的直观影响。 从局部放大图 7 a 中可 以看 出， 对岩体上每一 点有初始位移 ， 对 Z 处点初始位移约 0 ． 5 mm， 是 初始地应力作用下的压缩变形 。 f 处的点首先发 生负 向位移 ， 随后是 1 ／ 2和 0处 的点 ； 在 t 2 ． 5 ms 时岩块负向位移达到最大 ， 与 图 4应变能变 化过程所体现的特征一致。岩块不同点的稳定位移 相 同， 说明岩块初始地应力造成 的初始变形 随地应 力的释放发生了弹性回复。 鲁 g 毡 剐 她 图 7 模型不 l司点处位移时程 曲线 Fi g ．7 Di s p l a ce me nt c ur v e s o f d i f f e r e n t p o i n t s 总结 图 4 、 图 7， 可以看 出岩块位 移与应变能的 变化具有一致性 ， 能量 的积蓄阶段会挤压岩体 ， 岩体 主要产生负向位移； 能量的释放 阶段 引起位移的回 复 ， 而岩块整体位移的弹性 回复 ， 对应岩块初始应变 能的完全释放。爆破作用下 的节理岩体表现出明显 的“ 挤压蓄能” 和 “ 卸荷 回弹” 特征。 3 ． 2 不同爆炸荷载时的张开位移 对应 图 5的爆破荷载及地应力水平 ， 计算得岩 体节理张开位移如 图 8所示 。 吕 g 尉 社 图 8 不 同爆破荷载 F的节理 张开位移 F i g ．8 J o i n t o p e n i n g s u n d e r d i f f e r e n t b l a s t i n g l o a d s 从图 8可 以看出， 在地应力一定 时节理的张开 位移与爆破荷载峰值呈正相关， 爆破荷载越大， 节理 张开位移越大 。随爆破荷载 的增大 ， 这一影响更加 明显 。 3 ． 3 不同地应力水平时的张开位移 与图 7对应 ， 爆破荷载为 2 0 MP a时可得到不 同 地应力水平下 的节理张开位移如图9所示 。 从图 9可 以看出 ， 随地应力水平的增加 ， 节理张 开位移增大。究其原 因， 是 图 6所揭示的岩体积 蓄 的应变能与地应力水平关系 的体现， 高应变能造成 了大节理张开位移 。 6 爆破 2 0 1 8年 9月 吕 目 剐 牡 地应 力 ／ MP a 图 9 不同地应力水平下 的节理 张开位移 F i g ．9 J o i n t o p e n i n g s u n d e r d i f f e r e n t i n s i t u s t r e s s 4 结论 1 通过节理岩体爆破开挖理论模型推导 了爆 破荷载作用下节理岩体应变能变化 的理论解 ， 计算 结果表明节理岩体在爆破开挖过程中应变能密度表 现出明显的能量积蓄和释放特征， 引起岩体位移的 先挤压后回弹 ， 并导致岩体节理面张开。 2 不同爆破荷 载下 的应变 能变化规律表 明， 爆破荷载越大 ， 节理岩体积蓄的能量越高 ， 造成的节 理张开越大 ， 爆破荷载是节理张开的主导因素 。 3 初始地应力对岩体应变能 的影 响明显 ， 地 应力水平越高 ， 岩体应变能越大 ， 释放后造成的节理 张开越大。初始地应力的释放是造成岩块整体变形 的弹性 回复的原因。 必须说 明的是 ， 理论分析是基于部分假定条件 及一维波动理论进行的， 且未考虑岩体参数 、 地应力 卸载速率等因素的影 响， 实际的岩体爆破过程非常 复杂 ， 相关的内容仍有待进一步研究 。 [ 2 ] [ 2 ] 参考文献 R e f e r e n c e s 李正 刚． 二滩水电站地下厂房 系统洞室 围岩变形研 究 [ J ] ． 四川水力发电， 2 0 0 4 1 4 3 4 7 ． L I Z h e n g g a ng ． Er t an h y d r o p o we r s t a t i o n u n d e r g r o u n d p o we r - h o u s e s u r r o u n d i n g r o c k m a s s d e f o r m a t i o n r e s e a r c h『 J ] ． S i c h u a n Wa t e r P o w e r ， 2 O 0 4 1 4 3 4 7 ． i n C h i n e s e 江权， 侯靖 ， 冯夏庭 ， 等． 锦屏二级水 电站地下厂 房围岩局部不稳定问题的实时动态反馈分析与工程 调控研究[ J ] ． 岩石力学与工程学报 ， 2 0 0 8 9 1 8 9 9 一 l 9 07． J I A N G Q u a n ， H O U J i n g ， F E N G X i a ri n g ， e t a 1 ． D y n a m i c f e e d b a c k a n a l y s i s a n d e n g i n e e r i n g c o n t r o l o f s ur r o u nd i n g r o c k l o c a l i n s t a b i l i t y i n u n d e r g r o u n d p o w e r h o u s e o f J i n p i n g i i h y d r o p o w e r s t a t i o n [ J ] ． C h i n e s e J o u r n a l o f R e c k [ 3 ] [ 3 ] [ 4 ] [ 4 ] [ 5 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] [ 8 ] [ 9 ] [ 9 ] [ 1 0 ] [ 1 0 ] Me c h a n i c s a n d E n g i n e e ri n g ， 2 0 0 8 9 1 8 9 9 1 9 0 7 ． i n C h i n e s e 刘际飞， 璩世杰． 节理走向角度对爆炸应力波传播影 响的试验研究[ J ] ． 爆破 ， 2 0 1 4 ， 3 1 2 5 7 6 1 ， 6 6 ． L I U J i f e i ， Q U S h i - j i e ． E x p e ri m e n t a l s t u d y o f i n f l u e n c e o f j o i n t a n g l e o f s t r i k e o n e x p l o s i v e s t r e s s w a v e p r o p a g a t i o n [ J ] ． B l a s t i n g ， 2 0 1 4 ， 3 1 2 5 7 6 1 ， 6 6 ． i n C h i n e s e 张吉宏， 刘红岩． 节理性质对岩体爆破破坏模式影响 的数值分析[ J ] ． 爆破 ， 2 0 1 3 ， 3 0 2 7 9 8 3 ． Z H AN G J i h o n g ， L I U H o n g y a n ． N u me ri c a l a n a l y s i s o f e f f e c t 0 f j o i n t s p r o p e r t y o n r o c k m a s s b l a s t i n g f a i l u r e n x d e l [ J ] ． B l a s t i n g ， 2 0 1 3 ， 3 0 2 7 9 8 3 ． i n C h i n e s e 倪绍虎， 何世海 ， 陈益民， 等． 柱状节理玄武岩的破坏 模式、 破坏机制及工程对策[ J ] ． 岩石力学与工程学 报 ， 2 0 1 6 ， 3 5 S 1 3 0 6 4 3 0 7 5 ． NI S h a o h u， HE S h i h a l ， C HEN Yi Mi n g ， e t a 1 ． T h e f a i l u r e mo d e s， f a i l ur e me c h a n i s ms a nd c o u n t e r me a s u r e s o f c o l n m n a r i o i n t e d b a s a l t r o c k m a s s [ J ] ． C h i n e s e J o u r n a l o f R o c k Me c h a n i c s a n d E n g i n e e r i n g ， 2 0 1 6 ， 3 5 S I 3 0 6 4 3 0 7 5 ． i n C h i n e s e ABUOV M G， AI TAL I EV S M ， ERMEKOV T M ， e t a 1 ． S t u d i e s o f t h e e f f e c t o f dy n a mi c p r o c e s s e s d u rin g e x pl o s i v e b r e a k o u t u p o n t h e r o o f o f mi n i n g e x c a v a t i o n s [ J ] ． J o u r n a l o f Mi n i n g S c i e n c e ， 1 9 8 8， 2 4 6 5 8 1 - 5 9 0 ． TAO M ， LI X， L I D． Roc k f a i l u r e i nd u c e d by d y n a mi c H H l o a d i n g u n d e r 3 D s t r e s s s t a t e [ J ] ． T h e o r e t i c al＆A p p l i e d F r a c t u r e Me c h a n i c s ， 2 0 1 3 ， 6 5 3 4 7 - 5 4 ． i n C h i n e s e 黎立云， 谢和平， 鞠杨， 等． 岩石可释放应变能及耗 散能的实验研究[ J ] ． 工程力学 ， 2 0 1 1 3 3 5 4 0 ． I J l L i y u n， X I E He p i n， J U Ya n g ， e t a 1 ． E x p e r i me n t a l i n v e s t i g a t i o n s o f r e l e a s a b l e e ne r g y a n d d i s s i pa t i v e e ne r g y w i t h i n r o c k[ J ] ． E n g i n e e r i n g M e c h a n i c s ， 2 0 1 1 3 3 5 4 0 ． i n C h i n e s e 郭建强 ， 刘 新荣 ， 王 军保 ， 等 ． 基 于 弹性应 变能 的岩 石 强度准则[ J ] ． 岩土力