高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于脉冲激光的发射药参数快速检测方法

吴清潮 吴坚

引用本文:
Citation:

基于脉冲激光的发射药参数快速检测方法

    通讯作者: 吴坚, jxjawj@mail.xjtu.edu.cn
  • 中图分类号: TN24; O348.11

Rapid detection method of propellant parameters based on pulsed laser

    Corresponding author: WU Jian, jxjawj@mail.xjtu.edu.cn
  • CLC number: TN24; O348.11

  • 摘要: 为了验证基于脉冲激光作用的发射药能量性能参数快速测量方法的有效性,采用实验方法研究了脉冲激光作用下单、双、三基发射药和单晶冰糖的冲击波及等离子体特性,得到了激光参数对脉冲激光作用过程的影响规律以及单基发射药与单晶冰糖的冲击波、等离子体图像差异。结果表明,单基发射药的冲击波及等离子体膨胀特性主要受其化学反应释能的影响,而冰糖主要受激光辐照度的影响;采用支持向量机回归算法,建立了单、双、三基发射药激光作用下的冲击波特征速度与其火药力、爆热、爆温的线性定标模型,得到的决定系数R2值分别为0.9912、0.9998、0.9999。该方法对发射药火药力、爆热及爆温有较好的预测能力,为脉冲激光与含能材料作用的研究提供了参考。
  • 图 1  实验平台示意图

    Figure 1.  Schematic diagram of experimental platform

    图 2  发射药及单晶冰糖冲击波阴影图像

    Figure 2.  Shockwave shadow image of propellants and monocrystal rock sugar

    图 3  发射药及单晶冰糖等离子体演化图像

    Figure 3.  Plasma evolution image of propellants and monocrystal rock sugar

    图 4  不同激光参数条件下的单晶冰糖及单基发射药冲击波阴影图像

    Figure 4.  Shock wave shadow image of monocrystal rock sugar and single base propellant under different laser parameters

    图 5  不同激光参数条件下的单晶冰糖及单基发射药等离子体演化图像

    Figure 5.  Plasma evolution image of monocrystal rock sugar and single base propellant under different laser paremeters

    图 6  发射药冲击波速度及拟合曲线

    Figure 6.  Shockwave velocity and fitting curve of propellants

    图 7  冲击波传播距离修正流程

    Figure 7.  Correction process of shockwave propagation distance

    图 8  发射药冲击波速度及拟合曲线(修正后)

    Figure 8.  Shockwave velocity and fitting curve of propellants (after correction)

    图 9  火药力线性定标结果

    Figure 9.  Linear calibration result of explosive force

    图 10  爆热线性定标结果

    Figure 10.  Linear calibration result of explosive heat

    图 11  爆温线性定标结果

    Figure 11.  Linear calibration result of explosive temperature

    表 1  样品配方

    Table 1.  Sample formula

    sample formula(mass fraction)
    single base propellant nitrocotton(NC) 0.97(nitrogen content 0.126),additive 0.03
    double base propellant nitrocotton(NC) 0.55(nitrogen content 0.126),nitroglycerin(NG) 0.43,additive 0.02
    triple base propellant nitrocotton(NC) 0.25(nitrogen content 0.126),nitroglycerin(NG) 0.30,Hexogen(RDX) 0.45
    monocrystal rock sugar sucrose 1.00
    下载: 导出CSV

    表 2  发射药冲击波特征速度

    Table 2.  Shockwave characteristic velocity of propellants

    propellant shockwave characteristic velocity/(m·s-1) shockwave characteristic velocity(after correction)/(m·s-1)
    single base 1296 1047
    double base 1256 982.5
    triple base 1492 1140
    下载: 导出CSV

    表 3  发射药性能参数

    Table 3.  Propellant property parameters

    propellant explosive force/(kJ·kg-1) explosive heat/(kJ·kg-1) explosive temperature/K
    single base 1028 3876 2997
    double base 931.0 3188 2580
    triple base 1226 4812 3618
    下载: 导出CSV

    表 4  发射药性能参数预测

    Table 4.  Prediction of propellants property parameters

    propellant explosive force/(kJ·kg-1) explosive heat/(kJ·kg-1) explosive temperature/K
    measurement prediction error measurement prediction error calculation prediction error
    single base 1028 1000 2.7% 3876 3636 6.2% 2997 2857 4.7%
    double base 931.0 954.0 2.5% 3188 3388 6.3% 2580 2697 4.5%
    triple base 1226 1231 0.4% 4812 4853 0.9% 3618 3642 0.7%
    下载: 导出CSV

    表 5  发射药性能参数预测(修正后)

    Table 5.  Prediction of propellants property parameters (after correction)

    propellant explosive force/(kJ·kg-1) explosive heat/(kJ·kg-1) explosive temperature/K
    measurement prediction error measurement prediction error calculation prediction error
    single base 1028 1044 1.6% 3876 3863 0.3% 2997 3004 0.2%
    double base 931.0 921.0 1.1% 3188 3196 0.3% 2580 2576 0.2%
    triple base 1226 1219 0.6% 4812 4818 0.1% 3618 3615 0.1%
    下载: 导出CSV
  • [1] 李强, 闫光虎, 张玉成, 等. NG含量对改性单基发射药燃烧性能的影响[J]. 火炸药学报, 2012, 35(1): 73-76.

    LI Q, YAN G H, ZHANG Y Ch, et al. Effect of NG content on burning performance of modified sing-base gun propellant[J]. Chinese Journal of Explosives & Propellants, 2012, 35(1): 73-76(in Chin-ese). 
    [2]

    THOMAS J C, MORROW G R, DILLIER C A M, et al. Comprehensive study of ammonium perchlorate particle size/concentration effects on propellant combustion[J]. Journal of Propulsion and Power, 2019, 36(1): 1-6.
    [3] 郝海霞, 裴庆, 赵凤起, 等. 固体推进剂激光点火性能研究综述[J]. 含能材料, 2009, 17(4): 491-498.

    HAO H X, PEI Q, ZHAO F Q, et al. Summarization of laser ignition characteristics of solid propellant[J]. Chinese Journal of Energetic Materials, 2009, 17(4): 491-498 (in Chinese). 
    [4] 张陆, 王霆威, 王晓军, 等. 激光敏感型含能配合物类起爆药研究进展[J]. 含能材料, 2022, 30(4): 385-395.

    ZHANG L, WANG T W, WANG X J, et al. Review on laser sensitive energetic complex primary explosives[J]. Chinese Journal of Energetic Materials, 2022, 30(4): 385-395 (in Chinese). 
    [5] 刘彦汝, 孙杰, 金波, 等. 360 nm紫外激光辐照下HMX晶体的微观结构变化[J]. 含能材料, 2021, 29(12): 1208-1215.

    LIU Y R, SUN J, JIN B, et al. Microstructure changes of HMX crystals irradiated by 360 nm UV laser[J]. Chinese Journal of Energetic Materials, 2021, 29(12): 1208-1215 (in Chinese). 
    [6] 伍俊英, 刘嘉锡, 杨利军, 等. 不同频率飞秒激光脉冲序列加工炸药过程安全性的数值计算[J]. 含能材料, 2021, 29(3): 192-201.

    WU J Y, LIU J X, YANG L J, et al. Numerical calculation of the safety of processing explosives with femtosecond laser sequence with different frequencies[J]. Chinese Journal of Energetic Materials, 2021, 29(3): 192-201.
    [7]

    GOTTFRIED J L. Laser-induced plasma chemistry of the explosive RDX with various metallic nanoparticles[J]. Applied Optics, 2012, 51(7): B13-B21. doi: 10.1364/AO.51.000B13
    [8]

    KIMBLIN C, TRAINHAM R, CAPELLE G A, et al. Characterization of laser-induced plasmas as a complement to high-explosive large-scale detonations[J]. AIP Advances, 2017, 7(9): 095208. doi: 10.1063/1.4999793
    [9]

    HAUER M, FUNK D J, LIPPERT T, et al. Time resolved study of the laser ablation induced shockwave[J]. Thin Solid Films, 2004, 453/454(4): 584-588.
    [10]

    LUSSELL F C D, HARMON R S, MCNESBY K L, et al. Laser-induced breakdown spectroscopy analysis of energetic materials[J]. Applied Optics, 2003, 42(30): 6148-6152. doi: 10.1364/AO.42.006148
    [11]

    ZHANG Zh, WANG A, WU J, et al. Spatial confinement effects of bubbles produced by laser ablation in liquids[J]. AIP Advances, 2019, 9(12): 125048. doi: 10.1063/1.5127261
    [12] 葛一凡, 陆旭, 刘玉柱. 基于激光诱导击穿光谱和神经网络的蛋壳研究[J]. 激光技术, 2022, 46(4): 532-537.

    GE Y F, LU X, LIU Y Zh. Research on eggshell via laser-induced breakdown spectroscopy and neural network[J]. Laser Technology, 2022, 46(4): 532-537(in Chinese). 
    [13]

    GOTTFRIED J L. Influence of exothermic chemical reactions on laser-induced shock waves[J]. Physical Chemistry Chemical Physics, 2014, 16(39): 21452-21466. doi: 10.1039/C4CP02903H
    [14]

    GOTTFRIED J L. Laboratory-scale method for estimating explosive performance from laser-induced shock waves[J]. Propellants, Explosives, Pyrotechnics, 2015, 40(5): 674-681. doi: 10.1002/prep.201400302
    [15]

    COLLINS E S, GOTTFRIED J L. Laser-induced deflagration for the characterization of energetic materials[J]. Propellants, Explosives, Pyrotechnics, 2017, 42(6): 592-602. doi: 10.1002/prep.201700040
    [16]

    KALAM S A, MURTHY N L, MATHI P, et al. Correlation of molecular, atomic emissions with detonation parameters in femtosecond and nanosecond LIBS plasma of high energy materials[J]. Journal of Analytical Atomic Spectrometry, 2017, 32(8): 1535-1546. doi: 10.1039/C7JA00136C
    [17]

    BISS M M, BROWN K E, TILGER C F. Ultra-high fidelity laser-induced air shock from energetic materials[J]. Propellants, Explosives, Pyrotechnics, 2020, 45(3): 396-405. doi: 10.1002/prep.201900130
    [18] 王茜蒨, 赵宇, 卢小刚, 等. 激光诱导击穿光谱与拉曼光谱技术在危险物检测中的研究进展[J]. 光谱学与光谱分析, 2017, 37(8): 2430-2434.

    WANG X Q, ZHAO Y, LU X G, et al. Progress in laser induced breakdown spectroscopy and Raman spectroscopy for hazardous material detection[J]. Spectroscopy and Spectral Analysis, 2017, 37(8): 2430-2434 (in Chinese). 
    [19] 郭文灿, 郑贤旭, 张旭, 等. 含铝炸药在激光烧蚀下的发射光谱分布及瞬态温度测量[J]. 含能材料, 2018, 26(8): 671-676.

    GUO W C, ZHENG X X, ZHANG X, et al. Emission spectrum distribution and transient temperature measurement of aluminized explosives under laser ablation[J]. Chinese Journal of Energetic Materials, 2018, 26(8): 671-676 (in Chinese). 
    [20]

    VAPNIK V. The nature of statistical learning theory[M]. Berlin, Germany: Springer, 2000: 15-174.
  • [1] 黄宇姜银方金华黄利伟张仕培 . 环形激光诱导冲击波的传播及其对层裂的影响. 激光技术, 2013, 37(3): 301-305. doi: 10.7510/jgjs.issn.1001-3806.2013.03.007
    [2] 王飞姚红兵袁蓓于水生张永康 . 激光冲击强化的冲击波幅值和飞行时间的测量. 激光技术, 2010, 34(5): 694-696. doi: 10.3969/j.issn.1001-3806.2010.O5.032
    [3] 戴蜀娟刘富荣杨晓邹鸿承安承武陆冬生 . 激光冲击在LY12CZ中形成冲击波的研究. 激光技术, 1997, 21(6): 330-333.
    [4] 姚兵王甦曹雁刘迅曾传相 . 激光脉冲和放电脉冲产生冲击波的研究. 激光技术, 2000, 24(1): 12-14.
    [5] 陈建平倪晓武陆建卞保民王亚伟 . 用于探测激光等离子体冲击波的光纤传感器. 激光技术, 2001, 25(2): 85-90.
    [6] 张建花银群曹将栋 . 激光冲击铜薄膜的应力波传播特性模拟分析. 激光技术, 2016, 40(4): 601-605. doi: 10.7510/jgjs.issn.1001-3806.2016.04.030
    [7] 陈亮游利兵罗贤锋尹广玥方晓东 . 应用LIBS技术检测食盐中的镉. 激光技术, 2019, 43(1): 6-10. doi: 10.7510/jgjs.issn.1001-3806.2019.01.002
    [8] 倪晓武陈笑陆建 . 激光与液态物质相互作用机理的研究进展. 激光技术, 2002, 26(4): 258-261.
    [9] 余小敏陈慧敏李家泽 . 激光乳化白内障的新进展. 激光技术, 2003, 27(5): 400-402,406.
    [10] 倪晓武陆建贺安之王平秋马孜周九林 . 高功率激光对光学介质薄膜破坏机理的研究进展. 激光技术, 1994, 18(6): 348-352.
    [11] 田润妮王俊波邱荣周强蒋勇杨永佳 . 纳秒和飞秒激光烧蚀单晶硅的超快诊断. 激光技术, 2015, 39(6): 765-768. doi: 10.7510/jgjs.issn.1001-3806.2015.06.007
    [12] 陈子琪王新兵左都罗 . CO2激光诱导液滴射流等离子体的实验研究. 激光技术, 2016, 40(6): 888-891. doi: 10.7510/jgjs.issn.1001-3806.2016.06.023
    [13] 刘丽王声波吴鸿兴郭大浩廖培育 . 强激光诱导冲击波的实验研究. 激光技术, 2007, 31(2): 134-136.
    [14] 黄利伟龙昆黄宇郭海平丰琅周俊马睿佳 . 激光冲击波作用下3003H16铝合金板料层裂的研究. 激光技术, 2013, 37(2): 270-273. doi: 10.7510/jgjs.issn.1001-3806.2013.02.032
    [15] 王伟张永康鲁金忠殷苏民 . 基于激光冲击波三维无损打标的数值模拟. 激光技术, 2008, 32(1): 37-39,43.
    [16] 花银群吉光孙真真陈瑞芳 . 激光冲击波对ZnO压敏陶瓷电学性能的影响. 激光技术, 2011, 35(1): 133-136. doi: 10.3969/j.issn.1001-3806.2011.01.036
    [17] 叶亚云袁晓东向霞王海军晏良宏陈猛吕海兵贺少勃 . 激光冲击波清洗K9玻璃表面SiO2颗粒的研究. 激光技术, 2011, 35(2): 245-248. doi: 10.3969/j.issn.1001-3806.2011.02.028
    [18] 罗勇张华李月华肖敏徐健宁 . TIG焊快速制造激光视觉检测系统的研究. 激光技术, 2007, 31(4): 367-369.
    [19] 赵新才吴云峰方亮刘宁文 . 激光干涉测速中的小波快速算法研究. 激光技术, 2012, 36(5): 649-652. doi: 10.3969/j.issn.1001-3806.2012.05.019
    [20] 宋华峣王辉林曹泉泉张守宇秦正健 . 传递函数快速相关法的激光测速技术研究. 激光技术, 2020, 44(1): 61-65. doi: 10.7510/jgjs.issn.1001-3806.2020.01.011
  • 加载中
图(11) / 表(5)
计量
  • 文章访问数:  541
  • HTML全文浏览量:  343
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-02-23
  • 录用日期:  2023-04-17
  • 刊出日期:  2024-03-25

基于脉冲激光的发射药参数快速检测方法

    通讯作者: 吴坚, jxjawj@mail.xjtu.edu.cn
  • 1. 西安交通大学 电力设备电气绝缘国家重点实验室,西安 710049,中国
  • 2. 国网厦门供电公司,厦门 361001,中国

摘要: 为了验证基于脉冲激光作用的发射药能量性能参数快速测量方法的有效性,采用实验方法研究了脉冲激光作用下单、双、三基发射药和单晶冰糖的冲击波及等离子体特性,得到了激光参数对脉冲激光作用过程的影响规律以及单基发射药与单晶冰糖的冲击波、等离子体图像差异。结果表明,单基发射药的冲击波及等离子体膨胀特性主要受其化学反应释能的影响,而冰糖主要受激光辐照度的影响;采用支持向量机回归算法,建立了单、双、三基发射药激光作用下的冲击波特征速度与其火药力、爆热、爆温的线性定标模型,得到的决定系数R2值分别为0.9912、0.9998、0.9999。该方法对发射药火药力、爆热及爆温有较好的预测能力,为脉冲激光与含能材料作用的研究提供了参考。

English Abstract

    • 含能材料的反应性能研究与参数测试方法通常基于外场试验、爆炸容器实验开展[1-2],存在实验条件要求高、精细测量实施难等问题。同时,含能材料新配方、新材料在研发初期由于制备量有限而难以开展宏观反应性能测试。因此,迫切需要发展含能材料反应行为研究和评价的新方法。

      激光具有能量纯净、安全可控等优势,其与含能材料作用的研究得到广泛关注,主要包括激光点火[3]、激光起爆[4]、含能化合物激光辐照效应[5]、飞秒激光加工炸药[6]及纳秒脉冲激光作用含能材料[7]等研究。高功率纳秒脉冲激光作用于含能材料表面,引起材料瞬态相变和局部微爆,由于能量小、耗散快一般不会实现持续的燃烧爆炸,但其作用于含能材料形成的温度和压力环境,与含能材料实际应用时燃烧爆炸过程接近[8],进一步与先进的诊断测试方法结合[9-12],可成为精细研究含能材料点火等反应过程的新手段。

      美国陆军实验室的GOTTFRIED[7]基于纳秒脉冲激光与黑索金(cyclotrimethylenetrinitramine, RDX)的作用,证实了激光作用后烧蚀产生的含能材料会发生化学反应,含能材料普遍展现出比非爆炸材料更强的冲击波扩散和产物膨胀特征[13]。在上述研究的基础上,GOTTFRIED又进一步研究了脉冲激光与二硝基苯甲醚(dinitroanisole,DNAN)、三硝基甲苯(trinitrotoluene,TNT)、六硝基芪(hexanitrostilbene,HNS)、三氨基三硝基苯(triamino-trinitrobenzene,TATB)、硝基三唑酮(nitrotriazolone,NTO)、四硝酸季戊四醇酯(pentaerythritol tetranitrate,PETN)、RDX、奥克托今、六硝基六氮杂异伍兹烷(hexanitrohexaazaisowurtzitane,CL-20)等含能材料的作用特性,建立了激光作用后空气中冲击波特征速度与爆速的实验测量值,以及爆压、爆热等的理论计算值的定标曲线(laser-induced air shock from energetic materials,LASEM)方法,进而基于定标曲线测量了二氨基二硝基乙烯(diaminodinitroethy-lene, FOX-7)、纳米RDX、硝化纤维素等含能材料的爆速、爆热等参数,并估计了新型高氮炸药、新配方炸药的爆轰性能[14]。COLLINS等人[15]基于纳秒脉冲激光及高速彩色摄像机等诊断设备,对8种炸药的激光爆燃过程进行了温度、能量释放和发射特征的测量和分析,认为该实验室规模的研究方法适用于新型含能材料的开发。随后,KALAM等人[16]及BISS等人[17]也基于此研究方法进行了含能材料性能方面的研究。国内,WANG等人[18]研究分析了激光光诱导击穿光谱技术及喇曼光谱技术在爆炸物及危险物检测中的应用,提出将二者结合可以有效地提升检测准确度。GUO等人[19]采用脉冲激光烧蚀的手段激发了RDX基含铝炸药粉末的快速反应,结果表明,添加铝粉后炸药的光谱分布和反应行为存在明显区别,与炸药宏观性能有一定的对应性。

      由此可见,脉冲激光与含能材料的作用行为研究已成为一种深入揭示含能材料反应行为和机理的手段,特别是国外学者已基于激光作用下的冲击波速度、等离子体谱线强度比值等参数,建立了其与含能化合物宏观反应性能的定标曲线。但相关研究主要针对单一组分的含能化合物,对发射药等的研究仍非常有限,国内对该方法的研究更是非常缺乏。

      针对上述问题,本文作者研究了脉冲激光与单基发射药、双基发射药、三基发射药及单晶冰糖(非含能材料)的作用行为,分析了激光能量、光斑大小对3种发射药及冰糖的冲击波、等离子体特性的影响,并基于激光作用下的冲击波特性,建立了冲击波特征速度与单、双、三基发射药的火药力、爆热及爆温的定标关系,验证了脉冲激光作用参数用于测量发射药能量性能参数的有效性。

    • 本文中所用的样品为单基发射药、双基发射药、三基发射药和单晶冰糖,如表 1所示。单基发射药以石墨作为光泽剂,添加剂为2号中定剂(C15H16N2O1)。双基发射药的添加剂为钝感剂二氧化钛(TiO2)。单晶冰糖的成分为蔗糖,化学式为C12H22O11。本文中所用单基发射药及三基发射药为七孔粒状,双基发射药为管状。为了保证激光烧蚀面的平整,提高样品的统一度。对发射药样品进行了预处理,处理过程为:(a)沿轴线方向,并且避开孔结构对发射药进行纵切;(b)将纵切面在砂纸上轻轻打磨,使表面平整;(c)该纵切面作为激光烧蚀的平面;(d)每次实验进行前都用一发脉冲激光对样品表面杂质进行清理。

      表 1  样品配方

      Table 1.  Sample formula

      sample formula(mass fraction)
      single base propellant nitrocotton(NC) 0.97(nitrogen content 0.126),additive 0.03
      double base propellant nitrocotton(NC) 0.55(nitrogen content 0.126),nitroglycerin(NG) 0.43,additive 0.02
      triple base propellant nitrocotton(NC) 0.25(nitrogen content 0.126),nitroglycerin(NG) 0.30,Hexogen(RDX) 0.45
      monocrystal rock sugar sucrose 1.00
    • 基于脉冲激光作用的发射药反应行为及能量性能参数快速检测的实验平台如图 1所示。实验平台主要装置包括:激光器SGR-60,Beamtech(输出波长1064 nm,最大能量约3 J),输出激光与含能材料作用;激光器Dava-100,Beamtech(输出波长1064 nm,最大能量约100 mJ),用于搭建激光阴影诊断系统;增强型电荷耦合器件(intensified charge-coupled device,ICCD)iStar DH334T,Andor(像素1024×1024,最短曝光时间约2 ns);电荷耦合器件(charge-coupled device,CCD)EOS 750D,Canon(互补金属氧化物半导体图像感应器,约2420万pixel);激光能量计,表头NOVA Ⅱ,Ophir,探头PE50BF-DIF-C,Ophir(探测波长范围:0.19 μm~2.2 μm,能量范围:200 μJ~10 J,损伤阈值:4 J/cm2),用于监测激光能量。图中,L为平凸透镜;BS(beam splitter)为分束镜;R为反射镜;BE(beam expander)为激光扩束镜;EM(energy meter)为激光能量计;PD(photon detector)为光电探测器;OSC(oscilloscope)为示波器;AG(attenuator group)为能量衰减器。

      图  1  实验平台示意图

      Figure 1.  Schematic diagram of experimental platform

      脉冲激光(1064 nm)小部分能量被分光镜BS1反射至光电二极管用于时序监测,绝大部分能量透过分光镜BS1经衰减器衰减后再经第二分光镜BS2分光,部分光束(能量比例已知)反射至能量计进行能量检测,其余部分透过分光镜BS2后经透镜L1聚焦烧蚀样品,在样品表面形成等离子体并引起含能材料发生反应。等离子体的发光图像由ICCD相机搭配镜头直接测量。诊断激光束(1064 nm)由诊断激光器产生,经扩束后光斑覆盖待测区域,通过待测区域后经透镜聚焦,从而在装有1064 nm带通滤波片的CCD相机中成像,诊断激光束的传输同样需经过衰减和分光以确保安全和进行时序监测。

    • 基于上述实验平台,本节中研究发射药和单晶冰糖在不同激光能量和光斑大小条件下的冲击波及等离子体演化规律。激光能量和光斑直径条件包括:200 mJ/889 μm、400 mJ/889 μm和400 mJ/1474 μm。

    • 图 2为400 mJ/1474 μm激光参数下单、双、三基发射药及单晶冰糖的冲击波阴影图像,时间范围为1 μs~5 μs。由图 2可知,相同延时下发射药脉冲的冲击波传播距离显著大于非含能的单晶冰糖,这说明了“含能性”在脉冲激光作用过程中有显著影响。进一步分析,相同延时下三基发射药的冲击波传播距离最大,单基发射药次之,双基发射药最小,这说明了同样是含能的材料,在脉冲激光作用过程中仍存在差异性,该差异本质上与其性能参数有关。值得一提的是,冲击波传播距离不仅体现了平均速度上的区别,还包含了冲击波速度发展的趋势,通过参数化可进一步分析其差异性。

      图  2  发射药及单晶冰糖冲击波阴影图像

      Figure 2.  Shockwave shadow image of propellants and monocrystal rock sugar

      图 3为激光参数400 mJ/1474 μm条件下单、双、三基发射药及单晶冰糖的等离子体演化图像,时间范围为400 ns~2000 ns。由图 3可知,在激光参数条件相同的情况下,3种发射药的等离子体演化过程相似,等离子体形态呈扁条状,随着时间推移,等离子体膨胀并衰减。对比非含能的单晶冰糖,发射药等离子体厚度总体更大,膨胀速度更快。进一步分析,虽然材料的“含能性”增强了等离子体发光,但相同参数条件脉冲激光诱导的发射药等离子体并未表现出更强的持续性,结合上述冲击波特性推测,等离子体演化特性的特点可能与激光诱导发射药等离子体与化学反应区域之间的能量传递关系有关。

      图  3  发射药及单晶冰糖等离子体演化图像

      Figure 3.  Plasma evolution image of propellants and monocrystal rock sugar

      为了进一步研究脉冲激光作用发射药特性,明确作用过程以及冲击波、等离子体特性差异的机理,本文中以单基发射药为例,通过与非含能的单晶冰糖进行对比,研究了激光参数条件对脉冲激光作用过程的影响。

    • 图 4为3种激光参数下单晶冰糖(如图 4a图 4c图 4e所示)和单基发射药(如图 4b图 4d图 4f所示)的冲击波阴影图像,时间范围为脉冲激光起始后1 μs~5 μs。

      图  4  不同激光参数条件下的单晶冰糖及单基发射药冲击波阴影图像

      Figure 4.  Shock wave shadow image of monocrystal rock sugar and single base propellant under different laser parameters

      在激光阴影图像中,以同一时刻冲击波传播距离大小比较冲击波强弱。从图 4单基发射药的激光阴影图像可知,在激光参数为200 mJ/889 μm时,单基发射药的冲击波相比单晶冰糖更强,这可能是由于单基发射药作为含能材料,激光激发相变的物质后续存在一系列释能反应,对冲击波有较大增强作用。当激光参数为400 mJ/889 μm时,两者冲击波强度接近,这是由于此时较大激光功率密度诱导较强等离子体,冲击波强度主要受等离子体强度影响。当激光能量保持为400 mJ不变,激光光斑直径由889 μm增大至1474 μm时,单基发射药的冲击波速度呈现增大趋势,这与单晶冰糖的实验结果不同。激光光斑尺寸的增大使激光辐照区域增大,这导致了材料表面受激物质大大增加,而相变的含能材料物质在后续存在一系列显著的释能反应,反应释放的能量对冲击波有增强作用,而相比之下,释能反应对冲击波的增强作用,远大于激光功率密度下降导致等离子体减弱而造成的冲击波减弱作用,因此冲击波速度总体呈现增大趋势。值得一提的是,含能材料与非含能材料在脉冲激光作用过程冲击波特性上体现出来的重要区别,是基于脉冲激光作用参数评价含能材料性能这一方法可行性的重要依据。

      在等离子体演化图像中,通过发光强度和膨胀距离来比较等离子体的强度和衰减情况。图 5为3种激光参数条件下单晶冰糖的等离子体演化的图像(如图 5a图 5c图 5e所示)和单基发射药的等离子体演化图像(如图 5b图 5d图 5f所示),时间范围为400 ns~2000 ns。分析图 5单晶冰糖的等离子体演化图像可知,当激光光斑直径保持为889 μm不变,激光能量由200 mJ增大至400 mJ时,相同探测延时下等离子体膨胀的距离以及等离子体的尺寸都明显增大,并且等离子体形态在2维图像平面上由扁半圆形变为更加均匀的半圆形,这一现象是因为激光功率密度增大导致了激光诱导等离子体的显著增强。当激光能量保持为400 mJ不变,激光光斑直径由889 μm增大至1474 μm时,激光功率密度的减小导致了等离子体的显著减弱,并且在较大光斑的情况下,等离子体的形态在2维图像平面上呈扁条状,因此可以明显看出, 单晶冰糖在大光斑情况下,等离子体的膨胀距离明显小于小光斑的情况。

      图  5  不同激光参数条件下的单晶冰糖及单基发射药等离子体演化图像

      Figure 5.  Plasma evolution image of monocrystal rock sugar and single base propellant under different laser paremeters

      图 5中单基发射药的等离子体演化图像可知,在激光参数条件由200 mJ/889 μm变为400 mJ/889 μm时,等离子体演化特性的变化规律与单晶冰糖相似,而当激光参数条件为400 mJ/1474 μm时,单基发射药的激光诱导等离子体形态虽然也在2维图像平面上呈条状,但其膨胀距离却介于前两种激光参数条件下的膨胀距离之间,这是由于单基发射药作为含能材料,其受激物质存在后续释能化学反应,释放的能量对等离子体有一定的增强作用,但由于时间尺度太小(百纳秒),化学反应的进行程度很有限,反应释放的能量在一定程度上影响等离子体的强弱。因此,当激光功率密度较大时,相同时刻等离子体膨胀距离显然更大,如激光参数条件为400 mJ/889 μm(功率密度约为80 TW/m2)时;当激光功率密度接近,如激光参数条件为200 mJ/889 μm(功率密度约为40 TW/m2)和400 mJ/1474 μm(功率密度约为29 TW/m2)时,更大光斑尺寸引起更多受激物质后续的释能反应,在一定程度上增强了等离子体,因此后者的等离子体膨胀距离略大于前者。

      进一步分析对比两种样品的实验结果可以明显看出,在两种激光参数条件下,单基发射药的激光诱导等离子体膨胀距离均大于单晶冰糖,但等离子体衰减速度也明显大于单晶冰糖,即单基发射药的激光诱导等离子体起初强于单晶冰糖,其等离子体寿命明显更短。结合冲击波特性,分析如下:单基发射药的激光诱导等离子体起初强于单晶冰糖,这是由于激光直接激发的含能物质反应释放的能量增强等离子体,当直接激发的释能反应结束后,气相以及固相含能物质进行后续反应所需吸收的能量大部分将来自于等离子体,因此能量转移方向变为由等离子体转向含能物质反应区域,而当化学反应向流场释放大量能量时,冲击波被显著增强。值得一提的是,含能材料与非含能材料在等离子体演化特性上所体现出来的能量传递区别,也是基于脉冲激光作用参数评价含能材料性能这一方法可行性的重要依据。

    • 上述研究表明,较大激光光斑尺寸会引起更多质量发射药的反应,使得不同发射药在冲击波及等离子体特性上的差异更加显著。本节中在400 mJ/1474 μm激光参数条件下获得脉冲激光作用下3种发射药的冲击波传播距离,并结合支持向量机(support vector machines, SVM)机器学习方法,围绕单基、双基及三基发射药建立了脉冲激光作用下冲击波特征速度与火药力、爆热及爆温的定标关系。

    • 基于上述实验诊断平台,在激光能量400 mJ、光斑直径1474 μm条件下获取单、双、三基发射药脉冲激光作用过程冲击波阴影图像,并同步监测每发激光的能量。基于不同时刻的激光阴影图像,采用高斯滤波、Canny边缘提取等提取出冲击波前沿位置,从而得到包含1 μs~20 μs共20个散点的冲击波传播距离-时间数据。利用最小二乘法对散点数据进行4次多项式拟合后得到冲击波传播距离变化曲线,对该曲线求导即可得到冲击波速度随时间的变化。将冲击波速度-时间拟合曲线在零时刻取点,可得到该发射药样品的冲击波特征速度。LASEM技术是基于该冲击波特征速度建立与含能材料火药力、爆热及爆温等的定标关系[14]。需要指出的是,冲击波特征速度虽然在数学推导意义上为冲击波的起始速度,但其并无实际的物理意义,仅作为描述脉冲激光作用过程快速变化量的参数。

      图 6为3种发射药冲击波速度散点及其拟合曲线。由曲线y1y2y3的截距可得到特征冲击波速度:单基发射药(1296 m/s)、双基发射药(1256 m/s)、三基发射药(1492 m/s)。由图中误差棒分布可知,单基及三基发射药的冲击波速度抖动较小,而双基发射药在起始阶段(0 μs~5 μs)冲击波速度计算误差较大,从而对特征冲击波速度的截取造成影响。实际研究中,每次重复实验的脉冲激光能量存在微小抖动,对冲击波位置曲线及速度曲线的拟合存在显著影响。

      图  6  发射药冲击波速度及拟合曲线

      Figure 6.  Shockwave velocity and fitting curve of propellants

      为了尽量减小激光能量抖动对特征冲击波速度计算的影响,获取更为可靠、准确的定标基础数据,进一步基于每次实验对应激光能量参数的监测数据,通过SVM[20]的回归算法对获取的冲击波传播距离进行修正,修正流程如图 7所示。

      图  7  冲击波传播距离修正流程

      Figure 7.  Correction process of shockwave propagation distance

      流程中,通过固定时间参数t0为1 μs~20 μs来获取修正后的各时刻冲击波位置散点数据。激光能量参数E0取最大值与最小值的中间值,从而在算法上去除能量抖动影响。

      基于修正后的冲击波传播距离数据计算得到特征冲击波,图 8为经激光能量参数修正后的3种发射药冲击波速度散点及其拟合曲线。由曲线y1y2y3的截距可得到修正后的特征冲击波速度:单基发射药(1047 m/s)、双基发射药(982.5 m/s)、三基发射药(1140 m/s)。

      图  8  发射药冲击波速度及拟合曲线(修正后)

      Figure 8.  Shockwave velocity and fitting curve of propellants (after correction)

      综上所述,在400 mJ/1474 μm激光条件下获得的单、双、三基发射药冲击波特征速度如表 2所示。其在一定程度上体现了发射药的能量性能,下面进一步将其作为关键参数,建立其与发射药火药力、爆热及爆温的实验定标关系。

      表 2  发射药冲击波特征速度

      Table 2.  Shockwave characteristic velocity of propellants

      propellant shockwave characteristic velocity/(m·s-1) shockwave characteristic velocity(after correction)/(m·s-1)
      single base 1296 1047
      double base 1256 982.5
      triple base 1492 1140
    • 由单基发射药与单晶冰糖脉冲激光作用过程冲击波及等离子体特性研究可知,发射药脉冲激光作用冲击波速度与火药力、爆热及爆温的关联性在于:火药力衡量的气体做功能力、爆热及爆温衡量的能量释放及转化水平,对化学反应产物生成速率和气体运动速度产生影响,从而与脉冲激光作用过程的冲击波传播速度存在关联。本文中所采用的3种发射药的火药力、爆热及爆温参数如表 3所示,其中火药力、爆热为实测值,爆温为计算值,实测数据基于GJB770B-2005试验方法测得。

      表 3  发射药性能参数

      Table 3.  Propellant property parameters

      propellant explosive force/(kJ·kg-1) explosive heat/(kJ·kg-1) explosive temperature/K
      single base 1028 3876 2997
      double base 931.0 3188 2580
      triple base 1226 4812 3618

      基于表 2表 3数据,分别建立冲击波特征速度及修正的冲击波特征速度与发射药能量性能参数的线性定标模型,定标结果如图 9~图 11所示。未修正时,冲击波特征速度对于火药力、爆热及爆温的线性定标效果较好,决定系数R2分别为0.9718、0.9255、0.9380,但在95%的置信度下,都存在置信区间(confi-dence interval,CI)较大的问题。由图中红色点线可知,修正后的冲击波特征速度对于火药力、爆热及爆温的线性定标效果得到进一步提高,决定系数R2分别为0.9912、0.9998、0.9999,均优于未经激光能量修正的定标模型,并且在95%的置信度下,修正后的定标模型置信区间都有所减小,以爆热和爆温的定标关系最为显著。

      图  9  火药力线性定标结果

      Figure 9.  Linear calibration result of explosive force

      图  10  爆热线性定标结果

      Figure 10.  Linear calibration result of explosive heat

      图  11  爆温线性定标结果

      Figure 11.  Linear calibration result of explosive temperature

      表 4表 5为3种发射药火药力、爆热、爆温参数的预测值及误差。可见,基于冲击波特征速度建立的线性定标模型,对发射药性能参数的预测值误差均小于7%,在误差允许的条件下可实现火药力、爆热及爆温参数的有效预测。并且由表 5可知,经过激光能量参数修正,定标模型对火药力的预测误差减小至不高于2%,对爆热、爆温的预测误差减小至不高于0.5%,模型预测准确度显著提高。

      表 4  发射药性能参数预测

      Table 4.  Prediction of propellants property parameters

      propellant explosive force/(kJ·kg-1) explosive heat/(kJ·kg-1) explosive temperature/K
      measurement prediction error measurement prediction error calculation prediction error
      single base 1028 1000 2.7% 3876 3636 6.2% 2997 2857 4.7%
      double base 931.0 954.0 2.5% 3188 3388 6.3% 2580 2697 4.5%
      triple base 1226 1231 0.4% 4812 4853 0.9% 3618 3642 0.7%

      表 5  发射药性能参数预测(修正后)

      Table 5.  Prediction of propellants property parameters (after correction)

      propellant explosive force/(kJ·kg-1) explosive heat/(kJ·kg-1) explosive temperature/K
      measurement prediction error measurement prediction error calculation prediction error
      single base 1028 1044 1.6% 3876 3863 0.3% 2997 3004 0.2%
      double base 931.0 921.0 1.1% 3188 3196 0.3% 2580 2576 0.2%
      triple base 1226 1219 0.6% 4812 4818 0.1% 3618 3615 0.1%
    • 建立了脉冲激光与含能材料作用平台及冲击波、等离子体图像等瞬态诊断的实验平台,基于该平台研究了单基发射药和单晶冰糖在不同激光能量和光斑条件下的冲击波及等离子体演化规律。冲击波研究方面,相同激光能量下,对于单晶冰糖,更大的光斑导致冲击波减弱,这是由等离子体减弱引起,对于单基发射药,更大的光斑导致冲击波增强,这是由于含能物质的激发量增大引起的。等离子体图像方面,与单晶冰糖相比,单基发射药起始等离子体更强,但寿命更短,这是由于等离子体能量向含能材料反应区域转移引起的。可见含能材料与非含能材料脉冲激光作用过程冲击波及等离子体特性的主要影响因素不同,含能材料主要受含能物质化学反应的影响,非含能材料主要受等离子体强度影响。

      在激光参数条件为400 mJ/1474 μm条件下获取3种发射药冲击波传播距离,并同时监测每次激光能量。利用激光能量参数及SVM算法修正冲击波传播距离,基于LASEM方法建立的线性定标模型,对发射药火药力、爆热及爆温有较好的预测能力,R2分别达到0.9912、0.9998、0.9999,均优于未经激光能量修正的定标模型。研究验证了脉冲激光作用参数用于评估发射药能量性能参数的适用性。下一步将针对更多种类的含能材料开展研究。

参考文献 (20)

目录

    /

    返回文章
    返回