电磁能与炸药联合加载药型罩基础理论

豆剑豪;贾鑫;梁争峰;黄正祥;薛标

兵工学报 ›› 2023, Vol. 44 ›› Issue (6) : 1754-1763.

兵工学报 ›› 2023, Vol. 44 ›› Issue (6) : 1754-1763. DOI: 10.12382/bgxb.2022.0120

电磁能与炸药联合加载药型罩基础理论

作者信息 +

Basic Theory of Loading Liner With Electromagnetic Energy and Explosives

Author information +
本文二维码
文章历史 +

摘要

为提高聚能射流对炸药能量的利用率,结合电磁加载药型罩以及爆磁压缩发生器原理,提出一种电磁能与炸药联合加载药型罩的原理结构。联合加载分为两个阶段:第1阶段,电容器放电,在爆磁压缩发生器中形成初始磁通;第2阶段,炸药1爆轰,压缩磁场,电流增大,同时爆轰波传播至药型罩部分引爆炸药2,电磁力与爆轰压力共同压垮药型罩。结合炸药爆轰理论、传统PER射流成型理论以及电磁加载特性,建立电磁与炸药联合加载下药型罩射流成型的理论模型。利用理论模型计算ϕ56mm聚能装药在联合加载结构下的射流成型,并通过理论计算分析联合加载时序对联合加载射流成型的影响。计算结果表明:电磁与炸药联合加载能够增加成型射流的速度以及动能,提高聚能射流对炸药的利用率,并且优化联合加载时序能够进一步提高射流速度和动能;联合加载下射流动能相对传统聚能装药提升了34.5%,对炸药的能量利用率提高了36%。

Abstract

To improve the explosive’s energy utilization rate of the shaped charge jet, a principle structure in which electromagnetic energy and explosive are combined to load the liner is proposed considering the principles of the liner loaded by electromagnetic force and magnetic flux compression generator. Combined loading is divided into two stages: in the first stage, the capacitor discharge and the initial magnetic flux is formed in the magnetic flux compression generator; in the second stage, Explosive 1 detonates, which compresses the magnetic field and increases the current, and at the same time, the detonation wave propagates to the liner to detonate Explosive 2. The electromagnetic force and detonation pressure jointly crush the liner. Combined with the explosive detonation theory, traditional PER theory of jet formaton and electromagnetic loading characteristics, a theoretical model of jet formation of a liner under combined electromagnetic and explosive loading is established. Then, the jet formation of ϕ56mm shaped charge under the combined loading structure is calculated by the theoretical model, and the influence of the combined loading sequence on jet forming is analyzed by theoretical calculation. The computational results show that the combined electromagnetic and explosive loading can increase the velocity and kinetic energy of the shaped jet, improve the explosive utilization rate of the jet, and optimize the combined loading sequence, which will further improve the jet velocity and kinetic energy. The kinetic energy of the jet under combined loading is 34.5% higher than that of the traditional shaped charge, and the energy utilization rate of explosive is increased by 36%.

关键词

聚能射流 / 电磁能与炸药 / 爆磁压缩发生器 / 能量利用率 / 联合加载时序

Key words

shaped jet / electromagnetic energy and explosives / energy utilization / magnetic flux compression generator / combined loading sequence

引用本文

导出引用
豆剑豪 , 贾鑫 , 梁争峰 , 黄正祥 , 薛标. 电磁能与炸药联合加载药型罩基础理论. 兵工学报. 2023, 44(6): 1754-1763 https://doi.org/10.12382/bgxb.2022.0120
DOU Jianhao , JIA Xin , LIANG Zhengfeng , HUANG Zhengxiang , XUE Biao. Basic Theory of Loading Liner With Electromagnetic Energy and Explosives. Acta Armamentarii. 2023, 44(6): 1754-1763 https://doi.org/10.12382/bgxb.2022.0120
中图分类号: TJ410.1 (基础理论)   

基金

国家自然科学基金青年基金项目(11602110)

0 引言

破甲弹通过聚能装药爆炸后压垮药型罩形成的金属聚能射流来对装甲目标进行毁伤[1]。为了提升聚能装药的威力,国内外学者进行了大量研究,徐文龙等[2]提出一种超聚能结构,并通过数值模拟和试验进行了验证,超聚能结构下射流的有效质量及侵彻能力远大于传统聚能射流。张昊等[3]和苏成海等[4]对活性材料射流侵彻混凝土进行了研究,以提高射流侵彻开孔和后效效应。Kemmoukhe等[5-6]对改变药型罩结构来提高聚能装药威力进行了研究。郑宇等[7]对双层药型罩的射流形成进行了研究分析。杨丽等[8]对带隔板装药的杆式射流成型及侵彻特性进行研究,以提高射流威力。但根据聚能装药的作用原理,如图1所示(β为压垮角),聚能装药在爆轰过程中爆轰产物向四周飞散,非毁伤元驱动方向的装药几乎没有利用,很难将较多的爆轰能量集中于某一特定方向,同时圆柱部分装药主要用于传递爆轰和调整波形,炸药能量没有加载至药型罩上,因此现阶段传统聚能装药形成的射流对炸药能量利用率较低。现有国内外公开资料表明,聚能装药对钢靶的最大侵彻威力为 10倍装药直径[9],很难再通过提高炸药爆速、爆压或改变药型罩结构及材料的手段来进一步提升聚能装药的侵彻威力。
图1 传统聚能装药原理

Fig.1 Schematic diagram of traditional shaped charge

Full size|PPT slide

爆磁压缩发生器(FCG)是一种能够把炸药化学能转化成电磁能的脉冲能源装置[10-11],其基本概念是通过炸药爆轰驱动电枢压缩定子产生磁场,将化学能转化为电磁能,从而在负载中实现脉冲电流和电磁能量的放大。自前苏联的A.D.Caxapov和美国的J.L.Folwer提出FCG以来[12],至今被国内外大量研究。国内中国流体物理研究所和国防科技大学自上世纪60年代开始研究[11],已于2000年研制的8-6型FCG可在5.2μH电感性负载上产生 100kA 的脉冲电流,储能达25kJ,脉宽17μs。孙奇志等[13]于2009年设计的EMG-125型螺线型发生器可以在25nH的负载上输出大于3MA的电流,电磁能量放大50倍。夏明等[14]利用FCG作为初始能源,成功加载薄壁金属管变形。可见FCG技术相对成熟,能够实现对炸药能量的转化。
本文基于Dou等[15-17]以及Grace等[18]对电磁加载药型罩的研究以及FCG原理,提出一种电磁与炸药联合加载药型罩的原理结构,如图2所示。图2中,S1S2为闭合开关,R1L1分别为传输线路的电阻和电感,C为电容器,I为传输线路电流。该结构将原聚能装药的圆柱部分炸药设计为FCG结构,该部分炸药一方面用于压缩磁场放大电流,实现炸药能向电磁能的转化,另一方面传递爆轰波,引爆炸药2,爆轰加载药型罩。通过上述方法,将炸药1的部分能量也用于加载药型罩,提高成型射流的速度和威力以及对炸药能量的利用率。本文建立了炸药与电磁共同加载药型罩的简单理论模型,并对ϕ56mm装药结构进行联合加载设计,计算了炸药和电磁共同加载下的射流成型,并通过理论分析了联合加载时序对射流成型的影响。
图2 炸药与电磁联合加载原理示意图

Fig.2 Schematic diagram of the combined loading principle of explosives and electromagnetic force

Full size|PPT slide

1 联合加载下射流成型理论模型

1.1 加载原理及基本假设

炸药与电磁联合加载的电路与原理示意图如图2所示,炸药1为原聚能装药的圆柱部分装药。电磁的初始能源由电容器C提供,左端虚线框内为FCG结构。加载共分为两个阶段:
1)第1阶段为初始磁通形成阶段,闭合放电开关S1,电容器对整个回路放电,电路包括传输线路电阻R1和电感L1,FCG的等效电阻R2和电感L2,药型罩装置(负载)的电阻R3和电感L3。此时电流较小,电磁力无法压垮药型罩。当电流到达电流峰值后,闭合开关S2,进入第2阶段。
2)第2阶段为磁通压缩阶段,药型罩在电磁力及炸药的加载下压垮成型。此时电容器被短路,由FCG向负载供能,在闭合开关S2的同时控制炸药起爆,起爆点位于圆柱部底端轴线处。炸药1爆轰后驱动电枢向外运动,根据磁场冻结效应,极短时间内定子线圈与电枢之间的磁通量几乎不变,空间被压缩导致磁感应强度增大,电流被放大,放大的电流在药型罩上产生垂直于药型罩表面指向里的电磁力。同时爆轰波在炸药1中向右传播,引爆炸药2,炸药2 爆炸后,爆轰压力作用药型罩上与电磁力共同加载药型罩向轴线压垮。当电枢完全膨胀,FCG的电感不再减小时,电流增加到最大值。存储在电感中的电磁能继续向负载放电,此时的电路为RL电路,电流不断减小,电磁力与炸药爆轰压力继续共同加载药型罩,驱动药型罩在轴线处碰撞形成射流,最后电路断开,联合加载结束。
炸药爆轰会产生破坏效应,在试验或应用中应考虑爆轰的破坏,可在电路中设置隔爆装置,降低爆轰波的在电路器件中的传播,在电路被破坏前完成加载。本文只进行初步的理论研究,因此在建立理论模型前,先做出一些基本假设:
1) 在磁爆加载过程中,电路、炸药、磁场和力场互相不影响。
2) 电路加载过程中,药型罩装置与FCG之间的磁场变化相互独立,互不影响。
3) 在电路作用过程中,除了理论模型设计的关键参数外,其余各元件为理想元件。
4) 炸药能够理想传爆及爆轰,爆轰波以平面波的方式传播。
5) 不考虑炸药爆轰对电路的破坏。

1.2 炸药爆轰压力

传统炸药加载下的药型罩压垮理论已经相对成熟,这里采用有效装药的绝热压缩方法计算炸药爆轰压力[1]。假设炸药瞬时爆轰,且稀疏波沿装药表面的内法线方向向爆炸产物内部传播,爆轰产物是以稀疏波初始交界面为刚性边界进行定向膨胀的,药型罩向轴线压垮运动时,有效装药部分向内绝热膨胀做功。有效装药可以用装药剖面图的各个角平分线来确定,如图3所示,MN平分∠PMQ,△MQN围成的区域即为该装药结构下的有效装药,故为了不减小原有装药结构下形成的射流威力,在采用图2所示结构时,炸药2的装药长度不能小于PM。
图3 有效装药示意图

Fig.3 Schematic diagram of effective charge

Full size|PPT slide

根据有效装药绝热膨胀的假设,作用于药型罩微元表面的压力 pi[1]应满足:
pi=p0i(V0iVi)γ
(1)
式中:p0i为爆轰产物的初始压力;V0i为药型罩微元对应的有效装药的初始体积;Vi为爆轰产物所占体积;γ为绝热系数,一般的猛炸药可取3。瞬时爆轰后,爆轰产物的初始压力可近似取为
p0i= 12pCJ= 18ρeD2
(2)
式中:pCJ为炸药的CJ压力;ρe为装药密度;D为爆速。

1.3 FCG理论模型

在加载第1阶段,加载电路相当于RLC电路,电流可用式(3)[15]表达:
I(t)= U0ωLe-ζtsin(ωt)
(3)
式中:U0为电容器的初始电压;ω为电路振荡频率,ω= 1LC-(R2L)2,L为电路总电感,R为电路总电阻;ζ为衰减系数,ζ=-R/2L。第1阶段的电流到达峰值时,闭合开关S2,此时的电流记为I0,FCG的初始电感为L20,存储在FCG中的总磁通量为
ψ=I0L20=NΦ=IL2
(4)
式中:N为线圈匝数;Φ为单匝线圈的磁通量; L2[19]可表达为
L20πN2(rc2- ra2)f(l,rc,ra)K
(5)
μ0为真空磁导率,rc为定子线圈半径,ra为电枢半径,f为电感修正系数,
f= 1011l+8(rc-ra),l<rc-ra1011l+9(rc-ra),lrc-ra
(6)
l为线圈长度;K为计算动态电感时的体积修正系数,
K=1- ViV0i
(7)
假定磁感应强度在空间内均匀分布[13],因此空间内的磁感应强度B为
B= IL2NS
(8)
式中:S为定子线圈与电枢间空腔的剖面积。
在磁通压缩的第2阶段,炸药爆炸,电枢在爆轰压力作用下向外膨胀,此时作用在电枢上的有向外的爆轰压力以及向内的磁压,磁压力模型[20]
py= B22μ0cos θ
(9)
式中:θ为电枢膨胀角,即电枢向外膨胀时与轴线的夹角,如图4所示。在炸药爆轰的初始时刻,θ=0°。
图4 FCG作用过程

Fig.4 Process of FCG

Full size|PPT slide

电枢在爆轰压力和磁压共同作用下的运动方程为
mi dvidt=si(pdi-pyi)
(10)
式中:mi、vi和si分别为第i个电枢微元的质量、速度和面积;pdi和pyi分别为第i个电枢微元受到的爆轰压力和电磁力。
磁通压缩的第2阶段,开关S2闭合,电容器被短路,FCG作为供能源,此时的电路方程为
αI dL2dt+(αL2+L3) dIdt+RI=0
(11)
式中:α为磁通损耗系数,根据试验结果来拟合,一般取0.7~0.97。由于药型罩装置的电感远小于FCG电感,为了简化模型,不考虑药型罩的电感变化,且利用修正系数来近似替代电阻的热效应,动态电阻Rc=cR,c为修正系数[21]。将式(11)积分,即可得磁通压缩阶段的电流表达式:
I(t)= I0(αL20+L3)αL2(t)+L3exp(- 0tRcαL2(t')+L3dt')
(12)

1.4 压垮速度计算

记炸药1的长度为l1,炸药2的长度为l2,中间的间隔为l3,炸药1开始爆轰的时刻ts=(l1+l3)/D,此时药型罩可能已经在电磁力的作用下开始压垮,假定爆轰波速度无衰减,爆轰波追赶上压垮微元时,微元的运动距离为h1,此时对应的时刻即为爆轰压力与电磁力共同压垮药型罩的时刻:
te=ts+ x+h1D
(13)
式中:x为该微元距离顶端微元的轴向距离。药型罩微元的运动方程为
M dvdt= Fe-Fs,t<teFe+spd-Fs,tte
(14)
式中:M为药型罩微元的质量;v为压垮速度;Fe为电磁力,Fe0I2/4πro,ro为微元外径;Fs为塑性力,认为微元压垮在力的作用下开始压垮后Fs为定值,Fs= μ0Is24πr2,Is为微元的临界屈服电流[17];s为罩微元的表面积;pd为炸药爆轰压力。
将式(14)对时间积分,即可求得轴线处压垮速度v0。由于加载过程中电流与药型罩形态速度随时间而变化,不能通过直接积分求得压垮速度,需要利用压垮时各参数(电感L2,电流I,内径ri,外径ro等)对时间t的单值性,采用迭代的方法实现压垮速度的计算。单值是指计算中任意时刻,所有参数都只有唯一值。计算中以初始加载电压、初始电路参数为初始条件,计算得到第1阶段的峰值电流,再以电流I0、初始FCG参数及药型罩内外径为初始条件,得到微元在加载第2阶段和第3阶段的联合加载下的运动速度及位置等结果,具体的计算流程图如图5所示。
图5 计算联合加载下药型罩压垮速度流程图

Fig.5 Flow chart of calculating the collapse velocity of the liner under combined loading

Full size|PPT slide

图5中,CURRENT1、CURRENT2、EXPLOSIVEVELOCITY为自编的计算程序,根据式(1)~式(14)编写,输入初始的电压U0、电阻R0以及电感L0后,程序CURRENT1计算出加载第1阶段的峰值电流I0,输入至CURRENT2中,并且以初始FCG参数为初始输入条件,计算得到t时刻的电流值,进而求得作用在微元上的电磁力Fe。同时判断爆轰波是否传递至罩微元,若t>te,则由EXPLOSIVE计算出作用在罩微元上的爆轰压力pi,与电磁力Fe共同输入至VELOCITY程序中,求解t时刻微元新的位置(ri, ro)以及此时的压垮速度v0,反之则电磁力Fe单独加载药型罩。接着将得到的新的电阻电感以及微元位置反馈到下一步的计算中,作为下一个时间步的输入条件,依次循环直到计算结束。循环是否结束的判断依据是微元是否运动至轴线处(ri>0m)且压垮速度是否大于0m/s (v0>0m/s)。当ri≥0m 时,微元正常压垮至轴线处发生碰撞,输出此时刻的压垮速度即为在轴线处碰撞的压垮速度v0;当v0≤0m/s时,说明电磁力太小不能克服材料的屈服应力导致运动停止,计算结束,输出压垮速度 v0=0m/s, 以及此刻的微元位置。

1.5 射流成型理论

将坐标系建立在碰撞点Q上,从驻点坐标系来看药型罩的压垮过程,可以看见药型罩沿着母线方向向轴线压垮,记地球坐标系下的压垮速度为v0,驻点移动速度为v1,驻点坐标系下的压垮速度为v2,v0、v1和v2的关系如图6所示。图6中η为联合加载下药型罩微元的抛射角,可由其定义来计算。
图6 驻点坐标系下的速度关系

Fig.6 Velocity diagram in stationary point coordinate system

Full size|PPT slide

展示的速度之间的夹角由几何关系得到,各个速度之间的关系由正弦定理得到:
v1= cos(η+α-β)sinβv0v2= cos(α+η)sinβv0
(15)
根据伯努利方程以及电磁加载特性[16],联合加载下射流和杵体的表达式为
vj=v1+v2b= v0sinβ(cos(η+α-β)+λcos(α+η))vs=v1-v2a= v0sinβ(cos(η+α-β)-cos(α+η))
(16)
式中:λ= 2p2/ρv22+1,p2为在碰撞点作用在药型罩上的压力,p20I2/8π2 r22

2 联合加载下药型罩形成射流分析

以ϕ56mm装药为例,对联合加载进行计算分析。将原装药结构用于联合加载,ϕ56mm装药结构如图7[21]所示,外径为56mm,装药高度73.3mm,药型罩壁厚为1mm,锥角为60°,所使用的炸药为无壳8701,装药质量203g,装药密度1.713g/cm3,其爆速为7980m/s
图7 ϕ56mm聚能装药结构[21]

Fig.7 ϕ56mm shaped charge structure[21]

Full size|PPT slide

联合加载结构如图8所示,将基准装药分为炸药1和炸药2两部分,装药直径保持不变,都为56mm,对应长度l1=20mm、l2=53.3mm,药型罩结构保持一致,电枢直径ra=31mm,厚度为3mm,定子线圈的缠绕半径rc=50mm,选取导线直径 ds=2mm,因此线圈的匝数N=l1/ds=10。给定初始电路参数为加载电压U=20kV,电容C=50μF,连接导线电阻R1=0.1,电感L1=0.1μH,经计算FCG的电阻和电感分别为R2=2.2、L2=4.1μH,负载的电阻和电感分别为R3=0.01、L1=9.4nH。可以看出,相较于连接导线,FCG的电阻电感较大,负载的电阻电感可以忽略不计,在加载条件不变情况下,FCG的结构是影响加载电流的关键。
图8 联合加载装药结构

Fig.8 Charge structure under combined loading

Full size|PPT slide

动态电阻的修正系数取5.7[21],忽略磁通损耗,在上述加载结构以及条件下,第1阶段加载的电流峰值为66.5kA。到达电流峰值后,闭合开关S2短路电容器,同时引爆炸药,开始爆磁压缩阶段。记此时的时刻为0μs时刻,磁通压缩阶段的电流以及FCG的电感变化如图9所示。在炸药爆轰的初始阶段,电枢的膨胀速度较低,电感变化较慢且减小的幅度较小,故炸药能量对电磁能的转化较慢,在10μs前电流增加缓慢。随着FCG电感不断减小,电流增加速率变快,在14μs时,FCG电感减小到最小值时,电流增大到最大值,图9中所示电流峰值为 1739.6kA。 理想情况下,FCG电感减小到无限接近于0μH,随后开始加载第3阶段,FCG中存储的电磁能对电路放电,电流不断减小,直到加载结束。
图9 磁通压缩阶段电流以及FCG电感变化

Fig.9 Current and FCG inductance changes in the magnetic flux compression stage

Full size|PPT slide

取药型罩微元长度为1mm,共42个微元,由药型罩的顶端至底端依次编号为1~42,药型罩在炸药加载和电磁炸药联合加载下的压垮速度和形成的射流速度分别如图10图11所示,药型罩微元在联合加载下压垮至轴线的时间如图12所示。
图10 不同加载方式下压垮速度对比

Fig.10 Comparison of collapse velocities of different loading methods

Full size|PPT slide

图11 不同加载方式下射流速度对比

Fig.11 Comparison of jet velocities of different loading methods

Full size|PPT slide

图12 微元压垮至轴线时间

Fig.12 Time of elements collapsing to the axis

Full size|PPT slide

两种加载方式下的头部压垮速度和射流速度的变化趋势几乎一致,由图12所示,药型罩头部微元压垮至轴线的时刻在10μs之前,而由图9可见,在10μs之前电流较小,在压垮过程中作用在微元上的电磁力远小于爆轰力,电磁能未能完全作用于微元上,电磁力对头部微元的压垮速度的影响可忽略不计,故联合加载下头部微元的压垮速度和射流速度与炸药加载几乎一致。对于较靠后的微元,在压垮过程中,电磁力较大,在微元压垮至轴线前,大部分电磁能已经作用在微元上,故联合加载下的微元压垮速度和射流速度要高于炸药加载。
经计算,炸药加载下药型罩形成射流的动能为116kJ,8701炸药单位体积的能量为9.5J/mm3,故炸药总能量为1373kJ,对炸药能量的利用率为8.4%。在上述结构的联合加载下,射流动能为139kJ, 初始电容器的能量为10kJ,相较于炸药能量可以忽略不计,因此联合加载下射流动能提高了19.8%,对炸药的能量利用率为10.1%,相对提高了20.2%。

3 联合加载时序对射流成型的影响

电磁能与炸药的联合加载是涉及电场、磁场、炸药爆轰的复杂过程,药型罩结构、装药结构、初始加载能量、加载时序等因素都会影响射流的成型状态,为研究联合加载下射流的成型机制,需要对各影响因素进行分析,这里选择对联合加载时序对射流成型的影响进行分析。
在上述计算中,由于药型罩头部附近微元的压垮时间与电流峰值时间未联合,头部微元速度相较于传统加载无明显提升。为了提高联合加载下药型罩头部微元的压垮速度以及射流速度,需要对联合加载时序进行调整,本文提出两种方案,一是通过改变联合加载结构,减小电流上升时间,二是将炸药2延时爆轰。下面分别对这两种优化方案下的射流成型进行计算。

3.1 减小电流上升时间

由式(11)可知,电流上升速率与FCG电感变化相关,而电感变化速率,取决于炸药能量、爆速、电枢厚度以及电枢和线圈的半径,这里通过只改变线圈缠绕半径,即rc的大小来改变电流上升时间。计算算例以及电流计算结果如表1所示,表1中下标0表示加载0μs时刻的值,tu为电流上升时间,Imax为电流峰值,加载阶段的电流如图13所示。随着线圈缠绕半径的减小,FCG的初始电阻和电感都相对减小,在电容器电压和电容不变的情况下,第1阶段的电流峰值会随之增大。由于定子线圈的半径减小,电枢膨胀时间缩短,因此磁通压缩阶段的电流上升时间减小,但由于FCG初始电感的减小,导致FCG的电感在电路总电感中比重较小,因此放大后的电流峰值会相对减小。
表1 计算算例

Table 1 Numerical examples

编号 rc/mm R20/ L20/μH I0/kA tu/μs Imax/kA
s-1 50 2.2 4.1 66.5 14.3 1737
s-2 45 2.0 3.2 76.1 12.1 1585
s-3 40 1.8 2.2 91.5 9.8 1357
图13 不同算例的电流比较

Fig.13 Current comparison of different examples

Full size|PPT slide

随后计算了微元的压垮速度和射流速度,其中压垮速度以及第10、第20、第30和第40个微元的压垮速度对比如图14图15所示,算例s-0为炸药加载。由图14可以看出,3种算例下,联合加载下药型罩微元的头部的压垮速度与炸药加载几乎一致,而随着微元位置的后移,结合图15,在第10个微元处,算例s-3的压垮速度最先开始有明显的增大,随后在第15个微元处,算例s-2的压垮速度有明显增大。对于第20个微元,算例s-3的压垮速度最大而s-1的压垮速度最小,说明了随着电流上升时间的减小,电磁力对药型罩头部微元的加载作用增大,在微元压垮至轴线前,更多的电磁能作用于微元上,相较于算例s-0,电磁力使得更多药型罩微元的压垮速度增大。对于后续微元,如图15所示的第30 和第40个微元,压垮速度随着电流上升时间的减小先增加后减小,这是因为底端微元压垮至轴线前,绝大部分电磁能已经作用,而电流峰值的降低,导致作用于底端微元的电磁能量减少,故压垮速度降低。
图14 不同算例的压垮速度对比

Fig.14 Comparison of collapse velocities of different examples

Full size|PPT slide

图15 不同微元的压垮速度对比

Fig.15 Comparison of collapse velocities of different elements

Full size|PPT slide

图16展示了各微元形成的射流速度,图17为第10、第20、第30和第40个微元的射流速度对比,可以看出电流上升时间对射流速度的影响规律与压垮速度相同,由于压垮角的变化,射流速度对微元位置的变化相较于压垮速度出现了延迟,如第10个微元时,算例s-3的压垮速度相比于其他算例已经有明显的增大,但射流速度仍与其他算例几乎相同,在第13个微元时,才出现明显的增大。
图16 不同算例的射流速度对比

Fig.16 Comparison of jet velocities in different examples

Full size|PPT slide

图17 不同微元的射流速度对比

Fig.17 Comparison of jet velocities of different elements

Full size|PPT slide

算例s-2和s-3计算下的射流动能分别为 143kJ 和144kJ,3种算例下形成的射流动能都大于炸药加载下的射流动能,虽然随着电流上升时间的减小,药型罩顶端微元的速度有显著增加,但同时因为电流峰值的减小,相较于算例s-1,算例s-2和s-3的射流动能并没有显著增加。

3.2 炸药延时引爆

以算例s-1为基准,假定炸药2能延时起爆,由表1可知,算例s-2和算例s-3的电流上升时间相较与算例s-1分别减少了2.2μs和4.5μs,作为对比,假定炸药2分别相对于算例s-1延时起爆2.2μs和 4.5μs,记为s-2-1和s-3-1,相当于电流上升时间减小而电流峰值不变。计算得到的药型罩微元压垮速度和射流速度如图18图19所示。算例s-2-1和s-3-1的药型罩的压垮速度和射流速度随微元的变化趋势与其对应的算例s-2和s-3一致,但数值上较大,说明了炸药2延时起爆后,电流集中的时间区间与药型罩顶端微元的压垮时间吻合较好,使得电磁能更多地作用于靠近药型罩顶端的微元。加载算例s-2-1和s-3-1的电流峰值分别大于算例s-2和s-3,所以微元的压垮速度和射流速度均增大。
图18 压垮速度对比

Fig.18 Comparison of collapse velocities

Full size|PPT slide

图19 射流速度对比

Fig.19 Comparison of jet velocities

Full size|PPT slide

表2统计了所有算例下的射流动能以及对炸药的能量利用率,可以看出,相同炸药能量下,联合加载下形成的射流速度和动能都高于炸药加载,其中算例s-3-1,通过炸药延时起爆的方法,使得电流作用时间与炸药爆轰的作用时间联合更好,加载药型罩得到的射流动能最大(156kJ),相对于炸药加载,射流动能提高了34.5%,对炸药的能量利用率相对提高了36%。
表2 射流动能和能量利用率

Table 2 Kinetic energy and energy utilization rate of jet

算例 射流动能/kJ 炸药能量利用率/%
s-0 116 8.4
s-1 139 10.1
s-2 143 10.4
s-3 144 10.5
s-2-1 147 10.7
s-3-1 156 11.4

4 结论

本文基于FCG原理,将原聚能装药结构中的圆柱部分装药设计为FCG结构,炸药1能量转换为电磁能后,与炸药2共同加载药型罩,同时建立了联合加载药型罩的理论模型,并通过理论计算了ϕ56mm装药结构下联合加载的射流成型。得到主要结论如下:
1) 建立了电磁能与炸药联合加载药型罩的理论模型,理论计算结果表明,联合加载能够提高部分药型罩的压垮速度和射流速度,传统ϕ56mm聚能装药形成的射流动能为116kJ,联合加载下形成的射流动能为139kJ,提高了19.8%,对炸药的能量利用率提高了20.2%。
2) 电磁能与炸药能量作用时序对射流成型有较大影响,在本文算例下,当电流的上升时间减小或者炸药爆轰作用时间延后时,联合加载下药型罩头部微元的压垮速度和射流速度增加,电磁能作用更多的药型罩微元,药型罩形成的射流动能相应增加。在本文算例中,通过改变电流与炸药的联合作用时序,形成射流动能最大为156kJ,相较于传统炸药加载,提高了34.5%,对炸药的能量利用率提高了36%。

参考文献

[1]
黄正祥. 聚能装药理论与实践[M]. 北京: 北京理工大学出版社, 2014.
HUANG Z X. Theory and practice of shaped charge[M]. Beijing: Beijing University of Technology Press, 2014. (in Chinese)
[2]
徐文龙, 王成, 徐斌. 超聚能射流形成过程机理研究[J]. 兵工学报, 2018, 39(2):261-268.
摘要
提出一种超聚能射流形成的理论计算方法,给出了附加装置材料密度、厚度和药型罩锥角、密度、厚度与超聚能射流速度、有效质量之间的关系。基于非线性显式动力学分析软件AUTODYN,进行了超聚能射流形成过程数值模拟,其结果与理论计算结果基本一致。通过超聚能射流侵彻试验,证明试验结果与数值模拟结果吻合较好。研究表明:在一定范围内,随着药型罩锥角的增大,超聚能射流速度逐渐减小、质量逐渐增大;超聚能射流速度、有效质量及侵彻能力远大于传统聚能射流。
XU W L, WANG C, XU B. Investigation of hyper shaped charge jet formation theory[J]. Acta Armamentarii, 2018, 39(2):261-268. (in Chinese)
[3]
张昊, 王海福, 余庆波, 等. 活性射流侵彻钢筋混凝土靶后效超压特性[J]. 兵工学报, 2019, 40(7): 1365-1372.
摘要
为研究活性射流侵彻钢筋混凝土靶后效超压效应,采用实验、理论与数值仿真相结合的方法,开展了活性药型罩聚能装药作用钢筋混凝土靶后效超压特性研究。实验结果表明,在聚能装药结构一定条件下,低密度活性射流较高密度活性射流对钢筋混凝土靶的侵彻孔径更大,穿靶后在密闭空间内爆燃产生的冲击波峰值超压更高,且超压-时间曲线呈现出多峰现象。引入虚拟爆炸点的方法,建立了活性射流靶后爆燃超压与高能炸药爆炸等效分析模型,结合AUTODYN-3D有限元分析,得出模型预测的后效超压-时间曲线与实验结果吻合较好,揭示了活性射流在靶后密闭空间内爆燃超压传播规律及多峰现象。
ZHANG H, WANG H F, YU Q B, et al. Aftereffect overpressure of reactive jet perforating into reinforced concrete[J]. Acta Armamentarii, 2019, 40(7): 1365-1372. (in Chinese)
In order to investigate the aftereffect overpressure characteristics of reactive jet penetrating into reinforced concrete, the experimental, theoretical and numerical simulation methods are used to study the aftereffect overpressure characteristics of reactive liner shaped charge against reinforced concrete. Experimental results show that, for the same shaped charge configuration, the low density reactive jet produces a larger diameter penetration hole on the reinforced concrete and a higher peak overpressure in a closed space compared with the high density reactive jet, and the overpressure-time curves reveal multi-peaks phenomenon. An equivalent model for aftereffect deflagration overpressure of reactive jet and high explosive explosion is developed by using the virtual explosion point method.AUTODYN-3D finite element analysis shows that the predicted aftereffect overpressure-time curves are well agreement with the experimental results, and the propagation characteristics and multi-peaks phenomenon of reactive jet in the closed space are also revealed. Key
[4]
苏成海, 王海福, 谢剑文, 等. 活性射流作用混凝土靶侵彻与爆炸效应研究[J]. 兵工学报, 2019, 40(9): 1829-1835.
摘要
为分析活性药型罩配方、炸高对混凝土靶毁伤威力影响规律,开展了活性射流作用混凝土靶侵彻与爆炸联合毁伤效应研究。采用实验和数值模拟相结合的方法,对活性射流作用混凝土靶的典型毁伤模式进行探究,给出了配方与炸高对毁伤效应的影响特性。实验结果表明:活性射流作用下,混凝土靶呈现为显著的锥形爆坑和裂纹毁伤效应;气体产物量较高配方的活性射流对混凝土靶产生更强的毁伤效应;在1倍装药直径炸高下,炸坑直径可达10倍装药直径以上。基于有限元分析软件AUTODYN-3D平台开展了活性射流对混凝土靶侵彻与爆炸行为的数值模拟,揭示了炸高对毁伤效应的影响机理:随着炸高的增大,进入侵坑内部的活性材料随之减少,活性射流对混凝土的侵彻深度呈现先增大、后减小的趋势;当炸高为1倍装药直径时,活性射流动能侵彻与爆炸反应延迟匹配较好,侵彻与爆炸联合毁伤威力较强。数值模拟结果与实验结果吻合较好。
SU C H, WANG H F, XIE J W, et al. Penetration and damage effects of reactive material jet against concrete target[J]. Acta Armamentarii, 2019, 40(7): 1829-1835. (in Chinese)
[5]
KEMMOUKHE H, SAVI S, TERZI S, et al. Improvement of the shaped charge jet penetration capability by modifying the liner form using AUTODYN-2D[J]. Scientific Technical Review, 2019, 69(1):10-15.
[6]
KEMMOUKHE H, BURZIC Z, SAVIC S, et al. Improvement of shaped charge penetration capability and disturbation of the jet by explosive reactive armor[J]. Technical Gazette, 2019, 26(6):1658-1663.
[7]
郑宇, 王晓鸣, 李文彬, 等. 双层药型罩射流形成的理论建模与分析[J]. 火炸药学报, 2008, 31(3):10-14.
摘要
为研究乳化炸药基质含水量对其热分解特性和化学动力学参数的影响,制备了含水量分别为3.46%和12.27%的乳化基质。用DSCTG联用仪得到两试样在不同加热速率下的DSC、TGDTG图谱。通过对比外推起始分解温度、组成及图谱,研究了水含量对其热稳定性和热分解特性的影响。结果表明,无论含水率多少,在被加热初期,基质失重缓慢平稳;失重的主要原因是失水,这些水主要是游离于基质中的水和被加热时少量乳化微粒破乳后释放的水;低含水量基质外推起始分解温度明显低于高含水量。分解开始后,含水量高的乳化炸药基质放热速度和失
ZHENG Y, WANG X M, LI W B. Theoretical modeling and analysis on jet formation of double-layered conical liner[J]. Chinese Journal of Explosives & Propellants, 2008, 31(3):10-14. (in Chinese)
[8]
杨丽, 陈闯, 张健, 等. 带隔板装药的杆式射流成型试验及侵彻特性分析[J]. 兵工学报, 2016, 37(4):48-53.
YANG L, CHEN C, ZHANG J, et al. Analysis of formation experiment and penetration property of rod-like jet of shaped charge with foam[J]. Acta Armamentarii, 2016, 37(4):48-53. (in Chinese)
[9]
NAEEM K, HUSSAIN A, ABBAS S. A review of shaped charge variables for its optimum performance[J]. Engineering, Technology & Applied Science Research, 2019, 9(6): 4917-4924.
[10]
杜枢, 孙奇志, 刘伟. 圆盘型爆磁压缩发生器的数值模拟[J]. 强激光与粒子束, 2016, 28(8):085001.
DU S, SUN Z Q, LIU W. Numerical simulation of disk explosive magnetic generator[J]. High Power Laser and Particle Beam, 2016, 28(8): 085001. (in Chinese)
[11]
陈冬群, 曹胜光, 刘永贵, 等. 爆磁压缩发生器试验技术[J]. 实验技术与管理, 2002, 19(3):11-14.
CHEN D Q, CAO S G, LIU Y G, et al. Test technology of explosive magnetic compression generator[J]. Experimental Technology and Management, 2002, 19(3): 11-14. (in Chinese)
[12]
龚兴根, 孙奇志, 刘正芬, 等. 爆磁压缩技术综述[J]. 爆轰波与冲击波, 2003(3): 130-138.
GONG X G, SUN Q Z, LIU Z F, et al. Overview of explosive magnetic compression technology[J]. Detonation Wave and Shock Wave, 2003(3): 130-138. (in Chinese)
[13]
孙奇志, 孙承纬. 轴线起爆式螺线管型爆磁压缩发生器理论模型[J]. 强激光与粒子束, 2003, 15(4):385-390.
SUN Q Z, SUN C W. Theoretical model of axis initiating solenoid explosive magnetic compression generator[J]. Intense Laser and Particle Beam, 2003, 15(4): 385-390. (in Chinese)
[14]
夏明, 黄正祥, 顾晓辉, 等. 磁爆加载薄壁金属管的冲击变形实验研究[J]. 兵工学报, 2013, 34(3):257-262.
摘要
利用构建的磁爆加载实验平台,对脉冲磁动力与电爆炸耦合效应作用下薄壁金属管的冲击变形进行了研究,获得了加载过程的强流放电波形,以及不同结构和材料薄壁金属管的冲击变形。结果表明:以脉冲电容器组为能源时,RLC 电路过程能较好地反映加载时的强流放电特性;相同结构薄壁金属管的冲击变形主要受加载方式的影响,加载过程产生的焦耳热使材料显著升温,更易于变形,但高温将导致材料气化,无法形成磁聚毁伤元;磁聚毁伤元不宜采用低熔点或低延展性材料,相比铝材,紫铜更优,而截锥管比圆管结构更易于控制形成磁聚毁伤元。
XIA M, HUANG Z X, GU X H, et al. Experimental study on impact deformation of thin-walled metal tube loaded by magnetic explosion[J]. Acta Armamentarii, 2013, 34(3): 257-262. (in Chinese)
[15]
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.
[16]
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.
The use of a shaped liner driven by electromagnetic force is a new means of forming jets. To study the mechanism of jet formation driven by electromagnetic force, we considered the current skin effect and the characteristics of electromagnetic loading and established a coupling model of “Electric–Magnetic–Force” and the theoretical model of jet formation under electromagnetic force. The jet formation and penetration of conical and trumpet liners have been calculated. Then, a numerical simulation of liner collapse under electromagnetic force, jet generation, and the stretching motion were performed using an ANSYS multiphysics processor. The calculated jet velocity, jet shape, and depth of penetration were consistent with the experimental results, with a relative error of less than 10%. In addition, we calculated the jet formation of different curvature trumpet liners driven by the same loading condition and obtained the influence rule of the curvature of the liner on jet formation. Results show that the theoretical model and the ANSYS multiphysics numerical method can effectively calculate the jet formation of liners driven by electromagnetic force, and in a certain range, the greater the curvature of the liner is, the greater the jet velocity is. © 2020 The Authors
[17]
DOU J H, JIA X, HUANG Z X, et al. Theoretical and experimental study on jet formation and penetration of the liner loaded by electromagnetic force[J]. International Journal of Mechanical Sciences, 2020, 186:105883.
[18]
GRACE F, DEGNAN J, ROTH C, et al. Shaped charge jets driven by electromagnetic energy[C]//Proceedings of the 28th International Symposium of Ballistic. Atlanta, GE, US: International Ballistics Society, 2014.
[19]
APPELGREN P, BRENNING N, HURTIG T, et al. Modeling of a small helical magnetic flux-compression generator[J]. IEEE Transactions on Plasma Science, 2008, 36(5):2662-2672.
[20]
ALTGILBERS L L, BROWN M D J, GRISHNAEV I, 等. 磁通量压缩发生器[M]. 孙承纬, 周之奎, 译. 北京: 国防工业出版社, 2008.
ALTGILBERS L L, BROWN M D J, GRISHNAEV I, et al. Magnetocumulative generators[M]. SUN C W, ZHOU Z K, translated. Beijing: National Defense Industry Press, 2008. (in Chinese)
[21]
马彬. 强磁场与聚能装药耦合特性研究[D]. 南京: 南京理工大学, 2019.
MA B. Study on couplingcharacteristics of strong magnetic field and shaped charge[D]. Nanjing: Nanjing University of Science and Technology, 2019. (in Chinese)

文章所在专题

特种车辆理论与技术

61

Accesses

0

Citation

Detail

段落导航
相关文章

/