Deformation field in 316L stainless steel by single shot peening
-
摘要: 运用金属材料表面纳米化试验机对单个弹丸撞击316L不锈钢表面进行了撞击实验;采用激光共聚焦显微镜观察了弹坑的三维形貌,测量不同振动频率下弹坑的直径及离面位移;采用云纹干涉法对弹坑周围的面内应变场进行测量,并分析振动频率及撞击方式对弹坑尺寸、塑性应变大小以及塑性应变区范围的影响;采用有限元方法对单个弹丸垂直撞击试件表面的应变场进行数值模拟,与实验结果进行比较,分析了弹坑周围残余应力的分布。结果表明:随振动频率的增加,弹坑直径和离面位移都增加,频率在50~55Hz,弹坑直径有突变,离面位移和振动频率呈线性关系;振动频率越大,塑性应变越大,塑性应变分布范围均大于弹坑直径的2倍;同一振动频率下弹丸垂直撞击比倾斜撞击的塑性应变大,而塑性应变分布范围相差不大;面内残余应变场的数值模拟结果和实验结果吻合较好,最大误差小于10%。Abstract: The experiment of the single shot impacting the 316L stainless steel surface was carried out using a surface nano-crystallization testing machine. Three-dimensional morphology of the dimple was observed with a laser scanning confocal microscope, and the dimple's diameter and off-plane displacement in different vibration frequencies were also measured. The in-plane strain around the dimple was measured by moiré interferometry. The effect of the vibration frequency and the way of impacting on the dimple size, the plastic strain size and the plastic strain zone were also analyzed. In comparison with the experimental result, the strain field was simulated using the finite element method, and the distribution of the residual stress around the dimple was also analyzed. The result showed that the crater diameter and the off-plane displacement increases with the increase of the vibration frequency. When the frequency is from 50 to 55 Hz, the crater diameter experiences mutations. When the shot impacts the surface vertically, the greater the vibration frequency, the greater the plastic strain and plastic strain zone, and the plastic strain zone are two times larger than the crater diameter. The plastic strain by the vertical impact is slightly greater than the plastic strain by the oblique impact, but it has little effect on the plastic strain zone. The experimental U field strain is in fairly good agreement with the numerical simulation result, and the maximum error is less than 10%.
-
巷道掘进爆破中,一般只有一个自由面,所受夹制作用较大,通常会留有空孔为岩石破碎提供补偿空间[1]。而空孔数目、间距等参数是提高掘进进尺和掏槽取得成功的关键。目前一致认可的是补偿空间越大,掏槽效果越好,即增加空孔数量或空孔直径,然而大的补偿空间需要牺牲的是钻孔成本。那么,如何让有限的空孔实现最佳的爆破空间补偿效果,是掏槽爆破研究的热点问题。
就掏槽爆破参数优化问题,主要从数值模拟、现场实验等方面进行了研究。Li等[2]利用数值模拟手段验证了中空孔应力集中效应的存在;Zhang等[3]提出了掏槽爆破的一种先进的掏槽抛片(CCFT)概念;Sui等[4]研究了无限岩体在有自由面条件下裂纹扩展过程;单仁亮等[5]模拟研究了不同空孔孔径的直眼掏槽爆破效果,得出孔径对掏槽爆破的效果和振动均有重要影响;汪海波等[6]认为增大空孔直径可提高应力峰值和延长作用时间,空孔周围峰值为无空孔时的2.2倍;Li等[7]利用数值模拟手段研究得出两空孔间距越短,止裂作用越强,当两孔间距小于一定值时,可完全阻止外向裂纹;钟波波等[8]利用数值模拟方法,分析得出空孔有导向作用,且与空孔大小无明显相关;材料的非均匀性,对裂纹的扩展方式有显著影响;张召冉等[9]针对含空孔直眼掏槽参数确定没有可靠依据的问题,提出了动态碎胀系数的概念;Gao等[10]采用ANSYS/LS-DYNA软件对不同侧压力系数的复杂五孔掏槽爆破进行二维优化模拟研究,结果表明,设置延迟爆轰和孔径以及优化炮眼间距可以优化破岩效果;柴修伟等[11]利用数值模拟软件LS-DYNA模拟了不同数量的空孔情况下裂纹扩展情况。
鉴于以上研究成果,本文以在研的矿山巷道掏槽爆破掘进效率优化项目为依托,选择双大直径空孔(孔径100 mm)、侧向两孔起爆的掏槽方案,研究同等起爆条件下大直径空孔距离对掏槽槽腔断面大小的影响,以期确定最优的掏槽布孔方式,保证矿山巷道掘进进尺的同时,减少空孔钻凿数量。
1. 掏槽爆破的空孔效应
1.1 空孔的自由面作用
起爆后,冲击波随着传播的距离不断衰减,应力也不断下降,逐渐变为压缩波,当压缩波传到空孔孔壁会反射形成拉伸波,与应力波传到自由面的作用相同,因此空孔充当了自由面的作用,有利于岩石破碎。
根据动量守恒定律可推出[12]:
σ=ρcpvp (1) τ=ρcsvs (2) 式中:σ为纵波应力,MPa;τ为横波应力,MPa;ρ为岩石密度,kg/m3;cp为纵波波速,m/s;cs为横波波速,m/s;vp和vs分别为质点在波的传播方向和垂直与波的传播方向上的运动速度,m/s。
1.2 空孔的应力集中作用
力学计算和实验研究表明:当应力波传到空孔孔壁时,将产生反射拉应力波并在空孔附近形成拉应力区,最大拉应力峰值将出现在的两孔连心线方向上[13-14]。因此,空孔具有应力集中作用[6]。在某测点的径向压应力峰值为:
σrmax=p0(rr0)−a (3) 式中:σrmax为某点的爆炸径向应力,MPa;p0为孔壁初始应力峰值,MPa;r0为炮孔半径,m;r为某点到炮孔距离,m;a为应力波衰减指数。
切向拉应力峰值可通过径向压应力峰值得到:
σθmax=bσrmax (4) 式中:b为系数,在爆炸近区b≈1,b值与应力波传播距离成反比,并逐渐趋于稳定。
σθmax=μ1−μσrmaxb=μ1−μ (5) 式中:μ为岩石泊松比。
p0=14ρeD202ρmcpρeD0+ρmcp (6) 式中:ρe、ρm分别为炸药和岩石密度,kg/m3;D0为爆速,m/s;cp为纵波波速,m/s。
计算可得到某一点切向拉应力峰值σθmax,若切向拉应力峰值σθmax大于岩石的动态抗拉强度σθ,即满足岩石拉裂破坏条件,岩石破碎。
1.3 空孔的卸压作用
在复杂地压下,根据厚壁圆筒理论,爆炸与原岩应力场耦合为:
σ′rr=12(σh+σv)(1−r20r2)+12(σh−σv)(1−4r20r2+3r40r4)cos2θ+p0(r0r)2 (7) σ′rθ=12(σh+σv)(1+r20r2)−12(σh−σv)(1+3r40r4)cos[2(π2−θ)]−p0(r0r)2 (8) 式中:
σ′rr 、σ′rθ 分别为计算点的耦合径向和切向应力,MPa;σh、σv分别水平和垂直地应力,MPa;θ 为点和炮孔中心连线与水平方向的夹角,r为点到炮孔距离,m。若σh=σv=0(即无地压情况),则
σ′rr=p0(r0r)2 (9) σ′rθ=−p0(r0r)2 (10) 空孔使炮孔在径向的抵抗线要远小于轴向的抵抗线,因此,破碎岩体在爆生气体的推动作用下,优先向空孔方向移动,空孔提供的空间可以消掉一定的碎胀压力,有利于岩石抛掷出腔,即空孔的卸压作用。
由上文可知,大直径空孔直眼掏槽机理是应力集中、自由面、空孔的卸压作用共同作用的结果。当爆炸气体产生时,空孔起到自由面作用,为破碎岩体提供一定的碎胀空间,有利于破碎岩体在爆生气体的推动作用下抛出,形成槽腔,即实现空孔卸压作用。下文利用数值模拟手段研究双空孔间距对爆破效果的影响。
2. 数值分析模型的建立
2.1 模型尺寸及试验方案
(1) 补偿空间理论
补偿空间应满足矿岩碎胀后所需要的体积,不然就会发生挤死的现象,从而影响爆破效果。其应满足如下关系式:
S1×K≤S1+S2+S3 (11) 式中:S1为预爆岩体面积;S2为空孔面积;S3为装药孔面积;K 为岩石碎胀系数,K =1.5。
根据图1,可推导出装药孔与空孔距离的关系:
(LD+d2−πD28−πd28)K≤D+d2L+πD28+πd28 (12) L≤π4(D2+d2)(K+1)(D+d)(K−1) (13) 式中:L为空孔与装药孔的距离;d为装药孔直径;D为空孔直径。计算得到L<32.7 cm。
(2) 炮孔偏斜理论
为防止炮孔间贯穿,需要满足:
L>D+d2+2Hsina (14) 式中:H为炮孔深度,a为炮孔偏斜度。
现场使用KJ311型全液压掘进钻车,炮孔偏斜角度在1°以内,经计算得:L>18.9 cm。所以装药孔与空孔之间的距离范围为:18.9 cm<L<32.7 cm。
运用ANSYS/LSDYNA建立有限元二维模型,模型尺寸为4000 mm×4000 mm,边缘定义为无反射边界,装药炮孔直径为50 mm,大直径空孔直径为100 mm,装药方式为耦合装药。沿模型中心横向设置两个装药炮孔,同时起爆,孔距为40 cm,两空孔按模型中心竖向布置,空孔孔距用dv表示,通过上式计算装药孔与空孔之间孔距18.9 cm<L<32.7 cm。根据勾股定理原则可计算出dv最大为51.7 cm,本文将空孔孔距dv设置为15、25、35、45和55 cm,定义两空孔连心线的中点为P,各炮孔位置及模型网格划分如图2所示。
2.2 岩石模型
结合相关文献,本岩石模型选用更加适合岩体的爆炸冲击数值模拟的HJC本构模型[17-18]。根据大红山铜矿试验区域参数[19],确定HJC本构模型中的矿岩参数见表1。
表 1 大红山大理岩HJC本构模型参数Table 1. Dahongshan marble HJC constitutive model parametersρ0/(kg·m−3) ƒc/MPa A B C Smax G T D1 D2 2941 70.59 0.52 1.17 0.0163 4 22.27 7.68 0.036 1 pcrush/MPa µcrush plock/GPa µplock K1 K2 K3 EFmin N FS 23.65 0.00076 0.159 0.012 13 23 60 0.01 0.79 0.085 注:ρ0为岩石密度,ƒc为静态单轴抗压强度,T为岩石最大拉应力,A为无量纲粘性强度系数,B为无量纲压力硬化系数,C为应变率系数,N为压力硬化指数,Smax为最大无量纲的等效应力,D1、D2为损伤常数,EFmin为岩石最小塑性应变,µcrush为等效塑性应变增量,µplock塑性体积应变增量,K1、K2和K3分为压力常数,plock为压实后的静水压力,pcrush为弹性极限时静水压力值,G为剪切模量,FS失效参数。 2.3 炸药模型参数及状态方程
矿山作业使用的是1号岩石乳化炸药,选用*MAT_HIGH_EXPLOSIVE_BURN模型,采用EOS_JWL状态方程模拟炸药状态[20-21]
p1=A1(1−ωR1V)e−R2V+B1(1−ωR2V)e−R2V+ωE0V (15) 式中:p1为爆轰产物内部压力,V为相对体积,A1、B1、R1、R2和ω为炸药相关常数,E0为初始比内能。
本文采用的炸药模型参数见表2。
表 2 炸药的状态方程参数Table 2. HJC constitutive model parameters of explosiveρ1/(g·cm−3) D0/(m·s−1) A1/GPa B1/GPa R1 R2 ω 1.3 4500 214.4 0.182 4.2 0.9 0.15 注:D0为炸药爆速。 3. 数值计算结果及分析
3.1 爆破应力场分布情况
为了研究双空孔对直眼掏槽爆炸应力场的影响,模拟试验设置两个垂直掏槽孔装药起爆,炮孔直径50 mm,空孔直径100 mm,仅改变双空孔孔距,得到dv=15 cm、dv=25 cm、dv=35 cm、dv=45 cm和dv=55 cm五种情况下不同时刻的等效压力云图,见图3~图7,P点处应力时程曲线如图8所示。
分析压力云图可知:在应力波到达空孔前,五种模型压力云图的应力峰值基本一致;在应力波传播至空孔附近时,可明显观察出两空孔连线方向出现应力集中现象,应力波传播到中心P点处时两炮孔产生的应力波相互叠加,应力波在空孔壁形成反射拉伸波,随后出现不同的应力波叠加扰动,直至传播出视野。
通过观察不同dv模型P点处的应力时程曲线可以看出:受到的压应力峰值均为217 MPa,dv=15 cm、dv=25 cm、dv=35 cm、dv=45 cm 、dv=55 cm 时P点处拉应力峰值分别为197、180、165、160和156 MPa;即dv由15 cm增大到25 cm,拉应力峰值减小8.62%;dv由25 cm增大到35 cm,拉应力峰值减小8.3%;dv由35 cm增大到45 cm,拉应力峰值减小3.03%;dv由45 cm增大到55 cm,拉应力峰值减小0.25%。说明dv越小,P点受到的反射拉应力越大,这跟空孔自由面的距离有关。
3.2 岩石损伤破坏情况
只考虑双大直径空孔间距对爆破效果的影响时,空孔间距越小,P点处受到的拉应力峰值越大,中心岩体越容易破坏。为方便研究槽腔断面破坏情况,在k文件中添加*MAT_ADD_EROSION关键字[11],设定压力强度大于本文实际岩石力学参数,即单轴抗压、单轴抗拉、变形模量等时,模拟岩体的损伤与破坏,如图9~图13所示。
如图9所示,dv=15 cm时,空孔距离太近,补偿空间导向作用不明显;当 dv=25 cm时(图10),在t=0.05 ms时,装药孔周围岩石单元受压破坏,应力波刚刚达到空孔壁,在t=0.1 ms时应力波持续向外传播,装药孔周围岩石破坏向外扩展,空孔壁处应力波发生反射开始受拉破坏,出现空孔导向现象,使空孔与以破坏单元连接形成初步槽腔轮廓,在t=0.15 ms时,两装药孔与空孔之间完全贯穿,在t=0.2 ms时,槽腔轮廓面基本形成,直至0.25 ms时刻破坏终止。
当 dv=35 cm时(图11),炮孔周围裂纹发展更加充分,同时反射应力波对岩石的拉伸破坏减小,爆破后得到的槽腔轮廓面积较dv=25 cm时有所增大,中心区域虽有少量岩体,但基本贯通形成槽腔轮廓面。
当dv=45 cm时(图12),由0.25 ms时刻损伤破坏图可以看出两装药孔之间基本贯穿,但槽腔外断面仍存在大量未失效单元,说明掏槽区域仍有大量岩体存留,若应用于现场能会出现补偿空间不够的情况,不能保证掏槽效果;当dv=55 cm时(图13),空孔距离太大,爆生裂纹无法贯穿,无法形成有效槽腔。
岩石破碎区域若存在宽度为1单元的细长区域,即将该区域判定为裂隙,统计破碎区域面积时不计入统计。对模拟结果得到的槽腔外断面进行提取,即CAD等比例绘图计算,得到图14。根据图14可知,dv=15 cm模型爆破后形成槽腔断面宽度为63.03 cm,形成的槽腔断面积为0.1641 m2;dv=25 cm模型爆破后形成槽腔断面宽度为68.6 cm,形成的槽腔断面积为0.2116 m2;dv=35 cm模型爆破后形成的槽腔断面为67.4 cm,形成的槽腔断面面积为0.2436 m2;dv=45 cm模型爆破后槽腔断面内存在大量未失效岩石单元,面积为0.1740 m2,dv=55 cm模型爆破后形成槽腔断面宽度为62.27 cm,高度仅为16.65 cm,形成的槽腔断面积为0.0951 m2。可见双空孔方案中dv由15 cm增加到25 cm时,形成的槽腔断面增大了28.94%;dv由25 cm增加到35 cm时,槽腔断面增大15.1%;但当dv由35 cm增加到45 cm时,槽腔面积为0.1740 m2,形成的槽腔断面减小了17.8%;dv增大到55 cm时,槽腔面积为0.0951 m2,形成的槽腔断面减小了45.3%;即槽腔面积随空孔间距的增大先增加后减小,dv=35 cm时槽腔面积最大,确定为现场验证方案。
4. 爆破现场验证
4.1 炮孔布置
试验地点为大红山铜矿285中段、370水平,根据大红山铜矿已有资料显示,矿岩主要为大理岩,岩石力学参数见表1。装药孔直径为50 mm,空孔直径为100 mm,孔深为3.2 m。由于上述模拟得出dv=35 cm时形成的槽腔断面最大,故对dv=35 cm方案进行现场验证。掏槽区域炮孔布置方案如图2(c)所示。炸药为1号岩石乳化炸药,采用耦合装药,孔底起爆;所有炮孔均用炮泥填塞,炮泥长度为0.6 m。两个装药孔共需要1号岩石乳化炸药3.2 kg。
4.2 爆破成腔效果及分析
通过数值模拟计算分析,选择空孔间距dv=35 cm,两装药孔距离40 cm的掏槽方案进行现场试验,具体试验结果见表3,由于篇幅原因本文只对2号试验爆破后形成的槽腔进行分析,如图15所示。由于未进行全断面开挖爆破,实际上形成锥台型槽腔,测量槽腔底部断面宽度与高度并绘制槽腔断面轮廓,以此计算槽腔面积,并与模拟结果对比,如图16。现场试验槽腔断面宽度比模拟结果小4.0%,槽腔断面高度比模拟结果小3.4%,槽腔断面积比模拟结果小4.98%,多次现场试验与模拟结果误差均在5%以内。
表 3 现场试验结果Table 3. Field test results试验 高度/cm 宽度/cm 面积/m2 模拟面积误差/% 1 46.6 63.9 0.2491 2.26 2 42.5 64.7 0.2315 -4.98 3 42.1 66.3 0.2392 −1.8 4 51.5 62.7 0.2554 4.84 5. 结 论
以大红山铜矿地下某巷道为试验场地,并根据围岩力学性质和实际爆破条件进行数值模拟计算对比分析,得到如下结论:
(1) 两等间距掏槽孔同时起爆条件下,掏槽区域中心点的压应力峰值不变,反射拉应力峰值随着大直径空孔间距的增加而减小;
(2) 当空孔间距dv为15 cm、25 cm、35 cm、45 cm、55 cm时,模拟所得槽腔面积分别为0.1641 m2、0.2116 m2、0.2436 m2、0.1740 m2、0.0951 m2,表明存在一个最佳的空孔间距;
(3) 模拟预测的槽腔断面宽度、高度、断面积分别比现场试验结果大4.0%、3.4%和4.98%,多次现场试验与模拟结果误差均在5%以内,说本文预测爆破槽腔断面积的数值模拟方法具有一定可靠性。
因此,可进一步将此数值模拟方法从二维拓展至三维,研究不同掏槽方式下的成腔体积并选择性开展现场试验验证,以期找到预测地下巷道掏槽爆破成腔效果的数值方法。
-
-
[1] Marteau J, Bigerelle M, Mazeran P E, et al. Relation between roughness and processing conditions of AISI 316L stainless steel treated by ultrasonic shot peening[J]. Tribology International, 2015, 82:319-329. doi: 10.1016/j.triboint.2014.07.013 [2] Ganesh B K C, Sha W, Ramanaiah N, et al. Effect of shotpeening on sliding wear and tensile behavior of titanium implant alloys[J]. Materials and Design, 2014, 56(4):480-486. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=d2428b4f6fa62767ff66efd3f8f0b4cd [3] Benedetti M, Fontanari V, Santus C, et al. Notch fatigue behavior of shot peened high-strength aluminium alloys: Experiments and predictions using a critical distance method[J]. International Journal of Fatigue, 2010, 32(10):1600-1611. doi: 10.1016/j.ijfatigue.2010.02.012 [4] 栾伟玲, 涂善东.喷丸表面改性技术的研究进展[J].中国机械工程, 2005, 16(15):1405-1049. doi: 10.3321/j.issn:1004-132X.2005.15.023Luan Weiling, Tu Shandong. Recent trends on surface modification technology of shot peening[J]. China Mechanical Engineering, 2005, 16(15):1405-1049. doi: 10.3321/j.issn:1004-132X.2005.15.023 [5] Al-Obaid Y F. Shot peening mechanics: experimental and theoretical analysis[J]. Mechanics of Materials, 1995, 19(2/3):251-260. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0211961014/ [6] Menig R, Pintschovius L, Schulze V, et al. Depth profiles of macro residual stresses in thin shot peened steel plates determined by X-ray and neutron diffraction[J]. Scripta Materialia, 2001, 45(8):977-983. doi: 10.1016/S1359-6462(01)01063-6 [7] Xing Y M, Lu J. An experimental study of residual stress induced by ultrasonic shot peening[J]. Journal of Materials Processing Technology, 2004, 152(1):56-61. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f4518f2da0effa1755be7dca1bc56be4 [8] 张洪伟, 张以都, 吴琼.喷丸强化过程及冲击效应的数值模拟[J].金属学报, 2010, 46(1):111-117. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CAS201303040000302556Zhang Hongwei, Zhang Yidu, Wu Qiong. Numerical simulations of shot-peening process and impact effect[J]. Acta Metallurgica Sinica, 2010, 46(1):111-117. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CAS201303040000302556 [9] Taehyung K, Hyungyil L, Hong C H, et al. Effects of Rayleigh damping, friction and rate-dependency on 3D residual stress simulation of angled shot peening[J]. Materials and Design, 2013, 46(4):26-37. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1ab77b975ca1cdf02a33903920abb275 [10] Taehyung K, Hyungyil L, Minsoo K, et al. A 3D FE model for evaluation of peening residual stress under angled multi-shot impacts[J]. Surface & Coatings Technology, 2012, 206(19/20):3981-3988. http://www.sciencedirect.com/science/article/pii/S0257897212002599 [11] Sheng X F, Xia Q X, Cheng X Q, et al. Residual stress field induced by shot peening based on random-shots for 7075 aluminum alloy[J]. Transactions of Nonferrous Metals Society of China, 2012, 22:261-267. doi: 10.1016/S1003-6326(12)61717-8 [12] Mylonas G I, Labeas G. Numerical modeling of shot peening process and corresponding produces: Residual stress, surface roughness and cold work prediction[J]. Surface & Coatings Technology, 2011, 205(19):4480-4494. http://www.sciencedirect.com/science/article/pii/S0257897211002696 [13] Watanabe M, Kishimoto S, Xing Y M, et al. Evaluation of strain field around impacted particles by applying electron Moiré method[J]. Journal of Thermal Spray Technology, 2007, 16(5):940-946. doi: 10.1007/s11666-007-9129-1 [14] Schiffner K, Helling C D. Simulation of residual stresses by shot peening[J]. Computers & Structures, 1999, 72(1/2/3):329-340. http://www.sciencedirect.com/science/article/pii/S0045794999000127 [15] Meguid S A, Shagal G, Stranart J C, et al. Three-dimensional dynamic finite element analysis of shot-peening induced residual stresses[J]. Finite Elements in Analysis and Design, 1999, 31(3):179-191. doi: 10.1016/S0168-874X(98)00057-2 [16] Meo M, Vignjevic R. Finite element analysis of residual stress induced by shot peening process[J]. Advances in Engineering Software, 2003, 34(03):569-575. doi: 10.1016-j.clon.2010.02.005/ [17] Umbrello D, Saoubi R M, Outeiro J C. The influence of Johnson-Cook material constants on finite element simulation of machining of AISI 316L steel[J]. International Journal of Machine Tools & Manufacture, 2007, 47(3/4):462-470. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=da2b63f7ce1c3f04be8bf54f29588b1e [18] Wang J M, Liu F H, Yu F, et al. Shot peening simulation based on SPH method[J]. The International Journal of Advanced Manufacturing Technology, 2011, 56(5/6/7/8):571-578. http://d.old.wanfangdata.com.cn/Periodical/sdgydxxb201006013 期刊类型引用(11)
1. 高轩, 尹训江, 王帅帅, 苗龙刚, 郭亚斌. 铁路隧道机械化装药掏槽爆破参数的优化. 爆破器材. 2025(04) 百度学术
2. 国志雨. 双空孔菱形直眼掏槽力学机理及数值模拟研究与应用. 现代矿业. 2025(06) 百度学术
3. 杨成杰,张智宇,孟佳乐,曾晓昌. 基于762 mm直径空孔直眼掏槽爆破参数优化研究. 有色金属(中英文). 2025(01): 116-128 . 百度学术
4. 程兵,叶福,汪泉,郑强强,吕闹,高朋飞. 数值模拟技术在爆破安全课程教学中的应用. 实验室研究与探索. 2025(03): 88-93 . 百度学术
5. 王勉,王建国,马军,李进华,雷露刚,张伟. 小断面巷道精确延时逐孔起爆技术. 兵工学报. 2025(03): 245-256 . 百度学术
6. 徐钦明,陈叶群,黄昕旭,杨赛群,李洪伟. 延期时间对直眼掏槽爆破能量释放的影响规律研究. 煤矿爆破. 2025(01): 19-26 . 百度学术
7. 李曙光. 铁路隧道空孔直眼掏槽爆破数值模拟及现场应用. 地下空间与工程学报. 2025(02): 610-619 . 百度学术
8. 李祥龙,颜世骞,王建国,姚永鑫,黄原明. 分层一次爆破成井精确延时爆破参数研究. 高压物理学报. 2024(02): 175-187 . 百度学术
9. 赵艳伟,孙龙,孙进辉,李清林,陈浩,何丽华. 巷道掘进爆破参数优化设计及施工. 云南冶金. 2024(04): 39-44+58 . 百度学术
10. 陶子豪,李祥龙,胡启文,王建国. 掏槽爆破成腔空孔效应数值模拟研究与分析. 兵工学报. 2024(12): 4246-4258 . 百度学术
11. 刘树新,庄宇,戴谦君,安帅,曹飞. 盲天井中深孔掏槽爆破布孔优化及围岩损伤数值分析. 金属矿山. 2024(12): 44-51 . 百度学术
其他类型引用(4)
-