串联战斗部不同介质组合的隔爆能力

刘宏杰 王伟力 苗润 吴世永 王俊华

引用本文:
Citation:

串联战斗部不同介质组合的隔爆能力

    作者简介: 刘宏杰(1993-),男,硕士研究生,主要从事弹药设计与目标毁伤评估研究. E-mail: 1300202650@qq.com;
    通讯作者: 王伟力, 858769430@qq.com
  • 中图分类号: O382; TJ410.3

Explosive Interruption of Tandem Warhead with Different Multilayer Structures

    Corresponding author: WANG Weili, 858769430@qq.com ;
  • CLC number: O382; TJ410.3

  • 摘要: 为了有效提升串联切割战斗部隔爆结构衰减爆炸冲击波的性能,解决前级聚能装药结构与后级随进弹的匹配及隔爆问题,在前级切割器和后级随进弹之间加装隔爆结构,使用有限元分析软件ANSYS/LS-DYNA建立模型,进行不同组合结构隔爆性能的数值模拟,比较隔爆能力。模拟结果表明:前级装药起爆后,爆炸冲击波首先向后级随进弹头靠里区域汇聚,而不是向弹头尖端区域汇聚,因此可以适当减薄外层金属隔爆介质头部尖端区域;将外层金属由硬质钢改成铝时,后端壳体应力峰值的变化很小,故确定外层金属介质为铝;铝-聚脲的隔爆能力优于铝-泡沫铝结构,最终确定“软”隔爆介质为聚脲。通过调整铝和聚脲层的厚度,确定了最佳隔爆参数,能够满足实际应用。
  • 图 1  反射波与透射波Hugoniot关系计算曲线

    Figure 1.  Hugoniot curves of reflection and transmission wave

    图 2  串联战斗部及隔爆结构有限元模型

    Figure 2.  Finite element model of tandem warhead and explosive interruption structure

    图 3  工况I中装药压力及壳体应力分析

    Figure 3.  Loading pressure and shell stress of working condition I

    图 4  典型时刻外层钢介质的应力云图

    Figure 4.  Stress nephogram of steel at typical times

    图 5  工况II装药压力及壳体应力分析

    Figure 5.  Loading pressure and shell stress of working condition II

    图 6  工况III装药压力及壳体应力

    Figure 6.  Loading pressure and shell stress of working condition III

    图 7  工况IV装药压力及壳体应力

    Figure 7.  Loading pressure and shell stress of working condition IV

    图 8  工况V装药压力及壳体应力

    Figure 8.  Loading pressure and shell stress of working condition V

    图 9  后级随进弹加速度-时间变化曲线

    Figure 9.  Acceleration vs. time for the following stage’s incoming projectile

    图 10  后级随进弹速度-时间变化曲线

    Figure 10.  Velocity vs. time for the following stage’s incoming projectile

    图 11  串联战斗部及改进隔爆结构有限元模型

    Figure 11.  Finite element model of tandem warhead and improved explosive interruption structure

    图 12  工况I′装药压力及壳体应力

    Figure 12.  Loading pressure and shell stress of working condition I′

    图 13  工况II′装药压力及壳体应力分析

    Figure 13.  Loading pressure and shell stress of working condition II′

    图 14  工况III′装药压力及壳体应力分析

    Figure 14.  Loading pressure and shell stress of working condition III′

    图 15  工况IV′装药压力及壳体应力分析

    Figure 15.  Loading pressure and shell stress of working condition IV′

    图 16  后级随进弹加速度-时间变化曲线

    Figure 16.  Acceleration vs. time for the following stage’s incoming projectile

    图 17  后级随进弹加速度-时间曲线

    Figure 17.  Velocity vs. time for the following stage’s incoming projectile

    表 1  B炸药材料参数

    Table 1.  Material parameters of composition B

    ρ/(g·cm-3) D/(m·s-1) pCJ/GPa A/GPa B/GPa R1 R2 ω E0/GPa V0
    1.713 7 500 28.6 524.2 7.678 4.2 1.1 0.34 8.499 1.0
    下载: 导出CSV

    表 2  30CrMnSiNi2A钢和金属铝计算参数

    Table 2.  Material performance parameters of 30CrMnSiNi2A steel and aluminum

    Material AJC/MPa BJC/MPa n C m Tm/K T0/K S1 γ0 a
    Steel 1 280 420 0.30 0.030 1.00 1 793 294 1.490 2.17 0.46
    Al 265 426 0.34 0.015 1.00 775 294 1.345 2.13 0.10
    下载: 导出CSV

    表 3  不同工况及简便书写方式

    Table 3.  Different conditions and indications

    Condition Material
    A B C D
    Steel Polyurea Aluminum foam Aluminum
    Aluminum Polyurea Aluminum foam Aluminum
    Aluminum Aluminum foam Aluminum foam Aluminum
    Aluminum Aluminum foam Aluminum foam Aluminum foam
    Aluminum Polyurea Polyurea Polyurea
    下载: 导出CSV

    表 4  改进结构

    Table 4.  Improved structures

    Condition Material
    A B C
    Ⅰ′ Steel Aluminum Polyurea
    Ⅱ′ Steel Aluminum Aluminum foam
    Ⅲ′ Aluminum Aluminum Polyurea
    Ⅳ′ Aluminum Aluminum Polyurea(2 cm)
    下载: 导出CSV
  • [1] LEBLANC J, SHILLINGS C, GAUCH E, et al. Near field underwater explosion response of polyurea coated composite plates [J]. Experimental Mechanics, 2016, 56(4): 569–581. doi: 10.1007/s11340-015-0071-8
    [2] HARIS A, LEE H P, TAN V B C. An experimental study on shock wave mitigation capability of polyurea and shear thickening fluid based suspension pads [J]. Defence Technology, 2018, 14(1): 12–18. doi: 10.1016/j.dt.2017.08.004
    [3] DAI L, WU C, AN F, et al. Experimental investigation of polyurea-coated steel plates at underwater explosive loading [J]. Advances in Materials Science and Engineering, 2018: 1264276.
    [4] 曾必强, 姜春兰, 严翰新, 等. 串联攻坚战斗部前级爆轰场对随进弹随进影响分析 [J]. 兵工学报, 2010(Suppl 1): 162–166
    ZENG B Q, JIANG C L, YAN H X, et al. Analysis for effects of precursory detonation field on projectile following course in tandem warhead [J]. Acta Armamentarii, 2010(Suppl 1): 162–166
    [5] 姜夕博, 饶国宁, 徐森, 等. 冲击波在有机玻璃中衰减特性的数值模拟与实验研究 [J]. 南京理工大学学报(自然科学版), 2012, 36(6): 1059–1064 doi: 10.3969/j.issn.1005-9830.2012.06.028
    JIANG X B, RAO G N, XU S, et al. Numerical simulation and experimental research on shock wave attenuation properties in PMMA [J]. Journal of Nanjing University of Science and Technology, 2012, 36(6): 1059–1064 doi: 10.3969/j.issn.1005-9830.2012.06.028
    [6] 徐森, 刘大斌, 彭金华, 等. 药柱冲击波在有机玻璃中的衰减特性研究 [J]. 高压物理学报, 2010, 24(6): 431–437 doi: 10.11858/gywlxb.2010.06.005
    XU S, LIU D B, PENG J H, et al. Study on the shock wave attenuation of the booster charge in the PMMA gap [J]. Chinese Journal of High Pressure Physics, 2010, 24(6): 431–437 doi: 10.11858/gywlxb.2010.06.005
    [7] 侯海周, 彭金华, 胡毅亭. 爆炸冲击波在酚醛层压材料中衰减特性的实验研究 [J]. 火工品, 2016(2): 13–16 doi: 10.3969/j.issn.1003-1480.2016.02.005
    HOU H Z, PENG J H, HU Y T. Experimental study of shock wave attenuation properties in phenolic cotton fabric material [J]. Initiators & Pyrotechnics, 2016(2): 13–16 doi: 10.3969/j.issn.1003-1480.2016.02.005
    [8] TEDESCO J W, LANDIS D W. Wave propagation through layered systems [J]. Computers & Structures, 1989, 32(3/4): 625–638.
    [9] PETEL O E, JETTÉ F X, GOROSHIN S, et al. Blast wave attenuation through a composite of varying layer distribution [J]. Shock Waves, 2011, 21(3): 215–224. doi: 10.1007/s00193-010-0295-6
    [10] 石少卿, 刘仁辉, 汪敏. 钢板-泡沫铝-钢板新型复合结构降低爆炸冲击波性能研究 [J]. 振动与冲击, 2008, 27(4): 143–146 doi: 10.3969/j.issn.1000-3835.2008.04.037
    SHI S Q, LIU R H, WANG M. Shock wave reduction behavior of a new compound structure composed of a foam aluminum layer between two steel plates [J]. Journal of Vibration and Shock, 2008, 27(4): 143–146 doi: 10.3969/j.issn.1000-3835.2008.04.037
    [11] 董永香, 冯顺山, 李学林. 爆炸波在硬-软-硬三明治介质中传播特性的数值分析 [J]. 弹道学报, 2007, 19(1): 59–63 doi: 10.3969/j.issn.1004-499X.2007.01.017
    DONG Y X, FENG S S, LI X L. Numerical analysis of propagation characteristics of explosive wave in the hard-soft-hard sandwich media [J]. Journal of Ballistics, 2007, 19(1): 59–63 doi: 10.3969/j.issn.1004-499X.2007.01.017
    [12] 陈闯, 王晓鸣, 李文彬, 等. 多层介质阻抗匹配对隔爆效果的影响 [J]. 振动与冲击, 2014, 33(17): 105–110
    CHEN C, WANG X M, LI W B, et al. Influence of multilayered media impedance matching on explosion interruption effect [J]. Journal of Vibration and Shock, 2014, 33(17): 105–110
    [13] 北京工业学院八系. 爆炸及其作用 [M]. 北京: 国防工业出版社, 1994: 105-112.
    [14] 孙承纬. 爆炸物理学 [M]. 北京: 科学出版社, 2011: 205-207.
    [15] 宋博, 胡时胜, 王礼立. 分层材料的不同排列次序对透射冲击波强度的影响 [J]. 兵工学报, 2000, 21(3): 272–274 doi: 10.3321/j.issn:1000-1093.2000.03.021
    SONG B, HU S S, WANG L L. Influence on the transmitted intensity of shock wave through different tactic orders of layered materials [J]. Acta Armamentarii, 2000, 21(3): 272–274 doi: 10.3321/j.issn:1000-1093.2000.03.021
    [16] 汪文革, 杨世军, 韩永要, 等. 基于ANSYS/LS-DYNA的聚能射流侵彻装甲钢的有限元分析 [J]. 兵工自动化, 2008, 27(3): 39–41 doi: 10.3969/j.issn.1006-1576.2008.03.015
    WANG W G, YANG S J, HAN Y Y, et al. Finite element analysis of shaped charge jet penetrating into target based on ANSYS/LS-DYNA [J]. Ordnance Industry Automation, 2008, 27(3): 39–41 doi: 10.3969/j.issn.1006-1576.2008.03.015
    [17] TOGAMI T C, BAKER W E, FORRESTAL M J. A split Hopkinson bar technique to evaluate the performance of accelerometers [J]. Journal of Applied Mechanics, 1996, 63(2): 353–356. doi: 10.1115/1.2788872
    [18] FRANZEN R R, SCHNEIDEWIND P N. Observations concerning the penetration mechanics of tubular hypervelocity penetrators [J]. International Journal of Impact Engineering, 1991, 11(3): 289–303. doi: 10.1016/0734-743X(91)90040-M
    [19] WASMUND T L. New model to evaluate weapon effects and platform vulnerability: AJEM [J]. Wstiac Newsletter, 2001, 2: 1–3.
    [20] LEE W H, PAINTER J W. Material void-opening computation using particle method [J]. International Journal of Impact Engineering, 1999, 22(1): 1–22. doi: 10.1016/S0734-743X(98)00041-4
    [21] 甘云丹, 宋力, 杨黎明. 弹性体涂覆钢板抗冲击性能的数值模拟 [J]. 兵工学报, 2009(Suppl 2): 15–18
    GAN Y D, SONG L, YANG L M. Numerical simulation for anti-blast performances of steel plate coated with elastomer [J]. Acta Armamentarii, 2009(Suppl 2): 15–18
  • [1] 董永香黄晨光段祝平 . 多层介质对应力波传播特性影响分析. 高压物理学报, 2005, 19(1): 59-65 . doi: 10.11858/gywlxb.2005.01.011
    [2] 陈闯李伟兵王晓鸣李文彬吴巍 . 串联战斗部前级K装药结构的优化设计. 高压物理学报, 2014, 28(1): 73-78. doi: 10.11858/gywlxb.2014.01.012
    [3] 康彦龙蒋建伟王树有门建兵 . 面层混凝土对跑道内爆毁伤影响的仿真研究. 高压物理学报, 2012, 26(1): 76-82. doi: 10.11858/gywlxb.2012.01.011
    [4] 陈兴周兰伟李向东胡诤哲张瑞 . 破片式战斗部破片与冲击波相遇位置研究. 高压物理学报, 2018, 32(6): 065101-1-065101-9. doi: 10.11858/gywlxb.20180591
    [5] 徐浩铭顾文彬曾政王振雄赵长啸 . 延时起爆对串联EFP成型及其侵彻效应的影响. 高压物理学报, 2014, 28(1): 79-85. doi: 10.11858/gywlxb.2014.01.013
    [6] 杨佳喻寅蒋冬冬张明建李永宏刘福生贺红亮 . 低冲击应力下脆性介质电阻率的实时测量. 高压物理学报, 2012, 26(4): 389-394. doi: 10.11858/gywlxb.2012.04.005
    [7] 陈朗鲁建英冯长根 . 空气中激光支持爆轰波实验及理论分析. 高压物理学报, 2010, 24(2): 154-160 . doi: 10.11858/gywlxb.2010.02.012
    [8] 鲁建英陈朗伍俊英冯长根 . 激光支持等离子体爆轰波温度的实验测量. 高压物理学报, 2009, 23(2): 123-129 . doi: 10.11858/gywlxb.2009.02.008
    [9] 杨龙李平王刚华阚明先 . 固体套筒内爆非冲击压缩研究. 高压物理学报, 2016, 30(4): 344-352. doi: 10.11858/gywlxb.2016.04.012
    [10] 章冠人 . 冲击波作用下材料的极化率. 高压物理学报, 1990, 4(3): 161-166 . doi: 10.11858/gywlxb.1990.03.001
    [11] 金孝刚王弘 . 冲击波处理45钢的微结构研究. 高压物理学报, 1993, 7(4): 254-259 . doi: 10.11858/gywlxb.1993.04.003
    [12] 王洪波杨世全谢若泽钱立新卢永刚 . 混凝土靶体预开孔对弹丸侵彻性能的影响. 高压物理学报, 2015, 29(1): 69-74. doi: 10.11858/gywlxb.2015.01.012
    [13] 刘宏杰王伟力苗润吴世永 . 环形双锥罩聚能装药结构优化设计. 高压物理学报, 2018, 32(6): 065105-1-065105-9. doi: 10.11858/gywlxb.20180539
    [14] 章冠人 . 冲击波作用下铁电体的击穿. 高压物理学报, 1995, 9(1): 1-4 . doi: 10.11858/gywlxb.1995.01.001
    [15] 王同权张若棋汤文辉肖亚斌 . 强流质子束辐照形成的冲击波. 高压物理学报, 2003, 17(2): 117-121 . doi: 10.11858/gywlxb.2003.02.007
    [16] 章冠人 . 冲击波和分形集团的凝聚. 高压物理学报, 1997, 11(4): 241-244 . doi: 10.11858/gywlxb.1997.04.001
    [17] 谭华韩钧万贺红亮王小江 . 冲击波合成氮化硼的形貌及热稳定性. 高压物理学报, 1995, 9(1): 53-58 . doi: 10.11858/gywlxb.1995.01.009
    [18] 肖强强黄正祥顾晓辉 . 冲击波影响下的聚能射流侵彻扩孔方程. 高压物理学报, 2011, 25(4): 333-338 . doi: 10.11858/gywlxb.2011.04.008
    [19] 徐森刘大斌彭金华王建灵郭伟金朋刚贾宪政 . 药柱冲击波在有机玻璃中的衰减特性研究. 高压物理学报, 2010, 24(6): 431-437 . doi: 10.11858/gywlxb.2010.06.005
    [20] 陶俊林黎泽朋陈小伟 . 核爆聚变电站发电时的冲击波超压传播规律. 高压物理学报, 2012, 26(2): 221-226. doi: 10.11858/gywlxb.2012.02.016
  • 加载中
图(17)表(4)
计量
  • 文章访问数:  333
  • 阅读全文浏览量:  326
  • PDF下载量:  3
出版历程
  • 收稿日期:  2018-06-20
  • 录用日期:  2018-08-12
  • 网络出版日期:  2019-01-17
  • 刊出日期:  2019-02-01

串联战斗部不同介质组合的隔爆能力

    作者简介:刘宏杰(1993-),男,硕士研究生,主要从事弹药设计与目标毁伤评估研究. E-mail: 1300202650@qq.com
    通讯作者: 王伟力, 858769430@qq.com
  • 1. 海军航空大学岸防兵学院,山东 烟台 264001
  • 2. 海军工程大学兵器工程学院,湖北 武汉 430000
  • 3. 海军航空大学基础部,山东 烟台 264001
  • 4. 沈阳航空航天大学航空航天工程学部,辽宁 沈阳 110136

摘要: 为了有效提升串联切割战斗部隔爆结构衰减爆炸冲击波的性能,解决前级聚能装药结构与后级随进弹的匹配及隔爆问题,在前级切割器和后级随进弹之间加装隔爆结构,使用有限元分析软件ANSYS/LS-DYNA建立模型,进行不同组合结构隔爆性能的数值模拟,比较隔爆能力。模拟结果表明:前级装药起爆后,爆炸冲击波首先向后级随进弹头靠里区域汇聚,而不是向弹头尖端区域汇聚,因此可以适当减薄外层金属隔爆介质头部尖端区域;将外层金属由硬质钢改成铝时,后端壳体应力峰值的变化很小,故确定外层金属介质为铝;铝-聚脲的隔爆能力优于铝-泡沫铝结构,最终确定“软”隔爆介质为聚脲。通过调整铝和聚脲层的厚度,确定了最佳隔爆参数,能够满足实际应用。

English Abstract

  • 当前很多学者通过实验和数值仿真研究了聚脲衰减冲击波的特性[1-3]。曾必强等[4]通过数值模拟方法研究了不同初始速度的随进弹受前级装药爆炸冲击波影响产生的后抛、侧偏、翻转和力学响应规律,发现随进弹的后抛加速作用主要由前级装药爆炸产生的爆炸冲击波作用到弹体表面而产生,随进弹的侧偏现象则是作用在弹体的冲击波与冲击波后方的爆轰产物共同作用的结果。姜夕博等[5]、徐森等[6]研究了冲击波在有机玻璃中的衰减特性。侯海周等[7]的实验结果表明,药柱冲击波经酚醛隔板传播时,其峰值压力随着隔板厚度的增加呈指数衰减。然而,单层材料的隔爆能力十分有限,需要通过增加隔爆层厚度提升衰减冲击波能力,但是却占用了战斗部空间,增加了战斗部质量。近年来,许多学者对多层组合介质中冲击波的衰减特性开展了大量研究[8-10]。董永香等[11]通过分析不同介质层中和层界面的爆炸波形和幅值得出,在脆性材料中放置泡沫材料可以有效减小爆炸波应力峰值,降低自由表面的拉伸应力波幅值,软夹层在多层介质中能够起到明显的削波和增加加载时间的作用,中间夹层改变了爆炸波能量的分配。陈闯等[12]通过实验与数值模拟相结合的方法,将45钢、铝、有机玻璃进行不同组合,设计了顺序波阻抗梯度、逆序波阻抗梯度以及硬软硬3种多层隔爆结构,结果表明:逆序波阻抗梯度隔爆结构的输出冲击波压力最小,隔爆效果最佳。

    为了在有限的战斗部空间内尽可能提高隔爆效果并减轻隔爆体质量,本研究在后级随进弹头部壳体表面加装隔爆体,采用ANSYS/LS-DYNA有限元软件进行数值模拟,研究不同结构的隔爆特性,以期使后级随进弹能顺利随进并侵彻后续目标靶。

    • 高压下固体材料通常被认为是无黏性的可压缩流体,一般用介质Hugoniot关系曲线处理冲击波在两种不同介质上的反射和透射问题。定义介质的冲击阻抗为介质初始密度(固有)ρ0和冲击波速度DS的乘积,可近似用初始密度ρ0和声速C0的乘积表示。处理波从介质A传入介质B的最好方法是阻抗匹配法。碰撞界面的右行/左行冲击波满足连续方程和动量守恒。由界面连续条件可知:两介质界面上的粒子速度和压力相同。由一维冲击波理论可知,冲击波在两种界面处的反射波和透射波类型取决于两种介质的阻抗大小。当介质A的阻抗大于介质B时,反射波为稀疏波,反之反射波为冲击波[13-14]

      当冲击波传播至两层介质分界面处,可以利用介质的Hugoniot关系计算透射波和反射波,如图1所示,曲线1~曲线3分别表示不同介质的Hugoniot关系曲线,介质的阻抗越大,由原点出发与Hugoniot关系曲线上某点连线的斜率越大。图1反映了压力幅值为Pb的冲击波从介质A传入介质B的情形,分界面处的状态既在反射波Hugoniot曲线1′(曲线1′和曲线1呈镜像对称关系)上,也在介质B的右行透射波Hugoniot曲线上,即a点状态。同理,冲击波从介质A传入介质C时,状态从b点变为c点。

      图  1  反射波与透射波Hugoniot关系计算曲线

      Figure 1.  Hugoniot curves of reflection and transmission wave

      当固体处在高压条件下时,宋博等[15]提出用p=Au2分析冲击波在多层介质中的传播,其中A可用来表示材料的“相对软硬”程度,A值越大,材料越“硬”,反之,材料越“软”。假设冲击波初始状态为b(pb, ub),则曲线1′可以表示为

      则对于c

      b点状态为pb=A1ub2,联立(2)式,消去u,求得透射系数T2

      由(3)式可以得出冲击波传播n层介质后的透射系数

      其中

      由此可知对于相同层数(n)的相同介质来说,如果有某种顺序的排列结构,使得

      存在最大值,那么这种结构衰减冲击波的性能最好,即这种结构能最大限度地削弱透射波强度。

    • 第1节的分析仅考虑了不同介质的排列次序对衰减冲击波的影响,而未考虑每层介质材料本身对波的衰减及组合介质不同结构对衰减冲击波的影响。下面通过数值模拟方法首先确定最佳组合介质,以期最大限度地削弱透射波,同时兼顾质量和体积。

    • 为了便于建立计算模型,作以下假设:(1)空气、炸药、药型罩和装药壳体均为连续介质;(2)整个爆炸过程为绝热过程。串联战斗部及隔爆结构如图2所示。为了保证后级随进弹能顺利侵彻,隔爆体A选择金属介质,考虑到计算模型和爆炸载荷具有对称特性,为提高计算效率,建立1/4模型,施加对称边界条件。计算模型由炸药、空气、药型罩、壳体、后级随进弹和隔爆体组成,其中炸药、空气和药型罩定义为Euler网格,壳体、后级随进弹和隔爆体定义为Lagrange网格,隔爆体之间设置面面自动接触。选用多物质ALE(Arbitrary Lagrange-Euler)方法模拟分析前级爆轰场对后级随进弹的影响。起爆方式采用顶端环形起爆。

      图  2  串联战斗部及隔爆结构有限元模型

      Figure 2.  Finite element model of tandem warhead and explosive interruption structure

    • 炸药选用B炸药,采用MAT_HIGH-EXPLOSIVE-BURN模型和JWL状态方程[16-17]描述

      式中:p为爆轰压力;E为炸药比内能;v为相对比容;ABR1R2ω为参数,取值列于表1表1中:ρ为密度,D为B炸药爆速,pCJ为CJ面压力,V0为初始相对比容。

      ρ/(g·cm-3) D/(m·s-1) pCJ/GPa A/GPa B/GPa R1 R2 ω E0/GPa V0
      1.713 7 500 28.6 524.2 7.678 4.2 1.1 0.34 8.499 1.0

      表 1  B炸药材料参数

      Table 1.  Material parameters of composition B

      药型罩材料选用紫铜,采用MAT_Steinberg材料模型和Grüneisen状态方程[18]描述

      式中:E为单位体积内能;ρ0为材料初始密度;Cusup的交点,其中us为冲击波速度,up为粒子速度;S1S2S3us-up曲线的斜率;γ0为Grüneisen系数;a为对γ0的修正;μ=(1/V)-1,V为当前空气的相对体积。

      壳体选用30CrMnSiNi2A钢,用MAT_Johnson-Cook材料模型和Grüneisen状态方程描述;铝选用MAT_Johnson-Cook材料模型和Grüneisen状态方程描述,参数见表2,其中AJCBJCnm为Johnson-Cook材料模型参数,Tm为材料熔化温度,T0为常压下的熔化温度,a为对Grüneisen系数γ0的修正。空气采用MAT_Null材料模型和Linear-Polynomial状态方程描述。

      Material AJC/MPa BJC/MPa n C m Tm/K T0/K S1 γ0 a
      Steel 1 280 420 0.30 0.030 1.00 1 793 294 1.490 2.17 0.46
      Al 265 426 0.34 0.015 1.00 775 294 1.345 2.13 0.10

      表 2  30CrMnSiNi2A钢和金属铝计算参数

      Table 2.  Material performance parameters of 30CrMnSiNi2A steel and aluminum

      聚脲选用MAT_PIECEWISE_LINEAR_PLASTICITY材料模型[19]描述。该模型常用来模拟弹塑性材料,屈服准则由断裂应变定义,应力-应变曲线由屈服强度、剪切模量定义或者通过8组塑性应变和应力的数组定义,应变率的影响采用Cowper-Symbols模型,弹性体材料各参数为[20-21]: ρ=1 020 kg/m3,杨氏模量E=230 MPa,泊松比ν=0.4,初始屈服强度σ0=14 MPa。

    • 用B1、B2、B3和B4分别表示隔爆体A/B、B/C、C/D、D/后级弹壳体的界面,计算中两个界面处均采用面面自动接触条件。表3列出了5种工况。

      Condition Material
      A B C D
      Steel Polyurea Aluminum foam Aluminum
      Aluminum Polyurea Aluminum foam Aluminum
      Aluminum Aluminum foam Aluminum foam Aluminum
      Aluminum Aluminum foam Aluminum foam Aluminum foam
      Aluminum Polyurea Polyurea Polyurea

      表 3  不同工况及简便书写方式

      Table 3.  Different conditions and indications

    • 工程实践表明:具有软夹层的硬软硬多层介质具有较明显的抗冲击和抗爆炸能力,其中硬层抗击近距离的爆炸波作用,发挥其刚度与强度效应,而软层由多孔泡沫材料组成,波阻抗低,强度和刚度低,吸波能力强。如果外加爆炸波的持续时间足够长,软夹层不再发挥作用,软夹层下面的硬层将受到较大的冲击作用。

      这里假定几种材料的Hugoniot曲线分别为p=A1u(硬质钢),p=A2u(铝),p=A3u(泡沫铝),p=A4u (聚脲),并且有A1>A2>A3>A4

    • (1)工况I

      对“硬质钢-聚脲-泡沫铝-铝”结构,有

      装药起爆后,冲击波首先撞击复合隔爆结构中的硬质钢,该碰撞过程为非弹性碰撞。在该冲击速度下钢板发生弹塑性变形,历时很短,通常钢板发生局部的弹塑性变形,而其他区域受碰撞挤压作用影响较小。当冲击波传播至B1界面时,发生反射和透射,聚脲作为一种弹性涂覆体,其冲击阻抗小于硬质钢,界面B1处反射稀疏波,然而聚脲层的厚度仅为1 mm,衰减冲击波的能力有限。冲击波从介质B传入介质C时,在界面B2处反射冲击波,接下来讨论中间泡沫铝层的塑性变形和卸载冲击波强度的变化规律。泡沫铝在冲击载荷作用下产生递增硬化的变形即破坏波(压实波),其非线性应力-应变曲线可简化成3个阶段:弹性段、屈服塑性段、压实段,其力学性能除了与泡沫的泡孔结构有关外,还与基体材料性质、应变率和密度等有关,而且对静水压有很强的依赖性。应力波在泡沫铝中经历以下两个步骤:一是泡沫铝压实过程中的冲击波传播;二是泡沫铝压实后,应力波在包括外层硬质钢、聚脲涂覆层、压实泡沫铝以及后端铝介质的多次反射和透射。冲击波从介质C传入介质D时,在界面B3处反射冲击波,经过前面3层介质的衰减作用后,铝介质的冲击压力较低,铝介质作为一种软金属材料,在冲击载荷作用下易发生变形和破坏,因而能最大限度地吸收能量,最大程度地减小后级装药的压力峰值。

      为了检验工况I衰减冲击波的能力,查看随进弹壳体最大应力集中的区域以及装药压力峰值的区域和时刻。t=264 μs时,装药内部靠近轴线区域压力(见图3(a))达到最大值,为143.6 MPa;t=189 μs时,外侧壳体应力集中区域的应力峰值如图3(b)所示,可见该区域未发生塑性变形;应力峰值附近某单元应力-时间变化关系如图3(c)所示。

      图  3  工况I中装药压力及壳体应力分析

      Figure 3.  Loading pressure and shell stress of working condition I

      图4可知:t=30 μs时,外层钢介质头部尖端靠里区域首先发生弹性变形;随着冲击波在介质中继续传播,t=38 μs时,钢介质弹性变形区域向后级随进弹方向移动;t=45 μs时,应力集中区域为钢介质最薄处;t=51 μs时,钢介质最薄处应力峰值继续增加;t=60 μs时,应力峰值达到最大值,为1 130 MPa。前级装药采用中空结构,装药起爆后爆轰波向后运动并汇聚的过程中,首先作用到钢介质尖端靠里区域,考虑到节约质量成本,钢介质头部尖端区域可适当减薄。

      图  4  典型时刻外层钢介质的应力云图

      Figure 4.  Stress nephogram of steel at typical times

      (2)工况II

      与工况I相比,为了控制质量,外侧介质A由密度较大的硬质钢变为密度较小的铝。透射系数影响因子为

      比较Sn1Sn2的大小,得到Sn1>Sn2,可知工况I所示结构中衰减冲击波压力的性能优于工况II所示结构。

      装药起爆后,冲击波首先撞击复合隔爆结构中的金属铝,该碰撞过程为非弹性碰撞。金属铝的屈服强度小,在高速冲击载荷下铝板发生破坏,部分区域发生弹塑性变形。冲击波在聚脲层、泡沫铝和后端铝介质传播衰减分析方法与工况I类似。

      图5可知:t=120 μs时,装药尖端压力达到最大值,为421.2 MPa;相同时刻,外侧壳体应力出现最大值1 062.0 MPa,见图5(b);应力峰值附近某单元的应力-时间曲线如图5(c)所示,后级随进弹壳体未进入塑性变形段。硬软硬结构存在不利的一面:随着作用时间的增加,对于多层介质而言,在波动效应后期结构振动效应增强。

      图  5  工况II装药压力及壳体应力分析

      Figure 5.  Loading pressure and shell stress of working condition II

      (3)工况III

      工况III所示隔爆结构为铝-泡沫铝-铝硬软硬结构,与工况II相比,将介质B换为泡沫铝。含有低阻抗软夹层的多层介质能够减小应力幅值,增加波形脉宽(即增加作用时间),从而减小由于应力波峰超过材料强度而引起局部破坏的可能性。

      图6可知:t=117 μs时,装药尖端压力达到最大值363.5 MPa,与中间有聚脲层时对应的压力峰值时刻基本一致;t=174 μs时,外侧壳体应力峰值为1 062.0 MPa,峰值区域如图6(b)所示;应力峰值附近单元的应力-时间曲线如图6(c)所示,后级随进弹壳体未进入塑性变形段。

      图  6  工况III装药压力及壳体应力

      Figure 6.  Loading pressure and shell stress of working condition III

    • (1)工况IV

      隔爆介质为两层,外层为金属铝,里层为泡沫铝,透射系数影响因子为

      对比工况IV(见图7)和工况III:将厚端铝介质改为泡沫铝后,增加了泡沫铝的厚度,在爆炸冲击波作用下,泡沫铝内衬层能够吸收更多的能量,随进弹装药压力峰值比硬软硬结构中装药应力峰值出现得晚些;冲击波继续传播过程中,泡沫铝衰减冲击波的能力不再增加,最终泡沫铝被压实,密度增加,弹性模量也增加,紧挨软夹层的金属层受到较大的冲击作用,继续衰减冲击波。

      图  7  工况IV装药压力及壳体应力

      Figure 7.  Loading pressure and shell stress of working condition IV

      (2)工况V

      隔爆介质分层排列顺序为金属铝和弹性涂覆体聚脲,其透射系数影响因子为

      计算可得Sn5>Sn4

      考虑将隔爆结构简化为金属铝-聚脲结构(即工况V),对比工况IV(即将泡沫铝改为弹性介质聚脲),则弹性体厚度明显增加。工况V中装药压力及壳体应力见图8t=156 μs时,装药尖端压力达到最大值272.0 MPa,与工况IV中出现压力峰值时刻一致;t=192 μs时,外侧壳体应力峰值为1 050.0 MPa;应力峰值附近单元的应力-时间曲线显示,后级随进弹壳体未进入塑性变形段。聚脲作为一种喷涂材料,受制作工艺限制,喷涂有效厚度有限,该工况的聚脲层厚度为4 cm,远大于聚脲层的最大有效厚度。

      图  8  工况V装药压力及壳体应力

      Figure 8.  Loading pressure and shell stress of working condition V

      图9可知:采用工况II和工况III所示隔爆结构时,后级随进弹出现加速度峰值的时间最早,且峰值最大;工况I所示隔爆结构中后级弹的加速度峰值最小,工况IV和工况I出现加速度峰值的时间基本一致;工况V所示结构延缓了冲击波作用到后级随进弹的时间,出现加速度峰值的时间最晚。由图10可知,采用上述5种隔爆结构时,后级随进弹速度增加到一定值时不再增加。

      图  9  后级随进弹加速度-时间变化曲线

      Figure 9.  Acceleration vs. time for the following stage’s incoming projectile

      图  10  后级随进弹速度-时间变化曲线

      Figure 10.  Velocity vs. time for the following stage’s incoming projectile

    • 采用工况I的隔爆结构时,后端装药的压力峰值出现得最晚,且压力峰值最小,仅为143.6 MPa;受质量成本的限制,工况II中介质A由密度较大的硬质钢换成密度较小的金属铝,金属铝的屈服强度较低,爆炸冲击载荷下容易发生破坏,部分区域发生弹塑性变形,装药压力峰值区域为装药尖端区域,远早于A介质为硬质钢的情形,且峰值压力为421.2 MPa,远大于工况I中对应的峰值压力,但后级壳体应力峰值与工况I的应力峰值相近,均未进入塑性变形区;工况III与工况II相比,介质B由弹性涂覆体聚脲改为泡沫铝,即隔爆结构转变为“铝-泡沫铝-铝”硬软硬结构,由于介质B的厚度仅为1 mm,对冲击波反射和透射的影响不大,因此硬软硬结构出现装药压力峰值的时刻与工况II基本一致,装药峰值压力小于工况II的情形;工况IV和工况V都降低了隔爆层数,分别为铝-泡沫铝结构和铝-聚脲结构,对比透射系数影响因子,前者小于后者,因此得出铝-聚脲结构隔爆效果优于铝-泡沫铝结构,但工况V中聚脲层的厚度过大,无法实现实际应用。接下来将改进隔爆结构,使之既能充分提升隔爆性能,节约体积和质量成本,又能实现实际应用。

    • 由第3节分析可知,装药起爆后,冲击波首先作用到最外层隔爆体头部靠里区域。为了控制总质量,适当减薄金属层尖端区域的厚度,改进后隔爆结构如图11所示,保持后级弹尖端与前级聚能装药结构后端的距离不变,即图中d表4列出了4种改进结构。

      图  11  串联战斗部及改进隔爆结构有限元模型

      Figure 11.  Finite element model of tandem warhead and improved explosive interruption structure

      Condition Material
      A B C
      Ⅰ′ Steel Aluminum Polyurea
      Ⅱ′ Steel Aluminum Aluminum foam
      Ⅲ′ Aluminum Aluminum Polyurea
      Ⅳ′ Aluminum Aluminum Polyurea(2 cm)

      表 4  改进结构

      Table 4.  Improved structures

    • (1)工况I′

      工况I′下透射系数影响因子为

      图12可知:t=163 μs时,装药尖端压力达到最大值114.8 MPa;t=243 μs时,外侧壳体应力峰值为1 188.0 MPa;应力峰值附近单元的应力-时间曲线显示,后级随进弹壳体尖端区域进入塑性变形区。

      图  12  工况I′装药压力及壳体应力

      Figure 12.  Loading pressure and shell stress of working condition I′

      (2)工况II′

      工况II′下,透射系数影响因子为

      比较可得,工况I′对应的透射系数大于工况II′,说明工况I′所示隔爆结构的衰减冲击波能力大于工况II′。由图13可知:t=246 μs时,装药尖端压力达到最大值125.7 MPa;t=297 μs时,外侧壳体应力峰值为1 209.0 MPa;应力峰值附近单元的应力-时间曲线显示,后级随进弹壳体尖端区域进入塑性变形区。

      图  13  工况II′装药压力及壳体应力分析

      Figure 13.  Loading pressure and shell stress of working condition II′

    • (1)工况III′

      将最外层密度较大的硬质钢改为密度较小的金属铝,采用两层隔爆层,即金属铝-聚脲结构。由图14可知:t=159 μs时,装药内部尖端区域压力达到最大值170.8 MPa;t=195 μs时,外侧壳体应力峰值为1 005.0 MPa;应力峰值附近单元的应力-时间曲线显示,后级随进弹壳体未进入塑性变形区。

      图  14  工况III′装药压力及壳体应力分析

      Figure 14.  Loading pressure and shell stress of working condition III′

      (2)工况IV′

      考虑到聚脲喷涂的厚度有限,将3 cm聚脲层减薄至2 cm,研究其隔爆性能。由图15可知:t=144 μs时,装药内部尖端区域压力达到最大值296.3 MPa;t=177 μs时,外侧壳体应力峰值为1 035.0 MPa;应力峰值附近单元的应力-时间曲线显示,后级随进弹壳体未进入塑性变形区。

      图  15  工况IV′装药压力及壳体应力分析

      Figure 15.  Loading pressure and shell stress of working condition IV′

      图16图17可知:采用工况IV′所示隔爆结构时,后级随进弹出现加速度峰值的时间最早,且峰值最大;工况I′所示隔爆结构中,后级随进弹的加速度峰值最小;工况IV′和工况III′出现加速度峰值的时间基本一致;工况II′所示结构延缓了冲击波作用到后级随进弹的时间,出现加速度峰值的时间最晚。由图17可知,采用以上所述4种隔爆结构时,后级随进弹速度增加到一定值时将不再增加。

      图  16  后级随进弹加速度-时间变化曲线

      Figure 16.  Acceleration vs. time for the following stage’s incoming projectile

      图  17  后级随进弹加速度-时间曲线

      Figure 17.  Velocity vs. time for the following stage’s incoming projectile

    • 总结改进后的工况:采用工况I′对应的隔爆结构时,后端装药压力峰值最小,仅为114.8 MPa;受质量成本的限制,工况II′中将介质C换成泡沫铝,延缓了冲击波作用到达后级随进弹的时间。但是在工况I′和工况II′对应的隔爆结构中,由于外层硬质钢挤压变形,使得后级随进战斗部壳体头部发生变形,进入塑性变形区;工况IV′与工况III′相比,介质C对应弹性涂覆体聚脲的厚度减薄,装药压力较工况III′所示结构有所提高,但增幅不大,后端壳体应力峰值的变化不大。

    • 采用数值模拟方法研究不同介质组合衰减冲击波的能力,通过对比后级弹壳体应力峰值及装药压力峰值,考察不同介质组合的隔爆能力,得出以下结论:(1)前级中空结构装药起爆后,爆炸冲击波首先向后级随进弹头靠里区域汇聚,而不是向弹头尖端区域汇聚;(2)在“硬软硬”隔爆结构中,将外层介质A由密度较大的硬质钢改为金属铝后,后级弹壳体应力峰值的变化很小,故确定隔爆结构中外层金属介质为铝,从而大大节约质量成本;(3)“硬软”两层组合介质隔爆结构中,铝-聚脲的隔爆能力优于铝-泡沫铝结构,最终确定“软”隔爆介质为聚脲,但聚脲的喷涂厚度有限,需要进一步减薄聚脲层的厚度;(4)改进结构中,减薄了外层金属介质层尖端的厚度,若保持隔爆层厚度不变,通过增加金属铝的厚度降低聚脲层的厚度,最终确定聚脲层的厚度为2 cm,从而既节约质量和体积成本,又具有优良的隔爆性能,能够充分削弱作用在后级随进弹上的冲击波强度,使后级随进弹的加速度和速度较小。

参考文献 (21)

目录

    /

    返回文章
    返回