A modified reaction model of aluminum dust detonation
-
摘要: 铝粉反应模型是对悬浮铝粉尘气-固两相爆轰进行数值模拟研究的关键。通过考虑铝粉燃烧产物氧化铝(Al2O3)在高温下的分解吸热反应,改进了铝粉的扩散燃烧模型。将该模型嵌入到三维的气-固两相爆轰数值计算程序中,分别对铝粉/空气混合物以及铝粉/氧气混合物的爆轰进行了数值模拟,计算得到的稳定爆轰波速度与实验结果、文献值均吻合较好,误差小于5.5%,表明改进的铝粉反应模型适用于不同氧化气体氛围中铝粉尘爆轰的模拟计算。此外,对两相爆轰参数及爆轰流场的物理量分布进行分析,获得了铝粉反应模型对爆轰波结构的影响规律。Abstract: The reaction model of aluminum particles is the key to successfully simulate the two-phase detonation of aluminum suspensions. In this study, by considering the endothermic decomposition reaction of the aluminum oxide (Al2O3) product at high temperature, a diffusion combustion model for the aluminum particles was improved and was incorporated into the homemade numerical code for 3D simulation of gas-solid two-phase detonation. The numerical program is based on the theory of two-phase flows, both gaseous and solid phases are assumed to be continuous media with inter-phase transfer of mass, momentum and energy. The system of 3D governing equations is solved in Cartesian x-y-z coordinates using an Eulerian grid, the numerical simulation code uses an explicit finite difference scheme based on the space-time conservation element and solution element (CE/SE) method, and the fourth order Runge-Kutta method is used to solve the source terms of the governing equations. In addition, the stability is assured by the Courant-Friedrichs-Lewy (CFL) criterion. Program parallelization is realized based on the message-passing-interface (MPI) technique, and the reliability of the program is demonstrated by simulating the shock tube problem successfully. Based on the program and the improved reaction model for the aluminum particles, numerical simulations for detonations of Al/air mixtures and Al/O2 mixtures were performed, respectively, the simulated results of the steady detonation wave speeds are in agreement with the experimental results or the literature value, with the error of less than 5.5%, which demonstrate the validity of the improved reaction model for Al suspensions detonation in different oxidizing atmosphere. Moreover, the detonation parameters and the distributions of the physical quantities around the detonation wave are analyzed, and the influence law of the reaction model on the detonation wave structure is obtained.
-
Key words:
- aluminum particles /
- reaction model /
- two-phase detonation /
- reverse reaction
-
立井井筒开挖主要施工方式钻爆法中,随着立井井筒施工综合机械化水平提升,炮孔深度随之不断增加。炮孔深度超过5 m时,现有的孔内连续装药的形式难以满足高效掘进[1-2]。为此,本文针对孔内分段装药结构,开展爆炸应变场以及裂隙场研究。
针对条形柱状药包爆炸研究,颜事龙等[3]采用高速摄影对有机玻璃中条形药包破碎区分布特征进行研究,认为条型药包粉碎区可分为柱部和端部两个区域,端部的破坏效应比柱部弱,并且从中部到端部发展速度逐渐减弱。龚敏等[4-5]运用三维动态光弹性爆破模型实验系统,首次对条型药包端部、中间及两端同时起爆条件下物理模型应力场进行研究,发现两端同时起爆及中间起爆的应力场强度的分布比较均匀。郭洋等[6]以高速摄影和数值模拟为手段,研究不同位置起爆柱状药包爆生裂纹扩展和损伤破碎特征,得出正反向起爆时炮孔周围产生多条径向裂纹,爆生裂纹面经历了“密集波纹状→光滑平坦→少量径向微裂纹→密集径向裂纹”的逐步转变的四个典型过程。杨仁树等[7]结合动态焦散线、超动态应变测试系统和数值计算对中间起爆柱状药包的应力应变场演化规律进行研究,得出中间起爆可以显著改善沿炮孔轴线爆炸应力场的均匀程度,有效减少岩石爆破大块率的产生。向文飞等采用Starfield迭加法与动力有限元法,研究起爆点数量与起爆点位置对条形药包爆炸应力场影响规律,认为合理安排起爆点的数量及位置、缩小条形药包完成爆轰的时间有利于改善爆破效果[8]。李启月等[9]通过数值计算和现场试验,发现合理分段装药高度和空孔参数,能够有效提高深孔爆破一次成井爆破效果。杨国梁等[10]提出了三种轴向分段装药结构,分别对正向、反向和两端起爆时的应力波传播规律进行研究,认为增加分段装药上分层药量应力分布更加均匀并且平均压力最大。胡涛等[11]利用数值计算对比深孔爆破连续装药结构和分段装药结构有效应力、自由面有效应力峰值及孔壁压力,结果表明分段装药结构降低了压缩应力波和爆轰气体产物作用于孔壁的初始压力及周围点振动速度。杨仁树等[12-13]将分段装药应用于井巷掘进掏槽爆破,提出了先裂后掏的分段装药爆破核心,认为上段炸药爆破后下段爆破提供新的自由面及补偿空间,从而减低炮孔底部岩石夹制作用,增加岩石的破碎程度,改善了掏槽爆破效果,提高井巷掘进效率。上述研究成果多是集中在连续装药数值计算、模型试验及现场试验等,少见基于分段装药条件下全场应变以及三维裂隙场及损伤特性,因此开展分段装药爆炸应变场与裂隙场分布规律研究,为深井开挖爆破参数优化提供理论指导与科学依据。
本文针采用数字图像相关(digital image correlation, DIC)分析方法,探讨分段装药下全场应变的演化过程,分析不同分段装药比例下爆炸全场应变传播规律;采用大型工业电子计算机断层扫描仪(CT)研究手段,建立连续装药、分段装药爆破前后“岩石-爆炸裂纹”三维重构模型,获得岩石爆炸三维裂隙网络位置、形态的空间分布情况,分析爆后岩石三维裂隙场的体分形维数与损伤度,探索分段装药对岩体损伤破裂程度的影响。
1. 分段装药全场应变特征分布研究
1.1 数字图像相关法实验系统
数字图像相关法(digital image correlation, DIC)是通过待测介质受力变形前后的灰度矩阵进行相关计算,通过对测点位置受力变形前后的空间位置信息计算,获得试件分析点处表面位移与变形应变的光学方法[14]。
实验所采用的实验设备为超高速数字图像相关实验系统,如图1所示。该实验系统主要包含超高速摄像机、爆炸与摄像同步控制系统、补光与拍摄同步控制系统、VIC-2D分析软件。超高速相机采集幅数为180幅,每秒最高可拍500万幅图片,完全满足爆炸荷载应变监测过程,本实验采样频率为200 000 s−1。同步控制系统为自主开发的多通道电火花发爆器,能够精准控制起爆时间,同时可以与照明系统、高速相机相连,调控三者的启动顺序,满足多通道同时、微差爆破。
1.2 分段装药实验模型设计
本实验旨在研究孔内分段装药全场应变特征,分段掏槽的目的是缓解炮孔底部岩体所受的夹制力作用。考虑在工程实践中,当上分段装药比例过小时,上部岩体类似浅孔爆破,下部岩体还处于高夹制作用状态,不符合设计初衷;当上分段比例过大,小分段比例过小时,下分段炸药能量有限,不能很好地把炮孔底部岩体抛掷出来。为此,本实验设计上分段装药占比0.4、0.6两种方案。
试件采用聚碳酸酯(PC)板。试件尺寸为400 mm×300 mm×8 mm,柱状药包尺寸为100 mm×3 mm×4 mm,炸药为叠氮化铅,爆速为4600 m/s,密度为4.71×103 kg/m3,装药量为160 mg。设计上分段装药占比为0.4、0.6两种方案,第一种方案,上分段炸药长L1=40 mm,如图2(a)所示,第二种方案,上分段炸药长L1=60 mm,如图2(b)所示,装药实物如图2(c)所示。对于两分段微差起爆时间的确定,在实验过程中,发现如果微差时间过大,那么两分段相互作用明显减弱;微差时间过小,两分段爆炸应力场的相互影响增大,不利于分析各分段产生的爆炸应力应变状态,根据前期实验调试结果,确定两段间隔起爆时差为20 μs。沿炮孔上部10 mm处,间隔20 mm布置6个测点测点,E1对应下分段L2起爆端(0 mm),E6位于上分段顶端(100 mm),对比分析其全场应变的演化规律,采用孔底起爆的方式。
1.3 分段装药全场应变分析
图3~4分别为两种分段装药的全场应变,分别选取径向应变(εyy)和轴向应变(εxx)演化过程进行分析。图3(a)为分段装药上分段占比0.6的径向全场应变,由于炸药起爆和传爆过程在炮孔中完成,高速相机无法记录此过程,实验捕捉到的应变过程为爆炸应力波开始作用介质的时刻。t=30 μs时,炮孔上下两侧受爆炸应力波作用表现为压应变,此时上分段炸药对下分段炸药及介质没有产生影响;t=40 μs时,随上分段炸药爆轰的传递,压应变作用区域增大;t=50 μs时,下分段炸药起爆,炮孔两侧同样表现为压应变,这时上分段炸药爆轰传播完成;t=70 μs时,上分段炸药起爆时产生的压应变继续向外传播,同时在炮孔周围形成第二组压应变场,两组压应变场之间为拉应变场,下分段压应变场作用范围增大,两组应变场之间开始相互作用;t=80 μs时,上分段炮孔周围压应变场增大,下分段炮孔周围压应变场范围有所减小,这是由于下分段炸药起爆形成的爆炸应力波向上分段传递,造成上分段应变区域增大,说明两者在此区间内有应变场的相互叠加;t=100 μs时,全场应变的分布状态以上分段为主。
图3(b)为上分段占比0.6的轴向全场应变,t=30 μs时,炮孔左右两侧受爆炸应力波作用表现为压应变;t=40 μs时,随上分段炸药爆轰的传递,压应变作用范围逐渐扩大;t=50 μs时,下分段炸药起爆,炮孔两侧同样表现为压应变,下分段压应变场与上分段压应变场在连接处发生叠加;t=70 μs时,在各分段的两端形成压应变场,其余炮孔位置表现为拉应变,t=100 μs时,上分段产生的压应变在介质中继续传播,随后紧跟的是拉应变作用。炮孔中由于分段装药形式,产生了三组压应变作用区域,分别为炮孔两端位置以及上下分段的连接处。这就表明在分段装药结构下,下分段堵塞区域主要为压应变作用区域。
图4(a)为分段装药上分段占比0.4的径向全场应变,t=30 μs时,炮孔受爆炸应力波作用表现为压应变,此时上分段炸药对下分段炸药及介质没有产生影响;t=40 μs时,随着上分段炸药爆轰的传递,压应变作用区域增大;t=50 μs时,下分段炸药起爆,炮孔两侧表现为压应变,这时上分段炸药爆轰传播完成;t=70 μs时,下分段应变场压应变作用区域增大,上分段应变场作用强度随着应力波的向外传递而减小;t=80 μs时,炮孔应变场主要集中在下分段装药处,应变场的作用范围较上分段大,直至100 μs时,由于下分段爆炸应力波的传播,全炮孔周边应变场作用范围扩大,但主要集中在下半段装药位置。
图4(b)为分段装药上分段占比0.4的轴向全场应变,t=30 μs时,炮孔左右两侧受爆炸应力波作用表现为压应变;t=40 μs时,随上分段炸药爆轰的传递,压应变作用范围同样逐渐增大;t=50 μs时,下分段炸药起爆,炮孔两端表现为压应变,下分段压应变场与上分段压应变场发生叠加;t=70 μs时,在各分段的两端形成压应变场,其余炮孔位置表现为拉应变,此时下分段形成的压应变场作用强度最大;t=100 μs时,和上分段装药0.6相比,此时炮孔的两端表现为压应变,中间区域则为拉应变。介质在此分段比例下炮孔周边区域沿装药方向主要受拉伸应变的作用,更有利于介质的破坏。
图5为孔内分段装药测点径向应变时程曲线,针对分段装药形式,无论上下分段占比如何,这种装药形式改变了炸药的爆炸应力波传播形式。针对上分段占比0.6的情况,上分段爆破的介质产生的压应变作用大于下分段,同时下分段产生的爆炸作用对上分段会产生二次应变,这种作用形式会加大对介质的破坏能力。但由于下半段炸药对介质产生的应变强度小,在工程实践中掏槽区域往往下半段岩石受夹制作用更大,岩石更难抛掷,所需的爆破能量更多。从介质的应变曲线可以看出,对于上分段装药0.4的情况,下分段介质受爆炸作用应变峰值更大,更好满足工程实践中下半段岩体对爆炸能量的需求。与文献[15]中连续装药爆炸应变场结论进行对比,分段装药改变了连续装药对介质的全场应变形态,由原来对介质产生一次应变改变为两次应变。在满足第一段炸药对介质的破坏作用下,同时加大了第二段炸药对介质的作用效应。
2. 孔内分段装药三维裂隙场的分布
2.1 孔内分段岩石爆破实验方案
为了探究分段装药爆炸裂隙场,采用红砂岩作为实验材料,红砂岩的单轴抗压强度为68.30 MPa,单轴抗拉强度为5.95 MPa,密度为2.48×103 kg/m3,如图6所示,红砂岩试样直径为50 mm,高为150 mm。在试件中心轴线钻取直径为4 mm,高度为120 mm的炮孔。连续装药时,采用装药系数66.7%,即柱状药包装药高度为80 mm,封堵为40 mm。采用分段装药时,根据以上研究结果,着重分析上分段装药占比为0.4的情况,即上分段装药中,炸药为30 mm,封堵为20 mm,下分段装药中,炸药为50 mm,封堵为20 mm。炸药装药直径为3 mm的细管,炸药为叠氮化铅。分段装药时,上段先爆,下段后爆,两者间隔20 μs。
红砂岩样品放置在内径50 mm的被动围压装置内,装置如图7所示,被动围压装置分内外两层,外层为法兰;里层为内胆装置,用于放置试件;把试件放置在内胆后,上部用法兰盖拧紧。
2.2 分段爆破岩石CT扫描与三维重构
CT扫描仪器型号ACTIS300-320/225,扫描电压为280 kV,试样扫描区段为0~150 mm,每隔0.15 mm扫描一层,共扫描1000张图片,CT扫描输出图像的尺寸为55 mm×55 mm,像素数量为1024×1024。
图8为不同层位的扫描原图和灰度处理图像,连续装药表面没有产生爆炸裂隙,分段装药试件表面产生了爆炸裂隙,从直观上看,分段装药对岩体的破碎程度更大一些。这是由于深孔装药下,在装药系数以及炸药量一定的情况下,连续装药封堵段长,炸药爆炸能量对封堵段破坏小,造成炸药的能量利用效率降低。由于分段装药改变了炸药在炮孔中的分布形式,使得炸药在炮孔中分布的更为均匀。
图9和图10分别为两种装药形式下岩体的三维裂隙重构图,由图可知,爆炸裂纹主要沿径向扩展,轴向应力应变所形成的环向裂纹不明显,径向是岩石破坏的主要方向。连续装药结构下,爆炸裂纹没有贯穿试件整体,距孔口20 mm内爆炸裂隙较少,这也说明在工程实践中,由于装药深度的增大,非装药段岩体破碎效果差,更容易产生大块。分段装药结构下,由于提高了炸药的位置,使得上部分岩体能够更好地利用炸药爆炸的能量破碎岩石。
2.3 分形维数计算与分析
目前对于分形维数的计算方法有计盒维数、相似维数、信息维数等,其中计盒维数计算简单,能够直观反映介质所选区域目标的占有情况。计盒维数表达式为[16]:
lgNδk=−Dflgδk+b (1) 式中:
Nδk 为含有裂纹区域的盒子数目,Df 为区域裂纹场的分形维数,δk 为裂纹区域分割小方网格边长,b为拟合参数。计盒维数的计算方法为:建立一个边长为δk 的小立体盒;然后改变不同边长对应形成若干小盒子,计算覆盖集的小盒子数Nδk ,经过转换得到δk -Nδk 数据;最后取对数,采用最小二乘法求斜率,即为计盒维数。材料损伤度
ω 与分形维数的关系表达式为[17]:ω=Dt−D0Dmaxt−D0 (2) 式中:
Dt 为介质爆后内部造成损伤面积的分形维数;D0 为介质爆前内部初始损伤面积的分形维数;Dmaxt 为介质达到最大损伤面积时的分形维数,对于平面问题Dmaxt=2 ;对于三维问题Dmaxt=3 。图11为连续装药与分段装药爆后岩体的分形维数图,从图中可以看出:分段装药分形维数为2.5073,对岩体造成的损伤度为0.84;连续装药岩体的分形维数为2.0258,对岩体造成的损伤度为0.68;分段装药下岩石的损伤度提高23.5%。
由于岩石上分段装药长度为50 mm,以上分段50 mm、下分段100 mm分别作对比,分别重构出裂隙的三维分布图,如图12所示。上分段岩体爆炸裂隙差别最大,分别计算上下分段两种装药结构的分形维数和损伤度,如图13所示,连续装药上分段分形维数为1.6889,损伤度为0.56;分段装药上分段分形维数为2.4693,损伤度为0.82,比连续装药上分段提高了46.4%;连续装药下分段分形维数为2.6679,损伤度为0.89;分段装药上分段分形维数为2.5228,损伤度为0.84,下半段连续装药损伤度大于分段装药,这是由于连续装药炸药集中在下部,单位体积的炸药比分段装药大,单两者从对岩体的损伤角度来看差别不大。
3. 结 论
本文分析了孔内分段装药结构下爆炸全场应变传播规律,建立了爆后“岩石-爆炸裂隙”的三维重构模型,精准描述了爆炸裂纹位置与形态的空间分布情况,结合分形理论,得到岩石材料爆炸裂隙的分形维数与损伤度,定量研究岩石在爆炸作用下的三维裂隙场与损伤程度,得到以下结论:
(1) 分段装药改变了连续装药对介质的全场应变形态,由原来对介质产生一次应变改变为两次应变;通过装药比例的改变,上分段装药占比0.4时,下分段爆破对介质产生的压应变作用效应大于上分段,在满足第一段炸药对介质的破坏作用下,加大了第二段炸药对介质的作用效应,同时延长了介质受爆炸应力波作用时间;
(2) 连续装药结构下,爆炸裂纹没有贯穿试件整体,炮孔封堵段内的爆炸裂纹较少,更容易产生大块;分段装药结构下,由于提高了炸药的位置,使得上部分岩体能够更好的利用炸药爆炸的能量破碎岩石;
(3) 分段装药对岩体造成的损伤度为0.8401,连续装药对岩体造成的损伤度为0.6802,分段装药岩石的损伤度较连续装药提高了23.5%,其中上分段岩体两者的差异性分析显示,上分段50 mm、下分段100 mm清醒,上分段爆炸产生的裂隙差别最大,分段装药上分段损伤度比连续装药提高46.4%。
-
图 4 3.83 ms时刻铝粉/空气爆轰波阵面附近气体压力、气体温度、气体密度、气体轴向速度、铝粉体积分数以及铝粉轴向速度的分布
Figure 4. Distributions of gas pressure, gas temperature, gas density, gas velocity in the x direction, volume fraction of Al particles, and Al particles velocity in the x direction around the detonation wave for Al/air mixtures at t=3.83 ms
图 7 4.43 ms时刻铝粉/氧气爆轰波阵面附近气体压力、气体温度、气体密度、气体轴向速度、铝粉体积分数以及铝粉轴向速度的分布
Figure 7. Distributions of gas pressure, gas temperature, gas density, gas velocity in the x direction as well as volume fraction of Al particles, and Al particles velocity in the x direction around the detonation wave for Al/O2 mixtures at t=4.43 ms
-
[1] LIU X L, ZHANG Q. Influence of turbulent flow on the explosion parameters of micro- and nano-aluminum powder-air mixtures [J]. Journal of Hazardous Materials, 2015, 299: 603–617. DOI: 10.1016/j.jhazmat.2015.07.068. [2] VEYSSIERE B, KHASAINOV B A, BRIAND A. Investigation of detonation initiation in aluminium suspensions [J]. Shock Waves, 2008, 18(4): 307–315. DOI: 10.1007/s00193-008-0136-z. [3] FEDOROV A V, KHMEL T A. Numerical simulation of formation of cellular heterogeneous detonation of aluminum particles in oxygen [J]. Combustion, Explosion, and Shock Waves, 2005, 41(4): 435–448. DOI: 10.1007/s10573-005-0054-7. [4] BECKSTEAD M W. Correlating aluminum burning times [J]. Combustion, Explosion and Shock Waves, 2005, 41(5): 533–546. DOI: 10.1007/s10573-005-0067-2. [5] TANGUAY V, GOROSHIN S, HIGGINS A J, et al. Aluminum particle combustion in high-speed detonation products [J]. Combustion Science and Technology, 2009, 181(4): 670–693. DOI: 10.1080/00102200802643430. [6] BAZYN T, KRIER H, GLUMAC N. Evidence for the transition from the diffusion-limit in aluminum particle combustion [J]. Proceedings of the Combustion Institute, 2007, 31(2): 2021–2028. DOI: 10.1016/j.proci.2006.07.161. [7] GLORIAN J, GALLIER S, CATOIRE L. On the role of heterogeneous reactions in aluminum combustion [J]. Combustion and Flame, 2016, 168: 378–392. DOI: 10.1016/j.combustflame.2016.01.022. [8] ZHANG F, GERRARD K, RIPLEY R C. Reaction mechanism of aluminum-particle-air detonation [J]. Journal of Propulsion and Power, 2009, 25(4): 845–858. DOI: 10.2514/1.41707. [9] BRIAND A, VEYSSIERE B, KHASAINOV B A. Modelling of detonation cellular structure in aluminium suspensions [J]. Shock Waves, 2010, 20(6): 521–529. DOI: 10.1007/s00193-010-0288-5. [10] BALAKRISHNAN K. Diffusion- and kinetics-limited combustion of an explosively dispersed aluminum particle [J]. Journal of Propulsion and Power, 2014, 30(2): 522–526. DOI: 10.2514/1.B35059. [11] KWON Y S, GROMOV A A, ILYIN A P, et al. The mechanism of combustion of superfine aluminum powders [J]. Combustion and Flame, 2003, 133(4): 385–391. DOI: 10.1016/S0010-2180(03)00024-5. [12] NIGMATULIN R I. Methods used in mechanics of continuous media for a description of multiphase mixtures [J]. Journal of Applied Mathematics and Mechanics, 1970, 34(6): 1097–1112. DOI: 10.1016/0021-8928(70)90174-7. [13] HAYNES W M. Handbook of chemistry and physics [M]. Florida: CRC Press, 2014. [14] LEVITAS V I, PANTOYA M L, CHAUHAN G, et al. Effect of the alumina shell on the melting temperature depression for aluminum nanoparticles [J]. The Journal of Physical Chemistry C, 2009, 113(32): 14088–14096. DOI: 10.1021/jp902317m. [15] 洪滔, 秦承森. 铝颗粒激波点火机制初探 [J]. 爆炸与冲击, 2003, 23(4): 295–299.HONG T, QIN C S. Mechanism of shock wave ignition of aluminum particle [J]. Explosion and Shock Waves, 2003, 23(4): 295–299. [16] PRICE E W. Combustion of metalized propellants [M]// KUO K K. Fundamentals of Solid-Propellant Combustion. New York: American Institute of Aeronautics and Astronautics, 1984: 479−513. DOI: 10.2514/5.9781600865671.0479.0513. [17] STEINBERG T A, WILSON D B, BENZ F. The combustion phase of burning metals [J]. Combustion and Flame, 1992, 91(2): 200–208. DOI: 10.1016/0010-2180(92)90100-4. [18] GLASSMAN I. Combustion of metals revisited thermodynamically [C]// Proceedings of the Eastern States Section of the Combustion Institute. Princeton: The Combustion Institute, 1993: 17−26. [19] 沈维道, 童钧耕. 工程热力学[M]. 4版. 北京: 高等教育出版社, 2007.SHEN W D, TONG J G. Engineering thermodynamics [M]. 4th ed. Beijing: Higher Education Press, 2007. [20] TORO E F. Riemann solvers and numerical methods for fluid dynamics: a practical introduction [M]. Berlin Heidelberg: Springer, 2009. DOI: 10.1007/b79761. [21] TULIS A J, SELMAN J R. Detonation tube studies of aluminum particles dispersed in air [J]. Symposium (International) on Combustion, 1982, 19(1): 655–663. DOI: 10.1016/s0082-0784(82)80240-3. [22] LIU L J, ZHANG Q, SHEN S L, et al. Evaluation of detonation characteristics of aluminum/JP-10/air mixtures at stoichiometric concentrations [J]. Fuel, 2016, 169: 41–49. DOI: 10.1016/j.fuel.2015.11.090. [23] BENKIEWICZ K, HAYASHI K. Two-dimensional numerical simulations of multi-headed detonations in oxygen-aluminum mixtures using an adaptive mesh refinement [J]. Shock Waves, 2003, 12(5): 385–402. DOI: 10.1007/s00193-002-0169-7. 期刊类型引用(2)
1. 陈程,李全明. 地下矿山巷道掘进凿岩爆破技术参数优化研究进展与展望. 中国矿业. 2023(10): 103-110 . 百度学术
2. 费鸿禄,纪海楠,山杰. 露天台阶水介质间隔装药结构优选及对比试验研究. 黄金科学技术. 2023(06): 930-943 . 百度学术
其他类型引用(1)
-