Dynamic response mechanism of a rock-filling interfacial coupling body to blasting in it
-
摘要: 充填采矿法的充填体与矿岩体构成的界面耦合结构体,受采矿爆破影响会持续受到动力扰动,在充-岩界面耦合处易产生脱粘、裂隙扩展等行为,为井下生产带来安全隐患。采用ANSYS/LS-DYNA建立了充-岩界面耦合体模型,分析了爆破作用对界面耦合体结构的力学影响,获取了不同界面粗糙度、充填体养护龄期和起爆方式等因素对爆破裂隙扩展及应力波峰值应力的影响,探讨了爆破动力作用机理。结果表明:(1)爆破冲击在界面耦合体中存在拉、压和剪3种力学作用,且随着界面粗糙度的提高,界面受力呈先上升后下降趋势;(2)随着充填体养护时间增长,界面破坏逐步从受拉转化成剪切损伤;(3)同时起爆对耦合界面的损伤比逐孔起爆的小。Abstract: The interfacial coupling structure between the backfill and ore rock body in the filling mining method will be continuously disturbed by the influence of mining and blasting. In the filling-rock interfacial coupling, it is easy to produce the behaviors of debonding and fissure expansion, which will bring potential safety hazards to underground production. Because the field experiment is time-consuming and laborious, and it is difficult to observe the impact effect and the rock crack propagation process when the explosive is detonated, the simulation method was adopted for research. In the simulation, reasonable simplification is particularly important. According to the actual situation of blasting hole layout, the three rows of blasting holes that were detonated at one time were simplified into a single-row blasting hole model and an edge blasting hole model. According to the research results in related literatures, the coupling surfaces of the filling bodies and the ore rocks were simplified into three kinds (a flat interface, wavy interface and serrated interface). The three different shapes of the interfaces correspond to the different roughnesses of the interfaces, respectively. By considering that the hole arrangement method used in stope blasting is vertical hole arrangement, the holes are parallel, and at the same time, in order to improve the calculation of the software efficiency, simplified the three-dimensional model of the stope into a two-dimensional plane model. After a series of simplifications, a physical model for the filling-rock interface coupling body was proposed, and the corresponding geometric analysis model was established by using the ANSYS/LS-DYNA software And different material parameters were assigned to the different parts of the model, and the blasting effect was analyzed by software calculation. The mechanical influence of the interface coupling body structure was obtained, and the response law of different interface roughness, the curing age of the filling body and the blasting method on the blasting shock was obtained, and the mechanism of the blasting dynamics was discussed. The research results can reveal mechanical behaviors such as debonding and crack propagation at the coupling of filling-rock interface, and clarify the influence of different factors on the law of blasting shock response and the mechanism of blasting dynamics, which has certain guiding significance for downhole safety production. The results show as follows. (1) The blasting impact has three mechanical effects in the interface coupling body: tension, pressure and shear. When the stress wave passes through the interface, the peak acceleration of the monitoring point at the interface will increase due to different degrees of refraction. After passing through the coupling interface, the stress wave energy decays rapidly. (2) The impact of different interface roughness on blasting action is different. The joint roughness coefficient (JRC) represents the roughness of the interface coupling body. With the increase of the JRC value, the interface stress tends to rise first and then decline, but the overall damage decreases. (3) As the curing time of the backfill increases, the fracture range at the coupling interface shrinks, and the interface damage gradually changes from tensile damage to shear damage. (4) The damage of different detonation modes to the interface coupling body is different, and the damage of simultaneous detonation to the coupling interface is weaker than that of hole-by-hole detonation.
-
Key words:
- coupling structure /
- explosive impact /
- dynamic disturbance /
- crack propagation
-
动能穿甲弹的穿甲能力取决于材料强度、质量以及撞击速度等多种因素。在现代穿甲技术中,钨合金和贫铀合金是2种常用的弹芯材料,然而它们在应用过程中均存在不足:钨合金弹体在穿甲过程中因材料具有良好的塑性,弹头易形成类似“蘑菇头”的膨胀结构,这种结构会增大穿甲阻力,进而削弱穿甲效果[1];贫铀合金弹体则常发生“自锐”现象[2],即弹头前端保持尖锐,但贫铀材料伴随的放射性限制了其应用[3]。鉴于上述瓶颈,20余年来,军事领域一直致力于探索对人和环境无放射性危害且具有“自锐”特性的新型动能穿甲弹芯材料,其中钨纤维增强金属玻璃复合材料具有巨大潜力[4]。该复合材料弹体得益于其金属玻璃(metallic glass,MG)基体的高剪切敏感性[5-6],在侵彻/穿甲过程中弹头易呈现出尖锐构型,即“自锐”行为,其侵彻/穿甲能力得到明显提升[7-21]。
一般而言,侵彻/穿甲试验成本高昂,且试验难以捕捉撞击过程中弹靶的变形和破坏细节,也难以测量高速冲击条件下弹体的应力-应变情况。随着计算机技术和材料动态本构模型的发展,有限元模拟已成为侵彻/穿甲试验的良好补充,有助于更深入的机理研究[15-16, 22-24]。Li等[25-26]引入材料微结构自由体积参量,以准确描述金属玻璃基体特有的高强度和高剪切敏感性,进一步发展了热力耦合本构模型,并根据钨纤维增强金属玻璃复合材料的细观结构建立了弹体有限元模型,基于相关穿甲试验,开展了复合材料长杆弹在正侵彻条件下侵彻/穿甲钢靶的有限元模拟,分析了复合材料弹体的“自锐”机理,并系统地探讨了多个关键因素对复合材料弹体“自锐”行为的影响。
在实际工程应用中,动能穿甲弹受发射条件和装甲结构设计等因素的影响,对目标的打击多为斜撞击而非理想的正撞击。在斜撞击过程中,弹体受到靶板的非对称作用力,其弹道特征可能出现偏转甚至跳弹现象[27-28]。然而,对于钨纤维增强金属玻璃复合材料长杆弹斜撞击靶板过程中,弹体在非对称作用力下的变形和破坏特征及其“自锐”特性将可能与正撞击情形截然不同,试验[9]也显示,相较于正撞击,斜撞击条件下复合材料弹体的侵彻/穿甲性能降低,有必要开展进一步研究。
基于Li等[15-16, 22-23]、陈建良等[24]提出的有限元模拟方法,结合相关侵彻/穿甲试验,针对钨纤维增强金属玻璃复合材料和钨合金材料2种材质的长杆弹开展斜侵彻/穿甲钢靶的三维有限元模拟分析,研究2种材质弹体在斜侵彻过程中的弹靶变形和破坏特征,并对比复合材料弹体在正/斜侵彻条件下的“自锐”特征和弹道行为,分析不同撞击速度下撞击倾角对复合材料长杆弹“自锐”行为、侵彻/穿甲性能及相应弹道特征等的影响。
1. 有限元模型
1.1 几何模型
参考相关穿甲试验[9],采用弹道炮发射
∅ 8 mm×88 mm的复合材料弹体和钨合金弹体,正/斜侵彻50 mm厚的30CrMnMo钢靶。复合材料弹体的形貌如图1所示,弹体头部为半球形,弹体中增强钨纤维体积约占80%,直径约为0.3 mm。弹体斜侵彻/穿甲靶板的初始状态如图2所示,其中:D 为弹体直径;L 为弹体长度;v0 为弹体的初始速度;靶板倾角θ 为靶体表面法线方向与弹体初始速度方向的夹角;H 为靶板厚度(试验中,H =50 mm);S 为斜侵彻/穿甲时的靶板有效厚度,S=H/cosθ 。对于复合材料弹体正侵彻/穿甲靶板,二维数值模拟能够描述弹靶变形和破坏的主要特征[15-16, 22-24]。针对斜撞击条件,本文中,将开展三维模拟分析。考虑到弹靶结构和撞击条件的平面对称特征,为减小计算规模,数值模型将采用1/2弹靶结构。复合材料弹体的三维有限元模型如图3(a)所示,其中:MG matrix为金属玻璃基体,WF为钨纤维。弹体和靶板网格均采用八节点六面体单元进行划分,网格尺寸约为75 μm,网格之间以共节点的方式进行固接,弹身横截面的网格划分如图3(b)所示。Zhang等[29]的研究表明,钨纤维与金属玻璃基体之间的界面结合良好,扫描电子显微镜(scanning electron microscope,SEM)直观地显示了实际复合材料的横截面(图3(c))。由图3(a)~(c)可知,所构建的有限元模型能够准确描述复合材料的细观结构特征,验证了模型的准确性。
对钨合金弹体进行数值模拟时,将复合材料弹体模型中的基体材料替换为钨合金材料即可。
1.2 材料本构模型
1.2.1 金属玻璃材料
Li等[25-26]推导的三维修正热力耦合本构模型综合考虑了多种因素对金属玻璃变形和破坏的影响,如应变率、自由体积浓度、温度、应力以及静水应力等,从而能够准确描述金属玻璃材料在不同应变率和初始温度下的力学行为特性。试验弹体中,金属玻璃基体的本构模型参数如表1所示,其中:
˙ε 为应变率,Λc 为压缩条件下的静水应力敏感因子,Λt 为拉伸条件下的静水应力敏感因子。参量 符号 单位 数值 弹性模量 E GPa 96 泊松比 ν 0.36 密度 ρ kg/m3 6125 熔化温度 Tm K 993 玻璃转变温度 Tg K 625 初始温度 T0 K 300 比定容热容 cV J/(kg·K) 400 临界体积 v* m3 2.0×10−29 平均原子体积 Ω m3 2.5×10−29 原子振动频率 f s−1 1×1013 临界破坏自由体积浓度 ξc 0.065 初始自由体积浓度 ξ0 0.05 运动激活能 ΔGm eV ΔGm(˙ε) 几何因子 α 0.05 所需跃迁次数 nD 3 静水应力敏感因子 Λ 0.05(Λc) 0.35(Λt) 1.2.2 金属材料
30CrMnMo钢靶、钨合金弹体以及复合材料弹体中的增强钨纤维等均为常规的晶态合金材料。为了描述高应变率下合金材料的力学行为,综合Johnson-Cook本构模型[30]和累积损伤失效模型[31]来表征金属材料的力学性能,采用Grüneisen状态方程[32]来描述穿甲过程中金属材料的压力状态。钨合金和30CrMnMo钢的Johnson-Cook模型参数如表2所示,其中:
˙ε0 为参考应变率,Tr为参考温度,A为参考应变率条件下的屈服强度,B为应变硬化系数,C为应变率敏感系数,m为温度敏感系数,n为应变硬化指数,D1~D5为材料损伤参数,c0为体积声速,S1为vs-vp曲线斜率的系数(vs为冲击波速度,vp为粒子速度),γ0为Grüneisen系数,a为对γ0的一阶体积修正。其他材料模型参数可参考文献[15-16, 22-23]。表 2 金属材料的Johnson-Cook模型参数Table 2. Johnson-Cook model parameters of metallic materials材料 ρ/(kg·m−3) ν E/GPa ˙ε0/s−1 Tr/K Tm/K cV/(J·kg−1·K−1) 95W钨合金 17900 0.28 410 1 300 1752 134 30CrMnMo钢 7850 0.29 200 1 300 1793 477 材料 A/MPa B/MPa C m n D1 D2 95W钨合金 1650 450 0.016 1.00 0.12 3.00 0 30CrMnMo钢 1200 310 0.014 1.03 0.26 3.20 0 材料 D3 D4 D5 c0/(m·s−1) S1 γ0 a 95W钨合金 0 0 0 3850 1.44 1.58 0 30CrMnMo钢 0 0 0 4578 1.38 1.67 0.47 2. 模型验证及讨论
Li等[15]对文献[9]中的弹体正侵彻进行了模拟分析,本文中,对文献[9]中的2种斜侵彻试验进行数值模拟,试验和模拟结果如表3所示,其中:带“*”号数据是从文献中的图直接测量得到。
倾角θ=50°、撞击速度v0=
1235.1 m/s时,钨合金弹斜侵彻靶板的最终弹靶变形和破坏形貌如图4所示,其中:图4(a)~(c)分别为靶板纵截面、正面和背面的数值模拟结果,图4(d)~(e)分别为实际靶板正面和背面的最终形貌。可以看出,数值模拟结果中,靶板前端的椭圆形开孔和花瓣形撕裂形貌、靶板后端的隆起鼓包形貌等均与试验结果较为相符,侵彻深度和剩余弹长的模拟值与试验值也较一致(表3)。由于钨合金弹体具有良好的塑性,在侵彻过程中其头部发生镦粗变形,导致靶板内形成的弹孔直径相较于原始弹体直径明显增大,特别是在靶板的开坑区域,由于倾斜撞击作用,靶板的开孔口径显著增大,约为弹径的3倍,且靶孔呈现显著的非对称特征,下端孔口相对光滑而上端孔口发生显著的花瓣形撕裂(图4(a)~(b))。另外,受靶板两侧的挤压和摩擦作用,弹体发生质量侵蚀,即弹体碎片自弹头外边缘剥落并向后端流出,碎片的外流也表现出非对称特征,主要沿靶孔上表面流出(图4(a))。由于弹头镦粗变形导致较大的侵彻阻力,钨合金弹体未能穿透靶板,而是在靶板背面形成较大鼓包,其侵彻能力与钨合金弹在θ=0°(正侵彻)、v0=1076.2 m/s条件下的侵彻能力相近(表3);弹体侵蚀也较严重,剩余弹体长度仅为12.8 mm。特别地,由于斜侵彻(θ=50°),弹体在侵彻后期临近穿透靶板时,靶板背面施加于弹体的侵彻阻力并不对称,导致弹道轨迹偏转,同时弹体自身发生明显弯曲,弹头的镦粗变形也呈现出非对称特征(图4(a))。θ=50°、v0=
1263.9 m/s时,复合材料弹体斜侵彻钢靶的最终弹靶变形和破坏形貌如图5所示,结合表3中的侵彻深度和剩余弹长可知,数值模拟与试验结果符合较好。与图4中的钨合金弹斜侵彻相比,弹体在侵彻初期的镦粗变形程度减弱,进而导致靶板前端的开孔口径减小(图5(a)),上端孔口的花瓣形撕裂程度也减弱;在随后的侵彻过程中,复合材料弹体头部发生显著的剪切“自锐”行为,弹头表现为尖头构型,并导致靶板孔径减小,基本与弹径相当;在侵彻后期临近穿透靶板时,受靶板背面的非对称侵彻阻力作用,弹道轨迹偏转,弹头的锐化尖头构型也呈现出非对称特征(图5(a)),但弹体自身未弯曲。由于复合材料弹体的“自锐”行为,其受到的靶板阻力小于钨合金弹工况,因此,弹体侵彻性能显著提高,在相近的撞击速度下直接穿透靶板,同时弹体侵蚀程度也有所减小,剩余长度约为14.5 mm,比钨合金弹工况下的剩余弹长(12.8 mm)长。钨合金弹和复合材料弹在不同撞击倾角和速度下的残余弹体形貌如图6~7所示。可以看出,在正侵彻条件下,2种材质弹体的后端部分均较好地保持了初始形态,未出现明显变形;对于弹体头部,钨合金弹头镦粗为“蘑菇头”形状,临近弹头的弹身直径也明显增大(图6(a)~(b)和(d)),而复合材料弹头在“边缘层”[11, 14]区域锐化为尖头结构,且尖头结构后端区域的弹身与初始形貌基本相同(图7(a)~(b)和(d))。在斜侵彻条件下,因为侵彻后期的弹道轨迹偏转,2种弹体的后端区域均与靶板发生碰撞和摩擦作用,进而产生变形和破坏,钨合金弹体甚至发生了一定程度的弯曲(图6(c));对于弹体头部,钨合金弹头也镦粗为“蘑菇头”形状(图6(c)),复合材料弹头锐化为尖头构型(图7(c)),且二者均呈现出显著的非对称特征,与正侵彻工况明显不同,结合图4(a)和图5(a)可知,钨合金弹和复合材料弹头部的变形主要发生在与靶孔下表面接触的区域,即弹头的非对称变形和破坏形貌主要源于靶板作用的非对称特征。
3. 复合材料弹体斜侵彻/穿甲的“自锐”机理
本节进一步分析斜侵彻/穿甲条件下复合材料弹体的变形和破坏特征,详细讨论其独特的“自锐”行为及作用机制,并与钨合金弹进行对比。
钨合金弹(θ=50°,v0=
1235.1 m/s)斜侵彻50 mm厚的30CrMnMo钢靶时,弹靶内部的塑性应变发展历程如图8所示。可以看出,在斜侵彻的开坑阶段,由于钨合金具有较高的塑性,并且弹体受到靶板施加的非对称侵彻阻力,因此,弹头材料主要向阻力较小的上侧发生侵蚀和流动,初始的半球形弹头发生非对称镦粗变形,其上侧区域的镦粗程度显著大于下侧区域,导致靶板前端的靶孔直径远大于弹径,也呈现出显著的非对称特征(图8(a))。随后,弹头完全侵入靶板,此时弹体受到的靶板阻力较为对称,弹头的镦粗逐渐呈现出相对对称的“蘑菇头”形状,但形貌的对称程度较正侵彻/穿甲工况(图6(b))有所降低,该侵彻过程相对稳定,受靶板两侧的挤压作用,弹头材料从镦粗的“蘑菇头”两侧侵蚀剥落并向弹体后端流失,靶板孔径较初期开坑阶段明显减小(图8(b))。弹体临近靶板后端时,由于弹体前端的靶体材料厚度显著降低,使得靶板的变形模式转变为薄靶的整体响应,靶板背面开始出现鼓包变形;类似地,由于靶板倾斜放置,靶板后端未侵蚀区域的厚度不均匀,后端的鼓包变形不对称,靶体材料施加到弹体的侵彻阻力也不对称,弹道轨迹向侵彻阻力较小的下侧发生偏转,同时弹体发生一定程度的弯曲变形(图8(c)~(d)),弹头的镦粗变形再次呈现出非对称特征(图6(c))。相应地,弹体弯曲变形和弹道偏转进一步增大了弹体的侵彻阻力,且靶板倾斜放置增加了有效厚度,弹体最终未能穿透靶板,弹体侵蚀也更为严重,剩余弹长进一步减小(表3)。复合材料弹(θ=50°,v0=
1263.9 m/s)斜侵彻50 mm厚的30CrMnMo钢靶时弹靶内部的塑性应变发展历程如图9所示。在斜侵彻开坑阶段,同样受到靶板的非对称侵彻阻力,弹头材料也主要向上侧侵蚀和流动,弹头发生非对称的镦粗变形,导致靶板前端形成相对较大的孔径(图9(a)),但弹头的镦粗程度和靶板的孔径均明显小于钨合金弹工况(图8(a))。主要原因在于,高速冲击下,复合材料弹体中的金属玻璃基体具有极高的强度和剪切敏感性[15-16, 22-23, 28-29],在破坏之前可较好地约束钨纤维的变形,而在损伤破坏之后,其局域化剪切带或剪切裂纹对钨纤维产生强烈冲击,导致弹头前端形成一个“边缘层”[15-16],“边缘层”内的钨纤维受靶板的侧向挤压而迅速发生剪切破坏并从弹体剥落,因此,弹头的镦粗程度相对较弱。在随后的弹头完全侵入靶板的过程中,靶板的阻力逐渐对称,弹头的非对称变形和破坏特征减弱,弹头逐渐“自锐”为尖头结构(图9(b)),但形貌的对称程度较正侵彻/穿甲工况(图7(b))有所降低,相应地,侵彻过程中靶板的孔径逐渐与弹径相当。在弹体临近靶板后端时,靶板背面也发生不对称的鼓包变形,导致弹道轨迹向靶板下侧偏转,同时弹体发生弯曲变形,且弹头的“自锐”变形再次呈现非对称特征(图9(c)~(d)),由于弹体内金属玻璃基体的高强度和高硬度,弹体的弯曲程度较钨合金弹工况有所降低(图8(c))。由于弹头锐化导致侵彻阻力较低,弹体的侵彻能力提升,最终穿透靶板(图5(a)),弹体侵蚀程度严重,剩余弹长明显减小(表3),弹体最终的尖头形貌呈现出显著的非对称特征(图7(c))。图10为2种弹体斜侵彻后靶孔最终的轮廓形貌。可以看出:对于2种弹体,靶板的弹道形貌均表现出靶板前端面区域的开坑和临近后端面区域的弹道偏转特性;在整个侵彻过程中,钨合金弹形成的靶板孔径均大于复合材料弹形成的靶板孔径,前者约为后者的1.4倍,钨合金弹体的侵彻深度小于复合材料弹体的侵彻深度。相关差异也主要源于钨合金弹体的头部镦粗变形和复合材料弹体的“自锐”行为。
在弹体侵彻/穿甲过程中,弹体速度是靶板对弹体产生阻力的一个重要参量。图11给出了2种弹体在不同穿甲条件下的速度变化曲线。由图11和表3可知:撞击倾角一致时,侵彻前期,2种弹体的速度变化几乎一致;但在侵彻后期,钨合金弹体的速度下降明显快于复合材料弹体,这主要归因于复合材料弹的头部展现出的“自锐”特性,使得弹体在侵彻时所受阻力显著小于弹头呈镦粗状态的钨合金弹体。另外,相较于正侵彻/穿甲,斜侵彻/穿甲工况下弹体速度下降更快,这主要源于在斜侵彻/穿甲过程中弹头的非均匀变形和弹道弯曲特征。
综上所述,长杆弹斜撞击时,在开坑阶段,受到靶板的非对称阻力,弹头材料向阻力较小的一侧侵蚀和流动,进而导致弹头发生非对称的镦粗变形;在随后相对稳定的侵彻过程中,弹头的非对称变形和破坏特征逐渐减弱,且复合材料弹头开始发生“自锐”行为;在侵彻后期,弹道轨迹向靶板阻力较小的一侧偏转,同时弹体发生弯曲变形,弹头结构的非对称特征再次增强。相较于正撞击,斜侵彻/穿甲工况下,由于受力的非对称性,复合材料弹头的“自锐”特征减弱,易发生非对称变形和破坏,弹体侵蚀和弹道偏转明显加剧。在相同斜撞击工况下,复合材料弹的头部镦粗、弯曲变形和弹道偏转等均小于钨合金弹,其侵彻/穿甲性能仍高于钨合金弹。
4. 撞击倾角对复合材料弹体“自锐”行为和弹道特征的影响
如上所述,复合材料弹在斜侵彻/穿甲过程中将受到非对称的靶板阻力,进而导致与正侵彻/穿甲情形显著不同的弹头“自锐”特征和弹道特性。不同的撞击倾角下,弹体在开坑阶段和后期穿靶阶段受到靶板上下两侧所施加的阻力差异将有所不同,这可能使得弹靶的变形和破坏特征以及弹道偏转产生显著差别,最终导致不同的侵彻/穿甲性能。为了分析撞击倾角的影响,选取v0=857.5 m/s为典型中低速撞击速度,v0=
1263.9 m/s为典型高速撞击速度,模拟撞击倾角θ在0°~70°范围内复合材料弹体撞击钢靶的全过程。θ=30°时,2种典型速度下复合材料弹体撞击钢靶的弹靶变形和破坏历程如图12所示。由图12(a)~(c)可知,高速侵彻时,因基体材料的屈服强度较高,基体材料在破坏之前对钨纤维形成较强约束,使得在开坑阶段,弹头区域材料的侵蚀较小,基体材料破坏后其内部剪切带或剪切裂纹快速切断临近钨纤维,使其向后端流失,因此,弹头的镦粗变形并不明显(图12(a))。在后续侵彻过程中,随着靶板上下两侧的挤压作用逐渐对称,弹头的镦粗变形程度降低,结构的非对称性程度减弱(图12(b))。在后期弹体临近靶板后端时,由于靶板的倾斜放置,靶板后端的隆起鼓包呈非对称性,靶板两侧对弹头的挤压作用也不对称,使得弹头上侧表现为剪切破坏,而下侧钨纤维难以被切断,弹头又呈现出一定程度的非对称镦粗变形(图12(c))。此外,在高速侵彻过程中,侵彻初期弹体头部承受非对称阻力作用的时间较短,并且弹头镦粗程度较低,导致靶板阻力较小,因此,弹道的偏转程度较低,总体来说相对笔直(图12(c))。
如图12(d)~(f)所示,中低速侵彻时,金属玻璃基体对钨纤维的前期约束作用和后期冲击剪切作用减弱,使得弹头明显镦粗,同时弹头承受非对称阻力作用的时间较长,弹头的镦粗变形呈现出显著的非对称特征,且弹道也开始发生偏转(图12(d))。在后续侵彻过程中,弹头的镦粗程度和非对称特征减弱,但弹道持续偏转(图12(e))。在后期,弹头发生“自锐”行为,锐化的尖头构型呈现非对称特征,弹道继续偏转,弹体整体侵蚀程度严重(图12(f))。相较于高速侵彻,中低速侵彻时,复合材料弹体的侵彻能力显著降低。
θ取50°和70°、v0取
1263.9 和857.5 m/s时,复合材料弹体撞击钢靶的最终弹靶变形和破坏形貌如图13所示。由图13(b)可知,θ=50°、v0=857.5 m/s时,弹道严重偏转,弹体严重侵蚀破坏,未形成有效侵彻。θ=70°工况下,弹体的侵彻性能进一步衰减,v0=1263.9 m/s时,弹道严重偏转(图13(c)),v0=857.5 m/s时,弹体直接跳飞(图13(d))。图14给出了复合材料弹体以不同撞击速度和撞击倾角侵彻/穿甲钢靶后靶板孔洞的轮廓形貌,相应的剩余弹体形貌和弹体速度变化历程曲线如图15~16所示。θ=70°时,弹体未能形成有效侵彻甚至发生跳飞,因此,不列出相关结果。由图14(a)可知,v0=
1263.9 m/s时,随着撞击倾角的增大,靶板前端的开孔口径逐渐增大,孔洞的非对称特征也更加明显,表明随着撞击倾角的增大,弹头材料加剧向靶孔上侧侵蚀和流动,弹头的镦粗变形程度增强;在后续侵彻过程中,不同撞击倾角工况下靶板的孔径均逐渐减小,并最终与弹径一致,弹道相对笔直,表明弹体均发生了“自锐”行为;在侵彻后期弹体接近靶板后端面时,弹道轨迹向下偏转逐渐加剧,后端靶孔的非对称特征更加明显。v0=857.5 m/s(图14(b))时,靶板前端的开坑孔径随着倾角的增加而不断增大;在后续侵彻过程中,不同撞击倾角工况下的靶孔均逐渐变小,最终的靶孔底端呈现出尖锐的构型特征,构型的非对称特征随着倾角的增加而增强,弹道偏转也更为显著;当倾角达到50°以上时,弹体已无法形成有效侵彻。从图15(a)可以看出,在v0=
1263.9 m/s工况下,当θ<50°时,在侵彻后期临近穿透靶板时,弹体均受薄靶响应特征影响而发生一定的镦粗变形,且弹头的镦粗程度和构型的非对称特征均随着撞击倾角的增大而增强。当θ=50°时,受较大的靶板倾斜角度和弹体后期较低的穿透靶板速度影响,弹头又呈现出非对称的尖锐头型特征;此外,倾角增加使得靶板的有效厚度增大,随着倾角的增加,弹体侵蚀更为严重,剩余长度逐渐减小,且减小幅度更为明显。由图15(b)可知,v0=857.5 m/s工况下,随着撞击倾角的增大,弹体下侧面的侵蚀和弹体的弯曲变形程度加剧,剩余弹体头形的非对称特征也增强,但剩余弹长基本相同,主要是因为弹道偏转加剧,弹体的侵彻路程差距较小,弹体的侵蚀也基本相同。另外,对比图15(a)和图15(b)可知,在θ<50°范围内,相同撞击倾角工况下,撞击速度越低,弹体的侵蚀程度越严重,与雷波[9]的试验结果相同,这也源于金属玻璃基体强度和剪切敏感性的应变率相关性[15-16, 22-23, 25-26];当θ=50°时,由于v0=857.5 m/s时弹体未能形成有效侵彻,剩余弹长反而相对较长。由图16(a)可知,高速侵彻时,随着撞击倾角的增大,弹体速度的下降幅度增大,但变化特征基本一致。这主要是因为,在靶板倾斜放置产生非对称阻力的条件下,弹头仍能较好地发生“自锐”行为(图15(a)),且弹道相对笔直(图14(a))。由图16(b)可知,中低速侵彻时,相较于正侵彻,θ<40°范围内弹体速度下降更快,这是受弹头“自锐”能力降低、弹道偏转(图14(b))以及弹身与靶板孔壁相互作用(图15(b))等因素的综合影响;当θ≥40°时,随着撞击倾角的增大,弹体速度减速变小,这是因为弹道偏转严重,弹体主要在靶板的前端面附近穿行,受到的靶板总体阻力变小,此时弹体已难以形成有效侵彻。
综上所述,在侵彻初期开坑阶段,不同撞击速度下,随着撞击倾角的增大,由于弹体受力的非对称性,复合材料弹体头部的“自锐”行为更难发生;在侵彻后期弹体临近穿透靶板时,高速侵彻工况下,随着倾角的增大,弹道偏转程度以及弹头变形和破坏的非对称特征加剧,弹体侵蚀更严重。相应地,弹体的侵彻/穿甲性能随倾角增加而逐渐降低,弹体可形成有效撞击的速度阈值也逐渐提高,当倾角达到50°后,弹体在900 m/s速度范围内难以形成有效侵彻。总体来说,撞击倾角对复合材料弹体的侵彻/穿甲主要起到消极影响,实际工程应用中应尽量减小撞击倾角。
5. 结 论
基于细观有限元模拟,研究了钨纤维增强金属玻璃复合材料长杆弹对50 mm厚30CrMnMo钢靶的斜侵彻/穿甲性能,并与钨合金长杆弹进行对比分析,讨论了撞击倾角对复合材料弹体的“自锐”特性、穿甲性能以及相应弹道行为的影响,得到以下主要结论。
(1) 复合材料长杆弹对钢靶的斜侵彻/穿甲行为与正撞击行为完全不同。在弹体斜撞击开坑阶段,受靶板的非对称阻力作用,弹头发生非对称的镦粗变形;在后续侵彻过程中,弹头的非对称变形和破坏特征减弱,并发生“自锐”行为;在侵彻后期,非对称的靶板阻力使弹体弯曲,弹头结构又呈现非对称特征,弹道轨迹偏转。相较于正撞击,斜撞击时,复合材料弹体的“自锐”性能降低,弹体侵蚀和弹道偏转程度加剧,弹体速度加速下降,弹体的侵彻/穿甲性能降低,但仍高于钨合金弹的侵彻/穿甲性能。
(2) 随着撞击倾角的增大,在弹体初期开坑阶段弹头的“自锐”行为更难发生,在侵彻后期弹道的偏转以及弹头的非对称变形和破坏特征加剧,弹体侵蚀更严重,弹体侵彻/穿甲性能降低, 甚至难以形成有效侵彻。
总体来说,撞击倾角不利于复合材料弹体的侵彻/穿甲,在实际工程应用中,应尽量减小撞击倾角。
-
表 1 炮孔布置参数
Table 1. Parameters of blasting hole arrangement
布置方式 炸药密度/(kg·m−3) 孔径/mm 孔深/m 炮孔排距/m 炮孔间距/m 垂直中深孔 1060 90 8 2 2 表 2 耦合界面形态及对应节理粗糙度
Table 2. Coupling interface morphologies and the corresponding joint roughness coefficients
耦合界面类别 剖面线形态 cjr 平直形 0 波浪形 8.12 锯齿形 18.38 表 3 炸药材料及JWL状态方程参数
Table 3. Parameters for explosive materials and JWL equation of state
密度/(kg·m−3) 爆速/(km·s−1) A/GPa B/GPa R1 R2 ω E/GPa 1 060 4 220 0.2 4.5 1.1 0.35 4.2 表 4 岩石和充填体材料参数
Table 4. Parameters for rocks and filling materials
材料 密度/(kg·m−3) 泊松比 弹性模量/GPa 单轴抗压强度/GPa 岩石 2 551 0.25 25.00 100.00 7 d龄期充填体 2 180 0.31 0.92 2.10 28 d龄期充填体 2 200 0.24 2.20 4.17 表 5 不同粗糙度界面耦合体监测点1~4峰值拉应力
Table 5. Peak tensile stress at monitoring points 1−4 in the interface coupling bodies with different roughnesses
监测点编号 单排炮孔峰值拉应力/MPa 两帮炮孔峰值拉应力/MPa cjr=0 cjr=8 cjr=20 cjr=0 cjr=8 cjr=20 1 9.96 3.06 2.56 0.73 4.76 0.73 2 0 2.85 0.04 0.90 4.46 0.04 3 0 9.17×10−3 0.02 3.73×10−3 0.01 0.02 4 0 9.61×10−3 0.01 4.86×10−3 0.02 0.02 表 6 不同龄期界面耦合体监测点1~4峰值拉应力
Table 6. Peak tensile stress at monitoring points 1−4 in different-age interface coupling bodies
监测点编号 单排炮孔峰值拉应力/MPa 两帮炮孔峰值拉应力/MPa 7 d龄期 28 d龄期 7 d龄期 28 d龄期 1 9.96 3.12 0.73 0.70 2 1.15×10-3 2.10×10-3 0.90 0.97 3 1.43×10-3 0.03 3.73×10-3 0.56 4 9.56×10-4 0.02 4.86×10-3 0.26 表 7 不同起爆方式下界面耦合体监测点1~4峰值拉应力
Table 7. Peak tensile stress at monitoring points 1−4 in interfacial coupling bodies with different detonation modes
监测点编号 单排炮孔峰值拉应力/MPa 两帮炮孔峰值拉应力/MPa 同时起爆 逐孔起爆 同时起爆 逐孔起爆 1 9.96 2.66 0.73 14.64 2 0 0.09 0.90 2.34 3 0 0.05 3.73×10−3 0.19 4 0 0.04 4.86×10−3 0.16 -
[1] 张海波, 宋卫东. 评述国内外充填采矿技术发展现状 [J]. 中国矿业, 2009, 18(12): 59–62. DOI: 10.3969/j.issn.1004-4051.2009.12.018.ZHANG H B, SONG W D. Discussion on the current state of backfill mining from the domestic and foreign development [J]. China Mining Magazine, 2009, 18(12): 59–62. DOI: 10.3969/j.issn.1004-4051.2009.12.018. [2] 王湘桂, 唐开元. 矿山充填采矿法综述 [J]. 矿业快报, 2008, 24(12): 1–5.WANG X G, TANG K Y. Overview of cut and fill method for mines [J]. Express Information of Mining Industry, 2008, 24(12): 1–5. [3] 尹升华, 吴爱祥. 缓倾斜中厚矿体采矿方法现状及发展趋势 [J]. 金属矿山, 2007(12): 10–13. DOI: 10.3321/j.issn:1001-1250.2007.12.002.YIN S H, WU A X. Status quo of mining methods for gently inclined medium-thick orebodies and their development trend [J]. Metal Mine, 2007(12): 10–13. DOI: 10.3321/j.issn:1001-1250.2007.12.002. [4] 刘光生. 充填体与围岩接触成拱作用机理及强度模型研究[D]. 北京: 北京科技大学, 2017: 17−22.LIU G S. Required strength model of cemented backfill with research on arching mechanism considering backfill-rock interaction [D]. Beijing: University of Science and Technology Beijing, 2017: 17−22. [5] 韩斌, 王贤来, 肖卫国. 基于多元非线性回归的井下采场充填体强度预测及评价 [J]. 采矿与安全工程学报, 2012, 29(5): 714–718.HAN B, WANG X L, XIAO W G. Estimation and evaluation of backfill strength in underground stope based on multivariate nonlinear regression analysis [J]. Journal of Mining and Safety Engineering, 2012, 29(5): 714–718. [6] 贺桂成, 刘永, 丁德馨, 等. 废石胶结充填体强度特性及其应用研究 [J]. 采矿与安全工程学报, 2013, 30(1): 74–79.HE G C, LIU Y, DING D X, et al. Strength characteristic of cemented waste rock backfills and its application [J]. Journal of Mining and Safety Engineering, 2013, 30(1): 74–79. [7] 孙宁新, 雷明锋, 张运良, 等. 软弱夹层对爆炸应力波传播过程的影响研究 [J]. 振动与冲击, 2020, 39(16): 112–119, 147. DOI: 10.13465/j.cnki.jvs.2020.16.016.SUN N X, LEI M F, ZHANG Y L, et al. A study on the influence of weak interlayer on the propagation process of explosion stress wave [J]. Journal of Vibration and Shock, 2020, 39(16): 112–119, 147. DOI: 10.13465/j.cnki.jvs.2020.16.016. [8] 李夕兵. 论岩体软弱结构面对应力波传播的影响 [J]. 爆炸与冲击, 1993, 13(4): 334–342.LI X B. Influence of the structural weakness planes in rock mass on the propagation of stress waves [J]. Explosion and Shock Waves, 1993, 13(4): 334–342. [9] 李夕兵, 陈寿如. 应力波在层状矿岩结构中传播的新算法 [J]. 中南大学学报(自然科学版), 1993, 24(6): 738–742.LI X B, CHEN S R. A new method for investigating the propagation of stress waves through layered rockmass [J]. Journal of Central South University, 1993, 24(6): 738–742. [10] 李夕兵, 古德生, 赖海辉. 爆炸应力波遇夹层后的能量传递效果 [J]. 有色金属, 1993, 45(4): 1–6.LI X B, GU D S, LAI H H. Energy transmission effects of transient stress waves through sandwiches in rock mass [J]. Nonferrous Metals, 1993, 45(4): 1–6. [11] SCHOENBERG M. Elastic wave behavior across linear slip interfaces [J]. The Journal of the Acoustical Society of America, 1980, 68(5): 1516–1521. DOI: 10.1121/1.385077. [12] PYRAK L J, MYER L R, COOK N G W. Anisotropy in seismic velocities and amplitudes from multiple parallel fractures [J]. Journal of Geophysical Research: Solid Earth, 1990, 95(B7): 11345–11358. DOI: 10.1029/JB095iB07p11345. [13] SCHOENBERG M E, MUIR F. A calculus for finely layered anisotropic media [J]. Geophysics, 1989, 54(5): 581–589. DOI: 10.1190/1.1442685. [14] FAN L F, MA G W, LI J C. Nonlinear viscoelastic medium equivalence for stress wave propagation in a jointed rock mass [J]. International Journal of Rock Mechanics and Mining Sciences, 2012, 50: 11–18. DOI: 10.1016/j.ijrmms.2011.12.008. [15] LI J C, MA G W, ZHAO J. An equivalent viscoelastic model for rock mass with parallel joints [J]. Journal of Geophysical Research: Space Physics, 2010, 115(B3): B03305. DOI: 10.1029/2008jb006241. [16] MA G W, FAN L F, LI J C. Evaluation of equivalent medium methods for stress wave propagation in jointed rock mass [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2013, 37(7): 701–715. DOI: 10.1002/nag.1118. [17] 杨立云, 刘振坤, 周莹莹, 等. 爆炸应力波在含层理介质中传播规律的实验研究 [J]. 爆破, 2008, 35(2): 1–5, 11. DOI: 10.3963/j.issn.1001-487X.2018.02.001.YANG L Y, LIU Z K, ZHOU Y Y, et al. Study on propagation law of explosive stress wave in layered media [J]. Blasting, 2008, 35(2): 1–5, 11. DOI: 10.3963/j.issn.1001-487X.2018.02.001. [18] 杨仁树, 李炜煜, 方士正, 等. 层状复合岩体冲击动力学特性试验研究 [J]. 岩石力学与工程学报, 2019, 38(9): 1747–1757. DOI: 10.13722/j.cnki.jrme.2019.0021.YANG R S, LI W Y, FANG S Z, et al. Experimental study on impact dynamic characteristics of layered composite rocks [J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(9): 1747–1757. DOI: 10.13722/j.cnki.jrme.2019.0021. [19] 包会云. 层状岩体应力波传播特性及震源定位研究[D]. 沈阳: 东北大学, 2015: 45−61.BAO H Y. Research on propagation characteristics of stress wave propagation in layered rock and vibration source localization [D]. Shenyang: Northeastern University, 2015: 45−61. [20] 武仁杰, 李海波. SHPB冲击作用下层状千枚岩多尺度破坏机理研究 [J]. 爆炸与冲击, 2019, 39(8): 083106. DOI: 10.11883/bzycj-2019-0187.WU R J, LI H B. Multi-scale failure mechanism analysis of layered phyllite subject to impact loading [J]. Explosion and Shock Waves, 2019, 39(8): 083106. DOI: 10.11883/bzycj-2019-0187. [21] 刘婷婷, 李新平, 李海波, 等. 应力波在充填节理岩体中传播规律的数值研究 [J]. 岩石力学与工程学报, 2016, 35(S2): 3552–3560. DOI: 10.13722/j.cnki.jrme.2016.0829.LIU T T, LI X P, LI H B, et al. Numerical study on stress wave propagation across filled joints [J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(S2): 3552–3560. DOI: 10.13722/j.cnki.jrme.2016.0829. [22] 李秀虎, 郭连军, 潘博, 等. 节理岩体爆破数值模拟研究 [J]. 辽宁科技大学学报, 2017, 40(5): 390–395. DOI: 10.13988/j.ustl.2017.05.013.LI X H, GUO L J, PAN B, et al. Numerical simulation of blasting in jointed rock mass [J]. Journal of University of Science and Technology Liaoning, 2017, 40(5): 390–395. DOI: 10.13988/j.ustl.2017.05.013. [23] 朱鹏瑞, 宋卫东, 曹帅, 等. 爆破动载下胶结充填体的张拉力学响应机制 [J]. 采矿与安全工程学报, 2018, 35(3): 605–611. DOI: 10.13545/j.cnki.jmse.2018.03.022.ZHU P R, SONG W D, CAO S, et al. Tensile mechanical response mechanism of cemented backfills under blasting load [J]. Journal of Mining and Safety Engineering, 2018, 35(3): 605–611. DOI: 10.13545/j.cnki.jmse.2018.03.022. [24] 刘志祥, 李夕兵. 爆破动载下高阶段充填体稳定性研究 [J]. 矿冶工程, 2004, 24(3): 21–24. DOI: 10.3969/j.issn.0253-6099.2004.03.007.LIU Z X, LI X B. Research on stability of high-level backfill in blasting [J]. Mining and Metallurgical Engineering, 2004, 24(3): 21–24. DOI: 10.3969/j.issn.0253-6099.2004.03.007. [25] 王振昌. 二步回采围岩和充填体爆破损伤特征及安全阈值研究[D]. 福州: 福州大学, 2018: 10−26.WANG Z C. Study on blasting damage characteristics and safety threshold of surrounding rock and backfill in two-step mining [D]. Fuzhou: Fuzhou University, 2018: 10−26. [26] 李化, 黄润秋. 岩石结构面粗糙度系数JRC定量确定方法研究 [J]. 岩石力学与工程学报, 2014, 33(S2): 3489–3497. DOI: 10.13722/j.cnki.jrme.2014.s2.012.LI H, HUANG R Q. Method of quantitative determination of joint roughness coefficient [J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(S2): 3489–3497. DOI: 10.13722/j.cnki.jrme.2014.s2.012. [27] 王岐. 用伸长率R确定岩石节理粗糙度系数的研究[C]//地下工程经验交流会论文选集. 甘肃: 中国岩石力学与工程学会, 1982: 355−360.WANG Q. Study on determining the roughness coefficient of rock joints with elongation R [C] // Gansu: Meeting on Underground Engineering Experience Sharing, 1982: 355−360. [28] 杜时贵, 陈禹, 樊良本. JRC修正直边法的数学表达 [J]. 工程地质学报, 1996, 4(2): 36–43.DU S G, CHEN Y, FAN L B. Mathematical expression of JRC modified straight edge [J]. Journal of Engineering Geology, 1996, 4(2): 36–43. [29] 段军, 常诗源, 张鹏飞, 等. 双利铁矿台阶爆破空气间隔长度研究与应用 [J]. 中国矿业, 2020, 29(6): 155–158. DOI: 10.12075/j.issn.1004-4051.2020.06.006.DUAN J, CHANG S Y, ZHANG P F, et al. Study and application of bench blasting interval length in Shuangli mine [J]. China Mining Magazine, 2020, 29(6): 155–158. DOI: 10.12075/j.issn.1004-4051.2020.06.006. [30] 曹杨. 中深孔台阶爆破模拟与研究[D]. 太原: 中北大学, 2012: 27−47.CAO Y. Simulation and research on the medium-length hole bench blasting [D]. Taiyuan: North University of China, 2012: 27−47. [31] 季京晨. 花岗岩物理力学性质与宏微观力学特性研究[D]. 安徽淮南: 安徽理工大学, 2019: 27−28.JI J C. Study on the physical and mechanical properties of the granite and the micro-mechanical properties of the macro [D]. Huainan, Anhui: Anhui University of Science and Technology, 2019: 27−28. [32] 刘永涛. 不同龄期尾砂胶结充填体单轴压缩破坏的声发射及断口分形特征[D]. 武汉: 武汉科技大学, 2018: 21−26.LIU Y T. Acoustic emission and fractal characteristics of cemented tailings backfill in different age under uniaxial compression [D]. Wuhan: Wuhan University of Science and Technology, 2018: 21−26. [33] 蔡美峰, 何满潮, 刘东燕. 岩石力学与工程[M]. 2版. 北京: 科学出版社, 2013: 24−25.CAI M F, HE M C, LIU D Y. Rock mechanics and engineering [M]. 2nd ed. Beijing: Science Press, 2013: 24−25. -