Analysis on the blast resistance of polymer composite slabs under contact explosions
-
摘要: 高聚物材料具有成型快和膨胀性能好的特点,该材料与碎石和钢筋的复合结构应用于地基处理和城市道路脱空除险加固,具有明显优势。本文中,设计并制作了高聚物碎石板和钢筋高聚物板,开展了接触爆炸冲击下的试验研究,通过毁伤尺寸和所测冲击波数据探讨了2种板的毁伤特性。基于ANSYS/AUTODYN非线性显式有限元程序,建立了试验中毁伤更严重的钢筋高聚物板的接触爆炸全耦合模型,并通过与试验结果的对比,验证了所建耦合模型的准确性和适用性。参数化分析了钢筋高聚物板对炸药量和板厚的敏感性,并利用多参数非线性回归分析方法,提出了钢筋高聚物板迎爆面和背爆面破坏直径的预测公式。结果表明:接触爆炸作用下,高聚物碎石板的毁伤模式以接触部位的局部震塌冲切破坏为主,除此之还有一些毁伤裂纹;钢筋高聚物板的破坏模式主要是迎爆面爆坑毁伤、背爆面剥落损伤和中心冲切贯穿破坏。高聚物碎石板和钢筋高聚物板对爆炸冲击波都具有良好的衰减作用,有望应用到抗爆炸冲击防护领域。Abstract: Polymer materials have the characteristics of fast forming and good expansion performance. Their composite structures with gravel and reinforcement have obvious advantages in foundation treatment and urban road void removal and reinforcement. As reported in this paper, polymer gravel slabs and reinforced polymer slabs were designed and manufactured, and experimental study under contact explosion impact was carried out. The damage characteristics of the two kinds of slabs were investigated through the damage sizes and the measured shock wave data. Based on the ANSYS/AUTODYN nonlinear explicit finite element program, the damage modes and damage diameters of the reinforced polymer slabs were explored, and compared with the experimental results to verify the accuracy and applicability of the established finite element model. The sensitivity of the reinforced polymer slabs to explosive quantity and slab thickness was analyzed parametrically, and the prediction formula of the failure diameter of the top surface and bottom surface of the reinforced polymer slab was proposed by using a multi-parameter regression procedure. The results show that under the action of air contact explosion, the damage mode of the polymer gravel slab is mainly local collapse and perforation at the contact part. Under the impact load of contact explosion, punching and cutting explosion pits appear on the top surface of the slab, tensile failure and collapse area occur on the bottom surface, and a through failure hole is formed in the center of the slab, in addition to some damage cracks. Under the action of air contact explosion, the reinforced polymer slab mainly exhibits crater damage on the top surface, spalling damage on the bottom surface and central punching perforation damage. The reinforced polymer slab has a good attenuation effect on the explosion shock wave. The diffuse reflection effect of the closed bubble in the polymer structure on the shock wave can absorb more energy to alleviate the explosion shock wave, indicating that the polymer has the potential to be applied to the anti-explosion shock protection.
-
Key words:
- contact explosion /
- polymer composite slab /
- damage characteristics /
- dynamic response
-
动能穿甲弹的穿甲能力取决于材料强度、质量以及撞击速度等多种因素。在现代穿甲技术中,钨合金和贫铀合金是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. Test setup
试件 试件编号 试件尺寸 药量/g 高聚物碎石板 P1 50 cm×50 cm×8 cm 15 P2 50 cm×50 cm×8 cm 20 P3 50 cm×50 cm×8 cm 30 钢筋高聚物板 R1 50 cm×50 cm×6 cm 10 R2 50 cm×50 cm×6 cm 15 R3 50 cm×50 cm×6 cm 20 表 2 不同炸药量接触爆炸作用下不同厚度钢筋高聚物板的毁伤系数
Table 2. Damage coefficients of reinforced polymer slabs with different thicknesses subjected to contact explosions with different explosive quantities
w/g K w/g K t=6 cm t=8 cm t=10 cm t=6 cm t=8 cm t=10 cm 1 0.87 0.75 0.75 15 1.00 1.00 1.00 2 1.00 0.94 0.94 20 1.00 1.00 1.00 5 1.00 1.00 1.00 25 1.00 1.00 1.00 10 1.00 1.00 1.00 30 1.00 1.00 1.00 表 3 钢筋高聚物板破坏直径经验公式预测结果与数值模拟结果的比较
Table 3. Results of the damage diameters predicted by the empirical formulas compared with numerically-simulated ones for reinforced polymer slabs
板部位 w/g t/cm 破坏直径/cm 误差/% 数值模拟 经验公式预测 迎爆面 10 6 18.0 16.680 7.33 15 8 19.0 19.204 1.07 20 10 20.5 20.064 2.13 背爆面 10 6 21.0 19.255 8.31 15 8 15.5 16.454 6.15 20 10 14.0 15.113 7.95 -
[1] CHAUDHARI A, GITE V, RAJPUT S, et al. Development of eco-friendly polyurethane coatings based on neem oil polyetheramide [J]. Industrial Crops and Products, 2013, 50: 550–556. DOI: 10.1016/j.indcrop.2013.08.018. [2] CHEN J, YIN X J, WANG H, et al. Evaluation of durability and functional performance of porous polyurethane mixture in porous pavement [J]. Journal of Cleaner Production, 2018, 188: 12–19. DOI: 10.1016/j.jclepro.2018.03.297. [3] WANG F M, SHI M S, LI H J, et al. Experimental study on the anti-permeability properties of polymer grouting materials [J]. Advanced Materials Research, 2011, 284/285/286: 1952–1955. DOI: 10.4028/www.scientific.net/AMR.284-286.1952. [4] FANG H Y, SU Y J, DU X M, et al. Experimental and numerical investigation on repairing effect of polymer grouting for settlement of high-speed railway unballasted track [J]. Applied Sciences, 2019, 9(21): 4496. DOI: 10.3390/app9214496. [5] 刁龙. 高速公路软土路基不均匀沉降处治方法及效果评价 [J]. 交通世界, 2021(30): 106–107. DOI: 10.16248/j.cnki.11-3723/u.2021.30.049.DIAO L. Treatment method and effect evaluation of uneven settlement of expressway soft soil subgrade [J]. Transpoworld, 2021(30): 106–107. DOI: 10.16248/j.cnki.11-3723/u.2021.30.049. [6] 徐建国, 陈志豪, 王壬. 埋地排水管道高聚物注浆修复受力特性分析 [J]. 岩土工程学报, 2021, 43(1): 121–129. DOI: 10.11779/CJGE202101014.XU J G, CHEN Z H, WANG R. Mechanical characteristics of buried drainage pipes repaired by polymer grouting technology [J]. Chinese Journal of Geotechnical Engineering, 2021, 43(1): 121–129. DOI: 10.11779/CJGE202101014. [7] YATES B. Mechanics of cellular plastics [J]. Composites, 1982, 13(4): 362–363. DOI: 10.1016/0010-4361(82)90143-4. [8] 王海福, 冯顺山. 爆炸载荷下聚氨酯泡沫材料中冲击波压力特性 [J]. 爆炸与冲击, 1999, 19(1): 78–83.WANG H F, FENG S S. Properties of shock pressure caused by explosion loads in polyurethane foam [J]. Explosion and Shock Waves, 1999, 19(1): 78–83. [9] 徐全军, 唐丙文, 李裕春. 聚氨酯泡沫材料的隔爆作用 [J]. 解放军理工大学学报(自然科学版), 2007, 8(1): 77–81. DOI: 10.7666/j.issn.1009-3443.20070117.XU Q J, TANG B W, LI Y C. Polyurethane foam material’s gap effects [J]. Journal of PLA University of Science and Technology, 2007, 8(1): 77–81. DOI: 10.7666/j.issn.1009-3443.20070117. [10] LIU S C, ZHAO X H, FANG H Y, et al. Damage characteristics of polymer plates under the impact of the near-field and contact underwater explosion [J]. Advances in Materials Science and Engineering, 2021, 2021: 5957847. DOI: 10.1155/2021/5957847. [11] LIU Z D, ZHAO X H, FANG H Y, et al. Experimental study on the damage characteristics of polymer slabs subjected to air contact and close-in explosions [J]. Advances in Materials Science and Engineering, 2021, 2021: 2825062. DOI: 10.1155/2021/2825062. [12] 张勇. 聚氨酯泡沫铝复合结构抗爆吸能试验及数值模拟分析 [J]. 爆炸与冲击, 2022, 42(4): 045101. DOI: 10.11883/bzycj-2021-0182.ZHANG Y. Testing and numerical simulation of the antiknock energy absorption of polyurethane foam aluminum composite structure [J]. Explosion and Shock Waves, 2022, 42(4): 045101. DOI: 10.11883/bzycj-2021-0182. [13] GARGANO A, DONOUGH M, DAS R, et al. Damage to fibre-polymer laminates caused by surface contact explosive charges [J]. Composites Part B: Engineering, 2020, 197: 108162. DOI: 10.1016/j.compositesb.2020.108162. [14] KELLY M, ARORA H, DEAR J P. The comparison of various foam polymer types in composite sandwich panels subjected to full scale air blast loading [J]. Procedia Engineering, 2014, 88: 48–53. DOI: 10.1016/j.proeng.2014.11.125. [15] 曾祥, 刘彦, 许泽建, 等. 爆炸载荷作用下玻璃钢/硬质聚氨酯泡沫夹层结构抗冲击性能实验研究 [J]. 北京理工大学学报, 2021, 41(11): 1145–1153. DOI: 10.15918/j.tbit1001-0645.2021.036.ZENG X, LIU Y, XU Z J, et al. Experimental study on impact resistance of Glass Fiber Reinforced Plastic/Rigid Polyurethane Foam sandwich structures under air blast loading [J]. Transactions of Beijing Institute of Technology, 2021, 41(11): 1145–1153. DOI: 10.15918/j.tbit1001-0645.2021.036. [16] 邹广平, 孙杭其, 唱忠良, 等. 聚氨酯/钢夹芯结构爆炸载荷下动力学响应的数值模拟 [J]. 爆炸与冲击, 2015, 35(6): 907–912. DOI: 10.11883/1001-1455(2015)06-0907-06.ZOU G P, SUN H Q, CHANG Z L, et al. Numerical simulation on dynamic response of polyurethane/steel sandwich structure under blast loading [J]. Explosion and Shock Waves, 2015, 35(6): 907–912. DOI: 10.11883/1001-1455(2015)06-0907-06. [17] LI M J, DU M R, WANG F M, et al. Study on the mechanical properties of polyurethane (PU) grouting material of different geometric sizes under uniaxial compression [J]. Construction and Building Materials, 2020, 259: 119797. DOI: 10.1016/j.conbuildmat.2020.119797. [18] 刘志远. 高聚物注浆材料工程特性的试验研究 [D]. 郑州: 郑州大学, 2007: 22–30. [19] 石明生. 高聚物注桨材料特性与堤坝定向劈裂注桨机理研究 [D]. 辽宁大连: 大连理工大学, 2011: 22–61. [20] 王娟, 方宏远, 余自森, 等. 高聚物碎石混合料单轴受压性能试验研究 [J]. 建筑材料学报, 2019, 22(2): 320–326. DOI: 10.3969/j.issn.1007-9629.2019.02.024.WANG J, FANG H Y, YU Z S, et al. Experimental study on uniaxial compressive properties of polymer gravel mixtures [J]. Journal of Building Materials, 2019, 22(2): 320–326. DOI: 10.3969/j.issn.1007-9629.2019.02.024. [21] ZHAO X H, WANG G H, LU W B, et al. Damage features of RC slabs subjected to air and underwater contact explosions [J]. Ocean engineering, 2018, 147: 531–545. DOI: 10.1016/j.oceaneng.2017.11.007. [22] WANG Z Y, DU M R, FANG H Y, et al. Influence of different corrosion environments on mechanical properties of a roadbed rehabilitation polyurethane grouting material under uniaxial compression [J]. Construction and Building Materials, 2021, 301: 124092. DOI: 10.1016/J.CONBUILDMAT.2021.124092. [23] 王高辉. 极端荷载作用下混凝土重力坝的动态响应行为和损伤机理[D]. 天津: 天津大学, 2014: 163–167. [24] 廖真, 唐德高, 李治中, 等. 近地面空中爆炸马赫反射数值模拟研究 [J]. 振动与冲击, 2020, 39(5): 164–169,176. DOI: 10.13465/j.cnki.jvs.2020.05.022.LIAO Z, TANG D G, LI Z Z, et al. Numerical simulation for Mach reflection in air explosion near ground [J]. Journal of Vibration and Shock, 2020, 39(5): 164–169,176. DOI: 10.13465/j.cnki.jvs.2020.05.022. [25] 吴红波, 邢化岛, 缪志军, 等. 乳化炸药聚能射流侵彻靶板的数值仿真 [J]. 工程爆破, 2016, 22(1): 68–72. DOI: 10.3969/j.issn.1006-7051.2016.01.015.WU H B, XING H D, MIAO Z J, et al. Numerical simulation of shaped charge jet penetrating on target by emulsion explosive [J]. Engineering Blasting, 2016, 22(1): 68–72. DOI: 10.3969/j.issn.1006-7051.2016.01.015. -