榴弹机枪在现代战争中发展前景广阔，高初速、高精度、威力大、质量轻是未来榴弹机枪的发展趋势^[1-2].本文研究的榴弹机枪采用双药室喷管气流反推低后坐发射技术，具有弹丸初速高、威力大、后坐力小的特点.文献[3-4]对采用软后坐技术的大口径枪械进行研究，较常规发射技术，软后坐技术大幅降低了武器的后坐力，但为了得到高效减后坐效率带来首发预压困难、增加武器质量的问题.文献[5-6]提出直接在膛内开孔利用火药燃气减后坐技术降低武器后坐力，但文献[5-6]提出的单药室喷管减后坐技术存在下述问题：获得高效的减后坐效率，喷管导气孔开孔位置必须靠近膛底，这样导致弹丸初速衰减严重；若取喷管导气孔开孔位置靠近膛口，虽然弹丸初速小幅下降，但喷管武器减后坐效果不佳.文献[7-8]建立膨胀波枪炮发射过程的动力学模型，膨胀波减后坐技术有较高的减后坐效率、对弹丸初速影响较小，但结构复杂且膛尾打开时间对弹丸初速较敏感.文献[9]利用制退器降低大口径机枪后坐力，但获得高膛口制退效率必须具有较大装药量，本文研究的榴弹机枪装药量少，无法获得高效的制退效率.

为了克服上述研究不足，提出双药室喷管减后坐技术并应用于高初速榴弹发射器.为了进一步了解火药气流特征量在膛内和喷管内分布状况，建立双药室喷管数学模型.通过调整双药室结构参数，研究其对武器发射性能的影响规律，合理利用火药气体冲量，保证膛口初速几乎不变的前提下使武器获得高减后坐效率，降低武器后坐力.

1 双药室喷管减后坐计算模型 1.1 发射原理

基于双药室喷管减后坐技术的高初速榴弹机枪发射原理如图 1所示.喷管为一个可伸缩管，作战时将喷管拉伸至射手身后，携行时将喷管缩回.

减后坐原理：后药室发射药被点燃，当后药室气体压力达到破膜压力时，一方面火药燃气通过传火孔进入前药室，点燃前药室发射药，另一方面火药燃气通过喷管导气孔进入喷管，经收敛-扩张管道加速后高速冲出喷管尾部，产生喷管反推力.

1.2 气体数学模型及基本假设

图 2为双药室喷管气流数学计算模型，图中：p₁、T₁、ρ₁分别为后药室气体平均压力、平均温度、平均密度；p₂、T₂、ρ₂分别为收敛管道内气体压力、温度、密度；p₃、T₃、ρ₃分别为膛内气体平均压力、平均温度、平均密度；p₄、T₄、ρ₄分别为导气室内气体平均压力、平均温度、平均密度；p_d、T_d、ρ_d分别为膛内导气孔处气体压力、温度、密度；p_D、T_D、ρ_D分别为前药室底传火孔处气体压力、温度、密度；p_A、T_A、ρ_A分别为喷管气室内气体平均压力、平均温度、平均密度.为了研究喷管内气体特征参数在不同位置处变化规律，截取B~I这8个位置.气体数学模型包括后药室燃烧模型、前药室燃烧的内弹道模型、后药室喷管气流计算模型、导气室计算模型、气室之间流量方程.

基本假设：

1) 火药燃烧规律按照皮奥伯特定律^[10]，将实际燃烧过程进行简化，发射过程火药已燃尽，没有固相流动；

2) 内弹道采用修正经典内弹道模型求解，通过次要功系数简化弹丸旋转运动功、弹丸沿膛线运动摩擦功、火药气体运动功等复杂计算；

3) 研究后效期时膛内气流以准定常流处理，膛内任意点气流密度随时间的变化率与该点的密度成正比，膛内气流平均参数按膛轴中点选定；

4) 对于导气装置中气流参数沿空间变化不大的情况，假定气流参数只随时间变化，将导气室内气流参数在空间取平均值^[11]；

5) 由于喷管装置中具有较长的导气管，管内气体在管壁的摩擦、散热影响下，在管内各种复杂波系的作用下，气流参数沿管轴产生较大差别，喷管装置计算采用准一维非定常数学模型^[12].

1.2.1 后药室燃烧模型

后药室燃烧模型为

$ \left\{ \begin{array}{l} \frac{{{\rm{d}}{Z_1}}}{{{\rm{d}}t}} = \frac{{p_1^n}}{{{I_0}}};\\ \frac{{{\rm{d}}\psi }}{{{\rm{d}}t}} = \left\{ \begin{array}{l} \chi \left( {1 + 2\lambda {Z_1} + 3\mu Z_1^2} \right)\frac{{{\rm{d}}{Z_1}}}{{{\rm{d}}t}},0 \le {\psi _1} \le 1;\\ 0,{\psi _1} > 1; \end{array} \right.\\ \frac{{{\rm{d}}{p_1}}}{{{\rm{d}}t}} = \frac{1}{{{V_1}}}\left( {f{\omega _1}\frac{{{\rm{d}}\psi }}{{{\rm{d}}t}} - {e_{12}}{q_{12}} - {e_{13}}{q_{13}} - \theta \frac{{\delta {Q_1}}}{{{\rm{d}}t}}} \right);\\ {e_{12}} = \left\{ {\begin{array}{*{20}{c}} \begin{array}{l} {c_p}{T_1},\\ {c_p}{T_{\rm{A}}}, \end{array}&\begin{array}{l} {p_1} \ge {p_{\rm{A}}};\\ {p_{\rm{A}}} \ge {p_{\rm{1}}}; \end{array} \end{array}} \right.\\ {e_{13}} = \left\{ {\begin{array}{*{20}{c}} \begin{array}{l} {c_p}{T_1},\\ {c_p}{T_{\rm{D}}}, \end{array}&\begin{array}{l} {p_1} \ge {p_{\rm{D}}};\\ {p_{\rm{D}}} \ge {p_{\rm{1}}}; \end{array} \end{array}} \right.\\ \frac{{\delta {Q_1}}}{{{\rm{d}}t}} = \varepsilon {\rho _1}\left( {{T_1} - {T_{c0}}} \right){S_1};\\ {\rho _1} = \left\{ \begin{array}{l} \frac{{{\omega _1}{\psi _1} - \int_0^t {{q_{12}}{\rm{d}}t} - \int_0^t {{q_{13}}{\rm{d}}t} }}{{{V_1} - \frac{{{\omega _1}}}{\delta }\left( {1 - {\psi _1}} \right) - \alpha {\omega _1}{\psi _1}}},\;\;\;\;\;\;0 \le {\psi _1} \le 1;\\ \frac{{{\omega _1} - \int_0^t {{q_{12}}{\rm{d}}t} - \int_0^t {{q_{13}}{\rm{d}}t} }}{{{V_1} - \frac{{{\omega _1}}}{\delta }\left( {1 - {\psi _1}} \right) - \alpha {\omega _1}{\psi _1}}},\;\;\;\;\;\;{\psi _1} > 1. \end{array} \right. \end{array} \right. $

式中：Z₁、ψ₁、ω₁分别为后药室火药相对厚度、后药室火药燃烧去的百分比、后药室装药量；χ、λ、μ分别为火药药形系数；φ、I₀、n分别为计算次要功系数、压力全冲量、多变指数；θ=γ-1、γ分别为绝热指数；f、δ、α、c_p分别为火药力、火药密度、气体余容、火药气体比定压热容；V₁、S₁分别为后药室初始容积与初始表面积；ε=41.9(J·m)/(kg·s·(°))；T_c0为后药室壁面温度；e₁₂、q₁₂、e₁₃、q₁₃、δQ₁分别为后药室向喷管流入气体比能量、后药室流入喷管流量、后药室向前药室流入的气体比能量、后药室流入前药室流量、后药室对外界散热.

1.2.2 前药室燃烧的内弹道模型

前药室燃烧的内弹道模型为

$ \left\{ \begin{array}{l} \frac{{{\rm{d}}{Z_3}}}{{{\rm{d}}t}} = \frac{{p_3^n}}{{{I_0}}};\\ \frac{{{\rm{d}}{\psi _3}}}{{{\rm{d}}t}} = \left\{ \begin{array}{l} \chi \left( {1 + 2\lambda {Z_3} + 3\mu Z_3^2} \right)\frac{{{\rm{d}}{Z_3}}}{{{\rm{d}}t}},0 \le {\psi _3} \le 1;\\ 0,\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\psi _3} > 1; \end{array} \right.\\ \frac{{{\rm{d}}l}}{{{\rm{d}}t}} = v;\\ \frac{{{\rm{d}}v}}{{{\rm{d}}t}} = \frac{{{A_3}{p_3}}}{{\varphi m}};\\ \frac{{{\rm{d}}{p_3}}}{{{\rm{d}}t}} = \left\{ \begin{array}{l} \frac{1}{{{A_3}\left( {l + {l_{{\psi _3}}}} \right)}}\left( {f{\omega _3}\frac{{{\rm{d}}{\psi _3}}}{{{\rm{d}}t}} + {e_{13}}{q_{13}} + \cdots } \right.\\ \left. { - {e_{34}}{q_{34}} - v\theta \varphi m\frac{{{\rm{d}}v}}{{{\rm{d}}t}} - \theta \frac{{\delta {Q_3}}}{{{\rm{d}}t}}} \right),0 \le l \le {L_0};\\ - \frac{{\gamma {\mu _{\rm{k}}}{{\left( {A - C} \right)}^{\frac{{n + 1}}{{2\left( {n + 1} \right)}}}}{p_3}^{\frac{{3n - 1}}{{2n}}}}}{{{L_0}}}\sqrt {\frac{{\gamma p_{\rm{k}}^{\frac{1}{n}}}}{{{\rho _{\rm{k}}}}}} ,l > {L_0}; \end{array} \right.\\ {p_{\rm{d}}} = {\left[ {A - C{{\left( {\frac{{{L_{\rm{d}}}}}{{{L_0}}}} \right)}^2}} \right]^{\frac{n}{{n - 1}}}}{p_3};\\ {p_{\rm{d}}} = {\left[ {A - C{{\left( {\frac{{{L_{\rm{d}}}}}{{{L_0}}}} \right)}^2}} \right]^{\frac{1}{{n - 1}}}}{p_3};\\ {p_{\rm{D}}} = {A^{\frac{n}{{n - 1}}}}{p_3};\\ {p_{\rm{D}}} = {A^{\frac{1}{{n - 1}}}}{p_3};\\ {l_{{\psi _3}}} = \frac{1}{{{A_3}}}\left[ {{V_3} - \frac{{{\omega _3}}}{\delta }\left( {1 - {\psi _3}} \right) - \alpha {\omega _3}{\psi _3}} \right];\\ \frac{{\delta {Q_3}}}{{{\rm{d}}t}} = \varepsilon {\rho _3}\left( {{T_3} - {T_{{\rm{b0}}}}} \right)\left( {{S_3} + {\rm{ \mathsf{ π} }}{d_0}l} \right);\\ {\rho _3} = \left\{ \begin{array}{l} \frac{{{\omega _3}{\psi _3} + \int_0^t {{q_{13}}{\rm{d}}t} - \int_0^t {{q_{34}}{\rm{d}}t} }}{{{A_3}\left( {l + {l_{{\psi _3}}}} \right)}},0 \le {\psi _3} \le 1;\\ \frac{{{\omega _3} + \int_0^t {{q_{13}}{\rm{d}}t} - \int_0^t {{q_{34}}{\rm{d}}t} }}{{{A_3}\left( {l + {l_{{\psi _3}}}} \right)}},{\psi _3} > 1. \end{array} \right. \end{array} \right. $

式中：Z₃、ψ₃、ω₃分别为前药室火药相对厚度、前药室火药燃烧的百分比、前药室装药量；l_ψ3、V₃、A₃、S₃、L₀、L_d分别为前药室自由容积缩径长、前药室初始容积、线膛横截面积、膛内初始表面积、枪管长度、导气孔距膛底距离；l、v、d₀、m分别为弹丸的行程、速度、直径及质量；T_d、T_b0分别为膛内导气孔处气流温度、枪管壁面温度；p_k、ρ_k分别为膛口气流压力、膛口气流密度；e₃₄、q₃₄、δQ₃、μ_k分别为枪膛向导气室流入的气体比能量、导气孔处流入导气室的流量、枪管内气流对外界散热、膛口气流流量系数；p_k、ρ_k、A、C计算参考文献[13].

1.2.3 后药室喷管气流计算模型

后药室喷管气流计算模型为

$ \left\{ \begin{array}{l} \frac{{\partial \mathit{\boldsymbol{W}}}}{{\partial t}} + \frac{{\partial \mathit{\boldsymbol{F}}}}{{\partial x}} = \mathit{\boldsymbol{B}},\\ \mathit{\boldsymbol{W}} = \left[ \begin{array}{l} {\rho _2}\\ {\rho _2}{v_2}\\ {\rho _2}\left( {\frac{R}{{\gamma - 1}} + \frac{{v_2^2}}{2}} \right) \end{array} \right],\\ \mathit{\boldsymbol{F}} = \left[ \begin{array}{l} {\rho _2}{v_2}\\ {\rho _2}v_2^2 + {p_2}\\ {\rho _2}\left( {\frac{R}{{\gamma - 1}} + \frac{{v_2^2}}{2}} \right){v_2} + {p_2}{v_2} \end{array} \right],\\ \mathit{\boldsymbol{B}} = \left[ \begin{array}{l} - \frac{{{\rho _2}{v_2}}}{{{A_2}}}\frac{{{\rm{d}}{A_2}}}{{{\rm{d}}x}}\\ - \frac{{{\rho _2}v_2^2}}{{{A_2}}}\frac{{{\rm{d}}{A_2}}}{{{\rm{d}}x}} - {\rho _2}{Z_{\rm{f}}}\\ {\rho _2}\frac{{\delta {Q_2}}}{{{\rm{d}}t}} - \left[ {{\rho _2}\left( {\frac{R}{{\gamma - 1}} + \frac{{v_2^2}}{2}} \right){v_2} + {p_2}{v_2}} \right]\frac{1}{{{A_2}}}\frac{{{\rm{d}}{A_2}}}{{{\rm{d}}x}} \end{array} \right]. \end{array} \right. $

式中：v₂、A₂、Z_f分别为喷管内气流速度、管道截面积、表征管壁的摩擦；D_h、λ_f、T_n0分别为管道内径、气体摩擦系数、喷管管壁温度；δQ₂为管壁对管内气体的散热量；Z_f、δQ₂计算参考文献[14].

1.2.4 导气室计算模型

导气室计算模型为

$ \left\{ \begin{array}{l} \frac{{{\rm{d}}{\rho _4}}}{{{\rm{d}}t}} = \frac{{{q_{34}} - {q_{40}} - {\rho _4}{A_{\rm{h}}}{v_{\rm{h}}}}}{{{V_4} + {A_{\rm{h}}}{x_{\rm{h}}}}},\\ \frac{{{\rm{d}}{p_4}}}{{{\rm{d}}t}} = \frac{{\gamma - 1}}{{{V_4} + {A_{\rm{h}}}{x_{\rm{h}}}}}\left[ { - \frac{{\delta {Q_4}}}{{{\rm{d}}t}} + {e_{34}}{q_{34}} - {e_{40}}{q_{40}} - \frac{{\gamma {p_4}{A_{\rm{h}}}{v_{\rm{h}}}}}{{\gamma - 1}}} \right],\\ \frac{{\delta {Q_4}}}{{{\rm{d}}t}} = \varepsilon {\rho _4}\left( {{T_4} - {T_{{\rm{q0}}}}} \right)\left( {{S_4} + {\rm{ \mathsf{ π} }}{d_{\rm{h}}}{x_{\rm{h}}}} \right),\\ \frac{{{\rm{d}}{v_h}}}{{{\rm{d}}t}} = \frac{1}{{{m_{\rm{h}}}}}\left( {{p_4} - {p_0}} \right){A_{\rm{h}}} - {K_{\rm{c}}}{x_{\rm{h}}} - {P_{\rm{c}}} - {R_{\rm{c}}},\\ \frac{{{\rm{d}}{x_{\rm{h}}}}}{{{\rm{d}}t}} = {v_{\rm{h}}}. \end{array} \right. $

式中：p₀为大气压力；A_h、V₄、S₄、d_h分别为活塞横截面积、导气室初始容积、导气室初始表面积、活塞筒名义直径；v_h、x_h分别为活塞速度、活塞位移；m_h、M_c、K_c、P_c、R_c分别为活塞质量、机框组件质量、机框复进簧刚度、机框复进簧预压力、机框与活塞所受摩擦阻力；δQ₄、q₄₀、e₄₀分别为导气室气流对外界散热、导气室漏气流量、从导气室间隙漏出的单位质量气体所具有的能量.

1.2.5 气室之间流量方程

气室之间流量方程为

$ {q_{12}} = \left\{ \begin{array}{l} \frac{{{\rm{ \mathsf{ π} }}d_{\rm{p}}^2{\mu _{12}}}}{4}{\left( {\frac{2}{{\gamma + 1}}} \right)^{\frac{{\gamma + 1}}{{2\left( {\gamma - 1} \right)}}}}\sqrt {\gamma {p_1}{\rho _1}} ,{p_{\rm{A}}} \le \zeta {p_1};\\ \frac{{{\rm{ \mathsf{ π} }}d_{\rm{p}}^2{\mu _{12}}}}{4}\sqrt {\frac{{2\gamma }}{{\gamma - 1}}{p_1}{\rho _1}\left[ {{{\left( {\frac{{{p_{\rm{A}}}}}{{{p_{\rm{1}}}}}} \right)}^{\frac{2}{\gamma }}} - {{\left( {\frac{{{p_{\rm{A}}}}}{{{p_{\rm{1}}}}}} \right)}^{\frac{{\gamma + 1}}{\gamma }}}} \right]} ,\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\zeta {p_1} < {p_{\rm{A}}} < {p_1};\\ - \frac{{{\rm{ \mathsf{ π} }}d_{\rm{p}}^2{\mu _{12}}}}{4}\sqrt {\frac{{2\gamma }}{{\gamma - 1}}{p_{\rm{A}}}{\rho _{\rm{A}}}\left[ {{{\left( {\frac{{{p_{\rm{1}}}}}{{{p_{\rm{A}}}}}} \right)}^{\frac{2}{\gamma }}} - {{\left( {\frac{{{p_{\rm{1}}}}}{{{p_{\rm{A}}}}}} \right)}^{\frac{{\gamma + 1}}{\gamma }}}} \right]} ,\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\zeta {p_{\rm{A}}} < {p_{\rm{1}}} < {p_{\rm{A}}};\\ - \frac{{{\rm{ \mathsf{ π} }}d_{\rm{p}}^2{\mu _{12}}}}{4}{\left( {\frac{2}{{\gamma + 1}}} \right)^{\frac{{\gamma + 1}}{{2\left( {\gamma - 1} \right)}}}}\sqrt {\gamma {p_{\rm{A}}}{\rho _{\rm{A}}}} ,{p_{\rm{1}}} \le \zeta {p_{\rm{A}}}. \end{array} \right. $

$ {q_{13}} = \left\{ \begin{array}{l} \frac{{{\rm{ \mathsf{ π} }}d_{{\rm{13}}}^2{\mu _{13}}}}{4}{\left( {\frac{2}{{\gamma + 1}}} \right)^{\frac{{\gamma + 1}}{{2\left( {\gamma - 1} \right)}}}}\sqrt {\gamma {p_1}{\rho _1}} ,{p_{\rm{D}}} \le \zeta {p_1};\\ \frac{{{\rm{ \mathsf{ π} }}d_{{\rm{13}}}^2{\mu _{13}}}}{4}\sqrt {\frac{{2\gamma }}{{\gamma - 1}}{p_1}{\rho _1}\left[ {{{\left( {\frac{{{p_{\rm{D}}}}}{{{p_{\rm{1}}}}}} \right)}^{\frac{2}{\gamma }}} - {{\left( {\frac{{{p_{\rm{D}}}}}{{{p_{\rm{1}}}}}} \right)}^{\frac{{\gamma + 1}}{\gamma }}}} \right]} ,\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\zeta {p_1} < {p_{\rm{D}}} < {p_1};\\ - \frac{{{\rm{ \mathsf{ π} }}d_{{\rm{13}}}^2{\mu _{13}}}}{4}\sqrt {\frac{{2\gamma }}{{\gamma - 1}}{p_{\rm{D}}}{\rho _{\rm{D}}}\left[ {{{\left( {\frac{{{p_{\rm{1}}}}}{{{p_{\rm{D}}}}}} \right)}^{\frac{2}{\gamma }}} - {{\left( {\frac{{{p_{\rm{1}}}}}{{{p_{\rm{D}}}}}} \right)}^{\frac{{\gamma + 1}}{\gamma }}}} \right]} ,\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\zeta {p_{\rm{D}}} < {p_{\rm{1}}} < {p_{\rm{D}}};\\ - \frac{{{\rm{ \mathsf{ π} }}d_{{\rm{13}}}^2{\mu _{13}}}}{4}{\left( {\frac{2}{{\gamma + 1}}} \right)^{\frac{{\gamma + 1}}{{2\left( {\gamma - 1} \right)}}}}\sqrt {\gamma {p_{\rm{D}}}{\rho _{\rm{D}}}} ,{p_{\rm{1}}} \le \zeta {p_{\rm{D}}}. \end{array} \right. $

式中：ζ为临界压力比；μ₁₂、μ₁₃、μ₃₄、μ₄₀分别为喷管导气孔处流量系数、传火孔处流量系数、导气孔处流量系数、活塞处漏气系数；d_p、d₁₃、d分别为喷管导气孔直径、传火孔直径、导气孔直径；q₃₄、q₄₀计算参考文献[15].

1.3 减后坐效率及喷管反推力计算

计算武器减后坐效率是评测武器减后坐装置好坏的依据，喷管反推力则反映喷管气流对武器作用效果，是武器系统动力学仿真的基础载荷之一.

通过对比有、无双药室喷管装置的武器后坐动量，可将双药室武器的喷管减后坐效率定义为

$ \left\{ \begin{array}{l} {\eta _{{\rm{pg}}}} = \frac{{{I_0}\left( p \right) - \left[ {{I_{{\rm{pg}}}}\left( {{p_3}} \right) - {I_{{\rm{pg}}}}\left( {{F_{{\rm{pg}}}}} \right)} \right]}}{{{I_0}\left( p \right)}},\\ {I_0}\left( p \right) = \int_0^{{t_n}} {p\left( {{t_i}} \right){A_3}{\rm{d}}t} ,\\ {I_{{\rm{pg}}}}\left( {{p_3}} \right) = \int_0^{{t_n}} {{p_3}\left( {{t_i}} \right){A_3}{\rm{d}}t} ,\\ {I_{{\rm{pg}}}}\left( {{F_{{\rm{pg}}}}} \right) = \int_0^{{t_n}} {{F_{{\rm{pg}}}}\left( {{t_i}} \right){\rm{d}}t} . \end{array} \right. $

式中：I₀(p)为无双药室喷管武器膛内平均压力气体冲量；I_pg(p₃)为双药室喷管武器膛内平均压力气体冲量；I_pg(F_pg)为双药室喷管武器喷管反推力冲量.

喷管气流反推力随时间变化曲线通过气体数学模型及基本假设编程计算，程序中计算按下式计算:

$ \begin{array}{l} {F_{{\rm{pg}}}}\left( {{t_{i + 1}}} \right) = {F_{{\rm{pg}}}}\left( {{t_i}} \right) + {p_2}\left( {{t_i},{L_j}} \right) \cdot \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\left[ {{A_2}\left( j \right) - {A_2}\left( {j - 1} \right)} \right]. \end{array} $

式中：F_pg(t_i)为t_i时刻喷管气流反推力；p₂(t_i, L_j)为t_i时刻喷管截面位于L_j位置处喷管气流压力；A₂(j)为喷管截面位于L_j位置处喷管截面积.

2 模型验证及仿真结果分析 2.1 模型验证

表 1为双药室喷管武器预置参数.对气体数学模型及基本假设方程组通过Matlab进行编程计算，得到单发武器系统射击的后药室平均压力曲线、喷管气室压力曲线、膛内平均压力曲线、导气室压力曲线，如图 3所示.

由图 3可知：当t=(0~0.51) ms时后药室火药处于定容燃烧阶段，可看成密闭爆发器模型；当t₁=0.51 ms时后药室气体压力达到破膜压力10 MPa，后药室气体分别进入喷管气室、前药室；当t₂=1.02 ms时前药室气体压力达到弹丸启动压力30 MPa，弹丸开始挤进膛线，此刻后药室气体压力达到50.8 MPa，而喷管气室压力达到19.9 MPa，武器总体受力方向向后；当t₃=1.37 ms时前药室气体压力达到最大膛压79.9 MPa；当t₄=1.64 ms时后药室气体达到最大压力138.3 MPa；当t₅=1.73 ms时喷管气室达到最大压力71.1 MPa；t₆=1.89 ms时弹丸经过导气孔，膛内部分火药气体经导气孔进入导气室；当t₇=2.40 ms时导气室内压力达到最大值17.3 MPa；当t₈=3.39 ms时内弹道时期结束开始进入后效期，膛内平均压力曲线开始出现拐点，此时弹丸膛口速度为313.6 m/s；当t₉=8.50 ms时后药室气体压力降至0.5 MPa.

试验方法：将双药室喷管低后坐武器通过夹具安装于试验滑移架上，滑移架可在试验固定架座沿枪管轴线方向自由移动.试验过程时，枪身后坐带动滑移架尾端挤压电传感器.力传感信号通过DWE-43微型多功能数据采集器采集，由奥地利公司分析软件DEWESoft-DSA得出后坐力测试数据.通过天幕靶测得弹丸初速，膛压传感器测量后药室平均压力、喷管气室压力、前药室平均压力和导气室压力.图 4为天幕靶和高速摄影捕捉榴弹击发瞬间试验图.图 5为试验测得的武器后坐力F_hz曲线.

从表 2可以看出压力参数、弹丸初速和后坐力的仿真值与试验值相对误差在10%范围内，表明建立的双药室后喷装置数学模型是可信的，计算的后药室平均压力、喷管气室压力、前药室平均压力及弹丸初速是正确合理的.由图 5可知：B1处波峰为武器射击过程中最大后坐力(1 192.0 N)，主要由于弹丸击发后枪身后坐造成的；B2处波峰(533.2 N)主要由于机框后坐到位撞击机匣造成；A1处小波峰测试过程中干扰信号造成；A2、A3处波峰为武器射击过程中滑移架座自身振荡撞击压力传感器所致.

2.2 有、无双药室后喷装置武器发射性能对比

为了解采取双药室后喷减后坐装置武器与普通武器的发射性能的区别，分别对两种武器的内弹道及后效期进行数值计算，得到膛压曲线、导气室压力曲线和弹丸速度曲线，如图 6、7所示.前药室火药在后药室压力达到破膜压力(t₁=0.51 ms)时才开始燃烧，为了更直观区别两种武器压力曲线、速度曲线差别，将普通武器的膛压曲线、导气室压力曲线、速度曲线分别延时Δt=0.51 ms得到对应的压力曲线.

由图 6可知, 延迟普通武器的最大膛压值为76.6 MPa，略小于双药室武器的最大膛压值79.9 MPa，主要由于双药室武器后药室气体压力大于前药室压力，部分后药室火药燃气进入前药室；双药室武器膛压曲线达到最大膛压后衰减更平缓，结合图 3，主要由于在弹丸膛内运动期间后药室压力始终大于前药室压力，气流一直保持向前药室流入状态.双药室武器导气室压力曲线起始时刻早于普通武器的导气室压力曲线起始时刻，由于双药室武器膛内压力对弹丸做功多，弹丸速度大，故到达导气孔时间短.由图 7可知, 双药室武器弹丸膛口速度值为313.6 m/s，略大于普通武器弹丸膛口速度值309 m/s.将图 6中膛压曲线对时间积分，再与枪膛截面积相乘，得I₀(p)=67.59 N·s，I_pg(p₃)=68.43 N·s；将后药室装药量w₁=16 g对应的喷管反推力对时间积分，得I_pg(F_pg)=23.85 N·s.将上述结果代入喷管减后坐效率公式得η_pg=34%.可见采用双药室喷管减后坐技术的武器，对武器弹丸初速影响不大，制退效率高.

2.3 发射过程数值仿真结果及分析 2.3.1 双药室后喷装置气流参数随时间变化曲线

图 8为双药室导气管内流场状态参量在不同位置处随时间变化的曲线.整体来说，喷管气流温度、喷管气流速度受压强间断面影响较大，相应的变化曲线出现激增现象，而喷管气流压力曲线、喷管气流密度曲线随时间变化过度较平滑.

由图 8(a)可知，喷管截面距离喷管起始点B越近，喷管截面处达到气流压力峰值越大，气流压力峰值出现在t=(1.73~2.13) ms附近，而后迅速衰减；喷管截面位于最小截面D左侧(靠近喷管起始点)时，喷管截面距离喷管起始点越近，气流压力峰值出现时刻越早，压力波向喷管右侧传播；喷管截面位于最小截面D右侧(靠近喷管出口)时，距离喷管起始点越近，气流压力峰值越滞后，压力波向喷管左侧传播；喷管起始点在整个发射过程中气流压力变化剧烈，而喷管出口处I在整个发射过程中气流压力变化平缓.

由图 8(b)可知，喷管截面位于喷管起始点B时，气流温度变化曲线较平滑，主要在建模时将喷管入口处作为一个有初始容积的气室考虑，将B处气流参数作为喷管计算的初始边界条件，并没有考虑压强间断面的影响；喷管截面位于喷管起始点B与喷管最小截面D之间，越靠近喷管D处喷管气流温度峰值越高；喷管截面位于截面D与喷口I之间，越靠近喷管出口I处喷管气流温度峰值越低.

由图 8(c)可以看出喷管气流密度曲线变化趋势与喷管气流压力曲线变化趋势相同，喷管截面位置越靠近喷管喷口，管内气流密度越小.

由图 8(d)可知，喷管截面越靠近喷管喷口，管内气流速度峰值越大，且峰值出现的时间逐渐滞后，表明管内气流速度一直处于加速状态；在整个射击过程中，管内气流速度激增至最大值而后缓慢衰减.

2.3.2 双药室后喷装置气流参数随位置分布曲线

图 9为双药室后喷导气管内流场状态参量在喷管内分布曲线.

由图 9(a)可知，当t=(0.51~1.02) ms时，火药燃气进入喷管，喷管左侧压力高于右侧压力，压力波在喷管内逐渐向喷管口部传播，且随着时间增大，喷管内气体压力逐渐增大，但压力波未完全传至喷管口部；当t=(1.02~1.73) ms时，由于喷管气室气体压力逐渐增大直至达到压力峰值，喷管内呈现喷管两侧压力差增大趋势，当压力波传至喷管喉径D处后，喷管截面开始扩张，管内气流压力迅速衰减；当t=(1.73~3.39) ms时，随着时间增加，由于喷管气室压力降低，喷管两侧压力差逐渐缩小，当弹丸出膛口(t=3.39 ms)时喷管两侧压力差几乎为0；当t=(3.39~8.50) ms时，随着时间增加，喷管内气流压力总体下降.

由图 9(b)可知，当t=(0.51~1.02) ms时，火药燃气刚进入喷管气室，喷管内温度保持常温不变，当喷管气室气体进入喷管，靠近喷管起始位置喷管气流温度迅速升至温度峰值，且随着时间推移喷管内的气流温度逐渐升高；当t=(1.02~1.73) ms时，喷管气流温度在喷管内分布几乎保持不变，靠近喷管左侧的气流温度略高于右侧气流温度，喷管气流经过喷管最小截面D处后气流温度迅速下降，气流经扩张管后气流温度以一恒温向后坐方向喷出；当t=(1.73~8.50) ms时，气流温度在喷管内左侧温度高于右侧温度，且随时间增大喷管内气流温度逐渐降低，喷管两侧温度差变小.

由图 9(c)可以看出，当t=(0.51~1.02) ms时，当火药气体流入喷管，在管内气流密度左侧大于右侧，但压力波未传至处密度保持不变；当t=(1.02~ 1.73) ms时，喷管截面距喷管起始点越近，喷管截面处气流密度越大，当喷管气流经过截面D后密度迅速变小，且随时间增大进入喷管内火药燃气增多，喷管内气流密度增大；当t=(1.73~1.89) ms时，喷管内气流密度在喷管起始点处密度相同，喷管最小截面处气流密度随时间增大而增大；当t=(1.89~ 3.39) ms时，喷管最小截面处气流密度在上述时间段保持不变，喷管在起始点处喷管气流密度随时间增大而减小；当t=(3.39~8.50) ms时，喷管气流密度在最小截面左侧时几乎保持相同，喷管气流密度在最小截面右侧随喷管扩张急剧降低.

由图 9(d)可知，当t=(0.51~1.02) ms时，喷管气流还未冲出导气管，但气流速度随压力波到达喷管位置增大而增大；当t=(1.02~1.73) ms时，喷管气流在最小截面D处左侧时，管内气流速度缓慢增大，但经过D处进入扩张管后，管内气流速度骤增经扩张管加速后以恒速向后喷出；当t=(1.73~ 8.50) ms时，喷管内气流速度任一截面随时间增加而降低，喷管在此时间段内气流加速度减小.

3 双药室喷管参数对发射性能影响

喷管减后坐武器利用向后坐方向喷出火药气体动量降低后坐力，提高火药气体流量及增大火药气体速度可以增大喷管气流反推力(F_pg)、提高武器减后坐效率(η_pg).本文将从后药室装药量(ω₁)、喷管导气孔直径(d_p)、喷管内径(D_h)3个方面研究对武器发射性能的影响.

3.1 后药室装药量对发射性能影响分析

由图 10可知，随着ω₁增大，后药室火药燃气达到破膜压力时间缩短，越早进入喷管及前药室，F_pg作用时间越早；F_pg峰值随着ω₁增大而增大，对武器反后坐冲量越大.由表 3可知：增大ω₁，故后药室气体压力增大，流入喷管及前药室气体流量增大，F_pg冲量增大，前药室气体压力冲量增大，弹丸出膛口速度增大；ω₁=16 g的F_pg峰值是ω₁=12 g的1.48倍；ω₁由12 g增加至16 g，η_pg提升了32.1%，而v₀提升0.8%，可见增加ω₁对于前药室内弹道的改变可以忽略不计，但是大大提升了武器的η_pg.

3.2 喷管导气孔直径对发射性能影响分析

由图 11可知，随着d_p由6 mm增大至14 mm，F_pg峰值逐渐增大，结合表 5，F_pg峰值在d_p=12 mm时是d_p为6 mm时的3.16倍，增加幅度较大.由表 4可知，随着d_p增大，后药室火药流入喷管流量增大，F_pg冲量增大；当F_pg达到峰值后开始衰减时，d_p越小，后药室火药气体进入喷管流量变小，后药室压力越大，进入前药室气体流量增多，所以前药室气体冲量越大，v₀越高；d_p由6 mm增加2.3倍至14 mm，武器η_pg提升了42.9%，v₀下降了4.2%.

3.3 喷管内径对发射性能影响分析

由图 12可知，D_h越大，导气管内聚集的火药燃气的能力越大，单位时间内经喷管加速的火药燃气越多，F_pg越大，反后坐冲量越大，且F_pg峰值在D_h=(6~12) mm范围内较敏感.由表 5可知，当D_h由6 mm增加至14 mm，F_pg冲量增加了46%，前药室膛内气体压力冲量减小了2%，v₀降低了2.3%，武器η_pg提升了32.1%.

4 结论

1) 利用火药气体能量提出一种新型双药室减后坐装置，对设计的双药室减后坐方案进行实弹射击试验，v₀达到316.7 m/s时武器Max(F_pg)仅为1 192 N.

2) 将计算得到双药室喷管低后坐武器的Max(p₁)、Max(p_A)、Max(p₃)、v₀、Max(F_pg)与对应的试验结果进行对比，仿真值相对试验值误差在8.6%范围内，表明建立的计算模型是正确合理的.本文采用的双药室减后坐装置传火特性：在0.51 ms前，后药室处于定容燃烧；在0.51 ms时刻，p₁达到铝膜破膜压力10 MPa，后药室火药经过传火孔同时进入前药室和喷管气室；在1.02 ms时刻，p₃达到弹丸启动压力30 MPa，此时武器总体受力方向向后；在1.64 ms时在后药室位置处武器承受火药最大压力，压力峰值为138.3 MPa；火药燃烧直至后效期结束持续时间为8.5 ms左右.

3) 分析喷管气流参数随时间变化规律：温度、速度受压强间断面影响较大，相应的变化曲线出现激增现象，而压力、密度随时间变化过度较平滑；温度与速度曲线变化趋势相同，而压力与密度曲线变化趋势相同.分析喷管气流参数在喷管内分布状况：当压力波未传至喉径时，气流参数(压力、温度、密度)在喷管左侧时大于右侧，而气流速度在喷管右侧小于左侧；当压力波传至喉径时，喷管气流参数(压力、温度、密度)急剧衰减，而气流被加速从喷管口部高速喷出.

4) 增加ω₁、扩大d_p和D_h可以增大F_pg、提升η_pg.ω₁在(12~16) g范围内增加时，η_pg增幅显著；d_p和D_h在(6~14) mm范围内增大时，对提升η_pg越来越不显著.