基于平衡算法的无坐力炮不平衡冲量研究
分类：主编推荐 发布时间：2022-06-24 14:08 访问量：198
基于平衡算法的无坐力炮不平衡冲量研究
1 引言
无坐力炮是一种滑膛无反后坐装置的火炮，在火炮发射时与闭膛炮有明显区别，即有气流从尾部喷管流出以抵消膛内气体推动弹丸时所产生的后坐力，同时，还是会产生不平衡冲量。不平衡冲量对单兵无坐力炮武器性能的安全性有较大影响，因此对无坐力炮在发射过程中的不平衡冲量进行理论计算和试验测量是十分必要的。
为了研究无坐力炮的不平衡冲量，文献[1]通过建立准两相流内弹道模型得到了无坐力炮的不平衡冲量随时间的变化规律；文献[2]通过有限元软件分析了平衡体对平衡炮不平衡冲量的影响；文献[3]研究了带状变燃速发射药对无坐力炮不平衡冲量的影响；文献[4]对影响单兵武器发射的不平衡冲量进行了分析；文献[5]得出了75 mm无坐力炮不平衡冲量与喷管结构的关系；文献[6]通过自行设计的试验，测得了某无坐力炮的不平衡冲量。文献[7]建立了无坐力炮的一维两相流内弹道模型，但未对不平衡冲量进行研究。
本文利用经典无坐力炮平衡内弹道算法，建立了无坐力炮不平衡冲量的计算模型，并分别从弹丸挤进压力、堵片打开压力、火药弧厚、喷管面喉比等角度分析了初始设计参数对不平衡冲量的影响。通过该计算模型可以得到各设计参数对无坐力炮不平衡冲量的影响大小，可为无坐力炮不平衡冲量的优化提供了一定的参考。
2 数学物理模型
经典无坐力炮平衡内弹道算法将内弹道过程划分为前期、第一时期、第二时期这3个不同时期。这3个时期之间是相互连接的，前期的最终结果是第一时期的起始条件，第一时期的最终结果为第二时期的起始条件。按顺序分别作出各阶段的解，并在后两个过程中加入不平衡的计算。另外应用经典无坐力炮平衡内弹道算法还需作出如下简化和假设[8]：
1) 火药燃烧在平均压力下进行且服从正比燃速定律；
2) 在整个内弹道过程中，火药力f、余容α及绝热指数k均为常数；
3) 因点火药在整个装药能量中占比很小，将点火药质量折合成主装药质量考虑；
4) 弹丸挤进膛线和堵片打开是瞬时进行的，分别以挤进压力p0和堵片打开压力pom表示弹丸挤进和堵片打开。
2.1 无坐力炮内弹道方程
2.1.1 前期
在这一时期中，火药一直在容积为V的药室中燃烧，容积内的压力由点火药压力PB升高到弹丸起动压力P0[9]，求出此阶段结束后的火药形状尺寸参数ψ0、σ0、Z0。
(1)
(2)
(3)
得到的ψ0、σ0、Z0既是前期计算的最终结果，也是第一时期的起始条件。
2.1.2 第一时期计算
在这一时期中，弹丸达到挤进压力P0并开始向前运动而堵片未打开。在一般无坐力炮经典内弹道方程组中[7]，共有6个方程，有p、v、l、t、ψ、η及Z七个变量，在这7个变量中，只有ψ和Z的边界是已知的，即ψ的范围是ψ0到1，Z的范围是Z0到1，设x=Z-Z0作为自变量，结合一般无坐力炮经典内弹道方程组分别得出函数v、ψ、l、p的表达式：
1) 弹丸速度方程
(4)
2) 几何燃烧定律方程
ψ=ψ0+K1x+χλx2
(5)
式中K1=χσ0
3) 弹丸行程方程
(6)
式中
4) 压力方程
(7)
上述式中：v为弹丸运动速度；Z为火药相对已燃厚度；s为身管截面积；Ik为发射药的压力全冲量；φ′为第一时期次要功系数；m为弹丸质量；ψ为火药相对已燃质量；χ和λ为火药形状特征量；l为弹丸运动行程；p为膛内平均压力；f为火药力；ω为火药质量。
在这一时期中，当药室底部压力达到堵片打开压力pom时，记下此时刻的v1、ψ1、l1、p1并带入第二时期的计算。
2.1.3 第二时期计算
设第一时期结束诸元分别为Z1、v1、ψ1、l1、p1，取
为自变量，第二时期方程如下：
1) 弹丸速度方程
(8)
2) 几何燃烧定律方程
(9)
式中 K2=χ(1+2λZ1)
3) 气体流量方程
(10)
式中
4) 弹丸行程方程
(11)
式中
5) 压力方程
(12)
上述式中：φ″为第二时期次要功系数； φ2为流量损失系数；skp为喷管喉部面积；A为流量系数。
2.2 不平衡冲量的计算
不平衡冲量是一个过程量，但在本文1.1中的内弹道方程组中只能得到p-l、v-l之间的关系，所以在进行不平衡冲量计算之前，还应解出p-t、v-t之间的对应关系。将求出的v-l曲线做出1/v-l曲线进行如下的图解积分，计算与各l值相应的t值。由于当l接近于0时，1/v将趋近于无穷大，所以在起始段应采用近似方法来处理。先将t分成两段来处理，即
(13)
式中的t′是按下式来确定的：
(14)
其中v′/2代表v=0和v=v′之间的速度平均值，为了尽可能减小误差，必须在v-l曲线中取较小的l′值，当l′取得合适时，误差可以忽略不记。这样可以得出p-t、v-t之间的关系。
在第一时期中，弹底压力达到弹丸挤进压力p0并推动弹丸向前运动，则因弹丸运动所产生的后坐冲量为
(15)
式中:pd为第一时期弹底压力，s为身管横截面积，φ′为第一时期次要功系数。
在第二时期中，当药室底部压力达到堵片打开压力时pom时，这时炮身的不平衡冲量由弹丸运动和喷管流出的气体共同作用，设沿炮轴朝向炮尾的方向为合力的正方向，则炮身受力如图1所示。
图1 第二时期炮身受力示意图
Fig.1 Explosion diagram of the second period
当弹丸以速度v运动时，火药气体在药室压力p的作用下流经喷管的喉部截面skp后膨胀，并在通过喷管外部截面sa后具有压强pa，这里气体未充分膨胀，所以pa比一般大气压强要大得多，忽略外界大气压对其的影响，在喷管出口截面气体的流出速度为Ua。
将整个发射器看成是一个整体，根据内力作用下的系统各部分的动量等于外力冲量的力学原理对整个发射器进行分析[9]：
1) 内力作用下的动量：弹丸的动量为-m弹dv，经喷管流出的气体动量为
为火药气体每秒流出量，Ua为火药气体在喷管出口的流动速度。
2) 外力作用下的冲量：在喷管外部截面中的气体压力冲量为-pasadt。
根据内力与外力的分析，得到
(16)
火药气体对炮身的作用力可以表示为[9]：
(17)
则第二时期的不平衡冲量为
(18)
上述式中，φ″为第二时期的次要功系数；φ1为流速损失系数；φ2为流量修正系数； ξ为喷管推力系数，它是绝热指数K、skp/sa的函数，其数值详见文献[10]；sa为喷管出口截面积；skp为喷管喉部截面积；p总为膛内总静压。
3 试验及计算结果分析
3.1 试验方案
试验现场布置如图2所示，将发射器放置在自制的炮架上，使发射器与滑块相固连，滑块可在滑轨上沿膛线方向前后自由移动。通过位移传感器测得发射器前后运动位移随时间的变化，对位移数据进行一阶求导得到发射器沿轴线的运动速度，再与运动质量相乘得到不平衡冲量。试验的初始参数如表1、表2所示，并在相同的工况下对无坐力炮的发射过程进行3次不平衡冲量的测量。
图2 不平衡冲量试验布置示意图
Fig.2 Unbalance impulse test arrangement
表1 火炮构造诸元参数
Table 1 artillery structure elements
表2 火药装填参数
Table 2 Conditions of gunpowder loading
3.2 计算结果分析
将上述表格中的火炮初始参数带入经典无坐力炮平衡内弹道方程，通过数值仿真计算，不平衡冲量计算结果与试验结果对比如表3所示。结果表明，通过此平衡内弹道得到的不平衡冲量I与试验结果符合较好，相对误差均在10%以内，可用于指导工程设计。
表3 不平衡冲量计算结果与试验结果
Table 3 Comparison of calculated results and experimental resultsof unbalance impulse
4 初始设计参数对不平衡冲量的影响研究
在验证经典无坐力炮平衡内弹道算法合理可行之后，下面对某新型无坐力炮分别从弹丸挤进和堵片打开压力、火药弧厚、喷管面喉比等设计参数对其不平衡冲量进行分析。
4.1 弹丸挤进和堵片打开压力对不平衡冲量的影响
采取单一变量原则，在其他条件不变的情况下，根据2.1中的内弹道方程组分别计算当弹丸挤进压力po为 4 MPa、5 MPa、6 MPa、7 MPa、8 MPa时的p-t、p-l、v-t、v-l的关系，再根据式(15)、式(18)计算出这5种工况下的不平衡冲量I。计算结果如图3所示。
图3 不平衡冲量与弹丸挤进压力关系曲线
Fig.3 Relation diagram between projectile extrusion pressure and unbalance impulse
当弹丸挤进压力po为7 MPa，在其他条件不变的情况下，计算堵片打开压力pom分别取8 MPa、9 MPa、10 MPa、11 MPa、12 MPa这5种工况下的不平衡冲量I。计算结果如图4所示。
结合图3和图4可以看出，当弹丸挤进压力增大，逐渐接近于堵片打开压力时，不平衡冲量呈减小趋势；不平衡冲量随着堵片打开压力的增加而变大；弹丸挤进压力与堵片打开压力逐渐接近时，不平衡冲量呈减小趋势。
图4 不平衡冲量与堵片打开压力关系曲线
Fig.4 Relationship between unbalanced impulse and opening pressure of blockage
4.2 火药弧厚对不平衡冲量的影响
采取单一变量原则，在其他条件不变的情况下，根据2.1节中的内弹道方程组分别计算当火药弧厚分别为0.26 mm、0.28 mm、0.30 mm、0.32 mm、0.34 mm时的p-t、p-l、v-t、v-l 的关系，再根据式(15)、式(18)计算出这5种工况下的不平衡冲量I。计算结果如图5所示。
图5 不平衡冲量与火药弧厚关系曲线
Fig.5 Relation between unbalance impulse and thickness of gunpowder arc
由图5可以看出，火药弧厚越厚不平衡冲量越小。但由于火药弧厚对火炮膛压与炮口初速的影响较大，所以在满足膛压和初速指标的情况下，选用较厚的火药弧厚有利于不平衡冲量的下降。
4.3 喷管面喉比对不平衡冲量的影响
采取单一变量原则，在其他条件不变的情况下，根据2.1节中的内弹道方程组计算当喷管的面喉比分别为2、2.5、3.33、5、10时的p-t、p-l、v-t、v-l的关系，再根据式(15)、式(18)计算出这5种工况下的不平衡冲量I。计算结果如图6所示。
图6 不平衡冲量与喷管面喉比关系曲线
Fig.6 Relation diagram between unbalanced impulse and nozzle face throat ratio
从图6可以看出喷管面喉比的变化对不平衡冲量的影响很大，随着面喉比的增大，不平衡冲量由后坐变为前冲，且变化幅度较大，这说明选择合适的面喉比在无坐力炮的不平衡力初期设计和试验中起着重要作用。
5 结论
1) 喷管面喉比是影响无坐力炮平衡的重要参数，对冲量影响较大，在火炮的设计及试验时应重点考虑并进行调试。
2) 弹丸挤进压力和堵片打开压力对无坐力炮冲量大小起着调节作用，当两压力值接近相等时，有利于减小不平衡冲量，在火炮试验的中后期应重点关注。
3) 在满足膛压和炮口初速的设计指标前提下，选用较厚的火药弧厚有利于降低不平衡冲量。
对影响冲量的参数进行合理的匹配，可优化某新型无坐力炮的平衡性能，并可对其他身管武器的平衡设计提供参考。
[1] 张莺.无后坐炮推力及后坐冲量的计算[C]//中国兵工学会弹道学术讨论会论文集.桂林:1985：9.
Zhang Y.The calculation of thrust and impulse of recelless gun[C]//Ballistic Symposium of Chinese Society of Arm Engineering.Guilin:1985:9.
[2] 章勇，杨臻，吴臣.平衡炮内弹道不平衡冲量技术研究[J].火力与指挥控制,2019，44(02):102-107.
Zhang Y，Yang Z，Wu C.Research on unbalanced impulse technology of internal ballistics of balanced gun[J].Fire Control & Command Control，2019，44(02):102-107.
[3] 薛百文，马营，涂炯灿，等.带状变燃速发射药对无后坐力炮内弹道性能影响研究[J].火炮发射与控制学报，2020，41(03):6-11.
Xue B W，Ma Y，Tu J C.Research on the influence of ribon-shaped variable burning rate projector on the internal ballistic performance of non-recoil gun[J].Journal of Gun Launching and Control，2020，41(03):6-11.
[4] 陈志坚，林辉，张海凌.单兵武器发射器动不平衡冲量分析与计算[J].弹箭与制导学报，2008，28(06):164-165.
Chen Z J，Lin H，Zhang H L.Analysis and calculation of dynamic unbalance impulse of individual weapon launcher[J].Journal of Projectiles，Missiles and Guidance，2008，28(06):164-165.
[5] U.S.Army Materiel Command.Engineering design handbook:recoilless rifle weapon systems:AD-A023513[M].Alexandria，VA，US:National Technical Information Service，U.S.Department of Commerce，1976.
[6] 庞春桥，陶钢，李召，等.轻型无后坐力炮的动不平衡冲量特性[J].兵工学报，2020，41(12):24-31.
Pang C Q，Tao G，Li Z.Dynamic unbalance impulse characteristics of light non-recoil gun[J].Acta Armamentarii，2020，41(12):24-31.
[7] 付丽华，阮文俊，无后坐炮一维两相流内弹道参量分析[J].兵器装备工程学报，2019，40(07):51-55.
Fu L H，Ruan W J.One-dimensional two-phase flow internal ballistic parameters analysis of recarless gun[J].Journal of Ordnance and Equipment Engineering，2019，40(07):51-55.
[8] 金志明.枪炮内弹道学[M].北京:北京理工大学出版社，2004.
Jin Z M.Gun internal ballistics[M].Beijing:Beijing Institute of Technology Press，2004.
[9] 谢列伯梁可夫.M.E.身管武器和火药火箭内弹道学[M].北京:国防工业出版社，1965.
Serebelliangkov M E.Barreled weapon and gunpowder rocket internal ballistics[M].Beijing:National Defense Industry Press，1965.
[10] 张东妍.无座力炮设计手册[M].北京:国防工业出版社，1978.
Zhang D Y.Design manual of recoilless gun[M].Beijing:National Defense Industry Press，1978.
Study on unbalanced impulse of recoilless gun based on balance algorithm
Citation format:XUE Juyuan, TAO Gang, WANG Peng, et al.Study on unbalanced impulse of recoilless gun based on balance algorithm[J].Journal of Ordnance Equipment Engineering，2022，43(04)：1-6.