Modes and influencing factors of electromagnetically driven high velocity formed projectile
-
摘要: 为了研究电磁驱动药型罩形成高速成型弹丸的可行性及弹丸成型特性,基于CQ-7脉冲功率装置,开展了电磁驱动线性药型罩形成弹丸技术实验。结合激光多普勒测速技术,实现了电磁驱动药型罩形成弹丸的速度测量和侵彻铝靶验证。同时,基于流体动力学软件和相应电磁仿真模块,建立了电磁驱动弹丸成型的物理模型和数值模拟方法,模拟了弹丸成型和侵彻铝靶的动力学过程,利用实验结果验证了数值模拟方法的可靠性。在此基础上,研究了等壁厚球缺型药型罩的结构参数以及加载能量对弹丸成型参数的影响规律。结果表明:外曲率半径对弹丸的头部速度影响较小,而头部速度会随壁厚的减小和加载能量的增大显著增加;弹丸的长径比随外曲率半径和壁厚的减小、加载能量的增大呈逐渐增加的趋势。最后,利用数值模拟方法预测并验证了利用电磁驱动技术获得高速度和大质量成型弹丸的可行性。Abstract: To investigate the feasibility and characteristics of high-velocity formed projectile formation driven by electromagnetic loading, exploratory experiments of projectile formation by electromagnetically driven the linear liner were conducted using the pulsed power generator CQ-7. Photon Doppler velocimeter (PDV) was employed to measure the velocity of the electromagnetic-driven projectiles and validate their penetration into aluminum targets. A physical model and numerical simulation method for electromagnetic-driven projectile formation were established using fluid dynamics software and corresponding electromagnetic simulation modules. The changes in current density and magnetic pressure during the electromagnetic loading stage were studied and the dynamic processes of projectile formation and penetration into aluminum targets were simulated. The numerical simulation method was verified through the comparison between numerical results and experimental data. Based on this, the influences of liner configuration and loading energy on the projectile formation parameters of equal wall thickness hemispherical liner were explored. The results indicate that the outer curvature radius has a minor impact on the head velocity of the projectile, while the head velocity significantly increases with decreasing wall thickness and increasing loading energy. The aspect ratio of the projectile gradually increases with decreasing outer curvature radius and wall thickness, as well as increasing loading energy. The conversion between quasi-spherical and long rod-shaped projectile modes can be achieved by changing the structural parameters, and for the same structural parameter, the conversion between two modes can be achieved by controlling the loading energy. Finally, the feasibility of obtaining high-velocity and high-mass-formed projectiles using electromagnetic-driven technology was predicted using numerical simulation methods, and it can be figured out from the results that a projectile with a higher velocity and larger mass can be formed by increasing the loading energy and the sizes of the shaped liner, effectively breaking through the velocity limit of a traditional penetrator driven by explosive detonation.
-
金属在众多领域中的应用越来越广泛,其生产加工过程中的安全隐患日益突出。在冶金行业,铁水、钢水等金属液体在加工过程中,极易发生喷溅、泄露等事故[1-2]。这些高温熔融金属一旦与低温冷却剂相接触就会发生剧烈混合,引发蒸汽爆炸,造成不同程度的人员伤亡和财产损失。在这类沸腾液体扩展蒸汽爆炸的研究中,人们采用高速摄像系统、红外摄像仪等,观测水柱或单个水滴撞击高温熔融金属液表面[3-8]、熔融金属液柱或单个金属液滴落入冷却水中[9-13]的熔融物与冷却剂相互作用(fuel-coolant interaction, FCI)的过程。其中,关于低熔点锡与水的相互作用,也有以下实验研究。张荣金等[8]、林栋[13]在单个水滴撞击高温熔融锡液的实验中获知,锡液与水的反应剧烈程度随韦伯数的增大和熔融锡液温度的升高而增加;王骞[14]基于熔融锡液柱和冷却水接触反应现象获知,在液柱下落过程中,升高锡液温度和增大液柱直径会促进局部蒸汽爆炸的发生,同时通入一定惰性气体可减弱锡液表面的氧化程度,从而提升蒸汽爆炸威力;沈致远[15]通过调节炉温、水温和下落高度对锡液与水的相互作用机理进行研究得到,锡液细粒化程度与锡液温度成正比、与水温成反比,而下落高度的变化对金属细粒化程度影响不显著;李天舒等[16]在熔融锡液入水碎化实验研究中,展示了熔融金属热扩散系数对蒸汽爆炸的影响。此外,胡逊祥等[17]通过分析熔融金属锡与水作用过程中的运动图像和压力波动曲线提出,压力大小由进入锡液内部的液态冷却剂的质量和内部液态冷却剂所受加热膨胀程度决定。这些对于熔融锡液与水接触反应过程中的影响因素的探讨,为后续研究提供了一定理论指导,而在蒸汽爆炸过程中冲击波能量的定量分析尚有不足。
本文中,为了研究低熔点金属锡遇水爆炸机理,采用高速摄像机观测熔融锡液柱落入反应槽中并与槽内冷却水接触反应的过程;着重计算蒸汽爆炸过程中冲击波能量转换率,以期获得具有理论参考价值的结论。实验中,设计不同熔融锡液与水的质量比,以确保蒸汽爆炸,并了解金属与水的质量比对蒸汽爆炸过程的影响;同时,使熔融锡液温度远高于锡的熔点,以提高锡液的流动性[18],增加熔融锡液与水接触面积,使反应充分进行。在此基础上,选取中高熔点金属铝,进行在相同条件下的蒸汽爆炸实验,将熔融锡液遇水爆炸实验与周宁等[19]熔融铝液遇水爆炸实验作比较,进一步探究低熔点金属锡遇水爆炸反应特征。
1. 实 验
实验系统[19]主要包括高频熔融炉、高速摄像机(I-SPEED 3)、温度信号采集仪、高频信号采集仪、K型热电偶(WRNK-010,测温范围为0~1 300 ℃,相对误差为±0.75%)、压力传感器(KD2002-03X,测压范围为0~10 MPa,误差为±0.1 MPa)、反应槽和防护装置等,其中压力传感器呈90°安置在距离反应槽顶部边沿0.3 m处,热电偶放置于高频熔融炉腔内,如图1所示。实验环境温度为27 ℃,锡与水的质量比为5、10、15和20,其中锡的熔融态温度为1 050 ℃,反应结束时锡的温度为100 ℃。
实验过程如下。(1) 将一定量的干燥锡块放入高频熔融炉中,加热熔化至1 050 ℃。用调功器调整熔融炉的频率,控制炉内锡的熔融态温度;利用熔融炉中的转动装置可使锡受热均匀、充分熔化。(2) 采用遥控转动装置,将高频熔融炉平稳翻转一定角度,让熔融锡液以较均匀速度经导流管下落至反应槽中。(3) 采用高速摄像机,拍摄锡液与水的接触反应过程,温度测试系统和压力测试系统实时采集测点处的数据,获取的图像和数据传输并保存与电脑,用于后续分析。(4) 按步骤(1)~(3),在相同条件下,进行金属铝的对比实验。
2. 计算模型
2.1 熔融锡液内部储能
熔融锡液与水之间的快速热交换是蒸汽爆炸产生的动力,高温熔融锡液的总热量为:
Qr=mr(∫TnT0cdT+∫TTncprdT+γr) (1) 因金属锡的固态热容和液态热容随温度的变化不大[20],可将式(1)简化为:
Qr=mrc(Tn−T0)+mrcpr(T−Tn)+mrγr (2) 式中:mr为金属锡的质量,T0为室温,Tn为锡的熔点,T为锡的熔融态温度,c为锡的固态比热容,cpr为锡的液态比热容,γr为锡的熔化热。
2.2 冲击波能量
爆炸冲击波在空气中传播时,距离爆炸中心处的超压与该处TNT炸药药量对应。
利用压力传感器监测的距爆炸中心不同位置的冲击波压力,并基于不同爆炸环境下的Brode公式[21-22],可计算熔融锡液与水接触爆炸冲击波的TNT炸药药量:
Δp=0.096ˉr+0.143ˉr2+0.573ˉr3−0.00190.0098MPa≤Δp≤0.98 MPa (3) ˉr=3√wR (4) 式中:∆p为空气冲击波超压,R为距离爆炸中心的距离,
ˉr 为比例距离,w为TNT炸药药量。获得TNT炸药药量后,可转化为爆炸冲击波的能量为:
Qrb=QTNTw (5) 式中:Qrb为高温熔融锡液热量转化为爆炸能量的部分,QTNT=4 184 kJ/kg为TNT炸药爆炸的质量热量。
2.3 冲击波能量转化率
熔融锡液与水发生蒸汽爆炸时,冲击波能量的转化率为:
η=Qrb/Qr (6) 式中:η为熔融锡液热能转化为爆炸冲击波能量的转化率。
3. 实验结果分析
3.1 蒸汽爆炸过程
图2为n=5时(n为金属与水的质量比)高温熔融锡液与水的接触反应过程。高温熔融物与冷却剂的相互作用过程如下。 熔融锡液与水接触反应,在0~0.04 s时,为熔融锡液与水的粗混合阶段:锡液液柱受水力作用初步破碎成锡液液滴,并在水中分散开来;同时,低温水吸收熔融锡液的热量,开始升温、沸腾;系统内的熔融锡液与水逐步进入亚稳态状态,该阶段无明显反应特征。从0.06 s始,因锡液表面蒸汽膜塌陷,反应进入触发阶段,熔融锡液与水的直接接触使水槽上方及周边出现了大量锡液细粒和小范围的高温高压蒸汽团。随着高温熔融锡液持续不断地滴落至水槽,传热过程快速传播,高压蒸汽团内的压力不断叠加,反应进入膨胀阶段。在0.12 s时,系统发生第1次爆炸,在0.20 s时,紧接着发生第2次爆炸。爆炸过程中,锡液细粒化现象明显,飞溅加剧、范围加大,水槽上方及附近的蒸汽团以可见速度向周边迅速扩散,并伴有巨大响声和刺眼火光。2次爆炸结束后,水槽周边的锡液飞溅逐渐减弱,蒸汽团体积逐步减小,在1.04 s时蒸汽团几乎散尽。实验结束后,测得水槽内剩余锡块的质量,获得熔融锡液质量喷溅率高达33%。
图3为实验的冲击波压力曲线。随着高温熔融锡液与水的反应深入,测点压力在约0.12 s出现第1次突跃变化,反应槽上方高压蒸汽团向外扩散、锡液滴不断飞溅,冲击波能量迅速消耗,波阵面后压力急剧下降,在0.13~0.16 s出现负压区;在压力由负逐步回升至正并趋于平衡的过程中,在约0.2 s出现第2次压力突跃变化,极短时间内压力叠加,使第2次爆炸后的压力上下波动较大。熔融金属遇水爆炸达到一定规模时,往往伴随着水分解产生的氢气与空气中的氧气反应发生的化学爆炸,但水的热分解需在高温下触发,如温度达不到1 600 ℃,便不能分解成氢气和氧气[23];而实验中锡液加热的最高温度为1 050 ℃,无法使水受热分解,因此这2次爆炸均为物理爆炸。2次蒸汽爆炸产生原因,可能为:(1) 在第1次蒸汽爆炸发生后,随着熔融锡液持续落入水槽中与剩余的水反应,短时间内再次经历蒸汽爆炸发展的4个阶段[24],引发2次爆炸;(2) 在触发阶段,产生了两个尺寸不同的局部高压区,以致经后期传播、膨胀阶段后,在不同时间上发生了2次爆炸。
图4为n=5时相同条件下的熔融铝液与水的相互作用过程。在0.14 s时铝液明显碎化,在0.46 s时触发蒸汽爆炸,铝液剧烈飞溅和蒸汽团对外膨胀的过程持续了约0.32 s,铝液飞溅量为铝总质量的约15%。
与熔融铝液实验对比,熔融锡液遇水爆炸反应过程更剧烈、反应持续时间更长。究其原因,蒸汽爆炸反应的剧烈程度受金属碎化程度的影响,两者间呈正相关[12,24],因此细粒化程度较高的锡液反应更剧烈。在非稳态导热过程中,锡的固体比热容和液体比热容均大于铝的,锡的导热系数低于铝的;由比热容和导热系数构成的热扩散率可知,锡内部温度趋于一致的能力弱,表现为熔融锡液内部的热量向金属锡液表面传递的速度小,包裹金属表面的蒸汽膜更稳定,反应持续时间也就更长。
此外,由能量转换可知,熔融锡液遇水爆炸产生的冲击波能量来源于锡液内部热量,但这只是总热量的一部分。除仍以热的形式传递给周围环境或残留在锡液内部的能量,实验中观察到的强光为部分能量转化为光辐射能(见图2中0.12~0.28 s时);爆炸产生的巨大响声为部分能量转化为声能,而爆炸产物飞散为锡液内部热能转化为动能。
3.2 能量分析
3.2.1 不同质量比时锡遇水爆炸
当锡与水的质量比为5、10、15和20时,实验的熔融锡液内部储能、超压峰值、冲击波能量的转化率等,见表1。表中,mw为水的质量。冲击波能量转化率曲线如图5所示。
表 1 高温熔融锡液遇水的实验数据Table 1. Experimental data of high-temperature molten tin liquid contacted with watermr/kg mw/kg n Qr/kJ ∆p/MPa w/kg η/% 0.7031 0.140 5 199.36 0.22 0.0031 6.51 0.7224 0.144 204.84 0.24 0.0035 7.15 0.7490 0.150 212.38 0.25 0.0037 7.29 1.1230 0.225 318.43 0.28 0.0043 5.65 0.5944 0.059 10 168.54 0.27 0.0041 10.18 0.6108 0.061 173.19 0.29 0.0045 10.87 0.6355 0.063 180.20 0.30 0.0047 10.91 0.6452 0.065 182.95 0.30 0.0047 10.75 1.98 0.132 15 467.50 0.156 0.00348 3.11 2.01 0.134 473.67 0.272 0.00728 6.43 2.18 0.145 508.63 0.119 0.00237 1.95 2.26 0.151 525.08 0.242 0.00625 4.98 1.86 0.093 20 442.82 0.043 0.00048 0.45 2.00 0.100 471.61 0.072 0.00112 0.99 2.21 0.111 514.80 0.068 0.00102 0.83 2.28 0.114 529.20 0.065 0.00095 0.75 由实验数据,高温熔融锡液遇水爆炸过程中,0.45%~10.91%的热量转化为爆炸冲击波能量,质量比对爆炸过程中冲击波能量转化率影响较大。由图5可见,在n=10时冲击波能量转化率显然高于其他质量比的。对比n=15, 20时的冲击波能量转化率,可见:在相近锡质量时,因n=20的冷却水的质量占比少,反应槽中冷却水量在还未完全吸收熔融锡液内部储能时就已耗尽,所以熔融锡液与水之间的热传递总量减少,无法形成更大的蒸汽高压区,因此n=20的爆炸冲击波能量转化率较低。
此外,在相同质量比中,锡质量的增加一定程度上提高了超压(见表1中n=5, 10时),但两者并非正相关的增长关系,随着锡质量的持续增加,超压在一定质量范围内出现先增后减的趋势(见表1中n=15, 20时)。
3.2.2 锡/铝遇水爆炸
图6为熔融锡/熔融铝遇水爆炸实验的平均冲击波能量转化率随质量比的变化曲线。在实验中, 在n=5时,锡遇水爆炸的平均能量转化率较大;由FCI过程可见,冲击波能量转化率受反应剧烈程度和持续时间的影响,进而受与这2个量相关的金属碎化程度和金属热扩散率的影响。同时,在金属碎化程度和金属热扩散率的影响下,相比铝遇水爆炸,锡遇水爆炸的冲击波能量转化率峰值的位置,向质量比变小的方向偏移;在n≈10时,锡液内部储能转化为冲击波能量的占比最大,可达10.38%,而在n≈15时,铝遇水爆炸的冲击波能量转化率到达最大值8.43%。
2种金属遇水爆炸的平均冲击波能量转化率整体变化趋势一致,均随质量比增大先增后减。以两者曲线交点的质量比为界限,该点n≈12.69,
¯η ≈7.69;由实验,在n<12.69时,锡液遇水爆炸的平均冲击波能量转化率高于铝液遇水爆炸的冲击波能量转化率。4. 结 论
(1) 在n=5时,锡/铝遇水爆炸实验中,锡的热扩散率低,反应持续时间更长;熔融锡液的碎化程度高,锡液与水接触反应过程更剧烈。就能量转化而言,熔融锡液内部储能存在多种形式的耗散,有热能、光辐射能、声能、冲击波等;同时,熔融锡/熔融铝总热能的冲击波能量转化率受金属热扩散率和金属遇水后碎化程度的影响。
(2) 在相同质量比时,锡的质量增加一定程度上提高了超压峰值,但这种增长趋势不是持续的、绝对的,在 n=15, 20时,在一定质量范围超压出现了先增后减的变化。
(3) 在不同质量比时,熔融锡液遇水爆炸过程中0.45%~10.91%热能转化为冲击波能量;随着质量比的增加,冲击波能量转化率呈先增后减的趋势;在n=10时,冲击波能量转化率最大。
(4) 熔融锡/熔融铝遇水爆炸过程的平均能量转化率随着质量比的增大都呈先增后减的趋势;在n<12.69时,锡液遇水爆炸的冲击波能量转化率高于铝液遇水爆炸的冲击波能量转化率。
-
表 1 CQ-7的电参数
Table 1. Electrical parameters of CQ-7
电容/μF 电感/nH 电阻/mΩ 充电峰值电压 上升时间 20.48 4.12 3.35 ±65 kV 400~700 ns 表 2 实验条件
Table 2. Experimental conditions
实验 实验类型 充电电压/kV Shot-1 PDV测速 ±45 Shot-2 侵彻铝靶 ±45 Shot-3 PDV测速 ±55 Shot-4 侵彻铝靶 ±55 表 3 无氧铜的Burgess电导率模型参数
Table 3. Parameters of Burgess electrical resistivity model for Cu-OFHC
V0/(cm3·g−1) γ0 θm,0/eV LF/(kJ·mol−1) K 0.112 2.00 0.117 0.13 0.964 表 4 实验与数值模拟侵彻数据对比
Table 4. Comparison of the penetration data between experiments and numerical simulations
编号 方法 侵彻深度/mm 相对误差/% 开孔尺寸/mm 相对误差/% Shot-2 实验 26.5 7.1 14.2×16.3 4.2 计算 28.4 14.8 Shot-4 实验 32.3 2.8 17.1×16.7 −3.0 计算 33.2 16.6 -
[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]. 高压物理学报, 2006, 20(4): 429–433. DOI: 10.11858/gywlxb.2006.04.015.YANG J, JIANG J W, MEN J B. Numerical simulation for formation flight and penetration of sphericity EFP [J]. Chinese Journal of High Pressure Physics, 2006, 20(4): 429–433. DOI: 10.11858/gywlxb.2006.04.015. [3] 刘建青, 顾文彬, 徐浩铭, 等. 多点起爆装药结构参数对尾翼EFP成型的影响 [J]. 含能材料, 2014(5): 594–599. DOI: 10.3969/j.issn.1006-9941.2014.05.004.LIU J Q, GU W B, XU H M, et al. Effects of multi-point initiation charge configuration parameters on EFP with fins formation [J]. Chinese Journal of Energetic Materials, 2014(5): 594–599. DOI: 10.3969/j.issn.1006-9941.2014.05.004. [4] 郭莎, 任新联, 周涛, 等. 翻转成型大长径比爆炸成型弹丸的数值模拟 [J]. 科学技术与工程, 2019, 19(27): 272–276. DOI: 10.3969/j.issn.1671-1815.2019.27.039.GUO S, REN X L, ZHOU T, et al. Numerical simulation of large length diameter ratio overturn molding explosively formed penetrator [J]. Science Technology and Engineering, 2019, 19(27): 272–276. DOI: 10.3969/j.issn.1671-1815.2019.27.039. [5] 黄炫宁, 李伟兵, 程伟, 等. 锥弧结合罩形成长杆状密实EFP的可行性 [J]. 含能材料, 2019, 27(2): 90–96. DOI: 10.11943/CJEM2018051.HUANG X N, LI W B, CHENG W, et al. Feasibility of the formation of long rod-shaped compacted explosively formed penetrator by cone-arc liner [J]. Chinese Journal of Energetic Materials, 2019, 27(2): 90–96. DOI: 10.11943/CJEM2018051. [6] 王伟, 徐琳, 王玥兮, 等. 准球形EFP成形因素的正交优化设计与试验验证 [J]. 兵器材料科学与工程, 2020, 43(5): 22–25. DOI: 10.14024/j.cnki.1004-244x.20200513.001.WANG W, XU L, WANG Y X, et al. Orthogonal optimization design and experimental study on formation process of quasi-spheral explosively formed projectile [J]. Ordnance Material Science and Engineering, 2020, 43(5): 22–25. DOI: 10.14024/j.cnki.1004-244x.20200513.001. [7] 张雪朋, 刘亚昆, 伊建亚, 等. 复合装药包覆式活性侵彻体成型及侵彻研究 [J]. 兵器装备工程学报, 2021, 42(7): 1–5. DOI: 10.11809/bqzbgcxb2021.07.001.ZHANG X P, LIU Y K, YI J Y, et al. Study on formation and penetration of the wrapped reactive projectile formed by double-layer shaped charge [J]. Journal of Ordnance Equipment Engineering, 2021, 42(7): 1–5. DOI: 10.11809/bqzbgcxb2021.07.001. [8] 林加剑, 贾虎. 爆炸成型弹丸有效装药结构理论分析及试验研究 [J]. 弹箭与制导学报, 2015, 35(1): 59–62,67. DOI: 10.15892/j.cnki.djzdxb.2015.01.016.LIN J J, JIA H. Theoretical analysis and experimental research on the effective shaped charge with EFP [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2015, 35(1): 59–62,67. DOI: 10.15892/j.cnki.djzdxb.2015.01.016. [9] 张旭平. 电磁驱动实验技术及其加载下聚苯乙烯的动态行为研究 [D]. 绵阳: 中国工程物理研究院, 2019. [10] DEGNAN J H, BAKER W L, ALME M L, et al. Multimegajoule electromagnetic implosion of shaped solid-density liners [J]. Fusion Technology, 1995, 27(2): 115–123. DOI: 10.13182/FST95-A30368. [11] DEGNAN J H, TACCETTI J M, CAVAZOS T, et al. Implosion of solid liner for compression of field reversed configuration [J]. IEEE Transactions on Plasma Science, 2001, 29(1): 93–98. DOI: 10.1109/27.912947. [12] DOU J H, JIA X, HUANG Z X, et al. Theoretical and numerical simulation study on jet formation and penetration of different liner structures driven by electromagnetic pressure [J]. Defence Technology, 2021, 17(3): 846–858. DOI: 10.1016/j.dt.2020.05.016. [13] DOU J H, JIA X, HUANG Z X, et al. Theoretical study of the jet formation of a shaped charge liner driven by strong electromagnetic energy [J]. IEEE Transactions on Plasma Science, 2019, 47(12): 5283–5290. DOI: 10.1109/tps.2019.2951112. [14] 王桂吉, 罗斌强, 陈学秒, 等. 磁驱动平面准等熵加载装置、实验技术及应用研究新进展 [J]. 爆炸与冲击, 2021, 41(12): 121403. DOI: 10.11883/bzycj-2021-0119.WANG G J, LUO B Q, CHEN X M, et al. Recent progress on the experimental facilities, techniques and applications of magnetically driven quasi-isentropic compression [J]. Explosion and Shock Waves, 2021, 41(12): 121403. DOI: 10.11883/bzycj-2021-0119. [15] 王桂吉. 磁驱动等熵压缩和飞片加载技术和实验研究 [D]. 绵阳: 中国工程物理研究院, 2007.WANG G J. Research on magnetically driven isentropic compression and flyer plates [D]. Mianyang: China Academy of Engineering Physics, 2007. [16] 章征伟. 磁驱动固体套筒内爆理论与实验研究 [D]. 绵阳: 中国工程物理研究院, 2020.ZHANG Z W. Theoretic and experimental study on magnetically driven solid linerimplosion [D]. Mianyang: China Academy of Engineering Physics, 2020. [17] KNUDSON M D, LEMKE R W, HAYES D B, et al. Near-absolute Hugoniot measurements in aluminum to 500 GPa using a magnetically accelerated flyer plate technique [J]. Journal of Applied Physics, 2003, 94(7): 4420–4431. DOI: 10.1063/1.1604967. [18] GRACE F, DEGNAN J, ROTH C, et al. Shaped charge jets driven by electromagnetic energy [C]// Proceedings of the 28th International Symposium of Conference. Atlanta: International Ballistics Society, 2013: 15–26. [19] HUANG B Y, ZHANG X P, WANG G J, et al. Shaped charge liner collapse and jet formation by electromagnetic loading on high pulsed power generator [J]. IEEE Transactions on Plasma Science, 2023, 51(10): 3140–3151. DOI: 10.1109/TPS.2023.3320665. [20] CHEN X M, LUO B Q, ZHANG X P, et al. A compact pulsed power driver with precisely shaped current waveforms for magnetically driven loading experiments [J]. Review of Scientific Instruments, 2022, 93(8): 083910. DOI: 10.1063/5.0089939. [21] DOLAN D H. Extreme measurements with photonic Doppler velocimetry (PDV) [J]. Review of Scientific Instruments, 2020, 91(5): 051501. DOI: 10.1063/5.0004363. [22] ZELLNER M B, VUNNI G B. Photon Doppler velocimetry (PDV) characterization of shaped charge jet formation [J]. Procedia Engineering, 2013, 58: 88–97. DOI: 10.1016/j.proeng.2013.05.012. [23] L'EPLATTENIER P, COOK G, ASHCRAFT C, et al. Introduction of an electromagnetism module in LS-DYNA for coupled mechanical-thermal-electromagnetic simulations [J]. Steel Research International, 2009, 80(5): 351–358. DOI: 10.2374/SRI08SP152. [24] L'EPLATTENIER P, ÇALDICHOURY I. Recent developments in the electromagnetic module: a new 2D axi-symmetric EM solver [C]//Proceedings of the 10th European LS-DYNA Conference. Würzburg, German, 2015. [25] 张旭平, 赵剑衡, 谭福利, 等. 磁驱动飞片的三维数值模拟及分析 [J]. 高压物理学报, 2014, 28(4): 483–488. DOI: 10.11858/gywlxb.2014.04.015.ZHANG X P, ZHAO J H, TAN F L, et al. Three-dimensional numerical simulation and analysis of magnetically driven flyer plates [J]. Chinese Journal of High Pressure Physics, 2014, 28(4): 483–488. DOI: 10.11858/gywlxb.2014.04.015. [26] BURGESS T J. Electrical resistivity model of metals [C]//Proceedings of the 4th International Conference on Megagauss Magnetic-Field Generation and Related Topics. Santa Fe, NM, USA, 1986. -