Numerical simulation and experimental study on the damage of water partitioned structure by a shaped charge warhead with a combined charge liner
-
摘要: 为了研究组合药型罩聚能装药战斗部对含水复合结构的毁伤机理,基于LS-DYNA软件的任意拉格朗日-欧拉(arbitrary Lagrangian-Eulerian, ALE)流固耦合算法,对水下组合药型罩聚能装药战斗部侵彻体的形成以及穿靶过程开展研究,采用数值模拟等比例模型对水下组合药型罩聚能装药战斗部对靶板毁伤进行试验验证。研究结果表明,在偏心亚半球缺罩罩顶设计偏心亚半球形罩能够在侵彻体前端形成细长的杆式射流,可以增加整个侵彻体长度和头部侵彻体速度。在穿水和靶板过程中,利用头部杆式射流形成空腔帮助后续侵彻体低阻随进。对靶板毁伤过程的分析发现,与战斗部直接连接的第1层靶板将会受到侵彻体的高速冲击作用和爆炸波沿水介质传播过来的强冲击波联合作用,而随着水层厚度的增加,沿水中传播的爆炸冲击波强度会被迅速衰减,爆炸冲击波对后续靶板的作用变得不明显,主要为侵彻体的冲击作用。最后利用设计的组合药型罩结构开展了试验验证,对比分析了每层靶板的穿孔尺寸,试验结果与数值计算结果符合较好,最大误差小于15%。Abstract: In order to study the damage mechanism of the shaped charge warhead with a combined charge liner to the water containing composite structure, the formation and penetration process of the penetrator formed by the combined charge liner were studied based on the arbitrary Lagrangian-Euler (ALE) fluid structure coupling algorithm in the LS-DYNA. The damage of the shaped charge warhead with composite liner to the target was verified by experiments. A hemispherical liner eccentric to the axis was designed at the top of the original eccentric sub-hemispherical liner. The forming process of the penetrator, the response state of the water medium, the dynamic energy loss in the process of penetrating the target and the damage mechanism to the target were analyzed for the warhead with the combined liner. The results show that the design of the sub-hemispherical liner on the top of the eccentric sub-hemispherical liner can form a slender rod-like jet at the front of the penetrator, which can increase the whole length of the penetrator and the velocity of the head penetrator. In the process of the target, the head rod-like penetrators form a cavity to help the subsequent penetrators follow with low resistance. Through the analysis of the damage process to the target, it is found that the first layer of target directly connected with the warhead will be affected by both the high-speed impact of the penetrator and the strong shock wave transmitted by the explosion wave along the water medium. With the increase of the thickness of the water layer, the intensity of the explosion shock wave propagating along the water will be rapidly attenuated, and the effect of the explosion shock wave becomes less obvious to the subsequent target. The experimental verification was carried out by the warhead with composite liner structure. The perforation size of each target was compared and analyzed. The experimental results are in good agreement with the numerical simulation results, and the maximum error is within 15%.
-
随着武器装备防护能力的提高,对战斗部毁伤能力的要求进一步提升。聚能装药是为提高武器系统的定向威力而广泛采用的一种方式[1],其结构主要由主装药、药型罩、外壳以及波形发生器等组成。通过起爆主装药产生爆轰波压垮药型罩形成高速金属射流,金属射流的速度可达到9000 m/s甚至更高。药型罩是聚能装药战斗部的核心部件,也是影响其毁伤威力的关键因素之一。通过设计不同的药型罩结构,可以获得不同形态和速度的侵彻体,如爆炸成形弹丸(速度1.5~3 km)[2],杆式射流和射流(通常超过6 km/s,最高可达12 km/s)[3-5]。
水介质对侵彻体能量具有明显的衰减作用,因此提高射流或射弹在水中的毁伤威力是提高水下聚能战斗部毁伤效应的关键。王海福等[6]研究表明,杆式侵彻体径向尺寸适中,而且受拉伸会断裂形成多段侵彻体;前驱侵彻体可在水中形成较理想的空腔,在水中速度衰减最慢;偏心亚半球形罩可有效形成空腔随进效应。周方毅[7]研究表明,球锥组合药型罩形成的侵彻体能够形成兼具锥形罩(金属射流)和球缺罩(翻转弹丸)两种药型罩的侵彻体特点,使成型侵彻体质量增大、速度加快,形成高速杆式射流,非常有利于破坏目标。在此研究基础上,周方毅等[8]又提出一种双球缺组合药型罩,对于防护能力强的目标,双球缺组合药型罩比圆锥、球缺型组合药型罩破坏效果更好。张春辉等[9]设计了一种球缺罩与偏心亚半球形罩组合聚能战斗部,并通过改变小药型罩材质,使复合材质的组合药型罩与传统杆式射流相比侵彻相同厚度的水介质与靶板后,剩余动能提升了28.59%。李兵等[10]结合数值模拟和试验研究了半球形药型罩对双层圆柱壳结构的毁伤过程,结果表明水中聚能战斗部对双层圆柱壳结构的破坏载荷主要有金属射流、冲击波载荷及气泡载荷3种,金属射流穿透力强,造成结构局部小尺寸破口,而冲击波载荷及气泡载荷作用面积大,引起结构大面积破口及塑性凹陷。王玉等[11]对比分析了球缺罩聚能装药战斗部对水下单层潜艇靶板试验和数值模拟结果,发现射流头部距离和时间关系曲线、靶板变形扰度曲线、射流头部速度数值模拟计算与试验误差均在15%以内,在此基础上设计了圆锥+球缺组合药型罩和双球缺组合药型罩并展开数值模拟研究,研究结果表明双球缺组合药型罩质量分布均匀,射流形态好,有效射流转化率高。王长利等[12]通过实验以及数值模拟,开展了防雷舱结构在聚能装药水下爆炸作用下的毁伤研究,探讨冲击波在多介质结构中的传播规律及其对结构的毁伤机理,药型罩采用了变壁厚紫铜材料。
目前虽然组合药型罩聚能装药战斗部对含水复合结构的毁伤已有报道[6-15],但组合药型罩聚能侵彻体的形成、侵彻体对水下多层靶标毁伤机理以及侵彻体与水介质的作用规律仍需进一步研究,且组合药型罩聚能装药战斗部对于含水复合结构毁伤试验研究的报道相对较少,因此本文拟设计一种组合药型罩聚能装药战斗部,结合数值模拟和试验对含水复合结构多层靶的毁伤性能进行研究。
1. 数值模型
1.1 结构设计
根据报道[8],射流在水中由于过于拉伸,会造成射流分段且不连续,随进射流也会受到水干扰,而采用爆炸成形弹丸由于径向尺寸较大,与水接触面积较大能量衰减相比射流和杆式射流更快。因此,为了提高侵彻体的二次毁伤威力,采用偏心亚半球形罩和偏心亚半球缺罩形成组合药型罩。偏心亚半球形罩形成杆式射流,在侵彻体最前端侵彻水介质和靶板,为后续侵彻体提供瞬时空腔,偏心亚半球缺罩形成的侵彻体在空腔中低阻随进。
根据以上思路,设计了组合药型罩聚能装药战斗部,战斗部主要由3部分组成,分别是壳体、主装药和药型罩(见图1(a))。组合药型罩主要由大小药型罩构成,均采用紫铜材料,通过机加方式可形成连续结构。大药型罩采用偏心亚半球缺罩,小药型罩采用偏心亚半球形罩。主装药采用8701炸药,外壳采用铝材料。其中,小药型罩底部开口直径d1为20 mm,壁厚α为3.6 mm(即0.18d1),偏心距δ为6 mm(即0.3d1),炸高h1为175.4 mm(即8.77d1);大药型罩的装药直径d2为122 mm,炸高h2为147 mm(即1.2d2),外径R为68 mm(即0.56d2),内径r为76 mm(即0.62d2),外径与内径偏心距为17.5mm(即0.14d2);装药长径比l/d2为1.13。
另外,为研究增加的小药型罩(即偏心亚半球形罩)对侵彻体形成过程以及侵彻体与水介质和靶板作用过程的影响,设计了只含大药型罩(即偏心亚半球缺罩,见图1(b))的结构并进行了数值模拟。
1.2 数值模拟计算
1.2.1 计算模型
组合药型罩聚能战斗部对含水复合结构毁伤的数值模拟模型见图2,主要由战斗部、靶板、后效靶、水和空气组成。计算模型含有四层靶板,两层靶板的间距分别为172、119、172 mm,在第四层靶板后放置一层后效靶,与第四层靶板间隔100 mm。从第一层到后效靶的靶板厚度分别为28、12、12、28、14 mm。四层靶板置于水中,后效靶板置于空气中。战斗部直接与靶板贴合,通过主装药顶部起爆,靶板均采用45钢。采用LS-DYNA软件,计算模型采用二维轴对称结构,采用ALE流固耦合算法,网格尺寸0.1 cm,单位制为cm-g-μs。其中壳体和靶板采用拉格朗日网格,主装药和药型罩采用欧拉网格,采用*INITIAL_VOLUME_FRACTION_GEOMETRY关键字对欧拉网格进行材料填充,材料模型及本构模型使用见表1。偏心亚半球缺罩毁伤过程的数值模拟与组合药型罩相同,区别只是药型罩的结构。
表 1 材料本构模型及状态方程Table 1. Constitutive model and state equation of the materials材料 本构模型 状态方程 炸药 HIGH_EXPLOSIVE_BURN JWL 药型罩 STEINBERG GRÜNEISEN 壳体 JOHNSON_COOK GRÜNEISEN 水 NULL LINEAR_POLYNOMIAL 空气 NULL LINEAR_POLYNOMIAL 靶板 JOHNSON_COOK GRÜNEISEN 1.2.2 材料参数选取
炸药、药型罩、壳体、水、靶板采用的本构模型及状态方程见表1。药型罩采用STEINBERG本构模型,STEINBERG本构模型可用于模拟极高应变率下材料的变形,其参数选取参考文献[16]。
主装药采用8701炸药,其密度ρ=1.72 g/cm3,爆速为8425 m/s, pCJ=30.4 GPa,状态方程采用JWL方程:
p=AJWL|1−ωR1V|e−R1V+BJWL|1−ωR2V|e−R2V+ωE0V (1) 式中:p为冲击波压力;V为爆轰产物与初始炸药的体积比;AJWL、BJWL、R1、R2和
ω 为材料参数,E0为比内能,所有材料参数选取参考文献[17]。空气和水采用LINEAR_POLYNOMIAL描述流体的状态方程:
p=C0+C1u+C2u2+C3u3+(C4+C5u+C6u2)E (2) 式中:C0、C1、C2、C3、C4、C5、C6为常数,E为单位体积的内能,u为相对体积,其参数选取参考文献[18-19]。
靶板和壳体均采用JOHNSON_COOK本构模型描述高应变、应变率和高温下的材料本构方程,本构方程描述如下:
σ=(AJC+BJCεn)(1+CJCln˙ε)(1−|T−TrTm−Tr|m) (3) 式中:AJC、BJC、CCJ、n和m为材料参数;
ε 为塑性应变,˙ε 为塑性应变率;Tr为参考温度,Tm为熔点,其参数选取参考文献[20-21]。2. 组合药型罩聚能侵彻体穿靶过程数值模拟计算结果及分析
2.1 侵彻体形成过程
基于数值模拟计算结果,对组合药型罩形成侵彻体以及穿靶过程进行了分析。首先对比分析了偏心亚半球缺罩和组合药型罩形成聚能侵彻体并侵彻靶板过程。图3为侵彻第1层靶板时在爆轰波压垮作用下药型罩的侵彻体形成过程。分析发现,在20 μs时,组合药型罩小药型罩的内壁轴线位置附近最先获得加速(图3(a)红色区域),在50 μs时头部侵彻体发生了断裂,此时侵彻体头部到达第1层靶板位置(图3中未画第1层靶板),大药型罩此时并未完全形成稳定的侵彻体构型。在80 μs时,偏心亚半球缺罩形成了比较稳定的侵彻体构型并开始侵彻第1层靶板(图3(b))。此时组合药型罩中大药型罩的侵彻体构型也基本成型。
图4和图5分别为80 μs时偏心亚半球缺罩和组合药型罩形成的侵彻体外形及速度分布以及在轴线位置侵彻体速度随位置分布。分析侵彻体形态(见表2)发现,组合药型罩与偏心亚半球缺罩的侵彻体主要区别为前端和尾端(图5中的S1段和S3段),组合药型罩由于有小药型罩形成的杆式侵彻体在前端,因此整体更长(L1=1.64L2)。组合药型罩的S3段侵彻体头部已经到达靶板位置并处于侵彻靶板过程,若无靶板,组合药型罩形成的侵彻体会更长。由于小药型罩作用,组合药型罩侵彻体与偏心亚半球缺罩头部形状完全不同。通过对比偏心亚半球缺罩与组合药型罩外形尺寸以及轴线处的速度分布发现,组合药型罩S2段与偏心亚半球缺罩(L2段)外形尺寸非常接近,轴线处的速度分布也非常接近。组合药型罩中小药型罩产生了更多速度相对更高的侵彻体分布在侵彻体头部(S3段,80 μs时头部最高速度3267 m/s,在50 μs时最高速度达到4570 m/s,偏心亚半球缺罩80 μs时头部最高速度3005 m/s),同时也产生了一部分速度相比偏心亚半球缺罩速度更低的部分分布在尾部(S1段,80 μs最低速度为674 m/s,组合药型罩杵体最低速度1083 m/s)。组合药型罩头部出现速度急剧减小的情况(A点以后),主要原因为在80 μs时侵彻体头部与第1层靶板发生作用,由于质量较小而导致侵彻过程头部侵彻体速度迅速减小。
表 2 80 μs时组合药型罩和偏心亚半球缺罩的外形尺寸统计Table 2. Statistics of overall dimensions of combined liner and eccentric sub-hemispherical liner at 80 μsD1/mm D2/mm L1/mm S1/mm S2/mm 44 (0.36d2) 12 (0.10d2) 169 (1.39d2) 22 (0.18d2) 104 (0.85d2) S3/mm D3/mm D4/mm D5/mm L2/mm 43 18 46 (0.38d2) 12 (0.10d2) 103 (0.84d2) 通过以上对比分析发现,在偏心亚半球缺药型罩基础上增加的小药型罩能够形成速度相对更高的细长杆式射流分布于侵彻体前端,使整个侵彻体更长,而增加的小药型罩并不会明显改变原有侵彻体的速度分布和侵彻体构型,且可以使侵彻体头部速度更高。因此通过设计不同的小药型罩结构,可以在不显著改变原有侵彻体构型和速度分布情况下,增加整个侵彻体的长度和头部速度并改变头部侵彻体形态。
2.2 侵彻体穿过靶板过程
2.2.1 水介质响应状态
设计组合药型罩的主要目的是希望在侵彻水介质过程中,小药型罩形成侵彻体能够在前端为后续侵彻体开辟空腔通路,通过对侵彻体形成过程分析发现,在偏心亚半球缺药型罩顶部设计小药型罩能够实现。尽管如此,随后的侵彻靶板和水介质会对形成的分段侵彻体形成较大干扰,影响其穿靶过程。因此有必要对侵彻体穿靶过程水介质的响应状态进行分析。
图6为偏心亚半球缺罩和组合药型罩穿靶过程水介质的响应状态,其中Da和Db为射流的头部直径。分析结果发现组合药型罩和偏心亚半球缺罩在穿水过程中均能形成空腔随进效应。偏心亚半球缺罩在穿透第1层靶板后,侵彻体头部达到第2层靶板时,形成明显的空腔使后续侵彻体低阻随进,但侵彻体并没有断裂成多个部分,在215 μs时形成的空腔尺寸l2与侵彻体杵体的径向尺寸l1之比为2.25(l2=90 mm,l1=40 mm),即最大空腔径向尺寸为侵彻体最大径向尺寸的2.25倍,因此头部形成的较大空腔能够保证后续侵彻体在空腔中运动。在侵彻完成第2层靶板并到达第3层靶板时(345 μs),仍可观察到空腔随进效应,但形成的最大空腔径向尺寸变小(345 μs时头部随进侵彻体所处的径向最大空腔尺寸为56 mm,侵彻体径向尺寸为50 mm),与后续侵彻体尺寸接近。在侵彻体第四层靶板时(550 μs),已不能观察到空腔随进现象,整个侵彻体已汇聚到一起。在穿透第4层靶板后(810 μs),侵彻体头部径向尺寸超过轴向(长径比0.54),但最终侵彻体并未能穿透后效靶板(1000 μs)。对组合药型罩的穿靶过程分析发现,在760 μs时侵彻体已经完全穿透后效靶板。在155 μs时头部的侵彻体已经到达第2层靶板处,虽然侵彻体头部尺寸相对偏心亚半球缺罩头部尺寸更小,形成的空腔也相对更小(l3/l4=2),但能够观察到组合药型罩头部后面的侵彻体均处于空腔中。在400 μs时侵彻体由于速度梯度较大出现了分段(图6(b)中P1和P2为断裂点),分析断裂P1可以发现,头部侵彻体为后续侵彻体提供了非常理想的空腔环境,并在穿透完水介质后,头部侵彻体仍然可以为尾部侵彻体提前进行开孔。
王海福等[1]认为,水介质的扩展速度不仅反映空腔径向尺寸效应,同时也表明侵彻体水中侵彻过程消耗能量的多少。因此本文将对比分析在不同水层(每个水层的中间位置,距离轴线)、同一位置点侵彻体附近的水介质扩展速度(图7中vr表示侵彻体头部侵彻不同水层时使水介质产生的最大径向扩展速度)。通过对比发现,三个水层偏心亚半球缺罩水介质的最高速度均明显高于组合药型罩,表明在侵彻过程中偏心亚半球缺罩水介质获得的径向速度相对更高,也即偏心亚半球缺罩侵彻体在穿水过程中损失的能量相对更多。
通过对2种药型罩形成的侵彻体穿靶过程分析发现,侵彻体头部形状将直接影响组合药型罩在水中的侵彻性能以及二次毁伤效果。通过在偏心亚半球缺药型罩顶部设计小药型罩,小药型罩能够在爆轰波压垮作用下形成细长杆式射流,并在侵彻水介质和靶板过程中,前端杆式射流为后续侵彻体形成空腔通路。尽管如此,由于水介质和靶板对侵彻体很强的速度衰减作用并不断消耗头部侵彻体质量,导致在穿靶完成后,侵彻体只能保留尾部杵体的速度。
2.2.2 侵彻体穿过靶板过程速度以及动能衰减
为了进一步理解2种药型罩的穿靶过程,对2种药型罩形成的侵彻体头部位置及速度随时间变化、侵彻体动能以及侵彻体长径比随时间变化进行了分析,如图8所示。图8(a)中,曲线的斜率表征侵彻体头部速度的变化情况。分析发现,偏心亚半球缺罩头部侵彻体速度在穿靶过程中受到了较大的衰减,在1000 μs时,其斜率几乎为零,表明侵彻体并未能穿透后效靶板。而分析组合药型罩发现,在680 μs点处左右斜率不相同(左边斜率α小于右边斜率β),此时后效靶已完全穿透,可能原因为后端侵彻体与头部侵彻体汇聚导致速度突然增大,因此出现了头部侵彻体斜率突然增大的情况。分析图8(a)表明,对于穿靶过程,偏心亚半球缺罩侵彻体受到水的衰减作用更明显,通过对侵彻体动能随时间变化情况(图8(b))分析也能得到相同的结论。
对侵彻体的头部速度随位置变化曲线进行分析(见图9),发现在穿透后效靶板后剩余速度为1035 m/s,与杵体速度相近,说明在穿靶过程中速度梯度较大的侵彻体与杵体分离后在侵彻体前端不断被水介质和靶板消耗,速度衰减明显,虽然质量和速度均受到严重衰减,但杵体部分由于在空腔通道中低阻随进,因此在穿靶完成后,仍能保持(或接近)其最初的速度进一步完成毁伤作用,通过分析其构型发现(见图10),其径向尺寸Lx=24 mm,即0.2d2,轴向尺寸Ly=106 mm,长径比Ly/Lx=4.42。
分析结果表明,水介质对2种药型罩形成的侵彻体均有明显的速度衰减(偏心亚半球缺罩衰减了96%,组合药型罩速度衰减了80%)和动能衰减作用(偏心亚半球缺罩衰减了99%,组合药型罩速度衰减了92%),尽管对侵彻体的衰减作用很强,但由于组合药型罩其侵彻体头部有小药型罩形成的杆式射流开辟空腔通路,因此其速度和动能衰减趋势并没有球缺药型罩明显。研究结果表明,尽管穿靶过程水介质和靶板对侵彻体的有很强的衰减作用,通过将药型罩设计成组合结构,将前驱侵彻体设计成有利于穿水的杆式射流,可有效减弱水介质和靶板对侵彻体的衰减作用。
2.2.3 侵彻体对靶板毁伤机理分析
利用数值模拟结果研究了组合药型罩形成的侵彻体对靶板的毁伤机理。爆炸波在水中以压力波的形式传播,因此通过分析爆炸过程的压力波变化可对水中爆炸毁伤过程进行分析。
水中爆炸靶板可能受到侵彻体、爆炸冲击波、靶板应力波和气泡脉动等多种载荷的联合作用,通过对靶板表面的压力分析可以得出载荷对靶板的作用。第1层靶板由于与战斗部直接连接,因此将受到多种载荷联合作用。分析侵彻第1层靶板时在接近靶板表面炸高空腔(Q1点)和水中(Q2点)的压力发现(见图11),水中的峰值压力较空腔更早出现,且其峰值压力为673 MPa(85 μs),明显超过在炸高空腔中的靶板表面的峰值压力162 MPa(125 μs),表明爆炸形成的冲击波通过水介质传播并作用在第1层靶面,这种作用相比沿空气传播并作用在靶板上更快强度更高,但不如侵彻体作用的速度(聚能侵彻体在50 μs达到靶板)。因此分析结果表明,对于侵彻第1层靶板过程,首先爆炸形成的侵彻体高速冲击靶板,冲击产生的应力波在靶板中不断来回振荡并与冲击孔作用,这个过程会持续150 μs(即从侵彻体开始接触第1层靶面到尾部完全穿透所用时间),在此过程中,爆炸产生的冲击波通过水介质传播到靶板表面并与之发生作用,然后这种作用也会通过应力波方式最终作用在冲击孔处。
为了进一步分析爆炸沿水中传播的压力波对第一层靶板的作用,计算了当战斗部外部介质为空气时在相同位置处的压力随时间变化(即用空气替换包围在战斗部外部的水介质),分析发现(见图12),两种环境下Q1点处的峰值压力在接近(Q1点在空气介质中峰值压力及其达到时间分别为175 MPa和105 μs,在水介质中分别为162 MPa和125 μs)。但在战斗部外部P2点处,水介质时峰值压力明显高于空气介质处的峰值压力(P2点在空气介质中峰值压力和到达时间分别为2.71 MPa和120 μs,水介质中分别为673 MPa和85 μs),水介质下的峰值压力是空气介质下峰值压力的248倍,表明第1层靶板受到爆炸沿水介质传播冲击波的强烈作用。分析结果表明,当装药战斗部外部为水介质时,爆炸波一方面通过对药型罩的压垮形成金属侵彻体作用在靶板表面,也受到从水中传播过来的压力波作用,在两者共同作用下,导致了试验观察到的第1层靶板的撕裂。
试验结果表明,从第1层靶到第4层靶,受到毁伤程度逐渐减弱,表明侵彻体和冲击波对靶板作用强度被逐渐衰减,通过数值模拟方法可以计算出衰减的程度。分析在距离轴线40 mm附近各层靶板面的峰值压力值发现(见图13),在40 mm处爆炸产生的冲击波压力已经被衰减得非常小,最大不超过300 MPa,因此在距离40 mm处冲击波压力已经很难使靶板破坏或发生塑性变形,且随着距离增加,靶板表面的峰值压力明显减小。以上分析结果表明,从第2层靶开始,冲击波对结构作用已经不明显,主要是侵彻体对靶板的冲击作用以及靶板中应力波的作用。
以上分析结果表明,第1层靶在高速聚能侵彻体作用下冲击膨胀形成侵彻孔洞,并在孔洞附近形成压缩波向外传播,传播到靶板边缘反射成拉伸波,由于靶板径向尺寸相对较小,对孔洞形成拉伸效应明显,另一方面,由于距装药较近,爆炸形成的冲击波通过水介质传播到第1层靶面孔洞处,对靶板孔洞同样具有膨胀作用,因此第1层靶板将会受到侵彻体冲击作用和沿水介质传播的爆炸冲击波的联合作用。相较之下,第2层靶收到的冲击波作用的效应不明显,因此随后几层处于水中的靶板的毁伤效果并没有第1层强烈,随着水层厚度增加,水中的爆炸冲击波能量被逐渐衰减。
3. 组合药型罩聚能装药战斗部试验验证
3.1 试验设计
数值模拟计算结果表明,组合药型罩可以形成理想的空腔随进效应,因此对设计的结构展开试验验证。图14为设计的试验原理图,采用8号雷管中心起爆传爆药,传爆药再起爆主装药,端面采用磁环将整个战斗部吸附在靶板表面中心位置。采用数值模拟的结构布局尺寸以及靶板材料进行布置,水箱整体尺寸为300 mm×300 mm×950 mm,后效靶尺寸300 mm×300 mm×15 mm。整个水箱充满水介质。
3.2 试验结果
试验结果见图15,分析发现,四块靶板和后效靶均被完全穿透。第1块靶板完全碎裂,虽然仅找到两块碎片,但可以根据碎片评估侵彻体破孔尺寸,第2块和第3块靶板出现比较规则的圆孔,但整个靶板并未碎裂。第2块和第3块靶板在穿孔附近都出现了裂纹,但整体并未碎裂,第2块靶板比第3块靶板最长裂纹更长。第4块靶板上形成比较规则的圆孔,但孔附近并未出现裂纹,后效靶形成不规则圆孔,且靶板整体碎裂。结合数值模拟计算结果,表明第一层靶板受到了侵彻体的冲击作用以及爆炸冲击波联合作用,导致其沿对角线方向撕裂,而受到水介质对冲击波强烈的衰减作用,随后的第2层、第3层以及第四层靶板的冲击波的毁伤效果并不明显。
对靶板穿孔大小(入口孔径)进行了统计(见图16和表3),分析结果发现,组合药型罩数值模拟计算结果与试验结果相符,最大偏差为14.6%。
表 3 药型罩聚能装药穿孔尺寸数值模拟计算与试验结果对比Table 3. Comparison between numerical calculation and experimental results of perforation size靶板 偏心亚半球缺罩
破孔尺寸/mm组合药型罩破孔尺寸/mm 组合药型罩计算
结果偏差/%数值模拟结果 实验结果 第1层 49.6 46.6 50 −6.8 第2层 45.4 42.4 37 14.6 第3层 59.2 41.8 40 4.5 第4层 64.0 43.2 41 5.4 后效靶 未穿透 46.6 70×49 4.9 注:偏心亚半球缺罩破孔尺寸为数值模拟结果;组合药型罩计算结果偏差为数值模拟结果相对于试验结果的偏差。 4. 结 论
本文设计了一种组合药型罩结构并展开对含水复合结构的数值模拟研究,对组合药型罩形成侵彻体以及穿靶过程进行了分析,最后对设计的结构进行了试验验证,通过分析获得了以下结论。
(1) 设计了一种组合药型罩结构,在原有偏心亚半球缺药型罩结构的基础上,在外圆顶部轴线处设计偏心亚球缺罩,利用小药型罩在水下形成前驱杆式射流,辅助后续侵彻体在水中进行空腔随进运动,减弱水对后续侵彻体的作用。根据数值模拟结果表明设计的组合药型罩能形成理想空腔随进效果,对设计的组合药型罩展开试验验证,所有靶板均被完全穿透,靶板成孔尺寸数值模拟计算与试验结果符合较好,最大偏差不超过15%。
(2) 利用数值模拟方法对穿靶过程进行了分析,分析结果表明,在原有偏心亚半球缺罩基础上增加的小药型罩并不会明显影响原药型罩侵彻体的构型和速度,小药型罩在侵彻体头部形成细长的杆式射流,不仅增加了侵彻体总长,也增加了侵彻体头部的速度。通过对穿靶过程分析发现,侵彻体头部的形状将直接影响组合药型罩在水中的侵彻性能以及二次毁伤效果。通过将侵彻体头部设计成细长杆式射流,能形成理想的空腔随进效应。尽管如此,组合药型罩在完成水层的穿靶后,仅能保持其尾部杵体速度,水介质和靶板对侵彻体头部速度衰减作用非常明显。
(3) 通过数值模拟计算分析了组合药型罩穿靶机理。分析结果表明,第一层靶板不仅受到侵彻体高速冲击作用,还受到沿水传播的爆炸冲击波的强烈作用,而随着水层厚度增加,水中的爆炸冲击波能量被迅速衰减。对试验结果分析发现,由于第一层靶板受到了高速侵彻体冲击作用和爆炸冲击波联合作用,导致其沿对角线方向靶板的撕裂,而由于水介质对冲击波的衰减,导致其余靶板受爆炸冲击波作用并不明显。
-
表 1 材料本构模型及状态方程
Table 1. Constitutive model and state equation of the materials
材料 本构模型 状态方程 炸药 HIGH_EXPLOSIVE_BURN JWL 药型罩 STEINBERG GRÜNEISEN 壳体 JOHNSON_COOK GRÜNEISEN 水 NULL LINEAR_POLYNOMIAL 空气 NULL LINEAR_POLYNOMIAL 靶板 JOHNSON_COOK GRÜNEISEN 表 2 80 μs时组合药型罩和偏心亚半球缺罩的外形尺寸统计
Table 2. Statistics of overall dimensions of combined liner and eccentric sub-hemispherical liner at 80 μs
D1/mm D2/mm L1/mm S1/mm S2/mm 44 (0.36d2) 12 (0.10d2) 169 (1.39d2) 22 (0.18d2) 104 (0.85d2) S3/mm D3/mm D4/mm D5/mm L2/mm 43 18 46 (0.38d2) 12 (0.10d2) 103 (0.84d2) 表 3 药型罩聚能装药穿孔尺寸数值模拟计算与试验结果对比
Table 3. Comparison between numerical calculation and experimental results of perforation size
靶板 偏心亚半球缺罩
破孔尺寸/mm组合药型罩破孔尺寸/mm 组合药型罩计算
结果偏差/%数值模拟结果 实验结果 第1层 49.6 46.6 50 −6.8 第2层 45.4 42.4 37 14.6 第3层 59.2 41.8 40 4.5 第4层 64.0 43.2 41 5.4 后效靶 未穿透 46.6 70×49 4.9 注:偏心亚半球缺罩破孔尺寸为数值模拟结果;组合药型罩计算结果偏差为数值模拟结果相对于试验结果的偏差。 -
[1] 谭多望, 孙承纬. 成型装药研究新进展 [J]. 爆炸与冲击, 2008, 28(1): 50–56. DOI: 10.11883/1001-1455(2008)01-0050-07.TAN D W, SUN C W. Progress in studies on shaped charge [J]. Explosion and Shock Waves, 2008, 28(1): 50–56. DOI: 10.11883/1001-1455(2008)01-0050-07. [2] 梁争峰, 胡焕性. 爆炸成形弹丸技术现状与发展 [J]. 火炸药学报, 2004, 27(4): 21–25. DOI: 10.3969/j.issn.1007-7812.2004.04.006.LIANG Z F, HU H X. The current situation and future development direction of explosively formed projectile technology [J]. Chinese Journal of Explosives & Propellants, 2004, 27(4): 21–25. DOI: 10.3969/j.issn.1007-7812.2004.04.006. [3] 吴晗玲, 段卓平, 汪永庆. 杆式射流形成的数值模拟研究 [J]. 爆炸与冲击, 2006, 26(4): 328–332. DOI: 10.11883/1001-1455(2006)04-0328-05.WU H L, DUAN Z P, WANG Y Q. Simulation investigation of rod-like jets [J]. Explosion and Shock Waves, 2006, 26(4): 328–332. DOI: 10.11883/1001-1455(2006)04-0328-05. [4] 徐斌, 王成, 徐文龙. 高速杆式射流形成的数值模拟与实验研究 [J]. 北京理工大学学报, 2018, 38(4): 331–337. DOI: 10.15918/j.tbit1001-0645.2018.04.001.XU B, WANG C, XU W L. Numerical simulation and experiment investigation on hypervelocity jetting projectile charge formation [J]. Transactions of Beijing Institute of Technology, 2018, 38(4): 331–337. DOI: 10.15918/j.tbit1001-0645.2018.04.001. [5] 王利侠, 谷鸿平, 丁刚, 等. 聚能射流对带壳浇注PBX装药的撞击响应 [J]. 含能材料, 2015, 23(11): 1067–1072. DOI: 10.11943/j.issn.1006-9941.2015.11.006.WANG L X, GU H P, DING G, et al. Reaction characteristics for shelled cast-cured PBX explosive impacted by shaped charge jet [J]. Chinese Journal of Energetic Materials, 2015, 23(11): 1067–1072. DOI: 10.11943/j.issn.1006-9941.2015.11.006. [6] 王海福, 江增荣, 李向荣. 药型罩参数对聚能装药水下作用效应的影响 [J]. 北京理工大学学报, 2006, 26(5): 405–409. DOI: 10.3969/j.issn.1001-0645.2006.05.007.WANG H F, JIANG Z R, LI X R. Influences of liner parameters on the effects of shaped charge operating underwater [J]. Transactions of Beijing Institute of Technology, 2006, 26(5): 405–409. DOI: 10.3969/j.issn.1001-0645.2006.05.007. [7] 周方毅, 詹发民, 姜涛, 等. 一种组合药型罩聚能战斗部 [J]. 鱼雷技术, 2012, 20(5): 380–383; 400. DOI: 10.3969/j.issn.1673-1948.2012.05.014.ZHOU F Y, ZHAN F M, JIANG T, et al. An idea about shaped charge warhead with combined charge liner for torpedo [J]. Torpedo Technology, 2012, 20(5): 380–383; 400. DOI: 10.3969/j.issn.1673-1948.2012.05.014. [8] 周方毅, 黄雪峰, 詹发民, 等. 一种双球缺组合药型罩聚能鱼雷战斗部研究 [J]. 水下无人系统学报, 2017, 25(4): 278–281; 287. DOI: 10.11993/j.issn.2096-3920.2017.03.011.ZHOU F Y, ZHAN F M, JIANG T, et al. A shaped charge warhead with two spherical combined liners for torpedo [J]. Journal of Unmanned Undersea Systems, 2017, 25(4): 278–281; 287. DOI: 10.11993/j.issn.2096-3920.2017.03.011. [9] 张春辉, 张斐, 王志军, 等. 复合材质杆式射流侵彻水下目标的数值模拟 [J]. 爆破器材, 2019, 48(1): 8–14. DOI: 10.3969/j.issn.1001-8352.2019.01.002.ZHANG C H, ZHANG F, WANG Z J, et al. Numerical simulation of composite-material rod-like jet penetrating underwater targets [J]. Explosive Materials, 2019, 48(1): 8–14. DOI: 10.3969/j.issn.1001-8352.2019.01.002. [10] 李兵, 刘念念, 陈高杰, 等. 水中聚能战斗部毁伤双层圆柱壳的数值模拟与试验研究 [J]. 兵工学报, 2018, 39(1): 38–45. DOI: 10.3969/j.issn.1000-1093.2018.01.004.LI B, LIU N N, CHEN G J, et al. Numerical simulation and experimental research on damage of shaped charge warhead to double-layer columniform shell [J]. Acta Armamentarii, 2018, 39(1): 38–45. DOI: 10.3969/j.issn.1000-1093.2018.01.004. [11] 王玉, 卢熹, 张方方, 等. 反潜鱼雷战斗部对典型潜艇目标毁伤效应研究 [J]. 兵器装备工程学报, 2021, 42(12): 112–116. DOI: 10.11809/bqzbgcxb2021.12.016.WANG Y, LU X, ZHANG F F, et al. Damage effect of anti-Submarine torpedo warhead on typical submarine targets [J]. Journal of Ordnance Equipment Engineering, 2021, 42(12): 112–116. DOI: 10.11809/bqzbgcxb2021.12.016. [12] 王长利, 马坤, 周刚, 等. 防雷舱结构在聚能装药水下爆炸作用下的毁伤研究 [J]. 爆炸与冲击, 2018, 38(5): 1145–1154. DOI: 10.11883/bzycj-2017-0119.WANG C L, MA K, ZHOU G, et al. Damage effect of cabin near ship board under shaped charge exploding underwater [J]. Explosion and Shock Waves, 2018, 38(5): 1145–1154. DOI: 10.11883/bzycj-2017-0119. [13] 王长利, 周刚, 马坤, 等. 典型含水复合结构在聚能装药水下爆炸作用下的毁伤 [J]. 船舶力学, 2018, 22(8): 1001–1010. DOI: 10.3969/j.issn.1007-7294.2018.08.010.WANG C L, ZHOU G, MA K, et al. Damage anlysis of typical water partitioned structure under shaped charge underwater explosion [J]. Journal of Ship Mechanics, 2018, 22(8): 1001–1010. DOI: 10.3969/j.issn.1007-7294.2018.08.010. [14] 杨莉, 张庆明, 汪玉, 等. 反舰聚能战斗部装药结构研究 [J]. 兵工学报, 2009(S2): 154–158.YANG L, ZHANG Q M, WANG Y, et al. Research on shaped charge warhead of anti-ship missile [J]. Acta Armamentarii, 2009(S2): 154–158. [15] 傅磊, 王伟力, 李永胜, 等. 组合药型罩水介质中成型的数值仿真 [J]. 鱼雷技术, 2015, 23(5): 367–373. DOI: 10.11993/j.issn.1673-1948.2015.05.009.FU L, WANG W L, LI Y S, et al. Numerical simulation of combined liner formation in water [J]. Torpedo Technology, 2015, 23(5): 367–373. DOI: 10.11993/j.issn.1673-1948.2015.05.009. [16] 时党勇, 李裕春, 张胜民. 基于ANSYS/LS-DYNA 8.1进行显式动力分析 [M]. 北京: 清华大学出版社, 2005: 295–296. [17] 时党勇. 倾斜尾翼爆炸成型弹丸的数值模拟和外弹道计算 [C]// 第十届全国爆炸与安全技术会议论文集. 昆明: 2011, 286–292. [18] LEE S G, BAEK Y H, LEE I H, et al. Numerical simulation of 2D sloshing by using ALE2D technique of LS-DYNA and CCUP methods [C]// The Proceedings of the Twentieth (2010) International Offshore and Polar Engineering Conference (ISOPE-2010 Beijing). International Society of Offshore and Polar Engineers (ISOPE), 2010: 192–199. [19] MULLIN M J, O’TOOLE B J. Simulation of energy absorbing materials in blast loaded structures [C]// 8th International LS-DYNA Users Conference. 2004: 2–7. [20] 吴海军, 王可慧, 柯明, 等. 多点同步起爆条件下环形射流成型及侵彻过程的数值模拟 [J]. 现代应用物理, 2018, 9(2): 75–82. DOI: 10.12061/j.issn.2095-6223.2018.021002.WU H J, WANG K H, KE M, et al. Numerical simulation on formation and penetration processes of the annular jet with multi-points synchronous explosive circuit [J]. Modern Applied Physics, 2018, 9(2): 75–82. DOI: 10.12061/j.issn.2095-6223.2018.021002. [21] ROBBINS J R, DING J L, GUPTA Y M. Load spreading and penetration resistance of layered structures-a numerical study [J]. International Journal of Impact Engineering, 2004, 30(6): 593–615. DOI: 10.1016/j.ijimpeng.2003.08.001. 期刊类型引用(1)
1. 宗香华,王银,孔祥振,姜雅婷,孙留洋,袁俊成,杨涛春. 错位多次打击下UHPC靶体损伤破坏效应的数值模拟研究. 高压物理学报. 2024(03): 161-173 . 百度学术
其他类型引用(3)
-