CN110289053A - 发射药燃烧产生等离子体规律研究方法 - Google Patents
发射药燃烧产生等离子体规律研究方法 Download PDFInfo
- Publication number
- CN110289053A CN110289053A CN201910559092.XA CN201910559092A CN110289053A CN 110289053 A CN110289053 A CN 110289053A CN 201910559092 A CN201910559092 A CN 201910559092A CN 110289053 A CN110289053 A CN 110289053A
- Authority
- CN
- China
- Prior art keywords
- equation
- plasma
- combustion
- temperature
- powder
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Chemical & Material Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了发射药燃烧产生等离子体规律研究方法,包括以下步骤:S1:发射药燃烧生成等离子体参数的计算;S2:用平衡常数法计算燃气的化学成份;S3:结果分析,根据上面的计算方法,将方程组在MATLAB中编程,可以计算出等离子体参数之间的关系;S4:物理过程分析;S5:模型构建;S6:模型验证得出结论。本发明的研究方法能够得出电子密度和电离度都随着温度的增加而增大;相同温度下压力越大,电子密度越大,电离度越小,装药量和药室容积的不同改变燃气温度变化速率,相应的影响了等离子体电子密度的变化速率,火药力的不同改变了燃气最高温度,相应的改变了等离子体电子密度的最大值。
Description
技术领域
本发明涉及等离子体领域,具体为发射药燃烧产生等离子体规律研究方法。
背景技术
等离子体经过几十年的发展,已经发展成为包含天体等离子体、核聚变等离子体、低温等离子体等分支的独立学科,等离子体也被广泛应用于受控核聚变、磁流体发电、材料表面工程等领域,然而,在火炮方面,仅有外加等离子体被应用于电热化学炮的研究,对于发射药燃烧产生的等离子体用于提高火炮性能的研究还未涉及;发射药燃烧生成等离子体的规律,是研究利用发射药生成等离子体提高火炮性能的基础。
发明内容
本发明的目的在于提供发射药燃烧产生等离子体规律研究方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:发射药燃烧产生等离子体规律研究方法,包括以下步骤:
S1:发射药燃烧产生等离子体规律研究方法,对于化学反应中反应物A和B生成产物C和D,a、b、c、d分别是相应物质的物质的量;
其正向和逆向反应速度分别为:
V1=K1[A]a[B]b,V2=K2[C]c[D]d 2
其中,[A]、[B]、[C]、[D]分别表示各物质的浓度,当反应达到平衡状态时,即正逆反应的速率达到相等,即:
K1[A]a[B]b=K2[C]c[D]d 3
推导后可以得到该反应的化学平衡常数K:
根据质量守恒定律,任何参与反应的元素在反应中任何时刻的原子数目是不变的;
化学平衡常数法是以化学反应平衡和质量守恒定律为基础,通过联立各个反应的方程,组成方程组,求解各个组分,它的一个优点是能够直观的反映出化学反应的真实状态,另一个优点是方便调节,可以根据目标需要,改变反应体系中的不同反应;
S2:用平衡常数法计算燃气的化学成份,根据平衡常数法原理,做出几点假设:燃烧产物混合均匀;燃气组分处于化学平衡状态下;燃气等离子体服从理想气体状态方程,为了获得更高的等离子体密度,在发射药中加入碳酸钾,发射药和添加剂中一共含有C、H、O、N、K五种元素,对于燃烧的化学反应过程,可以用以下15个方程表示:
以上公式中Ki是各成份的化学平衡常数,K6K7K8K9K10K11的值可以通过以下公式计算得到;
发射药燃气生成等离子体的燃气成份计算中,一共考虑了21种成份,分别为CO2、CO、H2O、OH、O2、O、H2、H、N2、NO、N、N2O、NO2、K2CO3、K2O、KO、KOH、K、K+、OH-、e电子,按照上面的顺序相应的摩尔分数xi的下表依次为1,2,…21;
根据四个质量守恒定律,可以得到以下四个方程:
上面公式当中的MC MH MO MN MK分别表示反应物当中五种元素C、H、O、N、K的质量,它们的数值可以通过发射药当中添加剂的比例计算得到,如x千克发射药当中添加了质量百分比为n的种子,则
MC=12×(25.5243+1000×n/((1-x)×138)) 31
MH=33.1277 32
MO=16×(33.0757+3×1000×n/((1-x)×138)) 33
MN=14×9.3395 34
MK=78×1000×n/((1-x)×138) 35
根据道尔顿分压定律:
由于xi是摩尔分数,所以有:
pi=xip 37
由式37可以获得:
另外,对于带电粒子,由于电中性,假设第i组分带电量qi,则有:
这样,获得21个方程,组成了可以求解21种燃烧产物组分的闭合方程组,通过计算可以求解相关组分的参数;
燃气当中发生部分电离,其电子密度的计算公式为:
S3:结果分析,根据上面的计算方法,将方程组在MATLAB中编程,可以计算出等离子体参数之间的关系;
S4:物理过程分析,经典内弹道学是以热力学为理论依据,研究高压状态下各个火炮内弹道参数平均值的理论,经典内弹道学的数学模型是一阶常微分方程,由火药形状函数方程、能量方程、火药燃速方程、弹丸运动速度方程和弹丸速度与行程方程五个方程组成的方程组;
根据高压状态下火药燃烧情况的特点,做出以下假设:
假设无论主装药还是点火药都符合几何燃烧规律;
火药燃烧采用燃速指数公式,即
假定发射药燃烧和弹丸运动都是在平均压力下进行的;
射击时热散失比较复杂,而且难以描述,通常使用减小火药力或增大比热的方法进行修正,实际上是减少发射药量;
所有的次要功用系数来计算,由于射击过程中次要功主要是与动能有关,所有的次要功都可以用系数来加以考虑;
以弹丸挤进弹带的压力作为弹丸启动的条件;
膛内燃气符合诺贝尔——阿贝尔状态方程;
弹带密封良好,不存在漏气现象;
对于某种特定的火药来说,它的爆温是一定的,这和火药的性质有关,是一个常量,因此,定容状态下的气体状态方程如式41所示:
其中,Vψ为自由容积,随压力一起变化,使用火药力和装填密度表示的定容状态下的气体状态方程如式42所示:
从公式当中可以看出,压力随着火药燃烧百分比而逐渐变化,但是,燃烧温度一直是火药的爆温;
为了增加燃烧产物的热电离,本发明中在发射药里添加少量的碳酸钾,根据热电离理论,做出如下假设:
碳酸钾在高温环境中完全分解;
假定燃烧产物是均匀的;
不考虑热电离过程中出现概率极小的二次电离或者多重电离的情况;
S5:模型构建,火药形状函数方程如式43所示:
其中,χ、λ为火药***之前的形状特征量;Zk为燃烧结束时火药相对已燃厚度;χs、λs为火药***后的形状特征量;
燃烧速度方程如式44所示:
根据火药燃烧内弹道的基本假设,同时考虑次要功的影响,可由牛顿定律获得:
根据能量守恒定律,得到火药燃烧的能量方程如式47所示:
根据以上方程得到高压状态下内弹道方程组如式48所示:
由于火药燃气不断推动弹丸做功而损失能量,所以燃气温度不断降低,火药气体温度与体积的变化关系如式49所示:
p(Vψ+Sl)=ωψRT 49
通常以弹丸行程表示:
Sp(l+lψ)=ωψRT 50
温度是与压力p、弹丸行程l、火药已燃百分比ψ有关的函数,这些数据可以通过内弹道方程组求得;
通过燃气温度方程,可以获得火药燃气的燃温,由于火药燃气中生成等离子体的方式是热电离,所以可以应用萨哈方程计算燃气中电子密度;
高压状态下生成等离子体密度方程组由三部分组成:内弹道方程组、燃气温度方程和电子密度方程,仿真过程分为三步,首先根据内弹道方程计算出燃气压力、弹丸行程等与火药燃烧质量的关系,然后根据内弹道的计算结果,通过燃气温度方程计算出火药燃气的温度,最后再根据电子密度方程计算等离子体中电子密度;
S6:模型验证,由数值仿真获得的数据和试验数据进行对比,在弹丸质量m=0.39kg,装药量ω=0.136kg,药室容积V0=1.32×10-4m3的情况下,求得膛压最大值在374MPa,弹丸初速898m/s,与试验结果基本吻合,图像基本相近,验证了所建模型的正确性。
优选的,所述S1中参数aik ajk分别表示1mol气相物质和固相物质中某元素的原子数目,ngi nsj分别表示单位物质中气相组分和固相组分的摩尔数,Nk则表示某种元素的原子总数。
优选的,所述S2中K1K2K3K4K5K12K13K14K15的数值均为已知常数。
优选的,所述S3中电子密度和温度是正相关的关系,随着温度的增加,电子密度也在增加,而且当温度一样的时候,压力越大电子密度越大,这是因为随着温度的增加,电子获得的能量越多,电子越容易挣脱原子的束缚而成为自由电子。
优选的,所述S4中S是弹丸最大横截面积;p是火药燃气平均压力;m是弹丸的质量;l是弹丸行程;v是弹丸速度;是次要功系数。
优选的,所述S4中lψ是药室自由容积缩径长,ω是装药质量,θ是火药热力参数,ψ是火药已燃百分数,f是火药力,Δ是装填密度,ρp是火药密度,α是余容。
优选的,所述S5中T为热力学温度,Ei为离子的电离电位,go为原子基态的统计权重,gi为离子基态的统计权重,me为电子质量,ne为电子密度,ni为离子密度,n0为原子密度,对于碱金属来说,的值约为1,其他气体一般为2。
优选的,所述S5中的等离子体中电子密度
与现有技术相比,本发明的有益效果是:电子密度和电离度都随着温度的增加而增大;相同温度下压力越大,电子密度越大,电离度越小,装药量和药室容积的不同改变燃气温度变化速率,相应的影响了等离子体电子密度的变化速率,火药力的不同改变了燃气最高温度,相应的改变了等离子体电子密度的最大值。
附图说明
图1为本发明的研究方法的流程图;
表1为本发明发射药的配方组成;
图2为本发明的电子密度与温度的关系图;
图3为本发明的内弹道试验曲线与仿真曲线对比图;
图4为本发明的膛压时间曲线图;
图5为本发明的弹丸速度时间曲线图;
图6为本发明的膛压位移曲线图;
图7为本发明的弹丸速度位移曲线。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
请参照图1,本发明提供一种技术方案:发射药燃烧产生等离子体规律研究方法,包括以下步骤:
S1:如表1所示,发射药燃烧产生等离子体规律研究方法,对于化学反应中反应物A和B生成产物C和D,a、b、c、d分别是相应物质的物质的量;
其正向和逆向反应速度分别为:
V1=K1[A]a[B]b,V2=K2[C]c[D]d 2
其中,[A]、[B]、[C]、[D]分别表示各物质的浓度,当反应达到平衡状态时,即正逆反应的速率达到相等,即:
K1[A]a[B]b=K2[C]c[D]d 3
推导后可以得到该反应的化学平衡常数K:
根据质量守恒定律,任何参与反应的元素在反应中任何时刻的原子数目是不变的;
化学平衡常数法是以化学反应平衡和质量守恒定律为基础,通过联立各个反应的方程,组成方程组,求解各个组分,它的一个优点是能够直观的反映出化学反应的真实状态,另一个优点是方便调节,可以根据目标需要,改变反应体系中的不同反应;
参数aik ajk分别表示1mol气相物质和固相物质中某元素的原子数目,ngi nsj分别表示单位物质中气相组分和固相组分的摩尔数,Nk则表示某种元素的原子总数;
S2:如图2-4所示,用平衡常数法计算燃气的化学成份,根据平衡常数法原理,做出几点假设:燃烧产物混合均匀;燃气组分处于化学平衡状态下;燃气等离子体服从理想气体状态方程,为了获得更高的等离子体密度,在发射药中加入碳酸钾,发射药和添加剂中一共含有C、H、O、N、K五种元素,对于燃烧的化学反应过程,可以用以下15个方程表示:
以上公式中Ki是各成份的化学平衡常数,K2K3K4K5K12K13K1K14K15的数值均为已知常数,K6K7K8K9K10K11的值可以通过以下公式计算得到;
发射药燃气生成等离子体的燃气成份计算中,一共考虑了21种成份,分别为CO2、CO、H2O、OH、O2、O、H2、H、N2、NO、N、N2O、NO2、K2CO3、K2O、KO、KOH、K、K+、OH-、e电子,按照上面的顺序相应的摩尔分数xi的下表依次为1,2,…21;
根据四个质量守恒定律,可以得到以下四个方程:
上面公式当中的MC MH MO MN MK分别表示反应物当中五种元素C、H、O、N、K的质量,它们的数值可以通过发射药当中添加剂的比例计算得到,如x千克发射药当中添加了质量百分比为n的种子,则
MC=12×(25.5243+1000×n/((1-x)×138)) 31
MH=33.1277 32
MO=16×(33.0757+3×1000×n/((1-x)×138)) 33
MN=14×9.3395 34
MK=78×1000×n/((1-x)×138) 35
根据道尔顿分压定律:
由于xi是摩尔分数,所以有:
pi=xip 37
由式37可以获得:
另外,对于带电粒子,由于电中性,假设第i组分带电量qi,则有:
这样,获得21个方程,组成了可以求解21种燃烧产物组分的闭合方程组,通过计算可以求解相关组分的参数;
燃气当中发生部分电离,其电子密度的计算公式为:
S3:结果分析,根据上面的计算方法,将方程组在MATLAB中编程,可以计算出等离子体参数之间的关系,电子密度和温度是正相关的关系,随着温度的增加,电子密度也在增加,而且当温度一样的时候,压力越大电子密度越大,这是因为随着温度的增加,电子获得的能量越多,电子越容易挣脱原子的束缚而成为自由电子;
S4:如图5所示,物理过程分析,经典内弹道学是以热力学为理论依据,研究高压状态下各个火炮内弹道参数平均值的理论,经典内弹道学的数学模型是一阶常微分方程,由火药形状函数方程、能量方程、火药燃速方程、弹丸运动速度方程和弹丸速度与行程方程五个方程组成的方程组;
根据高压状态下火药燃烧情况的特点,做出以下假设:
假设无论主装药还是点火药都符合几何燃烧规律;
火药燃烧采用燃速指数公式,即
假定发射药燃烧和弹丸运动都是在平均压力下进行的;
射击时热散失比较复杂,而且难以描述,通常使用减小火药力或增大比热的方法进行修正,实际上是减少发射药量;
所有的次要功用系数来计算,由于射击过程中次要功主要是与动能有关,所有的次要功都可以用系数来加以考虑;
以弹丸挤进弹带的压力作为弹丸启动的条件;
膛内燃气符合诺贝尔——阿贝尔状态方程;
弹带密封良好,不存在漏气现象;
对于某种特定的火药来说,它的爆温是一定的,这和火药的性质有关,是一个常量,因此,定容状态下的气体状态方程如式41所示:
其中,Vψ为自由容积,随压力一起变化,使用火药力和装填密度表示的定容状态下的气体状态方程如式42所示:
从公式当中可以看出,压力随着火药燃烧百分比而逐渐变化,但是,燃烧温度一直是火药的爆温;
为了增加燃烧产物的热电离,本发明中在发射药里添加少量的碳酸钾,根据热电离理论,做出如下假设:
碳酸钾在高温环境中完全分解;
假定燃烧产物是均匀的;
不考虑热电离过程中出现概率极小的二次电离或者多重电离的情况;
S是弹丸最大横截面积;p是火药燃气平均压力;m是弹丸的质量;l是弹丸行程;v是弹丸速度;是次要功系数,lψ是药室自由容积缩径长,ω是装药质量,θ是火药热力参数,ψ是火药已燃百分数,f是火药力,Δ是装填密度,ρp是火药密度,α是余容;
S5:如图6-7所示,模型构建,火药形状函数方程如式43所示:
其中,χ、λ为火药***之前的形状特征量;Zk为燃烧结束时火药相对已燃厚度;χs、λs为火药***后的形状特征量;
燃烧速度方程如式44所示:
根据火药燃烧内弹道的基本假设,同时考虑次要功的影响,可由牛顿定律获得:
根据能量守恒定律,得到火药燃烧的能量方程如式47所示:
根据以上方程得到高压状态下内弹道方程组如式48所示:
由于火药燃气不断推动弹丸做功而损失能量,所以燃气温度不断降低,火药气体温度与体积的变化关系如式49所示:
p(Vψ+Sl)=ωψRT 49
通常以弹丸行程表示:
Sp(l+lψ)=ωψRT 50
温度是与压力p、弹丸行程l、火药已燃百分比ψ有关的函数,这些数据可以通过内弹道方程组求得;
通过燃气温度方程,可以获得火药燃气的燃温,由于火药燃气中生成等离子体的方式是热电离,所以可以应用萨哈方程计算燃气中电子密度;
高压状态下生成等离子体密度方程组由三部分组成:内弹道方程组、燃气温度方程和电子密度方程,仿真过程分为三步,首先根据内弹道方程计算出燃气压力、弹丸行程等与火药燃烧质量的关系,然后根据内弹道的计算结果,通过燃气温度方程计算出火药燃气的温度,最后再根据电子密度方程计算等离子体中电子密度,T为热力学温度,Ei为离子的电离电位,go为原子基态的统计权重,gi为离子基态的统计权重,me为电子质量,ne为电子密度,ni为离子密度,n0为原子密度,对于碱金属来说,的值约为1,其他气体一般为2;
S6:模型验证,由数值仿真获得的数据和试验数据进行对比,在弹丸质量m=0.39kg,装药量ω=0.136kg,药室容积V0=1.32×10-4m3的情况下,求得膛压最大值在374MPa,弹丸初速898m/s,与试验结果基本吻合,图像基本相近,验证了所建模型的正确性,等离子体中电子密度为
通过对上述实施例进行实验研究,本发明研究方法能够得出电子密度和电离度都随着温度的增加而增大;相同温度下压力越大,电子密度越大,电离度越小,装药量和药室容积的不同改变燃气温度变化速率,相应的影响了等离子体电子密度的变化速率,火药力的不同改变了燃气最高温度,相应的改变了等离子体电子密度的最大值。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (8)
1.发射药燃烧产生等离子体规律研究方法,其特征在于,包括以下步骤:
S1:发射药燃烧生成等离子体参数的计算,对于化学反应中反应物A和B生成产物C和D,a、b、c、d分别是相应物质的物质的量,
其正向和逆向反应速度分别为:
V1=K1[A]a[B]b,V2=K2[C]c[D]d 2
其中,[A]、[B]、[C]、[D]分别表示各物质的浓度,当反应达到平衡状态时,即正逆反应的速率达到相等,即:
K1[A]a[B]b=K2[C]c[D]d 3
推导后可以得到该反应的化学平衡常数K:
根据质量守恒定律,任何参与反应的元素在反应中任何时刻的原子数目是不变的,
化学平衡常数法是以化学反应平衡和质量守恒定律为基础,通过联立各个反应的方程,组成方程组,求解各个组分,它的一个优点是能够直观的反映出化学反应的真实状态,另一个优点是方便调节,可以根据目标需要,改变反应体系中的不同反应;
S2:用平衡常数法计算燃气的化学成份,根据平衡常数法原理,做出几点假设:燃烧产物混合均匀;燃气组分处于化学平衡状态下;燃气等离子体服从理想气体状态方程,为了获得更高的等离子体密度,在发射药中加入碳酸钾,发射药和添加剂中一共含有C、H、O、N、K五种元素,对于燃烧的化学反应过程,可以用以下15个方程表示:
以上公式中Ki是各成份的化学平衡常数,K6 K7 K8 K9 K10 K11的值可以通过以下公式计算得到;
发射药燃气生成等离子体的燃气成份计算中,一共考虑了21种成份,分别为CO2、CO、H2O、OH、O2、O、H2、H、N2、NO、N、N2O、NO2、K2CO3、K2O、KO、KOH、K、K+、OH-、e电子,按照上面的顺序相应的摩尔分数xi的下表依次为1,2,…21;
根据四个质量守恒定律,可以得到以下四个方程:
上面公式当中的MC MH MO MN MK分别表示反应物当中五种元素C、H、O、N、K的质量,它们的数值可以通过发射药当中添加剂的比例计算得到,如x千克发射药当中添加了质量百分比为n的种子,则
MC=12×(25.5243+1000×n/((1-x)×138)) 31
MH=33.1277 32
MO=16×(33.0757+3×1000×n/((1-x)×138)) 33
MN=14×9.3395 34
MK=78×1000×n/((1-x)×138) 35
根据道尔顿分压定律:
由于xi是摩尔分数,所以有:
pi=xip 37
由式37可以获得:
另外,对于带电粒子,由于电中性,假设第i组分带电量qi,则有:
这样,获得21个方程,组成了可以求解21种燃烧产物组分的闭合方程组,通过计算可以求解相关组分的参数;
燃气当中发生部分电离,其电子密度的计算公式为:
;
S3:结果分析,根据上面的计算方法,将方程组在MATLAB中编程,可以计算出等离子体参数之间的关系;
S4:物理过程分析,经典内弹道学是以热力学为理论依据,研究高压状态下各个火炮内弹道参数平均值的理论,经典内弹道学的数学模型是一阶常微分方程,由火药形状函数方程、能量方程、火药燃速方程、弹丸运动速度方程和弹丸速度与行程方程五个方程组成的方程组;
根据高压状态下火药燃烧情况的特点,做出以下假设:
假设无论主装药还是点火药都符合几何燃烧规律;
火药燃烧采用燃速指数公式,即
假定发射药燃烧和弹丸运动都是在平均压力下进行的;
射击时热散失比较复杂,而且难以描述,通常使用减小火药力或增大比热的方法进行修正,实际上是减少发射药量;
所有的次要功用系数来计算,由于射击过程中次要功主要是与动能有关,所有的次要功都可以用系数来加以考虑;
以弹丸挤进弹带的压力作为弹丸启动的条件;
膛内燃气符合诺贝尔——阿贝尔状态方程;
弹带密封良好,不存在漏气现象;
对于某种特定的火药来说,它的爆温是一定的,这和火药的性质有关,是一个常量,因此,定容状态下的气体状态方程如式41所示:
其中,Vψ为自由容积,随压力一起变化,使用火药力和装填密度表示的定容状态下的气体状态方程如式42所示:
从公式当中可以看出,压力随着火药燃烧百分比而逐渐变化,但是,燃烧温度一直是火药的爆温;
为了增加燃烧产物的热电离,本发明中在发射药里添加少量的碳酸钾,根据热电离理论,做出如下假设:
碳酸钾在高温环境中完全分解;
假定燃烧产物是均匀的;
不考虑热电离过程中出现概率极小的二次电离或者多重电离的情况;
S5:模型构建,火药形状函数方程如式43所示:
其中,χ、λ为火药***之前的形状特征量;Zk为燃烧结束时火药相对已燃厚度;χs、λs为火药***后的形状特征量;
燃烧速度方程如式44所示:
根据火药燃烧内弹道的基本假设,同时考虑次要功的影响,可由牛顿定律获得:
根据能量守恒定律,得到火药燃烧的能量方程如式47所示:
根据以上方程得到高压状态下内弹道方程组如式48所示:
由于火药燃气不断推动弹丸做功而损失能量,所以燃气温度不断降低,火药气体温度与体积的变化关系如式49所示:
p(Vψ+Sl)=ωψRT 49
通常以弹丸行程表示:
Sp(l+lψ)=ωψRT 50
温度是与压力p、弹丸行程l、火药已燃百分比ψ有关的函数,这些数据可以通过内弹道方程组求得;
通过燃气温度方程,可以获得火药燃气的燃温,由于火药燃气中生成等离子体的方式是热电离,所以可以应用萨哈方程计算燃气中电子密度;
高压状态下生成等离子体密度方程组由三部分组成:内弹道方程组、燃气温度方程和电子密度方程,仿真过程分为三步,首先根据内弹道方程计算出燃气压力、弹丸行程等与火药燃烧质量的关系,然后根据内弹道的计算结果,通过燃气温度方程计算出火药燃气的温度,最后再根据电子密度方程计算等离子体中电子密度;
S6:模型验证,由数值仿真获得的数据和试验数据进行对比,在弹丸质量m=0.39kg,装药量ω=0.136kg,药室容积V0=1.32×10-4m3的情况下,求得膛压最大值在374MPa,弹丸初速898m/s,与试验结果基本吻合,图像基本相近,验证了所建模型的正确性。
2.根据如权利要求1所述的发射药燃烧产生等离子体规律研究方法,其特征在于:所述S1中参数aik ajk分别表示1mol气相物质和固相物质中某元素的原子数目,ngi nsj分别表示单位物质中气相组分和固相组分的摩尔数,Nk则表示某种元素的原子总数。
3.根据如权利要求1所述的发射药燃烧产生等离子体规律研究方法,其特征在于:所述S2中K1 K2 K3 K4 K5 K12 K13K14 K15的数值均为已知常数。
4.根据如权利要求1所述的发射药燃烧产生等离子体规律研究方法,其特征在于:所述S3中电子密度和温度是正相关的关系,随着温度的增加,电子密度也在增加,而且当温度一样的时候,压力越大电子密度越大,这是因为随着温度的增加,电子获得的能量越多,电子越容易挣脱原子的束缚而成为自由电子。
5.根据如权利要求1所述的发射药燃烧产生等离子体规律研究方法,其特征在于:所述S4中S是弹丸最大横截面积;p是火药燃气平均压力;m是弹丸的质量;l是弹丸行程;v是弹丸速度;是次要功系数。
6.根据如权利要求1所述的发射药燃烧产生等离子体规律研究方法,其特征在于:所述S4中lψ是药室自由容积缩径长,ω是装药质量,θ是火药热力参数,ψ是火药已燃百分数,f是火药力,Δ是装填密度,ρp是火药密度,α是余容。
7.根据如权利要求1所述的发射药燃烧产生等离子体规律研究方法,其特征在于:所述S5中T为热力学温度,Ei为离子的电离电位,go为原子基态的统计权重,gi为离子基态的统计权重,me为电子质量,ne为电子密度,ni为离子密度,n0为原子密度,对于碱金属来说,的值约为1,其他气体一般为2。
8.根据如权利要求1所述的发射药燃烧产生等离子体规律研究方法,其特征在于:所述S5中的等离子体中电子密度
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910559092.XA CN110289053A (zh) | 2019-06-26 | 2019-06-26 | 发射药燃烧产生等离子体规律研究方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910559092.XA CN110289053A (zh) | 2019-06-26 | 2019-06-26 | 发射药燃烧产生等离子体规律研究方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110289053A true CN110289053A (zh) | 2019-09-27 |
Family
ID=68005941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910559092.XA Pending CN110289053A (zh) | 2019-06-26 | 2019-06-26 | 发射药燃烧产生等离子体规律研究方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110289053A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110765407A (zh) * | 2019-10-12 | 2020-02-07 | 中北大学 | 一种片状多层发射装药内弹道特性参数的计算方法 |
CN113643762A (zh) * | 2021-08-13 | 2021-11-12 | 中国人民解放军陆军装甲兵学院 | 计算电离种子影响的***、方法、装置及电子设备 |
CN115618171A (zh) * | 2022-06-06 | 2023-01-17 | 北京理工大学 | 一种基于同伦算法的推进剂燃烧平衡产物求解方法 |
CN116579259A (zh) * | 2023-04-24 | 2023-08-11 | 中国人民解放军陆军装甲兵学院 | 弹道三维瞬态流场建模及多物理场数值计算方法及装置 |
CN116794151A (zh) * | 2023-06-27 | 2023-09-22 | 中国人民解放军陆军装甲兵学院 | 常压条件下发射药燃烧产物电子密度的实验测试方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202797A (zh) * | 2016-07-21 | 2016-12-07 | 中国人民解放军装甲兵工程学院 | 一种火药燃烧时炮管内等离子体密度计算方法 |
CN107025325A (zh) * | 2017-01-03 | 2017-08-08 | 中国人民解放军装甲兵工程学院 | 一种使火药燃气生成等离子体的电离种子优化方法 |
CN107944150A (zh) * | 2017-11-27 | 2018-04-20 | 南京航空航天大学 | 内外流多组元化学反应流场一体化数值模拟模型 |
-
2019
- 2019-06-26 CN CN201910559092.XA patent/CN110289053A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202797A (zh) * | 2016-07-21 | 2016-12-07 | 中国人民解放军装甲兵工程学院 | 一种火药燃烧时炮管内等离子体密度计算方法 |
CN107025325A (zh) * | 2017-01-03 | 2017-08-08 | 中国人民解放军装甲兵工程学院 | 一种使火药燃气生成等离子体的电离种子优化方法 |
CN107944150A (zh) * | 2017-11-27 | 2018-04-20 | 南京航空航天大学 | 内外流多组元化学反应流场一体化数值模拟模型 |
Non-Patent Citations (6)
Title |
---|
TIKAI ZHENG等: "Numerical study of plasma assisted combustion for a rocket combustor using GCH4/GOX as propellants", 《JOURNAL OF PHYSICS: CONFERENCE SERIES》 * |
李俊: "火药燃气作用下电离种子产生等离子体研究", 《核聚变与等离子体物理》 * |
李晓刚 等: "火药燃烧等离子体电导率理论计算研究", 《火工品》 * |
毛保全 等: "含钾盐添加剂的火药燃烧产物导电特性研究", 《兵工学报》 * |
毛保全 等: "基于电离种子的发射药燃烧生成等离子体研究", 《火工品》 * |
钟孟春 等: "火药燃烧等离子体电导率影响因素研究", 《装甲兵工程学院学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110765407A (zh) * | 2019-10-12 | 2020-02-07 | 中北大学 | 一种片状多层发射装药内弹道特性参数的计算方法 |
CN110765407B (zh) * | 2019-10-12 | 2023-04-25 | 中北大学 | 一种片状多层发射装药内弹道特性参数的计算方法 |
CN113643762A (zh) * | 2021-08-13 | 2021-11-12 | 中国人民解放军陆军装甲兵学院 | 计算电离种子影响的***、方法、装置及电子设备 |
CN113643762B (zh) * | 2021-08-13 | 2024-05-31 | 中国人民解放军陆军装甲兵学院 | 计算电离种子影响的***、方法、装置及电子设备 |
CN115618171A (zh) * | 2022-06-06 | 2023-01-17 | 北京理工大学 | 一种基于同伦算法的推进剂燃烧平衡产物求解方法 |
CN115618171B (zh) * | 2022-06-06 | 2023-10-24 | 北京理工大学 | 一种基于同伦算法的推进剂燃烧平衡产物求解方法 |
CN116579259A (zh) * | 2023-04-24 | 2023-08-11 | 中国人民解放军陆军装甲兵学院 | 弹道三维瞬态流场建模及多物理场数值计算方法及装置 |
CN116579259B (zh) * | 2023-04-24 | 2024-02-09 | 中国人民解放军陆军装甲兵学院 | 弹道三维瞬态流场建模及多物理场数值计算方法及装置 |
CN116794151A (zh) * | 2023-06-27 | 2023-09-22 | 中国人民解放军陆军装甲兵学院 | 常压条件下发射药燃烧产物电子密度的实验测试方法 |
CN116794151B (zh) * | 2023-06-27 | 2024-02-02 | 中国人民解放军陆军装甲兵学院 | 常压条件下发射药燃烧产物电子密度的实验测试方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110289053A (zh) | 发射药燃烧产生等离子体规律研究方法 | |
CN110287599A (zh) | 一种发射药产生等离子体规律模型的构建方法 | |
CN106202797B (zh) | 一种火药燃烧时炮管内等离子体密度计算方法 | |
CN109446601B (zh) | 一种弹丸起始扰动的不确定优化方法 | |
CN107025325A (zh) | 一种使火药燃气生成等离子体的电离种子优化方法 | |
CN115795925B (zh) | 基于反应气体爆轰驱动的二级轻气炮炮口速度计算方法 | |
CN106908482A (zh) | 一种提高火药燃烧产物导电特性装置及其电导率计算方法 | |
Tran et al. | Analysis of gas flow losses in a gas-operated gun | |
Wang et al. | Influence of magnetically confined plasma on the muzzle velocity of gun projectile | |
CN111783268B (zh) | 燃烧轻气炮内弹道参数预估方法 | |
M. Rashad et al. | Two-phase flow interior ballistics model of naval large caliber guided projectile gun system | |
Rashad et al. | Interior ballistic two-phase flow model of guided-projectile gun system utilizing stick propellant charge | |
RU2500659C2 (ru) | Порох староверова - 2 | |
Lim | Predicting the accuracy of unguided artillery projectiles | |
RU180685U1 (ru) | Газопароимпульсная корректирующая двигательная установка для космических аппаратов | |
Kong et al. | Effect of modified method based on recoil motion on interior ballistic performance of ultralight high-low pressure artillery | |
Cao | Parametric cycle analysis of continuous rotating detonation ejector-augmented rocket engine | |
Ma et al. | Study on micro recoil mechanism of the weapon with a nozzle and two chambers separated by a partition | |
Adams et al. | Dynamic modeling and design of a bulk-loaded liquid monopropellant powered rifle | |
Deng et al. | Detonation Impact Characteristics of Nitrocellulose in Nail Gun | |
Ma et al. | Investigation of combustion process in a modular charge system with a primer gap | |
Chen et al. | Calculation and Analysis of Interior Ballistic Performance of Propellant Containing Ionized Seeds | |
Fan et al. | Study on theoretical calculation of quasi-static pressure for aluminized explosive in confined space | |
CN116150952A (zh) | 一种基于密闭爆发器试验的球扁药内弹道数值仿真方法 | |
Yang et al. | Design of a new multi-channel high-low pressure micro-ejection system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190927 |
|
RJ01 | Rejection of invention patent application after publication |