内爆加载下金属柱壳的冻结回收方法

罗渝松 李伟兵 陈志闯 王晓鸣 李文彬

引用本文:
Citation:

内爆加载下金属柱壳的冻结回收方法

    作者简介: 罗渝松(1995- ),男,博士研究生,njlgdxlys@163.com;
    通讯作者: 李伟兵, njustlwb@163.com
  • 中图分类号: O381;TJ410

A freezing recovery method for metallic cylinder shells under internal explosive loading

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

  • 摘要: 针对内爆炸载荷下膨胀态金属柱壳的回收问题,设计了冻结回收试验方法,实现了起爆后不同时刻金属柱壳的冻结回收。基于一体式壳体提出了3种改进结构,并分别对4种柱壳结构在内爆载荷下的膨胀断裂过程进行数值模拟。结果表明,两段式结构最有利于减小非起爆端对预期回收中间段壳体的影响。根据选定的最优壳体结构和金属柱壳在起爆后不同时刻的膨胀外形特征,设计与之匹配的冻结回收装置并进行冻结回收试验。试验结果表明,设计的冻结回收试验方法可以实现膨胀态金属柱壳的回收,回收壳体的轴向和径向尺寸与设计理想值符合较好,整体误差可控制在10%以内。
  • 图 1  4种壳体的结构方案

    Figure 1.  Structures of four shells

    图 2  一体式壳体的仿真计算模型(1~21为观测点)

    Figure 2.  The simulation calculation model for the integrated shell (1−21 are observed points)

    图 3  起爆后30 μs时刻4种壳体的径向位移

    Figure 3.  Radial displacements of four shells at 30 μs after detonation

    图 4  壳体断裂分区

    Figure 4.  Division of fracture area of the shell

    图 5  一体式壳体的膨胀变形响应压力云图

    Figure 5.  Variation of internal pressure of integrated shell with expansion and deformation process

    图 6  匹配不同时刻壳体膨胀外形设计的冻结回收装置结构尺寸

    Figure 6.  Structural dimensions of designed freezing recovery devices matching deformation responses of shells at different times

    图 7  3个回收时刻壳体径向变形响应

    Figure 7.  Radial deformation response of the shell at three recovery times

    图 8  3个时刻下的壳体回收半径

    Figure 8.  Recovery radii of the shell at three recovery times

    图 9  金属柱壳装药结构

    Figure 9.  Charge structures of the metal cylindrical shell

    图 10  冻结回收试验总体布局

    Figure 10.  Overall layout of the freeze recovery test

    图 11  23 µs时刻的壳体冻结回收试验结果

    Figure 11.  Shell obtained from freezing recovery test at 23 μs

    图 12  23 µs时刻中间段壳体的回收半径

    Figure 12.  Recovery radii of the middle section shell at 23 µs

  • [1] TAYLOR G I. Fragmentation of tubular bombs [M] // BATCHELOR O K. The Scientific Papers of Sir Geoffrey Ingram Taylor. Cambridge: Cambridge University Press, 1963: 387−390.
    [2] ZHU J J, LI W B, WANG X M, et al. Effect of tempering temperature on expansion and fracture mechanism of 40CrMnSiB steel cylinder shell [J]. International Journal of Impact Engineering, 2017, 107: 38–46. DOI: 10.1016/j.ijimpeng.2017.05.007.
    [3] 陈醇. 三代炸药爆炸加载下战斗部壳体形变规律及破片性能研究[D]. 南京: 南京理工大学, 2015: 14−22.
    [4] 朱建军, 李伟兵, 王晓鸣, 等. 回火温度对50SiMnVB钢壳体形成破片性能的影响 [J]. 兵工学报, 2015, 36(11): 2080–2086. DOI: 10.3969/j.issn.1000-1093.2015.11.009.
    ZHU J J, LI W B, WANG X M, et al. Effect of tempering temperature on the forming properties of fragments of 50SiMnVB steel shell [J]. Acta Armamentarii, 2015, 36(11): 2080–2086. DOI: 10.3969/j.issn.1000-1093.2015.11.009.
    [5] 朱如意, 韩保红, 赫万恒. 杀爆榴弹战斗部静态爆炸效能仿真研究 [J]. 计算机仿真, 2015, 32(11): 18–21. DOI: 10.3969/j.issn.1006-9348.2015.11.005.
    ZHU R Y, HAN B H, HE W H. Simulation study on static explosive efficiency of blast fragmentation warhead [J]. Computer Simulation, 2015, 32(11): 18–21. DOI: 10.3969/j.issn.1006-9348.2015.11.005.
    [6] SINGH M, SUNEJA H R, BOLA M S, et al. Dynamic tensile deformation and fracture of metal cylinders at high strain rates [J]. International Journal of Impact Engineering, 2002, 27(9): 939–954. DOI: 10.1016/S0734-743X(02)00002-7.
    [7] WANG X Y, WANG S S, MA F. Experimental study on the expansion of metal cylinders by detonation [J]. International Journal of Impact Engineering, 2018, 114: 147–152. DOI: 10.1016/j.ijimpeng.2017.12.017.
    [8] GOLD V M, BAKER E L. A model for fracture of explosively driven metal shells [J]. Engineering Fracture Mechanics, 2008, 75(2): 275–289. DOI: 10.1016/j.engfracmech.2007.02.025.
    [9] HUANG G Y, LI W, FENG S S. Axial distribution of fragment velocities from cylindrical casing under explosive loading [J]. International Journal of Impact Engineering, 2015, 76: 20–27. DOI: 10.1016/j.ijimpeng.2014.08.007.
    [10] GUO Z W, HUANG G Y, LIU C M, et al. Velocity axial distribution of fragments from non-cylindrical symmetry explosive-filled casing [J]. International Journal of Impact Engineering, 2018, 118: 1–10. DOI: 10.1016/j.ijimpeng.2018.03.011.
    [11] SUN G J, WANG X Y, GAO W L, et al. Expansion fracture behavior of metallic cylindrical shell caused by explosive detonation [J]. IOP Conference Series: Earth and Environmental Science, 2019, 267(4): 042124.
    [12] 禹富有, 董新龙, 俞鑫炉, 等. 不同填塞装药下金属柱壳断裂特性的实验研究 [J]. 兵工学报, 2019, 40(7): 1418–1424. DOI: 10.3969/j.issn.1000-1093.2019.07.011.
    YU F Y, DONG X L, YU X L, et al. Fracture characteristics of metal cylinder shells with different charges [J]. Acta Armamentarii, 2019, 40(7): 1418–1424. DOI: 10.3969/j.issn.1000-1093.2019.07.011.
    [13] 奥尔连科. 爆炸物理学: 下册 [M]. 孙承纬, 译. 北京: 科学出版社, 2011: 807−822.
    [14] 陈志闯, 李伟兵, 朱建军, 等. 40CrMnSiB钢圆柱壳体膨胀断裂中间状态回收试验研究 [J]. 兵工学报, 2018, 39(11): 2137–2144. DOI: 10.3969/j.issn.1000-1093.2018.11.007.
    CHEN Z C, LI W B, ZHU J J, et al. Recovery experiment study of cylindrical 40CrMnSiB steel shell in intermediate phase of expanding fracture processes [J]. Acta Armamentarii, 2018, 39(11): 2137–2144. DOI: 10.3969/j.issn.1000-1093.2018.11.007.
    [15] TAYLOR G I. Analysis of the explosion of a long cylindrical bomb detonated at one end [M] // BATCHELOR G K. The Scientific Papers of Sir Geoffrey Ingram Taylor: Ⅲ. Cambridge: Cambridge University Press, 1963: 277−286.
    [16] 任国武,郭昭亮,汤铁钢, 等. 高应变率加载下金属柱壳断裂的实验研究 [J]. 兵工学报, 2016, 37(1): 77–82. DOI: 10.3969/j.issn.1000-1093.2016.01.012.
    REN G W, GUO S L, TANG T G, et al. Experimental research on fracture of metal case under loading at high strain rate [J]. Acta Armamentarii, 2016, 37(1): 77–82. DOI: 10.3969/j.issn.1000-1093.2016.01.012.
    [17] 陈志闯. 基于冻结回收的金属壳体膨胀断裂机制研究[D]. 南京: 南京理工大学, 2019: 27−28.
    [18] 顾文彬, 胡亚峰, 徐浩铭, 等. 复合结构防爆罐抗爆特性的数值模拟 [J]. 含能材料, 2014, 22(3): 325–331. DOI: 10.3969/j.issn.1006-9941.2014.03.010.
    GU W B, HU Y F, XU H M, et al. Numerical simulation of blast resistant characteristics for the composite structure anti-explosion container [J]. Chinese Journal of Energetic Materials, 2014, 22(3): 325–331. DOI: 10.3969/j.issn.1006-9941.2014.03.010.
    [19] 汤铁钢, 李庆忠, 孙学林, 等. 45钢柱壳膨胀断裂的应变率效应 [J]. 爆炸与冲击, 2006, 26(2): 129–133. DOI: 10.11883/1001-1455(2006)02-0129-05.
    TANG T G, LI Q Z, SUN X L, et al. Strain-rate effects of expanding fracture of 45 steel cylinder shells driven by detonation [J]. Explosion and Shock Waves, 2006, 26(2): 129–133. DOI: 10.11883/1001-1455(2006)02-0129-05.
  • [1] 汤铁钢李庆忠孙学林孙占峰金山谷岩 . 45钢柱壳膨胀断裂的应变率效应. 爆炸与冲击, 2006, 26(2): 129-133. doi: 10.11883/1001-1455(2006)02-0129-05
    [2] 金山汤铁钢孙学林李庆忠 . 不同热处理条件下45钢柱壳的动态性能. 爆炸与冲击, 2006, 26(5): 423-428. doi: 10.11883/1001-1455(2006)05-0423-06
    [3] 刘龙飞周强 . 表面粗糙度对6061铝合金薄壁管冲击膨胀断裂性能的影响. 爆炸与冲击, 2018, 38(4): 749-758. doi: 10.11883/bzycj-2016-0389
    [4] 叶想平李英雷李英华 . 基于SHPB装置的膨胀圆柱管实验技术. 爆炸与冲击, 2012, 32(5): 528-534. doi: 10.11883/1001-1455(2012)05-0528-07
    [5] 任国武郭昭亮张世文汤铁钢金山胡海波 . 金属柱壳膨胀断裂的实验与数值模拟. 爆炸与冲击, 2015, 35(6): 895-900. doi: 10.11883/1001-1455(2015)06-0895-06
    [6] 胡海波汤铁钢胡八一王德生韩长生 . 金属柱壳在爆炸加载断裂中的单旋现象. 爆炸与冲击, 2004, 24(2): 97-107.
    [7] 高重阳施惠基姚振汉白春华 . 薄壁柱壳在内部爆炸载荷下膨胀断裂的研究. 爆炸与冲击, 2000, 20(2): 160-167.
    [8] 刘明涛汤铁钢胡海波2李庆忠胡秀章李永池 . 不同起爆方式下炸药驱动柱壳膨胀断裂的数值模拟. 爆炸与冲击, 2014, 34(4): 415-420. doi: 10.11883/1001-1455(2014)04-0415-06
    [9] 汤铁钢谷岩李庆忠华劲松孙学林 . 爆轰加载下金属柱壳膨胀破裂过程研究. 爆炸与冲击, 2003, 23(6): 529-533.
    [10] 任国武温上捷张茹郭昭亮汤铁钢 . 约束层对金属柱壳膨胀变形影响的数值模拟. 爆炸与冲击, 2017, 37(6): 946-951. doi: 10.11883/1001-1455(2017)06-0946-06
    [11] 郭昭亮范诚刘明涛任国武汤铁钢刘仓理 . 爆炸与电磁加载下无氧铜环、柱壳的断裂模式转变. 爆炸与冲击, 2017, 37(6): 1072-1079. doi: 10.11883/1001-1455(2017)06-1072-08
    [12] 郭昭亮任国武张世文汤铁钢刘仓理 . 膨胀柱壳恒定应变率的本构关系. 爆炸与冲击, 2016, 36(5): 583-589. doi: 10.11883/1001-1455(2016)05-0583-07
    [13] 李忠盛吴护林陈韵如易同斌李立潘复生 . 内爆炸载荷作用下7A55铝合金的动态性能及断裂行为. 爆炸与冲击, 2012, 32(2): 190-195. doi: 10.11883/1001-1455(2012)02-0190-06
    [14] 李永池董杰胡秀章范书群 . 内聚爆炸载荷下热塑性中厚球壳的变形和层裂. 爆炸与冲击, 2005, 25(4): 319-324. doi: 10.11883/1001-1455(2005)04-0319-06
    [15] 陈晓坤李海涛王秋红金永飞邓军张嬿妮 . 瓦斯爆炸载荷作用下矿用柱壳救生舱抗爆性分析及结构优化. 爆炸与冲击, 2018, 38(1): 124-132. doi: 10.11883/bzycj-2016-0248
    [16] 金山刘欣袁帅华劲松汤铁钢 . 爆轰加载下柱壳断裂时刻微小差异测量方法. 爆炸与冲击, 2015, 35(1): 130-134. doi: 10.11883/1001-1455(2015)01-0130-05
    [17] 李雪梅汪小松王鹏来卢敏贾路峰 . 内爆和柱壳条件下无氧铜的层裂特性. 爆炸与冲击, 2009, 29(2): 162-166. doi: 10.11883/1001-1455(2009)02-0162-05
    [18] 颜怡霞陈裕泽陈刚李思忠黄西成 . 激光辐照下充内压柱壳动态爆裂的数值模拟. 爆炸与冲击, 2004, 24(6): 487-492.
    [19] 张志彪黄风雷 . 内部爆炸加载下柱壳的环周分裂数. 爆炸与冲击, 2015, 35(5): 763-767. doi: 10.11883/1001-1455(2015)05-0763-05
    [20] 洪江波侯海量朱锡阚于龙 . 爆炸处理消除大型柱壳结构环焊缝残余应力实验研究. 爆炸与冲击, 2008, 28(3): 213-219. doi: 10.11883/1001-1455(2008)03-0213-07
  • 加载中
图(12)
计量
  • 文章访问数:  1556
  • HTML全文浏览量:  695
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-24
  • 录用日期:  2020-06-29
  • 网络出版日期:  2020-09-29
  • 刊出日期:  2020-10-05

内爆加载下金属柱壳的冻结回收方法

    作者简介:罗渝松(1995- ),男,博士研究生,njlgdxlys@163.com
    通讯作者: 李伟兵, njustlwb@163.com
  • 1. 南京理工大学智能弹药技术国防重点学科实验室,江苏 南京 210094
  • 2. 南京航空航天大学机电学院,江苏 南京 210016

摘要: 针对内爆炸载荷下膨胀态金属柱壳的回收问题,设计了冻结回收试验方法,实现了起爆后不同时刻金属柱壳的冻结回收。基于一体式壳体提出了3种改进结构,并分别对4种柱壳结构在内爆载荷下的膨胀断裂过程进行数值模拟。结果表明,两段式结构最有利于减小非起爆端对预期回收中间段壳体的影响。根据选定的最优壳体结构和金属柱壳在起爆后不同时刻的膨胀外形特征,设计与之匹配的冻结回收装置并进行冻结回收试验。试验结果表明,设计的冻结回收试验方法可以实现膨胀态金属柱壳的回收,回收壳体的轴向和径向尺寸与设计理想值符合较好,整体误差可控制在10%以内。

English Abstract

  • 在现有弹药战斗部的结构设计中,多数战斗部均采用柱壳装药结构,其破片质量分布、初速、飞散角等毁伤威力参数与金属柱壳的膨胀断裂过程直接相关[1-3]。国内外学者对内爆加载下金属柱壳的膨胀断裂过程采用了多种不同的方法进行研究,朱建军等[4]、朱如意等[5]数值模拟了內爆加载下金属圆柱壳的膨胀断裂过程,得到不同装药和不同结构下的壳体动态破碎规律。Singh等[6]、Wang等[7]利用高速分幅相机和光子多普勒测速系统研究了内爆载荷作用下圆柱形壳体的整个加速过程,得到壳体膨胀断裂参数与材料和结构之间的联系。Gold等[8]、Huang等[9]、Guo等[10]使用脉冲X光照相技术对不同结构柱壳装药下的壳体膨胀形态和破片飞散特征进行研究,建立了标准柱壳和非标准柱壳的破片初速沿轴向的分布修正公式。Sun等[11]运用高速摄影和扫描电镜分析技术分别对爆炸驱动下分别处于膨胀断裂状态和完全破碎状态的金属柱壳进行图像观测,分析表明金属柱壳断裂模式与炸药驱动能力和金属性能有关。禹富有等[12]进行了TA2钛合金柱壳在不同厚度药柱填塞下的爆炸实验,采用金相分析手段对回收破片进行破坏分析,发现柱壳的膨胀断裂模式与载荷脉宽和柱壳壁厚的比值相关。对于爆炸加载下壳体膨胀断裂过程的研究,已有的文献主要采取高速摄影和X光照相等手段来获取壳体的膨胀断裂图像,并对完全破碎后的壳体破片进行回收,辅以金相和扫描电镜等方式来完成分析[13-14]。但高速摄影和X光照相等拍摄手段获取的图像只能粗略地展示壳体外表面的膨胀断裂特征,对于内表面的情况和局部细节的体现则不足,而回收的破片则是完整壳体发生断裂失效后的结果,不能很好地体现壳体由初始膨胀变形发展至最终完全破碎的中间过程。为解决上述问题,需要在内爆炸载荷环境下回收处于膨胀断裂中间状态的金属壳体。

    本文中,针对内爆加下膨胀半径具有明显轴向效应的大长径比(L/d≥3)金属柱壳[15],利用AUTODYN软件数值模拟4种结构壳体的膨胀变形过程,统计、对比4种结构壳体自由膨胀的径向位移数值确定最优壳体结构。基于最优结构壳体在起爆后不同时刻的膨胀状态,匹配其外形设计冻结回收装置并进行回收过程的数值模拟和试验验证,对回收壳体的膨胀半径进行定量分析,以期为实现内爆加载下膨胀态金属柱壳的冻结回收提供参考。

    • 金属柱壳在爆轰产物的作用下首先是起爆端发生膨胀断裂,而非起爆端壳体的继续膨胀运动将不可避免地影响中间段壳体的回收效果[16-17]。为减小非起爆端对中间回收段壳体断裂形貌的影响,本节在一体式壳体结构的基础上进行改进,设计了两段式、刻槽式和单段式等3种柱壳结构,利用AUTODYN数值模拟软件对比研究4种壳体在内部爆炸驱动作用下的动态变形过程,分析起爆后各时刻壳体膨胀断裂形态的差异,并考虑到实际加工因素,确定用于进行冻结回收的最佳壳体结构。

    • 壳体的仿真计算方案是根据所设计的壳体结构而定的,除一体式结构外还包括两段式、刻槽式、单段式这3种改进结构,壳体改进遵循的原则为:结构改进后的壳体既要在预期回收段上与一体式结构壳体的膨胀外形相同,又要保证预期回收段以外的非起爆端部分对回收的影响要尽可能小。对于柱壳结构而言,减小非起爆端对预期回收段影响的具体做法是将壳体改动的位置设在变形和未变形的交界面处,使非起爆端壳体在內爆驱动作用下能与预期回收段壳体实现分离。根据壳体改进原则,在一体式结构壳体基础上设计了其余3种改进结构壳体,如图1(b)(d)所示。

      图  1  4种壳体的结构方案

      Figure 1.  Structures of four shells

      图1中一体式结构为完整圆柱壳体,其外表不进行任何削弱强度的机械加工,总长度为160 mm,其余3种结构则分别在距离起爆端120 mm位置处进行改进。两段式结构是在一体式结构基础上沿120 mm截面切断后粘接为一体,其中0~120 mm区间段包含了预期回收的壳体,120~160 mm区间段位于非起爆端,可降低稀疏波的提前传入对壳体膨胀产生的影响。刻槽式结构使用了刻槽工艺对壳体的机械强度进行部分削弱,使壳体在內爆驱动作用下能沿预设位置处断裂,实现前后两段壳体的分离。单段式壳体结构是考虑到当非起爆端的径向稀疏波对中间段壳体的影响不是很大时,在一体式壳体结构基础上保持药柱长度不变而直接缩短壳体长度,减少加工时间和成本。

      采用AUTODYN有限元分析软件建立上述4种柱壳装药结构的仿真计算模型,以一体式结构的参数设置为例来说明有限元模型的建立:炸药和壳体的模型长度均为160 mm,其中,炸药半径为20 mm,网格边长为0.5 mm壳体外壁半径为25 mm,厚度为5 mm,其网格相较于炸药进行适当加密,网格边长为0.2 mm网格类型均为Lagrange网格。为统计不同位置处的壳体膨胀半径,从壳体端部开始沿轴向每隔8 mm设置一个观测点,共计21个。起爆点设置在左端面中心,建立如图2所示的仿真计算模型,其余3种结构的模型建立方法相同。

      图  2  一体式壳体的仿真计算模型(1~21为观测点)

      Figure 2.  The simulation calculation model for the integrated shell (1−21 are observed points)

    • 仿真计算模型中的壳体材料为40CrMnSiB钢,热处理方式为860 ℃淬火2 h后再500 ℃回火2 h,內爆加载下材料压力和比体积之间的关系用Mie-Grüneisen状态方程进行描述,由于铁元素质量百分比相近,其状态方程参数选择用AUTODYN材料库中1006钢的参数近似代替。材料采用Johnson-Cook强度模型来描述其在变形过程中的强度变化,该模型定义的屈服应力表述如下:

      $\sigma = \left( {A + B{\varepsilon ^n}} \right)\left( {1 + C{\rm{ }}\ln{\rm{ }}{{\dot \varepsilon }^*}} \right)\left[ {1 - {{\left( {\frac{{T - {T_0}}}{{{T_{\rm m}} - T}}} \right)}^m}} \right]$

      式中:σ为流动应力,$ \varepsilon $为等效塑性应变,$ {\dot{\varepsilon }}^{*} $为无量纲塑性应变率,TT0Tm分别为材料温度、环境温度和材料熔点,ABnCm为材料特性常数。40CrMnSiB钢在上述热处理条件下的特性参数:A=1080 MPa,B=308 MPa,n=0.282,C=0.0175,m= 0.53,T0=293 K,Tm=1793 K。炸药选取为8701炸药,采用JWL状态方程进行描述:

      ${{{p}}_{\rm{e}}}{{ = }}{{{A}}_{\rm{e}}}\left( {1 - \frac{\omega }{{{{{R}}_1}{{V}}}}} \right){{\rm{e}}^{ - {{{R}}_1}{{V}}}} + {{{B}}_{\rm{e}}}\left( {1 - \frac{\omega }{{{{{R}}_2}{{V}}}}} \right){{\rm{e}}^{ - {{{R}}_2}{{V}}}} + \frac{{\omega {{E}}}}{{{V}}}$

      式中:pe为爆轰产物压力,AeBeR1R2ω为炸药的特性常数,V为相对体积,E为炸药单位体积的初始内能[18]。8701炸药的特性参数:Ae=542.4 GPa,Be=7.628 GPa,R1=4.2,R2=1.1,ω=0.3,E=9 GPa。

    • 图3给出了炸药起爆后30 μs时刻4种结构壳体的径向膨胀位移,从图中可以看出:从起爆端开始的0~96 mm区间内,4种壳体的径向位移并无明显差别;而在96~112 mm区间上,裸露药柱的壳体受到轴向稀疏波的影响较大,其径向位移明显减小。因此,为了减小非起爆端对中间段壳体冻结回收的影响,需将一体式壳体改为刻槽式或两段式壳体。由于两段式壳体可直接由一体式壳体经线切割加工后拼接而成,考虑到壳体实际加工的便捷性,确定最佳壳体结构为两段式。改进后的两段式结构既能保证回收段壳体的膨胀外形与一体式结构相同,又能减小非起爆端对回收段壳体的影响,因此满足改进设计的要求。

      图  3  起爆后30 μs时刻4种壳体的径向位移

      Figure 3.  Radial displacements of four shells at 30 μs after detonation

    • 根据采用前照明分幅照相技术拍摄到的金属柱壳膨胀断裂图像[19]可知,在柱壳最大膨胀半径与非起爆端初始膨胀位置之间,存在一个近似线性区,该区域内的壳体膨胀半径与轴向距离之间近似线性相关,如图4所示。区域内靠近起爆端的壳体已经发生断裂,而靠近非起爆端的部分则尚处在膨胀变形阶段,因此定义该区域为金属柱壳在内爆加载下的断裂过渡区,将断裂过渡区的膨胀态金属柱壳作为预期回收的目标。

      图  4  壳体断裂分区

      Figure 4.  Division of fracture area of the shell

      为获得起爆后不同时刻的膨胀态金属柱壳,根据一体式壳体自由膨胀形态随时间的变化,选定了3个回收时刻,结合3个时刻下壳体膨胀变形的初始位置和上一节中筛选出的最优壳体结构,确定了两段式结构的切割、粘接位置,匹配不同时刻下预期回收段壳体的膨胀外径设计了相应的冻结回收装置,数值模拟冻结回收过程并初步评估回收效果。

    • 冻结回收装置的内壁形状尺寸需根据金属柱壳在自由膨胀状态下的变形情况来确定。根据数值仿真的结果,可知起爆后0~23 µs内,炸药内部爆轰波经历了从起爆端向非起爆端传播的过程,壳体沿轴向随之发生不同程度的变形响应。为获取处于不同发展状态下的轴向裂纹传播段,选取起爆后13、16和23 µs等3个时刻的壳体进行冻结回收。图5反映了一体式壳体在内部爆炸载荷作用下于13、16、23 µs时刻的变形响应,其中8701药柱的尺寸为$\varnothing $40 mm×160 mm,40CrMnSiB钢柱壳厚度为5 mm。由图可知,13、16、23 µs时刻预期回收的中间段壳体分别介于高斯点5~13、6~16、7~21之间。

      图  5  一体式壳体的膨胀变形响应压力云图

      Figure 5.  Variation of internal pressure of integrated shell with expansion and deformation process

      为避免非起爆端壳体对中间段壳体冻结回收的影响,加工时将13、16 µs回收时刻的原壳体分别沿高斯点13、16所在横截面切断,然后粘接为一体,而23 µs时刻的膨胀态壳体从高斯点7到高斯点21呈近似线性膨胀,所以该回收时刻的原壳体加工时无需切断。

      进行冻结回收装置的设计时,将起爆端的冻结回收装置内壁设计成与壳体最大膨胀外径相切的圆柱形,以实现对该段壳体的径向约束,防止其在径向载荷的作用下将冻结段壳体带出回收装置;中间段壳体对应的冻结回收装置内壁则根据壳体膨胀外径与轴向位置的线性关系,设计成与壳体膨胀外形相匹配的圆台形;将非起爆端对应的回收装置内壁设计成扩口形状,利用内壁形状的突变在两段壳体的交界面处形成应力集中,加速非起爆端壳体的脱离。匹配13、16、23 µs时刻壳体膨胀外形设计的冻结回收装置如图6所示。

      图  6  匹配不同时刻壳体膨胀外形设计的冻结回收装置结构尺寸

      Figure 6.  Structural dimensions of designed freezing recovery devices matching deformation responses of shells at different times

    • 通过对柱壳装药在冻结回收装置约束作用下的静爆过程进行数值模拟,得到三个时刻下中间段壳体的径向膨胀位移随时间的变化情况,如图7所示。壳体从起爆端开始沿轴向逐渐发生膨胀变形,当自由膨胀进行到回收时刻时,从图中膨胀位移转折点可以看出,中间段壳体的轴向各点几乎是同时与回收装置内壁发生接触,在冻结回收装置对壳体的持续约束作用下,膨胀速度逐渐减小,最终壳体的膨胀运动完全停止,壳体回收的动态响应过程也说明了冻结回收装置内壁结构设计的合理性。

      图  7  3个回收时刻壳体径向变形响应

      Figure 7.  Radial deformation response of the shell at three recovery times

      图8给出了3个时刻下中间段壳体的最终回收半径与理想回收半径之间的关系,发现在壳体与冻结回收装置发生接触后,会继续进行不同程度的径向膨胀,引起冻结回收装置的塑性变形。其中,13、16、23 µs时刻的回收壳体在理想回收半径的基础上分别继续向外膨胀了4.7~8.2、4.8~8.8、4.0~5.4 mm,回收到的壳体分别更接近19.2、21.4、26.8 µs时刻的自由膨胀态。此外,不同轴向位置处壳体在与回收装置接触后的膨胀距离并不一致,这是由于起爆后爆轰产物向非起爆端流动堆积,使得越靠近非起爆端的冻结回收装置内壁所受到的压力载荷越大,造成壳体的偏移膨胀距离从起爆端到非起爆端呈现逐渐增大的趋势。

      图  8  3个时刻下的壳体回收半径

      Figure 8.  Recovery radii of the shell at three recovery times

    • 在对冻结回收过程进行数值模拟的基础上,为进一步验证冻结回收试验方法对于回收膨胀态金属柱壳的可行性与合理性,按照上述方案设计进行了内爆加载下金属柱壳的冻结回收试验。图9为金属柱壳装药结构的实物图,加工13和16 µs时刻回收的两段式结构壳体时,分别在距离起爆端96和120 mm位置处采取线切割的方法沿横截面切断,然后用胶水将两段壳体粘接成一体,23 µs时刻回收壳体则根据预期回收段壳体的膨胀断裂特点直接采用一体式结构,内部装填直径40 mm的8701药柱;壳体起爆端和非起爆端设计了螺纹段,以实现壳体与尼龙雷管座和端盖之间的稳固配合,其中尼龙雷管座匹配壳体和冻结回收装置的端口形状进行了定位设计,保证壳体在横向放置时的定位可靠。冻结回收试验的整体布置如图10所示。

      图  9  金属柱壳装药结构

      Figure 9.  Charge structures of the metal cylindrical shell

      图  10  冻结回收试验总体布局

      Figure 10.  Overall layout of the freeze recovery test

      按照设计的回收试验方案,对相应时刻处于膨胀断裂中间状态的壳体进行了冻结回收,这里选取23 µs时刻的壳体回收结果进行展示和量化分析。从图11中可以看出,设计的回收装置成功约束住了壳体的膨胀运动,得到的膨胀态壳体根据外形的差异可以分为柱面段和锥面段,分别与相应位置处对应的回收装置内壁形状相同,表明壳体的自由膨胀与回收阶段都在回收装置的轴向约束区内进行。将回收得到的壳体与图5所示的自由膨胀态对比,可知两者在预期回收段上的外形是相似的。

      图  11  23 µs时刻的壳体冻结回收试验结果

      Figure 11.  Shell obtained from freezing recovery test at 23 μs

      为进一步评估冻结回收试验的仿真设计与实际试验回收结果之间的差异,对23 µs时刻回收壳体的轴向和径向尺寸进行了测量。其中,壳体柱面段的轴向设计长度为40.78 mm,实际回收得到的长度为38.3 mm,比设计长度短2.48 mm,误差百分比为6.1%;锥面段轴向设计长度为119 mm,实际回收得到的长度为109 mm,比设计长度短10 mm,误差百分比为8.4%。图12反映了实际回收壳体的膨胀半径与其仿真值、设计时刻的理想值和26.8µs时刻的自由膨胀值之间的差异,可以看出几者在膨胀形态上是一致的。其中,壳体的实际回收半径略小于仿真回收半径,但整体吻合程度较高;由于冻结回收不可能瞬间完成,故壳体的实际回收半径要大于23 µs时刻的理想膨胀半径,轴向各点处的绝对误差在3.27~4.57 mm范围内,误差百分比小于10%,与26.8 µs时刻的自由膨胀值相近。评估结果表明,设计的冻结回收试验方法对于回收膨胀态壳体具有较好的效果,回收壳体的轴向和径向尺寸误差相较于其理想状态均可控制在10%以内。

      图  12  23 µs时刻中间段壳体的回收半径

      Figure 12.  Recovery radii of the middle section shell at 23 µs

      从冻结回收试验的结果不难看出,自由膨胀态的金属柱壳在与冻结回收装置接触后,其径向膨胀运动不可能立刻停止,在回收装置的约束作用下仍会继续向外膨胀一段距离。本文是从回收壳体的轴向和径向尺寸与设计理想值之间的近似程度出发,说明冻结回收方法设计的合理性,但对壳体与回收装置间的接触作用导致壳体在断裂形貌上发生的变化方面仍研究不足。为进一步评估回收过程的接触作用所带来的影响,后续可将回收得到的膨胀态壳体的表面失稳分布、裂纹扩展路径、断口形貌特征等断裂信息与自然破碎下的情况进行对比分析,通过比较断裂信息的异同来反映接触作用的影响程度。

    • 基于AUTODYN有限元分析软件对不同柱壳结构的壳体膨胀进行数值模拟,确定了应用于內爆加载下膨胀态金属柱壳回收的最优壳体结构。根据不同时刻下壳体的自由膨胀状态设计与之匹配的冻结回收装置,并进行回收过程的数值模拟和试验验证,最后对试验回收结果进行量化评估,得出结论如下:

      (1)两段式结构最有利于减小非起爆端对预期回收段壳体的影响,其膨胀特性最接于一体式壳体且加工方便,适用于进行内爆加载下膨胀态金属柱壳的冻结回收试验。

      (2)回收得到的壳体形态较为完整,壳体的轴向和径向尺寸与其设计理想值之间的误差均小于10%,在考虑回收装置塑性变形量的前提下,回收结果满足设计需求。

      (3)通过冻结回收试验得到的膨胀态金属柱壳可为研究内爆加载下的壳体膨胀断裂过程提供实物参考,且相比于常用回收破片的方式更能体现壳体膨胀断裂过程中的整体形态演变。

参考文献 (19)

目录

    /

    返回文章
    返回