CN106342303B - 一种固体火箭冲压发动机补燃室参数的计算方法 - Google Patents

一种固体火箭冲压发动机补燃室参数的计算方法

Info

Publication number
CN106342303B
CN106342303B CN200710081821.2A CN200710081821A CN106342303B CN 106342303 B CN106342303 B CN 106342303B CN 200710081821 A CN200710081821 A CN 200710081821A CN 106342303 B CN106342303 B CN 106342303B
Authority
CN
China
Prior art keywords
air
ramjet
temperature
calculate
centerdot
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.)
Expired - Fee Related
Application number
CN200710081821.2A
Other languages
English (en)
Inventor
曹军伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Airborne Missile Academy
Original Assignee
China Airborne Missile Academy
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Airborne Missile Academy filed Critical China Airborne Missile Academy
Priority to CN200710081821.2A priority Critical patent/CN106342303B/zh
Application granted granted Critical
Publication of CN106342303B publication Critical patent/CN106342303B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Of Engines (AREA)

Abstract

本发明属于固体火箭冲压发动机技术,涉及对固体火箭冲压发动机补燃室参数的计算方法的改进。其特征在于,计算的步骤如下:获取试验参数;计算试验状态下冲压补燃室温度Tbj;计算规定空燃比Ng下的冲压补燃室内的温度Tbj1;计算在来流空气总温发生变化下的补燃室内的温度Tbj2;计算冲压补燃室尾部压强Pbj。本发明基于冲压发动机工作过程中实际因素、计算速度快、计算结果准确。

Description

一种固体火箭冲压发动机补燃室参数的计算方法
技术领域
本发明属于固体火箭冲压发动机技术,涉及对固体火箭冲压发动机补燃室参数的计算方法的改进。
背景技术
在以往的固体火箭冲压发动机设计中,采用热力学计算的办法或者三维燃烧流动数值模拟的方法进行冲压发动机的补燃室参数计算。这两种方法存在的缺点有:一方面计算量太大、速度慢,满足不了在固体火箭冲压发动机优化设计和导弹变弹道多状态下大量快速计算的要求;另一方面没有考虑冲压发动机工作过程中实际因素的影响,计算结果与实际值相差较大。
发明内容
本发明的目的是:提供一种计算量小、计算速度快、基于冲压发动机工作过程中实际因素的固体火箭冲压发动机补燃室参数计算方法,该方法的计算结果与实际值较为一致。
本发明的技术方案是:一种固体火箭冲压发动机补燃室参数的计算方法,其特征在于,包括如下步骤:
1、获取试验参数;进行固体火箭冲压发动机试车,测量并记录以下参数:
qmas——实测空气质量流率,kg/s;
qmfs——实测燃气质量流率,kg/s;
Pbs——冲压补燃室实测压强,MPa;
Ns——实测的空燃比;
Tas——实测的来流进气总温,单位K;
2、计算试验状态下冲压补燃室温度Tbj;由流过冲压喷管的流量公式(1)计算出Tbj
T b j = ( m · P b s · A i q m a s + q m f s ) 2 - - - ( 1 )
式(1)中:m——密流系数;
At——冲压喷管喉部面积,cm2
3、计算规定空燃比Ng下的冲压补燃室内的温度Tbj1;设定空气来流总温不变,在规定空燃比Ng下,计算冲压补燃室内的温度Tbj1,用公式(2)进行计算:
T b j 1 = C P a ( N g - N s ) T a s + ( 1 + N s ) C P b · T b j C P a ( N g - N s ) + ( 1 + N s ) C P b - - - ( 2 )
式中:Cpa——空气定压比容;
Ng——规定的空燃比;其取值范围是:2~30;
Cpb——冲压补燃室燃气的定压比容;其取值范围是:1200~1800J/(kg·K);
4、计算在来流空气总温发生变化下的补燃室内的温度Tbj2;规定的来流空气总温为Tag,在规定的空燃比Ng下,冲压补燃室内的温度Tbj2由公式(3)计算:
T b j 2 = T b j 1 + C P a · N g ( T a g - T a s ) C P b ( 1 + N g ) - - - ( 3 )
式中:Tag——规定的来流空气总温,单位K;
5、计算冲压补燃室尾部压强Pbj;设定一个燃气质量流率qmfg,单位kg/s,在空燃比Ng下和空气来流总温Tag下,由公式(4)计算冲压补燃室尾部压强Pbj
P b j = ( q m f s + q m f s · N g ) q m f g T b j q m f s · m · A t - - - ( 4 )
步骤4中的Tbj2和5中的Pbj为所计算的固体火箭冲压发动机补燃室参数。
本发明的优点是:本发明基于冲压发动机工作过程中实际因素、计算速度快、计算结果准确,能够应用于固体火箭冲压发动机优化设计和飞行过程中变弹道多状态情况下固体火箭冲压发动机性能的快速准确计算,同时节省了导弹上有限的计算资源。
具体实施方式
下面对本发明做进一步详细说明。本发明的的计算方法,其特征在于,包括如下步骤:
1、获取试验参数;进行固体火箭冲压发动机试车,测量并记录以下参数:
qmas——实测空气质量流率,kg/s;
qmfs——实测燃气质量流率,kg/s;
Pbs——冲压补燃室实测压强,MPa;
Ns——实测的空燃比;
Tas——实测的来流进气总温,单位K;
2、计算试验状态下冲压补燃室温度Tbj;由流过冲压喷管的流量公式(1)计算出Tbj
T b j = ( m · P b s · A t q m a s + q m f s ) 2 - - - ( 1 )
式(1)中:m——密流系数;其计算方法参见《气体动力学基础》,潘锦珊,国防工业出版社,1980,6。
At——冲压喷管喉部面积,cm2
3、计算规定空燃比Ng下的冲压补燃室内的温度Tbj1;设定空气来流总温不变,在规定空燃比Ng下,计算冲压补燃室内的温度Tbj1,用公式(2)进行计算:
T b j 1 = C P a ( N g - N s ) T a s + ( 1 + N s ) C P b · T b j C P a ( N g - N s ) + ( 1 + N s ) C P b - - - ( 2 )
式中:Cpa——空气定压比容;
Ng——规定的空燃比;其取值范围是:2~30;
Cpb——冲压补燃室燃气的定压比容;其取值范围是:1200~1800J/(kg·K);
4、计算在来流空气总温发生变化下的补燃室内的温度Tbj2;规定的来流空气总温为Tag,在规定的空燃比Ng下,冲压补燃室内的温度Tbj2由公式(3)计算:
T b j 2 = T b j 1 + C P a · N g ( T a g - T a s ) C P b ( 1 + N g ) - - - ( 3 )
式中:Tag——规定的来流空气总温,单位K;
5、计算冲压补燃室尾部压强Pbj;设定一个燃气质量流率qmfg,单位kg/s,在空燃比Ng下和空气来流总温Tag下,由公式(4)计算冲压补燃室尾部压强Pbj
P b j = ( q m f s + q m f s · N g ) q m f g T b j q m f s · m · A t - - - ( 4 )
燃气质量流率qmfg的计算方法参见《固体火箭发动机设计》中有关固体火箭发动机燃气质量流率计算的内容,王元有等,国防工业出版社,1984,11。
步骤4中的Tbj2和5中的Pbj为所计算的固体火箭冲压发动机补燃室参数。
实施例。
通过固体火箭冲压发动机试车,测得的参数如下:实测的空气质量流率qmas为5.563kg/s,实测的燃气质量流率qmfs为0.589,kg/s,冲压补燃室实测压强Pbs为0.291MPa,实测的空燃比Ns为9.44,实测的来流进气总温Tas为573K。通过该台发动机的试验数据,要计算规定来流进气总温Tag为610K、规定空燃比Ng为12和燃气质量流率为0.62kg/s下的固体火箭冲压发动机补燃室参数Tbj2和Pbj
按照上述介绍的方法,并取m=0.04,Cpb为1600J/(kg·K),通过上述步骤的运算,得到固体火箭冲压发动机补燃室参数Tbj2为3109K和Pbj为0.361MPa。

Claims (1)

1.一种固体火箭冲压发动机补燃室参数的计算方法,其特征在于,包括如下步骤:
1.1、获取试验参数;进行固体火箭冲压发动机试车,测量并记录以下参数:
qmas——实测空气质量流率,kg/s;
qmfs——实测燃气质量流率,kg/s;
Pbs——冲压补燃室实测压强,MPa;
Ns——实测的空燃比;
Tas——实测的来流进气总温,K;
1.2、计算试验状态下冲压补燃室温度Tbj;由流过冲压喷管的流量公式(1)计算出Tbj
T b j = ( m · P b s · A t q m a s + q m f s ) 2 - - - ( 1 )
式(1)中:m——密流系数;
At——冲压喷管喉部面积,cm2
1.3、计算规定空燃比Ng下的冲压补燃室内的温度Tbj1;设定空气来流总温不变,在规定空燃比Ng下,计算冲压补燃室内的温度Tbj1,用公式(2)进行计算:
T b j 1 = C P a ( N g - N s ) T a s + ( 1 + N s ) C P b · T b j C P a ( N g - N s ) + ( 1 + N s ) C P b - - - ( 2 )
式中:Cpa——空气定压比热容;
Ng——规定的空燃比;其取值范围是:2~30;
Cpb——冲压补燃室燃气的定压比热容;其取值范围是:1200~1800J/(kg·K);
1.4、计算在来流空气总温发生变化下的补燃室内的温度Tbj2;规定的来流空气总温为Tag,在规定的空燃比Ng下,冲压补燃室内的温度Tbj2由公式(3)计算:
T b j 2 = T b j 1 + C P a · N g ( T a g - T a s ) C P b ( 1 + N g ) - - - ( 3 )
式中:Tag——规定的来流空气总温,单位K;
1.5、计算冲压补燃室尾部压强Pbj;给定一个燃气质量流率qmfg,单位kg/s,在空燃比Ng下和空气来流总温Tag下,由公式(4)计算冲压补燃室尾部压强Pbj
P b j = ( q m f s + q m f s · N g ) q m f g T b j q m f s · m · A t - - - ( 4 )
步骤1.4中的Tbj2和1.5中的Pbj为所计算的固体火箭冲压发动机补燃室参数。
CN200710081821.2A 2007-07-30 2007-07-30 一种固体火箭冲压发动机补燃室参数的计算方法 Expired - Fee Related CN106342303B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200710081821.2A CN106342303B (zh) 2007-07-30 2007-07-30 一种固体火箭冲压发动机补燃室参数的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200710081821.2A CN106342303B (zh) 2007-07-30 2007-07-30 一种固体火箭冲压发动机补燃室参数的计算方法

Publications (1)

Publication Number Publication Date
CN106342303B true CN106342303B (zh) 2011-06-22

Family

ID=57797825

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200710081821.2A Expired - Fee Related CN106342303B (zh) 2007-07-30 2007-07-30 一种固体火箭冲压发动机补燃室参数的计算方法

Country Status (1)

Country Link
CN (1) CN106342303B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109736970A (zh) * 2019-01-11 2019-05-10 中国人民解放军国防科技大学 一种含硼固体火箭冲压发动机补燃室特征长度设计方法
CN111058968A (zh) * 2019-12-12 2020-04-24 西安近代化学研究所 一种双燃烧室固体火箭发动机小燃烧室压强的计算方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109736970A (zh) * 2019-01-11 2019-05-10 中国人民解放军国防科技大学 一种含硼固体火箭冲压发动机补燃室特征长度设计方法
CN109736970B (zh) * 2019-01-11 2020-04-28 中国人民解放军国防科技大学 一种含硼固体火箭冲压发动机补燃室特征长度设计方法
CN111058968A (zh) * 2019-12-12 2020-04-24 西安近代化学研究所 一种双燃烧室固体火箭发动机小燃烧室压强的计算方法

Similar Documents

Publication Publication Date Title
CN103366078B (zh) 一种航空发动机轴承腔通风设计方法
CN101782029B (zh) 气-气喷嘴流量特性试验装置
CN106092420B (zh) 间接测量发动机有效推力的方法
CN104345118B (zh) 固体推进剂多靶线动态燃烧性能测试系统及方法
CN110309552B (zh) 一种考虑质量引射效应的飞行器湍流预测方法及系统
CN102656355A (zh) 内燃机的控制装置
CN111157248B (zh) 基于地面直连试验的冲压发动机及其燃烧室性能评估方法
CN204964196U (zh) 液体发动机喷嘴液流试验测试装置
CN103870683B (zh) 一种超燃冲压发动机燃烧室性能预估方法
CN110953090A (zh) 模拟等压多通道火箭发动机内两相流烧蚀环境的试验装置
CN103473396B (zh) 一种运载火箭低温贮箱换热过程仿真方法
CN106342303B (zh) 一种固体火箭冲压发动机补燃室参数的计算方法
CN105446167A (zh) 高超声速超燃冲压发动机实时模型、仿真方法
Forth et al. SCALING PARAMETERS IN FILM− COOLING
CN108644031A (zh) 一种固体火箭发动机绝热层烧蚀率测试方法
CN104318011B (zh) 一种基于实验与仿真互相耦合的真空羽流效应评估方法
Gu et al. A novel experimental method to the internal thrust of rocket-based combined-cycle engine
CN107655691A (zh) 一种飞行器进气道喘振锤击波试验装置及方法
CN105512405A (zh) 一种mpc引射喷管直径的优化设计方法
CN110442934A (zh) 一种考虑固体发动机尾喷流辐射的高精度气动热计算方法
Christensen et al. Measurement of Heat Transfer Inside a Channel Using External Infrared Thermography
CN103968907B (zh) 一种超临界态和气态碳氢燃料密流测量装置及方法
CN106556479A (zh) 一种热耗测量装置及其热耗测量方法
CN115906699B (zh) 超快速预测水或蒸汽管道破口处临界质量流速的方法
RU2586792C1 (ru) Способ определения коэффициента расхода газа через сопловой аппарат турбины двухконтурного газотурбинного двигателя

Legal Events

Date Code Title Description
GR03 Grant of secret patent right
DC01 Secret patent status has been lifted
DCSP Declassification of secret patent
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110622

Termination date: 20180730