非对称不耦合装药结构对预裂爆破效果的影响.pdf
第34卷 第3期 2017年9月 爆 破 BLASTING Vol. 34 No. 3 Sep. 2017 doi10. 3963/ j. issn. 1001 -487X. 2017. 03. 005 非对称不耦合装药结构对预裂爆破效果的影响* 李新平 1a, 陈萍萍1b, 罗 忆 1a, 黄俊红1b, 焦海棠2, 方延洪3, 詹剑霞3 (1.武汉理工大学a.道路桥梁与结构工程湖北省重点实验室; b.土木工程与建筑学院, 武汉430070;2.武汉地震工程研究院, 武汉430071; 3.中国葛洲坝集团三峡建设工程有限公司, 宜昌443002) 摘 要 为降低爆破开挖对保留区岩体的损伤, 采用数值模拟的方法, 研究不同装药结构对预裂爆破效果 的影响, 发现相比于中心单药卷和偏心双药卷, 偏心单药卷装药结构下的预裂缝成型效果最佳, 且其炮孔壁 上的爆炸荷载在保留区一侧较小, 在开挖区一侧较大, 能在降低爆破开挖对保留区岩体损伤的同时将爆炸能 量更多的作用于开挖区岩体。基于偏心单药卷装药结构, 在综合考虑预裂缝成型效果和保留区岩体损伤程 度的情况下, 通过数值模拟的方法研究了爆破参数对预裂爆破效果的影响, 结果表明 在不耦合系数为1. 68, 炮孔间距为70 cm的条件下, 预裂爆破效果最佳。 关键词 预裂爆破;装药结构;不耦合系数;炮孔间距;数值模拟 中图分类号 TD235 文献标识码 A 文章编号 1001 -487X(2017)03 -0025 -06 Impacts of Asymmetric Decoupling Charge Structure on Presplitting Blasting Effect LI Xin-ping1a,CHEN Ping-ping1b,LUO Yi1a,HUANG Jun-hong1b,JIAO Hai-tang2, FANG Yan-hong3,ZHAN Jian-xia3 (1. a. Hubei Key Laboratory of Roadway Bridge&Structure Engineering; b. School of Civil Engineering and Architecture,Wuhan University of Technology,Wuhan 430070,China;2. Wuhan Institute of Earthquake Engineering,Wuhan 430071,China; 3. China Three Gorges Construction Engineering Corporation Limited of Gezhouba Dam Group,Yichang 443002,China) Abstract To reduce the damage on the reserved rock mass by blasting excavation,the influence of different charging structures on blasting effect was studied through numerical simulation. It is found that the asymmetric char- ging structure with single cartridge formed the best presplitting crack,and caused smaller blasting load on the hole wall at the side of reserved area and bigger blasting load at the side of excavation area,which reduced the damage on reserved rock mass and takes more blasting energy on excavation area as well. On the basis of asymmetric charging structure with single cartridge,the influence on blasting effect of presplitting blasting parameters was also studied through numerical simulation,and the asymmetric charging structure with single cartridge reached the optimal blas- ting effect with decoupling coefficient 1. 68 and hole spacing 70 cm by comparing with presplitting crack forming and the damage degree on reserved rock mass. Key words presplitting blasting;charging structure;decoupling coefficient;hole spacing;numerical simula- tion 万方数据 收稿日期2017 -04 -04 作者简介李新平(1963 -) , 男, 湖北省武汉市, 教授、 博士, 从事爆 破控制技术与爆破效应的研究, (E-mail)xinpingli@ whut. edu. cn。 通讯作者陈萍萍(1993 -) , 女, 湖北省武汉市, 硕士, 从事爆破控制 技术与爆破效应的研究,(E-mail)pingchen@ whut. edu. cn。 基金项目湖北省自然科学基金重点项目(2015CFA136) ; 国家自然 科学基金资助项目(51274157) ; 国家自然科学基金青年 项目(51309183、51609183) 随着水利水电、 采矿、 交通运输工程的规模不断 增大, 工程上大量采用了爆破开挖技术。为了保证 施工质量、 降低事故发生的概率、 同时提高岩石的爆 破成型效果, 依据定向成缝爆破机理的预裂爆破技 术被广泛投入使用[ 1]。预裂爆破在对岩石进行爆 破成缝时, 不可避免地会对保留区岩体产生动力损 伤[ 2,3], 为了降低对保留区岩体的损伤, 在实际预裂 爆破工程中常常采用非对称不耦合装药结构, 即药 卷轴心和炮孔轴心不重合的装药结构, 又叫偏心不 耦合装药结构, 研究不对称装药结构对实际爆破工 程有重要意义[ 4-6]。 国内外学者在研究预裂爆破的非对称不耦合装 药结构发现, 采用非对称不耦合装药结构时, 耦合侧 的应力和应变峰值比非耦合侧高出8倍左右[ 7], 保 留区岩体的裂纹长度远小于开挖区岩体裂纹的长 度[ 4], 不耦合系数与线装药密度、 工程地质条件、 炮 孔孔径以及炮孔间距等参数相关[ 8], 且随着不耦合 系数的增大, 在炮孔壁上的应力波的峰值压力随之 减小, 炮孔壁的损伤程度也会由于峰值压力的降低 而减小[ 5,7,9]。非对称不耦合装药结构在爆破过程 中会在炮孔壁上产生一定的偏心效应, 在进行预裂 爆破装药结构设计时, 应该使较大的爆炸冲击力作 用在开挖区岩体上, 避免其作用在保留区岩体, 降低 对保留区岩体的损伤[ 10]。 除装药结构对预裂爆破效果有重要影响以外, 爆破参数也会对预裂爆破效果产生很大影响, 它主 要包括不耦合系数、 炮孔间距等。预裂爆破的参数 设计, 应该根据工程目的、 施工要求、 爆破对象以及 周围的环境对爆破安全的具体要求等, 综合分析并 慎重选出比较合理的爆破计算参数。预裂爆破参数 设计的基本要求是形成一条贯通的预裂缝, 将保留 区岩体与开挖区岩体充分隔离, 达到保护保留区岩 体, 并使开挖区岩体达到所需破碎程度的目的。 本文就不耦合装药结构对预裂爆破效果的影响 以及爆破参数的优化比选展开研究。 1 预裂爆破及其数值模拟的理论研究 1. 1 预裂爆破的技术设计 预裂爆破是通过形成一条贯通的预裂缝来降低 岩石的超、 欠挖量并且保护保留区岩体的一种爆破 方法。在进行技术设计时, 需要先确定预裂爆破的 爆破参数。如果将由工程爆破经验总结出的爆破参 数或者计算公式用于指导现场爆破施工, 由于不同 的施工现场的工程地质条件会存在一定的差异, 因 而获得的爆破效果不会很理想, 同时施工过程中会 存在一定隐患, 所以需要选择正确的爆破参数用于 爆破施工。通常需要合理的选择和确定线装药密 度、 炮孔间距以及不耦合系数等[ 11]。在重要的预裂 爆破工程中, 一般还应结合爆破工程的具体条件, 安 排小范围的局部试爆, 然后依据试爆结果来选择最 优的爆破参数。 工程爆破参数选定以后, 要进行炮孔的布孔设 计和装药结构设计, 当使用深孔装药结构时需要考 虑装药结构的轴向不耦合系数。总之, 要尽量使药 卷均匀的分布在炮孔中, 避免装药量过于集中。 1. 2 数值模拟参数的选取 (1)材料本构模型的选取 在爆炸冲击荷载作用下, 塑性随动强化本构模 型广泛被采用[ 12], 这是由于炸药爆炸后, 强大的爆 炸冲击力作用在炮孔壁上, 此时岩体内部的应力变 化一般是呈塑性的。之所以岩体内部的应力变化不 呈弹脆性, 是因为在爆炸冲击力作用的瞬间岩体实 际的动态屈服强度会远远高于在试验中得到的静态 屈服强度。 基于国外学者考珀与西蒙兹的研究基础可以建 立起塑性随动强化模型[ 13], 公式如下所示 σy=1 + ε J C 1 ] p (σ0+ βEpεp eff) (1) 式中σy为动态极限的屈服应力;ε为应变率; σ0为初始的屈服应力;C、P为考珀与西蒙兹应变率 参数;β为调节参数,β =0时公式是塑性随动强化模 型,β =1时公式是等向强化模型。 (2) 流固耦合算法 数值模拟采用多物质流固耦合算法, 其基本原 理是在所建模型中炸药和空气单元采用ALE算法, 岩石单元采用拉格朗日算法, 但是各单元之间的处 理方式采用流固耦合算法, 这样在定义烈性炸药材 料单元的时候, 就不会出现畸变, 与此同时还可模拟 出爆炸冲击荷载对炮孔壁周边岩石单元的加载过 程, 因此能在后处理器中直观的看出爆炸效果以及 62爆 破 2017年9月 万方数据 不同位置岩石单元所受应力的变化过程, 得到的效 果也更接近于实际工程爆破效果。 (3) 岩石介质在爆炸过程中的破坏判据 岩石介质在爆炸过程中的破坏判据一般情况下 有四种, 分别是岩石受张拉作用破坏判据、 岩石受剪 切作用破坏判据、 岩石发生滚动破坏判据、 岩石沿结 构面发生滑动破坏判据等。因为岩石的抗拉强度要 远小于其抗压强度, 且预裂爆破过程中预裂炮孔主 要受拉应力作用发生破坏, 故本次数值模拟选用岩 石张拉破坏判据来判断数值模拟中岩石单元是否发 生破坏。具体表现为当数值模型中岩石单元受到的 拉应力超过岩石的抗拉强度或爆炸冲击力形成的压 应力同切向拉应力的矢量方向互相垂直时, 岩石单 元会出现拉裂破坏。 2 数值模型的建立 2. 1 不同装药结构的数值模拟 为了研究不耦合装药结构对爆破效果的影响, 建立中心单药卷、 偏心单药卷和偏心双药卷三种不 耦合装药结构的数值模型来模拟自然条件下岩体受 到爆炸冲击荷载的变化过程。所建模型尺寸为 (X Y)140 cm 200 cm, 炮孔直径为42 mm, 药卷 直径为25 mm( 工程中药卷轴向采用空气间隔装药, 其数值模拟时等效为药卷直径为14 mm的连续装 药结构) , 炮孔间距为70 cm。三种不耦合装药结构 的示意图如图1所示。 图1 三种不耦合装药结构的数值模型 Fig. 1 Numerical models of three decoupling charge structures 2. 2 预裂爆破参数的优化比选 预裂爆破的成缝质量和预裂面的平整度是检验 预裂爆破效果的关键因素[ 14]。爆破参数会对预裂 爆破的效果产生很大影响, 故在对不同不耦合装药 结构进行数值模拟之后, 为了优化爆破参数, 本节以 偏心单药卷为基础研究不耦合系数和炮孔间距对预 裂爆破效果的影响。 (1) 不耦合系数的优化比选 数值模拟通过设置三组不同的装药不耦合系数 来研究不耦合系数对爆破效果的影响, 过程中所使 用的单位制、 算法、 岩石结构单元类型与炮孔直径等 均与上节相同, 炮孔间距为70 cm, 药卷直径分别设 置为18 mm、25 mm与32 mm, 通过计算可得其不耦 合系数分别为 当l1= 18 mm时,k1= 42 18 = 2. 33; 当 l2=25 mm时,k2= 42 25 =1. 68; 当l3=32 mm时,k3= 42 32 =1. 31。 (2) 炮孔间距的优化比选 数值模拟通过设置三组不同的炮孔间距来研究 炮孔间距对爆破效果的影响, 过程中所使用的单位 制、 算法、 岩石结构单元类型与炮孔直径等均与上节 相同,药卷直径取25 mm,炮孔间距分别设置为 50 cm、70 cm和90 cm。 2. 3 炮孔壁荷载的数值模拟 各工况的数值模型沿炮孔中轴线都是左右对称 的, 对其进行分析的时候, 在模型炮孔的左侧沿炮孔 壁从炮孔顶端(保留区)到炮孔底端(开挖区)等间 距选取10个岩石单元来分析爆破过程中炮孔壁荷 载的变化规律, 岩石单元的选取情况如下图2所示。 图2 研究炮孔壁荷载变化规律时选取的 岩石单元示意图 Fig. 2 The selected rock element for the changing rule research of hole wall loading 3 数值模拟结果及分析 3. 1 不同装药结构的数值模拟结果及分析 通过对比中心单药卷、 偏心单药卷和偏心双药 卷三种不耦合装药结构数值模拟的成缝效果可以发 现, 当炮孔直径和药卷直径一定时, 改变装药结构, 爆炸产生的效果会有明显差别。中心单药卷装药结 构可以形成一条贯通的预裂缝, 但是它会严重的损 伤保留区岩体, 无法达到保护保留区岩体的目的; 偏 心双药卷装药结构与中心单药卷装药结构的成缝效 果相似, 虽然沿两炮孔连线形成一条贯通的预裂缝, 但对保留区岩体造成很大的损伤; 偏心单药卷装药 结构不但可以沿炮孔连线形成一条贯通的预裂缝, 而且对保留区岩体的损伤很小, 从而提高了炸药的 72第34卷 第3期 李新平, 陈萍萍, 罗 忆, 等 非对称不耦合装药结构对预裂爆破效果的影响 万方数据 利用率。三种不耦合装药结构的数值模拟成缝效果 如图3所示。 图3 三种不耦合装药结构数值模拟的预裂缝成型图 Fig. 3 Presplitting crack forming of three decoupling charge structures 为了研究不同装药结构对炮孔壁上爆炸荷载的 影响, 分别作出三种装药结构炮孔壁上所选岩石单 元的爆炸荷载峰值变化曲线( 图4) 。从图中可以看 出, 中心单药卷装药结构炮孔壁上的爆炸荷载峰值 从炮孔顶部(保留区)沿炮孔壁到炮孔底部(开挖 区) 基本保持一致, 曲线整体比较平缓; 偏心双药卷 装药结构炮孔壁上的爆炸荷载峰值变化曲线从炮孔 顶部( 保留区)沿炮孔壁到炮孔底部(开挖区)先减 小后增大, 曲线整体起伏较大, 但炮孔顶部(保留 区) 和炮孔底部(开挖区)的爆炸荷载峰值大小相 近; 偏心单药卷装药结构炮孔壁上的爆炸荷载峰值 从炮孔顶部( 保留区) 沿炮孔壁逐渐增大, 在炮孔底 部( 开挖区)达到最大, 曲线整体起伏较大, 开挖区 一侧的爆炸荷载远大于保留区。 图4 不同装药结构炮孔壁上所选单元的 爆炸荷载峰值变化曲线 Fig. 4 Peak explosion loading curves of the selected elements 综合以上分析可以发现, 三种不耦合装药结构 中, 偏心单药卷装药结构不仅可以沿炮孔连线形成 一条贯通的预裂缝将保留区岩体和开挖区岩体隔 离, 而且其在开挖区一侧产生的爆炸荷载远大于保 留区, 既充分粉碎了开挖区岩体又达到了保护保留 区岩体的目的, 故三种不耦合装药结构中, 偏心单药 卷装药结构最佳。 3. 2 预裂爆破参数的优化比选结果及分析 通过上节中对三种不同不耦合装药结构进行的 数值模拟结果对比分析, 已经得出偏心单药卷为最 佳装药结构。本节以偏心单药卷装药结构为研究背 景, 对不同的不耦合系数与炮孔间距进行数值模拟, 研究爆破后的预裂缝成型效果以及爆炸荷载变化规 律, 从而达到优化预裂爆破参数以指导现场工程爆 破施工的目的。 3. 2. 1 不耦合系数的优化比选结果及分析 设置三组具有不同不耦合系数的工况, 令它们 的炮孔直径、 炮孔间距(a =70 cm) 以及其他参数均与 3.1节中相同, 通过改变药卷直径设计了三种不耦合 系数k1=2.33、k2=1. 68、k3=1. 31, 对比数值模拟结 果, 从而选出最优的不耦合系数。见图5、 图6。 图5 三种不耦合系数下的数值模拟预裂缝成型图 Fig. 5 Presplitting crack forming of three decoupling coefficients 通过对三种不耦合系数装药结构的数值模拟结 果进行对比, 可以得到以下结论 (1) 从图5三种不耦合系数下的预裂缝成型效 果可以看出, 不耦合系数k =2. 33时, 只在开挖区岩 体中形成一些裂缝, 岩石破碎程度很低, 炮孔连线上 没有形成贯通的预裂缝;k =1. 31时, 岩石的破碎程 度较高, 炮孔之间虽然形成了贯通的预裂缝, 但是它 对保留区岩体中也形成了很多裂缝, 对保留区岩体 的损伤较为严重;k =1. 68时, 在炮孔连线上形成了 贯通的预裂缝, 开挖区岩石破碎程度较高而保留区 岩体完整度良好, 对保留区岩体损伤很小, 其爆炸能 量主要作用于开挖区岩体, 达到了预裂爆破的根本 目的。 (2) 通过图6三种不耦合系数下炮孔壁上所选 82爆 破 2017年9月 万方数据 单元的爆炸荷载峰值变化曲线可知, 三种不耦合系 数下炮孔壁上的爆炸荷载峰值都是从炮孔顶部(保 留区) 沿炮孔壁逐渐增大, 在炮孔底部(开挖区)达 到最大。三条曲线走向基本一致, 不耦合系数k = 1. 68的爆炸荷载峰值变化曲线始终位于k =2. 33和 k =1. 31的两条曲线之间。 图6 三种不耦合系数下炮孔壁上所选单元的 爆炸荷载峰值变化曲线 Fig. 6 Peak explosion loading curves of the selected elements 综合以上分析可以得出, 三组不耦合系数的模拟 中,k =1.68时不仅可以沿炮孔连线形成一条贯通的 预裂缝将保留区岩体和开挖区岩体隔离, 而且其在充 分粉碎开挖区岩体的同时对保留区岩体损伤较小, 故 三组不耦合系数中,k =1.68时爆破效果最佳。 3. 2. 2 炮孔间距的优化比选结果及分析 设置三组具有不同炮孔间距的工况, 它们的炮 孔直径、 药卷直径(25 mm)以及其他参数均与3. 1 节中相同, 设计三组炮孔间距a =50 cm、a =70 cm、 a =90 cm, 通过对比数值模拟结果, 从而选出最优的 炮孔间距。见图7、 图8。 图7 不同炮孔间距下的数值模拟预裂缝成型图 Fig. 7 Presplitting crack forming of three hole spacings 通过对三种炮孔间距装药结构的数值模拟结果 进行对比研究和分析, 可以得到以下结论 (1) 从图7三种炮孔间距下的预裂缝成型图可 以看出, 在三种炮孔间距下, 炮孔之间均可形成贯通 的预裂缝,但是成缝效果有着很大的差别。a = 50 cm时, 炮孔两侧的预裂缝不在同一平面上, 无法 形成光滑平整的预裂面;a = 90 cm时, 所形成预裂 缝的宽度过窄, 这样在进行开挖区的爆破开挖时, 预 裂缝难以起到有效的隔振作用;a =70 cm时的预裂 缝成型效果相比前两者更好, 预裂缝既在一个平面 上, 又有一定的宽度, 这样不仅可以降低保留区岩体 的损伤, 保证岩体整体的稳定性, 还能够保证预裂面 的光滑平整。 图8 三种炮孔间距下炮孔壁上所选单元 的爆炸荷载峰值变化曲线 Fig. 8 Peak explosion loading curves of the selected elements (2) 从图8三种炮孔间距下炮孔壁上所选单元 的爆炸荷载峰值变化曲线中可以看出, 三种炮孔间 距下炮孔壁上的爆炸荷载峰值都是从炮孔顶部(保 留区) 往炮孔底部(开挖区)沿炮孔壁逐渐增大, 在 炮孔底部达到最大。三种炮孔间距下炮孔壁顶部 ( 保留区)的爆炸荷载相近, 三条曲线走向基本一 致,a =70 cm的爆炸荷载峰值沿着泡孔壁增长中始 终位于a =50 cm和a =90 cm的两条曲线之间。 三组炮孔间距的数值模拟中,a =70 cm时可以 沿炮孔连线形成一条在同一平面上、 且具有一定宽 度的、 贯通的预裂缝将保留区岩体和开挖区岩体隔 离, 而且其在充分粉碎开挖区岩体的同时对保留区 岩体损伤较小。故综合爆破效果和工程经济性来 看, 三组炮孔间距中,a =70 cm效果最佳。 4 结论 利用动力有限元分析软件对不同装药结构和预 裂爆破参数下的预裂缝成型效果及炮孔壁上的爆炸 荷载进行研究, 得到的主要结论如下 (1) 通过对比中心单药卷、 偏心单药卷和偏心 92第34卷 第3期 李新平, 陈萍萍, 罗 忆, 等 非对称不耦合装药结构对预裂爆破效果的影响 万方数据 双药卷三种不耦合装药结构下数值模拟的预裂爆破 效果发现, 偏心单药卷装药结构在保留区岩体一侧 的炮孔壁上产生的爆炸荷载明显低于开挖区, 其对 保留区岩体造成的损伤最小, 且预裂缝成型效果最 佳。( 2) 通过分别建立三组不同不耦合系数和炮孔 间距下的数值模型, 并对比其预裂爆破效果可以看 出, 当不耦合系数k =1. 68、 炮孔间距a =70 cm时, 预裂缝成型效果最佳。当不耦合系数过大时, 难以 形成贯通的预裂缝, 当不耦合系数过小时, 会对保留 区岩体造成较大的损伤; 当炮孔间距过大时, 同样难 以形成贯通的预裂缝, 而当炮孔间距过小时, 不仅不 经济, 而且形成的预裂缝不在一个平面上, 对保留区 岩体造成的损伤也较大, 难以达到预裂爆破的目的。 参考文献(References) [1] 张志呈.定向断裂控制爆破[M].重庆 重庆出版社, 2000. [2] 严 鹏, 邹玉君, 卢文波, 等.基于爆破振动监测的岩 石边坡开挖损伤区预测[J].岩石力学与工程学报, 2016,35(3) 538-548. [2] YAN Peng,ZOU Yu-jun,LU Wen-bo,et al. Predicting the damage zone of rock slopes under blasting excavation based on vibration monitoring[J]. Chinese Journal of Rock Mechanics and Engineering,2016,35(3) 538-548. (in Chinese) [3] 朱 强, 陈 明, 郑炳旭, 等.空气间隔装药预裂爆破 岩体损伤分布特征及控制技术[J].岩石力学与工程 学报,2016,35(S1) 2758-2765. [3] ZHU Qiang,CHEN Ming,ZHENG Bing-xu,et al. Distri- bution and control technology of rock damage induced by air-deck charge presplitting blasting[J]. Chinese Journal of Rock Mechanics and Engineering,2016,35(S1) 2758- 2765.(in Chinese) [4] 管少华, 蒲传金, 肖定军, 等.偏心不耦合装药爆破裂 纹扩展实验研究[J].爆破,2015,32(1) 16-21. [4] GUAN Shao-hua,PU Chuan-jin,XIAO Ding-jun,et al. Experimental study on crack propagation under eccentric decouple charge structure[J]. Blasting,2015,32(1) 16- 21.(in Chinese) [5] 田文高, 蒲传金, 陈晓玲.偏心不耦合装药爆破参数理 论分析[J].爆破,2008,25(2) 19-20. [5] TIAN Wen-gao,PU Chuan-jin,CHEN Xiao-ling. Theoreti- cal analysis on the parameters of eccentric decouple charge blasting[J]. Blasting,2008,25(2) 19-20.(in Chinese) [6] 张建华.偏心不耦合装药爆炸应力场研究及应用[J]. 工业安全与环保,2001,27(8) 20-23. [6] ZHANG Jian-hua. Stress-field research and application of eccentric decouple charge[J]. Industrial Safety and Envi- ronmental Protection,2001,27(8) 20-23.(in Chinese) [7] 张志呈, 史瑾瑾, 蒲传金, 等.偏心不耦合装药对岩石 损伤的试验研究[J].爆破,2006,23(4) 4-8. [7] ZHANG Zhi-cheng,SHI Jing-jing,PU Chuan-jin,et al. Experimental research of rock damage by eccentric decou- ple charge blasting[J]. Blasting,2006,23(4) 4-8.(in Chinese) [8] 罗 勇, 崔晓荣.工程爆破中装药不耦合系数的研究 [J].有色金属 矿山部分,2008,60(4) 39-43. [8] LUO Yong,CUI Xiao-rong. Study on the decoupling char- ging coefficient in engineering blasting[J]. Nonferrous Metals(Mining section) ,2008,60(4) 39-43.(in Chi- nese) [9] 张建华, 王玉杰, 梁 锐.偏心不耦合装药爆炸应力场 的动光弹研究[J].爆破,2001,18(1) 8-12. [9] ZHANG Jian-hua,WANG Yu-jie,LIANG Rui. Study on stress-field in eccentric decouple charge by dynamic pho- to-elastic experiments[J]. Blasting,2001,18(1) 8-12. (in Chinese) [10] 袁光英, 王庆国.不耦合装药爆破效应数值分析[J]. 陕西理工学院学报 自然科学版,2013(5) 28-31. [10] YUAN Guang-ying,WANG Qing-guo. Numerical analysis on the effect of decoupling charge blasting[J]. Journal of Shanxi University of Technology(Natural Science Edi- tion) ,2013(5) 28-31.(in Chinese) [11] HUANG T H,CHANG C S,YANG Z Y. Elastic module for fractured rock mass[J]. Rock Mechanics & Rock Engineering,1995,28(3) 135-144. [12] 罗 忆, 李新平, 董 千, 等.动静力组合作用下深埋 隧洞开挖卸荷诱发围岩动力损伤诱因研究[J].岩石 力学与工程学报,2015,34(S1) 3365-3371. [12] LUO Yi,LI Xin-ping,DONG Qian,et al. Damage cause study for static and dynamic combination of excavation and unloading in deep underground tunnels[J]. Chinese Journal of Rock Mechanics and Engineering,2015, 34(S1) 3365-3371.(in Chinese) [13] JONES N. Structural impact[M]. UKCambridge Uni- versity Press,1989. [14] 璩世杰, 刘际飞.节理角度对预裂爆破成缝效果的影 响研究[J].岩土力学,2015,36(1) 189-204. [14] QU Shi-jie,LIU Ji-fei. Numerical analysis of joint angle effect on cracking with presplit blasting[J]. Rock and Soil Mechanics,2015,36(1) 189-204.(in Chinese) 03爆 破 2017年9月 万方数据