内爆驱动式超高速发射技术的初步研究

王马法，HIGGINS Andrew J，焦德志，黄 洁，柳 森

（1. 中国空气动力研究与发展中心，四川 绵阳 621000；2. 麦吉尔大学工程力学系，魁北克 蒙特利尔 H3A2K6，加拿大）

随着空间科学研究的发展，为进一步开展空间碎片毁伤效应、防护结构以及超高速碰撞等研究工作，提出了10 km/s 以上的超高速发射需求。国内外常用的发射技术（设备）包括：二级/三级轻气炮、磁驱动飞片、聚能射流、等离子体加速器、激光驱动飞片装置等[1-6]。其中二级轻气炮是目前应用最广泛、技术最为成熟的发射设备。然而，受到炮体材料性能的限制，二级轻气炮的稳定发射速度基本都低于8 km/s，如果超过该速度，发射管将发生严重的烧蚀破坏[7]。虽然有些发射技术能够实现8 km/s 以上的超高速发射，但是这些发射技术要么难以控制弹丸的质量和形状，要么只能发射很薄的飞片状弹丸或质量很小的弹丸[8]，因此将其应用到超高速碰撞实验研究中受到了一定的限制。在20 世纪60 年代，内爆驱动式超高速发射技术已经实现了2 g 弹丸12.2 km/s 的发射速度[9-10]；进入21 世纪后，加拿大麦吉尔大学[1,11-13]又实现了15 g 弹丸7.5 km/s 左右的发射速度，这是一项发射质量和速度都超过二级轻气炮的超高速发射技术，能够满足空间碎片撞击的研究需求。

内爆驱动式超高速发射技术（国外有Explosive hypervelocity launcher[9]、Explosively driven gun[10]、Implosion-driven launcher[11-13]等不同说法，国内也有炸药爆轰驱动高速激波管发射技术[14]、炸药驱动枪[15]等不同名称，本文以Implosion-driven launcher 说法为准，简称内爆发射器）是一种利用炸药爆炸线性挤压压缩管内轻质气体，通过被压缩的轻质气体驱动弹丸的一种发射技术[11]，作用过程与二级轻气炮中活塞压缩高压轻质气体类似，不同的是其活塞速度为炸药爆速。国内一些文献中对该技术进行了相关介绍，但尚未见相关实验和仿真研究报道[14-15]。

为了掌握内爆发射技术，获得10 km/s 以上的发射能力，通过与加拿大麦吉尔大学的Higgins 教授合作，以8 mm 口径内爆发射装置为研究对象，对内爆发射器的发射能力和研究方法进行了初步探索，开展了内爆发射器内弹道分析和结构设计，并利用AUTODYN 软件对设计结构进行了初步仿真分析，最终针对典型状态开展了初步验证实验。

1 内爆发射器工作原理

内爆发射器主要由炸药线性压缩段、高压段和发射管组成，其中炸药线性压缩段包括外壳、装药、压缩管、驱动气体等，发射器的结构如图1 所示，作用过程如图2 所示。左端炸药起爆后，压缩管向内塌缩，形成锥形结构，当炸药爆轰波传播时，锥形结构向右运动并压缩驱动气体。该过程与二级轻气炮中活塞压缩高压轻质气体类似，因此该锥形结构可看成虚拟活塞，当冲击波到达弹丸底部时，高压驱动气体开始推动弹丸加速运动。

图1 内爆式超高速发射器结构示意图[11]Fig. 1 Structural diagram of implosion-driven hypervelocity launcher[11]

图2 内爆式超高速发射器工作原理[11]Fig. 2 Working diagram of implosion launcher after ignition[11]

2 内爆发射过程的仿真分析

设计了内径16 mm、长度500 mm 的压缩管和内径8 mm、长度300 mm 的发射管组成的内爆发射器，高压段采用锥段连接压缩管和发射管，锥段半锥角度为5°。采用AUTODYN 软件对该内爆发射器的工作过程进行了仿真计算。

2.1 仿真模型及方法

采用有限元中的流固耦合方法进行仿真。将模型分为Euler 模型和Lagrange 模型两个部分：Euler 模型包括填充的氦气、压缩管、装药以及周围空气介质；Lagrange 模型包括压缩段外壳、高压段、发射管以及弹丸。外壳、高压段、压缩管、发射管均采用4340 钢。本构模型采用Johnson Cook 模型，状态方程为Linear。装药采用爆速为7 km/s 左右的季戊四醇四硝酸酯（PETN），状态方程为JWL；驱动气体为氦气，采用理想气体状态方程；弹丸本构模型为Steinberg Guinan，状态方程为Shock。材料参数均选自AUTUDYN 数据库。

为了使弹丸能够约束初始状态的高压氦气，在弹丸与高压气体交界处增加一层Lagrange 网格。该网格与高压段和弹丸均采用共节点的方式固定，并将该部分的材料失效应变设置为0.5，当达到失效应变时自动删除该部分的网格，使其不会对弹丸发射造成过大影响。

2.2 仿真计算相关参数

针对弹丸材料和充气压力两个参数，共进行4 个状态的仿真计算，相关参数以及各状态下获得的弹丸出炮口速度见表1。

表 1 不同弹丸材料和充气压力下仿真计算参数Table 1 Simulation parameters of launchers with different projectile materials and filling pressure

2.3 仿真结果

2.3.1 内爆发射器工作过程

图3 给出了内爆式超高速发射装置发射铝弹丸的仿真结果，其中弹丸尺寸为 ∅8 mm × 4 mm，质量为0.55 g。零时刻从装药的左端面开始点火，装药爆炸后开始挤压压缩管形成锥形活塞，并在氦气（初始压力5 MPa）中形成压力达到500 MPa 左右的冲击波向右传播，外壳达到失效应变后被删除，如图3(b)所示，压力超过1 GPa 时显示为红色；在55 µs 左右，冲击波到达弹丸底部反射，形成峰值压力高达5 GPa的气体推动弹丸加速运动；70 µs 左右炸药爆轰结束；在108 µs 左右弹丸离开发射管，达到最高发射速度。

图3 8 mm 口径内爆式发射装置工作过程Fig. 3 Launch process of 8 mm caliber implosion launcher

2.3.2 不同初始充气压力下弹丸的发射速度

图4 给出了充气压力分别为4、5 和6 MPa 状态下铝弹丸的速度-时间历程。从速度-时间曲线上可以看到，弹丸加速过程有较为明显的二次冲击加载。初始冲击加载是由先导冲击波到达弹丸底部产生的；二次冲击加载是弹丸底部反射的冲击波到达左端后再次反射冲击波，该冲击波追赶上弹丸产生的。从表1 中的数据可知，4、5 和6 MPa 充气压力下弹丸的出炮口速度分别为8.62、9.17 和9.25 km/s。可以看到，炮口速度随充气压力的增加而增加，但增加幅度缩小。可能是由于当充气压力从5 MPa 增加至6 MPa 时，压缩管、高压段内壁的膨胀显著增加，削弱了驱动压力增加的增速效果。

图4 不同充气压力下铝弹丸的速度-时间历史Fig. 4 Velocity-time history of aluminum projectiles with different filling pressure

2.3.3 不同弹丸材料的发射速度

根据建立的仿真模型对尺寸为 ∅8 mm × 4 mm的铝合金弹丸和镁合金弹丸进行发射过程的计算，初始填充5 MPa 氦气，计算结果如图5 所示。0.55 g 铝合金弹丸的最终发射速度为9.17 km/s；0.37 g 镁合金弹丸的发射速度为10.66 km/s，比铝合金弹丸速度高1.49 km/s。

图5 不同材料弹丸速度-时间历史Fig. 5 Velocity-time history of projectiles with different materials

3 实验设计及结果

3.1 发射器结构

发射器总长度为914 mm，最大部位直径为81.4 mm，主要包括发射管、发射管护套、高压段、外壳、压缩管、弹丸和装药等。其中：发射管内径为8 mm，长度为304 mm；压缩管内径为16 mm，长度为495 mm。装药为柱壳装药，厚度约为5 mm，装药量约为300 g，采用两层以黑索金（RDX）为基底的片状挠性炸药包裹而成，挠性炸药的爆速为7.1 km/s。弹丸直径为8 mm，高度为4 mm。弹丸材料有两种，分别是铝合金和镁合金，其中铝合金弹丸质量约为0.55 g，镁合金弹丸质量约为0.37 g，安装完成后的发射器如图6 所示。

图6 内爆式超高速发射器实物Fig. 6 Image of implosion-driven hypervelocity launcher

3.2 实验装置总体布局

实验装置布局如图7 所示。实验装置主要包括发射器、测试系统、真空系统和防护系统等。发射器放置在厚壁防护筒内，防止爆炸产生的碎片向外飞散，损伤实验装置和设备。采用激光测速系统和靶网测速装置等测量弹丸的发射速度，并采用序列激光阴影成像仪拍摄模型的飞行姿态和形貌。实验中发射器压缩管内充5 MPa 氦气，从端面起爆装药。

图7 实验装置布局示意图Fig. 7 Layout of experimental equipment

3.3 实验结果与分析

共开展了7 次实验，其中铝合金弹丸实验3 次，镁合金弹丸实验4 次，具体实验参数如表2 所示。

表 2 内爆发射器实验相关参数Table 2 Parameters of the implosion-driven launchers in tests

3.3.1 回收发射器

实验后的发射器压缩管和外壳完全碎裂，高压段外壳明显膨胀变形并出现裂纹，回收的典型发射器装置如图8 所示。

图8 ILT08 实验发射器装置回收Fig. 8 Recycle launcher of test ILT08

3.3.2 铝合金弹丸速度测试结果

3 次铝合金弹丸实验结果如图9 所示，其中横坐标为测点位置离炮口的距离，纵坐标为测得的速度，各次实验测得的最大速度见表2。其中，实验中测得的弹丸最高出炮口速度为7.95 km/s，比仿真预测速度9.17 km/s 低1.22 km/s，相对偏差为15.3%。实验速度偏低可能是许多因素共同影响的结果，例如装置加工误差、装药松紧程度、发射管与弹丸有摩擦阻力等，而仿真状态则相对理想，没有考虑这些因素的影响。另外，实验中弹丸的速度也有1 km/s 左右的偏差，挠性炸药装填、装药爆轰的圆周同步性、发射管与弹丸的配合程度等问题都会对实际发射速度造成一定影响，导致重复实验出现一定的偏差。

图9 0.55 g 铝合金弹丸速度测试结果Fig. 9 Experimental muzzle velocities of 0.55 g aluminum projectile

3.3.3 镁合金弹丸速度测试结果

4 次镁合金有效实验获得的速度结果如图10所示，实验测得的最大速度见表2。4 次实验获得的最高发射速度与仿真预测速度10.66 km/s 非常接近，相对偏差为3.7%。镁合金弹丸在实验中的速度偏差为1 km/s左右，产生偏差的原因可能与铝合金弹丸相同。

图10 0.37 g 镁合金弹丸速度测试结果Fig. 10 Experimental muzzle velocities of 0.37 g magnesium projectile

3.3.4 弹丸成像结果

通过超高速八序列激光阴影成像仪拍摄弹丸在实验过程中的飞行情况，图11 为3 次实验拍摄的照片，其中ILT11、ILT12 为镁合金弹丸，ILT14 为铝合金弹丸。在图11 中可以看到，实验过程中的弹丸仍保持为一个整体弹丸，但存在一定程度变形，其中实验ILT12 的弹丸变形非常明显。发生变形的原因可能是材料强度不够，也可能是发射器结构存在一定的缺陷，有待进一步深入研究。

图11 弹丸序列激光阴影成像结果Fig. 11 Projectile photos shot by sequence laser shadowgraph imager

4 结 论

通过与加拿大麦吉尔大学合作，初步设计了8 mm 口径内爆式超高速发射器结构，采用AUTODYN 软件对内爆式超高速发射器内弹道进行数值模拟，所获得的发射器作用过程中弹底压力最高可达5 GPa，并计算了铝合金弹丸和镁合金弹丸的发射速度。最终通过开展验证实验，证明了设计的内爆发射技术能够将尺寸为 ∅8 mm × 4 mm、质量为0.37 g 的柱状镁合金弹丸发射到10.28 km/s，将尺寸为 ∅8 mm × 4 mm、质量为0.55 g 的柱状铝合金弹丸发射到7.95 km/s。实验获得的发射速度与数值仿真结果吻合较好。实验结果也证明了内爆式超高速发射技术的可行性，在空间碎片的撞击与防护方面具有较强的应用前景。利用超高速八序列激光阴影成像仪拍摄弹丸的飞行情况，获得了弹丸发射后的形状，发现发射后的弹丸是一个整体弹丸，但与原始柱形结构相比，存在不同程度的变形。未来将深入探索提高内爆发射器速度的方法，分析弹丸发生变形的原因，并在保持弹丸初始形状的发射器结构设计等方面开展进一步的研究工作。

感谢中国空气动力研究与发展中心的罗庆、龙耀、宋强、周毅、覃金贵、李俊玲、姜林、邹胜宇、李文光、廖强、刘晓龙、廖富强、丁建文等同事在实验过程中提供的无私帮助。