CN110516410B - 一种飞行器总体设计参数的确定方法及装置 - Google Patents
一种飞行器总体设计参数的确定方法及装置 Download PDFInfo
- Publication number
- CN110516410B CN110516410B CN201911001494.4A CN201911001494A CN110516410B CN 110516410 B CN110516410 B CN 110516410B CN 201911001494 A CN201911001494 A CN 201911001494A CN 110516410 B CN110516410 B CN 110516410B
- Authority
- CN
- China
- Prior art keywords
- model
- aircraft
- overall design
- parameter
- determining
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000000446 fuel Substances 0.000 claims abstract description 64
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 35
- 238000005457 optimization Methods 0.000 claims abstract description 35
- 230000002068 genetic effect Effects 0.000 claims abstract description 31
- 238000011084 recovery Methods 0.000 claims description 26
- 239000007789 gas Substances 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 18
- 230000003068 static effect Effects 0.000 claims description 17
- 238000010276 construction Methods 0.000 claims description 11
- 230000035945 sensitivity Effects 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 239000002737 fuel gas Substances 0.000 claims description 5
- 239000002245 particle Substances 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims 1
- 230000000875 corresponding effect Effects 0.000 description 24
- 239000007921 spray Substances 0.000 description 11
- 239000003380 propellant Substances 0.000 description 9
- 230000001808 coupling effect Effects 0.000 description 8
- 238000002485 combustion reaction Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000009467 reduction Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 4
- 230000002950 deficient Effects 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 239000001301 oxygen Substances 0.000 description 4
- 229910052760 oxygen Inorganic materials 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000000376 reactant Substances 0.000 description 3
- 229920002121 Hydroxyl-terminated polybutadiene Polymers 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 239000011777 magnesium Substances 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- SNAAJJQQZSMGQD-UHFFFAOYSA-N aluminum magnesium Chemical compound [Mg].[Al] SNAAJJQQZSMGQD-UHFFFAOYSA-N 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 229910052749 magnesium Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000007800 oxidant agent Substances 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明公开了一种飞行器总体设计参数的确定方法,包括:确定飞行器的各个组成部分的设计参数;根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解飞行器总体设计模型,得到飞行器总体设计参数;参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数。该方法能够考虑动力性能和弹道的耦合效应,提高飞行器的总体性能。本发明公开的一种飞行器总体设计参数的确定装置、设备及可读存储介质,也同样具有上述技术效果。
Description
技术领域
本发明涉及飞行器设计技术领域,更具体地说,涉及一种飞行器总体设计参数的确定方法、装置、设备及可读存储介质。
背景技术
飞行器设计中涉及到许多学科,从整体系统的角度来说,各个子系统间必然存在一定程度的相互作用和耦合效应。由于子系统均有各自独立的分析和设计工具,传统次序型设计往往无法综合考虑,从而没能充分考虑到各系统间的耦合因素。应用多学科设计优化解决上述整体系统的设计问题,既提高了设计效率,还改善了产品的性能。
在现有技术中,在进行多学科设计优化时,未考虑到动力参数和弹道参数之间的耦合效应,从而导致当飞行器的飞行空域和速域变化明显时,动力参数和弹道参数的耦合特征明显。当采用冲压发动机为飞行器提供动力时,若仍按照现有技术的设计方式设计总体参数,由于未考虑动力参数和弹道参数的耦合效应,会使得飞行器的总体性能降低。
因此,如何降低动力参数和弹道参数的耦合效应,提高飞行器的总体性能,是本领域技术人员需要解决的问题。
发明内容
本发明的目的在于提供一种飞行器总体设计参数的确定方法、装置、设备及可读存储介质,以降低动力参数和弹道参数的耦合效应,提高飞行器的总体性能。
为实现上述目的,本发明实施例提供了如下技术方案:
一种飞行器总体设计参数的确定方法,包括:
确定飞行器的各个组成部分和每个组成部分包括的设计参数,所述飞行器由冲压发动机提供动力;
根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;
根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数;
其中,所述参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,所述冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数;
其中,所述按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数,包括:
按照灵敏度分析所述飞行器总体设计模型,得到第一模型和第二模型;
当所述优化目标为射程最大时,利用遗传算法求解所述第二模型,得到第一最优解;
将所述第一最优解作为所述第一模型的初始迭代值,并利用遗传算法求解所述第一模型,得到第二最优解;
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
其中,所述飞行器质量模型的构建步骤包括:
确定所述飞行器的弹体组成部分和每个弹体组成部分的质量,并根据每个弹体组成部分的质量构建所述飞行器质量模型;
其中,所述飞行器气动模型的构建步骤包括:
采用部件组合法构建所述飞行器气动模型;所述飞行器气动模型为:
其中,所述弹道参数模型的构建步骤包括:
确定所述飞行器的弹道设计参数,并采用三自由度质点构建所述弹道参数模型,其中,所述弹道参数模型中的纵向运动方程组为:
其中,V表示飞行速度,m表示飞行器质量,F表示发动机推力,X表示气动阻力,Y表示气动升力,表示俯仰力矩,表示转动惯量,表示攻角,表示速度倾角,表示俯仰角,表示俯仰角速度,表示单位时间内燃料消耗量,表示飞行器射程,表示飞行器高度,t表示时间,g表示重力加速度。
其中,所述冲压发动机性能模型的构建步骤包括:
确定所述冲压发动机的组成部分,所述冲压发动机的组成部分包括:燃气发生器、补燃室和进气道,所述进气道为两个二元进气道;
根据所述冲压发动机的组成部分计算发动机推力、发动机比冲、所述空燃比和所述补燃室压力;
构建包括所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标的所述冲压发动机性能模型;所述冲压发动机性能模型为:
其中,F表示发动机推力,表示燃料流量,所述燃料流量为所述燃气发生器参数,v a 表示来流空气速度,v e 表示喷管出口速度,p e 表示喷管出口静压,p a 表示来流空气静压,A e 表示喷管出口面积;I s 表示发动机比冲;表示空燃比,表示空气流量;表示进气道总压恢复系数,所述进气道总压恢复系数为所述进气道参数,表示进气道出口总压,表示来流空气总压;表示补燃室压力,表示补燃室总压恢复系数。
其中,所述按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数,包括:
按照灵敏度分析所述飞行器总体设计模型,得到第一模型和第二模型;
当所述优化目标为射程最大时,利用遗传算法求解所述第二模型,得到第一最优解;
将所述第一最优解作为所述第一模型的初始迭代值,并利用遗传算法求解所述第一模型,得到第二最优解;
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
其中,所述根据所有参数模型确定飞行器总体设计模型,包括:
根据所述飞行器质量模型、所述飞行器气动模型、所述弹道参数模型和所述冲压发动机性能模型确定所述飞行器总体设计模型;
其中,所述飞行器总体设计模型的设计变量为:发动机动力参数和飞行器几何参数,弹道形式为高空-俯冲;
一种飞行器总体设计参数的确定装置,包括:
确定模块,用于确定飞行器的各个组成部分和每个组成部分包括的设计参数,所述飞行器由冲压发动机提供动力;
构建模块,用于根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;
执行模块,用于根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数;
其中,所述参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,所述冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数;
其中,所述执行模块具体用于:
按照灵敏度分析所述飞行器总体设计模型,得到第一模型和第二模型;
当所述优化目标为射程最大时,利用遗传算法求解所述第二模型,得到第一最优解;
将所述第一最优解作为所述第一模型的初始迭代值,并利用遗传算法求解所述第一模型,得到第二最优解;
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
一种飞行器总体设计参数的确定设备,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现上述任意一项所述的飞行器总体设计参数的确定方法的步骤。
一种可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现上述任意一项所述的飞行器总体设计参数的确定方法的步骤。
通过以上方案可知,本发明实施例提供的一种飞行器总体设计参数的确定方法,包括:确定飞行器的各个组成部分和每个组成部分包括的设计参数,所述飞行器由冲压发动机提供动力;根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数;其中,所述参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,所述冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数。其中,按照预设的优化目标和遗传算法求解飞行器总体设计模型,得到优化目标对应的飞行器总体设计参数,包括:按照灵敏度分析飞行器总体设计模型,得到第一模型和第二模型;当优化目标为射程最大时,利用遗传算法求解第二模型,得到第一最优解;将第一最优解作为第一模型的初始迭代值,并利用遗传算法求解第一模型,得到第二最优解;将第二最优解作为射程最大对应的飞行器总体设计参数。
可见,所述方法在确定飞行器总体设计模型时,使飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型均参与其中,且冲压发动机性能模型中的设计参数包括进气道参数,也就是说,该方法在确定飞行器总体设计模型时,考虑了所有模型之间的耦合效应。
具体的,本申请按照灵敏度分析飞行器总体设计模型,得到第一模型和第二模型;当所述优化目标为射程最大时,利用遗传算法求解第二模型,得到第一最优解;将第一最优解作为第一模型的初始迭代值,并利用遗传算法求解第一模型,得到第二最优解;将第二最优解作为射程最大对应的飞行器总体设计参数。如此则利用遗传算法求解得到了射程最大对应的总体设计参数。由于本申请提供的飞行器总体设计模型中涵盖了飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,因此遗传算法在求解过程中能够考虑动力参数和弹道参数的耦合效应,降低动力参数和弹道参数的耦合,从而提高飞行器的总体性能。
相应地,本发明实施例提供的一种飞行器总体设计参数的确定装置、设备及可读存储介质,也同样具有上述技术效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例公开的一种飞行器总体设计参数的确定方法流程图;
图2为本发明实施例公开的一种飞行器总体设计参数的确定装置示意图;
图3为本发明实施例公开的一种飞行器总体设计参数的确定设备示意图;
图4为本发明实施例公开的导弹总体设计参数的设计总流程示意图;
图5为本发明实施例公开的导弹弹体的各个组成部分示意图;
图6为本发明实施例公开的升力系数随马赫数变化示意图;
图7为本发明实施例公开的阻力系数随马赫数变化示意图;
图8为本发明实施例公开理论燃烧温度与理论比冲随空燃比变化曲线示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种飞行器总体设计参数的确定方法、装置、设备及可读存储介质,以降低动力参数和弹道参数的耦合效应,提高飞行器的总体性能。
参见图1,本发明实施例提供的一种飞行器总体设计参数的确定方法,包括:
S101、确定飞行器的各个组成部分和每个组成部分包括的设计参数,飞行器由冲压发动机提供动力;
S102、根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;
S103、根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解飞行器总体设计模型,得到优化目标对应的飞行器总体设计参数;
其中,参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数。
其中,所述飞行器质量模型的构建步骤包括:确定所述飞行器的弹体组成部分和每个弹体组成部分的质量,并根据每个弹体组成部分的质量构建所述飞行器质量模型。
其中,所述飞行器气动模型的构建步骤包括:
采用部件组合法构建所述飞行器气动模型;所述飞行器气动模型为:
其中,所述弹道参数模型的构建步骤包括:确定所述飞行器的弹道设计参数,并采用三自由度质点构建所述弹道参数模型,其中,所述弹道参数模型中的纵向运动方程组为:
其中,V表示飞行速度,m表示飞行器质量,F表示发动机推力,X表示气动阻力,Y表示气动升力,表示俯仰力矩,表示转动惯量,表示攻角,表示速度倾角,表示俯仰角,表示俯仰角速度,表示单位时间内燃料消耗量,表示飞行器射程,表示飞行器高度,t表示时间,g表示重力加速度。
其中,所述冲压发动机性能模型的构建步骤包括:确定所述冲压发动机的组成部分,所述冲压发动机的组成部分包括:燃气发生器、补燃室和进气道,所述进气道为两个二元进气道;根据所述冲压发动机的组成部分计算发动机推力、发动机比冲、所述空燃比和所述补燃室压力;根据所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标构建所述冲压发动机性能模型,也就是构建包括所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标的所述冲压发动机性能模型。
其中,所述冲压发动机性能模型为:
其中,F表示发动机推力,表示燃料流量,所述燃料流量为所述燃气发生器参数,v a 表示来流空气速度,v e 表示喷管出口速度,p e 表示喷管出口静压,p a 表示来流空气静压,A e 表示喷管出口面积;I s 表示发动机比冲;表示空燃比,表示空气流量;表示进气道总压恢复系数,所述进气道总压恢复系数为所述进气道参数,表示进气道出口总压,表示来流空气总压;表示补燃室压力,表示补燃室总压恢复系数。
需要说明的是,构建冲压发动机性能模型的目的为:基于冲压发动机的技术指标、空燃比、补燃室压力、进气道参数、燃气发生器参数等各方参数,求解发动机推力和发动机比冲。发动机推力和发动机比冲可表示冲压发动机的性能;发动机推力越大,发动机比冲越高,那么发动机性能就越高。发动机推力和发动机比冲的具体计算请参见下文描述的相应计算公式。发动机比冲即为理论比冲。
可见,本实施例提供了一种飞行器总体设计参数的确定方法,所述方法在确定飞行器总体设计模型时,使飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型均参与其中,且冲压发动机性能模型中的设计参数包括进气道参数,也就是说,该方法在确定飞行器总体设计模型时,考虑了所有模型之间的耦合效应,利用遗传算法求解得到了射程最大对应的总体设计参数。由于本申请提供的飞行器总体设计模型中涵盖了飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,因此遗传算法在求解过程中能够考虑动力参数和弹道参数的耦合效应,降低动力参数和弹道参数的耦合,从而提高飞行器的总体性能。
基于上述实施例,需要说明的是,所述按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数,包括:
按照灵敏度分析所述飞行器总体设计模型,得到第一模型和第二模型;
当所述优化目标为射程最大时,利用遗传算法求解所述第二模型,得到第一最优解;
将所述第一最优解作为所述第一模型的初始迭代值,并利用遗传算法求解所述第一模型,得到第二最优解;
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
需要说明的是,第一模型和第二模型的计算精度是相对而言的,第一模型的计算速度比第二模型的计算速度慢,但其精度比第二模型的精度高。也就是说,按照灵敏度分析飞行器总体设计模型,可得到具有不同计算精度的第一模型和第二模型。在具体计算时,首先用遗传算法求解较低精度的第二模型,得到第一最优解;进而将第一最优解作为较高精度的第一模型的初始迭代值,并利用遗传算法求解第一模型,从而得到第二最优解,那么第二最优解即为射程最大对应的飞行器总体设计参数。
基于上述实施例,需要说明的是,所述根据所有参数模型确定飞行器总体设计模型,包括:
根据所述飞行器质量模型、所述飞行器气动模型、所述弹道参数模型和所述冲压发动机性能模型确定所述飞行器总体设计模型;
其中,所述飞行器总体设计模型的设计变量为:发动机动力参数和飞行器几何参数,弹道形式为高空-俯冲;
下面对本发明实施例提供的一种飞行器总体设计参数的确定装置进行介绍,下文描述的一种飞行器总体设计参数的确定装置与上文描述的一种飞行器总体设计参数的确定方法可以相互参照。
参见图2,本发明实施例提供的一种飞行器总体设计参数的确定装置,包括:
确定模块201,用于确定飞行器的各个组成部分和每个组成部分包括的设计参数,所述飞行器由冲压发动机提供动力;
构建模块202,用于根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;
执行模块203,用于根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数;
其中,所述参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,所述冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数。
其中,所述构建模块具体用于:
确定所述飞行器的弹体组成部分和每个弹体组成部分的质量,并根据每个弹体组成部分的质量构建所述飞行器质量模型;
其中,所述构建模块具体用于:
采用部件组合法构建所述飞行器气动模型;所述飞行器气动模型为:
其中,所述构建模块具体用于:
确定所述飞行器的弹道设计参数,并采用三自由度质点构建所述弹道参数模型,其中,所述弹道参数模型中的纵向运动方程组为:
其中,V表示飞行速度,m表示飞行器质量,F表示发动机推力,X表示气动阻力,Y表示气动升力,表示俯仰力矩,表示转动惯量,表示攻角,表示速度倾角,表示俯仰角,表示俯仰角速度,表示单位时间内燃料消耗量,表示飞行器射程,表示飞行器高度,t表示时间,g表示重力加速度。
其中,所述构建模块具体用于:
确定所述冲压发动机的组成部分,所述冲压发动机的组成部分包括:燃气发生器、补燃室和进气道,所述进气道为两个二元进气道;
根据所述冲压发动机的组成部分计算发动机推力、发动机比冲、所述空燃比和所述补燃室压力;
构建包括所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标的所述冲压发动机性能模型;所述冲压发动机性能模型为:
其中,F表示发动机推力,表示燃料流量,所述燃料流量为所述燃气发生器参数,v a 表示来流空气速度,v e 表示喷管出口速度,p e 表示喷管出口静压,p a 表示来流空气静压,A e 表示喷管出口面积;I s 表示发动机比冲;表示空燃比,表示空气流量;表示进气道总压恢复系数,所述进气道总压恢复系数为所述进气道参数,表示进气道出口总压,表示来流空气总压;表示补燃室压力,表示补燃室总压恢复系数。
其中,所述执行模块具体用于:
按照灵敏度分析所述飞行器总体设计模型,得到第一模型和第二模型;
当所述优化目标为射程最大时,利用遗传算法求解所述第二模型,得到第一最优解;将所述第一最优解作为所述第一模型的初始迭代值,并利用遗传算法求解所述第一模型,得到第二最优解;
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
其中,所述执行模块具体用于:
根据所述飞行器质量模型、所述飞行器气动模型、所述弹道参数模型和所述冲压发动机性能模型确定所述飞行器总体设计模型;
其中,所述飞行器总体设计模型的设计变量为:发动机动力参数和飞行器几何参数,弹道形式为高空-俯冲;
可见,本实施例提供了一种飞行器总体设计参数的确定装置,包括:确定模块、构建模块以及执行模块。首先由确定模块确定飞行器的各个组成部分和每个组成部分包括的设计参数,所述飞行器由冲压发动机提供动力;然后构建模块根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;最后执行模块根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数;其中,所述参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,所述冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数。如此各个模块之间分工合作,各司其职,从而降低了动力参数和弹道参数的耦合效应,提高了飞行器的总体性能。
下面对本发明实施例提供的一种飞行器总体设计参数的确定设备进行介绍,下文描述的一种飞行器总体设计参数的确定设备与上文描述的一种飞行器总体设计参数的确定方法及装置可以相互参照。
参见图3,本发明实施例提供的一种飞行器总体设计参数的确定设备,包括:
存储器301,用于存储计算机程序;
处理器302,用于执行所述计算机程序时实现上述任意实施例所述的飞行器总体设计参数的确定方法的步骤。
下面对本发明实施例提供的一种可读存储介质进行介绍,下文描述的一种可读存储介质与上文描述的一种飞行器总体设计参数的确定方法、装置及设备可以相互参照。
一种可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如上述任意实施例所述的飞行器总体设计参数的确定方法的步骤。
按照本发明提供的飞行器总体设计参数的确定方法,同样可以设计导弹的总体设计参数。请参见图4,图4为导弹总体设计参数的设计总流程示意图。具体实现过程如下:
1、确定总体设计指标和约束条件。导弹的约束条件可以设置为总射程最大或质量最小。总体设计指标即冲压发动机的技术指标,具体如表1所示。
表1
序号 | 项目 | 单位 |
1 | 巡航高度 | m |
2 | 巡航攻角 | 度 |
3 | 巡航马赫数 | / |
4 | 巡航发动机推力 | N |
5 | 工作高度 | km |
6 | 来流马赫数 | / |
7 | 发动机最大外径 | mm |
8 | 工作时间 | s |
9 | 发动机质量 | kg |
10 | 飞行攻角及侧滑角 | 度 |
11 | 最大横向过载 | g |
12 | 最大纵向过载 | g |
2、确定导弹各个组成部分,以及各个组成部分的设计参数。导弹各个组成部分请参见表2。
表2
其中,冲压发动机的进气道参数如表3所示。
表3
序号 | 项目 | 单位 |
1 | 设计飞行高度 | m |
2 | 设计飞行攻角 | 度 |
3 | 设计飞行马赫数 | / |
4 | 进气道结构形式 | / |
5 | 设计空气流量 | kg/s |
6 | 最大飞行攻角及侧滑角 | 度 |
7 | 飞行马赫数工作范围 | / |
冲压发动机的性能指标如表4所示。
表4
项目 | 单位 |
设计点高度 | m |
设计点马赫数 | / |
设计点空燃比 | / |
设计点补燃室压力 | MPa |
设计点燃气流量 | kg/s |
设计点空气流量 | kg/s |
比冲效率 | / |
设计点理论比冲 | N·s/kg |
设计点比冲 | N·s/kg |
冲压喷管喉径 | mm |
冲压喷管出口直径 | mm |
理论额定推力 | N |
额定推力 | N |
导弹的弹体参数请参见表5。
表5
3、构建导弹质量模型、导弹气动模型、弹道参数模型和冲压发动机性能模型。
1)构建导弹质量模型:导弹的结构质量是将导弹几何模型进行一定厚度的抽壳后,对此壳体赋予密度等性质,进而计算出来其质量、质心参数。弹上设备质量设定为固定值,每一个设备的质量赋值为给定值,假设每一个设备都为一均质长方体,其中心位置即为设备质心位置。
2)构建导弹气动模型:采用部件组合法计算阻力系数和升力系数,并根据计算得到的阻力系数和升力系数构建导弹气动模型;阻力系数的计算公式为:,升力系数的计算公式为:;其中,,表示空气密度,S表示参考面积,q表示动压头,X表示气动阻力,Y表示气动升力,V表示飞行速度。
其中,升力系数随马赫数变化示意图请参见图6,阻力系数随马赫数变化示意图请参见图7。
3)构建弹道参数模型:确定导弹的弹道设计参数,并采用三自由度质点构建弹道参数模型,其中,弹道参数模型中的纵向运动方程组为:
其中,V表示飞行速度,m表示飞行器质量,F表示发动机推力,X表示气动阻力,Y表示气动升力,表示俯仰力矩,表示转动惯量,表示攻角,表示速度倾角,表示俯仰角,表示俯仰角速度,表示单位时间内燃料消耗量,表示飞行器射程,表示飞行器高度,t表示时间,g表示重力加速度。
4)构建冲压发动机性能模型:采用二个二元进气道;采用铝镁贫氧推进剂;采用非壅塞式燃气发生器;其中,冲压发动机主要由进气道、燃气发生器和补燃室(包括冲压喷管)组成。
其中,所述冲压发动机性能模型为:
其中,F表示发动机推力,表示燃料流量,所述燃料流量为所述燃气发生器参数,v a 表示来流空气速度,v e 表示喷管出口速度,p e 表示喷管出口静压,p a 表示来流空气静压,A e 表示喷管出口面积;I s 表示发动机比冲;表示空燃比,表示空气流量;表示进气道总压恢复系数,所述进气道总压恢复系数为所述进气道参数,表示进气道出口总压,表示来流空气总压;表示补燃室压力,表示补燃室总压恢复系数。
分析上式可知,理论比冲单位质量推进剂产生的推力越大,发动机性能越优,而理论比冲与飞行弹道参数(空气速度、空气压强)密切相关。耦合性突出体现在相同的发动机在不同弹道条件下性能差异极大。换句话说,欲使飞行器按期望的高度、速度飞行的距离最远,一定存在满足上述弹道条件的发动机设计参数。
具体的,冲压发动机性能模型的具体构成过程包括:
a、参数定义:
(1)定义空燃比,空燃比即为空气流量与燃料流量之比,用公式表示为:
(2)理论空气量(又称为当量比)定义为1kg燃料完全燃烧所需要的空气量,用公式表示为:
其中,注脚f表示燃料,注脚a表示空气;M a 为空气相对分子质量;M f 为燃料相对分子质量;L为元素数;b i为第i种元素的摩尔数;v i为第i种元素的原子价。
(3)余气系数定义为空燃比与理论空气量之比:
计算平衡状态下燃烧产物的组分必须首先计算反应物假定化学式。固体火箭冲压发动机由空气提供氧化剂,1kg由贫氧推进剂和空气共同组成的反应物的假定化学式由下式计算得到:
b、性能计算:
比冲既反映了推进剂能量的大小又反映了推进剂在燃烧过程中能量转换的效率。固冲发动机初步设计中,首先要在设计点处找到一组可用、较好的比冲数据,然后根据设计的要求,设计燃料流量和空气流量。燃料流量是设计燃气发生器的重要参数。空气流量是设计进气道的重要参数。
具体的,喷管出口速度v e 表示为:
其中,H c 为补燃室燃烧产物的焓;H e 为喷管出口燃烧产物的焓。
额定推力F由下式计算得到:
而比冲表示为:
故而,由上述公式可得:
其中,下标t表示喷管喉部,下标e表示喷管出口处。故为喷管喉部燃气流量,为喷管喉部的平衡因素,p t 表示喷管喉部静压,为通用气体常数,为喷管喉部的气相产物摩尔数,T t 为喷管喉部温度;为喷管出口处燃气流量,为喷管出口的平衡因素,p e 表示喷管出口静压,T e 为喷管出口温度,p e 表示喷管出口静压。
因此比冲可表示为:
c、冲压喷管设计:
d、进气道设计:
进气道设计计算主要是计算进气道总压恢复系数和流量系数。
(1)总压恢复系数:
进气道总压恢复系数为:
不变的条件下,当下降1%时,同样下降1%。而的下降会导致喷管出口总压减小,这就使冲压喷管完全膨胀所能达到的减小导致推力下降。进气道总压恢复系数下降1%而推力下降则在1%以上。因此,在设计进气道时应尽量提高总压恢复系数。
(2)流量系数:
进气道流量系数随飞行状态的不同改变。在发动机初步设计中期望进气道流量系数尽可能高,这样发动机在流量调节上会有更大的空间。流量系数直接关系到冲压喷管喉径的设计,具体关系可用下式表示:
给定补燃室内燃气理论特征速度,在燃气流量变化时,若喉径过大,则总压会有一个较大的减小以保证流量的不变;则正激波吸入进气道内,进气道在超临界状态下工作;同时,总压的下降导致推力的下降。若喉径过小,下游节流程度增加,总压升高,正激波推出进口,产生严重溢流。另外,流量的改变还会使燃料能量的释放发生改变,使得理论特征速度改变。给定其他参数,变大将使下游节流程度增加,进气道处于亚临界状态;变小将使下游节流程度减小,进气道处于超临界状态。亚临界工作状态会使进气道产生低频高振幅的“喘振”,可能引起进气道结构破坏、发动机停车甚至烧坏。超临界工作状态会使进气道产生高频低振幅的“轰鸣”,相对于“喘振”,“轰鸣”的危害较小。进气道流量系数关系到进气道乃至整个发动机工作的稳定性。
由流量系数可求得进气道入口面积:
按照上述内容,空燃比、补燃室压力、进气道参数、燃气发生器参数计算示例如下:
空燃比:若取高度H=10km,飞行马赫数Ma=2.5,补燃室压力0.5MPa,出口压力p e =13499.9Pa作为已知的设计点参数。推进剂选取22%HTPB、38% AP、18%Al和22%Mg组成的中能贫氧推进剂配方。通过热力计算得到理论燃烧温度与理论比冲随空燃比变化曲线请参见图8。
补燃室压力:若推进剂组分为22%HTPB、38% AP、18%Al和22%Mg,飞行马赫数为2.5,高度10km,则可估计来流总压约为1.48MPa。补燃室压力为0.85MPa时,应在空燃比3~13间选取,这里选取设计空燃比为12。
进气道参数:前体采用三波系压缩方案,压缩角分别为7,8,8度,总转角为23度;外壳角度分别为6,17度,总转角为23度;亚扩段采用圆弧曲线设计,出口处平均气流方向与轴线夹角为45度;隔道高度为12mm;内收缩比为1.03。
燃气发生器中的推进剂主要性能指标如下:燃气流量:≥1.0kg/s;工作时间:≥90s;最大纵向过载:20g;燃速:r=17 mm/s。
4、根据导弹质量模型、导弹气动模型、弹道参数模型和冲压发动机性能模型构建导弹总体设计模型,并按照1步骤设置的总体设计指标和约束条件进行求解。
其中,设计变量Z如表6所示。
在表6中,Z被替换为X。
表6
其中,通过求解器和灵敏度分析器工具进行求解,求解器包括了成熟的遗传算法、粒子群算法和进化算法,可得到导弹总体设计模型的最优解。灵敏度分析器包括正交试验设计等算法,可确定出不同精度的导弹总体设计模型,即低精度模型(第二模型)和高精度模型(第一模型)。
若以射程最大为优化目标,则求解过程可以为:
利用优化算法首先对低精度模型进行求解,得到其最优值作为高精度模型的初始迭代值,对高精度模型再次利用优化算法进行求解,得到最终的最优解。其中,低精度指模型的计算精度较低,运算速度较快,易于缩小求解空间;高精度模型的计算精度较高,但计算速度慢。本实例采用的优化算法为遗传算法。
请参见表7。其中低精度多学科优化迭代了262代,与初始射程参数相比,射程增加了189.42km,将低精度的最优值对应的参数作为中可信度多学科优化的迭代初值,经过189代的优化迭代,与初始射程参数相比,射程增加了8.23km。可见此求解过程的优化程度较佳。
表7
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (9)
1.一种飞行器总体设计参数的确定方法,其特征在于,包括:
确定飞行器的各个组成部分和每个组成部分包括的设计参数,所述飞行器由冲压发动机提供动力;
根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;
根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数;
其中,所述参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,所述冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数;
其中,所述按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数,包括:
按照灵敏度分析所述飞行器总体设计模型,得到第一模型和第二模型;
当所述优化目标为射程最大时,利用遗传算法求解所述第二模型,得到第一最优解;
将所述第一最优解作为所述第一模型的初始迭代值,并利用遗传算法求解所述第一模型,得到第二最优解;
将所述第二最优解作为射程最大对应的飞行器总体设计参数;
其中,所述第一模型的计算精度高于所述第二模型的计算精度。
5.根据权利要求1-4任意一项所述的飞行器总体设计参数的确定方法,其特征在于,所述冲压发动机性能模型的构建步骤包括:
确定所述冲压发动机的组成部分,所述冲压发动机的组成部分包括:燃气发生器、补燃室和进气道,所述进气道为两个二元进气道;
根据所述冲压发动机的组成部分计算发动机推力、发动机比冲、所述空燃比和所述补燃室压力;
构建包括所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标的所述冲压发动机性能模型;所述冲压发动机性能模型为:
7.一种飞行器总体设计参数的确定装置,其特征在于,包括:
确定模块,用于确定飞行器的各个组成部分和每个组成部分包括的设计参数,所述飞行器由冲压发动机提供动力;
构建模块,用于根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型;
执行模块,用于根据所有参数模型确定飞行器总体设计模型,并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型,得到所述优化目标对应的飞行器总体设计参数;
其中,所述参数模型至少包括:飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型,所述冲压发动机性能模型中的设计参数至少包括:空燃比、补燃室压力、进气道参数和燃气发生器参数;
其中,所述执行模块具体用于:
按照灵敏度分析所述飞行器总体设计模型,得到第一模型和第二模型;
当所述优化目标为射程最大时,利用遗传算法求解所述第二模型,得到第一最优解;
将所述第一最优解作为所述第一模型的初始迭代值,并利用遗传算法求解所述第一模型,得到第二最优解;
将所述第二最优解作为射程最大对应的飞行器总体设计参数;
其中,所述第一模型的计算精度高于所述第二模型的计算精度。
8.一种飞行器总体设计参数的确定设备,其特征在于,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如权利要求1-6任意一项所述的飞行器总体设计参数的确定方法的步骤。
9.一种可读存储介质,其特征在于,所述可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1-6任意一项所述的飞行器总体设计参数的确定方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911001494.4A CN110516410B (zh) | 2019-10-21 | 2019-10-21 | 一种飞行器总体设计参数的确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911001494.4A CN110516410B (zh) | 2019-10-21 | 2019-10-21 | 一种飞行器总体设计参数的确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110516410A CN110516410A (zh) | 2019-11-29 |
CN110516410B true CN110516410B (zh) | 2020-01-31 |
Family
ID=68634347
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911001494.4A Active CN110516410B (zh) | 2019-10-21 | 2019-10-21 | 一种飞行器总体设计参数的确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110516410B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111173644B (zh) * | 2020-01-06 | 2021-06-04 | 中国人民解放军国防科技大学 | 一种冲压发动机性能确定方法 |
CN112124620A (zh) * | 2020-09-13 | 2020-12-25 | 中国运载火箭技术研究院 | 应用于飞行器的伺服负载力矩计算方法、设备及存储介质 |
CN112362290B (zh) * | 2020-09-30 | 2021-08-03 | 成都飞机工业(集团)有限责任公司 | 一种机翼厚度公差对阻力系数影响的快速分析方法 |
CN112464364A (zh) * | 2020-11-19 | 2021-03-09 | 北京航空航天大学 | 飞行器设计框架生成方法及装置 |
CN113759956B (zh) * | 2020-12-14 | 2024-05-28 | 北京天兵科技有限公司 | 一种用于亚轨道飞行器的飞行弹道设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930550A (zh) * | 2016-04-01 | 2016-09-07 | 方洋旺 | 一种吸气式高超声速导弹助推-跳跃式弹道优化方法 |
CN107609307A (zh) * | 2017-10-10 | 2018-01-19 | 北京理工大学 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
CN110147521A (zh) * | 2019-04-25 | 2019-08-20 | 北京航空航天大学 | 一种高超声速飞行器跳跃滑翔弹道解析求解方法 |
-
2019
- 2019-10-21 CN CN201911001494.4A patent/CN110516410B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930550A (zh) * | 2016-04-01 | 2016-09-07 | 方洋旺 | 一种吸气式高超声速导弹助推-跳跃式弹道优化方法 |
CN107609307A (zh) * | 2017-10-10 | 2018-01-19 | 北京理工大学 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
CN110147521A (zh) * | 2019-04-25 | 2019-08-20 | 北京航空航天大学 | 一种高超声速飞行器跳跃滑翔弹道解析求解方法 |
Non-Patent Citations (3)
Title |
---|
Research on control-oriented coupling modeling for air-breathing hypersonic propulsion systems;DongZhang 等;《Aerospace Science and Technology》;20190131;第84卷;第143-157页 * |
Trajectory optimization for a ramjet-powered vehicle in ascent phase via the Gauss pseudospectral method;ShengboYang 等;《Aerospace Science and Technology》;20170831;第67卷;第88-95页 * |
冲压发动机为动力导弹爬升弹道优化;王华 等;《弹箭与制导学报》;20080630;第28卷(第3期);第185-188页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110516410A (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110516410B (zh) | 一种飞行器总体设计参数的确定方法及装置 | |
CN112613119B (zh) | 一种多几何参数可调的进/排/发一体化航空推进系统建模方法 | |
CN111339681A (zh) | 一种采用空气介质模拟发动机燃气介质喷流气动干扰效应的喷管出口参数匹配方法 | |
CN105975651B (zh) | 基于遗传粒子群多学科设计优化算法的导弹参数设计方法 | |
CN115688287A (zh) | 一种带引射喷管航空涡扇发动机的设计方法 | |
Gamble et al. | Improving off-design nozzle performance using fluidic injection | |
Naveen et al. | Design optimization of nozzle and second throat diffuser system for high altitude test using CFD | |
CN114934861A (zh) | 一种适用于低动压飞行的火箭增强型冲压发动机 | |
CN117993261A (zh) | 航空发动机整机喘振动态过程仿真计算方法 | |
Mishra et al. | Aerodynamic thrust characteristics and performance valuations for an aerospike nozzle-review | |
Hao et al. | A new design method for mode transition control law of variable cycle engine | |
CN114626146B (zh) | 一种用于组合动力系统的仿真方法 | |
Shi et al. | Numerical study of a boundary layer bleedfor a rocket-based combined-cycle inlet in ejector mode | |
CN115982943A (zh) | 一种超燃冲压发动机容积动力学一维建模方法 | |
Nordqvist et al. | Conceptual design of a turbofan engine for a supersonic business jet | |
Safdar et al. | Numerical Analysis of Afterburner Characteristics of a Low Bypass Ratio Turbofan Engine at Various Flight Conditions | |
CN115653784A (zh) | 水冲压发动机增压进水系统 | |
Christensen | Air turborocket/vehicle performance comparison | |
Patnaik et al. | Optimization for aircraft engines with regression and neural-network analysis approximators | |
Yu et al. | Experimental and numerical investigation on design parameters of slot nozzles | |
JP2002054503A (ja) | 亜音速機用高バイパス比・可変サイクルエンジン | |
Seidler et al. | Impact of Shape-Adaptive Fan Blading on the Performance of a Turbofan Engine | |
Ostrander et al. | Inlet analysis for ATR-powered supersonic missile | |
C¸ evik et al. | A robust design optimization of a mixed-flow compressor impeller | |
Kobayashi et al. | Optimal design of hypersonic turbojet engines for two-stage-to-orbit spaceplane |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |