膨胀波在火炮膛内传播规律的捕捉计算分析1.docx

膨胀波在火炮膛内传播规律的捕捉计算分析王颖泽，张小兵，袁亚雄(南京理工大学动力工程学院，南京210094)摘要：为了准确的获得膨胀波在火炮膛内的传播规律，根据膨胀波火炮的发射特点，分析了膨胀波火炮发射过程中膛内复杂波系的组成及传播特性；利用膨胀波在燃气介质中的传播特性，结合膨胀波火炮两相流内弹道方程，给出了膨胀波在火炮膛内传播速度的求解方法。通过对膛内压降传递规律的捕捉计算和试验比对，较为准确的获得了膨胀波在膛内的传播规律，为火炮系统膛内波动效应特性的研究提供了一种可行的求解方法。关键字：膨胀波；两相流；内弹道方程；数值计算中图分类号：TJ303.4CapturingCalculationAnalysisofPropagationLawaboutRarefactionWaveintheGunWANGYing-ze,ZHANGXiao-bing,YUANYan-xiong(CollegeofPowerEngineering,NUST,Nanjing210094,China)Abstract:Toobtaintherarefactionwavepropagationlawinthegun,theconstitutionofcomplexwavesystemintherarefactionwavegunanditspropagationcharacteristicsareanalyzedaccordingtotheemissioncharacteristicsoftherarefactionwavegun.Basedonthepropagationcharacteristicsoftherarefactionwave,thesolvingmethodaboutthepropagationvelocityoftherarefactionwaveinthegunisgivenbycombingwiththetwo-phaseflowinteriorballisticequationsoftherarefactionwavegun.Bythecapturingcalculationaboutthepropagationlawofthepressuredropandthecomparingwiththeexperimentalresults,thepropagationlawoftherarefactionwaveinthegunisobtainedexactly,andthenafeasiblemethodtoresearchthewaveeffectinthegunisgiven.Keywords:rarefactionwave;two-phaseflows;interiorballisticequation;numericalcalculation在火炮发射过程中，火药的瞬态燃烧及弹丸运动引起的弹后空间的变化，导致膛内产生强压缩波和膨胀波在膛底和弹底之间往复传播的复杂波动现象。对于传统火炮膛内复杂波动现象的相关研究由于膛内流动状态的复杂性及其对整个内弹道过程的宏观效应影响较小等原因而没有深入进行，只是在分析膛内异常压力产生机理的时候对于膛内压力波的传递现象进行了相关的分析”。基于新发射机理而研制的新型低后坐武器系统-膨胀波火炮系统34,通过在炮尾安装后喷装置。1火炮膛内复杂波系组成分析如图(Ia)所示，常规火炮系统发射过程中，伴随弹丸启动前后膛内燃烧流动状态的不同，膛内呈喷装置打开前后膛内流动状态的复杂变化，在膛内形成较传统火炮更加复杂的波动现象。如图(Ib)所示，动现象，同时由于弹丸不断向。(b)膨胀波火炮图1火炮膛内波系分布由表2可知，药室隔断，计算区域为惯性炮尾右端由药室-身管区域扩展为扩张喷管-药室-身管区。2膨胀波在膛内传递规律的捕捉计算考虑到膛内火药燃气不断生成补充的影响，弹丸传播规律决定了膨胀波火炮系统性能的发挥。2.1膨胀波传播特性分析由气体动的传播由于固相颗粒的存在而受到制约，难，必须借助描述膛内流动规律的内弹道方程组利用压降规律间接求解。2.2膨胀波火炮两相流内弹道方程根据膨胀波火炮得内弹道方程组为：。叫)+(*)=%+%H饵"g)+3"i)+噜"l(°)2(°)3(°)4(°)5(°)6(°)7(°)8(°)9(°)o(°)原布局19849571321984957132新设计72.5-70-1572.55570-表1凸轮曲线布局对比g、P分别代表气相和固相；e为空隙率,mc为发射药燃烧生成的气体源项；%为点火药生成的气体源项；£为相间阻力；.为发射药燃烧的化学潜能；Hgn为点火药比焰，&为相间热传导量；RP为颗粒间应力；其具体的求解公式可参考文献7。2.3 计算区域的划分根据膨胀波火炮发射过程分析可知，随着后喷装置的打开，炮发射过程分析可知，随着后喷装置的打开，扩张喷管与药室连通，膛内火药燃气流动区域发生变化，相应的计算求解区域也随之改变。表2凸轮曲线布局对比1(0)2(°)3(°)4(q)原布局19849571新设计72.5-70-通过图2所示膨胀波火炮内弹道计算区域划分简图可知，在后喷装置打开前，后喷部分与药室隔断，计算区域为惯性炮尾右端到弹底之间的药室-身管部分，药室-身管部分，药室-身管部分，断，计算区域为惯性炮尾右端药室-身管部分，当后喷装置打开后，扩张喷管与药室连通，计算区域左端面由惯性炮尾右端扩展至喷口，相应的计算区域由药室-身管区域扩展为扩张喷管-药室-身管区域。后喷装置打开前药室-身管图2计算区域划分简图2.4 数值差分格式针对上面建立的两相流内弹道方程选取具有二阶精度的MacCormack差分进行数值求解。将式改写为矢量的形式：则MacCormack差分格式为：预估计算Unj+1=UrMlx(F7"1-F/)+阳;(3)校正计算UT=t7f1-x(+1-¾1)+1(5)稳定性条件:At<CoAxu+c2.5 初始条件及边界条件处理2.5.1 初始条件处理把火炮开始射击时刻作为初始条件。此刻，气固两相速度均为零，压力为常压，温度为初始药温，密度由状态方程确定，空隙率由装填条件确定。2.5.2 惯性炮尾端部边界条件惯性炮尾运动前为固壁边界开前为运动边界，采用运动控制体法处理；喷管打开后，惯性炮尾端部边界条件失效。2.5.3 弹底边界条件弹丸运动前为固壁边界，采用镜面发射法处理；运动后为运动边界，采用运动控制体法处理。参考文献1金志明.火炮膛内压力波机理及其特征分析J南京理工大学学报,1992,61(1)：26-31JINGZhi-ming.Generationmechanismandcharacteristicanalysisofpressurewaveingun.JournalofNanjingUniversityofScienceandTechnologyJ,1992,61(l):26-31(inChinese)2金志明，袁亚雄，宋明.现代内弹道学M.北京：高等教育出版社，2003JINGZhi-ming,YUANYa-xiong,SONGMing.ModerninteriorballisticsM.Beijing:HigherEducationPublishingCompany,2003(inChinese)3 KatheE,DillonR.SonicrarefactionwavelowrecoilgunR.PaperAIAA-01-0743,2001,14 KatheE,DillonR,SopokS,etal.RarefactionwavegunpropulsionR.ADA389156,2001,45 WANGYing-ze,ZHANGXiao-bing,YUANYa-xiong.ComputationalanalysisaboutventingtimeofventdevicesofrarefactionwavegunJ.JournalofChinaOrdnance,2008,4(4):241-2466罗曼卢.气体动力学M.上海:上海交通大学出版社，1989LUOMan-luo.GasdynamicsM.Shanghai:PublishingCompanyofShanghaiJiaoTongUniversity,19897袁亚雄，张小兵.高温高压流体动力学基础M哈尔滨：哈尔滨工业大学出版社，2005YUANYa-xiong,ZHANGXiao-bing.MultiphasehydrokineticfoundationofhightemperatureandhighpressureM.Harbin:PublishingCompanyofHarbinInstituteofTechnology,2005(inChinese)8KatheE.RarefactionwavegunpropulsionD.NewYork:RensselaerPolytechnicInstitute,Troy,2002