膨胀波在火炮膛内传播规律的捕捉计算分析王颖泽，张小兵，袁亚雄(南京理工大学动力工程学院，南京210094)摘要：为了准确的获得膨胀波在火炮膛内的传播规律，根据膨胀波火炮的发射特点，分析了膨胀波火炮发射过程中膛内复杂波系的组成及传播特性；利用膨胀波在燃气介质中的传播特性，结合膨胀波火炮两相流内弹道方程，给出了膨胀波在火炮膛内传播速度的求解方法。通过对膛内压降传递规律的捕捉计算和试验比对，较为准确的获得了膨胀波在膛内的传播规律，为火炮系统膛内波动效应特性的研究提供了一种可行的求解方法。关键字：膨胀波；两相流；内弹道方程；数值计算：TJ303.4CapturingCalculationAnalysisofPropagationLawaboutRarefactionWaveintheGunWANGYing-ze,ZHANGXiao-bing,YUANYan-xiong(CollegeofPowerEngineering,NUST,Nan激ng210094,China)Abstract:Toobtaintherarefactionwavepropagationlawinthegun,theconstitutionofcomplexwavesystemintherarefactionwavegunanditspropagationcharacteristicsareanalyzedaccordingtotheemissioncharacteristicsoftherarefactionwavegun.Basedonthepropagationcharacteristicsoftherarefactionwave,thesolvingmethodaboutthepropagationvelocityoftherarefactionwaveinthegunisgivenbycombingwiththetwo-phaseflowinteriorballisticequationsoftherarefactionwavegun.Bythecapturingcalculationaboutthepropagationlawofthepressuredropandthecomparingwiththeexperimentalresults,thepropagationlawoftherarefactionwaveinthegunisobtainedexactly,andthenafeasiblemethodtoresearchthewaveeffectinthegunisgiven.Keywords:rarefactionwave;two-phaseflows;interiorballisticequation;numericalcalculation在火炮发射过程中，火药的瞬态燃烧及弹丸运动引起的弹后空间的变化，导致膛内产生强压缩波和膨胀波在膛底和弹底之间往复传播的复杂波动现象。对于传统火炮膛内复杂波动现象的相关研究由于膛内流动状态的复杂性及其对整个内弹道过程的宏观效应影响较小等原因而没有深入进行，只是在分析膛内异常压力产生机理的时候对于膛内压力波的传递现象进行了相关的分析[1,2]。基于新发射机理而研制的新型低后坐武器系统一膨胀波火炮系统[3,4]，通过在炮尾安装后喷装置。1火炮膛内复杂波系组成分析如图(1a)所示，常规火炮系统发射过程中，伴随弹丸启动前后膛内燃烧流动状态的不同，膛内呈喷装置打开前后膛内流动状态的复杂变化，在膛内形成较传统火炮更加复杂的波动现象。如图(1b)所示，动现象，同时由于弹丸不断向。(a)常规火炮(b)膨胀波火炮图1火炮膛内波系分布由表2可知，药室隔断，计算区域为惯性炮尾右端由药室-身管区域扩展为扩张喷管-药室-身管区。2膨胀波在膛内传递规律的捕捉计算考虑到膛内火药燃气不断生成补充的影响，弹丸传播规律决定了膨胀波火炮系统性能的发挥。2.1膨胀波传播特性分析由气体动的传播由于固相颗粒的存在而受到制约，难，必须借助描述膛内流动规律的内弹道方程组利用压降规律间接求解。2.2膨胀波火炮两相流内弹道方程根据膨胀波火炮得内弹道方程组为：(1)表1凸轮曲线布局对比θ1/（°）θ2/（°）θ3/（°）θ4/（°）θ5/（°）θ6/（°）θ7/（°）θ8/（°）θ9/（°）θ10/（°）原布局19849571321984957132新设计72.5－70－1572.55570－－g、p分别代表气相和固相；为空隙率，为发射药燃烧生成的气体源项；为点火药生成的气体源项；为相间阻力；为发射药燃烧的化学潜能；为点火药比焓，为相间热传导量；为颗粒间应力；其具体的求解公式可参考文献[7]。2.3计算区域的划分根据膨胀波火炮发射过程分析可知，随着后喷装置的打开，炮发射过程分析可知，随着后喷装置的打开，扩张喷管与药室连通，膛内火药燃气流动区域发生变化，相应的计算求解区域也随之改变。表2凸轮曲线布局对比θ1/（°）θ2/（°）θ3/（°）θ4（/（°）原布局19849571新设计72.5－70－通过图2所示膨胀波火炮内弹道计算区域划分简图可知，在后喷装置打开前，后喷部分与药室隔断，计算区域为惯性炮尾右端到弹底之间的药室-身管部分，药室-身管...