膨胀波在火炮膛内传播规律

《膨胀波在火炮膛内传播规律》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、膨胀波在火炮膛内传播规律的捕捉计算分析王颖泽，张小兵，袁亚雄(南京理工大学动力工程学院，南京210094)摘要：为了准确的获得膨胀波在火炮膛内的传播规律，根据膨胀波火炮的发射特点，分析了膨胀波火炮发射过程中膛内复杂波系的组成及传播特性；利用膨胀波在燃气介质中的传播特性，结合膨胀波火炮两相流内弹道方程，给出了膨胀波在火炮膛内传播速度的求解方法。通过对膛内压降传递规律的捕捉计算和试验比对，较为准确的获得了膨胀波在膛内的传播规律，为火炮系统膛内波动效应特性的研究提供了一种可行的求解方法。关键字：膨胀波；两相流；内弹道方程；数值计算在火炮发射过程中，火药的瞬态燃烧及弹丸运动引起的弹后

2、空间的变化，导致膛内产生强压缩波和膨胀波在膛底和弹底之间往复传播的复杂波动现象。对于传统火炮膛内复杂波动现象的相关研究由于膛内流动状态的复杂性及其对整个内弹道过程的宏观效应影响较小等原因而没有深入进行，只是在分析膛内异常压力产生机理的时候对于膛内压力波的传递现象进行了相关的分析[1,2]。1火炮膛内复杂波系组成分析如图(1a)所示，常规火炮系统发射过程中，伴随弹丸启动前后膛内燃烧流动状态的不同，膛内呈喷装置打开前后膛内流动状态的复杂变化，在膛内形成较传统火炮更加复杂的波动现象。如图(1b)所示，动现象，同时由于弹丸不断向。(a)常规火炮(b)膨胀波火炮图1火炮膛内波系分布由表

3、2可知，药室隔断，计算区域为惯性炮尾右端由药室-身管区域扩展为扩张喷管-药室-身管区。2膨胀波在膛内传递规律的捕捉计算考虑到膛内火药燃气不断生成补充的影响，弹丸传播规律决定了膨胀波火炮系统性能的发挥。2.1膨胀波传播特性分析由气体动的传播由于固相颗粒的存在而受到制约，难，必须借助描述膛内流动规律的内弹道方程组利用压降规律间接求解。2.2膨胀波火炮两相流内弹道方程根据膨胀波火炮得内弹道方程组为：(1)表1凸轮曲线布局对比θ1/（°）θ2/（°）θ3/（°）θ4/（°）θ5/（°）θ6/（°）θ7/（°）θ8/（°）θ9/（°）θ10/（°）原布局1984957132198495

4、7132新设计72.5－70－1572.55570－－g、p分别代表气相和固相；为空隙率，为发射药燃烧生成的气体源项；为点火药生成的气体源项；为相间阻力；为发射药燃烧的化学潜能；为点火药比焓，为相间热传导量；为颗粒间应力；其具体的求解公式可参考文献[7]。2.3计算区域的划分根据膨胀波火炮发射过程分析可知，随着后喷装置的打开，炮发射过程分析可知，随着后喷装置的打开，扩张喷管与药室连通，膛内火药燃气流动区域发生变化，相应的计算求解区域也随之改变。表2凸轮曲线布局对比θ1/（°）θ2/（°）θ3/（°）θ4（/（°）原布局19849571新设计72.5－70－通过图2所示膨胀波火

5、炮内弹道计算区域划分简图可知，在后喷装置打开前，后喷部分与药室隔断，计算区域为惯性炮尾右端到弹底之间的药室-身管部分，药室-身管部分，药室-身管部分，断，计算区域为惯性炮尾右端药室-身管部分，当后喷装置打开后，扩张喷管与药室连通，计算区域左端面由惯性炮尾右端扩展至喷口，相应的计算区域由药室-身管区域扩展为扩张喷管-药室-身管区域。图2计算区域划分简图2.4数值差分格式针对上面建立的两相流内弹道方程选取具有二阶精度的MacCormack差分进行数值求解。将式(1)改写为矢量的形式：(2)则MacCormack差分格式为：预估计算(3)校正计算(4)(5)稳定性条件：2.5初始条

6、件及边界条件处理2.5.1初始条件处理把火炮开始射击时刻作为初始条件。此刻，气固两相速度均为零，压力为常压，温度为初始药温，密度由状态方程确定，空隙率由装填条件确定。2.5.2惯性炮尾端部边界条件惯性炮尾运动前为固壁边界开前为运动边界，采用运动控制体法处理；喷管打开后，惯性炮尾端部边界条件失效。2.5.3弹底边界条件弹丸运动前为固壁边界，采用镜面发射法处理；运动后为运动边界，采用运动控制体法处理。后续省略参考文献[1]金志明.火炮膛内压力波机理及其特征分析[J].南京理工大学学报，1992，61(1)：26-31[2]金志明，袁亚雄，宋明.现代内弹道学[M].北京：高等教育出

7、版社，2003[3]KatheE,DillonR.Sonicrarefactionwavelowrecoilgun[R].PaperAIAA-01-0743,2001,1[4]KatheE,DillonR,SopokS,etal.Rarefactionwavegunpropulsion[R].ADA389156,2001,4[5]WANGYing-ze,ZHANGXiao-bing,YUANYa-xiong.Computationalanalysisaboutventingtimeofventdevic