一种燃气射流冲击液体水柱内弹道性能的仿真计算方法技术

本发明专利技术公开了一种燃气射流冲击液体水柱内弹道性能的仿真计算方法，包括如下步骤：预设仿真计算模型前提条件；基于步骤S1的前提条件建立气液物理作用过程模型；基于气液物理作用过程模型将整个内弹道过程划分为四个仿真计算阶段，所述四个仿真计算阶段分别为：燃气生成流量计算阶段、燃气流入圆柱管内状态计算阶段、圆柱管内液柱运动计算阶段和高低压室压力计算阶段；获取高低压室压力曲线。本发明专利技术可实现简洁、高效获得燃气射流冲击管内液柱平衡体现象的内弹道性能参数及推力特性，对于研究有燃气射流冲击管内液柱平衡体现象的工业产品、军事装备等系统结构优化设计具有重要意义。

【技术实现步骤摘要】
一种燃气射流冲击液体水柱内弹道性能的仿真计算方法
本专利技术涉及工业生产和军事领域，具体涉及一种燃气射流冲击液体水柱内弹道性能的仿真计算方法。
技术介绍
燃气射流冲击液柱平衡体的理论研究和试验系统中，包含了点火系统、燃气发生系统、喷管、以及圆柱管等组件，喷射过程以点火药被点燃为起点，燃气与液体全部从圆柱管中喷出后燃烧室内的压力降至大气压力为终点，期间包含了点火药的点燃、推进剂的燃烧、喷管破膜、火药燃气的流动及对液柱进行冲击后在强烈相互作用过程中共同从圆柱管中喷出等复杂而变化的过程。一般通过实验研究虽然可以准确地获得某些内弹道参数，但仍需耗费大量的时间和精力，且随着实验人员专业素质不稳定，会给实验获取数据带来不确定性误差。因此近年来，出现基于流体仿真计算来模拟该物理作用过程，从而仿真推算出内弹道参数，但由于燃气射流冲击液柱平衡体形态变化具有随机性，因此必须选择一个与实际物理作用过程、形变接近，又利于仿真计算的物理假设变化模型，才能使仿真计算合理有效。
技术实现思路
专利技术目的：为了克服现有技术中存在的不足，提供一种燃气射流冲击液体水柱内弹道性能的仿真计算方法，可准确有效的获取发射过程中各内弹道参数的变化情况，减少试验次数，节约试验成本，可广泛应用于大型装备制造、军事装备制造中涉及液体平衡发射过程的理论研究和验证模拟中，对采用液体平衡发射方式的制造装备的结构设计、内弹道过程的设计提供理论基础。技术方案：为实现上述目的，本专利技术提供一种燃气射流冲击液体水柱内弹道性能的仿真计算方法，包括如下步骤：S1：预设仿真计算模型前提条件；S2：基于步骤S1的前提条件建立气液物理作用过程模型；S3：基于气液物理作用过程模型将整个内弹道过程划分为四个仿真计算阶段，所述四个仿真计算阶段分别为：燃气生成流量计算阶段、燃气流入圆柱管内状态计算阶段、圆柱管内液柱运动计算阶段和高低压室压力计算阶段；S3：获取高低压室压力曲线。进一步的，所述步骤S1中仿真计算模型前提条件包括火药燃烧条件、燃气产生条件、燃气流动方式、平均压力模型和次要功处理方法五个条件。进一步的，所述步骤S2中气液物理作用过程模型分为气液未接触、气液冲击初始、气液作用推进、圆柱管内无液体四个阶段。进一步的，所述步骤S3中燃气生成流量的计算步骤如下：步骤S3.1-1：计算燃烧室内燃气发生量推进剂燃速方程：推进剂形状函数：式中Z为已燃厚度百分比，t为时间，u1为燃速系数，p为燃烧室内的平均压力，n为燃速指数，e1为推进剂初始弧厚的一半，Zk为推进剂分裂后碎粒全部燃完时的燃去相对厚度，如果ds为碎粒断面的内切圆半径，则Zk＝(e1+ds)/e1，ψ为火药已燃质量百分比，χ、λ、μ为推进剂形状特征量，χs、λs为推进剂在分裂时刻的形状特征量；燃烧室内燃气状态方程：式中f为推进剂火药力，ω为推进剂质量，τ＝T/Tz为燃烧室内的绝对温度与推进剂燃烧爆温的比值，η为燃烧室内燃气流入到圆柱管内的相对流出量，f1为点火药火药力，ω1为点火药质量，V0为燃烧室容积，ρp为推进剂密度，α为推进剂燃烧产生气体的余容，α1为点火药燃烧产生气体的余容；燃烧室系统内的能量方程：式中kn为推进剂燃烧产生燃气的比热比；步骤S3.1-2：计算燃烧室喷口燃气流出量燃烧室内燃气流出的相对流出量方程：式中为流量损耗系数，St为喷管喉部面积，pb为圆柱管内的平均压力。理想假设燃烧室喷口喷出的所有燃气均无损失喷入圆柱管口。进一步的，所述步骤S3中燃气流入圆柱管内状态的计算步骤如下：步骤S3.2-1：计算燃气管内状态圆柱管内燃气状态方程：式中τ1＝T1/Tz为圆柱管内的绝对温度与推进剂燃烧爆温的比值，η1为当液体水柱全部从圆柱管内喷出后圆柱管内燃气流入到外界大气环境的相对流出量，为次要功计算系数，m为液体水柱的质量，v为液体水柱的运动速度，Vψ为圆柱管内喷管出口与液体柱之间的初始自由容积，V1为燃气冲击液柱引起的自由容积的增加；圆柱管系统内的能量方程：步骤S3.2-2：计算燃气管内流出量圆柱管内燃气流出的相对流出量方程：式中Sd为圆柱管出口横截面积，Pa为外界环境大气压力。进一步的，所述步骤S3中圆柱管内液柱运动的计算过程如下：液柱运动方程：式中l为液柱运动位移，Pf为液柱运动阻压，Sn为燃气对液柱的有效推力面积，为液柱质量修正系数；空腔容积变化方程：式中r′为空腔顶部到内管口垂直距离，r为内管半径。进一步的，所述步骤S3中采用内弹道方程组对高低压室压力进行计算。进一步的，所述步骤S3中采用四阶龙格-库塔法对燃气射流冲击液柱过程的内弹道方程组进行求解，其具体计算方法如下：对于一阶微分方程组：设定变量i＝1,2,…,n，则上述微分方程组的初值为yi(t)＝yi0，进而四阶龙格—库塔法的公式可写成：在公式的右端h＝tk+1-tk为时间步长，括号内的各阶变量的赋值公式如下所示：通过以上四阶龙格-库塔法逐步迭代，得到一阶微分方程组式(12)的数值解，对于燃气射流冲击液柱平衡体的内弹道方程组而言，以时间t为自变量，当各变量被赋予初值后，调用四阶龙格-库塔法就可积分一个时间步长，重复逐步积分下去即可得到最终的计算结果。此计算结果为每增加一个时间步长后该时间点的火药燃烧百分比、高压室压力、低压室压力、液柱运动长度、液柱运动速度等数据，最后将所得到的数据通过作图即可获得在整个燃气射流作用液柱过程中各参量尤其是高压室压力、低压室压力随时间变化的曲线。本专利技术方法是针对燃气射流冲击圆柱管内液柱的物理作用过程中，对液柱在受到射流冲击后形状变形过程提出一种新的理想变形假设模型。即假设射流冲击管内液柱后，液柱呈现空腔，且空腔从早期呈现椭球形状，并假设此后空腔以头部为半球体，后面衔接为等直径的圆柱体形式在液体水柱内发展，液体水柱在燃气射流的推动下也不断向圆柱管外运动，当液体水柱靠近喷管出口的一端运动了与圆柱管长度相等的距离时，假设液体已全部从圆柱管中喷出。同时，假设该物理作用过程中，没有液体水汽化作用，液柱质量保持不变。有益效果：本专利技术与现有技术相比，具备如下优点：1、相比于将液柱平衡体假设为形状不变的类固体平衡体而言，本专利技术设计的物理假设模型更接近于实际工程应用情况，计算所得到的结果与试验结果更加符合，从而提高了仿真计算精度，对于研究有燃气射流冲击管内液柱平衡体现象的工业产品、军事装备等系统结构优化设计具有重要意义。2、可准确有效的获取发射过程中各内弹道参数的变化情况，减少试验次数，节约试验成本。附图说明图1为液体水柱的运动以及空腔的变化过程示意图，其中(a)为冲击早期液体水柱结构，(b)为冲击末期本文档来自技高网...

【技术保护点】
1.一种燃气射流冲击液体水柱内弹道性能的仿真计算方法，其特征在于，包括如下步骤：/nS1：预设仿真计算模型前提条件；/nS2：基于步骤S1的前提条件建立气液物理作用过程模型；/nS3：基于气液物理作用过程模型将整个内弹道过程划分为四个仿真计算阶段，所述四个仿真计算阶段分别为：燃气生成流量计算阶段、燃气流入圆柱管内状态计算阶段、圆柱管内液柱运动计算阶段和高低压室压力计算阶段；/nS3：获取高低压室压力曲线。/n

【技术特征摘要】
1.一种燃气射流冲击液体水柱内弹道性能的仿真计算方法，其特征在于，包括如下步骤：
S1：预设仿真计算模型前提条件；
S2：基于步骤S1的前提条件建立气液物理作用过程模型；
S3：基于气液物理作用过程模型将整个内弹道过程划分为四个仿真计算阶段，所述四个仿真计算阶段分别为：燃气生成流量计算阶段、燃气流入圆柱管内状态计算阶段、圆柱管内液柱运动计算阶段和高低压室压力计算阶段；
S3：获取高低压室压力曲线。


2.根据权利要求1所述的一种燃气射流冲击液体水柱内弹道性能的仿真计算方法，其特征在于，所述步骤S1中仿真计算模型前提条件包括火药燃烧条件、燃气产生条件、燃气流动方式、平均压力模型和次要功处理方法五个条件。


3.根据权利要求1所述的一种燃气射流冲击液体水柱内弹道性能的仿真计算方法，其特征在于，所述步骤S2中气液物理作用过程模型分为气液未接触、气液冲击初始、气液作用推进、圆柱管内无液体四个阶段。


4.根据权利要求1所述的一种燃气射流冲击液体水柱内弹道性能的仿真计算方法，其特征在于，所述步骤S3中燃气生成流量的计算步骤如下：
步骤S3.1-1：计算燃烧室内燃气发生量
推进剂燃速方程：



推进剂形状函数：



式中Z为已燃厚度百分比，t为时间，u1为燃速系数，p为燃烧室内的平均压力，n为燃速指数，e1为推进剂初始弧厚的一半，Zk为推进剂分裂后碎粒全部燃完时的燃去相对厚度，如果ds为碎粒断面的内切圆半径，则Zk＝(e1+ds)/e1，ψ为火药已燃质量百分比，χ、λ、μ为推进剂形状特征量，χs、λs为推进剂在分裂时刻的形状特征量；
燃烧室内燃气状态方程：



式中f为推进剂火药力，ω为推进剂质量，τ＝T/Tz为燃烧室内的绝对温度与推进剂燃烧爆温的比值，η为燃烧室内燃气流入到圆柱管内的相对流出量，f1为点火药火药力，ω1为点火药质量，V0为燃烧室容积，ρp为推进剂密度，α为推进剂燃烧产生气体的余容，α1为点火药燃烧产生气体的余容；
燃烧室系统内的能量方程：



式中kn为推进剂燃烧产生燃气的比热比；
步骤S3.1-2：计算燃烧室喷口燃气流出量
燃烧室内燃气流出的相对流出量方程：



式中为流量损耗系数，St为喷管喉部面积，pb为圆柱管内的平均压力。理想假设燃烧室喷口喷出的所有燃气均无损...

【专利技术属性】
技术研发人员：何淼，钱佳红，
申请(专利权)人：南京信息职业技术学院，
类型：发明
国别省市：江苏;32