A simple method of measuring impulse current of small high-voltage exploding device
-
摘要: 依据脉冲放电电路的等效电路及其微分方程,采用Levenberg-Marquarat算法对主电容放电电压测试波形数据进行衰减系数识别,从而获得模拟电流波形。该方法克服了分流器法和Rogowski线圈法等直接测量小型高压引爆装置冲击电流时,因附加电路引起的电流波形失真。MATLAB模拟结果表明,该方法得到的电流模拟波形与真实电流波形拟合精度高,可用于直列式引信电子安全与解除保险装置和低能冲击片雷管的优化匹配设计。
-
关键词:
- 爆炸力学 /
- Levenberg-Marquarat算法 /
- 放电电流模拟 /
- 小型高压引爆装置 /
- 衰减系数
Abstract: A simple method was developed for measuring the impulse current waveform produce by a small high-voltage exploding device. The model of impulse current was determined by attenuation coefficient. To estimate the attenuation coefficient from the actual discharge voltage curve, the Levenberg-Marquarat algorithm was applied based on the equivalent circuit of discharge circuit and its differential equations. Compared with the direct measuring methods such as using shunt or Rogowski coil, this method overcomes the distortion of impulse current waveform caused by additional measuring circuit. The results show that the simulation current waveform fits the actual current waveform well. The method can be used for optimum matching design of electronic safety, arming device of in-line fuse or low-energy slapper detonator. -
圆筒试验是用于评估炸药作功能力及确定爆轰产物状态方程的常用标准试验,为了确保试验数据的有效性或获取更充分的试验信息,一般要求铜管直径增至初始直径的3倍时,其壁面不能发生贯穿性破裂[1],这对铜管的材质性能及加工工艺均提出了较高的要求。铜管的材料一般选用经退火处理的无氧铜,其中,无氧铜的晶粒度作为铜管延展性能的重要影响因素之一,在材质选择和处理过程中受到密切关注,LANL实验室曾针对
∅ 25 mm标准铜管建议平均晶粒尺寸为25~30 μm。目前,无氧铜的晶粒细化主要是通过对棒料进行二次锻造实现,工艺较为复杂,尤其是对于大尺寸铜管所需的棒料,其加工难度较大,有时存在明显的粗晶或晶粒不均匀现象,那么这种无氧铜加工成的铜管是否能满足圆筒试验要求,还需要从两个方面进行分析,即晶粒度对铜管断裂过程和炸药驱动性能表征参量的影响。粗晶或晶粒不均匀一般会导致材料拉伸过程中的损伤增加[2-3],更易形成裂纹或断裂,但在铜管的膨胀过程中,材料性能与爆炸加载应力存在相互耦合作用,并共同决定着管壁的断裂表现[4],因此,判断这种铜管的断裂直径是否能达到初始直径的3倍,还需要相应的试验验证。而目前关于这方面的试验研究主要针对钢、铝合金等材质的管状结构[5-6]或不满足圆筒试验相似准则的铜管结构[7-10],其试验结果难以借鉴。另一方面,无氧铜的晶粒尺寸增大时,高应变率下的流动应力及应变能有时相对偏低[11],可能会造成炸药驱动性能表征参量(如,比动能、格尼速度等)发生明显的变化,但是否会对圆筒试验结果造成显著影响,还需要进行对比试验验证。
鉴于此,本研究采用两种具有不同晶粒度的软态无氧铜,加工成
∅ 25 mm圆筒试验用的标准铜管,通过高速扫描方法分析JO-159炸药爆轰加载下两种铜管的膨胀位移、比动能等曲线的差异;采用分幅相机观察铜管的断裂过程,分析两种铜管断裂直径、断裂方式等方面差异,以期能为圆筒试验用铜管的材质选择及加工工艺改进提供参考。1. 试 样
1.1 无氧铜管
采用中铝洛阳铜业有限公司生产的一号无氧铜TU1(GB5231-2001)棒料作为原材料,该材料经真空退火处理,其密度为8.94 g/cm3,其金相结构如图1(a)所示,可以看出,该无氧铜材料晶粒较粗,大部分的晶粒尺寸处于100~300 μm范围内,且形状不规则。作为对比,对该类棒料进行二次锻造加工,并再次进行退火处理,其材料的金相结构如图1(b)所示,其平均晶粒度为20~30 μm。将这两种棒料加工成
∅ 25 mm圆筒试验用标准铜管,其内径为25 mm,外径为30 mm,长度为300 mm,通过工业CT扫描,均未发现明显裂纹或缺陷。1.2 炸药
由于一般认为金属材料的断裂应变等参量受应变率影响较大,而在铜管膨胀过程中,其膨胀速率越高,则应变率越大,因此,本研究中选择在铜管内装填驱动能力较强的JO-159炸药[12]。该炸药的成分中HMX的质量分数为95%,将其压制成12个
∅ 25 mm×25 mm的药柱,密度为1.835 g/cm3。试验前,将药柱依次装入铜管后,在铜管的两个端面各粘接一个电探针,以便在试验中获取炸药在铜管中的平均爆速,然后在起爆端粘接一个∅ 25 mm起爆药柱。2. 铜管膨胀过程的测量及分析
2.1 两种铜管膨胀曲线的差异
表征炸药驱动性能的比动能、格尼速度等参量主要是根据铜管膨胀位移曲线计算所获,因此,这里首先通过狭缝扫描试验对比两种不同晶粒度无氧铜管的膨胀轨迹。试验主要由狭缝观测铜管某一横截面的直径变化,其布局如图2所示,对于
∅ 25 mm圆筒试验,则该横截面距离起爆端200 mm,采用SJZ-15型转镜式高速扫描相机进行观测,扫描速度设定为3 mm/μs。图3为粗晶无氧铜管膨胀过程的扫描试验底片,可通过判读黑白分界线的坐标获得铜管外表面的径向位移随时间的变化历程,其中,横坐标与相机扫描速度的比值即为时间值,纵坐标则对应着铜管外壁的径向位移。此外,根据扫描试验的底片有时也能看出管壁的断裂位置,这主要是由于溢出的产物具有较高的速度,使得图像的边界迹线出现一定程度的突跃,尤其是在管壁膨胀的早期和中期,该突跃现象尤为明显[7]。图3中的边界迹线较为光滑,这说明该时间段铜管断裂的可能性较小,因此其试验数据可认为是有效的。
通过对两种铜管的扫描试验底片进行数字化判读,可获得铜管外表面的径向位移
Δre 随时间t的变化曲线,如图4所示,可以看出,这两种铜管膨胀曲线的整体差异较小。为了深入对比铜管比动能等参量的差异,这里暂不讨论铜管的具体断裂位置,可先根据圆筒试验数据处理方法[13-14]进行分析。该数据处理方法的主要流程包括:(1)假定铜管的横截面积不变,将Δre-t 数据转换为铜管质量中心面的径向位移数据Δrm-t ,并按照特定的公式进行拟合处理;(2)计算铜管质量中心面的径向速度um=d(Δrm)/dt ;(3)结合铜管壁不同方向运动速度的几何关系,获得其质点速度us ,然后进一步计算出表征炸药驱动性能的铜管比动能E=u2s/2 等参量。图5为两种铜管比动能的变化曲线,可以看出,粗晶铜管的比动能相比细晶铜管略高。
Δre =6.0,12.5,19.0 mm三个位置(分别对应爆轰产物相对比容V=2.2,4.4,7.0的状态,可分别描述高压、中压及低压作用阶段的特征)的数据显示,在这三处,两种铜管的比动能相对偏差分别为0.64%、0.52%、0.55%,均处于较小的范围内。而对于炸药格尼速度√2Eg 的差异,结合装填比和V=7.0时的us 值,可计算出由细晶和粗晶铜管所获的该炸药格尼速度值分别为2.940和2.948 mm/μs。虽然该试验结果反映出相同炸药爆轰加载下这两种铜管的比动能及格尼速度差异较小,甚至可以忽略,但考虑到试验过程及数据判读时存在一定的误差,还不能直接确定两种铜管比动能的偏差范围,因此,这里通过分析铜管变形能与动能的比值以估算两种铜管试验数据的偏差上限。2.2 铜管变形能对炸药驱动性能表征参量的影响
常用于表征炸药驱动性能的铜管比动能、格尼速度等参量仅考虑爆轰转换的动能,未包括铜管的应变能、绝热压缩后的热能损耗等。铜材的晶粒尺寸增大后,铜管的变形能可能会随之减小[11],但若铜管的变形能与动能的比值较小,则两种铜管的试验数据不会有显著差异,更不会影响炸药驱动性能的比对与评价。
对于铜管应变能的计算,需要确定铜管流动应力
Y 的变化,这里采用Steinberg-Guinan模型[15-16]进行描述:Y=Y0(1+36ε)0.45G/G0≤Ymax (1) 式中:
Y0 为准静态屈服应力,Ymax 为最大流动应力,ε=ln(rm/rm0) 为铜管的真实应变,rm 为铜管质量中心面的半径(rm0 为其初始值);G 为剪切模量,G0 为其初始值,G/G0 与铜管的应变率相关。然而,由于铜管的真实应变率˙ε 接近105 s−1,且初期压力较高,导致铜材的温度快速上升,使得膨胀中后期的应变率效应会有所下降[1],因此,G/G0 与˙ε 的函数关系将较为复杂。为了便于计算,可根据炸药的猛度将G/G0 近似表示为一个常数值,即[15]G/G0≈0.368+0.151D (2) 式中:D为炸药的爆速,单位为 mm/μs。该试验所装填的JO-159药柱的爆速为8.76 mm/μs,则
G/G0 ≈1.69。由于目前在公开报道的文献中没有发现针对不同晶粒度退火无氧铜在高应变率拉伸过程的Y0 和Ymax 的值,因此,这里暂不考虑晶粒度的影响,仅采用文献[16]中的相应参数值(即Y0=0.049 GPa,Ymax=0.26 GPa)分析铜管应变能的变化过程。此外,铜管的应变范围较大,结合公式(1)可以计算出,当应变达到临界值ε∗ ≈0.325(即Δre ≈5 mm)时,流动应力Y 达到最大值,那么单位体积铜管的应变能计算式可写为Es=∫ε0Ydε={Y036GG0(1+36ε)1.451.45ε≤ε∗Y036GG0(1+36ε∗)1.451.45+Ymax(ε−ε∗)ε>ε∗ (3) 而单位质量铜管的变形能与动能之间的比值可表示为
η=EsEρCu (4) 式中:
ρCu 为铜管的密度。基于细晶铜管膨胀过程的扫描试验数据,并结合式(3)和(4)可分别计算出
Es-Δre 和η-Δre 曲线,如图6所示。从图中可以看出,随着铜管膨胀位移的增大,Es 和η 均呈递增趋势;在V=2.2,4.4,7.0三个特征位置,η 处于0.7%~1.5%范围,这说明铜管的变形能与动能的比值处于较小的范围,因此,可认为图5中两种晶粒度铜管的比动能偏差范围是基本合理的。当然,若炸药的作功能力有所降低,则η 值可能有所增大,以密度为1.63 g/cm3的TNT为例,V=7.0时的比动能约为JO-159的60%[12],则可以估算出,此时η ≈2.5%。考虑到两种晶粒度铜材在同种工况下的应变能相差20%~30%[11],则对于大多数高能炸药,两种铜管的比动能偏差可能不会超过1%。3. 铜管断裂特性及分析
3.1 两种晶粒度铜管断裂特性的差异
扫描试验中,虽然可以根据边界线是否存在“突跃点”来判断铜管的断裂位置,但仅适用于膨胀早期或中期,无法判断铜管在膨胀后期的具体断裂位置,更不可能判断出不同晶粒度铜管的具体断裂方式和断裂过程,难以为铜管材质选择和处理提供借鉴。因此,这里采用分幅观测试验分析两种晶粒度铜管的断裂过程。
该试验采用SJZ-15型转镜式高速分幅相机观测铜管的表面,并将铜管距离起爆端200 mm位置置于相机视场的中心,以便与狭缝扫描试验的结果进行对比。分幅观测试验的布局如图7所示,为了提升成像质量,试验时采用两发氩气弹从不同的方向对铜管进行同步照明,但布局时应控制氩气弹与铜管之间的距离,避免氩气弹的强光直接进入相机视场。试验时,分幅相机的摄影频率设置为5×105 s−1。
铜管膨胀时间一般为数十微秒,在膨胀的早、中期一般不会出现明显的裂纹,因此,重点关注膨胀的中后期。为了更明显地看出铜管在不同时刻的变化,这里取裂纹扩展过程几个典型时刻的图片,如图8所示。结合图中红色虚直线可以看出,在JO-159爆轰加载下,粗晶和细晶铜管被观测截面从开始膨胀至断裂分别约经历21 μs和28 μs。结合图4可以看出,21 μs时刻对应的膨胀位移
Δre 约为35 mm,而28 μs时刻已超过图4的显示范围,由于JO-159炸药为理想炸药,铜管膨胀后期的速度基本不再增加,即铜管膨胀后期的Δre-t 曲线近似为线性,因此将其进行相应地延伸则可确定出28 μs时刻对应的Δre 约为47 mm。从而可进一步计算出粗晶和细晶铜管的真实断裂应变分别约为0.77和0.93,所对应的铜管壁厚分别为0.69和0.56 mm。这表明细晶铜管的延展性较粗晶铜管有了显著提升,但粗晶铜管的断裂直径也达到了初始直径的3倍,满足圆筒试验的基本要求。对于两种铜管的断裂过程,虽然其裂纹的形成主要是由于铜管壁厚减少至一定程度后,变形过程不稳定,引起局部应力集中和温升,进而形成颈缩,但从图8中可以看出,两种铜管的裂纹扩展和分布有明显差异。图8(a)显示,当粗晶铜管接近断裂应变时,其表面出现较多随机分布的孔洞(如图中红色椭圆圈内所示),由于孔洞周围易出现应力集中,且孔洞之间的距离较近,使得这些孔洞没有出现明显的孤立增长趋势,而是在孔洞间快速形成大量交错的连接裂纹,最终产生形状各异的密集小碎片。而晶粒细化后,当铜管膨胀了24 μs时,其表面已出现了少量的孔洞(如图8(b)的红色虚直线上方所示);随着铜管的继续膨胀,该孔洞呈孤立增长趋势(如图中红色圆圈内所示),直至随后而至的裂纹将其贯穿;同时,大量裂纹主要沿铜管母线快速扩展,但不同区域的裂纹扩展速度差异较大(图中蓝色椭圆圈标示出了扩展速度较快的裂纹),最终形成相对较宽的长条状破片。考虑到铜管出现明显的孔洞后,会造成产物的溢出,使得后续的实验数据意义不大,因此,虽然细晶铜管的整体延展性明显提升,但孔洞出现前的有效膨胀时间不足24 μs,此时
Δre ≈40 mm,略大于粗晶铜管的最大有效膨胀距离。3.2 粗晶铜管装填不同炸药时断裂特性的差异
虽然粗晶铜管可满足作功能力较强的JO-159炸药的圆筒试验要求,但对于其他炸药是否适用,且断裂过程是否类似,还需要进一步试验验证,为此,这里对于装填TNT药柱(密度为1.58 g/cm3)的工况进行了分幅观测试验,以对比分析其断裂过程的差异。
从图9中可以看出,在TNT爆轰加载下,粗晶铜管被观测截面从开始膨胀至发生断裂约经历26 μs,结合扫描试验的数据可知[12],此时
Δre ≈32 mm,略小于装填JO-159炸药的工况,仍满足圆筒试验的要求。此外,对比图9和图8(a)可以发现,两种工况中,裂纹形成之前均首先产生了随机分布的孔洞,由于裂纹紧随其后,使得这些孔洞均未出现明显的孤立增长趋势。所不同的是,TNT加载时的孔洞数量明显减少,且孔洞间的距离相对增大,使得孔洞间交错裂纹的数量大幅度降低;此外,大部分裂纹主要沿铜管母线形成及扩展,并最终形成条状破片,但破片的宽度及长度均小于图8(b)所示工况。对于TNT加载时孔洞数量减少的原因,可尝试从粗晶粒无氧铜成核时断裂强度的变化规律进行分析,文献[17]认为铜的断裂强度σp 、真实应变率˙ε 、温度Tp 存在如下关系:σp=σ0˙εmf(Tp) (5) 式中:
σ0 为准静态断裂强度;˙ε =103~107 s−1时,m=0.2;f(Tp) 服从反幂指数规律。对于铜管断裂时的真实应变率差异,考虑到铜管膨胀后期的速度差异较小,可以
Δre =30 mm时为例进行计算,此时TNT和JO-159加载下,铜管质量中心面的速度um 分别约为1.40、1.84 mm/μs,则˙ε 分别为4.7×104、6.1×104 s−1。即TNT加载下,裂纹临近形成时,铜管的真实应变率˙ε 较JO-159低23%,但˙εm 的差异仅为5%左右,因此,真实应变率˙ε 的差异可能不是两种炸药加载下铜管断裂过程和破片形状发生显著变化的主要原因。铜管断裂前的温度主要由两部分组成,即加载初期,铜管所受的高压冲击后产生的温升,以及铜管膨胀过程中,应变能转换所产生的温升。对于冲击所致的温升,需首先计算铜管所受的初始载荷,这可结合爆轰产物及铜管的p(u)曲线进行计算。对于滑移爆轰加载方式,其爆轰产物的p(u)为:[1]
u=2Dγγ2−1[1−(ppCJ)γ−12γ] (6) 式中:TNT和JO-159的多方指数γ分别为3.11和2.97[12],爆速D分别为6.87、8.76 mm/μs;爆压
pCJ=ρ0D2/(γ+1) ,ρ0 为炸药密度。铜管的p(u)为:p=ρCuu(c0+λu) (7) 式中:
ρCu =8.93g/cm3,c0 =3.94 mm/μs,λ =1.489[1]。由公式(6)和(7)可计算出TNT和JO-159爆轰加载初始时刻,铜管所受的压力分别为14.5、26.5 GPa,根据文献[18]中铜的相关参数表可计算出冲击压缩导致的温升分别约为72、150 K,而卸载后的温升分别降为16、52 K。然后,假定铜管膨胀至Δre =30 mm的应变能全部转换为热量,结合式(3)则可计算出TNT和JO-159加载下铜管最终的温升分别为94、134 K。虽然两种炸药加载工况下铜管的温升差异约为40 K,但温度的绝对值较低,远未达到其退火温度,因此可能不足以造成断裂强度发生显著变化。综上分析认为,真实应变率和温升可能均不足以造成两种工况的断裂过程出现明显差异,那么最大的关注点则成为两种工况下铜管所受压力的差异。爆轰产物对铜管的作用力差异较大时,可能导致材料微观结构变化,进而影响其断裂行为,为此,采用扫描电镜对两种炸药爆轰加载下的铜管断口形貌进行了观测,如图10所示。
图10(a)显示,TNT加载下,铜管断口为韧窝和解理条纹的混合型断口,表现出一定的韧性断裂特征;而图10(b)显示,JO-159加载后,铜管的断口主要由解理条纹组成,未见韧窝出现,表明其脆性断裂特征明显。可见,铜管所受初始压力由14.5 GPa提升至26.5 GPa后,诱导其断裂模式发生了转变,进而导致其宏观断裂形态出现较大差异,但这方面的转变机理还有待于进一步研究。
从应用方面看,大多数炸药的猛度及作功能力处于TNT及JO-159的范围内,则粗晶铜管的断裂状况应介于图9和图8(a)之间,其最大有效膨胀位移也能大于30 mm,即断裂直径达到初始直径的3倍,基本满足圆筒试验要求。
4. 结 论
(1)由于标准铜管的变形能与动能的比值较小,晶粒尺寸为100~300 μm的粗晶无氧铜管与晶粒尺寸为20~30 μm的细晶无氧铜管的比动能相对偏差小于1%。
(2)JO-159爆轰加载下,细晶铜管的整体延展性能较粗晶铜管明显提升,其裂纹要沿母线扩展,形成条状破片;但裂纹形成之前,内部少量缺陷成核后,会形成孤立增长的孔洞,造成爆轰产物提前泄露,使得铜管的有效膨胀位移约为40 mm。
(3)粗晶铜管在JO-159和TNT加载下,其断裂直径均达到了初始直径的3倍,满足圆筒试验的基本要求。两种工况下的断裂过程类似,均首先出现随机分布的孔洞,但JO-159加载时孔洞数量较多,且距离较近,使得孔洞间快速形成大量交错的裂纹,最终产生密集小碎片;而TNT加载时,可能由于加载压力大幅度降低,使得孔洞数量大幅减少,且距离增大,裂纹主要沿母线扩展,形成条状破片。
-
表 1 衰减系数识别数值结果
Table 1. Numerical results of attenuation coefficient identification
步数 δ S(δ) pk P 0 2 000 000 1 171.61 0.004 16 1 2 000 010 1 171.53 10 0.004 16 2 2 000 030 1 171.36 20 0.004 16 3 2 000 070 1 171.03 40 0.004 16 4 2 000 150 1 170.36 80 0.004 16 5 2 000 310 1 169.03 160 0.004 16 6 2 000 630 1 166.37 320 0.004 15 7 2 001 270 1 161.06 640 0.004 14 8 2 002 550 1 150.47 1 280 0.004 12 9 2 005 110 1 129.46 2 560 0.004 08 10 2 010 230 1 088.07 5 120 0.004 00 11 2 020 470 1 007.77 10 240 0.003 84 12 2 040 950 857.06 20 480 0.003 52 13 2 081 910 594.38 40 960 0.002 90 14 2 163 830 217.92 81 920 0.001 71 15 2 291 189 0.074 463 127 259 3.05×10-5 16 2 293 600 9.72×10-9 2 442 1.10×10-8 表 2 模拟波形与实测波形拟合结果表
Table 2. The omparaison of simulation currents and measured currents
t/ns Isim/A Iact/A E/% 10 285.4 280 -1.93 20 551.1 540 -2.06 30 794.2 820 3.15 40 1 012.0 1 060 4.53 50 1 203.0 1 250 3.76 60 1 366.0 1 420 3.80 70 1 499.0 1 540 2.66 80 1 603.0 1 620 1.05 90 1 677.0 1 700 1.35 100 1 722.0 1 680 -2.50 110 1 738.0 1 700 -2.24 120 1 728.0 1 680 -2.86 130 1 692.0 1 640 -3.17 140 1 632.0 1 540 -5.97 150 1 551.0 1 460 -6.23 160 1 452.0 1 360 -6.76 -
[1] 陈景亮, 姚学玲, 孙伟.脉冲电流技术[M].西安: 西安交通大学出版社, 2008. [2] 戴建华, 李开成.基于Rogowski线圈的大电流测量[J].高电压技术, 2002, 28(1): 6-10.Dai Jian-hua, Li Kai-cheng. Heavy current measurement based on Rogowski coil[J]. High Voltage Engineering, 2002, 28(1): 6-10. [3] Eriksson J, Wedin P A, Gulliksson M E. Regularization methods for uniformly rank-deficient nonlinear least-squares problems[J]. Journal of Optimization Theory and Applications, 2005, 127(1): 1-26. doi: 10.1007%2Fs10957-005-6389-0 [4] Zhang J Z, Xue Y, Zhang K. A structured secant method based on a new quasi-newton equation for nonlinear least squares problems[J]. Bit Numerical Mathematics, 2003, 43(1): 217-229. doi: 10.1023/A:1023665409152 [5] 邹志伟.求解对称非线性方程组的一种修正共轭梯度法[D].长沙: 湖南大学, 2010. [6] 孙风建.基于新拟牛顿方程的非线性最小二乘的一类新算法[D].南京: 南京理工大学, 2007. [7] Fan Jin-yan. The modified Levenberg-Marquardt method for nonlinear equations with cubic convergence[J]. Mathematics of Computation, 2011, 81(277): 447-466. [8] 王新洲.非线性模型参数估计理论与应用[M].武汉: 武汉大学出版社, 2002. 期刊类型引用(1)
1. 张晓兵, 赵骥宇, 袁士伟. 水下冲击大电流间接测量方法研究. 科学技术与工程. 2016(24): 210-213 . 百度学术
其他类型引用(1)
-