精确延时对台阶爆破岩石破碎块度影响的数值模拟研究 .pdf

ｄｏｉ１０. ３９６９/ ｊ. ｉｓｓｎ. １００１-８３５２. ２０１６. ０３. ００３ 精确延时对台阶爆破岩石破碎块度 影响的数值模拟研究 ❋ 李顺波① 杨仁树① 杨 军② ①中国矿业大学北京力学与建筑工程学院北京，１０００８３ ② 北京理工大学爆炸科学与技术国家重点实验室北京，１０００８１ [摘 要] 为了研究精确延时间隔对岩石破碎的影响,利用有限元软件就不同延时间隔条件对岩石破碎的影响进 行了数值模拟研究 研究结果表明合理设定延时间隔利于台阶顶部岩石破碎块度控制,同时能够改善岩石破碎 块度分布情况 通过爆生气体和应力波关系的分析,计算得出的延时间隔和数值模拟结果较为接近；在本次数值 模拟中,延时间隔 ６ ｍｓ、抵抗线 ３ ｍ,即延时间隔和抵抗线关系为 ２ ｍｓ/ ｍ 情况下有利于岩石破碎 [关键词] 岩石破碎；精确延时；数值模拟 [分类号] ＴＤ８５３ 引言 随着数码电子雷管的诞生和使用,毫秒延时间 隔对爆破效果的影响越来越受到人们的重视,特别 是在爆破振动控制和岩石破碎这两个方面 数码电 子雷管具有延时精度高、延时间隔在 １ １ ０００ ｍｓ 范围内随意设定的优点,使得通过改变毫秒延时间 隔来控制岩石破碎块度成为了一种可能 国内外很多学者在这一方面进行了大量的尝 试,Ｓｔａｇｇ 等[１]通过试验得到了毫秒延时间隔 ３. ３ ｍｓ/ ｍ 时利于岩石的破碎；Ｏｔｔｅｒｎｅｓｓ 等[２]通过大量 的实践得出毫秒延时间隔 ３. ３ １３. ０ ｍｓ/ ｍ 时可以 使得大块率降低 １２％ ２０ ％；张乐等[３]通过对隆 芯 １ 号数码电子雷管在德兴铜矿工业试验中的爆破 振动监测,发现应用数码电子雷管可以改善矿岩的 破碎度,提高综合经济效益 Ｐａｌｅｙ 等[４]在美国 Ｒｅｄ Ｄｏｇ Ｍｉｎｅ 中,采用精确 延时电子雷管将延时间隔 ２５ ｍｓ 缩短到１７ ｍｓ；在普 通雷管的情况下,将岩石大块率增加 ２０ ％；若采用 电子雷管,大块率则降低 ３０ ％,这样能很好地改善 岩石破碎效果 Ｋｏｅｎｉｇ[５]利用精确延时电子雷管, 将块度小于 ５００ ｍｍ 的岩石由原来的 ５７％ 增加到 ７９％,同时块度大于 １ ０００ ｍｍ 的岩石由 ２７％ 减少 至 ５％ Ｌｅｗｉｓ[６]、ＭｃＫｉｎｓｔｒｙ[７]等通过研究发现,利 用精确延时电子雷管可以取得良好的破碎效果 Ｂａｔｔｅｎ 等[８]在 Ｎｅｗｌａｎｄｓ 煤矿中利用电子雷管改善 剥离岩体质量,提高了生产效率,取得了良好的经济 效果 Ｐｅｔｒｏｐｏｕｌｏｓ 等[９]通过室内的模型试验,探讨 了不同延时间隔对岩石破碎的影响 Ｓｃｈｉｌｌ[１０]通过 ＬＳ-ＤＹＮＡ 软件模拟了不同毫秒延时间隔情况下岩 石破碎形态 以上分析可以看出,精确延时间隔对岩石破碎 效果的改善已经受到越来越多学者的关注,同时通 过相关试验和现场实践,进一步验证了精确延时在 岩石块度改善方面具有优势 但是,目前对于毫秒 延时对岩石破碎影响的机理和毫秒延时间隔设定都 是不明确的 因此,本文从数值模拟角度出发,进一 步直观地揭示毫秒延时间隔和岩石破碎之间的本质 联系 １ 理论分析 炸药在岩体爆炸时产生冲击波,后续产生应力 波并向远处传播,其在时间和空间上的发展形态如 图 １ 中所示[１１] 其中,“ ＋ ”表示压缩波；“ － ”表示 拉伸波；ΛＷ表示波长；τＷ表示波传播 ΛＷ所需的时 间；σ 为应力；ｔ 为时间；Ｘ 为距离 采用一维拉格朗日描述应力波在时间和空间上 的关系,如图 ２ 所示 其中,ＳＥ为 Ｓ 波横波拉伸 波；ＳＦ为 Ｓ 波压缩波；ＰＥ为 Ｐ 波纵波拉伸波；ＰＦ 为 Ｐ 波压缩波 由于 Ｐ 波和 Ｓ 波传播速度的不同, 因此,在同一距离处二者已经发生了分离 从图 ２ 中看出了 Ｐ 波和 Ｓ 波相互作用的区域的形成 .１１.２０１６ 年 ６ 月 精确延时对台阶爆破岩石破碎块度影响的数值模拟研究 李顺波,等 ❋ 收稿日期２０１５-１０-２７ 作者简介李顺波１９８５ － ,男,博士,中国矿业大学北京博士后流动站从事岩石爆破理论方面的工作 Ｅ-ｍａｉｌｌｓｈｕｎｂｏ＠１２６. ｃｏｍ ａ波长 ｂ传播时间 图 １ 应力波在空间和时间上的形态示意图 Ｆｉｇ. １ Ｓｈａｐｅ ｏｆ ｓｔｒｅｓｓ ｗａｖｅ ｉｎ ｓｐａｃｅ ａｎｄ ｔｉｍｅ 图 ２ Ｐ 波与 Ｓ 波传播形态拉格朗日示意 Ｆｉｇ. ２ Ｌａｇｒａｎｇｅｓ ｃｈｅｍａｔｉｃ ｄｉａｇｒａｍ ｆｏｒ ｐｒｏｐａｇａｔｉｏｎ ｆｏｒｍ ｏｆ Ｐ ｗａｖｅ ａｎｄ Ｓ ｗａｖｅ 进一步考察应力波之间的相互作用,两孔产生 的应力波的一维拉格朗日表示方法见图 ３ １＃、２＃为 炮孔 １＃和炮孔 ２＃ 图３ 中可以看出 Ｓ 波和 Ｐ 波在拉伸波和压缩波 之间形成的不同作用区域 若二者之间存在延时间 隔,孔 ２＃产生的压缩波和拉伸波将沿坐标轴纵轴向 上移动,形成新的作用区域范围,如图 ３ｂ中所示 在 ｔ∗时刻 Ｘ∗处的 Ｏ∗点位于图 ３ａ中,受到孔 １＃ 产生的 Ｓ 波压缩波、孔 ２＃产生的 Ｐ 波压缩波的作 用；同样,在图 ３ｂ中,Ｏ∗点受到孔 １＃产生的 Ｓ 波 压缩波和孔 ２＃产生的 Ｓ 波拉伸波的作用 在同时 起爆情况下,Ｏ∗点受到压缩波的作用；在具有延时 间隔情况下,受到压缩波和拉伸波的作用,因此延时 间隔的存在,能够改变某一点的受力状态,从而能够 改变岩石的破碎效果 目前岩石爆破研究中,倾向于用爆生气体和应 力波相互作用的理论去解释岩石破碎机理 因此, 考虑到岩石裂纹和应力波的相互作用,同时设立图 ４ 中的坐标系 根据图 ４ 中的几何关系,可以得到 爆生气体到达 Ａ 点所需要的时间 ｔｇ＝ １ ＣｇｃｏｓθＸ ＝ ２Ｂ２＋ Ｌ２ ４ ＣｇＬ Ｘ１ ａ同时起爆 ｂ间隔时间为 Δｔ 图 ３ 两孔相互作用的一维拉格朗日示意图 Ｆｉｇ. ３ Ｏｎｅ ｄｉｍｅｎｓｉｏｎａｌ Ｌａｇｒａｎｇｅ ｓｃｈｅｍａｔｉｃ ｄｉａｇｒａｍ ｏｆ ｔｈｅ ｔｗｏ-ｈｏｌｅ ｉｎｔｅｒａｃｔｉｏｎ 图 ４ 爆生气体和应力波相遇示意图 Ｆｉｇ. ４ Ｓｃｈｅｍｅ ｄｉａｇｒａｍ ｒｅｐｒｅｓｅｎｔｉｎｇ ｔｈｅ ｉｎｔｅｒａｃｔｉｏｎ ｂｅｔｗｅｅｎ ｇａｓ ｆｌｏｗ ａｎｄ ｓｔｒｅｓｓ ｗａｖｅ 其中ｔｇ为爆生气体到达 Ａ 点所需的时间；Ｃｇ爆生气 体传播速度；Ｂ 为抵抗线；Ｌ 为孔间距 孔 ２＃产生的应力波到达 Ａ 点所需要的时间为 ｔｐ＝ ４Ｘ２Ｂ２ Ｌ２ ＋ Ｌ － Ｘ２ Ｃｐ ＋ Δｔ２ 其中Ｃｐ为纵波波速；Δｔ 为延时间隔 .２１. 爆 破 器 材 Ｅｘｐｌｏｓｉｖｅ Ｍａｔｅｒｉａｌｓ 第 ４５ 卷第 ３ 期 假设孔 １＃产生的爆生气体和孔 ２＃产生的应力 波在 Ａ 点相遇,在 Ｘ ＝ Ｌ/４ 时[１２]取的极值,此时 Δｔ 可以表示为 Δｔ ＝ Ｂ２＋ Ｌ２ ４ ２Ｃｇ － ９Ｌ２ １６ ＋ Ｂ２ ４ Ｃｐ ３ ２ 数值模拟 ２. １ 数值模拟模型 数码电子雷管在实现精确延时控制的同时,也 实现了逐孔起爆 因此,数值模拟中主要探讨相邻 两个炮孔由于延时间隔的改变造成的影响 由于受 到计算机计算能力限制,建立如图 ５ 所示的小型数 值模拟计算模型 模型高度为 １５. ０ ｍ,炮孔直径为 ２００ ｍｍ,台阶顶部宽度为 １０. ０ ｍ,台阶底部宽度为 １７. ５ ｍ,台阶高度为１０. ０ ｍ,台阶长１８. ０ ｍ,装药高 度 ７. ０ ｍ,填塞长度为 ３. ０ ｍ,两孔间距为 ６. ０ ｍ,抵 图 ５ 台阶模型图单位ｍ Ｆｉｇ. ５ Ｂｅｎｃｈ ｍｏｄｅｌ ｕｎｉｔ ｍ 抗线为 ３. ０ ｍ,台阶顶部和前部为自由边界,其他面 施加无反射边界 设定延时间隔分别为 ０、２、４、６、８ ｍｓ ５ 种情况 ２. ２ 数值模拟材料模型和参数 在 ＬＳ-ＤＹＮＡ 中,对炸药的描述一般采用 ＪＷＬ 状态方程 由于炸药种类、密度等参数不同,其对应 的 ＪＷＬ 状态方程的参数选取也有很大差别,爆轰压 力的状态方程 ｐ ＝ Ａ １ － ω Ｒ１Ｖ ｅ － Ｒ１Ｖ ＋ Ｂ １ － ω Ｒ２Ｖ ｅ － Ｒ２Ｖ ＋ ωＥ Ｖ ４ 其中ｐ 是爆轰压力；Ｅ 是炸药内能；Ｖ 是当前相对体 积；Ａ、Ｂ、Ｒ１、Ｒ２、ω 是 ＪＷＬ 状态方程参数 ２＃岩石乳 化炸药的 ＪＷＬ 参数,如表 １ 所示 岩石和填塞材料用 ＲＨＴ 模型进行描述,其强度 模型是通过引入 ３ 个极限面来实现初始的弹性屈 服面、失效面以及残余强度面 主要材料参数见表 ２ 其中Ｇ 为 Ｇｒｕｎｅｉｓｅｎ 常 数；Ｂ０、Ｂ１、Ｔ１、Ｔ２为状态方程参数；Ａ１、Ａ２、Ａ３为 Ｈａｇｏｎｉｏｔ 方程参数；Ａ、Ｎ 为屈服面参数；Ｑ０、Ｂ 为 Ｌｏｄｅ 角参数；Ｇｃ为压缩屈服面参数；ＧＴ为拉伸屈服 面参数；ＸＩ为剪切模量折减系数；Ｄ１、Ｄ２为损伤因 子；ＡＦ、ＮＦ为剩余面参数 ３ 个曲面解释了沿着不同的子午线强度的降 低,同时也解释了应变率效应 失效面,也就是岩石 最终强度,是由不同的材料参数组成,其中包括岩石 的压缩、拉伸、剪切强度 初始屈服面是有的用户所 表 １ 炸药状态参数 Ｔａｂ. １ Ｓｔａｔｕｓ ｐａｒａｍｅｔｅｒｓ ｏｆ ｔｈｅ ｅｘｐｌｏｓｉｖｅ 密度/ ｇ.ｍｍ －３ 爆速/ ｍｍ.ｍｓ －１ Ａ/ ＭＰａＢ/ ＭＰａωＲ１Ｒ２Ｅ/ Ｊ.ｍ －３ １ １０ －３ ３ ２００１ １０３３. １８ １０４０. ０４７ ９８３. ５４４. ２６８ １０９ 表 ２ 岩石材料参数 Ｔａｂ. ２ Ｐａｒａｍｅｔｅｒｓ ｏｆ ｒｏｃｋｓ 密度/ ｇ.ｍｍ －３ 弹性模量/ ＭＰａ侵蚀塑性应变Ｂ０Ｂ１Ｔ１/ ＭＰａＡＮ ２. ４１４ １０ －３ １６. ７ １０３２. ０１. ２２１. ２２３. ５２７ １０４１. ６０. ６１ 抗压强度/ ＭＰａ相对剪切强度相对拉伸强度Ｑ０ＢＴ２/ ＭＰａ参考压缩应变率/ ％参考拉伸应变率/ ％ ３５１. ００. ３０. ６８０ ５０. ０１０ ５０３ １０ －８ ３ １０ －９ 破坏压缩 应变率/ ％ 破坏拉伸 应变率/ ％ 压缩应变率 指数 拉伸应变率 指数 拉伸体积塑性 应变 ＧＣＧＴＸＩ ３ １０２２３ １０２２０. ０３２０. ０３６０. ００１０. ５３０. ７０. ５ Ｄ１Ｄ２最小损伤剩余应变ＡＦＮＦＧＡ１/ ＭＰａＡ２/ ＭＰａ ０. ０４１. ００. ０１１. ６０. ６１０３. ５２７ １０４３. ９５８ １０４ Ａ３/ ＭＰａ压溃压力/ ＭＰａ压实压力/ ＭＰａ多孔指数初始多孔性 ９. ０４ １０３２３. ３６ １０３３１. １８８ ４ .３１.２０１６ 年 ６ 月 精确延时对台阶爆破岩石破碎块度影响的数值模拟研究 李顺波,等 输入的一部分数据沿着拉伸和压缩子午线的失效 面加上另一段曲线使整个曲线在孔隙压溃压力 值处封闭组成 典型的加载状况,见图 ６ 图 ６ ＲＨＴ 模型的极限表面和加载情况 Ｆｉｇ. ６ Ｅｘｔｒｅｍｅ ｓｕｒｆａｃｅ ａｎｄ ｌｏａｄ ｏｆ ｔｈｅ ＲＨＴ ｍｏｄｅｌ 模型在达到初始屈服面之前为弹性,超过后认 为有塑性应变出现 通过初始屈服面和失效面的差 值的方法,把结合了岩石硬化特性的塑性应变用于 形成有效屈服面 同样,当应力达到失效面时,加入 参数化的损伤模型来控制损伤演化,这个损伤演化 也是由塑性应变驱动的 通过失效面和剩余强度面 的差值,损伤演化继而表现了失效应力后的极限面 对于完全损伤材料,没有子午线或者应变率关系,剪 切强度只是在限制条件即正压力的情况下下才 存在 ２. ３ 数值模拟结果 图７中的整个台阶的破碎形态图中,在不同延 时间隔情况下,炮孔周围的破碎形态呈现不同的变 化范围,特别是在右侧炮孔的根部和顶部的变化形 态显著地不同 在两孔中间的台阶坡面位置上,图 ７ｄ中形成了两条接近连接的破碎裂纹,有助于岩 石破碎形态的分割 其他几种延时间隔情况下,只 有在上部形成一条破碎裂纹,并形成多个“孤岛区 域”,这就容易造成岩石块度的不均匀,不利于整个 岩石块度分布形态 右侧炮孔的底部右下角形 成的延伸裂纹区域,在２ ｍｓ和６ ｍｓ情况下较长；同时 在 ８ ｍｓ 情况下,右侧炮孔底部形成一个对称区域, 可以看出不同的延时间隔对于整个炮孔区域的发展 具有重要的影响 图 ８ 给出了两炮孔连线中点处切面的损伤变化 情况,将损伤值大于 ０. ７ 的区域进行删除后形成图 ８ 中的区域范围 从图 ８ 中可以看出,随着延时间隔的增大,在台 阶上部范围内的预留面积存在先减小而后增大的过 程,延时间隔 ０ ｍｓ 时刻的上部面积最大,６ ｍｓ 时刻 的上部面积最小 在台阶中部范围内,延时间隔 ２ ｍｓ 时刻在后部和前部都形成较大的预留区域,其他 情况只是在前部范围内出现较大的预留区域 延时 间隔 ４ ｍｓ 和 ６ ｍｓ 两种情况下,在台阶中部的后部 位置的预留区域的面积变化较为均匀；在台阶的中 后部位置,延时间隔 ０ ｍｓ 时刻的预留区域面积较小 [图 ８ａ中圆圈位置],存在一些过度破碎的情况 在台阶的底部,随着延时间隔的增大,根底的变化存 在较小的差异 ａ０ ｍｓ ｂ２ ｍｓ ｃ４ ｍｓ ｄ６ ｍｓ ｅ８ ｍｓ 图 ７ 台阶整体破碎形态图 Ｆｉｇ. ７ Ｂｅｎｃｈ ｆｒａｇｍｅｎｔａｔｉｏｎ .４１. 爆 破 器 材 Ｅｘｐｌｏｓｉｖｅ Ｍａｔｅｒｉａｌｓ 第 ４５ 卷第 ３ 期 ａ０ ｍｓ ｂ２ ｍｓ ｃ４ ｍｓ ｄ６ ｍｓ ｅ８ ｍｓ 图 ８ 两炮孔中点处切面 Ｆｉｇ. ８ Ｓｅｃｔｉｏｎ ｉｎ ｔｈｅ ｍｉｄｄｌｅ ｏｆ ｔｗｏ ｂｌａｓｔ ｈｏｌｅｓ 为了更直观地考察延时间隔对岩石破碎的影 响,将图 ８ 的预留区域的面积进行了统计,两炮孔中 点切面剩余的面积如图 ９ 所示 图 ９ 两炮孔中点切面剩余的面积 Ｆｉｇ. ９ Ｒｅｍａｉｎｉｎｇ ａｒｅａ ｏｆ ｔｈｅ ｓｅｃｔｉｏｎ ｉｎ ｔｈｅ ｍｉｄｄｌｅ ｏｆ ｔｗｏ ｂｌａｓｔ ｈｏｌｅｓ 从图 ９ 可以看出,延时间隔 ０ ｍｓ 时刻的预留面 积最小,延时间隔 ４ ｍｓ 时刻的预留面积最大,延时 间隔 ２ ｍｓ 时的预留面积次之,延时间隔 ６ ｍｓ 和 ８ ｍｓ 时的预留面积差别不是很大 在图 ９ 中,预留面 积最小值与最大值相差 １０％ 从而可以看出,合理 的毫秒延时间隔可以较好地改善岩石破碎情况,同 时大幅降低岩石的破碎块度比例范围,对于台阶上 部的岩石破碎情况有一定的改善作用,从而降低大 块的产生 根据文献[１２] 给出的爆生气体速度 Ｃｇ＝３００ ｍ/ ｓ,Ｃｐ＝４ ０００ ｍ/ ｓ,利用公式３计算出 延时间隔约为 ５. ７７ ｍｓ,这和数值模拟在 ６ ｍｓ 的情 况较为吻合 图１０给出了两孔顶部连线中点处图 ５ 中 Ｅ 点的损伤值随着时间变化关系 从图１０可以看 出,具有延时间隔的情况下,随着两孔间的延时间隔 增加,损伤值变化在前期较为一致；同时,４ ｍｓ和 ６ ｍｓ 两种情况的损伤变化曲线出现了重叠,说明二者 的损伤值随着时间的变化差别较小 图 １０ 两孔顶部连线中点处损伤曲线 Ｆｉｇ. １０ Ｄａｍａｇｅ ｃｕｒｖｅｓ ａｔ ｔｈｅ ｍｉｄｄｌｅ ｐｏｉｎｔ ｏｆ ｔｈｅ ｔｗｏ ｈｏｌｅｓ ｏｎ ｔｈｅ ｔｏｐ ｏｆ ｔｈｅ ｌｉｎｅ ３ 结论 利用一维拉格朗日波形理论,直观地揭示了延 时间隔能够改变岩石中某点的受力状态,并通过对 不同延时间隔情况下的台阶爆破进行数值模拟,进 一步明确了延时间隔对于岩石破碎具有较大的影 响 合理设置延时间隔对改善台阶上部岩石破碎情 况具有重要的影响,同时可以避免岩石过度破碎,进 一步优化岩石破碎块度比例范围 对在文中所采用参数的情况下的数值模拟结果 进行分析,同时,通过爆生气体和应力波关系计算出 的延时间隔和数值模拟得出的结果较为接近,有利 于岩石破碎的延时间隔为 ６ ｍｓ,其与抵抗线的关系 为 ２ ｍｓ/ ｍ 参 考 文 献 [１] ＳＴＡＧＧ Ｍ Ｓ,ＲＨＯＬＬ Ｓ Ａ. Ｅｆｆｅｃｔｓ ｏｆ ａｃｃｕｒａｔｅ ｄｅｌａｙｓ ｏｎ ｆｒａｇｍｅｎｔａｔｉｏｎ ｆｏｒ ｓｉｎｇｌｅ-ｒｏｗ ｂｌａｓｔｉｎｇ ｉｎ ａ ６. ７ ｍ ２２ ｆｔ ｂｅｎｃｈ[Ｃ] / / Ｆｒａｇｂｌａｓｔ ２. Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ２ｎｄ Ｉｎｔｅｒｎａ- ｔｉｏｎａｌ Ｓｙｍｐｏｓｉｕｍ ｏｎ Ｒｏｃｋ Ｆｒａｇｍｅｎｔａｔｉｏｎ ｂｙ Ｂｌａｓｔｉｎｇ. ＳＥＭ, Ｂｅｔｈｅｌ ＣＴ,１９８７２１０-２３０. [２] ＯＴＴＥＲＮＥＳＳ Ｒ Ｅ, ＳＴＡＧＧ Ｍ Ｓ,ＲＨＯＬＬ Ｓ Ａ. Ｃｏｒｒｅｌａ- .５１.２０１６ 年 ６ 月 精确延时对台阶爆破岩石破碎块度影响的数值模拟研究 李顺波,等 ｔｉｏｎ ｏｆ ｓｈｏｃｋ ｄｅｓｉｇｎ ｐａｒａｍｅｔｅｒｓ ｔｏ ｆｒａｇｍｅｎｔａｔｉｏｎ[Ｃ] / / Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ７ｔｈ ＩＳＥＥ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｅｘｐｌｏｓｉｖｅｓ ａｎｄ Ｂｌａｓｔｉｎｇ Ｒｅｓｅａｒｃｈ. Ｌａｓ Ｖｅｇａｓ,Ｎｅｖａｄａ, １９９１１７９- １８１. [３] 张乐,颜景龙,李风国,等. 隆芯 １ 号数码电子雷管在 露天采矿中的应用[Ｊ]. 工程爆破,２０１０,１６４７３- ７７. ＺＨＡＮＧ Ｌ,ＹＡＮ Ｊ Ｌ, ＬＩ Ｆ Ｇ, ｅｔ ａｌ. Ａｐｐｌｉｃａｔｉｏｎ ｏｆ Ｌｕｘ Ｎｏ. １ ｄｉｇｉｔａｌ ｅｌｅｃｔｒｏｎｉｃ ｄｅｔｏｎａｔｏｒ ｉｎ ａｎ ｏｐｅｎ ｐｉｔ ｍｉｎｅ[Ｊ]. Ｅｎｇｉｎｅｅｒｉｎｇ Ｂｌａｓｔｉｎｇ, ２０１０,１６４７３-７７. [４] ＰＡＬＥＹ Ｎ,ＡＬＡＳＫＡ Ｔ. Ｔｅｓｔｉｎｇ ｅｌｅｃｔｒｏｎｉｃ ｄｅｔｏｎａｔｏｒｓ ｔｏ ｉｎｃｒｅａｓｅ ＳＡＧ ｍｉｌｌ ｔｈｒｏｕｇｈ ｐｕｔ ａｔ ｔｈｅ Ｒｅｄ Ｄｏｇ Ｍｉｎｅ [Ｃ] / / Ｐｒｏｃｅｅｄｉｎｇ ｏｆ ｔｈｅ ３６ｔｈ Ａｎｎｕａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｅｘｐｌｏｓｉｖｅｓ ａｎｄ Ｂｌａｓｔｉｎｇ Ｔｅｃｈｎｉｑｕｅ. Ｓａｎ Ｄｉｅｇｏ, Ｃａｌｉｆｏｒ- ｎｉａ, ２０１０３０２-３１５. [５] ＫＯＥＮＩＧ Ｒ. Ｆｉｖｅ ｙｅａｒｓ ｅｘｐｅｒｉｅｎｃｅ ｗｉｔｈ ｔｈｅ ｄｙｎａｔｒｏｎｉｃ ｅｌｅｃｔｒｏｎｉｃ ｆｉｒｉｎｇ ｓｙｓｔｅｍ [Ｃ] / / Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ２４ｔｈ Ａｎｎｕａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｅｘｐｌｏｓｉｖｅｓ ａｎｄ Ｂｌａｓｔｉｎｇ Ｔｅｃｈ- ｎｉｑｕｅ. Ｎｅｗ Ｏｒｌｅａｎｓ, １９９８ ３５９-３６４. [６] ＬＥＷＩＳ Ｎ, ＰＥＲＥＩＲＡ Ｐ. Ｏｐｅｒａｔｉｎｇ ｉｍｐｒｏｖｅｍｅｎｔｓ ａｔ ｖｕｌｃａｎ ｍａｔｅｒｉａｌｓ ｍｃｃｏｏｋ ｑｕａｒｒｙ ｕｓｉｎｇ ｅｌｅｃｔｒｏｎｉｃ ｄｅｔｏｎａ- ｔｏｒｓ[Ｃ] / / Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ２９ｔｈ Ａｎｎｕａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｅｘｐｌｏｓｉｖｅｓ ａｎｄ Ｂｌａｓｔｉｎｇ Ｔｅｃｈｎｉｑｕｅ. Ｎａｓｈｖｉｌｌｅ, ２００３ １-１４. [７] ＭＣＫＩＮＳＴＲＹ Ｒ, ＦＬＯＹＤ Ｊ, ＢＡＲＴＬＥＹ Ｄ. Ｅｌｅｃｔｒｏｎｉｃ ｄｅｔｏｎａｔｏｒ ｐｅｒｆｏｒｍａｎｃｅ ｅｖａｌｕａｔｉｏｎ [Ｃ] / / Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ２８ｔｈ Ａｎｎｕａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｅｘｐｌｏｓｉｖｅｓ ａｎｄ Ｂｌａｓｔｉｎｇ Ｔｅｃｈｎｉｑｕｅ. Ｌａｓ Ｖｅｇａｓ, ２００２１-２０. [８] ＢＡＴＴＥＮ Ｃ, ＷＩＬＬＡＮ Ｌ, ＳＰＡＲＫＥＳ Ｍ. Ｐｒｏｄｕｃｔｉｖｅ ｐａｒｔｎｅｒｓｈｉｐ-ａｄｖａｎｃｅｄ ｂｌａｓｔｉｎｇ ｔｅｃｈｎｏｌｏｇｙ ａｔ Ｎｅｗｌａｎｄｓ Ｃｏａｌ Ｍｉｎｅ [ Ｃ ] / / Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ３７ｔｈ Ａｎｎｕａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｅｘｐｌｏｓｉｖｅｓ ａｎｄ Ｂｌａｓｔｉｎｇ Ｔｅｃｈｎｉｑｕｅ. Ｓａｎ Ｄｉｅｇｏ, Ｃａｌｉｆｏｒｎｉａ, ２０１０１２０-１３２. [９] ＰＥＴＲＯＰＯＵＬＯＳ Ｎ, ＪＯＨＡＮＳＳＯＮ Ｄ,ＯＵＣＨＴＥＲＬＯＮＹ Ｆ. Ｆｒａｇｍｅｎｔａｔｉｏｎ ｕｎｄｅｒ ｄｉｆｆｅｒｅｎｔ ｃｏｎｆｉｎｅｍｅｎｔ ｃｏｎｄｉｔｉｏｎｓ ａｎｄ ｔｈｅ ｂｕｒｄｅｎ ｂｅｈａｖｉｏｒｓｍａｌｌ ｓｃａｌｅ ｔｅｓｔｓ [ Ｍ ] / / ＢＬＡＮＣＯ Ｊ Ａ Ｓ,ＳＩＮＧＨ Ａ Ｋ. Ｍｅａｓｕｒｅｍｅｎｔ ａｎｄ Ａｎａｌｙｓｉｓ ｏｆ Ｂｌａｓｔ Ｆｒａｇｍｅｎｔａｔｉｏｎ. Ｂｏｃａ Ｒａｔｏｎ, ＦＬ ＣＲＣ Ｐｒｅｓｓ, ２０１２６１. [１０] ＳＣＨＩＬＬ Ｍ. Ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ｓｉｍｕｌａｔｉｏｎｓ ｏｆ ｂｌａｓｔｉｎｇ ａｎｄ ｔｈｅ ｅｆｆｅｃｔｓ ｏｆ ｐｒｅｃｉｓｅ ｉｎｉｔｉａｔｉｏｎ ｏｎ ｆｒａｇｍｅｎｔａｔｉｏｎ[Ｒ]. ＤＹＮＡ ｍｏｒｅ ｒｅｐｏｒｔ,２０１１. [１１] ＲＯＳＳＭＡＮＩＴＨ Ｈ Ｐ. Ｔｈｅ ｕｓｅ ｏｆ ｌａｇｒａｎｇｅ ｄｉａｇｒａｍｓ ｉｎ ｐｒｅｃｉｓｅ ｉｎｉｔｉａｔｉｏｎ ｂｌａｓｔｉｎｇ. Ｐａｒｔ Ｉｔｗｏ ｉｎｔｅｒａｃｔｉｎｇ ｂｌａｓｔ- ｈｏｌｅｓ[Ｊ]. ＦｒａｇｂｌａｓｔＩｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｆｏｒ Ｂｌａｓｔｉｎｇ ａｎｄ Ｆｒａｇｍｅｎｔａｔｉｏｎ,２００２,６１１０４-１３６. [１２] ＣＨＯ Ｓ Ｈ, ＫＡＮＥＫＯ Ｋ. Ｒｏｃｋ ｆｒａｇｍｅｎｔａｔｉｏｎ ｃｏｎｔｒｏｌ ｉｎ ｂｌａｓｔｉｎｇ[Ｊ]. Ｍａｔｅｒｉａｌｓ Ｔｒａｎｓａｃｔｉｏｎｓ,２００４,４５１７２２- １７３０ Ｎｕｍｅｒｉｃａｌ Ｓｉｍｕｌａｔｉｏｎ ｏｆ ｔｈｅ Ｉｍｐａｃｔ ｏｆ Ｐｒｅｃｉｓｅ Ｔｉｍｅ Ｄｅｌａｙ ｏｎ Ｒｏｃｋ Ｆｒａｇｍｅｎｔａｔｉｏｎ ｉｎ Ｂｅｎｃｈ Ｂｌａｓｔｉｎｇ ＬＩ Ｓｈｕｎｂｏ①, ＹＡＮＧ Ｒｅｎｓｈｕ①,ＹＡＮＧ Ｊｕｎ② ①Ｓｃｈｏｏｌ ｏｆ Ｍｅｃｈａｎｉｃ ＆ Ｃｉｖｉｌ Ｅｎｇｉｎｅｅｒｉｎｇ, Ｃｈｉｎａ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｍｉｎｉｎｇ ＆ Ｔｅｃｈｎｏｌｏｇｙ Ｂｅｉｊｉｎｇ Ｂｅｉｊｉｎｇ, １０００８３ ②Ｓｔａｔｅ Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｅｘｐｌｏｓｉｏｎ Ｓｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ, Ｂｅｉｊｉｎｇ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｔｅｃｈｎｏｌｏｇｙ Ｂｅｉｊｉｎｇ, １０００８１ [ＡＢＳＴＲＡＣＴ] Ｉｎ ｏｒｄｅｒ ｔｏ ｓｔｕｄｙ ｅｆｆｅｃｔｓ ｏｆ ｔｈｅ ｐｒｅｃｉｓｅ ｔｉｍｅ ｄｅｌａｙ ｉｎｔｅｒｖａｌ ｏｎ ｒｏｃｋ ｆｒａｇｍｅｎｔａｔｉｏｎ, ｔｈｅ ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ｓｏｆｔ- ｗａｒｅ ｗａｓ ｕｓｅｄ ｔｏ ｓｉｍｕｌａｔｅ ｔｈｅ ｒｏｃｋ ｆｒａｇｍｅｎｔａｔｉｏｎ ｕｎｄｅｒ ｄｉｆｆｅｒｅｎｔ ｔｉｍｅ ｄｅｌａｙ ｉｎｔｅｒｖａｌｓ. Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｒｅａｓｏｎａｂｌｅ ｄｅｌａｙ ｔｉｍｅ ｓｅｔｔｉｎｇ ｉｓ ａｄｖａｎｔａｇｅｏｕｓ ｔｏ ｃｏｎｔｒｏｌ ｔｈｅ ｒｏｃｋ ｆｒａｇｍｅｎｔａｔｉｏｎ ｏｎ ｔｏｐ, ａｎｄ ｃａｎ ｉｍｐｒｏｖｅ ｔｈｅ ｆｒａｇｍｅｎｔａｔｉｏｎ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｒｏｃｋ. Ｔｈｒｏｕｇｈ ｔｈｅ ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｂｅｔｗｅｅｎ ｓｔｒｅｓｓ ｗａｖｅ ａｎｄ ｄｅｔｏｎａｔｉｏｎ ｇａｓ, ｉｔ ｉｓ ｆｏｕｎｄ ｔｈａｔ ｔｈｅ ｒｅｓｕｌｔｓ ｏｆ ｄｅｌａｙ ｉｎｔｅｒｖａｌ ｔｉｍｅ ｉｓ ｃｌｏｓｅ ｔｏ ｔｈｏｓｅ ｏｆ ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ. Ｔｈｅ ｄｅｌａｙ ｔｉｍｅ ｉｎｔｅｒｖａｌ ｉｓ ６ ｍｓ,ａｎｄ ｔｈｅ ｒｅｓｉｓｔａｎｃｅ ｌｉｎｅ ｉｓ ３ ｍ. Ａｎｄ ｗｈｅｎ ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｂｅｔｗｅｅｎ ｄｅｌａｙ ｔｉｍｅ ｉｎｔｅｒｖａｌ ａｎｄ ｒｅｓｉｓｔａｎｃｅ ｉｓ ２ ｍｓ/ ｍ, ｉｔ ｉｓ ｃｏｎｄｕｃｉｖｅ ｔｏ ｔｈｅ ｒｏｃｋ ｆｒａｇｍｅｎｔａｔｉｏｎ ｕｎｄｅｒ ｔｈｉｓ ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ. [ＫＥＹ ＷＯＲＤＳ] ｒｏｃｋ ｆｒａｇｍｅｎｔａｔｉｏｎ； ｐｒｅｃｉｓｅ ｔｉｍｅ ｄｅｌａｙ； ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ .６１. 爆 破 器 材 Ｅｘｐｌｏｓｉｖｅ Ｍａｔｅｒｉａｌｓ 第 ４５ 卷第 ３ 期