CN115238396B - 三维气动外形反设计方法及装置 - Google Patents
三维气动外形反设计方法及装置 Download PDFInfo
- Publication number
- CN115238396B CN115238396B CN202211093692.XA CN202211093692A CN115238396B CN 115238396 B CN115238396 B CN 115238396B CN 202211093692 A CN202211093692 A CN 202211093692A CN 115238396 B CN115238396 B CN 115238396B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- target
- deformation control
- pressure distribution
- control point
- 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
- 238000013461 design Methods 0.000 title claims abstract description 169
- 238000000034 method Methods 0.000 title claims abstract description 139
- 238000009826 distribution Methods 0.000 claims abstract description 138
- 230000009467 reduction Effects 0.000 claims abstract description 63
- 238000004422 calculation algorithm Methods 0.000 claims description 23
- 238000012549 training Methods 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 239000012530 fluid Substances 0.000 claims description 6
- 230000001502 supplementing effect Effects 0.000 claims description 5
- 239000013598 vector Substances 0.000 claims description 5
- 238000004590 computer program Methods 0.000 description 9
- 238000013507 mapping Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000005457 optimization Methods 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 4
- 238000012938 design process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000010801 machine learning Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013079 data visualisation Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
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
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Processing Or Creating Images (AREA)
Abstract
本申请公开了一种三维气动外形反设计方法及装置,通过获取每个样本点的参数化变形控制点和气动外形;根据样本点的参数化变形控制点和气动外形,建立高维数据;采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;获取目标参数化变形控制点,对目标参数化变形控制点进行压力系数插值,获得与目标参数化变形控制点对应的目标压力分布值;根据目标压力分布值,通过POD降维模型,确定目标压力分布下的设计变量值;根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,利用本申请实施例可以直接对三维气动外形进行反设计。
Description
技术领域
本申请属于流体力学技术领域,尤其涉及一种三维气动外形反设计方法及装置。
背景技术
气动外形设计是飞行器设计中极为重要的组成部分,其在飞行器设计的整个阶段有着主导作用,气动外形直接决定了整个飞行器的飞行性能与飞行品质。气动外形和气动性能之间存在复杂的非线性关系且受多种约束限制,使得气动外形设计成为整个飞行器设计流程中最为困难的环节之一,往往需要贯穿整个飞行器设计过程。现代飞行器的高载荷、高效率、高性能、高可靠性和经济性等发展趋势,使得传统的飞行器设计方法越来越难以满足设计需求,随着计算机技术的高速发展,基于计算流体力学(CFD)和数值优化方法的飞行器气动外形数值设计技术取得了跨越式的发展。这些数值设计技术主要分为两类:(1)将空气动力学正问题与优化算法结合的气动外形优化设计方法;(2)根据空气动力学反问题建立的气动外形设计方法。
反设计方法是在设计人员丰富的气动设计经验基础上,根据设计要求针对气动外形提出合理的压力分布或速度分布,然后通过改变气动外形,使气动外形表面的压力分布或速度分布吻合目标压力分布或速度分布。该方法设计效率高,近些年得到了研究人员的广泛关注,并在翼型和机翼的气动外形优化设计中取得了相当大的成果。现行的反设计方法,一般需要设计人员具有相当的设计经验,且目标压力分布较为真实合理,否则难以收敛,不能达到预期设计效果。
随着机器学习技术不断发展,其数据降维思想和方法已被广泛引入到与飞行器设计相关的各个领域,来提升高维数据下的计算效率,或高维数据可视化过程。目前,基于机器学习降维方法和神经网络的反设计方法,多以二维翼型为研究对象,难以直接对三维飞行器进行反设计。
发明内容
本发明意在提供一种三维气动外形反设计方法及装置,以解决现有技术中存在的不足,本发明要解决的技术问题通过以下技术方案来实现。
第一个方面,本发明实施例提供一种三维气动外形反设计方法,所述方法包括:
获取每个样本点的参数化变形控制点和气动外形;
根据所述样本点的参数化变形控制点和所述气动外形,建立高维数据;
采用ISOMAP非线性降维方法,将所述高维数据降到低维流形空间的低维数据;
在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;
根据所述目标压力分布值,通过POD降维模型,确定所述目标压力值分布下的设计变量值;
根据所述目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形。
可选地,所述采用ISOMAP非线性降维方法,将所述高维数据降到低维流形空间的低维数据,包括:
获取N个D维数据的样本;
根据所述样本中任意两点,确定领域图;
在所述领域图中,根据Dijkstra算法或Floyd算法,计算所述领域图中所有点对的最短距离,并将得到的最短距离的线段补充在所述领域图中,得到重构的领域图;
采用MDS算法对所述重构的领域图进行降维处理,得到低维流形空间的低维数据。
可选地,所述根据所述样本中任意两点,确定领域图,包括:
计算样本中任意两点i和j的距离;
判断i和j是否相邻,若i和j相邻,则将i和j两点连接起来,计算两点之间的长度;
根据满足条件的所有点,确定领域图。
可选地,所述采用MDS算法对所述重构的领域图进行降维处理,得到低维流形空间的低维数据,包括:
根据所述重构领域图,确定内积矩阵;
计算所述内积矩阵的特征值和特征向量;
根据所述特征向量,确定原始数据集在低维流形上的坐标。
可选地,所述在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值,包括:
在低维流形空间中,获取目标参数化变形控制点;
根据气动外形的流形结构和所述目标参数化变形控制点的流形空间坐标值进行RBF插值,获得所述目标参数化变形控制点的目标压力分布值,其中,所述气动外形的流形结构为气动外形在低维空间中的坐标值。
可选地,所述根据所述目标压力分布值,通过POD降维模型,确定所述目标压力值分布下的设计变量值,包括:
将所述目标参数化变形控制点的目标压力分布值和对应的设计变量进行组合,获得训练数据集;
将所述训练数据集与所述目标参数化变形控制点的目标压力分布值进行组合,确定缺失目标设计变量的数据集;
通过POD降维模型迭代出符合参数化变形控制点对应的目标压力分布的设计变量值。
可选地,所述根据所述目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,包括:
根据所述目标压力值分布下的设计变量值,利用FFD参数化方法对原始气动外形进行参数化处理,得到目标气动外形;
在通过网格变形方法获得目标气动外形对应的计算网格;
利用CFD方法,确定所述目标气动外形下的压力分布值,并根据所述压力分布值,确定与所述设计变量值对应的压力分布;
若所述压力分布满足预设条件,则得到相应的三维气动外形。
第二个方面,本发明实施例提供一种三维气动外形反设计装置,所述装置包括:
获取模块,用于获取每个样本点的参数化变形控制点和气动外形;
建立模块,用于根据所述样本点的参数化变形控制点和所述气动外形,建立高维数据;
降维模块,用于采用ISOMAP非线性降维方法,将所述高维数据降到低维流形空间的低维数据;
插值模块,用于在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;
确定模块,用于根据所述目标压力分布值,通过POD降维模型,确定所述目标压力值分布下的设计变量值;
设计模块,用于根据所述目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形。
可选地,所述降维模块用于:
获取N个D维数据的样本;
根据所述样本中任意两点,确定领域图;
在所述领域图中,根据Dijkstra算法或Floyd算法,计算所述领域图中所有点对的最短距离,并将得到的最短距离的线段补充在所述领域图中,得到重构的领域图;
采用MDS算法对所述重构的领域图进行降维处理,得到低维流形空间的低维数据。
可选地,所述降维模块用于:
计算样本中任意两点i和j的距离;
判断i和j是否相邻,若i和j相邻,则将i和j两点连接起来,计算两点之间的长度;
根据满足条件的所有点,确定领域图。
可选地,所述降维模块用于:
根据所述重构领域图,确定内积矩阵;
计算所述内积矩阵的特征值和特征向量;
根据所述特征向量,确定原始数据集在低维流形上的坐标。
可选地,所述插值模块用于:
在低维流形空间中,获取目标参数化变形控制点;
根据气动外形的流形结构和所述目标参数化变形控制点的流形空间坐标值进行RBF插值,获得所述目标参数化变形控制点的目标压力分布值,其中,所述气动外形的流形结构为气动外形在低维空间中的坐标值。
可选地,所述确定模块用于:
将所述目标参数化变形控制点的目标压力分布值和对应的设计变量进行组合,获得训练数据集;
将所述训练数据集与所述目标参数化变形控制点的目标压力分布值进行组合,确定缺失目标设计变量的数据集;
通过POD降维模型迭代出符合参数化变形控制点对应的目标压力分布的设计变量值。
可选地,所述设计模块用于:
根据所述目标压力值分布下的设计变量值,利用FFD参数化方法对原始气动外形进行参数化处理,得到目标气动外形;
在通过网格变形方法获得目标气动外形对应的计算网格;
利用CFD方法,确定所述目标气动外形下的压力分布值,并根据所述压力分布值,确定与所述设计变量值对应的压力分布;
若所述压力分布满足预设条件,则得到相应的三维气动外形。
第三个方面,本发明实施例提供一种终端设备,包括:至少一个处理器和存储器;
所述存储器存储计算机程序;所述至少一个处理器执行所述存储器存储的计算机程序,以实现第一个方面提供的三维气动外形反设计方法。
第四个方面,本发明实施例提供一种计算机可读存储介质,该计算机可读存储介质中存储有计算机程序,所述计算机程序被执行时实现第一个方面提供的三维气动外形反设计方法。
本发明实施例包括以下优点:
本发明实施例提供的三维气动外形反设计方法及装置,通过获取每个样本点的参数化变形控制点和气动外形;根据样本点的参数化变形控制点和气动外形,建立高维数据;采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,利用 ISOMAP方法的特点,将三维气动外形降维到二维流形空间,并结合POD模型,提出了基于等距特征映射的三维气动外形反设计方法,直接对三维气动外形进行反设计。
附图说明
为了更清楚地说明本申请实施例或现有的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请一实施例中一种三维气动外形反设计方法的流程示意图;
图2为本申请一实施例中基于流形学习的三维机翼反设计流程;
图3为本申请一实施例中M6机翼及设计点;
图4为本申请一实施例中机翼上翼面低维流形结构和设计点分布图;
图5为本申请一实施例中低维流形空间的设计点压力插值;
图6为本申请一实施例中M6机翼反设计结果;
图7是本发明的一种三维气动外形反设计装置实施例的结构框图;
图8是本发明的一种终端设备的结构示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合具体实施例及相应的附图对本申请的技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
等距特征映射(ISOMAP)作为机器学习中流形学习领域的降维模型,能实现高维数据的非线性降维,被广泛应用于数据分析及可视化分析。ISOMAP通过测地距离的方式将高维空间的拓扑关系映射到流形空间,保证了数据集中任意两点的距离在降维前后不变,从而捕捉数据集所在低维流形的全局结构。
本发明一实施例提供一种三维气动外形反设计方法,用于对三维机翼进行反设计。本实施例的执行主体为三维气动外形反设计装置,设置在计算机上。
参照图1,示出了本发明的一种三维气动外形反设计方法实施例的步骤流程图,该方法具体可以包括如下步骤:
S101、获取每个样本点的参数化变形控制点和气动外形;
具体地,利用拉丁超立方方法获得样本点;
S102、根据样本点的参数化变形控制点和气动外形,建立高维数据;
将每个样本下的气动外形和参数化变形控制点组合,构成高维非线性数据集,其中每个样本下气动外形的所有网格节点坐标,是每个样本下参数化变形控制点的三维坐标。
S103、采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;
根据ISOMAP非线性降维方法获得每个样本高维数据集的低维流形空间。
S104、在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;
在低维流形空间中根据气动外形的流形结构(即气动外形在流形空间中的坐标值)和参数化变形控制点的流形空间坐标值进行RBF插值,获得每个参数化变形控制点的压力分布值。
S105、根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;
S106、根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形。
本申请实施例通过ISOMAP非线性降维方法,将每个样本的参数化变形控制点三维坐标和气动外形三维坐标组成高维数据,并进行降维,获得其在低维流形空间中的二维结构,在低维流形空间对参数化变形控制点进行压力系数插值,获得参数化变形控制点对应的压力分布。将插值获得的参数化变形控制点压力分布与设计变量组合,构成样本数据集,并通过POD降维模型捕捉目标压力分布下的设计变量值。根据获得的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,完成三维机翼的反设计。
本发明实施例提供的三维气动外形反设计方法,通过获取每个样本点的参数化变形控制点和气动外形;根据样本点的参数化变形控制点和气动外形,建立高维数据;采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,利用 ISOMAP方法的特点,将三维气动外形降维到二维流形空间,并结合POD模型,提出了基于等距特征映射的三维气动外形反设计方法,直接对三维气动外形进行反设计。
本发明又一实施例对上述实施例提供的三维气动外形反设计方法做进一步补充说明。
可选地,采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据,包括:
获取N个D维数据的样本;
根据样本中任意两点,确定领域图;
在领域图中,根据Dijkstra算法或Floyd算法,计算领域图中所有点对的最短距离,并将得到的最短距离的线段补充在领域图中,得到重构的领域图;
采用MDS算法对重构的领域图进行降维处理,得到低维流形空间的低维数据。
可选地,根据样本中任意两点,确定领域图,包括:
计算样本中任意两点i和j的距离;
判断i和j是否相邻,若i和j相邻,则将i和j两点连接起来,计算两点之间的长度;
根据满足条件的所有点,确定领域图。
可选地,采用MDS算法对重构的领域图进行降维处理,得到低维流形空间的低维数据,包括:
根据重构领域图,确定内积矩阵;
计算内积矩阵的特征值和特征向量;
根据特征向量,确定原始数据集在低维流形上的坐标。
可选地,在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值,包括:
在低维流形空间中,根据气动外形的流形结构和参数化变形控制点的流形空间坐标值进行RBF插值,获得每个参数化变形控制点的压力分布值,其中,气动外形的流形结构为气动外形在低维空间中的坐标值。
可选地,根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值,包括:
将参数化变形控制点的压力分布值和对应的设计变量进行组合,获得训练数据集;
将训练数据集与参数化变形控制点的压力分布值进行组合,确定缺失目标设计变量的数据集;
通过POD降维模型迭代出符合参数化变形控制点对应的目标压力分布的设计变量值。
可选地,根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,包括:
根据设计变量值,利用FFD参数化方法对原始气动外形进行参数化处理,得到目标气动外形;
在通过网格变形方法获得目标气动外形对应的计算网格;
利用CFD方法,确定目标气动外形下的压力分布值,并根据压力分布值,确定与设计变量值对应的压力分布;
若压力分布满足预设条件,则得到相应的三维气动外形。
图2为本申请一实施例中基于流形学习的三维机翼反设计流程;以跨声速M6机翼(图3所示),在低维流形空间的反设计为具体实例进行说明。
该流程包括:
1)气动外形和设计点组成的高维数据集;
2)ISOMAP降维并差值参数化变形控制点对应的压力分布;
3)获得设计变量和压力分布组成的数据集;
4)降维迭代;
5)最优设计结果;
6)进行参数化处理;
7)进行网格变形;
8)CFD校验;
9)对校验的结果进行判断,是否满足设计要求;
10)若不满足设计要求,则设计压力分布,且更新训练数据集。
本发明解决其技术问题所采用的技术方案包括以下步骤:
① 利用拉丁超立方方法获得样本点;
② 将每个样本下的气动外形Si(xi, yi, zi)和参数化变形控制点Xi(xi, yi, zi)组合,构成高维非线性数据集Φi=[Si, Xi]T,其中Si每个样本下气动外形的所有网格节点坐标,Xi是每个样本下参数化变形控制点的三维坐标。
③ 根据ISOMAP非线性降维方法获得每个样本高维数据集Φi的低维流形空间Ψi,其中ISOMAP的具体降维流程如下:
b、根据Dijkstra算法或Floyd算法计算两点之间的最短距离。即在领域图G中,计算所有点对之间的最短路径,以重新填充步骤(1)中所构建的领域图G。具体为查看所有点对i—j,是否存在一条更短的路径i—k—j进行替代,即:若,则,否则
c、利用重构的领域图G结合MDS算法进行降维,即将领域图G代入式(1)求得内积矩阵B,然后计算内积矩阵B的特征值与特征向量,根据特征向量得到原始数据集在低维流形上的坐标。
本发明的有益效果是:(1)通过流形学习的非线性算法,构造三维气动外形的低维流形结构,并将设计目标转换到参数化变形控制点的压力分布上,使得三维气动外形反设计方法更加工程实用。(2)将三维气动外形的压力分布数据降维与设计变量相同的维度,利用POD方法在正交空间中进行反设计,进一步降低了设计维度,极大的提高了设计效率。(3)在非线性降维方法构造的低维流形空间中插值参数化变形控制点的压力分布,确保了设计点与气动外形的空间非线性,且参数化变形控制点的压力分布插值结果具有较高精度,保证了方法的可靠性。(4)无需严格定义设计目标的物理意义,使得反设计方法更加灵活性和鲁棒性。
对M6机翼进行反设计,计算状态为马赫数Ma=0.83950,雷诺数Re=1.172×107,攻角α=3.06,湍流模型采用SA模型,温度T=288.15。采用FFD参数化方法对M6机翼气动外形进行参数化,设计变量为FFD控制体上翼面110个控制点的z方向位移,如图3所示。具体实施步骤为:
步骤1:采用拉丁超立方抽样方法获得200个样本点,通过FFD参数化方法对M6机翼进行参数化,通过CFD求解器计算得到M6机翼表面的压力分布和气动响应。
步骤2:将每一个样本点下的M6机翼上翼面的三维气动外形Si(xi, yi, zi)和110个参数化变形控制点 Xi(xi, yi, zi)组合,构建每个样本点的高维非线性数据集Φi=[Si,Xi]T;
步骤3:采用ISOMAP非线性降维方法对步骤2中所得的每个样本点的高维非线性数据集 Φi,进行非线性降维,获得M6机翼上翼面二维流形结构和110个参数化变形控制点在流形空间中的位置,如图4所示,图4为本申请一实施例中机翼上翼面低维流形结构和设计点分布图,图中星号为二维流形空间中110个参数化变形控制点的分布,圆点为M6机翼上翼面在流形空间中的流形结构和压力分布情况;在图4中的横纵坐标没有具体的物理意义,仅仅是显示的刻度值,图中的星号表示每个点插值的压力系数大小。
步骤4:根据步骤3中所得的上翼面流形结构(即上翼面网格节点在流形空间中的坐标和其压力分布),及流形空间中分布的参数化变形控制点,采用RBF插值方法,计算出参数化变形控制点的压力分布,如图5所示,图5为本申请一实施例中低维流形空间的设计点压力插值;其中,横纵坐标没有具体的物理意义,仅仅是显示的刻度值,图中的星号表示每个点插值的压力系数大小,而且图4中的Isomap_1和图5中的Isomap_2为一种图的标识,并不具体实际物理意义。
步骤5:将每个样本的设计变量和步骤4中所得的参数化变形控制点的压力分布进行组合,得到数据集Q;
步骤6:将目标压力分布与步骤5中所得数据集Q得到缺失目标设计变量的数据集Y;
步骤8:根据步骤7所得的设计变量D T ,利用FFD参数化方法对初始气动外形Sorg进行参数化,得到目标气动外形S T ;在通过网格变形方法获得新气动外形对应的计算网格。
图6给出了M6机翼的反设计结果。从图6中可知,本申请所提供的基于等距特征映射方法的三维气动外形反设计方法,能有效完成三维气动外形的反设计,且将设计目标定义在参数化变形控制点的压力分布上,使得方法工程实用的同时具有较高的设计效率。此外,设计目标定义灵活,不需要严格定义三维气动外形压力分布的物理意义,极大地提升了反设计方法的鲁棒性,使得三维气动外形的设计变的更加灵活,使用更加方便。
需要说明的是,对于方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明实施例并不受所描述的动作顺序的限制,因为依据本发明实施例,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作并不一定是本发明实施例所必须的。
本发明实施例提供的三维气动外形反设计方法,通过获取每个样本点的参数化变形控制点和气动外形;根据样本点的参数化变形控制点和气动外形,建立高维数据;采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,利用 ISOMAP方法的特点,将三维气动外形降维到二维流形空间,并结合POD模型,提出了基于等距特征映射的三维气动外形反设计方法,直接对三维气动外形进行反设计。
本发明另一实施例提供一种三维气动外形反设计装置,用于执行上述实施例提供的三维气动外形反设计方法。
参照图7,示出了本发明的一种三维气动外形反设计装置实施例的结构框图,该装置具体可以包括如下模块:获取模块701、建立模块702、降维模块703、插值模块704、确定模块705和设计模块706,其中:
获取模块701用于获取每个样本点的参数化变形控制点和气动外形;
建立模块702用于根据样本点的参数化变形控制点和气动外形,建立高维数据;
降维模块703用于采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;
插值模块704用于在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;
确定模块705用于根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;
设计模块706用于根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形。
本发明实施例提供的三维气动外形反设计装置,通过获取每个样本点的参数化变形控制点和气动外形;根据样本点的参数化变形控制点和气动外形,建立高维数据;采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,利用 ISOMAP方法的特点,将三维气动外形降维到二维流形空间,并结合POD模型,提出了基于等距特征映射的三维气动外形反设计方法,直接对三维气动外形进行反设计。
本发明又一实施例对上述实施例提供的三维气动外形反设计装置做进一步补充说明。
可选地,降维模块用于:
获取N个D维数据的样本;
根据样本中任意两点,确定领域图;
在领域图中,根据Dijkstra算法或Floyd算法,计算领域图中所有点对的最短距离,并将得到的最短距离的线段补充在领域图中,得到重构的领域图;
采用MDS算法对重构的领域图进行降维处理,得到低维流形空间的低维数据。
可选地,降维模块用于:
计算样本中任意两点i和j的距离;
判断i和j是否相邻,若i和j相邻,则将i和j两点连接起来,计算两点之间的长度;
根据满足条件的所有点,确定领域图。
可选地,降维模块用于:
根据重构领域图,确定内积矩阵;
计算内积矩阵的特征值和特征向量;
根据特征向量,确定原始数据集在低维流形上的坐标。
可选地,插值模块用于:
在低维流形空间中,根据气动外形的流形结构和参数化变形控制点的流形空间坐标值进行RBF插值,获得每个参数化变形控制点的压力分布值,其中,气动外形的流形结构为气动外形在低维空间中的坐标值。
可选地,确定模块用于:
将参数化变形控制点的压力分布值和对应的设计变量进行组合,获得训练数据集;
将训练数据集与参数化变形控制点的压力分布值进行组合,确定缺失目标设计变量的数据集;
通过POD降维模型迭代出符合参数化变形控制点对应的目标压力分布的设计变量值。
可选地,设计模块用于:
根据设计变量值,利用FFD参数化方法对原始气动外形进行参数化处理,得到目标气动外形;
在通过网格变形方法获得目标气动外形对应的计算网格;
利用CFD方法,确定目标气动外形下的压力分布值,并根据压力分布值,确定与设计变量值对应的压力分布;
若压力分布满足预设条件,则得到相应的三维气动外形。
对于装置实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本发明实施例提供的三维气动外形反设计装置,通过获取每个样本点的参数化变形控制点和气动外形;根据样本点的参数化变形控制点和气动外形,建立高维数据;采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,利用 ISOMAP方法的特点,将三维气动外形降维到二维流形空间,并结合POD模型,提出了基于等距特征映射的三维气动外形反设计方法,直接对三维气动外形进行反设计。
本发明再一实施例提供一种终端设备,用于执行上述实施例提供的三维气动外形反设计方法。
图8是本发明的一种终端设备的结构示意图,如图8所示,该终端设备包括:至少一个处理器801和存储器802;
存储器存储计算机程序;至少一个处理器执行存储器存储的计算机程序,以实现上述实施例提供的三维气动外形反设计方法。
本实施例提供的终端设备,通过获取每个样本点的参数化变形控制点和气动外形;根据样本点的参数化变形控制点和气动外形,建立高维数据;采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,利用 ISOMAP方法的特点,将三维气动外形降维到二维流形空间,并结合POD模型,提出了基于等距特征映射的三维气动外形反设计方法,直接对三维气动外形进行反设计。
本申请又一实施例提供一种计算机可读存储介质,该计算机可读存储介质中存储有计算机程序,计算机程序被执行时实现上述任一实施例提供的三维气动外形反设计方法。
根据本实施例的计算机可读存储介质,通过获取每个样本点的参数化变形控制点和气动外形;根据样本点的参数化变形控制点和气动外形,建立高维数据;采用ISOMAP非线性降维方法,将高维数据降到低维流形空间的低维数据;在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;根据目标压力分布值,通过POD降维模型,确定目标压力值分布下的设计变量值;根据目标压力值分布下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,利用 ISOMAP方法的特点,将三维气动外形降维到二维流形空间,并结合POD模型,提出了基于等距特征映射的三维气动外形反设计方法,直接对三维气动外形进行反设计。
应该指出,上述详细说明都是示例性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语均具有与本申请所属技术领域的普通技术人员的通常理解所相同的含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式。此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的术语在适当情况下可以互换,以便这里描述的本申请的实施方式能够以除了在这里图示或描述的那些以外的顺序实施。
此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含。例如,包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
为了便于描述,在这里可以使用空间相对术语,如“在……之上”、“在……上方”、“在……上表面”、“上面的”等,用来描述如在图中所示的一个器件或特征与其他器件或特征的空间位置关系。应当理解的是,空间相对术语旨在包含除了器件在图中所描述的方位之外的在使用或操作中的不同方位。例如,如果附图中的器件被倒置,则描述为“在其他器件或构造上方”或“在其他器件或构造之上”的器件之后将被定位为“在其他器件或构造下方”或“在其他器件或构造之下”。因而,示例性术语“在……上方”可以包括“在……上方”和“在……下方”两种方位。该器件也可以其他不同方式定位,如旋转90度或处于其他方位,并且对这里所使用的空间相对描述作出相应解释。
在上面详细的说明中,参考了附图,附图形成本文的一部分。在附图中,类似的符号典型地确定类似的部件,除非上下文以其他方式指明。在详细的说明书、附图及权利要求书中所描述的图示说明的实施方案不意味是限制性的。在不脱离本文所呈现的主题的精神或范围下,其他实施方案可以被使用,并且可以作其他改变。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种三维气动外形反设计方法,其特征在于,所述方法包括:
获取每个样本点的参数化变形控制点和气动外形;
根据所述样本点的参数化变形控制点和所述气动外形,建立高维数据;
采用ISOMAP非线性降维方法,将所述高维数据降到低维流形空间的低维数据;
在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;具体包括:
在低维流形空间中,获取目标参数化变形控制点:
根据气功外形的流形结构和所述目标参数化变形控制点的流形空间坐标值进行RBF插值,获得所述目标参数化变形控制点的目标压力分布值,其中,所述气动外形的流形结构为气功外形在低维空间中的坐标值;
根据所述目标压力分布值,通过POD降维模型,确定所述目标压力分布值下的设计变量值;包括:
将所述目标参数化变形控制点的目标压力分布值和对应的设计变量进行组合,获得训练数据集;
将所述训练数据集与所述目标参数化变形控制点的目标压力分布值进行组合,确定缺失目标设计变量的数据集;
通过POD阵维模型迭代出符合参数化变形控制点对应的目标压力分布的设计变量值;
根据所述目标压力分布值下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形。
2.根据权利要求1所述的方法,其特征在于,所述采用ISOMAP非线性降维方法,将所述高维数据降到低维流形空间的低维数据,包括:
获取N个D维数据的样本;
根据所述样本中任意两点,确定领域图;
在所述领域图中,根据Dijkstra算法或Floyd算法,计算所述领域图中所有点对的最短距离,并将得到的最短距离的线段补充在所述领域图中,得到重构的领域图;
采用MDS算法对所述重构的领域图进行降维处理,得到低维流形空间的低维数据。
3.根据权利要求2所述的方法,其特征在于,所述根据所述样本中任意两点,确定领域图,包括:
计算样本中任意两点i和j的距离;
判断i和j是否相邻,若i和j相邻,则将i和j两点连接起来,计算两点之间的长度;
根据满足条件的所有点,确定领域图。
4.根据权利要求2所述的方法,其特征在于,所述采用MDS算法对重构的领域图进行降维处理,得到低维流形空间的低维数据,包括:
根据所述重构领域图,确定内积矩阵;
计算所述内积矩阵的特征值和特征向量;
根据所述特征向量,确定原始数据集在低维流形上的坐标。
5.根据权利要求4所述的方法,其特征在于,所述根据所述目标压力分布值下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形,包括:
根据所述目标压力分布值下的设计变量值,利用FFD参数化方法对原始气动外形进行参数化处理,得到目标气动外形;
在通过网格变形方法获得目标气动外形对应的计算网格;
利用CFD方法,确定所述目标气动外形下的压力分布值,并根据所述压力分布值,确定与所述设计变量值对应的压力分布;
若所述压力分布满足预设条件,则得到相应的三维气动外形。
6.一种三维气动外形反设计装置,其特征在于,所述装置包括:
获取模块,用于获取每个样本点的参数化变形控制点和气动外形;
建立模块,用于根据所述样本点的参数化变形控制点和所述气动外形,建立高维数据;
降维模块,用于采用ISOMAP非线性降维方法,将所述高维数据降到低维流形空间的低维数据;
插值模块,用于在低维流形空间中,获取目标参数化变形控制点,对所述目标参数化变形控制点进行压力系数插值,获得与所述目标参数化变形控制点对应的目标压力分布值;具体包括:
在低维流形空间中,获取目标参数化变形控制点:
根据气功外形的流形结构和所述目标参数化变形控制点的流形空间坐标值进行RBF插值,获得所述目标参数化变形控制点的目标压力分布值,其中,所述气动外形的流形结构为气功外形在低维空间中的坐标值;
确定模块,用于根据所述目标压力分布值,通过POD降维模型,确定所述目标压力分布值下的设计变量值;包括:
将所述目标参数化变形控制点的目标压力分布值和对应的设计变量进行组合,获得训练数据集;
将所述训练数据集与所述目标参数化变形控制点的目标压力分布值进行组合,确定缺失目标设计变量的数据集;
通过POD阵维模型迭代出符合参数化变形控制点对应的目标压力分布的设计变量值;
设计模块,用于根据所述目标压力分布值下的设计变量值,通过FFD方法对原始气动外形进行参数化,得到相应的三维气动外形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211093692.XA CN115238396B (zh) | 2022-09-08 | 2022-09-08 | 三维气动外形反设计方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211093692.XA CN115238396B (zh) | 2022-09-08 | 2022-09-08 | 三维气动外形反设计方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115238396A CN115238396A (zh) | 2022-10-25 |
CN115238396B true CN115238396B (zh) | 2022-12-30 |
Family
ID=83681370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211093692.XA Active CN115238396B (zh) | 2022-09-08 | 2022-09-08 | 三维气动外形反设计方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115238396B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109899248A (zh) * | 2019-03-04 | 2019-06-18 | 天津工业大学 | 一种基于多项式拟合的水平轴风力机叶片霜冰冰形形成方法 |
CN114355777A (zh) * | 2022-01-06 | 2022-04-15 | 浙江大学 | 一种基于分布式压力传感器与分段姿态控制的动态滑翔方法及*** |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101655990B (zh) * | 2009-06-25 | 2012-06-06 | 浙江大学 | 一种基于非线性流形学习的三维人体运动合成方法 |
EP3209803A4 (en) * | 2014-10-21 | 2018-06-13 | Ubiome, Inc. | Method and system for microbiome-derived diagnostics and therapeutics |
CN106960000A (zh) * | 2017-01-23 | 2017-07-18 | 中山大学 | 一种基于测地线保持的非线性数据降维方法 |
CN107066741B (zh) * | 2017-04-18 | 2020-03-10 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于数据挖掘的气动外形优化设计方法 |
CN111340685A (zh) * | 2020-02-14 | 2020-06-26 | 中国地质大学(武汉) | 一种用于遥感数据处理的流形降维方法 |
CN111553034B (zh) * | 2020-04-28 | 2022-07-08 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于插值方法的自由变形参数化方法 |
CN111783223B (zh) * | 2020-06-17 | 2023-11-10 | 中国空气动力研究与发展中心 | 一种面向超声速飞行器远场声爆抑制的近场伴随优化方法 |
CN112016167B (zh) * | 2020-10-22 | 2021-01-29 | 中国人民解放军国防科技大学 | 基于仿真和优化耦合的飞行器气动外形设计方法及*** |
CN113268818A (zh) * | 2021-07-19 | 2021-08-17 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于生成拓扑映射的气动全局优化方法、存储介质及终端 |
CN114077771A (zh) * | 2022-01-07 | 2022-02-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 气动外形反设计方法 |
CN114036650B (zh) * | 2022-01-07 | 2022-03-15 | 中国空气动力研究与发展中心计算空气动力研究所 | 非轴对称自然层流短舱反设计方法及自然层流短舱 |
CN114154275B (zh) * | 2022-02-10 | 2022-04-12 | 北京航空航天大学 | 基于最佳负荷分布模型优化的低压涡轮叶型气动设计方法 |
CN114925539B (zh) * | 2022-06-01 | 2022-12-30 | 中国科学院计算机网络信息中心 | 声爆波反演方法、装置、设备及存储介质 |
-
2022
- 2022-09-08 CN CN202211093692.XA patent/CN115238396B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109899248A (zh) * | 2019-03-04 | 2019-06-18 | 天津工业大学 | 一种基于多项式拟合的水平轴风力机叶片霜冰冰形形成方法 |
CN114355777A (zh) * | 2022-01-06 | 2022-04-15 | 浙江大学 | 一种基于分布式压力传感器与分段姿态控制的动态滑翔方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN115238396A (zh) | 2022-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106384384B (zh) | 一种三维产品模型的形状优化方法 | |
CN101013449A (zh) | 利用交叉模型交叉模态的结构模型修正方法 | |
CN112862956B (zh) | 基于HRBFs的人体与服装模型碰撞检测和处理方法 | |
Chen et al. | Deep generative model for efficient 3D airfoil parameterization and generation | |
CN110334450B (zh) | 一种多块结构网格生成中物面投影错误的修复方法 | |
Lin et al. | Vertex-ball spring smoothing: an efficient method for unstructured dynamic hybrid meshes | |
Xiong et al. | A point cloud deep neural network metamodel method for aerodynamic prediction | |
Li et al. | An efficient multiple point selection study for mesh deformation using radial basis functions | |
CN115238396B (zh) | 三维气动外形反设计方法及装置 | |
WO2022046889A1 (en) | Flex representation in computer aided design and computer aided engineering | |
Strofylas et al. | An agglomeration strategy for accelerating RBF-based mesh deformation | |
CN116342836B (zh) | 一种基于深度隐式细节表征的人体建模方法 | |
Cao et al. | Computation of the medial axis of planar domains based on saddle point programming | |
Wang et al. | Kirchhoff–Love shell representation and analysis using triangle configuration b-splines | |
Poole et al. | Optimal domain element shapes for free-form aerodynamic shape control | |
Han et al. | Adjoint method for aerodynamic shape improvement in comparison with surface pressure gradient method | |
Wei et al. | Latent Representation of CFD Meshes and Application to Airfoil Aerodynamics | |
CA3009499A1 (en) | Composite design direction | |
Xiao et al. | Three-dimensional aerodynamic shape inverse design based on ISOMAP | |
CN112182739A (zh) | 一种飞行器结构非概率可信可靠性拓扑优化设计方法 | |
Bhat et al. | Constraint-Based Parameterization and Disentanglement of Aerodynamic Shapes Using Deep Generative Models | |
Peng et al. | Hexahedral Mesh Quality Improvement with Geometric Constraints | |
CN113626893B (zh) | 一种基于隐式几何模型的计算机力学分析数值模拟方法 | |
CN116579269B (zh) | 多面体网格原位分布式并行边界层物理量获取方法和装置 | |
CN117291123B (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 |