阅读说明：本技术 一种渗透工艺球扁发射药的燃速计算方法 (Burning rate calculation method for penetration process ball flat propellant powder ) 是由 张江波 梁泰鑫 赵树森 肖飞 刘威 刘幼平 肖忠良 于 2021-08-29 设计创作，主要内容包括：本发明公开了一种渗透工艺球扁发射药的燃速计算方法,适用于渗透工艺球扁发射药燃速的计算。本发明采用将液体组分渗透到固体中的不同工艺方法统一采用渗透分布函数进行表征,去除了工艺参数的影响,能够实现任意多种渗透组分发射药的燃烧特性表征；通过采用渗透分布函数结合均质发射药组分组合成新配方的方式直接进行能量计算,能够准确地获得球扁发射药不同分布位置处的火药力；采用将渗透层和均质层分为两种药形函数的计算方法,能够较好的描述非完全平行层燃烧的物理过程。本发明解决了渗透工艺球扁发射药的燃速精确计算问题,提供了一种渗透工艺球扁发射药燃速计算的技术手段。(The invention discloses a method for calculating the burning rate of a penetration process ball flat propellant powder, which is suitable for calculating the burning rate of the penetration process ball flat propellant powder. The method adopts different process methods for permeating the liquid components into the solid to uniformly characterize by adopting the permeation distribution function, thereby removing the influence of process parameters and realizing the combustion characteristic characterization of any multiple permeation component propellant powder; energy calculation is directly carried out by combining a penetration distribution function with a homogeneous propellant powder to form a new formula, so that the propellant powder at different distribution positions of the flat propellant powder can be accurately obtained; by adopting a calculation method of dividing the permeable layer and the homogeneous layer into two types of medicine shape functions, the physical process of combustion of the incomplete parallel layer can be better described. The invention solves the problem of accurate calculation of the burning rate of the flat propellant powder of the infiltration process ball and provides a technical means for calculating the burning rate of the flat propellant powder of the infiltration process ball.)

一种渗透工艺球扁发射药的燃速计算方法

技术领域

本发明涉及发射药技术领域，尤其涉及发射药燃烧速度的计算，具体涉及一种针对采用渗透工艺制备的球扁发射药，能够实现其燃速的计算。

背景技术

球扁发射药是具有椭圆体形状的一种发射药，主要应用于小口径武器或者需要低膛压速燃装药的武器。为了降低球扁发射药的起始燃烧压力、增加其能量或燃烧渐增性，通常需要对球扁发射药进行后处理，一般采用深钝感或浸渍高能组分的方法，这些方法的共同特点就是采用液态功能试剂对成型后的球扁发射药进行渗透，而后进行烘干等后续的工艺，采用这种渗透工艺的球扁发射药具有能量组分分布不均匀，功能试剂从球扁发射药表层到内部按一定的规律分布，通过对功能试剂品种及分布规律的调整可以有效地对球扁发射药的燃烧性能进行调控，进而提高武器的技术性能。其中，渗透工艺球扁发射药的燃烧性能决定了其装药应用特性，因此，开展其燃速的研究具有重大意义。现在的渗透工艺球扁发射药的燃速计算方法主要有：(1)将其不同组分物质的特性进行平均，将渗透工艺球扁发射药作为宏观的均一物质进行处理，这明显与实际物质的组成不符，无法得到较准确的燃烧过程参数；(2)对每一组分采用与渗透深度相关的近似函数表示，总组分的参数认为是每一组分参数的加和，该方法尽管考虑了渗透物质的分布，但是采用参数加和的方法无法得到较可靠的结果，如不同组分的氧平衡不同，且火药力也不同，将两种一正一负氧平衡的组分混合后，采用火药力加和方法计算两种单独组分获得的火药力比实际火药力小，这是因为在实际燃烧过程中，两种物质将同时参加燃烧，氧将得到有效利用，并不会出现一种组分燃烧产物为正、另一种组分燃烧产物为负的情况；(3)将发射药的药形参数视作一个不变的参数，从开始燃烧到燃烧完毕均符合均值物质燃烧的平行层规律，这显然与实际情况不同，渗透工艺球扁发射药随着弧厚的变化，组分和燃速均在变化。因此，现有计算方法无法较好的获得渗透工艺球扁发射药的燃速。

从目前检索的技术资料中，尚未见具有能够考虑组分渗透深度分布与发射药参数性能及药形变化的渗透工艺球扁发射药的燃速计算方法。

发明内容

为了解决现有技术存在的问题，本发明提供一种考虑组分渗透深度分布与发射药参数性能及药形变化的渗透工艺球扁发射药的燃速计算方法。

本发明提供的计算方法，步骤如下：

(1)取渗透工艺球扁发射药厚度的一半，坐标轴为由外层向内层为x的正方向，渗透组分的特性与x的关系为

H₁(x)＝f₁(x)

H₂(x)＝f₂(x)

x为弧厚方向上的坐标变量，由外层向内层为x正向，H₁(x)为第一种渗透组分的含量，f₁(x)为第一种渗透组分的分布函数，H₂(x)为第二种渗透组分的含量，f₂(X)为第二种渗透组分的分布函数；

(2)不同位置处火药力的计算

f(x)＝power(H，H₁(x)，H₂(x))

f(x)为不同位置处的火药力，power()为火药力计算函数，根据组分的含量进行计算，H为均质基体发射药的各组分含量，即未渗透前的发射药配方；

(3)发射药燃烧的已燃百分数计算

ψ(x)＝χZ(x)(1+λZ(x))

χ＝2+β₁ x≤x₀

λ＝-(1+2β₁)/(2+β₁)x≤x₀

α₁＝e₁/b，β₁＝e₁/c

χ＝2+2β₂ x＞x₀

α₂＝(e₁-x₀)/(b-x₀)，β₂＝(e₁-x₀)/(c-x₀)

x₁₀＝(f₁(x)＝0，x)

x₂₀＝(f₂(x)＝0，x)

x₀＝max(x₁₀，x₂₀)

ψ(x)为不同位置处的已燃百分数，Z(x)为不同位置处的已燃相对厚度，χ为第一项药形系数，λ为第二项药形系数，e₁为发射药弧厚的一半，b为球扁发射药宽度的一半，c为球扁发射药长度的一半，α₁为药形变化前第一项中间参数，β₁为药形变化前第二项中间参数，α₂为药形变化后第一项中间参数，β₂为药形变化后第二项中间参数，x₁₀为第一种渗透组分含量为零的位置，x₂₀为第二种渗透组分含量为零的位置，x₀为发射药渗透层与均质层的边界位置，max()为取两者数值中最大值的函数；

(4)由密闭爆发器实验所测数据计算发射药已燃百分数

Δ为装填密度，ρ为发射药的密度，p′为所测得的燃气压力，p_B为密闭爆发器中的点火压力，α_f为发射药的余容；

(5)计算发射药的已燃相对厚度，联立第(3)步和第(4)步，求出发射药的已燃相对厚度Z(x)与p′的关系

Z(x)＝f_zp(p′)

f_zp()为函数关系表达式；

(6)计算发射药的燃速

t为所测得的对应于p′的时间值，结合第(5)步，即可计算得到渗透工艺球扁发射药的燃速；

(7)重复第(2)到第(6)步，计算不同压力下渗透工艺球扁发射药的燃速。

因此，基于渗透组分的实际分布规律，采用渗透分布函数与发射药特性参数及药形函数紧密关联的计算方法，能够获得较符合实际燃烧过程的渗透工艺球扁发射药燃烧特性数据，较好的对渗透工艺球扁发射药的燃速进行表征。

本发明提供的一种渗透工艺球扁发射药的燃速计算方法，与现有技术相比，有益效果体现在以下几个方面：

(1)将液体组分渗透到固体中的不同工艺方法统一采用渗透分布函数进行表征，去除了工艺参数的影响，能够实现任意多种渗透组分发射药的燃烧特性表征，并且工艺复杂度和组分的多少不增加计算的复杂度。与现有技术相比，实现了各类渗透工艺的统一化，如浸渍工艺、钝感工艺均可用渗透分布函数表征，从而拓展了计算范围，降低了计算难度。

(2)通过采用渗透分布函数结合均质发射药组分组合成新配方的方式直接进行能量计算，能够准确地获得球扁发射药不同分布位置处的火药力，更好的与发射药物理组成相符，较准确的表达发射药的燃烧性能参数。与现有技术相比，实现了与渗透分布函数对应的发射药能量参数的求解，为燃速的准确计算提供了技术支撑，填补了现有技术空白。

(3)采用将渗透层和均质层分为两种药形函数的计算方法，能够较好的描述非完全平行层燃烧的物理过程，进一步提高了计算的精度。与现有技术相比，适用范围更广，能够应用于各类具有非均质结构的球扁发射药。

因此，本发明解决了渗透工艺球扁发射药的燃速精确计算问题，提供了一种非均质球扁发射药燃速计算的技术手段。

附图说明

图1是渗透工艺球扁发射药的示意图。

图2是图1的A向剖视图。

图3是渗透工艺球扁发射药燃速计算流程图。

具体实施方式

下面结合附图对本发明的具体实施方式作进一步的详细说明。

如图1、图2所示，一种渗透工艺球扁发射药1的整体结构呈现椭圆形，短轴长6mm，长轴长10mm，厚度0.5mm；渗透层2的厚度为0.15mm，第一种渗透物质为邻苯二甲酸二丁酯，第二种渗透物质为聚己二酸丙二醇酯，均质层3为含氮量12.6％的硝化棉，厚度为0.2mm。

根据图3所示的渗透工艺球扁发射药燃速计算步骤如下：

(1)取渗透工艺球扁发射药的弧厚的一半(0.25mm)，坐标轴为由外层向内层为x的正方向，渗透组分的特性与x的关系为

H₁(x)＝0.06-7×(x-0.05)²

H₂(x)＝0.05-6.8×(x-0.04)²

(2)不同位置处火药力的计算，计算第一个点，x＝0处的火药力

f(0)＝power(H，H₁(0)，H₂(0))

求得火药力为967.25kJ/kg；

(3)发射药燃烧的已燃百分数计算

x₀＝max(0.15，0.12)＝0.15

α₁＝0.083，β₁＝0.05

χ＝2.05 x≤0.15

λ＝-0.54 x≤0.15

α₂＝0.035，β₂＝0.021

χ＝2.042 x＞0.15

λ＝-0.041 x＞0.15

ψ(x)＝χZ(x)(1+λZ(x))

在不同位置处计算得到了ψ(x)与Z(x)之间的关系；

(4)由密闭爆发器实验所测数据计算发射药已燃百分数，取密闭爆发器中30MPa的数据作为第一个计算数据点，装填密度为0.2g/cm³，发射药密度为1.6g/cm³，点火压力取10MPa，发射药余容为0.001m³，x取厚度的10％即0.025，根据

求得发射药的已燃百分数为0.09；

(5)计算发射药的已燃相对厚度，由第(3)步得到

0.09＝2.05×Z(0.025)×(1-0.54×Z(0.025))

求得Z(0.025)＝0.045；

(6)计算发射药的燃速，根据

重复第(2)到第(6)步对数据点逐次求解，即可获得发射药的燃速。计算获得的特定压力点处渗透工艺球扁发射药燃速如下：

渗透工艺球扁发射药的燃速

序号	压力(MPa)	燃速(mm/s)
			1	30	28.34
2	60	56.50
			3	90	82.79
4	120	105.67
			5	150	128.22
6	180	148.25