飞行器结构表面被高空核爆产生的软X射线辐照，瞬时沉积大量能量，导致结构表层材料发生熔化、气化等现象。这些熔化或气化的表层材料背离结构面向外喷射，导致结构表面受到物质喷射时的反作用冲量作用，引起结构动态响应^[1]。在缺少大型X射线模拟源的情况下，化爆模拟技术是目前研究空间结构在软X射线辐照下动态响应的重要技术手段^[2]。

核爆软X射线辐照在飞行器表面形成的冲量载荷具有典型的余弦分布特征，余弦分布载荷的化爆模拟技术所使用的加载源一般有柔爆索、光敏炸药和薄片炸药等。Lindberg等^[3-4]和赵国民等^[5]发展了柔爆索加载模拟技术，按照一定的优化设计实现了近似余弦分布的爆炸载荷。20世纪70年代，Benham等^[6]建立了光敏炸药加载模拟技术，通过在全尺寸再入弹头上喷涂一定厚度的光敏炸药来模拟它在大气层外软X射线作用下的动态响应。Rivera等^[7-9]和Covert等^[10]进一步发展了光敏炸药驱动飞片加载壳体的实验技术，同时对光敏炸药的喷涂系统^[11]和测试系统^[12]进行了升级改进。近几年，光敏炸药加载技术在炸药合成技术、喷涂技术、起爆技术和冲量标定技术方面也取得了重要进展^[13-17]。薄片炸药加载技术在20世纪70年代发展起来^[18]，此技术把炸药条按照一定规则离散并放置在壳体表面一定距离处，在炸药条和壳体之间放置泡沫橡胶等软材料对爆炸载荷进行匀化，以实现在壳体表面施加近似余弦分布的脉冲载荷^[19-21]。薄片炸药由于加工简便、安全性高、起爆可靠，已经成为飞行器抗软X射线辐照考核的重要加载手段。

近年来，随着抗X射线防护技术的发展，相同辐照环境下，软X射线在新型空间飞行器结构表面的沉积量大幅降低，这对飞行器抗软X射线辐照考核技术提出了新的要求，即高均匀度、高同步性和低比冲量。对于传统炸药条加载，它的条状排布、滑移爆轰和橡胶匀化的加载方式，难以满足新型飞行器考核的需求，原因包括：(1)炸药条排布间距越小，总条数越多，均匀性越好，但比冲量越大，难以同时满足高均匀度和低比冲量的要求；(2)一端起爆、滑移爆轰加载，导致结构不同部位加载时间不同步，难以满足高同步性的要求；(3)在壳体质量较小或刚度较低时，采用泡沫橡胶进行载荷匀化，引入的附加质量容易引起结构响应的失真；(4)飞行器结构较复杂时，条状布药缺乏设计灵活性。

本文中，为适应新型空间飞行器结构考核的复杂构型、高均匀度、高同步性和低比冲量载荷设计要求，针对传统炸药条加载方式的不足，提出用十字形超细药条离散群同步起爆实现超低比冲量的加载方式；基于自研的太安基薄片炸药，制作最小截面为0.50 mm×0.33 mm的系列炸药条，并压制组装精密的十字形离散片炸药组件；通过平面多点同步起爆离散片炸药实验，验证该加载方式的可靠性和同步性。通过离散片炸药爆炸流场数值模拟，分析平面布药条件下比冲量空间分布特征和传播演化规律，以期为基于薄片炸药离散群的余弦载荷设计提供参考。

1.   实验方法

1.1   十字形超细药条组件制备

为实现基于炸药条的低比冲量加载，需减小药条的宽度，但药条宽度越小，可稳定传爆的难度越大。通过改进压药和裁剪工艺，制作了可稳定爆轰的厚度为0.33 mm、长度为10 mm、宽度分别为0.5、0.75、1.0、1.25 mm的炸药条，如图1所示。每2根药条与1根直径为0.5 mm的银皮柔爆索（线装药密度为0.09 g/m）通过模具压制成十字形片炸药组件，制作过程和实物如图2所示：首先将柔爆索与金属支架固定，然后将2根短炸药条垂直交叉，通过模具压紧在金属支架上，从而实现柔爆索起爆端与片炸药中心的无缝连接。

图  1  不同宽度片炸药条

Figure  1.  Sheet explosive rods with different widths

下载: 全尺寸图片 幻灯片

图  2  十字形片炸药组件及其制作过程

Figure  2.  Cross sheet explosive components and their manufacturing process

下载: 全尺寸图片 幻灯片

1.2   分布式布药方法

为了实现给定空间构型加载，将十字形炸药布设在特定形状的定位板上，定位板的形状决定了载荷的初始空间构型。对于平面加载，定位板形状为平面壳；对圆柱壳进行加载时，定位板形状为圆柱壳。定位板上预制与片炸药分布相同的定位孔，将柔爆索穿过定位孔，并将片炸药粘贴在定位板上，即可实现片炸药的空间定位。

使用平面加载构型，共使用21支十字形片炸药组件，布药间距d=25 mm，如图3所示。若以截面尺寸相同的0.5 mm宽炸药条按滑移爆轰加载方式条状排布，则相同布药面密度下的条状布药间距为15.625 mm。为对比2种布药方式的均匀度，用能够连续排布并铺满整个二维平面的正方形将布药平面分割成M个布药单元，将可使每个布药单元内的含药量都相同的最小正方形称作最小布药单元。因此，可以用单位长度内的最小布药单元数来衡量布药均匀度。如图4所示，十字药条排布情况下，最小布药单元边长为8.8388 mm，布药均匀度为113.1 m⁻¹；条状排布情况下的最小布药单元边长与布药间距相同，即15.625 mm，布药均匀度为64 m⁻¹。因此，相比条状布药，十字布药的空间均匀度提高了76.7%。

图  3  十字形片炸药离散起爆阵列

Figure  3.  Cross sheet explosive array for discrete detonation

下载: 全尺寸图片 幻灯片

图  4  十字形布药与条状布药的均匀度对比

Figure  4.  Comparison of the evenness between cross and strip distributions

下载: 全尺寸图片 幻灯片

为提高加载的时间同步性，分布式布药采用多点同步起爆的方式。21根柔爆索汇聚至1个一入多出的传爆接头，并由1根直径1.0 mm的铅皮柔爆索（线装药密度为0.2 g/m）同步起爆，降低了传统滑移爆轰加载一端起爆的时间不同步性。该传爆模式可逐级扩展，最终实现由1根导爆索逐级同步起爆成百上千个薄片炸药组件，从而应用于大型结构的加载。

1.3   测试系统

主要测试包括基于电探针的起爆同步性测试、起爆过程的高速摄影以及加载比冲量测试，测试控制系统如图5～6所示。由于爆炸在微秒级的瞬间完成，测试系统需高度的时间同步。

图  5  测试系统原理图

Figure  5.  Schematic diagram of test system

下载: 全尺寸图片 幻灯片

图  6  测试系统实物图

Figure  6.  Physical photo of test system

下载: 全尺寸图片 幻灯片

图  7  起爆阵列的安装方式

Figure  7.  Installation mode of the explosive array

下载: 全尺寸图片 幻灯片

由同步起爆器发出3路ns级同步控制信号：第1路用于起爆雷管，并依次引爆1根一级柔爆索、21根二级柔爆索、21片十字形片炸药；第2路用于触发一号示波器并记录3路探针信号，其中1路探针贴于雷管外壁，用于触发高速相机，其余2路探针贴于片炸药阵列两对角处，用于衡量起爆同步性；第3路用于触发二号示波器并记录冲量摆测量系统的输出信号。

高速相机为自研八通道分幅高速相机，用于捕捉爆炸瞬间图像。高速相机各通道曝光时间为0.5 μs，最多记录时长为4.0 μs。高速相机由雷管起爆导通电探针产生的电信号触发，雷管起爆时刻距离薄片炸药起爆时刻间隔着柔爆索传爆时间。2级柔爆索总长度为300 mm，以爆速6.7 km/s估算，传爆总时间约44.7 μs。因此，设置高速相机触发后延迟44 μs曝光。

冲量摆测量系统^[22-23]用于测定爆炸产生的比冲量，摆锤放置于距片炸药垂直距离40 mm处，两摆锤平行放置，一个正对药片阵列中心，另一个中心偏左30 mm。

图3所示的定位板通过顶角处的4个安装孔固定于一块中心开孔的立板中心，立板尺寸约2 m×2 m，中心开孔尺寸略小于定位板尺寸。立板的一侧分布起爆组件，包括雷管、导爆索；立板的另一侧分布测试系统（见图6～7），包括比冲量测试系统、高速摄影测试系统和同步探针测试系统等。立板垂直固定于地面，用于隔绝雷管和柔爆索的爆炸能量对比冲量测试的影响。

2.   实验结果及分析

2.1   多点起爆同步性

实验中21根十字形片炸药组件全部起爆，起爆率100%。如图8所示，出于超低比冲量的设计要求，所用药量很小，起爆能量极低，金属支架的十字支撑端（见图2）仅在铝合金材质的定位板上产生了微米级极浅压痕，部分位置由于与定位板贴合不够紧密，压痕甚至不甚明显。

图  8  实验后定位板样貌

Figure  8.  Appearance of the positioning plate after experiment

下载: 全尺寸图片 幻灯片

图9为高速相机捕捉到的起爆瞬间，亮斑即为冲击波发光，光斑阵列中的圆形及条状阴影为冲量摆锤及摆杆遮光所致。在第3个通道曝光之间内，亮斑开始出现；在第4个通道曝光时间内，21个亮斑完全出现。各通道的曝光时间均为0.5 μs，说明起爆的不同步性小于1 μs。左下角的黑色阴影为固定探针所用胶带的影响。

图  9  高速摄影捕捉爆炸瞬间

Figure  9.  Explosion moment captured by high-speed photography

下载: 全尺寸图片 幻灯片

探针信号给出了系统各事件的发生时间，如图10所示。以起爆器发出的触发信号为零时，探针1的信号时刻为7 578.96 μs，即雷管起爆延迟时间，与所用雷管性能相符；探针2与探针1信号时刻相差45.02 μs，即两级柔爆索传爆所需总时间，与事先估算的44.7 μs基本一致；直接表征起爆同步性的探针2与探针3的信号时差为0.38 μs，与高速摄影同步性判读结果映证。导致起爆同步性误差的原因包括：(1)柔爆索与片炸药连接点情况，影响起爆过程；(2)柔爆索长度误差及柔爆索爆速的不稳定性，影响从雷管到片炸药的传爆过程。

图  10  探针信号

Figure  10.  Signals of probes

下载: 全尺寸图片 幻灯片

2.2   比冲量测试结果

由动量矩定理、机械能守恒定律，可推出冲量摆的比冲量计算公式为：

式中：I为待测比冲量，S为载荷作用面积，L为载荷对转轴的力臂，J为冲量摆系统对转轴的转动惯量，T₀为冲量摆系统固有周期，θ_m为冲量摆最大摆角，m和l分别为冲量摆系统质量和质心距，g为重力加速度。

本系统摆锤测量面直径约为15 mm，中心和左侧摆锤质量分别为68.9和69.1 g，摆杆长200 mm，摆杆质量为2.3 g。图11为实测冲量摆系统输出的电压信号，满量程±2.5 V对应±180°摆角。综合考虑系统各部件转动惯量，计算得出中心冲量摆所测比冲量为79.9 Pa·s，左侧冲量摆所测比冲量为61.4 Pa·s，左侧冲量较小是由边侧稀疏效应造成的。测试误差可由下式计算：

图  11  冲量摆系统输出信号

Figure  11.  Output voltage signals of the impulse pendulum system

下载: 全尺寸图片 幻灯片

通过简化，认为转动惯量J、摆杆长度L和摆锤面积S引入的不确定度可以忽略，则式(2)可以写为：

所用编码器为12位绝对编码器，其角度分辨率为0.00153°，时间分辨率为0.01 s；2个冲量摆测得的最大摆角为8.65°，冲量摆周期为0.9 s，则$\left| {\dfrac{{{\text{d}}I}}{I}} \right| {\text{＜}} 3.1$%，该误差忽略了系统阻尼的影响。

3.   数值模型及验证

实验测试可以给出有限区域的比冲量信息，但很难给出比冲量在空间上的分布及其传播演化规律，更丰富的信息可通过数值计算获得。为此，依据实验建立了如图12所示的数值模型。由于实验中定位板仅产生了微米级的极浅压痕，加载过程可近似为理想反射加载；从实测结果看，冲量摆的固有周期为0.9 s，远长于爆炸加载过程（约几十微秒），且金属靶亦未发生塑性变形，应力波在靶体中的传播过程可忽略。因此，模型忽略定位板和金属靶中的应力波效应，将其简化为刚体，采用Lagrange单元描述；炸药和空气采用ALE单元，以共节点形式实现物质传输，并与Lagrange单元在空间位置有重合，以便程序自动搜索流固耦合边界；空气域四周为无反射边界。

图  12  数值计算模型

Figure  12.  The numerical model according to the experiment

下载: 全尺寸图片 幻灯片

空气采用多项式状态方程描述：

式中：p为空气压力，γ为比热比，ρ和ρ₀分别为密度和初始密度，E为空气比内能。取值分别为：γ=1.4，ρ₀=1.29 kg/m³，E=2.5 MPa。

炸药爆轰产物采用JWL状态方程描述：

式中：p为爆轰产物压力，V为爆轰产物相对比容，E为爆轰产物比内能，A、B、R₁、R₂、ω为常数。参数取值参考手册UCRL-52997中PETN 1.77^[24]，分别为：A=617 GPa, B=17 GPa, R₁=4.4, R₂=1.2, ω=0.25；初始比内能为10.1 GPa；初始相对比容为1。炸药其余参数取值为：密度ρ=1770 kg/m³, 爆速D=8300 m/s, CJ爆轰压力p_CJ=33.5 GPa。

模型的布局和尺寸与实验一致，并做了2处简化：

(1)简化了冲量摆系统。将它等效为等尺寸的金属质量块，以质量块速度信息计算比冲量。

(2)简化了片炸药尺寸。由于片炸药短边尺寸仅有0.5 mm，则片炸药与周边网格尺寸需小于0.5 mm，这将使模型的单元数量变得非常巨大而降低计算效率，因此，在此将炸药短边尺寸变为1 mm，同时调整药片厚度为0.174 mm，使总药量相同。

在爆炸加载问题中，网格尺寸对计算结果往往有较大影响，且所需的网格尺寸与炸药当量关系较大^[25-26]。为此，建立了网格尺寸分别为2.0、1.0和0.5 mm的模型（见图13），炸药厚度方向均设置两层网格。计算直至质量块速度v趋于恒定，则质量块局部比冲量为I=ρvl，ρ、l分别为质量块的密度和长度。

图  13  三种网格尺寸

Figure  13.  Three types of mesh sizes

下载: 全尺寸图片 幻灯片

3种网格尺寸下的总ALE单元数分别为564万、75.3万、7.73万，随网格尺寸减小，计算耗时几何增长，而计算结果越来越趋近实测值（见图14）。当网格尺寸为2.0 mm时，计算结果比实测值大30%～40%；当网格尺寸为1.0 mm时，误差小于15%；当网格尺寸为0.5 mm时，误差小于5%。由此可知，网格尺寸偏大会使计算比冲量偏大，在该尺寸薄片炸药计算中，应当选取0.5 mm的网格尺寸。

图  14  不同网格尺寸的计算结果与实测值的对比

Figure  14.  Comparison of experimental and numerical results under various mesh sizes

下载: 全尺寸图片 幻灯片

4.   比冲量特性分析

4.1   数值模型和评价方法

为探究比冲量在空间的分布和传播演化规律，建立了如图15所示的数值模型，图中黑色网格标记部分为代表性压力监测单元，模型网格尺寸为0.5 mm，布药间距d分别为20、25和30 mm，其余参数与图12中模型一致。由于药片布设具有周期性，在忽略边侧稀疏波影响的前提下，代表性压力监测单元可表征布药阵列的比冲量匀化过程。

图  15  比冲量演化规律研究数值模型

Figure  15.  The numerical model for studying the evolution of specific impulse

下载: 全尺寸图片 幻灯片

不同z向距离处的峰值比冲量I(z)和比冲量均匀度偏差U(z)分别为：

式中：N为z层（单元中心z坐标为z）监测单元总数，I_i(z)为z层第i个监测单元的峰值比冲量。U(z)越小，代表均匀度越高。

4.2   比冲量匀化特性

图16为布药间距为25 mm时不同z层峰值比冲量的分布。在同一层监测单元内，峰值比冲量分布关于2个对角线具有较好的对称性，这与布药方式的几何对称性一致。z较小时（以z=2 mm为例），峰值比冲量分布与初始布药位置相对应，药片附近的峰值比冲量较大，其余位置峰值比冲量较小；随z的增大（以z=15 mm为例），由于冲击波在传播过程中的相互碰撞、叠加，峰值比冲量发生了重新分布，呈现中心大、四周小的形态；随z继续增大（以z=30 mm为例），冲击波在多次碰撞、叠加后进一步匀化，峰值比冲量的分布已较均匀。

图  16  布药间距25 mm时比冲量空间分布及演化规律

Figure  16.  Spatial distribution and evolution of specific impulse with a spacing of 25 mm

下载: 全尺寸图片 幻灯片

各层峰值比冲量的计算结果如图17(a)所示，考虑匀化过程与布药间距有关，横坐标采用无量纲量z´=z/d表示；峰值比冲量与布药面密度α相关，纵坐标以I´=I/α^1/3表示。此时，不同布药面密度的峰值比冲量演化曲线几乎重合。这表明，在该周期性离散布药条件下，峰值比冲量由布药面密度决定，峰值比冲量随传播距离的演化过程由布药间距决定。

图  17  不同药片间距下比冲量演化规律

Figure  17.  Evolution of specific impulse under different spacings of sheet explosive

下载: 全尺寸图片 幻灯片

各层峰值比冲量的均匀度偏差如图17(b)所示，横坐标同样采用无量纲量z´=z/d表示。从均匀度偏差的演化特征看，匀化过程可大致分为3个阶段：第1阶段为扩散段，距离起爆点较近，受离散布药的影响，均匀度偏差较大，主要通过冲击波自由扩散匀化，匀化速度最快；第2阶段为叠加段，约起始于0.1倍布药间距附近，通过冲击波叠加、碰撞进行匀化，此时均匀度偏差仍然较高，匀化速度慢于扩散段；第3阶段为均匀段，起始于0.8倍布药间距附近，峰值比冲量的分布已达到较高的均匀度，均匀度偏差小于10%。回看图17(a)，峰值比冲量的演化过程亦可对应分为3个阶段，峰值比冲量在扩散段迅速降低，在叠加段逐渐增大，在均匀段逐渐趋于稳定。

4.3   药片形状对匀化过程的影响

如前所述，峰值比冲量与布药面密度相关，匀化过程与布药间距相关。若要实现给定比冲量加载，在相同的布药面密度下，布设点数越多，匀化越快，但同时会带来制作难度增高、熄爆点出现概率增大的问题。因此，主要讨论在相同布药面密度和布设点数下，药片形状对匀化过程的影响，此时单个药片的表面积不变（药片厚度恒定）。

十字形药片的特点是在4个方向上各有一个分支，因此，从分支数量和分支长度2个维度来考虑药片形状的改变。如图18所示，一种方法是增加分支数量，减小分支长度，极限情况是，分支数无限多，形成圆形药片。为方便划分网格，数值模型中将药片设计为方形。另一种方法是减少分支数量，增大分支长度，极限情况是，分支数只有一个，形成短条形药片。图19为3种类型药片的数值模型：方形、十字形和短条形。3种模型的布药面密度相同，图中黑色网格标记处为代表性压力监测单元。

图  18  不同形状药片的设计思路

Figure  18.  Design ideas for sheet explosives with different shapes

下载: 全尺寸图片 幻灯片

图  19  不同药片形状的数值模型

Figure  19.  Numerical models of different shapes of sheet explosive

下载: 全尺寸图片 幻灯片

如图20所示，与十字形药片类似，方形和短条形药片阵列的峰值比冲量和均匀度偏差演化过程亦可大致分为3个阶段，此为阵列式布药的共性特征。由于布药面密度相同，三者的峰值比冲量曲线较接近，但匀化过程差异较大。扩散段，3种阵列初始均匀度偏差基本相同，由于形状的差异，方形药片阵列最快进入叠加段，短条形最慢，十字形处于二者之间。叠加段，方形均匀度偏差最大；短条形均匀度偏差最小，但起伏较大，进入均匀段的速度与方形相当，均在1倍布药间距附近；十字形均匀度偏差起伏较小，最快进入均匀段，约在0.8倍布药间距处。均匀段，短条形和方形均匀度偏差相当，曲线几乎重合，十字形均匀度偏差最小。以进入均匀段的速度和均匀段的匀化程度来看，十字形最优，方形和短条形接近。

图  20  不同形状药片阵列的峰值比冲量和均匀度偏差演化过程

Figure  20.  Evolution of peak specific impulse and uniformity deviation of sheet explosive arrays with different shapes

下载: 全尺寸图片 幻灯片

同理，进一步开展与传统条状排布滑移爆轰加载形式的对比分析。在保持布药面密度、药条宽度和厚度相同的前提下，该加载方式的布药间距为15.625 mm。如图21(a)所示，在药条一端起爆，选取紧邻y=0平面的一层单元作为代表性压力监测单元。不考虑滑移爆轰所带来的时间不同步性，如图21(b)所示，条形布药匀化过程与短条形较类似，但由于布药间距较大，均匀度偏差起伏较大，进入均匀段较晚，且偏差值高于短条形和十字形。由此可见，相比传统条状排布滑移爆轰加载，使用十字形药条离散群同步加载，不仅提高了时间同步性，而且提高了载荷均匀性。

图  21  条状布药滑移爆轰加载数值模型和比冲量匀化过程

Figure  21.  The numerical model and homogenization process of specific impulse under sliding detonation loading of explosive rods

下载: 全尺寸图片 幻灯片

5.   结　论

使用十字形超细药条离散群同步起爆的方法，实现了高同步性、高均匀度、低比冲量平面加载。验证实验的结果表明：(1)同步性方面，各十字形片炸药起爆同步性误差小于1 μs；(2)均匀度方面，相比炸药条滑移爆轰加载模式，布药空间均匀度提高了76.7%；(3)比冲量方面，实测比冲量低至几十帕秒量级。在此基础上，保证了起爆的可靠性，起爆率达到100%。

为进一步明确该加载方式下载荷比冲量演化特性，建立了离散薄片炸药加载数值模型，模型可靠性得到了比冲量实测结果的验证；选取了代表性压力监测单元进行比冲量匀化过程分析，得出以下结论。

(1)在周期性离散布药条件下，峰值比冲量由布药面密度决定，峰值比冲量的演化过程由布药间距决定。

(2)比冲量匀化过程可分为3个阶段，即扩散段、叠加段和均匀段，扩散段通过冲击波自由扩散匀化，叠加段通过冲击波叠加、碰撞匀化，最终进入比冲量分布较为均匀的均匀段。

(3)方形和短条形药片阵列进入均匀段所需的匀化距离约为1倍布药间距，而十字形药片阵列仅需约0.8倍布药间距，且在均匀段的匀化程度更高，因而在平面布药加载情况下，十字形药片具有较大优势。

(4)离散炸药群同步起爆的加载方式，相比条状布药滑移爆轰加载，不仅提高了载荷同步性，而且提高了载荷均匀度；利用片炸药和结构之间的空气层对载荷进行匀化，可避免橡胶附加质量过大带来的结构响应失真。

本文中仅以平面同步加载作为原理性验证，若进一步改变各薄片炸药所对应柔爆索的长度、改变定位板的形状，可实现特定空间构型、特定时序的爆炸加载，从而应用于复杂形状、复杂工况下的爆炸加载考核。