CN110516410B - 一种飞行器总体设计参数的确定方法及装置 - Google Patents

一种飞行器总体设计参数的确定方法及装置 Download PDF

Info

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
parameters
parameter
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
Application number
CN201911001494.4A
Other languages
English (en)
Other versions
CN110516410A (zh
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.)
Xingshan Zhilian Technology Shenzhen Co ltd
Original Assignee
Hunan Deyakun Creative Technology Co Ltd
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 Hunan Deyakun Creative Technology Co Ltd filed Critical Hunan Deyakun Creative Technology Co Ltd
Priority to CN201911001494.4A priority Critical patent/CN110516410B/zh
Publication of CN110516410A publication Critical patent/CN110516410A/zh
Application granted granted Critical
Publication of CN110516410B publication Critical patent/CN110516410B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

一种飞行器总体设计参数的确定方法及装置
技术领域
本发明涉及飞行器设计技术领域，更具体地说，涉及一种飞行器总体设计参数的确定方法、装置、设备及可读存储介质。
背景技术
飞行器设计中涉及到许多学科，从整体系统的角度来说，各个子系统间必然存在一定程度的相互作用和耦合效应。由于子系统均有各自独立的分析和设计工具，传统次序型设计往往无法综合考虑，从而没能充分考虑到各系统间的耦合因素。应用多学科设计优化解决上述整体系统的设计问题，既提高了设计效率，还改善了产品的性能。
在现有技术中，在进行多学科设计优化时，未考虑到动力参数和弹道参数之间的耦合效应，从而导致当飞行器的飞行空域和速域变化明显时，动力参数和弹道参数的耦合特征明显。当采用冲压发动机为飞行器提供动力时，若仍按照现有技术的设计方式设计总体参数，由于未考虑动力参数和弹道参数的耦合效应，会使得飞行器的总体性能降低。
因此，如何降低动力参数和弹道参数的耦合效应，提高飞行器的总体性能，是本领域技术人员需要解决的问题。
发明内容
本发明的目的在于提供一种飞行器总体设计参数的确定方法、装置、设备及可读存储介质，以降低动力参数和弹道参数的耦合效应，提高飞行器的总体性能。
为实现上述目的，本发明实施例提供了如下技术方案：
一种飞行器总体设计参数的确定方法，包括：
确定飞行器的各个组成部分和每个组成部分包括的设计参数，所述飞行器由冲压发动机提供动力；
根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型；
根据所有参数模型确定飞行器总体设计模型，并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数；
其中，所述参数模型至少包括：飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，所述冲压发动机性能模型中的设计参数至少包括：空燃比、补燃室压力、进气道参数和燃气发生器参数；
其中，所述按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数，包括：
按照灵敏度分析所述飞行器总体设计模型，得到第一模型和第二模型；
当所述优化目标为射程最大时，利用遗传算法求解所述第二模型，得到第一最优解；
将所述第一最优解作为所述第一模型的初始迭代值，并利用遗传算法求解所述第一模型，得到第二最优解；
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
其中，所述飞行器质量模型的构建步骤包括：
确定所述飞行器的弹体组成部分和每个弹体组成部分的质量，并根据每个弹体组成部分的质量构建所述飞行器质量模型；
所述飞行器质量模型为：
Figure 267271DEST_PATH_IMAGE001
，
Figure 22738DEST_PATH_IMAGE002
表示飞行器的质量，表示采用等效厚度的壳体模型的机体质量，
Figure 430902DEST_PATH_IMAGE004
表示发动机的质量，
Figure 390899DEST_PATH_IMAGE005
表示有效载荷质量，
Figure 950056DEST_PATH_IMAGE006
表示弹翼或机翼的质量。
其中，所述飞行器气动模型的构建步骤包括：
采用部件组合法构建所述飞行器气动模型；所述飞行器气动模型为：
Figure 662797DEST_PATH_IMAGE007
其中，表示阻力系数；
Figure 17872DEST_PATH_IMAGE009
表示升力系数，表示空气密度，S表示参考面积，q表示动压头，X表示气动阻力，Y表示气动升力，V表示飞行速度。
其中，所述弹道参数模型的构建步骤包括：
确定所述飞行器的弹道设计参数，并采用三自由度质点构建所述弹道参数模型，其中，所述弹道参数模型中的纵向运动方程组为：
Figure 964280DEST_PATH_IMAGE011
其中，V表示飞行速度，m表示飞行器质量，F表示发动机推力，X表示气动阻力，Y表示气动升力，
Figure 856012DEST_PATH_IMAGE012
表示俯仰力矩，
Figure 977552DEST_PATH_IMAGE013
表示转动惯量，
Figure 878512DEST_PATH_IMAGE014
表示攻角，表示速度倾角，
Figure 379212DEST_PATH_IMAGE016
表示俯仰角，
Figure 519206DEST_PATH_IMAGE017
表示俯仰角速度，
Figure 223857DEST_PATH_IMAGE018
表示单位时间内燃料消耗量，
Figure 513499DEST_PATH_IMAGE019
表示飞行器射程，
Figure 12614DEST_PATH_IMAGE020
表示飞行器高度，t表示时间，g表示重力加速度。
其中，所述冲压发动机性能模型的构建步骤包括：
确定所述冲压发动机的组成部分，所述冲压发动机的组成部分包括：燃气发生器、补燃室和进气道，所述进气道为两个二元进气道；
根据所述冲压发动机的组成部分计算发动机推力、发动机比冲、所述空燃比和所述补燃室压力；
构建包括所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标的所述冲压发动机性能模型；所述冲压发动机性能模型为：
Figure 905484DEST_PATH_IMAGE021
其中，F表示发动机推力，表示燃料流量，所述燃料流量为所述燃气发生器参数，v a 表示来流空气速度，v e 表示喷管出口速度，p e 表示喷管出口静压，p a 表示来流空气静压，A e 表示喷管出口面积；I s 表示发动机比冲；
Figure 295325DEST_PATH_IMAGE023
表示空燃比，
Figure 168603DEST_PATH_IMAGE024
表示空气流量；
Figure 283189DEST_PATH_IMAGE025
表示进气道总压恢复系数，所述进气道总压恢复系数为所述进气道参数，
Figure 64064DEST_PATH_IMAGE026
表示进气道出口总压，表示来流空气总压；表示补燃室压力，
Figure 39607DEST_PATH_IMAGE029
表示补燃室总压恢复系数。
其中，所述按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数，包括：
按照灵敏度分析所述飞行器总体设计模型，得到第一模型和第二模型；
当所述优化目标为射程最大时，利用遗传算法求解所述第二模型，得到第一最优解；
将所述第一最优解作为所述第一模型的初始迭代值，并利用遗传算法求解所述第一模型，得到第二最优解；
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
其中，所述根据所有参数模型确定飞行器总体设计模型，包括：
根据所述飞行器质量模型、所述飞行器气动模型、所述弹道参数模型和所述冲压发动机性能模型确定所述飞行器总体设计模型；
其中，所述飞行器总体设计模型的设计变量为：发动机动力参数和飞行器几何参数，弹道形式为高空-俯冲；
所述飞行器总体设计模型的约束条件为：
Figure 827434DEST_PATH_IMAGE030
，
Figure 198373DEST_PATH_IMAGE031
和
Figure 210191DEST_PATH_IMAGE032
，m表示飞行器质量，L表示飞行器长度，Z指代所有设计变量。
一种飞行器总体设计参数的确定装置，包括：
确定模块，用于确定飞行器的各个组成部分和每个组成部分包括的设计参数，所述飞行器由冲压发动机提供动力；
构建模块，用于根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型；
执行模块，用于根据所有参数模型确定飞行器总体设计模型，并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数；
其中，所述参数模型至少包括：飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，所述冲压发动机性能模型中的设计参数至少包括：空燃比、补燃室压力、进气道参数和燃气发生器参数；
其中，所述执行模块具体用于：
按照灵敏度分析所述飞行器总体设计模型，得到第一模型和第二模型；
当所述优化目标为射程最大时，利用遗传算法求解所述第二模型，得到第一最优解；
将所述第一最优解作为所述第一模型的初始迭代值，并利用遗传算法求解所述第一模型，得到第二最优解；
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
一种飞行器总体设计参数的确定设备，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现上述任意一项所述的飞行器总体设计参数的确定方法的步骤。
一种可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述任意一项所述的飞行器总体设计参数的确定方法的步骤。
通过以上方案可知，本发明实施例提供的一种飞行器总体设计参数的确定方法，包括：确定飞行器的各个组成部分和每个组成部分包括的设计参数，所述飞行器由冲压发动机提供动力；根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型；根据所有参数模型确定飞行器总体设计模型，并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数；其中，所述参数模型至少包括：飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，所述冲压发动机性能模型中的设计参数至少包括：空燃比、补燃室压力、进气道参数和燃气发生器参数。其中，按照预设的优化目标和遗传算法求解飞行器总体设计模型，得到优化目标对应的飞行器总体设计参数，包括：按照灵敏度分析飞行器总体设计模型，得到第一模型和第二模型；当优化目标为射程最大时，利用遗传算法求解第二模型，得到第一最优解；将第一最优解作为第一模型的初始迭代值，并利用遗传算法求解第一模型，得到第二最优解；将第二最优解作为射程最大对应的飞行器总体设计参数。
可见，所述方法在确定飞行器总体设计模型时，使飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型均参与其中，且冲压发动机性能模型中的设计参数包括进气道参数，也就是说，该方法在确定飞行器总体设计模型时，考虑了所有模型之间的耦合效应。
具体的，本申请按照灵敏度分析飞行器总体设计模型，得到第一模型和第二模型；当所述优化目标为射程最大时，利用遗传算法求解第二模型，得到第一最优解；将第一最优解作为第一模型的初始迭代值，并利用遗传算法求解第一模型，得到第二最优解；将第二最优解作为射程最大对应的飞行器总体设计参数。如此则利用遗传算法求解得到了射程最大对应的总体设计参数。由于本申请提供的飞行器总体设计模型中涵盖了飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，因此遗传算法在求解过程中能够考虑动力参数和弹道参数的耦合效应，降低动力参数和弹道参数的耦合，从而提高飞行器的总体性能。
相应地，本发明实施例提供的一种飞行器总体设计参数的确定装置、设备及可读存储介质，也同样具有上述技术效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例公开的一种飞行器总体设计参数的确定方法流程图；
图2为本发明实施例公开的一种飞行器总体设计参数的确定装置示意图；
图3为本发明实施例公开的一种飞行器总体设计参数的确定设备示意图；
图4为本发明实施例公开的导弹总体设计参数的设计总流程示意图；
图5为本发明实施例公开的导弹弹体的各个组成部分示意图；
图6为本发明实施例公开的升力系数随马赫数变化示意图；
图7为本发明实施例公开的阻力系数随马赫数变化示意图；
图8为本发明实施例公开理论燃烧温度与理论比冲随空燃比变化曲线示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例公开了一种飞行器总体设计参数的确定方法、装置、设备及可读存储介质，以降低动力参数和弹道参数的耦合效应，提高飞行器的总体性能。
参见图1，本发明实施例提供的一种飞行器总体设计参数的确定方法，包括：
S101、确定飞行器的各个组成部分和每个组成部分包括的设计参数，飞行器由冲压发动机提供动力；
S102、根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型；
S103、根据所有参数模型确定飞行器总体设计模型，并按照预设的优化目标和遗传算法求解飞行器总体设计模型，得到优化目标对应的飞行器总体设计参数；
其中，参数模型至少包括：飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，冲压发动机性能模型中的设计参数至少包括：空燃比、补燃室压力、进气道参数和燃气发生器参数。
其中，所述飞行器质量模型的构建步骤包括：确定所述飞行器的弹体组成部分和每个弹体组成部分的质量，并根据每个弹体组成部分的质量构建所述飞行器质量模型。
其中，所述飞行器质量模型为：
Figure 784523DEST_PATH_IMAGE001
，
Figure 172779DEST_PATH_IMAGE002
表示飞行器的质量，
Figure 398224DEST_PATH_IMAGE003
表示采用等效厚度的壳体模型的机体质量，
Figure 846523DEST_PATH_IMAGE004
表示发动机的质量，表示有效载荷质量，
Figure 120606DEST_PATH_IMAGE006
表示弹翼或机翼的质量。
其中，所述飞行器气动模型的构建步骤包括：
采用部件组合法构建所述飞行器气动模型；所述飞行器气动模型为：
Figure 200557DEST_PATH_IMAGE007
其中，表示阻力系数；表示升力系数，
Figure 302002DEST_PATH_IMAGE010
表示空气密度，S表示参考面积，q表示动压头，X表示气动阻力，Y表示气动升力，V表示飞行速度。
其中，所述弹道参数模型的构建步骤包括：确定所述飞行器的弹道设计参数，并采用三自由度质点构建所述弹道参数模型，其中，所述弹道参数模型中的纵向运动方程组为：
其中，V表示飞行速度，m表示飞行器质量，F表示发动机推力，X表示气动阻力，Y表示气动升力，
Figure 26562DEST_PATH_IMAGE012
表示俯仰力矩，表示转动惯量，
Figure 596532DEST_PATH_IMAGE014
表示攻角，
Figure 651075DEST_PATH_IMAGE015
表示速度倾角，
Figure 815340DEST_PATH_IMAGE016
表示俯仰角，
Figure 322545DEST_PATH_IMAGE017
表示俯仰角速度，
Figure 659986DEST_PATH_IMAGE018
表示单位时间内燃料消耗量，
Figure 585347DEST_PATH_IMAGE019
表示飞行器射程，
Figure 451672DEST_PATH_IMAGE020
表示飞行器高度，t表示时间，g表示重力加速度。
其中，所述冲压发动机性能模型的构建步骤包括：确定所述冲压发动机的组成部分，所述冲压发动机的组成部分包括：燃气发生器、补燃室和进气道，所述进气道为两个二元进气道；根据所述冲压发动机的组成部分计算发动机推力、发动机比冲、所述空燃比和所述补燃室压力；根据所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标构建所述冲压发动机性能模型，也就是构建包括所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标的所述冲压发动机性能模型。
其中，所述冲压发动机性能模型为：
Figure 915015DEST_PATH_IMAGE021
其中，F表示发动机推力，
Figure 524987DEST_PATH_IMAGE022
表示燃料流量，所述燃料流量为所述燃气发生器参数，v a 表示来流空气速度，v e 表示喷管出口速度，p e 表示喷管出口静压，p a 表示来流空气静压，A e 表示喷管出口面积；I s 表示发动机比冲；
Figure 554123DEST_PATH_IMAGE023
表示空燃比，
Figure 339152DEST_PATH_IMAGE024
表示空气流量；
Figure 352108DEST_PATH_IMAGE025
表示进气道总压恢复系数，所述进气道总压恢复系数为所述进气道参数，表示进气道出口总压，
Figure 649414DEST_PATH_IMAGE027
表示来流空气总压；
Figure 326383DEST_PATH_IMAGE028
表示补燃室压力，
Figure 46208DEST_PATH_IMAGE029
表示补燃室总压恢复系数。
需要说明的是，构建冲压发动机性能模型的目的为：基于冲压发动机的技术指标、空燃比、补燃室压力、进气道参数、燃气发生器参数等各方参数，求解发动机推力和发动机比冲。发动机推力和发动机比冲可表示冲压发动机的性能；发动机推力越大，发动机比冲越高，那么发动机性能就越高。发动机推力和发动机比冲的具体计算请参见下文描述的相应计算公式。发动机比冲即为理论比冲。
可见，本实施例提供了一种飞行器总体设计参数的确定方法，所述方法在确定飞行器总体设计模型时，使飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型均参与其中，且冲压发动机性能模型中的设计参数包括进气道参数，也就是说，该方法在确定飞行器总体设计模型时，考虑了所有模型之间的耦合效应，利用遗传算法求解得到了射程最大对应的总体设计参数。由于本申请提供的飞行器总体设计模型中涵盖了飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，因此遗传算法在求解过程中能够考虑动力参数和弹道参数的耦合效应，降低动力参数和弹道参数的耦合，从而提高飞行器的总体性能。
基于上述实施例，需要说明的是，所述按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数，包括：
按照灵敏度分析所述飞行器总体设计模型，得到第一模型和第二模型；
当所述优化目标为射程最大时，利用遗传算法求解所述第二模型，得到第一最优解；
将所述第一最优解作为所述第一模型的初始迭代值，并利用遗传算法求解所述第一模型，得到第二最优解；
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
需要说明的是，第一模型和第二模型的计算精度是相对而言的，第一模型的计算速度比第二模型的计算速度慢，但其精度比第二模型的精度高。也就是说，按照灵敏度分析飞行器总体设计模型，可得到具有不同计算精度的第一模型和第二模型。在具体计算时，首先用遗传算法求解较低精度的第二模型，得到第一最优解；进而将第一最优解作为较高精度的第一模型的初始迭代值，并利用遗传算法求解第一模型，从而得到第二最优解，那么第二最优解即为射程最大对应的飞行器总体设计参数。
基于上述实施例，需要说明的是，所述根据所有参数模型确定飞行器总体设计模型，包括：
根据所述飞行器质量模型、所述飞行器气动模型、所述弹道参数模型和所述冲压发动机性能模型确定所述飞行器总体设计模型；
其中，所述飞行器总体设计模型的设计变量为：发动机动力参数和飞行器几何参数，弹道形式为高空-俯冲；
所述飞行器总体设计模型的约束条件为：
Figure 997984DEST_PATH_IMAGE030
，和
Figure 380741DEST_PATH_IMAGE032
，m表示飞行器质量，L表示飞行器长度，Z指代所有设计变量。
下面对本发明实施例提供的一种飞行器总体设计参数的确定装置进行介绍，下文描述的一种飞行器总体设计参数的确定装置与上文描述的一种飞行器总体设计参数的确定方法可以相互参照。
参见图2，本发明实施例提供的一种飞行器总体设计参数的确定装置，包括：
确定模块201，用于确定飞行器的各个组成部分和每个组成部分包括的设计参数，所述飞行器由冲压发动机提供动力；
构建模块202，用于根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型；
执行模块203，用于根据所有参数模型确定飞行器总体设计模型，并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数；
其中，所述参数模型至少包括：飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，所述冲压发动机性能模型中的设计参数至少包括：空燃比、补燃室压力、进气道参数和燃气发生器参数。
其中，所述构建模块具体用于：
确定所述飞行器的弹体组成部分和每个弹体组成部分的质量，并根据每个弹体组成部分的质量构建所述飞行器质量模型；
所述飞行器质量模型为：
Figure 853441DEST_PATH_IMAGE001
，
Figure 608908DEST_PATH_IMAGE002
表示飞行器的质量，表示采用等效厚度的壳体模型的机体质量，
Figure 282652DEST_PATH_IMAGE004
表示发动机的质量，表示有效载荷质量，
Figure 473910DEST_PATH_IMAGE006
表示弹翼或机翼的质量。
其中，所述构建模块具体用于：
采用部件组合法构建所述飞行器气动模型；所述飞行器气动模型为：
Figure 452230DEST_PATH_IMAGE007
其中，
Figure 438640DEST_PATH_IMAGE008
表示阻力系数；
Figure 869622DEST_PATH_IMAGE009
表示升力系数，
Figure 449114DEST_PATH_IMAGE010
表示空气密度，S表示参考面积，q表示动压头，X表示气动阻力，Y表示气动升力，V表示飞行速度。
其中，所述构建模块具体用于：
确定所述飞行器的弹道设计参数，并采用三自由度质点构建所述弹道参数模型，其中，所述弹道参数模型中的纵向运动方程组为：
Figure 813100DEST_PATH_IMAGE011
其中，V表示飞行速度，m表示飞行器质量，F表示发动机推力，X表示气动阻力，Y表示气动升力，
Figure 704832DEST_PATH_IMAGE012
表示俯仰力矩，
Figure 639421DEST_PATH_IMAGE013
表示转动惯量，表示攻角，
Figure 962135DEST_PATH_IMAGE015
表示速度倾角，
Figure 290348DEST_PATH_IMAGE016
表示俯仰角，
Figure 633605DEST_PATH_IMAGE017
表示俯仰角速度，
Figure 88988DEST_PATH_IMAGE018
表示单位时间内燃料消耗量，表示飞行器射程，
Figure 395522DEST_PATH_IMAGE020
表示飞行器高度，t表示时间，g表示重力加速度。
其中，所述构建模块具体用于：
确定所述冲压发动机的组成部分，所述冲压发动机的组成部分包括：燃气发生器、补燃室和进气道，所述进气道为两个二元进气道；
根据所述冲压发动机的组成部分计算发动机推力、发动机比冲、所述空燃比和所述补燃室压力；
构建包括所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标的所述冲压发动机性能模型；所述冲压发动机性能模型为：
Figure 22812DEST_PATH_IMAGE021
其中，F表示发动机推力，
Figure 750728DEST_PATH_IMAGE022
表示燃料流量，所述燃料流量为所述燃气发生器参数，v a 表示来流空气速度，v e 表示喷管出口速度，p e 表示喷管出口静压，p a 表示来流空气静压，A e 表示喷管出口面积；I s 表示发动机比冲；
Figure 412653DEST_PATH_IMAGE023
表示空燃比，
Figure 82669DEST_PATH_IMAGE024
表示空气流量；
Figure 197256DEST_PATH_IMAGE025
表示进气道总压恢复系数，所述进气道总压恢复系数为所述进气道参数，
Figure 181392DEST_PATH_IMAGE026
表示进气道出口总压，
Figure 242365DEST_PATH_IMAGE027
表示来流空气总压；
Figure 286544DEST_PATH_IMAGE028
表示补燃室压力，
Figure 685164DEST_PATH_IMAGE029
表示补燃室总压恢复系数。
其中，所述执行模块具体用于：
按照灵敏度分析所述飞行器总体设计模型，得到第一模型和第二模型；
当所述优化目标为射程最大时，利用遗传算法求解所述第二模型，得到第一最优解；将所述第一最优解作为所述第一模型的初始迭代值，并利用遗传算法求解所述第一模型，得到第二最优解；
将所述第二最优解作为射程最大对应的飞行器总体设计参数。
其中，所述执行模块具体用于：
根据所述飞行器质量模型、所述飞行器气动模型、所述弹道参数模型和所述冲压发动机性能模型确定所述飞行器总体设计模型；
其中，所述飞行器总体设计模型的设计变量为：发动机动力参数和飞行器几何参数，弹道形式为高空-俯冲；
所述飞行器总体设计模型的约束条件为：
Figure 20462DEST_PATH_IMAGE030
，
Figure 594663DEST_PATH_IMAGE031
和
Figure 606481DEST_PATH_IMAGE032
，m表示飞行器质量，L表示飞行器长度，Z指代所有设计变量。
可见，本实施例提供了一种飞行器总体设计参数的确定装置，包括：确定模块、构建模块以及执行模块。首先由确定模块确定飞行器的各个组成部分和每个组成部分包括的设计参数，所述飞行器由冲压发动机提供动力；然后构建模块根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型；最后执行模块根据所有参数模型确定飞行器总体设计模型，并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数；其中，所述参数模型至少包括：飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，所述冲压发动机性能模型中的设计参数至少包括：空燃比、补燃室压力、进气道参数和燃气发生器参数。如此各个模块之间分工合作，各司其职，从而降低了动力参数和弹道参数的耦合效应，提高了飞行器的总体性能。
下面对本发明实施例提供的一种飞行器总体设计参数的确定设备进行介绍，下文描述的一种飞行器总体设计参数的确定设备与上文描述的一种飞行器总体设计参数的确定方法及装置可以相互参照。
参见图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
Figure 430080DEST_PATH_IMAGE033
其中，冲压发动机的进气道参数如表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）构建导弹质量模型：导弹的结构质量是将导弹几何模型进行一定厚度的抽壳后，对此壳体赋予密度等性质，进而计算出来其质量、质心参数。弹上设备质量设定为固定值，每一个设备的质量赋值为给定值，假设每一个设备都为一均质长方体，其中心位置即为设备质心位置。
请参见图5，图5为导弹弹体的各个组成部分示意图。所述飞行器质量模型为：，
Figure 242813DEST_PATH_IMAGE002
表示飞行器的质量，
Figure 22550DEST_PATH_IMAGE003
表示采用等效厚度的壳体模型的机体质量，
Figure 948918DEST_PATH_IMAGE004
表示发动机的质量，
Figure 28869DEST_PATH_IMAGE005
表示有效载荷质量，
Figure 664381DEST_PATH_IMAGE006
表示弹翼或机翼的质量。
2）构建导弹气动模型：采用部件组合法计算阻力系数和升力系数，并根据计算得到的阻力系数和升力系数构建导弹气动模型；阻力系数的计算公式为：
Figure 462573DEST_PATH_IMAGE035
，升力系数的计算公式为：
Figure 927052DEST_PATH_IMAGE036
；其中，，
Figure 120453DEST_PATH_IMAGE010
表示空气密度，S表示参考面积，q表示动压头，X表示气动阻力，Y表示气动升力，V表示飞行速度。
其中，升力系数随马赫数变化示意图请参见图6，阻力系数随马赫数变化示意图请参见图7。
3）构建弹道参数模型：确定导弹的弹道设计参数，并采用三自由度质点构建弹道参数模型，其中，弹道参数模型中的纵向运动方程组为：
Figure 888164DEST_PATH_IMAGE011
其中，V表示飞行速度，m表示飞行器质量，F表示发动机推力，X表示气动阻力，Y表示气动升力，
Figure 687493DEST_PATH_IMAGE012
表示俯仰力矩，表示转动惯量，
Figure 437460DEST_PATH_IMAGE014
表示攻角，
Figure 960977DEST_PATH_IMAGE015
表示速度倾角，表示俯仰角，
Figure 145150DEST_PATH_IMAGE017
表示俯仰角速度，
Figure 277054DEST_PATH_IMAGE018
表示单位时间内燃料消耗量，
Figure 537134DEST_PATH_IMAGE019
表示飞行器射程，
Figure 163419DEST_PATH_IMAGE020
表示飞行器高度，t表示时间，g表示重力加速度。
4）构建冲压发动机性能模型：采用二个二元进气道；采用铝镁贫氧推进剂；采用非壅塞式燃气发生器；其中，冲压发动机主要由进气道、燃气发生器和补燃室（包括冲压喷管）组成。
其中，所述冲压发动机性能模型为：
Figure 192555DEST_PATH_IMAGE021
其中，F表示发动机推力，表示燃料流量，所述燃料流量为所述燃气发生器参数，v a 表示来流空气速度，v e 表示喷管出口速度，p e 表示喷管出口静压，p a 表示来流空气静压，A e 表示喷管出口面积；I s 表示发动机比冲；
Figure 977157DEST_PATH_IMAGE023
表示空燃比，
Figure 594083DEST_PATH_IMAGE024
表示空气流量；
Figure 228458DEST_PATH_IMAGE025
表示进气道总压恢复系数，所述进气道总压恢复系数为所述进气道参数，
Figure 436586DEST_PATH_IMAGE026
表示进气道出口总压，
Figure 405679DEST_PATH_IMAGE027
表示来流空气总压；表示补燃室压力，
Figure 108985DEST_PATH_IMAGE029
表示补燃室总压恢复系数。
分析上式可知，理论比冲单位质量推进剂产生的推力越大，发动机性能越优，而理论比冲与飞行弹道参数（空气速度、空气压强）密切相关。耦合性突出体现在相同的发动机在不同弹道条件下性能差异极大。换句话说，欲使飞行器按期望的高度、速度飞行的距离最远，一定存在满足上述弹道条件的发动机设计参数。
具体的，冲压发动机性能模型的具体构成过程包括：
a、参数定义：
（1）定义空燃比，空燃比即为空气流量与燃料流量之比，用公式表示为：
Figure 475561DEST_PATH_IMAGE039
（2）理论空气量（又称为当量比）定义为1kg燃料完全燃烧所需要的空气量，用公式表示为：
其中，注脚f表示燃料，注脚a表示空气；M a 为空气相对分子质量；M f 为燃料相对分子质量；L为元素数；b i为第i种元素的摩尔数；v i为第i种元素的原子价。
（3）余气系数定义为空燃比与理论空气量之比：
Figure 292525DEST_PATH_IMAGE041
计算平衡状态下燃烧产物的组分必须首先计算反应物假定化学式。固体火箭冲压发动机由空气提供氧化剂，1kg由贫氧推进剂和空气共同组成的反应物的假定化学式由下式计算得到：
Figure 858766DEST_PATH_IMAGE042
其中，
Figure 536872DEST_PATH_IMAGE043
为贫氧推进剂和空气中都含有的元素；
Figure 96030DEST_PATH_IMAGE044
为贫氧推进剂中独有的元素；
Figure 808771DEST_PATH_IMAGE045
为空气中独有的元素；为元素
Figure 242474DEST_PATH_IMAGE047
的摩尔数；
Figure 543005DEST_PATH_IMAGE048
为元素
Figure 110253DEST_PATH_IMAGE047
的相对分子质量，为元素
Figure 670996DEST_PATH_IMAGE047
在空气中的百分比。反应物假定化学式的一般形式为：。
b、性能计算：
比冲既反映了推进剂能量的大小又反映了推进剂在燃烧过程中能量转换的效率。固冲发动机初步设计中，首先要在设计点处找到一组可用、较好的比冲数据，然后根据设计的要求，设计燃料流量和空气流量。燃料流量是设计燃气发生器的重要参数。空气流量是设计进气道的重要参数。
具体的，喷管出口速度v e 表示为：
Figure 524868DEST_PATH_IMAGE051
其中，H c 为补燃室燃烧产物的焓；H e 为喷管出口燃烧产物的焓。
额定推力F由下式计算得到：
Figure 853081DEST_PATH_IMAGE052
其中，p e 为喷管出口静压；p a 为来流空气静压；v a 表示来流空气速度，v e 表示喷管出口速度；A e 为喷管出口面积；
Figure 196338DEST_PATH_IMAGE053
为燃料流量；
Figure 383212DEST_PATH_IMAGE054
为空燃比。
而比冲表示为：
Figure 925052DEST_PATH_IMAGE055
特征速度
Figure 424166DEST_PATH_IMAGE056
表示为：
其中，
Figure 310531DEST_PATH_IMAGE058
为补燃室内燃气总流量；
Figure 972456DEST_PATH_IMAGE059
为补燃室总压；A t 为喉部面积。
故而，由上述公式可得：
其中，下标t表示喷管喉部，下标e表示喷管出口处。故
Figure 491479DEST_PATH_IMAGE061
为喷管喉部燃气流量，
Figure 819824DEST_PATH_IMAGE062
为喷管喉部的平衡因素，p t 表示喷管喉部静压，
Figure 70677DEST_PATH_IMAGE063
为通用气体常数，
Figure 911594DEST_PATH_IMAGE064
为喷管喉部的气相产物摩尔数，T t 为喷管喉部温度；
Figure 998629DEST_PATH_IMAGE065
为喷管出口处燃气流量，
Figure 583194DEST_PATH_IMAGE066
为喷管出口的平衡因素，p e 表示喷管出口静压，T e 为喷管出口温度，p e 表示喷管出口静压。
进一步的，
Figure 954133DEST_PATH_IMAGE067
。
因此比冲可表示为：
c、冲压喷管设计：
冲压喷管的主要设计参数是喉部几何面积和膨胀比
Figure 191189DEST_PATH_IMAGE070
。假设喷管是理想喷管，喷管收敛段无总压损失，则壅塞条件下喷管喉部的有效面积为：
Figure 682213DEST_PATH_IMAGE071
其中，
Figure 864933DEST_PATH_IMAGE022
为燃料流量；
Figure 926561DEST_PATH_IMAGE072
为空燃比；
Figure 56191DEST_PATH_IMAGE073
为补燃室燃烧效率；
Figure 401721DEST_PATH_IMAGE074
为补燃室燃气理论特征速度；
Figure 755342DEST_PATH_IMAGE028
为补燃室总压。
燃气在流过喷管喉部时会产生“颈缩”效应，引入冲压喷管喉部流量系数
Figure 553534DEST_PATH_IMAGE075
来考虑其影响。因此，得到喉部几何面积：
Figure 299904DEST_PATH_IMAGE076
喷管出口静压
Figure 499942DEST_PATH_IMAGE077
可以作为设计值，假设冲压喷管出口气流完全膨胀则
Figure 227726DEST_PATH_IMAGE078
。
由和气动函数
Figure 46963DEST_PATH_IMAGE079
可求得出口速度系数
Figure 852240DEST_PATH_IMAGE080
，则有：
Figure 813242DEST_PATH_IMAGE081
再由气动函数求得喷管膨胀比
Figure 657887DEST_PATH_IMAGE070
，则有：
Figure 338178DEST_PATH_IMAGE083
其中，若给定的设计值是喷管膨胀比
Figure 470082DEST_PATH_IMAGE070
，则可反过来计算出口静压
Figure 667845DEST_PATH_IMAGE077
。
d、进气道设计：
进气道设计计算主要是计算进气道总压恢复系数和流量系数。
（1）总压恢复系数：
进气道总压恢复系数为：
Figure 340135DEST_PATH_IMAGE084
其中，
Figure 854424DEST_PATH_IMAGE029
为补燃室总压恢复系数；
Figure 157229DEST_PATH_IMAGE026
为进气道出口总压；
Figure 639026DEST_PATH_IMAGE029
表示补燃室总压恢复系数。
Figure 52690DEST_PATH_IMAGE029
不变的条件下，当
Figure 687065DEST_PATH_IMAGE025
下降1%时，
Figure 364034DEST_PATH_IMAGE028
同样下降1%。而的下降会导致喷管出口总压减小，这就使冲压喷管完全膨胀所能达到的
Figure 288630DEST_PATH_IMAGE086
减小导致推力下降。进气道总压恢复系数下降1%而推力下降则在1%以上。因此，在设计进气道时应尽量提高总压恢复系数。
在
Figure 683971DEST_PATH_IMAGE027
、
Figure 405939DEST_PATH_IMAGE029
确定的条件下，进气道总压恢复系数随着补燃室总压设计值正相关。但总压恢复系数不能超过进气道的临界总压恢复系数，否则进气道不能与补燃室匹配工作。
（2）流量系数：
进气道流量系数随飞行状态的不同改变。在发动机初步设计中期望进气道流量系数尽可能高，这样发动机在流量调节上会有更大的空间。流量系数直接关系到冲压喷管喉径的设计，具体关系可用下式表示：
Figure 161406DEST_PATH_IMAGE087
其中，
Figure 222902DEST_PATH_IMAGE024
进入补燃室的空气流量；
Figure 772833DEST_PATH_IMAGE022
为燃料流量；为进气道流量系数，
Figure 289057DEST_PATH_IMAGE089
为进入进气道入口的空气流量；
Figure 267377DEST_PATH_IMAGE074
为补燃室燃气理论特征速度。
给定补燃室内燃气理论特征速度，在燃气流量变化时，若喉径过大，则总压
Figure 966660DEST_PATH_IMAGE028
会有一个较大的减小以保证流量的不变；则正激波吸入进气道内，进气道在超临界状态下工作；同时，总压的下降导致推力的下降。若喉径过小，下游节流程度增加，总压
Figure 63929DEST_PATH_IMAGE028
升高，正激波推出进口，产生严重溢流。另外，流量的改变还会使燃料能量的释放发生改变，使得理论特征速度
Figure 100018DEST_PATH_IMAGE074
改变。给定其他参数，
Figure 8062DEST_PATH_IMAGE074
变大将使下游节流程度增加，进气道处于亚临界状态；
Figure 660761DEST_PATH_IMAGE074
变小将使下游节流程度减小，进气道处于超临界状态。亚临界工作状态会使进气道产生低频高振幅的“喘振”，可能引起进气道结构破坏、发动机停车甚至烧坏。超临界工作状态会使进气道产生高频低振幅的“轰鸣”，相对于“喘振”，“轰鸣”的危害较小。进气道流量系数关系到进气道乃至整个发动机工作的稳定性。
由流量系数可求得进气道入口面积：
Figure 827300DEST_PATH_IMAGE090
为了使进气道能稳定的工作，在设计上一般应使进气道在接近临界的超临界工况下工作，以超临界裕度
Figure 514633DEST_PATH_IMAGE091
来定量描述，定义：
Figure 593579DEST_PATH_IMAGE092
其中，
Figure 936835DEST_PATH_IMAGE093
为进气道临界总压恢复系数，
Figure 375907DEST_PATH_IMAGE025
表示进气道总压恢复系数。
按照上述内容，空燃比、补燃室压力、进气道参数、燃气发生器参数计算示例如下：
空燃比：若取高度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
Figure 183326DEST_PATH_IMAGE094
其中，通过求解器和灵敏度分析器工具进行求解，求解器包括了成熟的遗传算法、粒子群算法和进化算法，可得到导弹总体设计模型的最优解。灵敏度分析器包括正交试验设计等算法，可确定出不同精度的导弹总体设计模型，即低精度模型（第二模型）和高精度模型（第一模型）。
若以射程最大为优化目标，则求解过程可以为：
利用优化算法首先对低精度模型进行求解，得到其最优值作为高精度模型的初始迭代值，对高精度模型再次利用优化算法进行求解，得到最终的最优解。其中，低精度指模型的计算精度较低，运算速度较快，易于缩小求解空间；高精度模型的计算精度较高，但计算速度慢。本实例采用的优化算法为遗传算法。
请参见表7。其中低精度多学科优化迭代了262代，与初始射程参数相比，射程增加了189.42km，将低精度的最优值对应的参数作为中可信度多学科优化的迭代初值，经过189代的优化迭代，与初始射程参数相比，射程增加了8.23km。可见此求解过程的优化程度较佳。
表7
Figure 416861DEST_PATH_IMAGE095
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (9)

1.一种飞行器总体设计参数的确定方法，其特征在于，包括：
确定飞行器的各个组成部分和每个组成部分包括的设计参数，所述飞行器由冲压发动机提供动力；
根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型；
根据所有参数模型确定飞行器总体设计模型，并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数；
其中，所述参数模型至少包括：飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，所述冲压发动机性能模型中的设计参数至少包括：空燃比、补燃室压力、进气道参数和燃气发生器参数；
其中，所述按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数，包括：
按照灵敏度分析所述飞行器总体设计模型，得到第一模型和第二模型；
当所述优化目标为射程最大时，利用遗传算法求解所述第二模型，得到第一最优解；
将所述第一最优解作为所述第一模型的初始迭代值，并利用遗传算法求解所述第一模型，得到第二最优解；
将所述第二最优解作为射程最大对应的飞行器总体设计参数；
其中，所述第一模型的计算精度高于所述第二模型的计算精度。
2.根据权利要求1所述的飞行器总体设计参数的确定方法，其特征在于，所述飞行器质量模型的构建步骤包括：
确定所述飞行器的弹体组成部分和每个弹体组成部分的质量，并根据每个弹体组成部分的质量构建所述飞行器质量模型；
所述飞行器质量模型为：
Figure 899105DEST_PATH_IMAGE001
，
Figure 347404DEST_PATH_IMAGE002
表示飞行器的质量，表示采用等效厚度的壳体模型的机体质量，
Figure 787930DEST_PATH_IMAGE004
表示发动机的质量，表示有效载荷质量，
Figure 237814DEST_PATH_IMAGE006
表示弹翼或机翼的质量。
3.根据权利要求2所述的飞行器总体设计参数的确定方法，其特征在于，所述飞行器气动模型的构建步骤包括：
采用部件组合法构建所述飞行器气动模型；所述飞行器气动模型为：
Figure 36006DEST_PATH_IMAGE007
其中，
Figure 31643DEST_PATH_IMAGE008
表示阻力系数；
Figure 231681DEST_PATH_IMAGE009
表示升力系数，
Figure 506935DEST_PATH_IMAGE010
表示空气密度，S表示参考面积，q表示动压头，X表示气动阻力，Y表示气动升力，V表示飞行速度。
4.根据权利要求3所述的飞行器总体设计参数的确定方法，其特征在于，所述弹道参数模型的构建步骤包括：
确定所述飞行器的弹道设计参数，并采用三自由度质点构建所述弹道参数模型，其中，所述弹道参数模型中的纵向运动方程组为：
Figure 792423DEST_PATH_IMAGE011
其中，V表示飞行速度，m表示飞行器质量，F表示发动机推力，X表示气动阻力，Y表示气动升力，
Figure 60593DEST_PATH_IMAGE012
表示俯仰力矩，
Figure 318399DEST_PATH_IMAGE013
表示转动惯量，
Figure 279402DEST_PATH_IMAGE014
表示攻角，
Figure 799989DEST_PATH_IMAGE015
表示速度倾角，
Figure 606271DEST_PATH_IMAGE016
表示俯仰角，表示俯仰角速度，
Figure 912804DEST_PATH_IMAGE018
表示单位时间内燃料消耗量，表示飞行器射程，表示飞行器高度，t表示时间，g表示重力加速度。
5.根据权利要求1-4任意一项所述的飞行器总体设计参数的确定方法，其特征在于，所述冲压发动机性能模型的构建步骤包括：
确定所述冲压发动机的组成部分，所述冲压发动机的组成部分包括：燃气发生器、补燃室和进气道，所述进气道为两个二元进气道；
根据所述冲压发动机的组成部分计算发动机推力、发动机比冲、所述空燃比和所述补燃室压力；
构建包括所述发动机推力、所述发动机比冲、所述空燃比、所述补燃室压力、所述进气道参数、所述燃气发生器参数和所述冲压发动机的技术指标的所述冲压发动机性能模型；所述冲压发动机性能模型为：
Figure 923617DEST_PATH_IMAGE019
其中，F表示发动机推力，
Figure 533590DEST_PATH_IMAGE020
表示燃料流量，所述燃料流量为所述燃气发生器参数，va表示来流空气速度，ve表示喷管出口速度，pe表示喷管出口静压，pa表示来流空气静压，Ae表示喷管出口面积；Is表示发动机比冲；
Figure 765988DEST_PATH_IMAGE021
表示空燃比，
Figure 68793DEST_PATH_IMAGE022
表示空气流量；
Figure 816169DEST_PATH_IMAGE023
表示进气道总压恢复系数，所述进气道总压恢复系数为所述进气道参数，
Figure 980566DEST_PATH_IMAGE024
表示进气道出口总压，
Figure 67470DEST_PATH_IMAGE025
表示来流空气总压；表示补燃室压力，
Figure 510270DEST_PATH_IMAGE027
表示补燃室总压恢复系数。
6.根据权利要求5所述的飞行器总体设计参数的确定方法，其特征在于，所述根据所有参数模型确定飞行器总体设计模型，包括：
根据所述飞行器质量模型、所述飞行器气动模型、所述弹道参数模型和所述冲压发动机性能模型确定所述飞行器总体设计模型；
其中，所述飞行器总体设计模型的设计变量为：发动机动力参数和飞行器几何参数，弹道形式为高空-俯冲；
所述飞行器总体设计模型的约束条件为：
Figure 727625DEST_PATH_IMAGE028
，
Figure 482085DEST_PATH_IMAGE029
和
Figure 329955DEST_PATH_IMAGE030
，m表示飞行器质量，L表示飞行器长度，Z指代所有设计变量。
7.一种飞行器总体设计参数的确定装置，其特征在于，包括：
确定模块，用于确定飞行器的各个组成部分和每个组成部分包括的设计参数，所述飞行器由冲压发动机提供动力；
构建模块，用于根据每个组成部分包括的设计参数分别构建与每个组成部分对应的参数模型；
执行模块，用于根据所有参数模型确定飞行器总体设计模型，并按照预设的优化目标和遗传算法求解所述飞行器总体设计模型，得到所述优化目标对应的飞行器总体设计参数；
其中，所述参数模型至少包括：飞行器质量模型、飞行器气动模型、弹道参数模型和冲压发动机性能模型，所述冲压发动机性能模型中的设计参数至少包括：空燃比、补燃室压力、进气道参数和燃气发生器参数；
其中，所述执行模块具体用于：
按照灵敏度分析所述飞行器总体设计模型，得到第一模型和第二模型；
当所述优化目标为射程最大时，利用遗传算法求解所述第二模型，得到第一最优解；
将所述第一最优解作为所述第一模型的初始迭代值，并利用遗传算法求解所述第一模型，得到第二最优解；
将所述第二最优解作为射程最大对应的飞行器总体设计参数；
其中，所述第一模型的计算精度高于所述第二模型的计算精度。
8.一种飞行器总体设计参数的确定设备，其特征在于，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如权利要求1-6任意一项所述的飞行器总体设计参数的确定方法的步骤。
9.一种可读存储介质，其特征在于，所述可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1-6任意一项所述的飞行器总体设计参数的确定方法的步骤。
CN201911001494.4A 2019-10-21 2019-10-21 一种飞行器总体设计参数的确定方法及装置 Active CN110516410B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 北京航空航天大学 一种高超声速飞行器跳跃滑翔弹道解析求解方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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) 一种飞行器总体设计参数的确定方法及装置
CN107575309B (zh) 一种高性能矩形双通道外并联tbcc进气道及设计方法
US20220398354A1 (en) Modeling method for integrated intake/exhaust/engine aero propulsion system with multiple geometric parameters adjustable
CN109506517B (zh) 一种带约束的中制导弹道优化方法
US9206729B2 (en) Throttleable exhaust venturi
CN112613119A (zh) 一种多几何参数可调的进/排/发一体化航空推进系统建模方法
CN105446167B (zh) 高超声速超燃冲压发动机实时模型、仿真方法
Sreenath et al. Design and analysis of contour bell nozzle and comparison with dual bell nozzle
CN115688287A (zh) 一种带引射喷管航空涡扇发动机的设计方法
CN117606295A (zh) 变掠翼制导火箭弹弹道多目标优化方法
CN110758730B (zh) 一种高超声速飞行器及其弹道设计
Hao et al. A new design method for mode transition control law of variable cycle engine
CN117932793A (zh) 采用耦合迭代优化的制导火箭弹内外弹道联合设计方法
CN118607104A (zh) 一种暴露在高空羽流中箭体的三维干扰力计算方法
Safdar et al. Numerical Analysis of Afterburner Characteristics of a Low Bypass Ratio Turbofan Engine at Various Flight Conditions
Gounko et al. Gasdynamic design of a two-dimensional supersonic inlet with the increased flow rate factor
JP2002221092A (ja) 可変ジェットノズル付き高バイパス比ターボファン・エンジン
Yuhendri et al. Impact of Cowl Deflections on Scramjet Intake Aerodynamics with Shock-on-Lip Condition
CN113742847B (zh) 用于固液动力火箭飞行器的多学科设计方法和系统
Saha et al. CFD prediction of ramjet intake characteristics at angle of attack
Bhandarkar et al. CFD Driven Aero-Propulsive Design of a Ducted Ramjet Missile
Krishna et al. Anomalies in the Flow over Projectile with Wrap-around Fins
Branch Design of Streamline-Traced Inlets and Their Integration with Rotating Detonation Engines
Kobayashi et al. Optimal design of hypersonic turbojet engines for two-stage-to-orbit spaceplane
Mizener et al. Preliminary installed performance of rotating detonation engines onto waverider configurations

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
TR01 Transfer of patent right

Effective date of registration: 20250109

Address after: D706, Zhongshan Yijing Garden, No. 3866 Nanshan Avenue, Nanshan District, Shenzhen City, Guangdong Province, 518000

Patentee after: Ming Qiang

Country or region after: China

Address before: Room 903, building B, Yongtong Jiayuan, 303 Sany Avenue, Sifangping street, Kaifu District, Changsha, Hunan 410000

Patentee before: HUNAN DEYA KUNCHUANG TECHNOLOGY Co.,Ltd.

Country or region before: China

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20250319

Address after: B105-36 Phoenix Road, Fubao Community, Fubao Street, Futian District, Shenzhen City, Guangdong Province, 518000. Ritong International Logistics Phase I and II, 3rd Floor, 328

Patentee after: Xingshan Zhilian Technology (Shenzhen) Co.,Ltd.

Country or region after: China

Address before: D706, Zhongshan Yijing Garden, No. 3866 Nanshan Avenue, Nanshan District, Shenzhen City, Guangdong Province, 518000

Patentee before: Ming Qiang

Country or region before: China

TR01 Transfer of patent right