钝感双基发射药老化迁移机理及动力学过程

梁昊; 丁亚军; 李世影; 赵先正; 肖忠良; LIANG Hao; DING Ya-jun; LI Shi-ying; ZHAO Xian-zheng; XIAO Zhong-liang

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

钝感双基发射药老化迁移机理及动力学过程 PDF

- ORCID：
梁昊 ^1,2
✉
- ORCID：
丁亚军 ^1,2
- ORCID：
李世影 ^1,2
- ORCID：
赵先正 ^1,2
- ORCID：
肖忠良 ^1,2
✉

1. 南京理工大学化工学院，江苏南京 210094； 2. 特种能源材料教育部重点实验室（南京理工大学），江苏南京 210094

中图分类号： TJ55； O63

最近更新：2022-03-03

DOI：10.11943/CJEM2021007

目录contents

摘要

为研究钝感双基发射药老化迁移过程及影响因素，利用钝感剂的正态分布特性及Fickian第二扩散定律的推论，建立了钝感发射药老化迁移导致钝感剂浓度分布及燃烧性能变化的理论模型。对6种不同硝化甘油（NG）及邻苯二甲酸二丁酯（DBP）含量的发射药进行不同温度下加速老化实验，采用密闭爆发器测试了不同老化时间后发射药的燃烧性能，分析了发射药老化迁移的动力学过程及热力学影响。理论分析表明，钝感球形发射药燃气生成猛度积分值的增长率随老化时间呈现出较好的线性增长模型。密闭爆发器结果表明，6种发射药加速老化后的动态活度最大值升高百分比随老化时间均呈较强的线性关系，与理论研究结果一致。建立的钝感发射药老化迁移热力学方程与试验结果基本吻合，迁移驱动因子对数与温度的倒数呈现出较为明显的线性负相关。发射药中NG与DBP总含量高于15%时，发射药的迁移驱动因子（动态活度最大值的百分比增长速率）更易受温度影响从而加剧钝感剂迁移导致钝感失效；而NG和DBP的总含量不高于15%时，钝感剂迁移导致失效的速率较为缓慢，更有利于发射药贮存后的使用和射击安全。

图文摘要

The functional invalidation rate caused by the deterrent migration is faster than the rate of chemical stability invalidation in the aging process of deterred DB propellants. The aging migration model of deterred DB propellants was established， and the relevant parameters were acquired by experiments.

关键词

钝感; 双基发射药; 老化迁移; 扩散原理; 动态活度

1 引言

轻武器是一类应用最广、数量最多的武器种类，其水平是军事强国的特征标志之一^［

1］。目前，世界各国轻武器发射装药广泛采用钝感双基球形（扁）药，该药具有装填密度高、流散性好、弹道优良等特点。但是钝感剂为极易扩散的小分子物质，可能随着武器的使用地域、温差、湿度等不同而发生巨大变化，给武器的使用安全性带来了严重的威胁^{［参考文献 2

百度学术}2］。随着武器信息化、人性化的发展，对武器的烟、焰、残渣、精度、使用安全性、尤其是弹道稳定性等方面提出了更高的要求^{［参考文献 3

百度学术}3］。

化学安定性对发射药的制造、贮存和使用具有十分重要的意义^［

4］。发射药的初始热分解的部分产物，通常对发射药的后续分解具有催化作用，使得发射药热分解反应加速，贮存安定性变差。实际生产中通过添加二苯胺（DPA）、二甲基二苯脲（C₂）等化学安定剂来提高发射药的化学安定性，原理是安定剂能够和发射药的热分解产物发生反应，从而抑制发射药的热分解自催化作用^{［参考文献 5-6}5-6］。

近年来在测试一些贮存老化后钝感双基发射药的性能时发现，虽然发射药基体中的安定剂消耗量低于50%阈值^［

7］，即贮存安定性满足要求，但在射击时出现了严重的超压现象，显著影响了轻武器的精度要求及使用安全性。这种现象产生的根本原因是，钝感发射药为非均质体系，钝感剂在钝感发射药表面呈一定的梯度分布，老化过程中在浓度差和温度驱使下发射药表面的钝感剂逐渐向内部扩散，导致发射药表面阻燃效果下降，起始燃烧活度变大，膛压升高。

钝感剂迁移对发射药的贮存安全性并没有实质影响，但钝感剂在发射药中的分布改变，就会使得发射药燃烧过程异常，在内弹道上表现为膛压增高、初速降低，既影响作战能力、又对武器系统带来安全隐患。因此，研究钝感发射药老化迁移的机理及动力学过程意义重大。

在钝感剂分布以及扩散迁移方面，国内外均进行了大量的理论及实验研究。起初染色法、同位素法、电感耦合等离子发射光谱法等被用于钝感剂的分布及扩散规律研究，获取到了一定的钝感剂分布深度等信息^［

8-10］。后来逐渐发展了显微红外光谱及激光显微共焦拉曼光谱测定钝感剂浓度分布技术，该方法不破坏发射药的原始状态，测试结果准确性大大提高，在非均质发射药微区定量分析中显示了独特优势^{［参考文献 11-14}11-14］。

虽然钝感剂在发射药表面的浓度分布可由多种方式检测获取，但通常轻武器用发射药尺寸小、装药样本量大，这就导致一方面测试误差较大，另一方面少样本量发射药粒的浓度分布变化很难与大样本量发射装药的燃烧性能建立对应关系。密闭爆发器是一种测试发射药静态燃烧性能的超高耐压装置^［

15］，结果稳定、可重复性好，可用于钝感双基发射药老化迁移性质的表征与评价。

本研究从钝感剂在发射药中的扩散机理入手，建立了钝感发射药老化迁移导致钝感剂浓度分布及燃烧性能变化的数学模型。对6种发射药进行不同温度和不同时间的加速老化实验，采用密闭爆发器测试了老化发射药的燃烧性能，分析了发射药老化迁移的动力学过程及热力学影响，以期为钝感发射药迁移评价和钝感发射药迁移抑制研究提供理论基础。

2 钝感双基发射药老化迁移机理

2.1 钝感剂扩散规律

发射药钝感处理是指用不含能或者能量较低的钝感剂渗入发射药，在扩散作用下到达一定的深度，从而使发射药实现渐增燃烧。起始时，钝感剂以较高浓度分布在发射药表面，随着老化进行，在浓度差和温度共同驱动下，钝感剂由发射药表面逐渐向内迁移，这是发射药老化过程中的物理变化层面。综合热分解化学变化层面的发射药安定性问题，钝感发射药老化失效途径及原理如图1所示。

图1 钝感发射药老化失效途径及原理

Fig.1 Aging pathway and principle of deterred propellant

钝感剂在发射药中的扩散规律通常遵循Fickian第二扩散定律^［

16］，即非稳态扩散模型：

\frac{\partial C}{\partial t} = D \frac{\partial^{2} C}{\partial x^{2}}

（1）

式中，x为所在处到表面的距离（深度），μm；t为扩散时间，s；C为钝感剂质量浓度，mg·μm^-3。

对于钝感发射药而言，起初钝感剂分布在发射药颗粒表面，之后逐渐向内扩散，符合典型非稳态扩散的第二类情况。可采用正态分布N（μ，σ）的右支描述，标准正态分布的形式为：

f (x) = \frac{1}{\sqrt[]{2 π} σ} e^{- \frac{{(x - μ)}^{2}}{2 σ^{2}}}

（2）

\int_{0}^{+ \infty} f (x) d x = \frac{1}{2}

（3）

对于钝感发射药而言，μ=0；钝感剂在发射药中的分布记为C_x；保证分布函数积分为钝感剂总含量m，mg。则钝感剂分布方程为：

C_{x} = \frac{2 m}{\sqrt[]{2 π} σ} e^{- \frac{x^{2}}{2 σ^{2}}} (x \geq 0)

（4）

在正态分布中，3σ准则表明数值分布在（μ-3σ，μ+3σ）中的概率为0.9974，可以认为，函数的取值几乎全部集中在（μ-3σ，μ+3σ）区间内。因此，在钝感发射药中，可以近似认为钝感剂分布深度d为半个置信区间，即：

d = 3 σ

（5）

进一步，钝感剂分布函数为：

C_{x} = \frac{6 m}{\sqrt[]{2 π} d} e^{- \frac{{(3 x)}^{2}}{2 d^{2}}} (x \geq 0)

（6）

直接应用Fickian第二定律解决实际问题很困难，但是利用Fickian定律的推论十分简单、有效。其中一条推论为，扩散时间t与扩散距离d的平方成正比^［

16］。假定初始发射药钝感剂分布深度为d₀，记等效初始时刻为t₀，则：

d_{0} = k \sqrt[]{D t_{0}}

（7）

式中，k为相关系数；D为系统的扩散系数，μm²·s^-1。

那么扩散t时间后对应的扩散深度d为：

d = k \sqrt[]{D (t_{0} + t)} = k \sqrt[]{(d_{0} {/ k)}^{2} + D t}

（8）

将扩散距离d与扩散时间t的关系带入钝感剂分布函数，则：

C_{x, t} = \frac{6 m}{\sqrt[]{2 π} d_{0} \sqrt[]{1 + {(k / d_{0})}^{2} D t}} e^{- \frac{{[(3 / d_{0}) x]}^{2}}{2 [1 + (k / d_{0})^{2} D t]}}

（9）

进一步，令：

t^{'} = (k / d_{0})^{2} D t

（10）

x^{'} = (3 / d_{0}) x

（11）

则：

\frac{d_{0}}{6 m} C_{x^{'}, t^{'}} = \frac{1}{\sqrt[]{2 π} \sqrt[]{1 + t^{'}}} e^{- \frac{{x^{'}}^{2}}{2 (1 + t^{'})}}

（12）

由式（12）可知，不同老化时间t′对应的发射药钝感剂浓度分布如图2所示。初始阶段钝感剂主要分布在发射药表层，随着老化进行钝感剂逐渐向内扩散，到达的深度更深，发射药表层的钝感剂浓度下降，钝感剂分布的梯度性质逐渐减弱。

图2 发射药中钝感剂分布随老化时间的变化关系

Fig.2 Relationship between the deterrent distribution and aging time

2.2 钝感剂迁移对燃烧性能的影响

2.2.1 球形药的形状函数

钝感发射药的燃烧性能主要受两方面影响：一是本身的药型尺寸，它决定了发射药燃烧过程中的燃烧面积变化规律（几何燃烧定律）^［

17］；二是钝感剂分布，它决定了发射药燃烧过程中燃速系数的变化规律^{［参考文献 18

百度学术}18］。

发射药的燃面变化一般通过形状函数来描述^［

19］：

Z = \frac{e}{e_{1}}

（13）

ψ = χ Z (1 + λ Z + μ Z^{2})

（14）

式中，Z为发射药已燃厚度百分比；2e为燃烧层厚度，mm；2e₁为弧厚，mm；Ψ为发射药已燃质量百分比；χ、λ、μ是只与发射药药型相关的特性参数。

对于球形药，χ=3，λ=-1，μ=1/3，则球形药的形状函数为：

ψ = 3 Z (1 - Z + \frac{1}{3} Z^{2}) = 1 - {(1 - Z)}^{3}

（15）

\frac{d Ψ}{d Z} = 3 {(1 - Z)}^{2}

（16）

2.2.2 燃速系数与钝感剂浓度的关系

发射药燃烧过程中燃烧面沿法向向内的推进速度称为燃烧速度，内弹道计算广泛采用指数函数式来描述，即：

\dot{r} = u_{1} p^{n}

（17）

式中，p为燃烧压力，MPa；u₁为燃速系数，cm·s^-1·MPa^-1；n为燃速压力指数，通常n取1。

对于钝感发射药而言，不同深度发射药的钝感剂含量不同，因此其能量和燃速均有差异。通常实验条件下发射药的燃速系数随钝感剂浓度呈现出较好的线性下降趋势^［

3］，因此它们的关系可表示为：

u_{1} = u_{1_0} - k^{'} C_{x, t^{'}} > 0

（18）

其中，u_{1_0}为钝感剂含量为0时发射药的燃速系数，k′为燃速系数的浓度影响因子。

2.2.3 球形药燃气生成猛度的变化规律

钝感发射药与均质发射药的区别在于，其燃烧过程受本身的形状函数和燃速系数随钝感剂浓度变化共同影响，一般采用燃气生成猛度Γ函数来表示发射药的燃烧规律性。

以球形药为例，对钝感剂迁移过程中钝感剂分布变化进行理论分析。燃气生成猛度Γ可表示为：

Γ = \frac{1}{p} \frac{d Ψ}{d t} = \frac{1}{p} \frac{d Ψ}{d Z} \frac{d Ζ}{d t}

（19）

其中：

\frac{d Ζ}{d t} = \frac{1}{e_{1}} \frac{d e}{d t}

（20）

\frac{d e}{d t} = \dot{r} = u_{1} p^{n}

（21）

综上（9）、（16）以及（18）~（21）式，Γ可转化为：

\begin{array}{l} Γ = \frac{1}{e_{1}} \frac{d Ψ}{d Z} u_{1} \\ = \frac{1}{e_{1}} 3 {(1 - Z)}^{2} (u_{1_0} - \frac{6 m k^{'}}{d_{0}} \frac{1}{\sqrt[]{2 π} \sqrt[]{1 + t^{'}}} e^{- \frac{{[(3 / d_{0}) x]}^{2}}{2 (1 + t^{'})}}) \end{array}

（22）

此外，位置x也是这里所描述的燃烧层厚度（2e）的一半，即：

Z = \frac{e}{e_{1}} = \frac{x}{e_{1}}

（23）

因此，可以得到Γ与Z的关系式：

\frac{e_{1}}{3 u_{1_0}} Γ_{Z, t^{'}} = {(1 - Z)}^{2} (1 - \frac{6 m k^{'}}{d_{0} u_{1_0}} \frac{1}{\sqrt[]{2 π} \sqrt[]{1 + t^{'}}} e^{- \frac{9 {(e_{1} / d_{0})}^{2} Z^{2}}{2 (1 + t^{'})}})

（24）

分别令：

p = \frac{6 m k^{'}}{d_{0} u_{1_0}}

（25）

q = \frac{e_{1}}{d_{0}}

（26）

定义p为钝感强度，它代表了钝感剂对于发射药基体燃速系数的综合抑制强度；q则是初始钝感相对深度的倒数。

同时，将球形药Ψ‑Z关系式（15）带入式（24），于是Γ与Ψ的关系可最终表示为：

\frac{e_{1}}{3 u_{1_0}} Γ_{ψ, t^{'}} = {(1 - ψ)}^{\frac{2}{3}} (1 - p \frac{1}{\sqrt[]{2 π} \sqrt[]{1 + t^{'}}} e^{- \frac{9 q^{2} {(1 - \sqrt[3]{1 - ψ})}^{2}}{2 (1 + t^{'})}})

（27）

基于以上物理模型和数学推导，取不同钝感条件（p，q）参数，绘制了发射药的燃气生成猛度曲线随老化时间的变化关系，如图3所示。由图3a、图3b、图3c可知，随着p增大，钝感剂对发射药基体燃速系数的综合抑制强度增加，发射药燃烧起始动态活度更低；由图3d、图3e、图3f可知，随着q增大（钝感深度越浅），发射药达到动态活度最大值对应的已燃质量百分比越靠前。此外，钝感发射药老化前后燃气生成猛度曲线的增减性分两种情况：一种是老化迁移程度较低时，发射药燃气生成猛度曲线先微弱降低、再急剧升高、最后逐渐下降；第二种是老化迁移程度较高时，发射药燃气生成猛度曲线单调递减。随老化进行，发射药初期燃气生成猛度逐渐上升，燃气生成猛度抑制作用范围后延，这是发射药表层钝感剂向内迁移的原因和作用结果。

a. p=1.6， q=3.0

b. p=2.0， q=3.0

c. p=2.4， q=3.0

d. p=2.0， q=2.0

e. p=2.0， q=4.0

f. p=2.0， q=6.0

图3 不同钝感条件发射药燃气生成猛度曲线随老化时间的变化关系

Fig.3 Relationship between the dynamic vivacity curves and aging time of propellant under different deterred conditions

钝感发射药具有典型的不同时点火特性，且钝感发射药实际燃烧过程中燃气生成猛度是大样本发射药粒不同时点火的共同影响结果。为此，可以按式（28）对发射药的Γ进行Ψ在0~0.5范围内的积分，记为I。并且按式（29）计算燃气生成猛度积分值的增长率R_I：

I (t^{'}) = \int_{0}^{0.5} \frac{e_{1}}{3 u_{1_0}} Γ_{ψ, t^{'}} d ψ = \frac{e_{1}}{3 u_{1_0}} \int_{0}^{0.5} Γ_{ψ, t^{'}} d ψ

（28）

R_{I} (t^{'}) = \frac{I_{t^{'}} - I_{0}}{I_{0}} \times 100 %

（29）

钝感条件（p=2.0，q=3）发射药燃气生成猛度曲线积分值及其增长率随老化时间t′的变化关系如图4所示。图4可以看出，燃气生成猛度积分值的增长率随老化时间呈现出较好的线性关系。

图4 发射药燃气生成猛度积分随老化时间的线性增长

Fig.4 Linear increasing of the integral of dynamic vivacity of propellants with aging time

综合以上钝感剂扩散规律、钝感剂浓度与发射药燃速系数关系、发射药形状函数以及燃气生成猛度的变化规律，可以看出，老化过程中钝感剂扩散迁移对发射药的燃烧性质产生了显著的规律性的影响。随着老化进行，钝感剂向内迁移，钝感剂浓度梯度分布性质减弱，初始燃烧阶段发射药燃速变快，燃气生成猛度升高。同样体现在发射装药的内弹道试验中，钝感剂的迁移导致发射药初始燃烧过快，发射最大膛压升高，初速不稳定，出现“弹道漂移”现象，严重影响武器使用安全。

基于发射药燃气生成猛度积分值的增长率随老化时间呈现出较好的线性关系，可以采用发射药密闭爆发器试验所得发射药的动态活度随老化时间的变化来评价钝感发射药的老化迁移程度，从理论上是可行的，同时能够避免直接进行射击试验可能出现的危险。

2.3 老化迁移模型及特征参量

2.3.1 动力学方程

发射药的动态活度通常用L表示，进一步记每组老化样品不同老化天数t_aging的动态活度升高百分比为R_L（t_aging），则钝感双基发射药老化迁移的动力学方程可表述为：

R_{L} (t_{a g i n g}) = k_{d} \cdot t_{a g i n g}

（30）

其中，k_d为R_L（t_aging）-t_aging直线的斜率，定义为迁移驱动因子。

迁移驱动因子k_d代表发射药老化过程中动态活度最大值的百分比增长速率，受环境温度、发射药配方、残余溶剂含量以及钝感剂种类等多种因素共同影响。

2.3.2 热力学方程

温度是影响钝感剂在发射药中扩散速率的关键因素，表现在迁移驱动因子k_d的差异。迁移驱动因子k_d与老化环境温度T的关系可表述为：

k_{d} = c_{d} e^{- l_{d} \cdot \frac{1}{T}}

（31）

式中，c_d为驱动系数；l_d为驱动阈能，J·g^-1。

采用对数法，可将老化迁移的热力学方程转化为：

l n (k_{d}) = l n (c_{d}) - l_{d} \cdot \frac{1}{T}

（32）

3 钝感发射药老化实验及燃烧性能测试

3.1 加速老化实验条件

实验选取了不同硝化甘油（NG）含量、不同钝感剂邻苯二甲酸二丁酯（DBP）用量的6种发射药，在55、65、75 ℃和85 ℃多个温度下进行0~50 d不等的加速老化试验，具体见表1。

表1 实验样品及老化条件

Table 1 Experimental samples and aging conditions

number	NG / %	DBP / %	T / ℃	t / d
1	5	5	65， 75	0~50
2	10	3.5	55， 75， 85
3	10	5	55， 75， 85
4	10	6.5	55， 65， 75， 85
5	15	5	55， 65， 75， 85
6	20	5	55， 65， 75， 85

每种样品在每个温度点和每个老化时间的待老化试样称取（100±0.1） g，分别装入9 cm×14 cm长方形的铝箔/PE复合密封袋中，如图5所示。再放入相应温度的水浴烘箱进行加速老化实验。

图5 装入铝箔/PE复合密封袋的发射药老化实验样品

Fig.5 Aging experimental samples of propellants in aluminum foil / PE composite sealing bags

3.2 密闭爆发器实验

采用50 mL密闭爆发器研究老化前后发射药的燃烧性能。发射药装填密度为0.2 g·mL^-1，装药量以实际标定爆发器容积乘以装填密度来计算，点火药为0.55 g C型硝化棉。密闭爆发器实验在常温下进行，即发射药初温为常温。

实验中通过压力传感器记录密闭爆发器内腔压力（p）随发射药燃烧时间（t）的关系，再计算得到发射药燃烧过程中动态活度的变化：

L (t) = \frac{d p (t) / d t}{p (t) \cdot p_{m}}

（33）

进一步，取动态活度曲线（L‑t）上的最大值，记为L_m：

L_{m} = M a x \{L (t)\}

（34）

对于同一个未老化发射药样品，记老化t_aging时间后动态活度最大值为L_m（t_aging），未老化样品的动态活度最大值为L_m（t₀）。老化t时间后发射药动态活度最大值的升高百分比，记为R_L（t_aging），即：

R_{L} (t_{a g i n g}) = \frac{L_{m} (t_{a g i n g}) - L_{m} (t_{0})}{L_{m} (t_{0})} \times 100 %

（35）

4 结果与讨论

4.1 迁移失效的动力学过程

为研究钝感发射药老化过程中迁移失效的动力学过程，分析了前述6种发射药加速老化后的动态活度最大值升高百分比R_L（t_aging），随老化时间t_aging的关系，如图6所示。可以看出，钝感发射药老化过程中R_L（t_aging）随t_aging呈较强的线性关系，这与前述的理论研究结果是一致，说明该理论模型构建合理、能够清楚地表述钝感发射药老化过程中燃烧性能的变化规律。

a. 5NG， 5DBP

b. 10NG， 3.5DBP

c. 10NG， 5DBP

d. 10NG， 6.5DBP

e. 15NG， 5DBP

f. 20NG， 5DBP

图6 动态活度最大值升高百分比与老化时间的关系

Fig.6 Relationship between R_L（t_aging） and ageing time

进一步，对每组老化样品不同老化天数t_aging的动态活度升高百分比R_L（t_aging）进行截距为零的线性拟合，6种发射药不同温度下拟合方程的迁移驱动因子k_d及拟合度r²如表2中所示。

表2 6种发射药不同温度老化的拟合参数

Table 2 Fitting parameters of the six propellants at different temperature

No.	NG / %	DBP / %	T / ℃	k_d / %·d^-1	r²
1	5	5	65	0.1103	0.9784
2	5	5	75	0.2681	0.9949
3	10	3.5	55	0.4792	0.9941
4			75	0.6644	0.9872
5			85	0.9538	0.9983
6	10	5	55	0.3129	0.9896
7			75	0.7200	1.0000
8			85	0.7655	0.9466
9	10	6.5	55	0.1617	0.9946
10			65	0.4656	0.9989
11			75	1.4650	1.0000
12			85	2.1234	0.973
13	15	5	55	0.5823	0.9951
14			65	0.8056	0.8701
15			75	4.6978	0.9988
16			85	4.7266	0.8873
17	20	5	55	1.8673	0.9837
18			65	6.5239	0.9332
19			75	7.7533	0.8981
20			85	11.7920	0.9535

可以看出，各组老化样品的拟合度基本在0.98以上，说明钝感发射药的动态活度升高百分比R_L（t_aging）与老化时间t_aging线性度较高。进一步分析可知，k_d受NG含量、DBP用量以及老化温度的影响较显著。总体上看NG与DBP含量越高、老化温度越高，k_d越大。75 ℃下，NG含量由5%提高到20%时，k_d提高27.9倍；85 ℃下，DBP用量由3.5%提高到6.5%时，k_d提高1.2倍。

4.2 迁移失效的热力学影响

前述6种发射药k_d随T的关系如图7所示，并采用式（32）钝感发射药老化迁移热力学方程进行拟合，得到的驱动系数c_d、驱动阈能l_d以及拟合度r²如表3所示。

a. 5NG， 5DBP

b. 10NG， 3.5DBP

c. 10NG， 5DBP

d. 10NG， 6.5DBP

e. 15NG， 5DBP

f. 20NG， 5DBP

图7 6种发射药的老化驱动因子随温度变化的拟合曲线

Fig. 7 Fitting curves of the aging driving factors of six propellants with the change of temperature

表3 6种发射药老化迁移热力学方程拟合参数

Table 3 Fitting parameters of thermodynamic equation for the aging migration of six propellants

samples	c_d	l_d	r²
5NG 5DBP	2.961×10¹²	1.046×10⁴	1
10NG 3.5DBP	1.144×10³	2.562×10³	0.9432
10NG 5DBP	2.648×10⁴	3.709×10³	0.9384
10NG 6.5DBP	1.271×10¹³	1.047×10⁴	0.9700
15NG 5DBP	1.790×10¹²	9.470×10³	0.8609
20NG 5DBP	2.090×10⁹	6.755×10³	0.8761

从图7和表3可以看出，上述k_d与T的热力学方程基本符合实验结果，ln（k_d）与1/T呈现出较为明显的线性负相关性，理论模型假设合理。此外，当NG和DBP总含量超过20%时，老化迁移的拟合度降低到0.9以下，这表明当发射药体系内小分子总量高时，发射药体系的扩散迁移及其受温度影响的规律性下降，这与前述k_d随NG和DBP含量增大的规律是一致的。

4.3 增塑剂与钝感剂总含量的影响

将6种发射药的拟合特性参数代入老化迁移热力学方程中，绘制了6种发射药的k_d随T的变化关系，如图8所示。可以明显看出，温度对发射药老化迁移驱动因子具有显著影响，控制贮存温度有利于减缓发射药的老化迁移失效速率。此外，当NG和DBP的总含量高于15%时，发射药k_d受温度影响增长迅速，高温下更容易因钝感剂迁移而失效；当NG和DBP的总含量不高于15%时，发射药k_d受温度影响增长较为缓慢，钝感剂迁移导致的失效速率增长也较为缓慢，更有利于贮存后的使用和射击安全。

图8 老化驱动因子随温度的变化关系

Fig.8 Relationship between the aging driving factor of propellants and temperature

4.4 常温环境迁移失效规律

根据迁移失效热力学方程，可以外推计算出常温下发射药的老化迁移驱动因子。此外假定动态活度升高百分比R_L（t_aging）=8.64%作为比较依据，则发射药的迁移失效时间如式（36）：

t_{i n v a l i d a t i o n} = \frac{8.64}{k_{d}} = \frac{8.64}{c_{d} e^{- l_{d} \cdot \frac{1}{T}}}

（36）

绘制了前述3种相同DBP用量、不同NG含量发射药的迁移失效时间与老化温度的变化关系，并且计算了常温下3种发射药的有效使用寿命，如图9中所示。从图中可以看出，双基钝感发射药受迁移影响的贮存寿命，在钝感剂含量相同时，随NG含量增加发射药的贮存年限大大缩短，这是因为NG增加了发射药体系的塑性，增加了DBP在发射药基体内的扩散速率；双基钝感发射药受迁移影响的贮存寿命，随老化环境的温度升高而大幅降低；经估算可得：（5NG 5DBP）样品受迁移影响的常温有效使用寿命为13.8年，（15NG 5DBP）样品为10.1月，（20NG 5DBP）样品为28.8天。

图9 3种不同NG含量发射药有效使用寿命随老化温度的变化

Fig.9 Chang of effective service life of three NG content propellants with aging temperature

5 结论

利用钝感剂的正态分布特性及Fickian第二扩散定律的推论，建立了钝感发射药老化迁移理论模型，分析了钝感剂迁移导致钝感剂浓度分布以及燃烧性能变化的影响规律。对6种不同配方的发射药进行了加速老化实验，采用密闭爆发器研究了老化迁移发射药的燃烧性能，分析了发射药老化迁移的动力学过程及热力学影响，主要结论如下：

（1）钝感球形发射药燃气生成猛度积分值的增长率随老化时间呈现出较好的线性关系，从理论上证明了大样本发射药粒燃烧性能随老化迁移的影响过程。

（2）6种发射药加速老化后的动态活度最大值升高百分比随老化时间均呈较强的线性关系，与建立的老化迁移动力学模型一致；发射药中NG与DBP含量越高、老化温度越高，迁移驱动因子k_d越大，相同时间发射药老化迁移程度越大，动态活度升高百分比越大。

（3）建立的钝感发射药老化迁移热力学方程与试验结果基本吻合，ln（k_d）与1/T呈现出较为明显的线性负相关性。

（4）发射药中NG和DBP的总含量高于15%时，发射药的迁移驱动因子k_d更易受温度影响，也就是高温下更容易迁移失效；而NG和DBP的总含量不高于15%时，钝感剂迁移导致失效的速率较为缓慢，更有利于发射药贮存后的使用和射击安全。

参考文献

肖忠良. 火炸药导论［M］. 北京：国防工业出版社， 2019. [百度学术]

XIAO Zhong‑liang. Introductory theory of propellants and explosives［M］. Beijing： National Defense Industry Press， 2019. [百度学术]

Manning T， Wyckoff J， Klingaman K， et al. Enhanced propellant performance via environmentally friendly curable surface coating［J］. Defence Technology， 2017， 13（3）： 131-142. [百度学术]

Li S， Tao Z， Ding Y， et al. Gradient denitration strategy eliminates phthalates associated potential hazards during gun propellant production and application［J］. Propellants， Explosives， Pyrotechnics， 2020 ，45（7）： 1156-1167. [百度学术]

Zhao X， Rui X， Wang Y， et al. Propellant shelf‑life prediction under temperature and relative humidity conditions based on DPA kinetics［J］. Journal of Energetic Materials，2019，37（4）： 407-419. [百度学术]

宋力骞，刘大斌，钱华，等. 某三基发射药的老化寿命评估［J］. 火炸药学报， 2018， 41（6）： 105-109. [百度学术]

SONG Li‑qian， LIU Da‑bin， QIAN Hua， et al. Aging life evaluation of a tri‑base gun propellant［J］. Chinese Journal of Explosives and Propellants， 2018， 41（6）： 105-109. [百度学术]

周静，丁黎，张俊林，等. 基于近红外技术的火药安定期评价新方法［J］. 兵工学报， 2018， 39（3）： 463-467. [百度学术]

ZHOU Jing， DING Li， ZHANG Jun‑lin， et al. New approach for evaluating the stability period of propellant based on near‑infrared spectroscopy［J］. Acta Armamentarii， 2018， 39（3）： 463-467. [百度学术]

国防科学技术工业委员会. GJB 770B-2005火药试验方法［S］. 北京：国防科工委军标出版发行部， 2005. [百度学术]

COSTIND. GJB 770B-2005 Test method of propellant［S］. Beijing： Publication Department of COSTIND， 2005. [百度学术]

应三九，罗运军，徐复铭. 无机钝感剂在深钝感发射药中的应用研究［J］. 兵工学报， 1995， 1（1）： 10-14. [百度学术]

YING San‑jiu， LUO Yun‑jun， XU Fu‑ming. Study on the application of inorganic deterrent agent in deep deterring propellant［J］. Acta Armamentarii， 1995， 1（1）： 10-14. [百度学术]

Fong C W， Cooke C. Diffusion of deterrents into a nitrocellulose‐based small arms propellant. The effect of deterrent structure and solvent［J］. Journal of Applied Polymer Science， 1982， 27（8）： 2827-2832. [百度学术]

Brodman BW，Devine MP，Finch RW， et al. Autoradiographic determination of the di‑n‑butyl phthalate concentration profile in a nitrocellulose matrix［J］. Journal of Applied Polymer Science， 1974， 18（12）： 3739-3744. [百度学术]

Vogelsanger B， Ossola B， Brönnimann E. The diffusion of deterrents into propellants observed by FTIR microspectroscopy‐quantification of the diffusion process［J］. Propellants， explosives， pyrotechnics， 1996 ，21（6）： 330-336. [百度学术]

Louden J D， Kelly J， Phillipson J. Raman microspectroscopic determination of the methyl centralite （N，N'‐dimethyl‐N，N'‐diphenylurea） concentration profile in a nitrocellulose extruded mono‐perforated small arms propellant［J］. Journal of Raman spectroscopy， 1987， 18（2）： 137-140. [百度学术]

Trewartha S， Shapter J， Gibson C T， et al. Determination of deterrent profiles in nitrocellulose propellant grains using confocal raman microscopy［J］. Propellants， Explosives， Pyrotechnics， 2011， 36（5）： 451-458. [百度学术]

Roduit B， Brandsch R， Folly P， et al. Simulation of deterrent diffusion in double base propellant under different temperature profiles［C］//International Annual Conference‑Fraunhofer‑Institut fur Chemische Technologie ICT， 2010. [百度学术]

赵小亮，张小兵. 管状发射药内孔侵蚀燃烧与流动特性［J］. 含能材料， 2019， 27（3）： 202-209. [百度学术]

ZHAO Xiao‑liang， ZHANG Xiao‑bing. Internal perforation erosive burning and flow characteristics of tubular propellant［J］， Chinese Journal of Energetic Materials（Hanneng Cailiao）， 2019， 27（3）： 202-209. [百度学术]

沈青. 高分子物理化学［M］. 北京：科学出版社， 2016. [百度学术]

SHEN Qing. Polymer physical chemistry Ⅰ［M］. Beijing： Science Press， 2016. [百度学术]

王泽山，何卫东，徐复铭. 火炮发射装药设计原理与技术［M］. 北京：北京理工大学出版社， 2014. [百度学术]

WANG Ze‑shan， HE We‑dong， XU Fu‑ming. Principle and Technique for Gun Propellant Charge Design［M］. Beijing： Beijing Institute of Technology Press， 2014. [百度学术]

易文俊，汤瑞峰，邱佩蓉. 钝感发射药的内弹道特性研究［J］. 南京理工大学学报， 1997，21（02）：133‑136. [百度学术]

YI Wen‑jun， TANG Rui‑feng， QIU Pei‑rong. A Study on interior ballistic characteristics of passivation propellant［J］. Journal of Nanjing University of Science and Technology， 1997，21（2）： 133-136. [百度学术]

张小兵，金志明. 枪炮内弹道学［M］. 北京：北京理工大学出版社， 2014. [百度学术]

ZHANG Xiao‑bing，JIN Zhi‑ming. Interior Ballistics of Guns［M］. Beijing： Beijing Institute of Technology Press， 2014. [百度学术]

首页

期刊简介

编委会

作者专区

审稿专区

道德声明

联系我们

进入EMF

钝感双基发射药老化迁移机理及动力学过程 PDF

摘要

图文摘要

关键词

1 引言

2 钝感双基发射药老化迁移机理

2.1 钝感剂扩散规律

2.2 钝感剂迁移对燃烧性能的影响

2.3 老化迁移模型及特征参量

3 钝感发射药老化实验及燃烧性能测试

3.1 加速老化实验条件

3.2 密闭爆发器实验

4 结果与讨论

4.1 迁移失效的动力学过程

4.2 迁移失效的热力学影响

4.3 增塑剂与钝感剂总含量的影响

4.4 常温环境迁移失效规律

5 结论

参考文献

钝感双基发射药老化迁移机理及动力学过程 PDF

摘要

图文摘要

关键词

1 引 言

2 钝感双基发射药老化迁移机理

2.1 钝感剂扩散规律

2.2 钝感剂迁移对燃烧性能的影响

2.3 老化迁移模型及特征参量

3 钝感发射药老化实验及燃烧性能测试

3.1 加速老化实验条件

3.2 密闭爆发器实验

4 结果与讨论

4.1 迁移失效的动力学过程

4.2 迁移失效的热力学影响

4.3 增塑剂与钝感剂总含量的影响

4.4 常温环境迁移失效规律

5 结 论

参考文献

1 引言

5 结论