CN111611646B - 一种飞行器前体气动构型的快速设计方法 - Google Patents
一种飞行器前体气动构型的快速设计方法 Download PDFInfo
- Publication number
- CN111611646B CN111611646B CN202010267497.9A CN202010267497A CN111611646B CN 111611646 B CN111611646 B CN 111611646B CN 202010267497 A CN202010267497 A CN 202010267497A CN 111611646 B CN111611646 B CN 111611646B
- Authority
- CN
- China
- Prior art keywords
- precursor
- design
- aircraft
- flow field
- base vector
- 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
- 239000002243 precursor Substances 0.000 title claims abstract description 89
- 238000013461 design Methods 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims abstract description 30
- 239000013598 vector Substances 0.000 claims abstract description 32
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 13
- 238000004088 simulation Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 239000012530 fluid Substances 0.000 claims description 6
- 230000035939 shock Effects 0.000 claims description 4
- 238000012938 design process Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 5
- 230000006978 adaptation Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种飞行器前体气动构型的快速设计方法。首先,根据参数化的前体设计方法设计出大量几何样本;其次,对几何样本数据进行矩阵特征值正交分解,获得其基矢量;将高维数据样本集在基矢量空间正交投影,获得基矢量的系数;最后,对基矢量系数进行插值,求得新设计变量的基矢量的系数,从而获得前体的气动型面参数;同时,对所获得前体构型样本流场仿真结果进行矩阵本征正交分解,获得基矢量,并将样本空间向基矢量空间进行投影,获得基矢量的系数;采用径向基函数对基矢量的系数进行插值重构,该方法基于矩阵本征正交分解理论建立的设计平台,能加快飞行器概念设计阶段的速度,减小计算时间和设计成本。
Description
技术领域
本发明涉及飞行器前体设计技术领域,主要涉及一种飞行器前体气动构型的快速设计方法。
背景技术
高超声速技术是一项重要的颠覆性技术,飞行器前体的设计是其关键技术之一。高超声速技术是继莱特兄弟的飞机、喷气推进之后的人类航空史上的第三次革命。军事方面,高超声速导弹可以深刻地未来战争的游戏规则;民用方面,高超声速技术若用于民机,则可将从洲际航行的民航飞机从13小时缩减为1个小时。
高超声速飞行器设计不同于常规的超声速飞行器,高超声速飞行器的前体通常采用乘波体构型,整个飞行器就“骑”在激波上,从而能最大限度地获得较高的升阻比。
然而,乘波体构型的容积效率通常较低,使得飞行器没有太多的有效容积。飞行器设计过程中,设计师通常希望在飞行器内部放置雷达天线、武器弹药,而且在可能的情况下,这些有效载荷越多越好。民用方面,若有效容积效率低,则未来高超声速客机的载客数量会减小,影响航空公司的经济收入。因此,工业界一直希望能有一种技术,使得高超声速飞行器的兼具乘波体的高升阻比,同时容积效率也较高的问题。
另一方面,升力体方案被广泛应用于导弹或者飞机的设计,飞行器的容积效率通常能够比较高。那么能否将两者融合起来,使得飞行器设计既具有乘波体的高升阻比,还具有升力体的大容积效率?基于此需求而提出本发明。
同时,过去飞行器前体的设计者属于专业人员,专业人员必须学习空气动力学和飞行器设计之后,才能根据个人经验,采用乘波体或升力体方案进行飞行器前体的设计。通常,这种设计是一种经验设计,对设计者的经验和能力要求很高。设计者具备上述能力之后,才能初步进行几何设计,然后对设计出来的前体划分网格、进行计算流体力学的数值模拟,然后获得飞行器前体的流场及其性能。如果性能不好,还得重新进行设计,再进行网格划分和流场计算、后处理等工作,如此不断循环,花费的计算代价和设计时间非常漫长。
能否自动、快速地设计出大容积率、高升阻比的飞行器前体?现代数学,尤其是线性代数理论的迅速发展给这种设计带来可能。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明采用线性代数的矩阵本征正交分解理论,结合乘波体和升力体的设计理论,提出一种快速、精确的飞行器前体设计方法。其设计时间在数秒时间量级内。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种飞行器前体气动构型的快速设计方法,包括如下步骤:
步骤1、建立基于空气动力学理论的前体设计模型;
所述前体下壁面采用密切锥设计方法,在密切面内通过特征线理论求解Taylor-Maccoll方程,获得轴对称基准流场;在该基准流场内进行流线追踪,获得前体的下壁面几何构型;
所述前体上壁面构型采用冯卡门曲线设计方法;
所述前体上壁面出口型线采用类型函数设计方法;
设计过程中,通过调整前体半宽度D、类型函数指数n、来流条件、上壁面基准高度H、上表面出口型线方程、前体前缘型线方程和前体激波出口型线方程,可得到不同气动构型的飞行器前体,获得所需几何前体样本集合;对所述几何前体样本集合分别开展计算流体力学数值模拟,获得几何前体样本构型的流场样本集合;
步骤2、对步骤1中所述几何前体样本集合{x(M),y(M),zi(M)}T,进行本征正交分解,获得其特征值λ1>λ2>…λP及其对应的特征向量;其中M为样本数目;根据特征值大小所占比重,选取总和占全部特征值99%的特征值,用于获取设计参数空间内相互正交的基矢量,如下所示:
步骤3、根据步骤2中选取的特征值,确定投影空间的基矢量系数,在投影空间内,前体几何型面可以采用特征值所对应的特征向量进行线性加权,如下所示:
式中,Φi为特征向量,通过矩阵分解获得;ai为系数;求解ai时约束为所有样本的型面函数S(x,y,z)与预测结果的误差最小;即满足如下约束条件:
步骤4、采用伽辽金投影方法获得系数ai,采用径向基函数来对新设计变量进行插值,获得新的系数,快速确定所需设计的前体气动几何型面;
步骤5、参照步骤2-步骤5、对各个前体型面样本的流场结果进行本征正交分解,提取其99%的特征,进行伽辽金投影,获得简化的设计参数空间的基矢量及其系数;最后,针对新的设计变量,采用径向基函数进行系数计算,从而快速获得所需设计点前体型面对应的流场。
进一步地,所述步骤1中在调整各项设计参数时,采用优化拉丁超立方试验设计,使得设计参数均匀分布,获得设计参数的数据集合;根据设计参数的数据集合,采用步骤1的方法获得一系列不同气动构型的几何前体样本集合;并且通过对几何前体样本集合进行计算流体力学模拟,获得这些前体样本集合的流场样本集合。
有益效果:
1、经过前期样本数据处理,构建好几何快速设计平台和流场与性能快速预估平台之后,进行具体的飞行器前体设计时能够快速得到所需前体几何构型和流场及性能数据,极大缩短设计和计算时间,减小计算设计成本。
2、降低了设计难度和对设计者专业要求,设计者仅需根据设计需求,改变不同设计变量即可获得前体构型及其流场与性能数据。
3、本发明对于未来飞行器精细化、自动化和快速化设计具有良好的应用潜能。
附图说明
图1为本发明提供的前体气动型面快速设计流程图;
图2为本发明提供的前体流场快速预测流程图;
图3为飞行器前体设计示意图;
图4为前体密切面设计示意图;
图5为本发明提供的快速设计前体几何构型图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1-2所示的一种飞行器前体气动构型的快速设计方法,包括如下步骤:
步骤1、建立基于空气动力学理论的前体设计模型;
前体下壁面采用密切锥设计方法,在密切面内通过特征线理论求解Taylor-Maccoll方程,获得轴对称基准流场;在该基准流场内进行流线追踪,获得前体的下壁面几何构型;
前体上壁面构型采用冯卡门曲线设计方法;
前体上壁面出口型线采用类型函数设计方法;
设计过程中,通过调整前体半宽度D、类型函数指数n、来流条件、上壁面基准高度H、上表面出口型线(USEC,Upper Surface Exit Curve)方程、前体前缘型线(FCT,ForebodyCapture Tube)方程和前体激波出口型线(ICC,Inlet Capture Curve)方程,可得到不同气动构型的飞行器前体,获得所需几何前体样本集合;对所述几何前体样本集合分别开展计算流体力学数值模拟,获得几何前体样本构型的流场样本集合。
调整各项设计参数时,采用优化拉丁超立方试验设计,使得设计参数均匀分布,获得设计参数的数据集合;根据设计参数的数据集合,采用步骤1的方法获得一系列不同气动构型的几何前体样本集合;并且通过对几何前体样本集合进行计算流体力学模拟,获得这些前体样本集合的流场样本集合。
步骤2、对步骤1中所述几何前体样本集合{x(M),y(M),zi(M)}T,进行本征正交分解,获得其特征值λ1>λ2>…λP及其对应的特征向量;其中M为样本数目;根据特征值大小所占比重,选取总和占全部特征值99%的特征值,用于获取设计参数空间内相互正交的基矢量,如下所示:
步骤3、根据步骤2中选取的特征值,确定投影空间的基矢量系数,在投影空间内,前体几何型面可以采用特征值所对应的特征向量进行线性加权,如下所示:
式中,Φi为特征向量,通过矩阵分解获得;ai为系数;求解ai时约束为所有样本的型面函数S(x,y,z)与预测结果的误差最小;即满足如下约束条件:
步骤4、采用伽辽金投影方法获得系数ai,采用径向基函数来对新设计变量进行插值,获得新的系数,快速确定所需设计的前体气动几何型面。
步骤5、参照步骤2-步骤5、对各个前体型面样本的流场结果进行本征正交分解,提取其99%的特征,进行伽辽金投影,获得简化的设计参数空间的基矢量及其系数;最后,针对新的设计变量,采用径向基函数进行系数计算,从而快速获得所需设计点前体型面对应的流场。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (2)
1.一种飞行器前体气动构型的快速设计方法,其特征在于,包括如下步骤:
步骤1、建立基于空气动力学理论的前体设计模型;
前体下壁面采用密切锥设计方法,在密切面内通过特征线理论求解Taylor-Maccoll方程,获得轴对称基准流场;在基准流场内进行流线追踪,获得前体下壁面构型;
前体上壁面构型采用冯卡门曲线设计方法;
前体上壁面出口型线采用类型函数设计方法;
设计过程中,通过调整前体半宽度D、类型函数指数n、来流条件、上壁面基准高度H、上表面出口型线方程、前体前缘型线方程和前体激波出口型线方程,得到不同气动构型的飞行器前体,从而获得所需的前体几何结构样本集合;对所述前体几何结构样本集合分别开展计算流体力学数值模拟,获得不同前体构型对应的流场样本集合;
步骤2、对步骤1中的前体几何结构样本集合进行本征值正交分解,获得对应特征值λ1>λ2>…λP及对应的特征向量;根据特征值大小所占比重,选取总和占全部特征值之和99%的特征值,用于获取设计参数空间内相互正交的基矢量,如下所示:
步骤3、根据步骤2中选取的特征值,确定投影空间的基矢量系数,在投影空间内,前体几何型面采用特征值所对应的特征向量进行线性加权,如下所示:
式中,Φi为特征向量,通过矩阵分解获得;ai为系数;求解ai时约束为所有样本的型面函数S(x,y,z)与预测结果的误差最小;即满足如下约束条件:
步骤4、对于新的前体几何结构,采用径向基函数来对新设计变量进行插值,获得新的基矢量系数,从而快速确定所需设计的前体气动几何型面;
步骤5、参照步骤2-步骤5,对各个前体几何结构样本的流场结果进行本征值正交分解,提取99%的特征,进行伽辽金投影,获得简化的设计参数空间的基矢量及系数;最后,针对新的设计变量,采用径向基函数进行系数计算,从而快速获得所需设计的前体几何结构对应的流场。
2.根据权利要求1所述的一种飞行器前体气动构型的快速设计方法,其特征在于,所述步骤1中在选取各项设计参数时,采用优化拉丁超立方进行采样设计,使得设计参数均匀分布,获得设计参数的数据集合;根据设计参数的数据集合,采用步骤1的方法获得一系列不同气动构型的前体几何结构样本集合;并且通过对前体几何结构样本集合进行流体力学模拟,获得前体几何结构样本集合的流场样本集合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010267497.9A CN111611646B (zh) | 2020-04-08 | 2020-04-08 | 一种飞行器前体气动构型的快速设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010267497.9A CN111611646B (zh) | 2020-04-08 | 2020-04-08 | 一种飞行器前体气动构型的快速设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111611646A CN111611646A (zh) | 2020-09-01 |
CN111611646B true CN111611646B (zh) | 2024-04-12 |
Family
ID=72201697
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010267497.9A Active CN111611646B (zh) | 2020-04-08 | 2020-04-08 | 一种飞行器前体气动构型的快速设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111611646B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111240209B (zh) | 2020-03-16 | 2020-10-09 | 广东工业大学 | 构型动型控型优型联动响应的自适应组态方法及系统 |
CN113127994A (zh) * | 2021-04-29 | 2021-07-16 | 复旦大学 | 考虑三维因素的短舱进气道设计几何参数化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103902782A (zh) * | 2014-04-11 | 2014-07-02 | 北京理工大学 | 基于pod和代理模型的高超声速气动热模型降阶方法 |
CN105975645A (zh) * | 2016-02-26 | 2016-09-28 | 西北工业大学 | 一种基于多步的含激波区域飞行器流场快速计算方法 |
CN107220403A (zh) * | 2017-04-20 | 2017-09-29 | 南京航空航天大学 | 飞行器弹性模态的控制关联建模方法 |
CN110059410A (zh) * | 2019-04-19 | 2019-07-26 | 中国人民解放军国防科技大学 | 一种基于b样条的腹部进气高超声速前体优化设计方法 |
-
2020
- 2020-04-08 CN CN202010267497.9A patent/CN111611646B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103902782A (zh) * | 2014-04-11 | 2014-07-02 | 北京理工大学 | 基于pod和代理模型的高超声速气动热模型降阶方法 |
CN105975645A (zh) * | 2016-02-26 | 2016-09-28 | 西北工业大学 | 一种基于多步的含激波区域飞行器流场快速计算方法 |
CN107220403A (zh) * | 2017-04-20 | 2017-09-29 | 南京航空航天大学 | 飞行器弹性模态的控制关联建模方法 |
CN110059410A (zh) * | 2019-04-19 | 2019-07-26 | 中国人民解放军国防科技大学 | 一种基于b样条的腹部进气高超声速前体优化设计方法 |
Non-Patent Citations (2)
Title |
---|
基于本征正交分解和代理模型的高超声速气动热模型降阶研究;陈鑫等;《航空学报》;第462-472页 * |
基于本征正交分解的气动外形设计空间重构方法研究;刘南等;《西北工业大学学报》;第171-177页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111611646A (zh) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111079228B (zh) | 一种基于流场预测的气动外形优化方法 | |
CN111611646B (zh) | 一种飞行器前体气动构型的快速设计方法 | |
Mengistu et al. | Aerodynamic optimization of turbomachinery blades using evolutionary methods and ANN-based surrogate models | |
Fike et al. | Optimization with gradient and hessian information calculated using hyper-dual numbers | |
US6606612B1 (en) | Method for constructing composite response surfaces by combining neural networks with other interpolation or estimation techniques | |
CN107515994B (zh) | 一种自主水下航行器的壳体结构多保真度设计优化方法 | |
CN113297686B (zh) | 固体火箭发动机数据融合设计方法、装置、设备及介质 | |
CN104750948B (zh) | 一种处理飞行器设计中多极值多约束问题的优化方法 | |
CN113673027B (zh) | 一种基于代理模型的高超声速飞行器气动载荷优化设计方法 | |
Dennis et al. | Multi-objective optimization of turbomachinery cascades for minimum loss, maximum loading, and maximum gap-to-chord ratio | |
CN111898212A (zh) | 基于BezierGAN和贝叶斯优化的叶轮机械型线设计优化方法 | |
Heidebrecht et al. | Parametric geometry and CFD process for turbofan nacelles | |
CN108446481B (zh) | 一种基于灰狼群智能优化的低阻低噪便携式auv外形设计方法 | |
Liou et al. | Challenges and progress in aerodynamic design of hybrid wingbody aircraft with embedded engines | |
CN114662209A (zh) | 一种碟形水下直升机外形的cfd仿真优化及验证方法 | |
Mengistu et al. | Aerodynamic shape optimization of turbine blades using a design-parameter-based shape representation | |
CN108664742B (zh) | 短舱参数化设计的多目标优化方法 | |
Zheng et al. | Dynamic space reduction optimization framework and its application in hull form optimization | |
Sasaki et al. | Navier-Stokes optimization of supersonic wings with four design objectives using evolutionary algorithm | |
CN115577453A (zh) | 一种变距螺旋桨气动设计方法、装置及系统 | |
Yi et al. | A fast forecast method based on high and low fidelity surrogate models for strength and stability of stiffened cylindrical shell with variable ribs | |
CN115358167A (zh) | 一种考虑发动机参数的飞发一体气动伴随优化设计方法 | |
Popov et al. | Estimation of the Effect of 3D Grid Parameters on the Simulation of the Working Process of Axial Turbines. | |
CN109214056B (zh) | 一种基于流动物理的气动优化设计变量选取方法 | |
CN113283021A (zh) | 一种确定后涵道引射器可调面积的方法 |
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 |