CN115994410B - 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 - Google Patents
基于八叉树细化四面体网格的飞行器仿真驱动设计方法 Download PDFInfo
- Publication number
- CN115994410B CN115994410B CN202310283034.5A CN202310283034A CN115994410B CN 115994410 B CN115994410 B CN 115994410B CN 202310283034 A CN202310283034 A CN 202310283034A CN 115994410 B CN115994410 B CN 115994410B
- Authority
- CN
- China
- Prior art keywords
- grid
- aircraft
- geometric
- tetrahedron
- vertex
- 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 48
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004088 simulation Methods 0.000 title claims abstract description 29
- 238000005520 cutting process Methods 0.000 claims description 8
- 238000007670 refining Methods 0.000 claims description 7
- 238000001727 in vivo Methods 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000013499 data model Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000004898 kneading Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012216 screening Methods 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
- Processing Or Creating Images (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
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.3、对当前已有四面体网格作八叉树细化网格***,方法如图3所示,以一个四面体为例,在当前四面体每个面上连接四面体边的中点,可以切割得到顶点处四个子四面体和一个核心几何体,连接核心几何体的两个顶点,可再***得到四个子四面体,共得到八个子四面体,形状与原四面体完全相同。细化后得到四面体网格单元,其中/>包含四面体四个顶点/>,/>分别为四面体网格单元/>的坐标值。每次对当前所有四面体作一次八叉树细化网格***,增加一次八叉树网格细化层数。
3.4、判断当前***得到的四面体大小是否达到预设尺寸,或细化层数是否达到预设上限,若达到,转至3.5;若未达到转至3.3。
3.5、搜索出与三维几何边界相交的四面体网格单元、与几何边界相交单元的相邻单元以及几何体内部单元,保留这些单元并删除其他单元,得到接近飞行器及外流场几何外形的初始八叉树网格,包含筛选后的四面体网格单元/>。
3.6.2、根据几何参数曲面得到已有网格中与几何边相交的四面体网格单元,计算这些四面体网格单元每个顶点与其相对网格面的距离,再乘以预设系数得到每个四面体顶点的最大可移动距离为
此步骤执行后的效果如图5所示,左图为包含几何点的网格及相邻网格示意图,右图为网格顶点移动到几何点的网格示意图。
3.6.3、根据几何参数曲面得到已有网格中与几何面相交的四面体网格单元,对比四面体顶点与几何面的距离和此顶点的最大可移动距离/>,若顶点与几何面的距离小于此点的最大可移动距离,更改此顶点坐标为投影点坐标,否则坐标不变。
此步骤执行后的效果如图6所示,(a)图为被几何面横穿的四面体网格,(b)为符合条件的网格顶点投影到几何面后的效果图。
此步骤执行后的效果如图7所示,四面体网格边与几何面相交的情况在交换边操作后消除了四面体被几何面分割的现象,且在几何面上生成了新的三角形网格;
3.6.6、删除包含,/>,/>元素的四面体网格单元,在相交位置处建立新的/>连接关系,更新网格/>。图8所示结果为飞行器外流场三维几何模型的表面网格生成结果,经过网格与三维几何边界的拟合操作,在处理体网格边界的同时也得到了三角形表面网格。
本申请通过对比网格节点最大可移动距离来判断网格节点是否要移动的算法,最大可移动距离由网格节点与四面体内相对网格面的距离得到,这种方法可以有效调整网格在几何边界处的形状,并可以避免点移动导致的网格错误。再通过几何点拟合、几何线拟合、几何面拟合来修正几何边界与网格的交叉,实现网格在几何边界处的调整。这种方式不仅简化了传统网格生成算法在边界处的处理步骤,并且可以使边界处网格形状更接近于内部,提升了网格质量,并且本申请无需逐点计算,计算量需求小,减少了得到三维几何模型的体网格和面网格的计算量,计算效率更高,鲁棒性好,并且四面体形状更接近于正四面体,网格质量更好,从而缩短了飞行器性能指标的计算过程,进而减少了飞行器仿真驱动的设计周期。
步骤106,将体网格和面网格的信息构建格式网格文件,利用CFD数值仿真软件对格式网格文件进行计算,得到飞行器的多个性能指标;性能指标包括飞行器升力系数和飞行器阻力系数。
步骤108,根据预先设置的飞行器设计要求对飞行器的多个性能指标进行评估,若性能指标达到设计要求则输出设计方案。
在飞行器几何外形设计后,对输入的几何数模作四面体网格生成,在CFD计算后评估并反馈结果,形成完整的飞行器仿真驱动设计方案。
上述基于八叉树细化四面体网格的飞行器仿真驱动设计方法中,本申请在飞行器几何外形设计后,能够对输入的几何数模作四面体网格生成,将体网格和面网格的信息构建格式网格文件,利用CFD数值仿真软件对格式网格文件进行计算,得到飞行器的多个性能指标,然后对性能指标进行评估并反馈结果,形成完整的飞行器仿真驱动设计方案,通过八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理,减少了得到飞行器的三维几何模型的体网格和面网格的计算量和人工交互,计算效率更高,鲁棒性好,并且四面体形状更接近于正四面体,网格质量更好,从而减少了飞行器仿真驱动的设计周期。
在其中一个实施例中,若性能指标未达到设计要求则由设计人员对飞行器的几何形状重新进行模型调整并进行性能指标评估,直至性能指标达到设计要求输出设计方案。
在其中一个实施例中,长方体包围盒的远场长度为4倍机身长度,高度为2倍机身长度,宽度为1.5倍机身长度;利用八叉树细化四面体网格对飞行器的三维几何模型进行网格化处理,得到三维几何模型的体网格和面网格,包括:
获取三维几何模型的坐标范围,根据三维几何模型的坐标范围等比例更改初始根四面体尺寸至根网格可将几何模型完全包围为止,得到更改后的初始根网格四面体;
对更改后的初始根网格四面体作八叉树细化网格***,判断当前***得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限,直至***得到最终的四面体网格单元集;
从四面体网格单元集中搜索出与三维几何边界相交的四面体网格单元、与几何边界相交单元的相邻单元以及几何体内部单元,保留这些单元并删除四面体网格单元集中其他单元,得到接近飞行器及外流场几何外形的初始八叉树网格;
根据飞行器的几何参数曲面对初始八叉树网格进行移动网格节点、交换边和切割网格单元,使得网格与三维几何表面完全贴合,得到三维几何模型的体网格和面网格。
在其中一个实施例中,对更改后的初始根网格四面体作八叉树细化网格***,直至得到达到预先设置尺寸的四面体,包括:
在当前四面体每个面上连接四面体边的中点进行切割,得到顶点处四个子四面体和一个核心几何体;连接核心几何体的两个顶点进行***,再得到四个子四面体;子四面体形状与原四面体完全相同;子四面体表示为,其中/>包含子四面体四个顶点,/>分别为子四面体网格单元/>的坐标值/>;每次对当前所有子四面体作一次八叉树细化网格***,增加一次八叉树网格细化层数;
判断当前***得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限,直至***得到最终的四面体网格单元集。
在其中一个实施例中,根据飞行器的几何参数曲面对初始八叉树网格进行移动网格节点、交换边和切割网格单元,使得网格与三维几何表面完全贴合,得到三维几何模型的体网格和面网格,包括:
根据飞行器的几何参数曲面确定初始八叉树网格中包含几何点的四面体网格单元,对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算,得到网格质量平均值;
取每个四面体网格单元顶点中的移动后的网格质量平均值最大的顶点作为最佳可移动顶点,将最佳可移动顶点坐标修改为内部几何点坐标值;
根据飞行器的几何参数曲面确定初始八叉树网格中与飞行器的几何边相交的四面体网格单元,计算与几何边相交的四面体网格单元中每个顶点与顶点相对网格面的距离再乘以预设系数,得到每个四面体顶点的最大可移动距离;
根据几何参数曲面得到初始八叉树网格中与飞行器的几何面相交的四面体网格单元,对比四面体顶点与几何面的距离和此顶点的最大可移动距离/>,若顶点与几何面的距离小于此点的最大可移动距离,更改此顶点坐标为投影点坐标,否则坐标不变。
在其中一个实施例中,对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算,得到网格质量的平均值,包括:
对四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算,得到网格质量的平均值为
应该理解的是,虽然图1的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (6)
1.一种基于八叉树细化四面体网格的飞行器仿真驱动设计方法,其特征在于,所述方法包括:
对飞行器的几何形状进行设计,得到飞行器的几何模型文件;所述几何模型文件包括飞机部件的几何参数曲面;所述飞机部件包括机身、机翼和尾翼;所述几何参数曲面包括曲面的阶数、控制点数、节点向量和控制点坐标;
从所述几何模型文件中获取飞行器的几何参数曲面,根据所述几何参数曲面以飞行器为中心建立远场的长方体包围盒,得到飞行器的三维几何模型;
利用八叉树细化四面体网格对所述飞行器的三维几何模型进行网格化处理,得到三维几何模型的体网格和面网格;
将所述体网格和面网格的信息构建格式网格文件,利用CFD数值仿真软件对所述格式网格文件进行计算,得到飞行器的多个性能指标;所述性能指标包括飞行器升力系数和飞行器阻力系数;
根据预先设置的飞行器设计要求对所述飞行器的多个性能指标进行评估,若性能指标达到设计要求则输出设计方案;
利用八叉树细化四面体网格对所述飞行器的三维几何模型进行网格化处理,得到三维几何模型的体网格和面网格,包括:
定义一个固定形状的初始根网格四面体,其中根网格四面体的四个顶点坐标为:(1,0,0),(-1,0,0),(0,1,1),(0,-1,1);
获取三维几何模型的坐标范围,根据所述三维几何模型的坐标范围等比例更改初始根四面体尺寸至根网格可将几何模型完全包围为止,得到更改后的根网格四面体;
对所述更改后的根网格四面体作八叉树细化网格***,判断当前***得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限,直至***得到最终的四面体网格单元集;
从所述四面体网格单元集中搜索出与三维几何边界相交的四面体网格单元、与几何边界相交单元的相邻单元以及几何体内部单元,保留这些单元并删除四面体网格单元集中其他单元,得到接近飞行器及外流场几何外形的初始八叉树网格;
根据飞行器的几何参数曲面对所述初始八叉树网格进行移动网格节点、交换边和切割网格单元,使得网格与三维几何表面完全贴合,得到三维几何模型的体网格和面网格;
对所述更改后的根网格四面体作八叉树细化网格***,直至得到达到预先设置尺寸的四面体,包括:
在当前四面体每个面上连接四面体边的中点进行切割,得到顶点处四个子四面体和一个核心几何体;连接所述核心几何体的两个顶点进行***,再得到四个子四面体;所述子四面体形状与原四面体完全相同;所述子四面体表示为,其中/>包含子四面体四个顶点/>,/>分别为子四面体网格单元/>的坐标值;每次对当前所有子四面体作一次八叉树细化网格***,增加一次八叉树网格细化层数;
判断当前***得到的四面体大小是否达到预设尺寸或细化层数是否达到预设上限,直至***得到最终的四面体网格单元集。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
若性能指标未达到设计要求则由设计人员对飞行器的几何形状重新进行模型调整并进行性能指标评估,直至性能指标达到设计要求输出设计方案。
3.根据权利要求1所述的方法,其特征在于,所述长方体包围盒的远场长度为4倍机身长度,高度为2倍机身长度,宽度为1.5倍机身长度。
4.根据权利要求1所述的方法,其特征在于,根据飞行器的几何参数曲面对所述初始八叉树网格进行移动网格节点、交换边和切割网格单元,使得网格与三维几何表面完全贴合,得到三维几何模型的体网格和面网格,包括:
根据飞行器的几何参数曲面确定所述初始八叉树网格中包含几何点的四面体网格单元,对所述四面体网格单元中每个顶点移动后所有相邻网格单元的网格质量进行计算,得到网格质量平均值;
取每个四面体网格单元顶点中的移动后的网格质量平均值最大的顶点作为最佳可移动顶点,将所述最佳可移动顶点坐标修改为内部几何点坐标值;
根据飞行器的几何参数曲面确定所述初始八叉树网格中与飞行器的几何边相交的四面体网格单元,计算所述与几何边相交的四面体网格单元中每个顶点与顶点相对网格面的距离再乘以预设系数,得到每个四面体顶点的最大可移动距离;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310283034.5A CN115994410B (zh) | 2023-03-22 | 2023-03-22 | 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310283034.5A CN115994410B (zh) | 2023-03-22 | 2023-03-22 | 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115994410A CN115994410A (zh) | 2023-04-21 |
CN115994410B true CN115994410B (zh) | 2023-05-30 |
Family
ID=85992401
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310283034.5A Active CN115994410B (zh) | 2023-03-22 | 2023-03-22 | 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115994410B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116628899B (zh) * | 2023-07-25 | 2023-09-26 | 中国人民解放军国防科技大学 | 基于自动生成边界层网格的喷管气动外形设计方法及设备 |
CN117708904B (zh) * | 2023-11-08 | 2024-07-23 | 浙江大学 | 面向半导体数值模拟的混合单元网格生成方法 |
CN117473655B (zh) * | 2023-12-27 | 2024-03-15 | 中国人民解放军国防科技大学 | 基于边坍缩网格优化的飞行器仿真驱动设计方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102332179A (zh) * | 2010-09-20 | 2012-01-25 | 董福田 | 三维模型数据化简、渐进传输方法及装置 |
CN110457860A (zh) * | 2019-08-22 | 2019-11-15 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 网格的自动生成、装置、设备和存储介质 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7671858B1 (en) * | 2005-09-06 | 2010-03-02 | Sandia Corporation | Unconstrained paving and plastering method for generating finite element meshes |
EP2528042B1 (de) * | 2011-05-23 | 2019-12-18 | Deutsche Telekom AG | Verfahren und Vorrichtung zum Re-Meshing von 3D-Polygonmodellen |
CN102567594B (zh) * | 2012-02-20 | 2014-01-15 | 上海海洋大学 | 一种近岸岛礁型人工鱼礁群流场仿真建模的方法 |
CN104866648B (zh) * | 2015-04-21 | 2018-10-30 | 中车青岛四方机车车辆股份有限公司 | 一种高速行驶交通工具仿真模型前处理模板及处理方法 |
WO2019156947A1 (en) * | 2018-02-07 | 2019-08-15 | Siemens Corporation | System and method for design and manufacture of steady lattice structures |
US20200410147A1 (en) * | 2019-06-28 | 2020-12-31 | Viettel Group | Aerodynamic derivatives calculation method for flight vehicle |
CN110610050B (zh) * | 2019-09-18 | 2022-11-08 | 中国人民解放军国防科技大学 | 基于改进径向基函数变形算法的翼型气动减阻方法 |
CN112380795A (zh) * | 2020-12-09 | 2021-02-19 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种基于预置模型库的试车台流场数值仿真*** |
CN113379929B (zh) * | 2021-05-25 | 2023-02-28 | 西南交通大学 | 一种基于物理仿真的骨组织修复虚拟现实解决方法 |
CN115797600B (zh) * | 2023-02-08 | 2023-06-16 | 安世亚太科技股份有限公司 | 一种三维几何模型包面网格生成*** |
-
2023
- 2023-03-22 CN CN202310283034.5A patent/CN115994410B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102332179A (zh) * | 2010-09-20 | 2012-01-25 | 董福田 | 三维模型数据化简、渐进传输方法及装置 |
CN110457860A (zh) * | 2019-08-22 | 2019-11-15 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 网格的自动生成、装置、设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN115994410A (zh) | 2023-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115994410B (zh) | 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 | |
EP3092626B1 (en) | Method for creating three dimensional lattice structures in computer-aided design models for additive manufacturing | |
CN108647370B (zh) | 基于双环迭代的无人直升机气动外形优化设计方法 | |
CN110069800B (zh) | 具有光滑边界表达的三维结构拓扑优化设计方法及设备 | |
CN111581776B (zh) | 一种基于几何重建模型的等几何分析方法 | |
JPWO2007083602A1 (ja) | 補間処理方法、補間処理装置、形状評価方法、および形状評価装置 | |
US20150127301A1 (en) | Updating A CAD Model To Reflect Global Or Local Shape Changes | |
CN110689620B (zh) | 一种多层次优化的网格曲面离散样条曲线设计方法 | |
Zuo et al. | Isogeometric analysis for CSG models | |
CN104156546A (zh) | 基于t样条的汽车覆盖件模具的形面再设计方法 | |
Wang et al. | From computer-aided design (CAD) toward human-aided design (HAD): an isogeometric topology optimization approach | |
CN115481497A (zh) | 一种基于特征框架的体参数化建模方法 | |
CN115688276A (zh) | 一种基于离散伴随方法的飞行器外形自动化优化方法、***、设备、介质 | |
CN117473655B (zh) | 基于边坍缩网格优化的飞行器仿真驱动设计方法和装置 | |
US20100153074A1 (en) | Design support apparatus | |
CN108763668B (zh) | 基于细分技术与边界替换的齿轮模型区域参数化方法 | |
CN110532670B (zh) | 一种适分析样条实体模型构建方法及*** | |
CN114119628B (zh) | 一种基于特征模板的点云精确分割方法 | |
CN114741798A (zh) | 一种考虑电磁与力学性能的电机转子结构拓扑优化方法 | |
Shang et al. | Effective re-parameterization and GA based knot structure optimization for high quality T-spline surface fitting | |
CN115659480A (zh) | 飞机参数化建模和自动优化集成设计*** | |
CN112749438B (zh) | 基于二维虚拟断面的三维大坝坝体结构模型构建方法 | |
Navangul et al. | A vertex translation algorithm for adaptive modification of STL file in layered manufacturing | |
CN110765506B (zh) | 实体模型的多分辨率等几何拓扑优化方法 | |
Zhu et al. | A simple automatic hexahedron mesh generation and polyhedral smoothed finite element method for mechanics problems |
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 |