基于离散元方法的急倾斜极薄矿脉机械化开采爆破参数仿真优化.pdf

第 卷第期 年月有色金属(矿山部分)N o n f e r r o u sM e t a l s(M i n i n gS e c t i o n)V o l ,N o M a y d o i:/j i s s n 基于离散元方法的急倾斜极薄矿脉机械化开采爆破参数仿真优化熊有为,朱建国,刘福春,张军,陈伟(长沙有色冶金设计研究院有限公司,长沙 ;中南大学 资源与安全工程学院,长沙 ;广西桂华成有限责任公司,广西 贺州 )摘要:急倾斜极薄矿脉留矿法开采爆破空间规整平齐是实现小型采矿机器人高效通行和稳定作业的关键因素.为了提高狭小空间爆破光面控制效果,采用离散元方法和黏结颗粒模型耦合仿真方法,开展了不同炮孔参数条件下急倾斜极薄矿脉爆破模拟,以爆破空间轮廓规整性系数为评价指标,对急倾斜极薄矿脉炮孔布置方式和参数进行了优化.研究结果表明:同等排距和孔距条件下,矩形方式布孔的爆破规整性系数最高,矩形间隔空孔布孔方式爆破规整性系数次之,之字形炮孔布置方式爆破效果最差;结合爆破经济性指标和采场断面规整性系数综合考虑,对于急倾斜极薄矿脉机械化开采,推荐采用矩形间隔空孔装药的爆破方式,炮孔排距 m,孔距 m.关键词:急倾斜极薄矿脉;机械化开采;留矿法;爆破;离散元仿真中图分类号:T D 文献标志码:A文章编号:()S i m u l a t i o na n do p t i m i z a t i o no fb l a s t i n gp a r a m e t e r s f o rm e c h a n i z e dm i n i n go f s t e e p l y i n c l i n e dt h i nv e i n sb a s e do nd i s c r e t ee l e m e n tm e t h o dX I ONGY o u w e i,Z HUJ i a n g u o,L I UF u c h u n,Z HANGJ u n,CHE N W e i(C I N FE n g i n e e r i n gC o,L t d,C h a n g s h a ,C h i n a;S c h o o l o fR e s o u r c ea n dS a f e t yE n g i n e e r i n g,C e n t r a lS o u t hU n i v e r s i t y,C h a n g s h a ,C h i n a;G u a n g x iG u i h u a c h e n gC o,L t d,H e z h o uG u a n g x i ,C h i n a)A b s t r a c t:I t i s t h ek e yf a c t o r t or e a l i z e t h eh i g he f f i c i e n c ya n ds t a b l eo p e r a t i o no f s m a l l s i z e dm i n i n gr o b o t s t h a tt h eb l a s t i n gs p a c ei sr e g u l a ra n de v e ni ns t e e pa n de x t r e m e l yt h i nv e i ns h r i n k a g em i n i n g T oi m p r o v et h es m o o t hs u r f a c ec o n t r o le f f e c to fb l a s t i n gi nn a r r o ws p a c e,t h ed i s c r e t ee l e m e n tm e t h o d(D EM)a n dt h ec o u p l e ds i m u l a t i o nm e t h o do fb o n d e dp a r t i c l e m o d e la r ea d o p t e d,t h eb l a s t i n ge f f e c to fs t e e p l yi n c l i n e de x t r e m e l yt h i nv e i n su n d e rd i f f e r e n tb l a s t i n gp a r a m e t e r sw a ss t u d i e d,a n dt h eb l a s t i n gl a y o u ta n dp a r a m e t e r sw e r eo p t i m i z e db yt a k i n gt h er e g u l a r i t yo fb l a s t i n gs p a t i a lp r o f i l ea st h ee v a l u a t i o ni n d e x T h er e s u l t ss h o wt h a t,u n d e rt h es a m er o wa n dh o l es p a c i n g,t h er e g u l a r i t yc o e f f i c i e n to ft h er e c t a n g u l a rh o l ea r r a n g e m e n t i st h eh i g h e s t,a n dt h a to ft h er e c t a n g l ea n di n t e r v a l e m p t yh o l e i s t h es e c o n d C o m b i n e dw i t ht h ee c o n o m i c i n d e xo fb l a s t i n ga n dt h er e g u l a r i t yc o e f f i c i e n t o f t h es t o p es e c t i o n,t h eb l a s t i n g m e t h o do fr e c t a n g l ea n di n t e r v a le m p t yh o l ec h a r g ei sr e c o mm e n d e d,t h eb e s tm e s hp a r a m e t e r so f t h eh o l e s i s m m K e yw o r d s:s t e e pe x t r e m e l y t h i nv e i n;m e c h a n i z e dm i n i n g;s h r i n k a g em e t h o d;b l a s t i n g;d i s c r e t e e l e m e n t s i m u l a t i o n收稿日期:基金项目:湖南省科技创新计划资助项目(R C )作者简介:熊有为(),男,高级工程师,博士,采矿工程专业,主要从事金属非金属矿山安全高效开采技术研究,E m a i l:x y w c i n f c o m.通信作者:刘福春(),男,正高级工程师,采矿工程专业,主要从事金属非金属矿山安全高效开采技术研究,E m a i l:l f c s i n a c o m.长期以来,金属矿山急倾斜极薄矿脉主要采用浅孔留矿法或削壁充填法进行开采 .该类方法以人工凿岩方式分层上行浅孔回采,有利于控制矿石损失和贫化,但由于机械化程度极低,导致劳动强度大,作业效率低,安全风险高,严重制约矿山安全第期熊有为等:基于离散元方法的急倾斜极薄矿脉机械化开采爆破参数仿真优化高效生产.随着机器人技术的不断发展,小型采矿机器人进入采场作业成为可能 .然而,急倾斜极薄矿脉开采幅宽极窄,在受限窄小空间内爆破夹制作用显著,若爆破参数控制不当,极易造成超挖或欠挖,导致采场边帮凹凸不平,空间轮廓规整性差,影响机械化设备的高效通行和稳定作业 .在矿岩爆破参数优化的研究工作中,数值模拟是一种行之有效的研究手段 .安龙等 采用L S D YNA数值仿真方法,分析了不同爆破参数下薄矿脉中深孔爆破围岩裂纹演化规律,优化了薄矿体中深孔爆破参数;马春德等 利用非线性动力分析软件AN S Y S/L S D YNA对连续段毫秒微差爆破、跳段毫秒微差爆破和瞬间爆破方式进行了数值模拟对比分析,提出了一种浅孔多段毫秒微差爆破技术.然而,以往针对急倾斜薄矿体爆破方面的研究,主要是以控制围岩扰动、减少矿石贫化、提高爆破效率为目标进行爆破参数优化.对于急倾斜极薄矿脉机械化开采,如何通过合理的炮孔布置以提高爆后空间轮廓规整性是保障机械化设备充分发挥效能的关键因素,目前对于该方面的研究甚少,难以有效指导现场生产.本文通过采用离散元和黏结颗粒模型耦合仿真方法,研究不同炮孔参数条件下急倾斜极薄矿脉爆破效果,以爆破空间轮廓规整性系数为评价指标,优化急倾斜极薄矿脉精细化控制爆破参数,为急倾斜极薄矿脉机械化开采提供技术支撑.离散元方法基本原理离散元方法(D EM)基本原理是将系统划分为许多小体积或小颗粒,并对每个小颗粒进行运动方程的求解,通过模拟颗粒之间以及颗粒与外界的相互作用力,从而模拟整个系统的行为.离散元中描述颗粒之间的相互作用力的模型为H e r t z M i n d l i n模型,模型作用力由切向与法向两部分组成.FFnFt()F(kn ni jnv ni j)(kt ti jtv ti j)()式中,F为颗粒间相互作用力,N;Fn与Ft分别为法向和切向作用力,N;kn与kt为法向与切向弹性系数,P a/m;ni j与 ti j为颗粒法向与切向重叠量面积,m;n与t为法向与切向阻尼系数,Nms;ni j与 ti j为法向与切向相对速度,m/s.现实中的固体颗粒一般为非球形,而离散元方法采用球形颗粒提升计算效率,因此需要补偿颗粒形状所致额外的滚动阻力.此处采用E l a s t i c P l a s t i cS p r i n g D a s h b o a r d模 型 补 偿 颗 粒 的 滚 动扭矩.Mkr,t tMkr,tkrrt;|Mkr,tkrrt|rR|Fn|rRFnMkr,tkrrtMkr,tkrrt o t h e r w i s e()式中,Mkr,tt与Mkr,t为颗粒滚动扭矩,二者相差一个时间步长t;kr与r为颗粒的滚动摩擦系数和静摩擦系数;r为颗粒的相对角速度,r a d/s;R为颗粒等效半径,m.离散元爆破仿真模型及参数设置 离散元爆破仿真模型适用于爆破过程模拟及岩石破碎过程的离散元数值模拟一般采用黏结颗粒模型(B o n d e dP a r t i c l eM e t h o d,B P M).B P M模型是P o t y o d y和C u n d a l l 为模拟岩石破碎而提出的,该方法不需要复杂的本构方程,只需设置黏结颗粒模型的微观属性参数去满足岩石宏观变化结果,因此被广泛应用到岩土工程、岩石力学、采矿工程等领域中.B PM模型的基本思想是将一群装配好的球形颗粒集合体,在各球的黏结点处通过添加平行键,形成可破碎的矿石颗粒模型.因此,它将矿石颗粒模型看成是由颗粒和平行键组成,分别由刚性颗粒单元和键来模拟,如图所示.颗粒间的力位移行为由法向刚度kn、切向刚度k、摩擦系数三个参数描述,平行键的力位移行为由单位面积上的法向、切向的刚度kn和kt,抗拉、抗压强度c和c,键的半径倍数五个参数描述.其中黏结键半径为两颗粒半径的m i n(R(A),R(B).两颗粒间的黏结键为弹性关系,它与颗粒自身的力位移为平行存在,称为平行键,因此颗粒间即使存在平行键也不能阻止颗粒间滑动.两者不同之处,在于平行键可传递力和力矩,而颗粒之间只传递力.平行键的总受力及扭矩由式()()定义,Fb ni jknAn()Fb ti jktAt()Tb ni jktJn()Tb ti jknJt()式()()中,kn和kt是键的法向和切向刚度;n、t、n、t分别是法向和切向的位移和旋有色金属(矿山部分)第 卷图黏结颗粒模型及平行键的力位移行为F i g T h e f o r c e d i s p l a c e m e n tb e h a v i o ro fb o n d e dp a r t i c l em o d e l a n dp a r a l l e l b o n d转.当平行键的张力或剪切力满足如式()所示的条件时,键断裂.m i nFb ni jATb ni jI,Fb ti jATb ti jJb()式()中,b为平行键的断裂阈值.离散元仿真参数及模型设置爆破过程中炸药对岩体所产生的破坏由动态冲击载荷及爆破后气体的准静压作用两部分组成.通过采用B PM模型,可模拟岩体在爆破压力下的破碎过程,而如何在选定的爆孔区域施加动态冲击载荷及爆破后的准静压表征目前无相应的离散元模型.因此,本文中爆破动、静压的施加通过在爆破瞬间增加爆孔内颗粒粒径的方式来仿真炸药所产生的瞬时能量.如图所示,建立与矿山机械化留矿法采场实际尺寸一致的离散元仿真模型.其中,采场作业空间(近似为一巷道)宽度为 m,与水平方向的夹角为 o;拟上向浅孔爆破岩体的厚度为 m.除采场作业空间外,都设置为 mm粒径的球形离散元颗粒.基于B PM模型的岩石模拟模型需定义多项数值模拟参数,如表所示,本文采用文献报导中描述类似岩体的B PM模型参数进行数值模拟.图离散元模拟爆破的截面及炮孔布置方式示意图F i g S e c t i o no fb l a s t i n g s i m u l a t e db yd i s c r e t e e l e m e n tm e t h o da n dt h ea r r a n g e m e n t o fb l a s th o l e s表离散元数值模拟参数设置T a b l eP a r a m e t e r s e t t i n go fd i s c r e t ee l e m e n tn u m e r i c a l s i m u l a t i o n离散元参数离散元颗粒参数值密度/(k gm)剪切刚度/P a 泊松比 静摩擦系数 恢复系数 滚动摩擦系数 颗粒黏结模型参数法向刚度/(Nm)切向刚度/(Nm)法向应力阈值/MP a 切向应力阈值/MP a 时间步长/s e针对数值模拟的设置流程,在仿真初始化后,图中的爆破相岩体首先被生成.根据不同的布孔方式,相应的爆孔位置中的颗粒粒径在起爆瞬间按照预先选定的值进行设置,产生爆破效果.爆破后,追踪并保存遗留岩体在选定位置的截面进行超挖与欠挖的评估.以实际爆破断面轮廓面积S与爆破理想断面轮廓面积S的比值作为评价爆破轮廓断面规整性的指标,即“爆破轮廓规整性系数S规”.S规S/S()第期熊有为等:基于离散元方法的急倾斜极薄矿脉机械化开采爆破参数仿真优化如图所示的离散元爆破数值模拟的示意结果,可通过选取特定的岩体界面来评估当前布孔方式及参数的超挖及欠挖的表现.图离散元模拟爆破的结果F i g T h e r e s u l t so fb l a s t i n g s i m u l a t i o nb yd i s c r e t e e l e m e n tm e t h o d(a)O v e r b r e a ka n du n d e r b r e a ka s s e s s m e n t;(b)S i m u l a t i o no fb l a s t i n g i nt h ed i r e c t i o no f s t o p e s t r i k e仿真结果分析 之字形布孔方式数值模拟结果采用之字形布孔方式时,通过采用不同的排距及孔距参数进行组合模拟,结果如图所示.结果表明,在同一排距工况下,逐步增加孔距,爆破后空间断面不断增大,欠挖量逐渐减小,爆破轮廓规整性系数不断提高,但总体均表现为欠挖现象,爆破轮廓规整性系数最大为 (图).在选定孔距工况下,逐渐增加排距,爆破后的空间轮廓的断面逐渐缩小,即欠挖量增加,这一现象是由于炮孔排距的增加导致单位体积矿石的炸药用量降低,爆破不充分.较多的欠挖轮廓将可能对机械化设备的通行效率和安全稳定性造成影响.图之字形布孔方式下排距为 m及不同孔距的数值模拟结果F i g N u m e r i c a l s i m u l a t i o nr e s u l t so f z i g z a gh o l ea r r a n g e m e n tw i t h mb o t t o mr o ws p a c i n ga n dd i f f e r e n th o l e s p a c i n g s图之字形布孔不同炮孔参数爆破规整性系数F i g T h eb l a s t i n gr e g u l a r i t yc o e f f i c i e n t o fd i f f e r e n tp a r a m e t e r so f z i g z a gh o l e s 矩形布孔方式数值模拟结果采用矩形布孔方式时,通过采用不同的排距及孔距参数进行组合模拟,结果如图所示.结果表明,在同一排距工况下,逐步增加孔距,爆破后壁面逐渐从欠挖转变为超挖,爆破轮廓规整性系数达到 .造成这一现象的原因是炮孔孔距增大,使炮孔更加接近理想轮廓边界,爆破冲击波对采场边界破坏性增加.在一定孔距工况下,逐渐增加排距,爆破后的壁面超挖程度有所降低(欠挖程度有所升高),如图所示.相比之字形,在同等排距及孔距设置条件下,矩形布孔方式的爆破轮廓规整性系数明显提升.从图中可以看出,当排距为 m,孔距为 m时,矩形布孔爆破轮廓规整性系数可达 ,此时爆破后的空间轮廓规整性较好,但还需结合炸药单耗量和爆破块度均匀性指标,对其参数合理性进行综合评价.有色金属(矿山部分)第 卷图矩形布孔方式下排距为 m及不同孔距的数值模拟结果F i g N u m e r i c a l s i m u l a t i o nr e s u l t so f mb o t t o mr o ws p a c i n ga n dd i f f e r e n th o l e s p a c i n g i nr e c t a n g u l a rh o l e l a y o u t图矩形布孔不同炮孔参数爆破规整性系数F i g B l a s t i n gr e g u l a r i t yc o e f f i c i e n t o fd i f f e r e n tp a r a m e t e r so f r e c t a n g u l a rh o l e s 矩形间隔空孔布孔方式数值模拟结果采用矩形间隔空孔的布孔方式,通过采用不同的排距及孔距参数组合进行模拟,结果如图所示.结果表明,固定排距条件下增加孔距,或固定孔距条件下增加排距,爆破轮廓规整性系数均未出现图矩形间隔空孔布孔,排距为 m不同孔距的数值模拟结果F i g N u m e r i c a l s i m u l a t i o no f mr o ws p a c i n gw i t hd i f f e r e n th o l e s p a c i n go f r e c t a n g u l a r a n d i n t e r v a l o f e m p t yh o l e s线性 变 化(图);在 炮 孔 排 距 为 m、孔 距 为 m时,爆破轮廓规整性系数达到 ;其次为炮孔排距为 m、孔距为 m时,爆破轮廓规整性系数为 ,其他工况条件下爆破轮廓规整性系数均不高.图矩形间隔空孔布孔不同炮孔参数爆破规整性系数F i g B l a s t i n gr e g u l a r i t yc o e f f i c i e n t o fd i f f e r e n tp a r a m e t e r so f r e c t a n g u l a ra n d i n t e r v a l o f e m p t yh o l e s图 排距 m不同布孔方式条件下爆破规整性系数F i g T h e r e g u l a r i t yc o e f f i c i e n t o fb l a s t i n gw i t h md i s t a n c ea n dd i f f e r e n th o l ea r r a n g e m e n t如图 所示,同等排距和孔距条件下,矩形方式布孔的爆破规整性系数均为最高;当排距为 m和 m时,矩形间隔空孔比之字形布孔方式的爆破规整性系数要高,其爆破效果更好;当排距为 m和 m时,矩形间隔空孔比之字形布孔方式的爆破规整性系数相差不大.对于急倾斜极薄矿脉,矩形布孔方式由于炮孔装药密集,爆破炸药单耗较高.结合爆破经济性指标和采场断面规整性系数综合考虑,对于急倾斜极薄矿脉机械化开采,推荐采用矩形间隔空孔装药的爆破方式,炮孔排距 m,孔距 m.第期熊有为等:基于离散元方法的急倾斜极薄矿脉机械化开采爆破参数仿真优化图 排距 m不同布孔方式条件下爆破规整性系数F i g T h er e g u l a r i t yc o e f f i c i e n t o fb l a s t i n gw i t h md i s t a n c ea n dd i f f e r e n th o l ea r r a n g e m e n t图 排距 m不同布孔方式条件下爆破规整性系数F i g T h er e g u l a r i t yc o e f f i c i e n t o fb l a s t i n gw i t h md i s t a n c ea n dd i f f e r e n th o l ea r r a n g e m e n t图 排距 m不同布孔方式条件下爆破规整性系数F i g T h er e g u l a r i t yc o e f f i c i e n t o fb l a s t i n gw i t h md i s t a n c ea n dd i f f e r e n th o l ea r r a n g e m e n t结论)采用离散元方法和黏结颗粒模型B PM耦合,实现了急倾斜极薄矿脉不同炮孔参数条件下的爆破仿真;提出以爆破规整性系数作为爆破质量评价指标,适应机械化留矿法小型采矿机器人高效通行的工艺特点.)在同等排距和孔距条件下,矩形方式布孔的爆破 规 整 性 系 数 均 为 最 高;当 排 距 为 m和 m时,矩形间隔空孔比之字形布孔方式的爆破规整性系数明显要高,其爆破效果更好;当排距为 m和 m时,矩形间隔空孔比之字形布孔方式的爆破规整性系数相差不大.)结合爆破经济性指标和采场断面规整性系数综合考虑,对于急倾斜极薄矿脉机械化开采,推荐采用矩 形间 隔 空 孔 装 药 的 爆 破 方 式,炮 孔 排 距 m,孔距 m.参考文献周佳琦,任玉东黄金矿山急倾斜薄矿脉开采工艺及发展趋势J黄金,():Z HOUJ i a q i,R E NY u d o n g M i n i n g t e c h n o l o g ya n dd e v e l o p m e n t t r e n do fs t e e p l yi n c l i n e dv e i n s i ng o l dm i n e sJG o l d,():张庆旭,刘东锐,李士超典型急倾斜薄矿脉浅孔回采工艺优化J有色金属(矿山部分),():,Z HAN G Q i n g x u,L I U D o n g r u i,L IS h i c h a o O p t i m i z a t i o no fs h a l l o wh o l es t o p i n gt e c h n o l o g y i nt y p i c a l s t e e pt h i nv e i n sJN o n f e r r o u sM e t a l s(M i n i n gS e c t i o n),():,郑相悦,关贤欢,荣帅急倾斜薄矿脉采矿方法 优选模型研究J化工矿物与加工,():Z HE N GX i a n g y u e,GUAN X i a n h u a n,R ON GS h u a i S t u d yo no p t i m i z a t i o nm o d e lo fm i n i n gm e t h o df o rs t e e p l yi n c l i n e dt h i no r ev e i nJI n d u s t r i a lM i n e r a l s&P r o c e s s i n g,():贾小东,韩国生急倾斜薄矿脉探矿与回采工艺实践J甘肃冶金,():J I A X i a o d o n g,HAN G u o s h e n g P r a c t i c eo fp r o s p e c t i n ga n dm i n i n gt e c h n o l o g yf o r s t e e p l y i n c l i n e da n dt h i nv e i nJ G a n s uM e t a l l u r g y,():胡际平国外小型地下采掘设备的发展与应用J国外金属矿山,():HUJ i p i n g D e v e l o p m e n ta n da p p l i c a t i o no f s m a l lu n d e r g r o u n dm i n i n ge q u i p m e n ta b r o a dJ M i n i n gE n g i n e e r i n g,():刘翔宇论采矿设备的自动化J科技创新导报,():L I U X i a n g y u O nt h ea u t o m a t i o no fm i n i n ge q u i p m e n tJS c i e n c ea n dT e c h n o l o g yI n n o v a t i o nH e r a l d,():张芫涛,林卫星,欧任泽,等深埋复杂倾斜薄矿体高效回采方案研究J矿业研究与开发,():Z HAN G Y u a n t a o,L I N W e i x i n g,OU R e n z e,e ta l S t u d yo ne f f i c i e n tm i n i n gs c h e m ef o rd e e p b u r i e dc o m p l e xi n c l i n e dt h i no r e b o d yJ M i n i n gR e s e a r c ha n dD e v e l o p m e n t,():有色金属(矿山部分)第 卷郑迪,陈才贤,杨福斗,等某钨矿近水平薄矿体采矿方法选择及开采保障措施J矿业研究与开发,():Z HE N GD i,C HE N C a i x i a n,YAN GF u d o u,e ta l S e l e c t i o no fm i n i n gm e t h o da n d m i n i n gs a f e g u a r d m e a s u r e sf o ran e a r l y h o r i z o n t a l t h i no r e b o d y i nat u n g s t e no r eJ M i n i n gR e s e a r c ha n dD e v e l o p m e n t,():张莉鹏缓倾斜至倾斜薄矿体的进路掏槽式采矿法J现代矿业,():Z HAN G L i p e n g T u n n e lc u t t i n g m i n i n g m e t h o df o r m i n i n gl o w a n g l ed i pt od i pt h i no r eb o d i e sJ M o d e r nM i n i n g,():廖九 波,范 文 录急 倾 斜 极 薄 矿 脉 机 械 化 开 采 新 技 术 及 应用J矿业研究与开发,():L I A OJ i u b o,F AN W e n l u An e w m a c h i n e dm i n i n gt e c h n o l o g yo fs t e e p l y i n c l i n e d a n d e x t r e m e l y t h i n o r e v e i n a n d i t sa p p l i c a t i o nJM i n i n g R e s e a r c h a n d D e v e l o p m e n t,():柏杨,周亚博,崔越,等急倾斜薄矿体脉外采准上向水平分层充填采矿法研究J矿业研究与开发,():B O Y a n g,Z HOU Y a b o,C U I Y u e,e ta l S t u d yo nu p w a r dh o r i z o n t a l c u t a n df i l lm i n i n gm e t h o dw i t hm i n i n gp r e p a r a t i o no u to fv e i n si ns t e e p l y i n c l i n e dt h i n o r eb o d yJM i n i n gR e s e a r c ha n dD e v e l o p m e n t,():A L E X AN D E RPN,C L AUD I OTV,L U I SAG,e t a l P r e s p l i tb l a s t i n g d e s i g n t o r e d u c e c o s t s a n d i m p r o v e s a f e t y i nu n d e r g r o u n dm i n i n gJ E n e r g yR e p o r t s,():L IXD,L I UK W,YAN GJC,e t a l N u m e r i c a l s t u d yo nb l a s t i n d u c e df r a g m e n t a t i o ni n d e e pr o c k m a s sJ I n t e r n a t i o n a lJ o u r n a l o fI m p a c tE n g i n e e r i n g,():D O I:/j i j i m p e n g L IG W,D E N G G Z,MAJG N u m e r i c a lm o d e l l i n go ft h er e s p o n s eo fc e m e n t e dp a s t eb a c k f i l lu n d e rt h eb l a s t i n go fa na d j a c e n to r es t o p eJ C o n s t r u c t i o na n dB u i l d i n g M a t e r i a l s,():D O I:/j c o n b u i l d m a t 安龙,张家华,李元辉,等急倾斜薄矿脉夹制作用下中深孔爆破模拟与参数优化J东北大学学报(自然科学版),():ANL o n g,Z HAN GJ i a h u a,L IY u a n h u i,e ta l S i m u l a t i o na n dp a r a m e t e r so p t i m i z a t i o n o f m e d i u m d e e p h o l e b l a s t i n gi ns t e e p l yn a r r o w v e i n m i n eu n d e rc o n f i n e m e n tJ J o u r n a lo fN o r t h e a s t e r nU n i v e r s i t y(N a t u r a l S c i e n c e),():ANL,S UO R I N E N IFT,XUS,e ta l Af e a s i b i l i t ys t u d yo nc o n f i n e m e n te f f e c to n b l a s t i n g p e r f o r m a n c ei n n a r r o w v e i nm i n i n gt h r o u g hn u m e r i c a lm o d e l i n gJ I n t e r n a t i o n a lJ o u r n a lo fR o c kM e c h a n i c sa n dM i n i n gS c i e n c e s,():马春德,王业顺,胡顺喜,等毫秒微差爆破技术在深部缓倾斜薄矿体安全强采中的应用研究J矿冶工程,():MAC h u n d e,WAN GY e s h u n,HUS h u n x i,e t a l A p p l i c a t i o no fm i l l i s e c o n d d e l a yb l a s t i n g t e c h n o l o g ys a f ea n d i n t e n s i v em i n i n go fg e n t l y d i p p i n gt h i no r eb o d i e sJ M i n i n ga n dM e t a l l u r g i c a lE n g i n e e r i n g,():C HE N W,R O B E R T SA,KA T T E R F E L D A,e ta l M o d e l l i n gt h e s t a b i l i t yo f i r o no r eb u l kc a r g o e sd u r i n gm a r i n et r a n s p o r tJP o w d e rT e c h n o l o g y,:熊有为,刘福春,陈伟,等基于离散元与光滑粒子流体力学耦合的高浓度尾砂浓密过程仿真J金属矿山,():X I ON GY o u w e i,L I UF u c h u n,C HE N W e i,e ta l S i m u l a t i o no ft h i c k e n i n gp r o c e s sf o rh i g hs o l i d sc o n c e n t r a t i o nt a i l i n g su s i n gc o u p l e dD EM S P H m e t h o dJ M e t a l M i n e,():WE N S R I C H C,KA T T E R F E L D AR o l l i n g f r i c t i o n a s at e c h n i q u ef o r m o d e l l i n