双模成型装药战斗部药型罩结构的正交设计

王彦盛 李伟兵 黄炫宁 王晓鸣

引用本文:
Citation:

双模成型装药战斗部药型罩结构的正交设计

    作者简介: 王彦盛(1998-),男,本科,主要从事弹药研究. E-mail:514270553@qq.com;
    通讯作者: 李伟兵, njustlwb@163.com
  • 中图分类号: O385; TJ410.33

Orthogonal Design of the Liner Structure in Dual-Mode Charge Warhead

    Corresponding author: LI Weibing, njustlwb@163.com ;
  • CLC number: O385; TJ410.33

  • 摘要: 针对双模战斗部小型化设计及其应用问题,利用LS-DYNA仿真软件,研究了双模战斗部药型罩结构参数(药型罩锥角、药型罩壁厚)对双模毁伤元成型性能的影响规律,揭示了各结构参数对双模毁伤元成型性能的控制规律:随着药型罩锥角和壁厚的增大,双模毁伤元的头部速度下降明显。结合正交设计方法,确定了药型罩壁厚是决定两毁伤元头部速度差的主要因素,锥角是决定各毁伤元头部速度的主要因素。得到了双模毁伤元成型性能均较佳的药型罩结构参数组合:药型罩锥角为80°,药型罩上端壁厚为5.0 mm,药型罩下端壁厚为4.0 mm,药型罩倒角弧度半径为10.0 mm。为验证模拟结果,进行了X射线成像试验,数值模拟结果与试验结果吻合较好。研究结果可为双模战斗部的进一步优化设计提供参考依据。
  • 图 1  成型装药结构示意图

    Figure 1.  Schematic of shaped charge structure

    图 2  JPC(a)和JET(b)的成型图

    Figure 2.  JPC (a) and JET (b) molding charts

    图 3  毁伤元头部速度随药型罩锥角的变化曲线(a)及毁伤元成型图像(b)

    Figure 3.  Change curve of the head velocity of the damaged element with the cone angle of the liner (a) and the corresponding damage molding chart (b)

    图 4  毁伤元头部速度与药型罩壁厚关系曲线及毁伤元成型图

    Figure 4.  Relationship between head velocity of damage element and liner thickness and the corresponding damage molding charts

    图 5  毁伤元头部速度与药型罩变壁厚关系曲线及毁伤元成型图

    Figure 5.  Relationship between head velocity of damage element and varied thickness of liner and the corresponding damage molding charts

    图 6  双模毁伤元各指标随因素水平的变化曲线

    Figure 6.  Change curves of each index of bimodal damage element with the level of factors

    图 7  毁伤元成型形态仿真(上)与试验结果(下)的对比

    Figure 7.  Comparison between simulation results (above) and X-ray pictures (below) of penetrators

    表 1  材料参数及计算模型

    Table 1.  Material parameters and calculation models

    ComponentMaterialρ/(g·cm−3Material modelEquation of state
    LinerCopper8.960Johnson_CookGrüneisen
    Explosive87011.695High_Explosive_BurnJWL
    Air1.250 × 10−3NullGrüneisen
    下载: 导出CSV

    表 2  药型罩结构参数仿真方案

    Table 2.  Simulation scheme of the structural parameters of the liner

    ProjectDe/mmL/mmα/(°)h/mmh1/mmh2/mmR/mmInitiation mode
    1100.00.90De70−1000.051De0.10De1, 2
    2100.00.90De85, 900.030De−0.070De0.10De1, 2
    3100.00.90De850.035De−0.055De0.055De0.10De1, 2
    950.050De0.050De−0.030De
    下载: 导出CSV

    表 3  正交试验的各因素水平

    Table 3.  Factor levels in orthogonal test

    LevelFactor
    α/(°)h1h2R
    1800.03De0.03De0
    2850.04De0.04De0.05De
    3900.05De0.05De0.10De
    4950.06De0.06De0.15De
    51000.07De0.07De0.20De
    下载: 导出CSV

    表 4  正交阵列各方案的计算结果

    Table 4.  Calculation results of orthogonal array schemes

    ProjectLαLh1Lh2LRvtip1/(m·s−1)vtip2/(m·s−1)Δv/(m·s−1)
    111115 2118 2893 078
    212224 7988 6833 885
    313334 4367 3842 948
    414444 1266 5092 383
    515553 8364 9921 156
    621234 6598 3613 702
    722344 3206 9002 580
    823454 0005 3531 353
    924513 7416 7413 000
    1025124 4497 1932 744
    1131354 2205 9821 762
    1232413 9326 9443 012
    1333523 6486 6292 981
    1434134 3666 4412 075
    1535243 9805 3561 376
    1641423 8817 2743 393
    1742533 5826 0352 453
    1843144 3215 7681 447
    1944253 8944 9571 063
    2045313 5706 1402 570
    2151543 5195 4391 920
    2252154 2855 5861 301
    2353213 6486 3092 661
    2454323 5105 9052 395
    2555433 2175 1441 927
    下载: 导出CSV

    表 5  各指标极差

    Table 5.  Range of indicators

    FactorSvtip1/(m·s−1)Svtip2/(m·s−1)SΔv/(m·s−1)
    α845.61 494.8649.2
    h1487.61 304.0816.4
    h2861.2766.0408.4
    R36.81 762.81 752.6
    下载: 导出CSV
  • [1] 郭美芳, 范宁军. 多模式战斗部与起爆技术分析研究 [J]. 探测与控制学报, 2005, 27(1): 31–34. doi: 10.3969/j.issn.1008-1194.2005.01.009
    GUO M F, FAN N J. The study on a multimode warhead and the initiation technology [J]. Journal of Detection and Control, 2005, 27(1): 31–34. doi: 10.3969/j.issn.1008-1194.2005.01.009
    [2] 谭多望, 孙承纬. 成型装药研究新进展 [J]. 爆炸与冲击, 2008, 28(1): 50–56. doi: 10.3321/j.issn:1001-1455.2008.01.009
    TAN D W, SUN C W. Progress in studies on shaped charge [J]. Explosion and Shock Waves, 2008, 28(1): 50–56. doi: 10.3321/j.issn:1001-1455.2008.01.009
    [3] BENDER D, FONG R, NG W, et al. Dual mode warhead technology for future smart munitions [C]//19th International Symposium on Ballistics. Interlaken, Switzerland: Casino Kursaal Congress Center, 2001: 679–684.
    [4] ARNOLD W, GRASWALD M, ROTTENKOLBER E. A novel technology for switchable modes warheads [C]//26th International Symposium on Ballistics. Miami, Florida, USA, 2011.
    [5] LI W B, WANG X, LI W B, et al. Research on the optimum length-diameter ratio of the charge of a multimode warhead [J]. Shock Waves, 2012, 22(3): 265–274. doi: 10.1007/s00193-012-0365-z
    [6] 纪冲, 龙源, 余道强, 等. 切割式双模战斗部毁伤元成型及侵彻钢靶特性研究 [J]. 高压物理学报, 2012, 26(5): 508–516. doi: 10.11858/gywlxb.2012.05.005
    JI C, LONG Y, YU D Q, et al. Experimental and numerical study on the formation and penetration properties of dual-mode warhead [J]. Chinese Journal of High Pressure Physics, 2012, 26(5): 508–516. doi: 10.11858/gywlxb.2012.05.005
    [7] 杨亚东, 陈智刚. 起爆方式对杆式射流形成影响的数值模拟研究 [J]. 弹箭与制导学报, 2008, 28(5): 89–91. doi: 10.3969/j.issn.1673-9728.2008.05.027
    YANG Y D, CHEN Z G. Numerical simulation of jetting penetrator charge formation under different initiation ways [J]. Journal of Projectiles, Rocket, Missiles and Guidance, 2008, 28(5): 89–91. doi: 10.3969/j.issn.1673-9728.2008.05.027
    [8] 樊菲. 实现杆流与射流转换的研究 [D]. 南京: 南京理工大学, 2012.
    FAN F. Research on achieving the conversion of JPC and JET [D]. Nanjing: Nanjing University of Science and Technology, 2012.
    [9] 李伟兵, 樊菲, 王晓鸣, 等. 杆式射流与射流转换的双模战斗部优化设计 [J]. 兵工学报, 2013, 34(12): 1500–1505. doi: 10.3969/j.issn.1000-1093.2013.12.003
    LI W B, FAN F, WANG X M, et al. Optimal design of dual mode warhead for jetting penetrator and jet conversion [J]. Acta Armamentarii, 2013, 34(12): 1500–1505. doi: 10.3969/j.issn.1000-1093.2013.12.003
    [10] 樊雪飞, 李伟兵, 王晓鸣, 等. 爆轰驱动钽药型罩形成双模毁伤元仿真与试验研究 [J]. 兵工学报, 2017, 38(10): 1918–1925. doi: 10.3969/j.issn.1000-1093.2017.10.006
    FAN X F, LI W B, WANG X M, et al. Simulation and experimental study of tantalum liner to form dual-mode damage element by detonation [J]. Acta Armamentraii, 2017, 38(10): 1918–1925. doi: 10.3969/j.issn.1000-1093.2017.10.006
    [11] FU J P, CHEN Z G, HOU X C, et al. Simulation and experimental investigation of jetting penetrator charge at large stand-off distance [J]. Defence Technology, 2013, 9(2): 117–124.
    [12] LI W B, WANG X M, LI W B. The effect of annular multi-point initiation on the formation and penetration of an explosively formed penetrator [J]. International Journal of Impact Engineering, 2010, 37(4): 414–424. doi: 10.1016/j.ijimpeng.2009.08.008
  • [1] 刘宏杰王伟力苗润吴世永 . 环形双锥罩聚能装药结构优化设计. 高压物理学报, 2018, 32(6): 065105-1-065105-9. doi: 10.11858/gywlxb.20180539
    [2] 于川王伟陈浩虞德水谢刚张振涛 . 小口径药型罩爆炸成型弹丸设计与多层钢靶侵彻实验. 高压物理学报, 2014, 28(1): 69-72. doi: 10.11858/gywlxb.2014.01.011
    [3] 郭腾飞李伟兵李文彬洪晓文 . 钽罩结构参数对EFP成型及侵彻性能的控制. 高压物理学报, 2018, 32(3): 035104-1-035104-8. doi: 10.11858/gywlxb.20170667
    [4] 王铁福王雷阮文俊赵同虎 . 药型罩材料的晶粒度对射流性能的影响. 高压物理学报, 1996, 10(4): 291-298 . doi: 10.11858/gywlxb.1996.04.009
    [5] 赵长啸龙源纪冲李裕春谢全民路亮 . 多点起爆下药型罩表面压力分布规律研究. 高压物理学报, 2013, 27(1): 83-89. doi: 10.11858/gywlxb.2013.01.012
    [6] 高永宏庞瑞强王凤英顾晓辉 . 材料颗粒度对粉末药型罩特性的影响. 高压物理学报, 2014, 28(4): 455-460. doi: 10.11858/gywlxb.2014.04.011
    [7] 郝彦君刘金旭李建崇冯伟刘旭真初东山李树奎 . 钛合金药型罩对纯铜及45钢靶的破甲特性. 高压物理学报, 2017, 31(5): 535-540. doi: 10.11858/gywlxb.2017.05.005
    [8] 龚柏林李明初哲吴海军柯明 . 贫铀合金药型罩聚能破甲性能实验研究. 高压物理学报, 2018, 32(3): 035102-1-035102-6. doi: 10.11858/gywlxb.20170647
    [9] 郑宇王晓鸣李文彬李伟兵 . 曲率半径对双层球缺罩形成串联爆炸成型弹丸的影响. 高压物理学报, 2009, 23(3): 229-235 . doi: 10.11858/gywlxb.2009.03.011
    [10] 唐蜜柏劲松李平张展冀 . 超高速撞击中影响碎片云形状因素分析. 高压物理学报, 2007, 21(4): 425-432 . doi: 10.11858/gywlxb.2007.04.016
    [11] 陈闯李伟兵王晓鸣李文彬吴巍 . 串联战斗部前级K装药结构的优化设计. 高压物理学报, 2014, 28(1): 73-78. doi: 10.11858/gywlxb.2014.01.012
    [12] 邓爱东张华杨显俊 . 电磁驱动产生超强磁场的参数优化设计. 高压物理学报, 2015, 29(2): 123-128. doi: 10.11858/gywlxb.2015.02.006
    [13] 龚柏林卢芳云李翔宇 . 可变形战斗部弹体变形型面研究. 高压物理学报, 2010, 24(2): 102-106 . doi: 10.11858/gywlxb.2010.02.004
    [14] 纪冲龙源余道强周翔张洋溢 . 切割式双模战斗部毁伤元成型及侵彻钢靶特性研究. 高压物理学报, 2012, 26(5): 508-516. doi: 10.11858/gywlxb.2012.05.005
    [15] 段建杨黔龙周刚初哲田亚军张颖 . 串联战斗部前级聚能装药和隔爆结构设计与实验研究. 高压物理学报, 2006, 20(2): 202-206 . doi: 10.11858/gywlxb.2006.02.015
    [16] 李如江沈兆武王凤英 . 药型罩初始孔隙度对聚能射流行为的影响. 高压物理学报, 2008, 22(4): 445-448 . doi: 10.11858/gywlxb.2008.04.019
    [17] 孙建军李如江杨玥万清华张明孙淼李杨 . 辅助药型罩材料对线型聚能射流性能影响的数值仿真. 高压物理学报, 2018, 32(6): 065106-1-065106-7. doi: 10.11858/gywlxb.20180542
    [18] 王铁福赵国虎阮文俊王雷 . 贫铀药型罩材料中微量杂质对聚能射流性能的影响. 高压物理学报, 1996, 10(2): 102-106 . doi: 10.11858/gywlxb.1996.02.004
    [19] 庞勇于川桂毓林 . 球缺药型罩爆炸成型弹丸数值模拟. 高压物理学报, 2005, 19(1): 86-92 . doi: 10.11858/gywlxb.2005.01.015
    [20] 刘志跃翟俊昭 . 铝铜药型罩射流与侵彻数值模拟. 高压物理学报, 2019, 33(6): 064107-1-064107-8. doi: 10.11858/gywlxb.20190728
  • 加载中
图(7)表(5)
计量
  • 文章访问数:  515
  • 阅读全文浏览量:  380
  • PDF下载量:  1
出版历程
  • 收稿日期:  2020-04-07
  • 录用日期:  2020-04-26
  • 网络出版日期:  2020-06-25
  • 刊出日期:  2020-08-05

双模成型装药战斗部药型罩结构的正交设计

    作者简介:王彦盛(1998-),男,本科,主要从事弹药研究. E-mail:514270553@qq.com
    通讯作者: 李伟兵, njustlwb@163.com
  • 南京理工大学智能弹药技术国防重点学科实验室,江苏 南京 210094

摘要: 针对双模战斗部小型化设计及其应用问题,利用LS-DYNA仿真软件,研究了双模战斗部药型罩结构参数(药型罩锥角、药型罩壁厚)对双模毁伤元成型性能的影响规律,揭示了各结构参数对双模毁伤元成型性能的控制规律:随着药型罩锥角和壁厚的增大,双模毁伤元的头部速度下降明显。结合正交设计方法,确定了药型罩壁厚是决定两毁伤元头部速度差的主要因素,锥角是决定各毁伤元头部速度的主要因素。得到了双模毁伤元成型性能均较佳的药型罩结构参数组合:药型罩锥角为80°,药型罩上端壁厚为5.0 mm,药型罩下端壁厚为4.0 mm,药型罩倒角弧度半径为10.0 mm。为验证模拟结果,进行了X射线成像试验,数值模拟结果与试验结果吻合较好。研究结果可为双模战斗部的进一步优化设计提供参考依据。

English Abstract

  • 为了适应现代战场复杂多变的作战环境,武器系统正朝着信息化、精确化和多功能化发展。这就要求武器系统能够自动识别目标类型,从而自适应地在多种模态中选择特定作战模式,以此打击多种目标。在众多弹药系统中,多模战斗部[1]是当前战斗部发展的主要方向之一。它是在同一成型装药[2]的基础上采用不同起爆方式来实现射流(JET)、杆式射流(JPC)、爆炸成型弹丸(EFP)或破片等多种毁伤模式的转换[3]。2001年,Bender等[3]率先开展了多模技术的基础性研究。近年来,Arnold等[4]发展了一种新型切换模式战斗部技术。对此我国很多学者也开展了大量研究,例如:Li等[5]研究了聚能装药长径比对多模毁伤元成型参数的影响;纪冲等[6]对同一成型装药结构实现EFP和多爆炸成型弹丸(MEFP)的双模转换战斗部开展了研究;杨亚东等[7]、樊菲[8]、李伟兵等[9]对于同一成型装药结构实现JPC与JET的双模转换进行了研究。然而,对于成型装药药型罩上、下壁厚变化对JPC和JET转换的影响,以及将药型罩上、下壁厚与其他药型罩结构参数进行匹配设计的研究工作,国内外较少见诸报道。

    本研究利用LS-DYNA有限元软件对JPC与JET的转换过程进行数值模拟,分析药型罩锥角和壁厚对JPC与JET转化的影响,并对JPC与JET转换的双模战斗部成型装药结构参数进行正交设计,对数值模拟结果进行分析,结合成型装药结构优化指标(毁伤元头部速度、毁伤元成型状态),获得最佳药型罩结构参数,并进行双模毁伤元X射线成像试验验证,为双模转换战斗部结构的进一步优化设计提供参考。

    • 本研究的成型装药结构如图1所示,其装药口径De为100.0 mm,药型罩采用大锥角罩,顶部倒角,其结构参数包括倒角弧度半径R、壁厚h、上壁厚h1、下壁厚h2。通过改变药型罩结构参数设计方案,分析不同药型罩结构参数对聚能毁伤元成型的影响规律。成型装药为船尾形,装药高度L = 0.90De。根据成型装药的相关研究结果[10]可知:装药顶端中心单点起爆(起爆方式1)形成JPC毁伤元;装药船尾底端环起爆(起爆方式2)形成JET毁伤元,其中起爆环距离装药顶点56.5 mm。JPC头部速度一般为3~5 km/s,JPC成型仿真结果见图2(a);JET头部速度一般为5~10 km/s,JET成型仿真结果见图2(b)[11]

      图  1  成型装药结构示意图

      Figure 1.  Schematic of shaped charge structure

      图  2  JPC(a)和JET(b)的成型图

      Figure 2.  JPC (a) and JET (b) molding charts

      聚能装药的作用过程是多物质相互作用的大变形运动过程,在金属射流形成过程中炸药和药型罩材料会发生愈来愈剧烈的变形,采用Lagrange方法难以准确模拟。因此,本研究采用LS-DYNA的多物质ALE(Arbitrary Lagrange-Euler)方法和运动网格法,模拟JET和JPC的形成、延展和断裂。对于多物质ALE方法,除了聚能装置外,还需要建立足以覆盖整个射流范围的空气网格,并且在模型边界施加压力流出边界条件,避免压力在边界上反射。计算中各部分的材料参数及计算模型[9, 12]表1所示,其中ρ为密度。

      ComponentMaterialρ/(g·cm−3Material modelEquation of state
      LinerCopper8.960Johnson_CookGrüneisen
      Explosive87011.695High_Explosive_BurnJWL
      Air1.250 × 10−3NullGrüneisen

      表 1  材料参数及计算模型

      Table 1.  Material parameters and calculation models

    • 药型罩是影响聚能毁伤元成型效果和成型质量的最关键部件。药型罩结构参数直接影响毁伤元成型的形态。为了寻找双模战斗部正交设计中正交表所需的药型罩结构参数的取值范围,首先研究药型罩结构参数对双模毁伤元成型的影响规律。仿真计算方案如表2所示,分别分析大锥角罩的锥角α、壁厚h、上壁厚h1和下壁厚h2对毁伤元转换的影响(研究变壁厚时,h为变量,故变壁厚方案中h未取值)。

      ProjectDe/mmL/mmα/(°)h/mmh1/mmh2/mmR/mmInitiation mode
      1100.00.90De70−1000.051De0.10De1, 2
      2100.00.90De85, 900.030De−0.070De0.10De1, 2
      3100.00.90De850.035De−0.055De0.055De0.10De1, 2
      950.050De0.050De−0.030De

      表 2  药型罩结构参数仿真方案

      Table 2.  Simulation scheme of the structural parameters of the liner

    • 针对方案1,研究70°~100°锥角下药型罩的成型情况(每隔5°为一个工况)。选取毁伤元成型时刻为100 μs时的结果,获得了毁伤元头部速度vtip随药型罩锥角的变化曲线,如图3(a)所示。以药型罩锥角α = 90°工况为例,毁伤元成型图见图3(b)。从图3(a)中可以得出:在装药顶端中心单点起爆条件下,药型罩锥角α越大,毁伤元头部速度越小,降幅为1317 m/s;在装药船尾底端环起爆条件下,毁伤元头部速度也随α的增大而减小,降幅为3280 m/s;随着α的增大,两种起爆方式下毁伤元的头部速度差Δv越来越小,从70°时的3879 m/s减小到100°时的1916 m/s。

      图  3  毁伤元头部速度随药型罩锥角的变化曲线(a)及毁伤元成型图像(b)

      Figure 3.  Change curve of the head velocity of the damaged element with the cone angle of the liner (a) and the corresponding damage molding chart (b)

      JPC头部速度为3~5 km/s,结合图3得出,对于α < 80°的情况可不作进一步研究。随着α的增大,两种起爆方式下的毁伤元头部越来越接近,不利于实现JPC与JET两种模态的转换,在此对α > 100°的情况也不作进一步研究。

    • 针对方案2,研究0.030De~0.070De药型罩壁厚下药型罩的成型情况(每隔1 mm为一个工况)。为减少计算量,只讨论85°和90°锥角下药型罩的壁厚变化。选取毁伤元成型时刻为100 μs时的结果,获得了毁伤元头部速度与药型罩壁厚的关系曲线,如图4所示,其中右侧插图为85°和90°锥角下药型罩壁厚h = 5.0 mm时的毁伤元成型图。

      图  4  毁伤元头部速度与药型罩壁厚关系曲线及毁伤元成型图

      Figure 4.  Relationship between head velocity of damage element and liner thickness and the corresponding damage molding charts

      图4可以看出,药型罩壁厚越大,两种起爆方式下毁伤元头部速度差Δv越小,并逐渐趋于平稳。同时,可以观察到85°锥角下的头部速度差Δv变化较剧烈。对于85°和90°锥角的药型罩,当药型罩壁厚为0.050De时,速度差Δv趋于平稳。

    • 针对方案3,研究在上壁厚和下壁厚其中一个量不变而另一量改变时药型罩的成型情况。为减少工作量,只研究锥角为95°、h2 = 5.0 mm条件下上壁厚在0.030De~0.050De范围(每隔0.5 mm为一个工况)和锥角为85°、h1 = 5.5 mm条件下下壁厚在0.035De~0.055De范围(每隔0.5 mm为一个工况)的头部速度变化规律。选取毁伤元成型时刻为100 μs时的结果,分析得到毁伤元头部速度与药型罩变壁厚关系曲线,如图5所示。以α = 95°、h1 = 5.0 mm、h2 = 5.0 mm以及α = 85°、h1 = 5.5 mm、h2 = 5.5 mm两种工况为例,计算得到双模毁伤元成型图,见图5右侧插图。在上壁厚不变的条件下,当下壁厚从0.050De降至0.030De时,Δv降幅为122 m/s,下壁厚减小使两种起爆方式下毁伤元的头部速度差Δv降低。在下壁厚不变的条件下,当上壁厚从0.055De降至0.035De时,Δv增幅为688 m/s,上壁厚减小使两种起爆方式下毁伤元的头部速度差Δv增大,且上壁厚较下壁厚对毁伤元头部速度差的影响更大。

      图  5  毁伤元头部速度与药型罩变壁厚关系曲线及毁伤元成型图

      Figure 5.  Relationship between head velocity of damage element and varied thickness of liner and the corresponding damage molding charts

    • 正交设计是通过正交表安排多因素试验、利用统计数学原理进行数据分析的一种科学方法,符合“以尽量少的试验获得有效、足够的信息”的试验设计原则。正交表能够在因素变化范围内均衡抽样,使每次试验都具有较强的代表性。由于正交表具备均衡分散的特点,保证了全面试验的某些要求,这些试验能够较好或更好地达到试验目的。如在4因素5水平的条件下,根据L25正交表进行25次仿真计算,即可完成组合计算54 = 625次试验的内容,方便找到各个因素对最终目标的影响规律。

    • 在双模战斗部装药高度及各部分材料确定的条件下,选择药型罩锥角α、药型罩上壁厚h1、药型罩下壁厚h2、药型罩倒角弧度半径R[7]这4个结构参数作为正交设计的4个因素,每个因素取5个水平,参与正交计算,每个因素及其对应的水平见表3

      LevelFactor
      α/(°)h1h2R
      1800.03De0.03De0
      2850.04De0.04De0.05De
      3900.05De0.05De0.10De
      4950.06De0.06De0.15De
      51000.07De0.07De0.20De

      表 3  正交试验的各因素水平

      Table 3.  Factor levels in orthogonal test

      将装药顶端单点起爆与装药船尾底端环起爆两种起爆方式下毁伤元成型后的头部速度vtip1vtip2,以及两种起爆方式下毁伤元成型后的头部速度差Δv作为评价指标。由于JPC与JET的头部速度范围较大(JET的头部速度一般可以达到5~10 km/s,JPC的头部速度在3~5 km/s之间),且两种模态的头部速度差越大越好,即JET的头部速度尽可能大,而JPC的头部速度尽可能接近3 km/s,因此将形成最佳转换的药型罩结构条件约束为较小的JPC头部速度及较大的JET头部速度。各组合的计算结果如表4所示,其中LαLh1Lh2LR分别为αh1h2R的水平数。

      ProjectLαLh1Lh2LRvtip1/(m·s−1)vtip2/(m·s−1)Δv/(m·s−1)
      111115 2118 2893 078
      212224 7988 6833 885
      313334 4367 3842 948
      414444 1266 5092 383
      515553 8364 9921 156
      621234 6598 3613 702
      722344 3206 9002 580
      823454 0005 3531 353
      924513 7416 7413 000
      1025124 4497 1932 744
      1131354 2205 9821 762
      1232413 9326 9443 012
      1333523 6486 6292 981
      1434134 3666 4412 075
      1535243 9805 3561 376
      1641423 8817 2743 393
      1742533 5826 0352 453
      1843144 3215 7681 447
      1944253 8944 9571 063
      2045313 5706 1402 570
      2151543 5195 4391 920
      2252154 2855 5861 301
      2353213 6486 3092 661
      2454323 5105 9052 395
      2555433 2175 1441 927

      表 4  正交阵列各方案的计算结果

      Table 4.  Calculation results of orthogonal array schemes

      利用极差分析方法对25次仿真计算结果进行分析。将各列水平数相同的结果相加,记为K,5个水平的计算结果分别为K1K2K3K4K5,分别除以5,优化出每个因素的水平;将5个水平中的最大值减去最小值,得到极差S,通过S可以得到各因素对指标影响的主次顺序。表5列出了各因素影响下对应各个指标的极差S,即Svtip1Svtip2SΔv

      FactorSvtip1/(m·s−1)Svtip2/(m·s−1)SΔv/(m·s−1)
      α845.61 494.8649.2
      h1487.61 304.0816.4
      h2861.2766.0408.4
      R36.81 762.81 752.6

      表 5  各指标极差

      Table 5.  Range of indicators

      表5可以看出各因素对各个指标影响的主次顺序:对vtip1的影响因素由主到次依次为h2αh1R,对vtip2的影响因素由主到次依次为Rαh1h2,对Δv的影响因素由主到次依次为Rh1αh2,其中药型罩上下壁厚h1h2是影响各个指标的重要因素。为了分析每个因素中各个水平对3个指标的影响情况,计算得到了头部速度vtip1vtip2和头部速度差Δv随不同因素水平的变化关系,如图6所示。图6清晰地反映了各个因素对3个评价因素的影响规律,从中可以得到不同因素对同一指标的影响差别。

      图  6  双模毁伤元各指标随因素水平的变化曲线

      Figure 6.  Change curves of each index of bimodal damage element with the level of factors

      根据表5图6的结果,进行各因素的最优水平组合。h1是影响Δv的第2重要因素,h1越小,Δv越大。h2是影响JPC头部速度的最主要因素,h2越大,JPC毁伤元头部速度越小,对JET头部速度以及Δv的影响较小。双模战斗部装药结构的要求:在射流尽可能不断裂的情况下,JET的头部速度越大越好,JPC和JET的头部速度差尽可能大。综合各指标要求,选择h1 = 0.05Deh2 = 0.04De,使JET既有较大的头部速度,又有较大的头部速度差。药型罩倒角弧度半径R是影响JET头部速度的最主要因素,但对JPC头部速度几乎没有影响。从仿真结果来看,R越小,JET头部速度越大,头部速度差越大,侵彻体越长。从成型效果图来看,R越小,装药船尾底端环起爆形成的JET侵彻体越细,杵体越大,不利于侵彻。综合考虑后,选择R = 0.10De。药型罩锥角α是影响两种毁伤元头部速度和头部速度差的重要因素,α越小,毁伤元头部速度越大,毁伤元头部速度差也越大。权衡考虑,选择α = 80°。

    • 基于正交设计计算结果,选取药型罩结构参数:α = 80°,h1 = 0.05Deh2 = 0.04DeR = 0.10De。由于此工况不在正交表中,因此需要按照新优化方案进行计算,结果如图7所示。基于图1所示成型装药结构,采用优化方案,开展X射线成像试验。试验仪器主要包括成型装药、托弹架、X射线管和脉冲X光机、底片、底片保护盒等。试验过程中,将成型装药固定在一定高度的托弹架上,通过控制X射线管的出光时间,便可在底片上得到毁伤元的X射线成像照片。根据所拍摄的时间及不同毁伤元来布置靶块及炸高筒,每次靶块及炸高筒的总高度由仿真结果确定,保证毁伤元在拍摄时间通过底片盒。

      图  7  毁伤元成型形态仿真(上)与试验结果(下)的对比

      Figure 7.  Comparison between simulation results (above) and X-ray pictures (below) of penetrators

      试验获得了两个时刻装药顶端中心单点起爆和装药船尾底端环形多点起爆下毁伤元成型照片,选取对应时刻的仿真计算结果进行对比,如图7所示。考虑到JPC和JET的长度较长,需要用两张底片来获得完整的毁伤元成型形状,而试验中只有3个底片盒,为此在第1时刻选取两张底片拍摄完整的毁伤元形态,第2时刻用1张底片拍摄毁伤元部分形态。因为第1时刻射流尾部未完全进入底片,因此只能将两个时刻头部平均速度与仿真结果进行对比。

      对毁伤元X射线照片进行数字化处理,得到了优化方案形成的毁伤元两个时刻头部速度的平均值。其中,试验获得的JPC毁伤元头部平均速度为3901 m/s,仿真计算得到的头部速度为4006 m/s;试验获得的JPC毁伤元在100 μs时的长度为353.3 mm,仿真值为351.5 mm。因底片问题,只获得了两个时刻JET毁伤元头部平均速度,为5638 m/s,相应的仿真结果为5925 m/s。从毁伤元成型结果的试验和仿真对比来看,X射线拍摄的毁伤元成型形态与仿真结果具有较好的一致性,JPC毁伤元成型参数的试验与仿真结果的相对偏差不超过5%,JET毁伤元头部平均速度的试验结果与仿真结果的相对偏差不超过10%。

    • (1) 通过改变起爆方式,设计并实现了双模成型装药战斗部中JPC与JET的双模转换,其中装药顶端中心单点起爆时形成了JPC,装药船尾底端环起爆时形成了JET。获得了药型罩结构参数对双模毁伤元成型的影响规律:随着药型罩锥角的增大,两种起爆方式下形成的毁伤元头部速度差Δv逐渐减小;随着药型罩壁厚增大,Δv变小,且逐渐趋于平稳;药型罩上壁厚不变、下壁厚变小时,Δv基本不变;药型罩下壁厚不变、上壁厚变小时,Δv变大。由此找出了双模毁伤元成型较佳时各参数的取值范围:药型罩锥角α为80°~100°,药型罩上壁厚h1为0.030De~0.070De,药型罩下壁厚h2为0.030De~0.070De

      (2) 通过正交设计方法,得到了药型罩各结构参数在两种起爆方式下对毁伤元头部速度和头部速度差影响的主次顺序,其中:对装药顶端中心单点起爆下毁伤元头部速度的影响因素由主到次依次为下壁厚、锥角、上壁厚、倒角弧度半径,对装药船尾底端环起爆下毁伤元头部速度的影响因素由主到次依次为倒角弧度半径、锥角、上壁厚、下壁厚,对Δv的影响因素由主到次依次为倒角弧度半径、上壁厚、锥角、下壁厚。

      (3) 通过正交设计得到了双模毁伤元成型性能较佳的药型罩结构参数组合,即药型罩锥角为80°,药型罩上端壁厚为5.0 mm,药型罩下端壁厚为4.0 mm,药型罩倒角弧度半径为10.0 mm。在此装药结构下,JPC和JET的成型效果都较好,试验结果与仿真计算结果较吻合。

参考文献 (12)

目录

    /

    返回文章
    返回