CN115994410A - 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 - Google Patents

Abstract

本申请涉及一种基于八叉树细化四面体网格的飞行器仿真驱动设计方法。所述方法包括：根据几何参数曲面以飞行器为中心建立远场的长方体包围盒，得到飞行器的三维几何模型；利用八叉树细化四面体网格对所述飞行器的三维几何模型进行网格化处理，得到三维几何模型的体网格和面网格；将所述体网格和面网格的信息构建格式网格文件，利用CFD数值仿真软件对所述格式网格文件进行计算，得到飞行器的多个性能指标；所述性能指标包括飞行器升力系数和飞行器阻力系数；根据预先设置的飞行器设计要求对所述飞行器的多个性能指标进行评估，若性能指标达到设计要求则输出设计方案。采用本方法能够缩短并优化飞行器设计的仿真周期。

Description

基于八叉树细化四面体网格的飞行器仿真驱动设计方法

技术领域

本申请涉及航空飞行器的仿真驱动技术领域，特别是涉及一种基于八叉树细化四面体网格的飞行器仿真驱动设计方法。

背景技术

随着仿真驱动技术的发展，数值仿真技术已经成为飞行器气动外形设计过程中性能分析的重要手段，在飞行器外形设计完成后，一般会根据设计的几何模型增加远场包围，进而生成可供CFD数值仿真软件计算的网格文件，在运用有限体积法或有限元法等CFD数值仿真方法计算后评估并优化飞行器外形设计方案。其中，CFD数值仿真计算的一个关键前提是需要通过网格生成技术首先将飞行器模型对应的计算域进行离散化表示，目的是将给定飞行器外流场区域划分成有限个基本几何单元，一般将三维几何区域划分为四面体或六面体单元的网格。在飞行器外流场计算区域转化成网格表示后，CFD方法在每个网格单元上存储温度、速度等物理量，然后利用物理守恒方程求解给定工况下的物理解,即各网格单元上的物理量，并通过相关计算获得飞行器升力系数、阻力系数等性能指标。因此，生成网格的单元尺寸分布和单元质量对仿真结果影响巨大，从而影响飞行器外流场设计方案的性能。此外网格生成方法的鲁棒性也往往决定着网格生成的成败，从而极大影响着飞行器设计方案的仿真周期。

然而，飞行器驱动设计方案的优化依赖CFD计算结果，作为CFD数值仿真计算的一个关键前提，通常的非结构四面体网格生成方法大多依照先生成表面网格，再基于表面网格生成体网格。因此若几何数模存在缺陷就很难生成封闭的表面网格，导致体网格生成失败，从而使得生成的飞行器驱动设计方案不准确。目前的飞行器驱动设计方案的Delaunay网格剖分法需要逐个插入新的节点到当前已有的网格，阵面推进法从边界出发逐层向区域内部生成网格。这两种方法生成三维空间的四面体网格都需要首先生成边界表面网格，在内部再逐点计算出新增网格。面对大型几何模型的高精度仿真往往会出现的三维网格生成计算量过大的情况，在生成效率方面难以有较大突破，飞行器仿真驱动设计周期长。

发明内容

基于此，有必要针对上述技术问题，提供一种能够缩短并优化飞行器设计仿真周期的基于八叉树细化四面体网格的飞行器仿真驱动设计方法。

一种基于八叉树细化四面体网格的飞行器仿真驱动设计方法，所述方法包括：

对飞行器的几何形状进行设计，得到飞行器的几何模型文件；几何模型文件包括飞机部件的几何参数曲面；飞机部件包括机身、机翼和尾翼；几何参数曲面包括曲面的阶数、控制点数、节点向量和控制点坐标；

从几何模型文件中获取飞行器的几何参数曲面，根据几何参数曲面以飞行器为中心建立远场的长方体包围盒，得到飞行器的三维几何模型；

利用八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理，得到三维几何模型的体网格和面网格；

将体网格和面网格的信息构建格式网格文件，利用CFD数值仿真软件对格式网格文件进行计算，得到飞行器的多个性能指标；性能指标包括飞行器升力系数和飞行器阻力系数；

根据预先设置的飞行器设计要求对飞行器的多个性能指标进行评估，若性能指标达到设计要求则输出设计方案。

在其中一个实施例中，若性能指标未达到设计要求则由设计人员对飞行器的几何形状重新进行模型调整并进行性能指标评估，直至性能指标达到设计要求输出设计方案。

在其中一个实施例中，长方体包围盒的远场长度为4倍机身长度，高度为2倍机身长度，宽度为1.5倍机身长度；利用八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理，得到三维几何模型的体网格和面网格，包括：

定义一个固定形状的初始根网格四面体，其中初始根网格四面体的四个顶点坐标为：；

获取三维几何模型的坐标范围，根据三维几何模型的坐标范围等比例更改初始根四面体尺寸至根网格可将几何模型完全包围为止，得到更改后的初始根网格四面体；

对更改后的初始根网格四面体作八叉树细化网格分裂，判断当前分裂得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限，直至分裂得到最终的四面体网格单元集；

从四面体网格单元集中搜索出与三维几何边界相交的四面体网格单元、与几何边界相交单元的相邻单元以及几何体内部单元，保留这些单元并删除四面体网格单元集中其他单元，得到接近飞行器及外流场几何外形的初始八叉树网格；

根据飞行器的几何参数曲面对初始八叉树网格进行移动网格节点、交换边和切割网格单元，使得网格与三维几何表面完全贴合，得到三维几何模型的体网格和面网格。

在其中一个实施例中，对更改后的初始根网格四面体作八叉树细化网格分裂，直至得到达到预先设置尺寸的四面体，包括：

在当前四面体每个面上连接四面体边的中点进行切割，得到顶点处四个子四面体和一个核心几何体；连接核心几何体的两个顶点进行分裂，再得到四个子四面体；子四面体形状与原四面体完全相同；子四面体表示为，其中包含子四面体四个顶点，分别为子四面体网格单元的坐标值；每次对当前所有子四面体作一次八叉树细化网格分裂，增加一次八叉树网格细化层数；

判断当前分裂得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限，直至分裂得到最终的四面体网格单元集。

在其中一个实施例中，根据飞行器的几何参数曲面对初始八叉树网格进行移动网格节点、交换边和切割网格单元，使得网格与三维几何表面完全贴合，得到三维几何模型的体网格和面网格，包括：

根据飞行器的几何参数曲面确定初始八叉树网格中包含几何点的四面体网格单元，对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量平均值；

取每个四面体网格单元顶点中的移动后的网格质量平均值最大的顶点作为最佳可移动顶点，将最佳可移动顶点坐标修改为内部几何点坐标值；

根据飞行器的几何参数曲面确定初始八叉树网格中与飞行器的几何边相交的四面体网格单元，计算与几何边相交的四面体网格单元中每个顶点与顶点相对网格面的距离再乘以预设系数，得到每个四面体顶点的最大可移动距离；

获取飞行器的几何边的离散点，通过对比离散点与四面体顶点的距离得到每个顶点在几何边上的投影点，再计算出四个四面体顶点到几何边的距离；

对比四面体顶点与几何边的距离和四面体顶点的最大可移动距离，若四面体顶点与几何边的距离小于四面体顶点的最大可移动距离，更改四面体顶点坐标为投影点坐标值上，否则坐标不变；

根据几何参数曲面得到初始八叉树网格中与飞行器的几何面相交的四面体网格单元，对比四面体顶点与几何面的距离和此顶点的最大可移动距离，若顶点与几何面的距离小于此点的最大可移动距离，更改此顶点坐标为投影点坐标，否则坐标不变。

在其中一个实施例中，若四面体网格与几何面相交，则更改四面体关联的顶点，删除相交的四面体形成空腔，在空腔内使同一几何面上的顶点建立连接关系，形成新的四面体簇；

在体网格中根据几何参数曲面和初始八叉树网格得到所有与飞行器的几何线或几何面相交的四面体网格边、四面体网格面和四面体网格体。

定义分割边、分割面和分割体，并记录相交位置，删除包含，，元素的四面体网格单元，在相交位置处建立新的连接关系，更新四面体簇网格，得到体网格；

删除飞行器几何体外的网格单元，提取四面体网格单元与飞行器几何体贴合的三角形面单元形成三维模型的面网格。

在其中一个实施例中，对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量的平均值，包括：

对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量的平均值为

；

其中，为相邻网格数；为相邻网格编号；为相邻四面体网格体积；为相邻四面体外接圆体积。

上述基于八叉树细化四面体网格的飞行器仿真驱动设计方法，本申请在飞行器几何外形设计后，能够对输入的几何数模作四面体网格生成，将体网格和面网格的信息构建格式网格文件，利用CFD数值仿真软件对格式网格文件进行计算，得到飞行器的多个性能指标，然后对性能指标进行评估并反馈结果，形成完整的飞行器仿真驱动设计方案，通过八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理，减少了得到飞行器三维几何模型的体网格和面网格的计算量和人工交互，计算效率更高，鲁棒性好，并且四面体形状更接近于正四面体，网格质量更好，从而减少了飞行器仿真驱动的设计周期。

附图说明

图1为一个实施例中一种基于八叉树细化四面体网格的飞行器仿真驱动设计方法的流程示意图；

图2为一个实施例中飞行器的几何模型图；

图3为一个实施例中面网格生成的流程图；

图4为一个实施例中八叉树根网格分裂方式示意图；

图5为一个实施例中网格顶点移动到几何点示意图；（a）为网格顶点移动前示意图，（b）为网格顶点移动后示意图；

图6为一个实施例中网格顶点投影到几何面示意图；（a）为被几何面横穿的四面体网格示意图，（b）为符合条件的网格顶点投影到几何面后的效果图；

图7为一个实施例中网格边交换示意图；（a）为网格边交换前示意图，（b）为网格边交换后示意图；

图8为一个实施例中飞行器经八叉树细化网格生成后的表面网格图；

图9为一个实施例中飞行器经八叉树细化网格生成后的体网格图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。

在一个实施例中，如图1所示，提供了一种基于八四面体网格的飞行器仿真驱动设计方法，包括以下步骤：

步骤102，对飞行器的几何形状进行设计，得到飞行器的几何模型文件；几何模型文件包括飞机部件的几何参数曲面；飞机部件包括机身、机翼和尾翼；几何参数曲面包括曲面的阶数、控制点数、节点向量和控制点坐标；从几何模型文件中获取飞行器的几何参数曲面，根据几何参数曲面以飞行器为中心建立远场的长方体包围盒，得到飞行器的三维几何模型。

根据几何参数曲面以飞行器为中心建立远场的长方体包围盒，得到飞行器的三维几何模型如图2所示，

步骤104，利用八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理，得到三维几何模型的体网格和面网格。

如图3和图4所示，利用八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理主要包括以下步骤：

3.1、定义一个固定形状的初始根网格四面体，四个顶点坐标为：(1,0,0),(-1,0,0),(0,1,1),(0,-1,1)。

3.2、获取三维几何模型M的最大坐标和最小坐标，保持初始初始根网格四面体形状不变，根据其坐标范围等比例更改根四面体尺寸至根网格可将几何模型完全包围为止。更新四面体的顶点坐标。

3.3、对当前已有四面体网格作八叉树细化网格分裂，方法如图3所示，以一个四面体为例，在当前四面体每个面上连接四面体边的中点，可以切割得到顶点处四个子四面体和一个核心几何体，连接核心几何体的两个顶点，可再分裂得到四个子四面体，共得到八个子四面体，形状与原四面体完全相同。细化后得到四面体网格单元，其中包含四面体四个顶点，分别为四面体网格单元的坐标值。每次对当前所有四面体作一次八叉树细化网格分裂，增加一次八叉树网格细化层数。

3.4、判断当前分裂得到的四面体大小是否达到预设尺寸，或细化层数是否达到预设上限，若达到，转至3.5；若未达到转至3.3。

3.5、搜索出与三维几何边界相交的四面体网格单元、与几何边界相交单元的相邻单元以及几何体内部单元，保留这些单元并删除其他单元，得到接近飞行器及外流场几何外形的初始八叉树网格，包含筛选后的四面体网格单元。

3.6、通过移动网格节点、交换边、切割网格单元三个步骤调整几何边界上的网格，使得网格与三维几何表面完全贴合，初始八叉树网格调整后得到网格。具体步骤为：

3.6.1、根据几何参数曲面得到已有网格中包含几何点的四面体网格单元，计算四面体每个顶点移动后所有相邻网格单元的网格质量的平均值，这里取四面体体积占其外接圆体积的比例为网格质量标准。

取每个四面体网格单元顶点中的移动后值最大的顶点作为最佳可移动顶点。将这一顶点坐标修改为内部几何点坐标值。

3.6.2、根据几何参数曲面得到已有网格中与几何边相交的四面体网格单元，计算这些四面体网格单元每个顶点与其相对网格面的距离，再乘以预设系数得到每个四面体顶点的最大可移动距离为

；

式中，为四面体顶点最大可移动距离；为预设的距离系数，默认设为0.35，可根据计算需求调整；为顶点到四面体相对网格面的垂线距离。

获取几何边的离散点，通过对比离散点与四面体顶点的距离得到每个顶点在几何边上的投影点，进而计算出四个四面体顶点到几何边的距离；

对比四面体顶点与几何边的距离和此顶点的最大可移动距离，若顶点与几何边的距离小于此点的最大可移动距离，更改此顶点坐标为投影点坐标值上，否则坐标不变。

此步骤执行后的效果如图5所示，左图为包含几何点的网格及相邻网格示意图，右图为网格顶点移动到几何点的网格示意图。

3.6.3、根据几何参数曲面得到已有网格中与几何面相交的四面体网格单元，对比四面体顶点与几何面的距离和此顶点的最大可移动距离，若顶点与几何面的距离小于此点的最大可移动距离，更改此顶点坐标为投影点坐标，否则坐标不变。

此步骤执行后的效果如图6所示，（a）图为被几何面横穿的四面体网格，（b）为符合条件的网格顶点投影到几何面后的效果图。

3.6.4、针对四面体网格与几何面相交的情况，更改四面体关联的顶点。方法为删除相交的四面体形成空腔，在空腔内使同一几何面上的顶点建立连接关系，形成新的四面体簇，更新得到新的网格。

此步骤执行后的效果如图7所示，四面体网格边与几何面相交的情况在交换边操作后消除了四面体被几何面分割的现象，且在几何面上生成了新的三角形网格；

3.6.5、在网格中根据几何参数曲面和已有网格得到所有与几何线或几何面相交的四面体网格边，网格面和网格体。定义为分割边、分割面和分割体，并记录相交位置；

3.6.6、删除包含，，元素的四面体网格单元，在相交位置处建立新的连接关系，更新网格。图8所示结果为飞行器外流场三维几何模型的表面网格生成结果，经过网格与三维几何边界的拟合操作，在处理体网格边界的同时也得到了三角形表面网格。

3.7、删除几何体外的网格单元，用体单元与几何体贴合的三角网格形成三维模型的面网格。

本申请通过对比网格节点最大可移动距离来判断网格节点是否要移动的算法，最大可移动距离由网格节点与四面体内相对网格面的距离得到，这种方法可以有效调整网格在几何边界处的形状，并可以避免点移动导致的网格错误。再通过几何点拟合、几何线拟合、几何面拟合来修正几何边界与网格的交叉，实现网格在几何边界处的调整。这种方式不仅简化了传统网格生成算法在边界处的处理步骤，并且可以使边界处网格形状更接近于内部，提升了网格质量，并且本申请无需逐点计算，计算量需求小，减少了得到三维几何模型的体网格和面网格的计算量，计算效率更高，鲁棒性好，并且四面体形状更接近于正四面体，网格质量更好，从而缩短了飞行器性能指标的计算过程，进而减少了飞行器仿真驱动的设计周期。

步骤106，将体网格和面网格的信息构建格式网格文件，利用CFD数值仿真软件对格式网格文件进行计算，得到飞行器的多个性能指标；性能指标包括飞行器升力系数和飞行器阻力系数。

将体网格和面网格构建成cgns格式的网格文件，图9所示结果为删除模型外部网格后的三维体网格生成结果。CFD软件导入上一步获得的cgns格式网格文件并计算获得飞行器升力系数、阻力系数等性能指标。

步骤108，根据预先设置的飞行器设计要求对飞行器的多个性能指标进行评估，若性能指标达到设计要求则输出设计方案。

在飞行器几何外形设计后，对输入的几何数模作四面体网格生成，在CFD计算后评估并反馈结果，形成完整的飞行器仿真驱动设计方案。

上述基于八叉树细化四面体网格的飞行器仿真驱动设计方法中，本申请在飞行器几何外形设计后，能够对输入的几何数模作四面体网格生成，将体网格和面网格的信息构建格式网格文件，利用CFD数值仿真软件对格式网格文件进行计算，得到飞行器的多个性能指标，然后对性能指标进行评估并反馈结果，形成完整的飞行器仿真驱动设计方案，通过八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理，减少了得到飞行器的三维几何模型的体网格和面网格的计算量和人工交互，计算效率更高，鲁棒性好，并且四面体形状更接近于正四面体，网格质量更好，从而减少了飞行器仿真驱动的设计周期。

在其中一个实施例中，若性能指标未达到设计要求则由设计人员对飞行器的几何形状重新进行模型调整并进行性能指标评估，直至性能指标达到设计要求输出设计方案。

在其中一个实施例中，长方体包围盒的远场长度为4倍机身长度，高度为2倍机身长度，宽度为1.5倍机身长度；利用八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理，得到三维几何模型的体网格和面网格，包括：

定义一个固定形状的初始根网格四面体，其中初始根网格四面体的四个顶点坐标为：；

获取三维几何模型的坐标范围，根据三维几何模型的坐标范围等比例更改初始根四面体尺寸至根网格可将几何模型完全包围为止，得到更改后的初始根网格四面体；

对更改后的初始根网格四面体作八叉树细化网格分裂，判断当前分裂得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限，直至分裂得到最终的四面体网格单元集；

从四面体网格单元集中搜索出与三维几何边界相交的四面体网格单元、与几何边界相交单元的相邻单元以及几何体内部单元，保留这些单元并删除四面体网格单元集中其他单元，得到接近飞行器及外流场几何外形的初始八叉树网格；

根据飞行器的几何参数曲面对初始八叉树网格进行移动网格节点、交换边和切割网格单元，使得网格与三维几何表面完全贴合，得到三维几何模型的体网格和面网格。

在其中一个实施例中，对更改后的初始根网格四面体作八叉树细化网格分裂，直至得到达到预先设置尺寸的四面体，包括：

在当前四面体每个面上连接四面体边的中点进行切割，得到顶点处四个子四面体和一个核心几何体；连接核心几何体的两个顶点进行分裂，再得到四个子四面体；子四面体形状与原四面体完全相同；子四面体表示为，其中包含子四面体四个顶点，分别为子四面体网格单元的坐标值；每次对当前所有子四面体作一次八叉树细化网格分裂，增加一次八叉树网格细化层数；

判断当前分裂得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限，直至分裂得到最终的四面体网格单元集。

在其中一个实施例中，根据飞行器的几何参数曲面对初始八叉树网格进行移动网格节点、交换边和切割网格单元，使得网格与三维几何表面完全贴合，得到三维几何模型的体网格和面网格，包括：

根据飞行器的几何参数曲面确定初始八叉树网格中包含几何点的四面体网格单元，对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量平均值；

取每个四面体网格单元顶点中的移动后的网格质量平均值最大的顶点作为最佳可移动顶点，将最佳可移动顶点坐标修改为内部几何点坐标值；

根据飞行器的几何参数曲面确定初始八叉树网格中与飞行器的几何边相交的四面体网格单元，计算与几何边相交的四面体网格单元中每个顶点与顶点相对网格面的距离再乘以预设系数，得到每个四面体顶点的最大可移动距离；

获取飞行器几何边的离散点，通过对比离散点与四面体顶点的距离得到每个顶点在几何边上的投影点，再计算出四个四面体顶点到几何边的距离；

对比四面体顶点与几何边的距离和四面体顶点的最大可移动距离，若四面体顶点与几何边的距离小于四面体顶点的最大可移动距离，更改四面体顶点坐标为投影点坐标值上，否则坐标不变；

根据几何参数曲面得到初始八叉树网格中与飞行器的几何面相交的四面体网格单元，对比四面体顶点与几何面的距离和此顶点的最大可移动距离，若顶点与几何面的距离小于此点的最大可移动距离，更改此顶点坐标为投影点坐标，否则坐标不变。

在其中一个实施例中，若四面体网格与几何面相交，则更改四面体关联的顶点，删除相交的四面体形成空腔，在空腔内使同一几何面上的顶点建立连接关系，形成新的四面体簇；

在体网格中根据几何参数曲面和初始八叉树网格得到所有与飞行器的几何线或几何面相交的四面体网格边、四面体网格面和四面体网格体。

定义分割边、分割面和分割体，并记录相交位置，删除包含，，元素的四面体网格单元，在相交位置处建立新的连接关系，更新四面体簇网格，得到体网格；

删除飞行器几何体外的网格单元，提取四面体网格单元与飞行器几何体贴合的三角形面单元形成三维模型的面网格。

在其中一个实施例中，对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量的平均值，包括：

对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量的平均值为

；

其中，为相邻网格数；为相邻网格编号；为相邻四面体网格体积；为相邻四面体外接圆体积。

应该理解的是，虽然图1的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。

以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。

Claims (7)

1.一种基于八叉树细化四面体网格的飞行器仿真驱动设计方法，其特征在于，所述方法包括：

对飞行器的几何形状进行设计，得到飞行器的几何模型文件；所述几何模型文件包括飞机部件的几何参数曲面；所述飞机部件包括机身、机翼和尾翼；所述几何参数曲面包括曲面的阶数、控制点数、节点向量和控制点坐标；

从所述几何模型文件中获取飞行器的几何参数曲面，根据所述几何参数曲面以飞行器为中心建立远场的长方体包围盒，得到飞行器的三维几何模型；

利用八叉树细化四面体网格对所述飞行器的三维几何模型进行网格化处理，得到三维几何模型的体网格和面网格；

将所述体网格和面网格的信息构建格式网格文件，利用CFD数值仿真软件对所述格式网格文件进行计算，得到飞行器的多个性能指标；所述性能指标包括飞行器升力系数和飞行器阻力系数；

根据预先设置的飞行器设计要求对所述飞行器的多个性能指标进行评估，若性能指标达到设计要求则输出设计方案。

2.根据权利要求1所述的方法，其特征在于，所述方法还包括：

若性能指标未达到设计要求则由设计人员对飞行器的几何形状重新进行模型调整并进行性能指标评估，直至性能指标达到设计要求输出设计方案。

3.根据权利要求1所述的方法，其特征在于，所述长方体包围盒的远场长度为4倍机身长度，高度为2倍机身长度，宽度为1.5倍机身长度；利用八叉树细化四面体网格对所述飞行器的三维几何模型进行网格化处理，得到三维几何模型的体网格和面网格，包括：

定义一个固定形状的初始根网格四面体，其中初始根网格四面体的四个顶点坐标为：；

获取三维几何模型的坐标范围，根据所述三维几何模型的坐标范围等比例更改初始根四面体尺寸至根网格可将几何模型完全包围为止，得到更改后的初始根网格四面体；

对所述更改后的初始根网格四面体作八叉树细化网格分裂，判断当前分裂得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限，直至分裂得到最终的四面体网格单元集；

从所述四面体网格单元集中搜索出与三维几何边界相交的四面体网格单元、与几何边界相交单元的相邻单元以及几何体内部单元，保留这些单元并删除四面体网格单元集中其他单元，得到接近飞行器及外流场几何外形的初始八叉树网格；

根据飞行器的几何参数曲面对所述初始八叉树网格进行移动网格节点、交换边和切割网格单元，使得网格与三维几何表面完全贴合，得到三维几何模型的体网格和面网格。

4.根据权利要求3所述的方法，其特征在于，对所述更改后的初始根网格四面体作八叉树细化网格分裂，直至得到达到预先设置尺寸的四面体，包括：

在当前四面体每个面上连接四面体边的中点进行切割，得到顶点处四个子四面体和一个核心几何体；连接所述核心几何体的两个顶点进行分裂，再得到四个子四面体；所述子四面体形状与原四面体完全相同；所述子四面体表示为，其中包含子四面体四个顶点，分别为子四面体网格单元的坐标值；每次对当前所有子四面体作一次八叉树细化网格分裂，增加一次八叉树网格细化层数；

判断当前分裂得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限，直至分裂得到最终的四面体网格单元集。

5.根据权利要求3所述的方法，其特征在于，根据飞行器的几何参数曲面对所述初始八叉树网格进行移动网格节点、交换边和切割网格单元，使得网格与三维几何表面完全贴合，得到三维几何模型的体网格和面网格，包括：

根据飞行器的几何参数曲面确定所述初始八叉树网格中包含几何点的四面体网格单元，对所述四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量平均值；

取每个四面体网格单元顶点中的移动后的网格质量平均值最大的顶点作为最佳可移动顶点，将所述最佳可移动顶点坐标修改为内部几何点坐标值；

根据飞行器的几何参数曲面确定所述初始八叉树网格中与飞行器的几何边相交的四面体网格单元，计算所述与几何边相交的四面体网格单元中每个顶点与顶点相对网格面的距离再乘以预设系数，得到每个四面体顶点的最大可移动距离；

获取飞行器的几何边的离散点，通过对比离散点与四面体顶点的距离得到每个顶点在几何边上的投影点，再计算出四个四面体顶点到几何边的距离；

对比四面体顶点与几何边的距离和四面体顶点的最大可移动距离，若四面体顶点与几何边的距离小于四面体顶点的最大可移动距离，更改四面体顶点坐标为投影点坐标值上，否则坐标不变；

根据几何参数曲面得到初始八叉树网格中与飞行器的几何面相交的四面体网格单元，对比四面体顶点与几何面的距离和此顶点的最大可移动距离，若顶点与所述几何面的距离小于此点的最大可移动距离，更改此顶点坐标为投影点坐标，否则坐标不变。

6.根据权利要求5所述的方法，其特征在于，所述方法还包括：

若四面体网格与几何面相交，则更改四面体关联的顶点，删除相交的四面体形成空腔，在空腔内使同一几何面上的顶点建立连接关系，形成新的四面体簇；

在体网格中根据几何参数曲面和初始八叉树网格得到所有与飞行器的几何线或几何面相交的四面体网格边、四面体网格面和四面体网格体；

定义分割边、分割面和分割体，并记录相交位置，删除包含，，元素的四面体网格单元，在相交位置处建立新的连接关系，更新四面体簇网格，得到体网格；

删除飞行器几何体外的网格单元，提取四面体网格单元与飞行器几何体贴合的三角形面单元形成三维模型的面网格。

7.根据权利要求5所述的方法，其特征在于，对所述四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量的平均值，包括：

对所述四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算，得到网格质量的平均值为

；

其中，为相邻网格数；为相邻网格编号；为相邻四面体网格体积；为相邻四面体外接圆体积。

Images