CN115270378A - 一种弓形激波外场网格的生成方法 - Google Patents

一种弓形激波外场网格的生成方法 Download PDF

Info

Publication number
CN115270378A
CN115270378A CN202211189574.9A CN202211189574A CN115270378A CN 115270378 A CN115270378 A CN 115270378A CN 202211189574 A CN202211189574 A CN 202211189574A CN 115270378 A CN115270378 A CN 115270378A
Authority
CN
China
Prior art keywords
grid
shock wave
target
point
node
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.)
Granted
Application number
CN202211189574.9A
Other languages
English (en)
Other versions
CN115270378B (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.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
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 Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202211189574.9A priority Critical patent/CN115270378B/zh
Publication of CN115270378A publication Critical patent/CN115270378A/zh
Application granted granted Critical
Publication of CN115270378B publication Critical patent/CN115270378B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Algebra (AREA)
  • Probability & Statistics with Applications (AREA)
  • Evolutionary Biology (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请公开了一种弓形激波外场网格的生成方法，在生成弓形激波的外场网格的过程中，在获取到模型网格文件的基础上，对模型网格文件表示出的目标模型是全模还是半模进行区分，无论目标模型的模态是全模还是半模，均能够高效地生成外场网格。进一步地，本申请中的弓形激波外场网格的生成方法从参考线、截面两个方面着手构建弓形激波的外场网格，有利于保障得到的外场网格的网格准确度和网格精度。并且，生成外场网格的过程无需过多的人工干预，有利于降低生成外场网格的门槛、降低人工成本。

Description

一种弓形激波外场网格的生成方法
技术领域
本申请属于气动数据研究领域，特别涉及一种弓形激波外场网格的生成方法。
背景技术
目标物（例如飞行器）在运动的过程中，由目标物的运动产生的能量传递给了目标物周边的空气，由此，在目标物的头部会形成一道弓形激波（也可以称作“头激波”）。弓形激波对目标物的运动会造成一定的影响，例如，气体通过激波后温度和压力大大增加而速度降低，激波层内高温气体，将以传导和辐射的方式把热量传递给目标物。同时，气体和目标物的表面之间有摩擦，不仅会使目标物减速，而且会产生热量并传递给周边的空气。可见，对弓形激波的研究是对目标物进行气动分析的较为重要的环节之一。飞行器（flightvehicle）是在大气层内或大气层外空间（太空）飞行的器械。飞行器分为3类：航空器、航天器、火箭和导弹。在大气层内飞行的称为航空器，如气球、飞艇、飞机等。它们靠空气的静浮力或空气相对运动产生的空气动力升空飞行。在太空飞行的称为航天器，如人造地球卫星、载人飞船、空间探测器、航天飞机等。它们在运载火箭的推动下获得必要的速度进入太空，然后依靠惯性做与天体类似的轨道运动。计算流体动力学方法现已成为飞行器气动研究的主流手段之一。
在现有的针对目标物、弓形激波网格生成软件中，网格生成均采用交互式的生成方式。即使是对于自动化程度比较高的非结构网格生成，在生成过程中也需要大量的人机实时交互，才能保证所生成的网格具有较好的质量，以满足后续的CFD流场求解和分析。也就是说，现有的网格生成方法存在大量的人机实时交互，导致网格生成非常耗时，由于人工干预较多也增加了人工成本。因此急需一种更加高效的网格生成方法。
发明内容
为了解决现有技术的不足，本申请提供了一种弓形激波外场网格的生成方法，在生成弓形激波的外场网格的过程中，在获取到模型网格文件的基础上，对模型网格文件表示出的目标模型是全模还是半模进行区分，无论目标模型的模态是全模还是半模，均能够高效地生成外场网格。进一步地，本申请中的弓形激波外场网格的生成方法从参考线、截面两个方面着手构建弓形激波的外场网格，有利于保障得到的外场网格的网格准确度和网格精度。并且，生成外场网格的过程无需过多的人工干预，有利于降低生成外场网格的门槛、降低人工成本。本发明中的方法可以应用于与飞机配合或者装到飞机上的设备、飞行衣、降落伞、动力装置或推进传动装置在飞机中的配置或安装。本发明中的方法应用于一般的控制或调价系统、这种系统的功能单元、用于这种系统或单元的监视或测试装置。本发明中的方法还可以应用于宇宙航行，及其所用的飞行器或设备。
本申请所要达到的技术效果通过以下方案实现：
第一方面，本发明提供一种弓形激波外场网格的生成方法，所述方法包括：
获取模型网格文件，其中，所述模型网格文件表示出目标模型在运动过程中能够在周围产生弓形激波，所述目标模型的表面包含初始拉伸网格面和表面网格；所述弓形激波在垂直于运动方向上的截面为参考面；
根据所述目标模型运动过程中处于所述运动方向的最前端的网格单元包含在不同方向上的节点数量，确定所述目标模型的模态，其中，所述模态是全模、半模中的一种；
根据所述目标模型的模态、以及预设的构造参数，识别出表征初始弓形激波的表面的网格单元，其中，构造参数是用于表征初始弓形激波的至少部分尺寸的参数；
对所述初始拉伸网格面的网格单元上的至少部分节点进行差值处理，得到参考线；
确定所述初始弓形激波的表面的网格单元的节点，相较于所述参考线的比例值，其中，所述比例值表示出所述节点和与之相邻的节点之间的距离，与所述节点在所述初始弓形激波的表面的网格单元的位置的相对关系；
根据所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线上点，作为第一目标点；
对所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线与所述目标模型的表面之间区域的点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点；
确定所述初始弓形激波的表面的网格单元的节点与所述最前端构成的向量与所述截面的夹角；
根据所述夹角，确定所述初始弓形激波的表面的网格单元的节点在所述截面上对应的点，作为第三目标点；
根据所述第一目标点、所述第二目标点、以及所述第三目标点，构建弓形激波的外场网格。
在本发明一个可选的实施例中，根据所述目标模型运动过程中处于所述运动方向的最前端的网格单元包含在不同方向上的节点数量，确定所述目标模型的模态，包括：
从所述表面网格的各网格单元中，确定出距所述最前端距离最小、边长为指定尺寸的网格单元阵列；
将所述网格单元阵列在第一方向上的节点数量确定为第一数量，并且，将所述网格单元阵列在第二方向上的节点数量确定为第二数量，其中，所述第一方向、第二方向与所述运动方向两两垂直；
若所述第一数量等于所述第二数量，则所述目标模型是全模；若所述第一数量不等于所述第二数量，则所述目标模型是半模。
在本发明一个可选的实施例中，所述构造参数包括：球头半径、第一距离HeadDis、第二距离TailDis、初始弓形激波的侧面高度、以及所述截面的短轴，所述第一距离是所述最前端距所述初始弓形激波与所述运动方向所在轴线的交点之间的距离，所述第二距离是所述初始弓形激波的尾部距所述目标模型的表面的距离，其中，根据所述目标模型的模态、以及预设的构造参数，识别出表征初始弓形激波的表面的网格单元之前，所述方法还包括以下至少一项：
对所述网格单元阵列中的节点进行球形拟合，将所述球形拟合的结果表示出的球形半径，作为所述球头半径；
将所述球头半径的第一系数倍，作为所述第一距离；其中，所述第一系数是预设值；
将所述目标模型在所述运动方向上的尺寸，作为第一尺寸；将所述第一尺寸的第二系数倍，作为所述第二距离；其中，所述第二系数是预设值；
将所述第一尺寸的第三系数倍，作为所述初始弓形激波的侧面高度；其中，所述第三系数是预设值；
所述侧面高度与所述截面的短轴相等或不相等。
在本发明一个可选的实施例中，根据所述目标模型的模态、以及预设的构造参数，识别出表征初始弓形激波的表面的网格单元，包括：
根据所述目标模型的模态、以及预设的构造参数，识别出待处理弓形激波；
根据网格单元属性和连接性，从所述待处理弓形激波中去除缝隙填充面和/或所述待处理弓形激波的尾部网格面，得到初始弓形激波。
在本发明一个可选的实施例中，确定所述初始弓形激波的表面的网格单元的节点，相较于所述参考线的比例值ratio，包括：
将所述初始弓形激波上任意一个节点，确定为目标节点pt；
确定所述目标节点的第一向量
Figure 700058DEST_PATH_IMAGE001
、第二向量
Figure 745374DEST_PATH_IMAGE002
、以及第三向量
Figure 202900DEST_PATH_IMAGE003
；
Figure 751693DEST_PATH_IMAGE004
式中，
Figure 839735DEST_PATH_IMAGE005
是所述参考线在i方向上节点，i是所述初始弓形激波的网格单元的排列方向之一；j是所述初始弓形激波的网格单元除i以外的排列方向之一；
根据所述第一向量
Figure 919818DEST_PATH_IMAGE006
、第二向量
Figure 853138DEST_PATH_IMAGE007
、以及第三向量
Figure 584334DEST_PATH_IMAGE008
确定所述目标节点的比例值ratio；
Figure 843277DEST_PATH_IMAGE009
式中，
Figure 332027DEST_PATH_IMAGE010
，
Figure 147668DEST_PATH_IMAGE011
是所述目标节点
Figure 671053DEST_PATH_IMAGE012
在所述第一方向上的坐标，
Figure 835318DEST_PATH_IMAGE013
是所述目标节点
Figure 139260DEST_PATH_IMAGE012
在所述第二方向上的坐标。
在本发明一个可选的实施例中，根据所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线上点，作为第一目标点，包括：
确定所述初始弓形激波的节点中离散点的占比s；
Figure 679963DEST_PATH_IMAGE014
针对所述目标节点，若所述目标节点的比例值与所述占比相等，则将所述目标节点确定为第一目标点。
在本发明一个可选的实施例中，对所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线与所述目标模型的表面之间区域的点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点，包括：
针对所述目标节点，若所述目标节点的比例值小于所述占比，则将所述目标节点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点。
第二方面，本发明提供一种弓形激波外场网格的生成装置，所述装置用于实现第一方面之任一所述方法。
第三方面，本申请实施例还提供一种电子设备，包括：
处理器；以及
被安排成存储计算机可执行指令的存储器，所述可执行指令在被执行时使所述处理器执行所述第一方面中的方法。
第四方面，本申请实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储一个或多个程序，所述一个或多个程序当被包括多个应用程序的电子设备执行时，使得所述电子设备执行所述第一方面中的方法。
技术效果描述：交底书记载根据导入网格可自动区分半模和全模，生成不同的外场网格；可精确控制弓形激波外场的形状，实时调整；
根据导入网格自动构造激波面，用户可更改关键形状控制参数。
附图说明
为了更清楚地说明本申请实施例或现有的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本申请一实施例中一种弓形激波外场网格的生成方法流程示意图；
图2为本申请一实施例中目标模型及针对目标模型的表面进行拉伸得到初始弓形激波示意图；
图3为本申请一实施例中涉及的网格单元阵列示意图；
图4为本申请一实施例中涉及的参考线示意图；
图5为本申请一实施例中得到的弓形激波在一视角下的示意图；
图6为本申请一实施例中得到的弓形激波在另一视角下的示意图；
图7为本申请一实施例中一种弓形激波外场网格的生成装置结构示意图；
图8为本申请一实施例中一个实施例电子设备的结构示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合具体实施例及相应的附图对本申请的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
下面通过具体实施方式结合附图对本发明作进一步详细说明。其中不同实施方式中类似元件采用了相关联的类似的元件标号。在以下的实施方式中，很多细节描述是为了使得本申请能被更好的理解。然而，本领域技术人员可以毫不费力的认识到，其中部分特征在不同情况下是可以省略的，或者可以由其它元件、材料、方法所替代。在某些情况下，本申请相关的一些操作并没有在说明书中显示或者描述，这是为了避免本申请的核心部分被过多的描述所淹没，而对于本领域技术人员而言，详细描述这些相关操作并不是必要的，他们根据说明书中的描述以及本领域的一般技术知识即可完整了解相关操作。
另外，说明书中所描述的特点、操作或者特征可以以任意适当的方式结合形成各种实施方式。同时，方法描述中的各步骤或者动作也可以按照本领域技术人员所能显而易见的方式进行顺序调换或调整。因此，说明书和附图中的各种顺序只是为了清楚描述某一个实施例，并不意味着是必须的顺序，除非另有说明其中某个顺序是必须遵循的。
本文中为部件所编序号本身，例如“第一”、“第二”等，仅用于区分所描述的对象，不具有任何顺序或技术含义。而本申请所说“连接”、“联接”，如无特别说明，均包括直接和间接连接(联接)。
以下结合附图，详细说明本申请各实施例提供的技术方案。
本发明中的弓形激波外场网格的生成方法包括以下步骤：
S100：获取模型网格文件。
本发明中的模型网格文件是针对目标体生成的文件，可选地为数模绘制网格。本发明对模型网格文件的格式不做具体限制，在本发明一个可选的实施例中，模型网格文件是网格的CAD文件。本发明涉及的目标体包括但不限于飞行器。
本发明中的目标模型是针对目标物形成的模型。在可选的实施例中，目标模型是以文件的形式进行表征的，在某些情况下，目标模型可以具有可视化的图形结构。所述模型网格文件表示出目标模型在运动过程中能够在周围产生弓形激波，所述目标模型的表面包含初始拉伸网格面strechDoms和表面网格surDoms，初始拉伸网格面strechDoms和表面网格surDoms是通过网格属性信息进行识别。所述弓形激波在垂直于所述运动方向上的截面为参考面。
在本发明一个可选的实施例中，本发明中的网格单元是四边形；在其他可选的实施例中，网格单元还可以是三角形。为便于说明，在下文的各示例中，以网格单元均为四边形为例进行说明。可选地，本发明中涉及的截面是圆形或椭圆形。为便于说明，在下文的各示例中，以截面是椭圆形为例进行说明。
模型网格文件表示出的目标模型示例性地如图2所示的包含了目标模型和弓形激波的图形。在图2所示的图形内容中，surfDoms部分在示意图的DPE弧形区域内，由于它不参与拉伸计算，故在图中未画出。初始的StrechDoms是包括DPE弧线的面、DCE直线的面，如果是半模，还包括PC上的面。其中，A与B点对称，特征长度为L，激波截面为椭圆，球头半径为r，P点到O点距离大于0.1r，尾部大于3L。
可选地，本步骤的过程具体可以是导入模型网格文件，如果目标模型的网格不在坐标系上，则调整坐标系（例如笛卡尔坐标系），目标模型的网格调整到指定轴上，下文均假设为x轴。判断目标模型的网格是否满足弓形激波外场生成条件，若满足以下条件生成失败，需要重新生成模型网格文件：没有表面网格、没有空间网格、球头拟合失败。
S102：根据所述目标模型运动过程中处于所述运动方向的最前端的网格单元包含在不同方向上的节点数量，确定所述目标模型的模态。
本发明中涉及的模态是指全模和半模。模态即为目标模型的状态，是受到模型的建模过程的影响的结果。本步骤即为本案的执行主体通过自动的方式判断目标模型的模态是全模还是半模的过程。该过程无需人工干预，能够自动地实现，有利于降低人工成本。
在本发明一个可选的实施例中，确定目标模型的模态的过程可以是：从所述表面网格的各网格单元中，确定出距所述最前端距离最小、边长为指定尺寸的网格单元阵列。将所述网格单元阵列在第一方向上的节点数量确定为第一数量，并且，将所述网格单元阵列在第二方向上的节点数量确定为第二数量，其中，所述第一方向、第二方向与所述运动方向两两垂直。若所述第一数量等于所述第二数量，则所述目标模型是全模；若所述第一数量不等于所述第二数量，则所述目标模型是半模。在后续步骤中，针对模态为半模的目标模型，则将对称面上的网格去掉，不拉伸外场。
可选地，本发明中的目标模型是建立在笛卡尔坐标系中的，运动方向与笛卡尔坐标系的x轴平行，且目标模型的运动轨迹与x轴至少部分地重合。下文以该笛卡尔坐标系所示的方位对本发明中的方法进行说明，则第一方向即为y轴指示的方向、第二方向即为z轴指示的方向。
其中，指定尺寸是以网格单元一边的长度为步长的指定数量（例如，第一数量、第二数量）倍得到的。指定数量是与网格单元的尺寸相关的量（例如，负相关）。在指定数量等于6的情况下，网格单元阵列示例性地如图3所示。在图3所示的示例中，第一数量和第二数量均等于6。指定数量可以由人工设定。
具体地，根据网格属性信息分类出表面网格surDoms和初始的外场拉伸网格strechDoms。计算表面网格大小GridSize。获取表面网格中X值最小单元，找到周围6*6区域内的顶点SpherePts并记录X最小单元在Y方向和Z方向的6个区域的节点数目。采用SpherePts数据进行球拟合，得到球头半径SphereR。根据节点数目数模区分是半模还是全模，如果相等则为全模，不相等则为半模，后续假设半模情况为Y轴，全模也适应。
S104：根据所述目标模型的模态、以及预设的构造参数，识别出表征初始弓形激波的表面的网格单元。
本发明中的构造参数可以是在本步骤之前确定出的，也可以是在本步骤确定出。在一个可选的实施例中，构造参数可以包括：球头半径、第一距离（HeadDis）、第二距离（TailDis）、初始弓形激波的侧面高度、以及所述截面的短轴。其中，第一距离是所述最前端距所述初始弓形激波与所述运动方向所在轴线的交点之间的距离，所述第二距离是所述初始弓形激波的尾部距所述目标模型的表面的距离。在其他可选的实施例中，还可以采用其他构造参数，在此不一一列举。构造参数是用于表征初始弓形激波的至少部分尺寸的参数。
外场网格实际是对输入的目标模型的网格拓扑的放大，但输入如果有凸起的部分变为抛物线的适应位置边长形状。可输入构造参数的不同的参数值，改变外场的形状。
在一个可选的实施例中，识别出初始弓形激波的表面的网格单元的过程可以是：根据所述目标模型的模态、以及预设的构造参数，生成待处理弓形激波；根据网格单元属性和连接性，从所述待处理弓形激波中去除缝隙填充面（缝隙填充的块的网面关联两个块）和/或所述待处理弓形激波的尾部网格面，得到初始弓形激波。
本发明中的过程是基于参考线和截面构建外场网格的过程，下面将对参考线、截面进行分别地说明。
具体地，根据尾部面在一个近似平面的属性，先根据最大X的网格面，自动调用角度选，去掉strechDoms中的尾部网格面。根据网格块属性和连接性去掉缝隙填充的面（缝隙填充的块的网格面关联两个块）。计算空间网格大小SpaceSize。计算默认椭球参数。根据参数构造抛物线方程和椭圆方程。
从StrechDoms中识别去掉尾部的网格面。然后从StrechDoms中去掉缝隙填充的网格面。然后，采用以下计算空间参数值大小。
SpaceX = maxX-minX
式中，
Figure 859185DEST_PATH_IMAGE015
，
Figure 459931DEST_PATH_IMAGE016
，
Figure 923273DEST_PATH_IMAGE017
，X的取值为 StrechDoms[i][j][k]的x坐标值，Y、Z的取值类似。k表示第几个外场网格面，i表示此网格面的i向索引，j表示此网格面的j向索引。GridSize是表面网格大小。
之后，采用以下计算
Figure 329984DEST_PATH_IMAGE018
。根据最小X值确定弓形激波顶点，弓形激波尾点取为最大X值。
Figure 827961DEST_PATH_IMAGE019
S106：对所述初始拉伸网格面的网格单元上的至少部分节点进行差值处理，得到参考线。
由于，从侧面看，弓形激波的面是一个二维的抛物线，则本发明中的参考线（refCn）是二维的线条、是抛物线，以PTi表征参考线上的节点（可选地，参考线上的节点数量设定为101，参考线的长度为Len），参考线示例性地，如图4所示。参考线可以通过以下公式表征：
Figure 334029DEST_PATH_IMAGE020
具体地，获取StrechDoms上的点。采用插值方式计算对应X坐标的最远距离MaxDis的Y、Z值，为保证精度越靠近球头得位置，取点越密集，目的是保留球头。参考线的X球头方向顶点取为最小X点，保证映射可靠。根据这些有序的离散点构造参考线refCn。
构建外场网格的过程可以概括为：
S1：获取StrechDoms线的去重端点。
首先，计算当前点在二维参考线refCn上的比例的值。
然后，根据比例值计算到抛物线上的坐标。
之后，根据YZ方向偏移角度计算椭圆上的点。对称面法向上点的值进行校正，如果原来点在对称正平面，算出的点也在对称正平面，否则在对称负平面。计算这些点在弓形激波上的点。
S2：根据当前点和弓形激波上的点构造三次样条曲线，并设置分布点数。
S3：根据底层线和弓形激波上的点构造外框线，得到弓形激波的外形。
构建外场网格的过程具体如下：
S108：确定所述初始弓形激波的表面的网格单元的节点，相较于所述参考线的比例值。
本发明中的比例值表示出所述节点和与之相邻的节点之间的距离，与所述节点在所述初始弓形激波的表面的网格单元的位置的相对关系。
由于初始弓形激波的表面的网格单元的节点不唯一，则需要针对每个节点确定出其比例值。不同节点的比例值的取值有可能相同、也有可能不同。
S110：根据所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线上点，作为第一目标点。
本发明中的目标点是外场网格上的节点，除本步骤中的第一目标点以外，目标点还包括后续步骤中涉及的第二目标点和第三目标点。第一目标点是初始弓形激波的表面的网格单元上与外场网格共用的节点。
S112：对所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线与所述目标模型的表面之间区域的点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点。
第二目标点是通过对初始弓形激波的表面的网格单元上、且位于外场网格所属区间上的节点进行调整得到的。
第一目标点和第二目标点都是位于外场网格对应于弓形激波的侧壁上的节点。
S114：确定所述初始弓形激波的表面的网格单元的节点与所述最前端构成的向量与所述截面的夹角。
由于初始弓形激波的表面的网格单元的节点不唯一，则需对弓形激波的表面的网格单元的每个节点执行本步骤。在下文中，该夹角可以通过angle进行表示。
从正面看，激波面是一个椭圆（也可退化为圆），截面可以通过以下公式表征：
Figure 363296DEST_PATH_IMAGE021
S116：根据所述夹角，确定所述初始弓形激波的表面的网格单元的节点在所述截面上对应的点，作为第三目标点。
S118：根据所述第一目标点、所述第二目标点、以及所述第三目标点，构建弓形激波的外场网格。
可选地，可以对第一目标点、所述第二目标点、以及所述第三目标点构造三次样条曲线，并设置分布点数。然后，根据底层线和激波上的点构造外框线，我们就得到弓形激波的外形。最后装配面装配块得到外场网格。得到的弓形激波示例性地如图5和图6所示。
得到外场网格之后，即可执行基于外场网格的目标体气动分析。
通过本发明中的方法，在生成弓形激波的外场网格的过程中，无论目标模型的模态是全模还是半模，均能够高效地生成外场网格。并且，生成外场网格的过程无需过多的人工干预，有利于降低生成外场网格的门槛、降低人工成本。此外，通过本发明中的方法得到的外场网格其精度、准确性也较高。
现就本发明中的技术方案进行进一步地说明。
1、关于构造参数的确定。
（1）针对球头半径R：对所述网格单元阵列中的节点进行球形拟合，将所述球形拟合的结果表示出的球形半径，作为所述球头半径。具体地，计算表面网格大小GridSize。获取表面网格中X值最小单元，找到周围6*6区域内的顶点SpherePts并记录X最小单元在Y方向和Z方向的6个区域的点数目。可选地，可以通过以下公式（一）根据对称面上的靠近外场的网格面计算对称轴。
Figure 245801DEST_PATH_IMAGE022
公式（一）
式中，i表示此网格面的i向索引，j表示此网格面的j向索引。如果sumY>sumZ，则是Z向对称，否则Y向对称，记录对称信息。
Figure 660602DEST_PATH_IMAGE023
是索引为i、j的节点Pt在y轴上的坐标。
对6*6的网格单元阵列内的节点进行球拟合，得到球头半径R。后续假设半模情况为Y轴，全模也适应。
（2）针对第一距离
Figure 337571DEST_PATH_IMAGE024
：获取第一系数（示例性地，第一系数可以是经验值，例如，第一系数等于0.5）。将所述球头半径的第一系数倍，作为所述第一距离；其中，所述第一系数是预设值。可选地，可以通过以下公式（二）计算第一距离。
Figure 775506DEST_PATH_IMAGE025
公式（二）
（3）针对第二距离TailDis：获取第二系数（示例性地，第二系数可以是经验值，例如，第二系数等于4）。将所述目标模型在所述运动方向上的尺寸
Figure 274751DEST_PATH_IMAGE026
，作为第一尺寸；将所述第一尺寸的第二系数倍，作为所述第二距离；其中，所述第二系数是预设值。可选地，可以通过以下公式（三）计算第二距离。
Figure 481742DEST_PATH_IMAGE027
公式（三）
（4）针对初始弓形激波的侧面高度SideDis：获取第三系数（示例性地，第三系数可以是经验值，例如，第三系数等于0.5）。将所述第一尺寸的第三系数倍，作为所述初始弓形激波的侧面高度。可选地，可以通过以下公式（四）计算初始弓形激波的侧面高度。
Figure 657508DEST_PATH_IMAGE028
公式（四）
（5）针对截面的短轴EllipseR：所述侧面高度与所述截面的短轴相等或不相等。可选地，可以通过以下公式（五）计算截面的短轴。
Figure 582739DEST_PATH_IMAGE029
公式（五）
2、关于计算比例值ratio。
将所述初始弓形激波上任意一个节点，确定为目标节点
Figure 807047DEST_PATH_IMAGE030
，坐标为（x,y,z）。先根据对称面确定和坐标值确定是上半激波面还是下半激波面（如：关于Z轴对称，y大于0为上半激波面，小于等于0为下半激波面，后续所有的假设都是基于Z轴对称的上半激波面展开，其他三种情况类似计算方法）。计算对称方向上的基线值的公式如下：
Figure 947172DEST_PATH_IMAGE031
然后，确定所述目标节点的第一向量
Figure 231523DEST_PATH_IMAGE001
、第二向量
Figure 440788DEST_PATH_IMAGE002
、以及第三向量
Figure 468786DEST_PATH_IMAGE003
。各个向量可以通过以下公式（六）计算得到。Pt1,Pt2…… Ptn为二维参考线上的离散点，N为最大的离散点数索引，startI为比例值所在的参考线段的离散点起始坐标索引值。节点可以通过其在笛卡尔坐标系下的坐标进行表征，则公式（六）中的计算可以是基于坐标进行的计算。
Figure 650369DEST_PATH_IMAGE032
公式（六）
式中，
Figure 918671DEST_PATH_IMAGE033
是所述参考线在i方向上节点，i是所述初始弓形激波的网格单元的排列方向之一；j是所述初始弓形激波的网格单元除i以外的排列方向之一；
之后，根据所述第一向量
Figure 818493DEST_PATH_IMAGE034
、第二向量
Figure 446921DEST_PATH_IMAGE035
、以及第三向量
Figure 483010DEST_PATH_IMAGE036
确定所述目标节点的比例值ratio。可选地，可以通过以下公式（七）计算得到比例值。初始弓形激波上的节点不唯一，则可以针对每个节点分别计算其比例值。
Figure 109163DEST_PATH_IMAGE009
公式（七）
式中，
Figure 43752DEST_PATH_IMAGE010
，
Figure 413554DEST_PATH_IMAGE011
是所述目标节点
Figure 569729DEST_PATH_IMAGE012
在所述第一方向上的坐标，
Figure 429100DEST_PATH_IMAGE013
是所述目标节点
Figure 37936DEST_PATH_IMAGE012
在所述第二方向上的坐标。
确定比例值的过程的伪代码示例性地，如下：
For i From 0 to N-1 step 1
StartI ←i-1
If(pt和pti的距离等于0) then
Ratio ← 1
结束循环
End if
v1←pti-pti-1
v2←pt-pti-1
If(v1*v2<0) then
Ratio ← 0
继续循序
End if
v3←pt-pti
If(v1*v3>0) then
startI←I;
继续查找
Else
Ratio ← v1*v2/(v1*v1)
结束循序
End if
End for
If(startI == N-1) then
Ratio ←1
sumLen←0
for i From 1 to startI step 1
sumLen ←sumLen+pti-1和pti的距离
End for
SumLen = ratio *ptstartI-1到ptstartI的距离
Ratio←sumLen/Len
3、关于确定第一目标点。
首先，如果ratio的取值是0或1，则目标节点为首点或尾点，如果是首点
Figure 758899DEST_PATH_IMAGE037
值为0，坐标为（minX+
Figure 769580DEST_PATH_IMAGE038
,0,0）如果是尾点
Figure 737536DEST_PATH_IMAGE039
,坐标值为(maxX,
Figure 895985DEST_PATH_IMAGE040
,0),然后结束。如果ratio的取值不是0或1，则遍历计算当前离散点的占比s。可选地，可以通过以下公式（八）计算占比s。
Figure 607589DEST_PATH_IMAGE041
公式（八）
然后，针对所述目标节点，若所述目标节点的比例值与所述占比相等，当前
Figure 738356DEST_PATH_IMAGE042
点为抛物线上坐标，则将所述目标节点确定为第一目标点。
若所述目标节点的比例值大于所述占比，则重新确定目标节点，按b更新占比s。
4、关于确定第二目标点。
针对所述目标节点，若所述目标节点的比例值小于所述占比，则将所述目标节点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点。
5、关于确定第三目标点。
本发明一个可选的实施例中，根据y、z方向偏移角度计算截面上的节点。
首先，确定目标向量v1为（0,0,1），当前目标节点的向量v2为
Figure 690263DEST_PATH_IMAGE043
。然后，计算v1与v2的夹角angle。之后，通过以下公式（九）计算椭圆的长轴和短轴。
Figure 8112DEST_PATH_IMAGE044
公式（九）
之后，通过以下公式（十），根据长轴和短轴及角度得到椭圆界面上的坐标点。
Figure 523407DEST_PATH_IMAGE045
公式（十）
在本发明一个可选的实施例中，此后，对称面法向上点的值进行校正，如果原来点在对称正平面，算出的点也在对称正平面，否则在对称负平面。具体地，如果目标节点是正平面的，就将目标节点设为正平面的值，如关于y轴对称的网格，目标节点Z值大于0，则算出的值Z值校正为正值。
基于同样的思路，本发明实施例还提供了对应于图1所示部分过程的一种弓形激波外场网格的生成装置。
如图7所示，所述一种弓形激波外场网格的生成装置，该装置可以包括以下单元中的一个或多个：
文件获取模块700，配置为：获取模型网格文件，其中，所述模型网格文件表示出目标模型在运动过程中能够在周围产生弓形激波，所述目标模型的表面包含初始拉伸网格面和表面网格；所述弓形激波在垂直于所述运动方向上的截面为参考面。
模态判断模块702，配置为：根据所述目标模型运动过程中处于所述运动方向的最前端的网格单元包含在不同方向上的节点数量，确定所述目标模型的模态，其中，所述模态是全模、半模中的一种。
初始弓形激波生成模块704，配置为：根据所述目标模型的模态、以及预设的构造参数，识别出表征初始弓形激波的表面的网格单元，其中，构造参数是用于表征初始弓形激波的至少部分尺寸的参数。
参考线生成模块706，配置为：对所述初始拉伸网格面的网格单元上的至少部分节点进行差值处理，得到参考线。
比例值确定模块708，配置为：确定所述初始弓形激波的表面的网格单元的节点，相较于所述参考线的比例值，其中，所述比例值表示出所述节点和与之相邻的节点之间的距离，与所述节点在所述初始弓形激波的表面的网格单元的位置的相对关系。
第一目标点确定模块710，配置为：根据所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线上点，作为第一目标点。
第二目标点确定模块712，配置为：对所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线与所述目标模型的表面之间区域的点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点。
夹角确定模块714，配置为：确定所述初始弓形激波的表面的网格单元的节点与所述最前端构成的向量与所述截面的夹角。
第三目标点确定模块716，配置为：根据所述夹角，确定所述初始弓形激波的表面的网格单元的节点在所述截面上对应的点，作为第三目标点。
构建模块718，配置为：根据所述第一目标点、所述第二目标点、以及所述第三目标点，构建弓形激波的外场网格。
在本发明一个可选的实施例中，模态判断模块702，具体配置为：从所述表面网格的各网格单元中，确定出距所述最前端距离最小、边长为指定尺寸的网格单元阵列；将所述网格单元阵列在第一方向上的节点数量确定为第一数量，并且，将所述网格单元阵列在第二方向上的节点数量确定为第二数量，其中，所述第一方向、第二方向与所述运动方向两两垂直；若所述第一数量等于所述第二数量，则所述目标模型是全模；若所述第一数量不等于所述第二数量，则所述目标模型是半模。
在本发明一个可选的实施例中，构造参数包括：球头半径、第一距离HeadDis、第二距离TailDis、初始弓形激波的侧面高度、以及所述截面的短轴，所述第一距离是所述最前端距所述初始弓形激波与所述运动方向所在轴线的交点之间的距离，所述第二距离是所述初始弓形激波的尾部距所述目标模型的表面的距离。
在本发明一个可选的实施例中，所述装置还包括参数获取模块，配置为：对所述网格单元阵列中的节点进行球形拟合，将所述球形拟合的结果表示出的球形半径，作为所述球头半径。
在本发明一个可选的实施例中，参数获取模块还配置为：将所述球头半径的第一系数倍，作为所述第一距离；其中，所述第一系数是预设值。
在本发明一个可选的实施例中，参数获取模块还配置为：将所述目标模型在所述运动方向上的尺寸，作为第一尺寸；将所述第一尺寸的第二系数倍，作为所述第二距离；其中，所述第二系数是预设值。
在本发明一个可选的实施例中，参数获取模块还配置为：将所述第一尺寸的第三系数倍，作为所述初始弓形激波的侧面高度；其中，所述第三系数是预设值。
在本发明一个可选的实施例中，所述侧面高度与所述截面的短轴相等。
在本发明一个可选的实施例中，初始弓形激波生成模块704具体配置为：根据所述目标模型的模态、以及预设的构造参数，生成待处理弓形激波；根据网格单元属性和连接性，从所述待处理弓形激波中去除缝隙填充面和/或所述待处理弓形激波的尾部网格面，得到初始弓形激波。
在本发明一个可选的实施例中，比例值确定模块708具体配置为：将所述初始弓形激波上任意一个节点，确定为目标节点
Figure 570997DEST_PATH_IMAGE046
；
确定所述目标节点的第一向量
Figure 880756DEST_PATH_IMAGE006
、第二向量
Figure 685901DEST_PATH_IMAGE007
、以及第三向量
Figure 552356DEST_PATH_IMAGE008
；
Figure 392137DEST_PATH_IMAGE032
式中，
Figure 872796DEST_PATH_IMAGE047
是所述参考线在i方向上节点，i是所述初始弓形激波的网格单元的排列方向之一；j是所述初始弓形激波的网格单元除i以外的排列方向之一；
根据所述第一向量
Figure 493134DEST_PATH_IMAGE006
、第二向量
Figure 84652DEST_PATH_IMAGE007
、以及第三向量
Figure 44518DEST_PATH_IMAGE008
确定所述目标节点的比例值ratio；
Figure 509128DEST_PATH_IMAGE009
式中，
Figure 288865DEST_PATH_IMAGE010
，
Figure 684075DEST_PATH_IMAGE011
是所述目标节点
Figure 560764DEST_PATH_IMAGE012
在所述第一方向上的坐标，
Figure 383226DEST_PATH_IMAGE013
是所述目标节点
Figure 650260DEST_PATH_IMAGE012
在所述第二方向上的坐标。
在本发明一个可选的实施例中，第一目标点确定模块710具体配置为：确定所述初始弓形激波的节点中离散点的占比s；
Figure 662209DEST_PATH_IMAGE014
针对所述目标节点，若所述目标节点的比例值与所述占比相等，则将所述目标节点确定为第一目标点。
在本发明一个可选的实施例中，第二目标点确定模块712具体配置为：针对所述目标节点，若所述目标节点的比例值小于所述占比，则将所述目标节点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点。
图8是本申请的一个实施例电子设备的结构示意图。请参考图8，在硬件层面，该电子设备包括处理器，可选地还包括内部总线、网络接口、存储器。其中，存储器可能包含内存，例如高速随机存取存储器(Random-Access Memory，RAM)，也可能还包括非易失性存储器（non-volatile memory），例如至少1个磁盘存储器等。当然，该电子设备还可能包括其他业务所需要的硬件。
处理器、网络接口和存储器可以通过内部总线相互连接，该内部总线可以是ISA(Industry Standard Architecture，工业标准体系结构）总线、PCI(PeripheralComponent Interconnect，外设部件互连标准)总线或EISA(Extended Industry StandardArchitecture，扩展工业标准结构）总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图8中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
存储器，用于存放程序。具体地，程序可以包括程序代码，所述程序代码包括计算机操作指令。存储器可以包括内存和非易失性存储器，并向处理器提供指令和数据。
处理器从非易失性存储器中读取对应的计算机程序到内存中然后运行，在逻辑层面上形成一种弓形激波外场网格的生成方法。处理器，执行存储器所存放的程序，并具体用于执行前述任意一种弓形激波外场网格的生成方法。
上述如本申请图1所示实施例揭示的一种弓形激波外场网格的生成方法可以应用于处理器（即，本发明中的删除控制模块）中，或者由处理器实现。处理器可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器，包括中央处理器（Central Processing Unit，CPU）、网络处理器（Network Processor，NP）等；还可以是数字信号处理器（Digital Signal Processor，DSP）、专用集成电路（ApplicationSpecific Integrated Circuit，ASIC）、现场可编程门阵列（Field－Programmable GateArray，FPGA）或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。
该电子设备还可执行图1中一种弓形激波外场网格的生成方法，并实现图1所示实施例的功能，本申请实施例在此不再赘述。
本申请实施例还提出了一种计算机可读存储介质，该计算机可读存储介质存储一个或多个程序，该一个或多个程序包括指令，该指令当被包括多个应用程序的电子设备执行时，能够使该电子设备执行图1所示实施例中一种弓形激波外场网格的生成方法执行的方法，并具体用于执行前述的任意一种弓形激波外场网格的生成方法。
本领域内的技术人员应明白，本申请的实施例可提供为方法、方法、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备（方法）、和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器 (CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器 (RAM) 和/或非易失性内存等形式，如只读存储器 (ROM) 或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存 (PRAM)、静态随机存取存储器 (SRAM)、动态随机存取存储器 (DRAM)、其他类型的随机存取存储器 (RAM)、只读存储器 (ROM)、电可删除可编程只读存储器 (EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘 (DVD) 或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体 (transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、方法或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (10)

1.一种弓形激波外场网格的生成方法，其特征在于，所述方法包括：
获取模型网格文件，其中，所述模型网格文件表示出目标模型在运动过程中能够在周围产生弓形激波，所述目标模型的表面包含初始拉伸网格面和表面网格；所述弓形激波在垂直于运动方向上的截面为参考面；以及
根据所述目标模型运动过程中处于所述运动方向的最前端的网格单元包含在不同方向上的节点数量，确定所述目标模型的模态，其中，所述模态是全模、半模中的一种；
根据所述目标模型的模态、以及预设的构造参数，识别出表征初始弓形激波的表面的网格单元，其中，构造参数是用于表征初始弓形激波的至少部分尺寸的参数；
对所述初始拉伸网格面的网格单元上的至少部分节点进行差值处理，得到参考线；
确定所述初始弓形激波的表面的网格单元的节点，相较于所述参考线的比例值，其中，所述比例值表示出所述节点和与之相邻的节点之间的距离，与所述节点在所述初始弓形激波的表面的网格单元的位置的相对关系；
根据所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线上点，作为第一目标点；
对所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线与所述目标模型的表面之间区域的点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点；
确定所述初始弓形激波的表面的网格单元的节点与所述最前端构成的向量与所述截面的夹角；
根据所述夹角，确定所述初始弓形激波的表面的网格单元的节点在所述截面上对应的点，作为第三目标点；
根据所述第一目标点、所述第二目标点、以及所述第三目标点，构建弓形激波的外场网格。
2.如权利要求1所述的方法，其特征在于，根据所述目标模型运动过程中处于所述运动方向的最前端的网格单元包含在不同方向上的节点数量，确定所述目标模型的模态，包括：
从所述表面网格的各网格单元中，确定出距所述最前端距离最小、边长为指定尺寸的网格单元阵列；
将所述网格单元阵列在第一方向上的节点数量确定为第一数量，并且，将所述网格单元阵列在第二方向上的节点数量确定为第二数量，其中，所述第一方向、第二方向与所述运动方向两两垂直；
若所述第一数量等于所述第二数量，则所述目标模型是全模；若所述第一数量不等于所述第二数量，则所述目标模型是半模。
3.如权利要求2所述的方法，其特征在于，所述构造参数包括：球头半径、第一距离、第二距离、初始弓形激波的侧面高度、以及所述截面的短轴，所述第一距离是所述最前端距所述初始弓形激波与所述运动方向所在轴线的交点之间的距离，所述第二距离是所述初始弓形激波的尾部距所述目标模型的表面的距离，其中，根据所述目标模型的模态、以及预设的构造参数，识别出表征初始弓形激波的表面的网格单元之前，所述方法还包括以下至少一项：
对所述网格单元阵列中的节点进行球形拟合，将所述球形拟合的结果表示出的球形半径，作为所述球头半径；
将所述球头半径的第一系数倍，作为所述第一距离；其中，所述第一系数是预设值；
将所述目标模型在所述运动方向上的尺寸，作为第一尺寸；将所述第一尺寸的第二系数倍，作为所述第二距离；其中，所述第二系数是预设值；
将所述第一尺寸的第三系数倍，作为所述初始弓形激波的侧面高度；其中，所述第三系数是预设值；
所述侧面高度与所述截面的短轴相等或不相等。
4.如权利要求1所述的方法，其特征在于，根据所述目标模型的模态、以及预设的构造参数，识别出表征初始弓形激波的表面的网格单元，包括：
根据所述目标模型的模态、以及预设的构造参数，生成待处理弓形激波；
根据网格单元属性和连接性，从所述待处理弓形激波中去除缝隙填充面和/或所述待处理弓形激波的尾部网格面，得到初始弓形激波。
5.如权利要求2所述的方法，其特征在于，确定所述初始弓形激波的表面的网格单元的节点，相较于所述参考线的比例值ratio，包括：
将所述初始弓形激波上任意一个节点，确定为目标节点pt；
确定所述目标节点的第一向量
Figure 130143DEST_PATH_IMAGE001
、第二向量
Figure 576167DEST_PATH_IMAGE002
、以及第三向量
Figure 769251DEST_PATH_IMAGE003
；
Figure 908109DEST_PATH_IMAGE004
式中，
Figure 491537DEST_PATH_IMAGE005
是所述参考线在i方向上节点，i是所述初始弓形激波的网格单元的排列方向之一；j是所述初始弓形激波的网格单元除i以外的排列方向之一；
根据所述第一向量
Figure 554302DEST_PATH_IMAGE001
、第二向量
Figure 539575DEST_PATH_IMAGE002
、以及第三向量
Figure 849334DEST_PATH_IMAGE003
确定所述目标节点的比例值ratio；
Figure 982375DEST_PATH_IMAGE006
式中，
Figure 35782DEST_PATH_IMAGE007
，
Figure 875562DEST_PATH_IMAGE008
是所述目标节点
Figure 169271DEST_PATH_IMAGE009
在所述第一方向上的坐标，
Figure 461712DEST_PATH_IMAGE010
是所述目标节点
Figure 318810DEST_PATH_IMAGE009
在所述第二方向上的坐标。
6.如权利要求5所述的方法，其特征在于，根据所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线上点，作为第一目标点，包括：
确定所述初始弓形激波的节点中离散点的占比s；
Figure 75413DEST_PATH_IMAGE011
针对所述目标节点，若所述目标节点的比例值与所述占比相等，则将所述目标节点确定为第一目标点。
7.如权利要求6所述的方法，其特征在于，对所述比例值确定出所述初始弓形激波的表面的网格单元的节点中位于所述参考线与所述目标模型的表面之间区域的点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点，包括：
针对所述目标节点，若所述目标节点的比例值小于所述占比，则将所述目标节点进行二分法迭代计算，得到其位于所述参考线上点，作为第二目标点。
8.一种弓形激波外场网格的生成装置，其特征在于，所述装置用于实现权利要求1~7之任一所述方法。
9.一种电子设备，包括：
处理器；以及
被安排成存储计算机可执行指令的存储器，所述可执行指令在被执行时使所述处理器执行所述权利要求1~7之任一所述方法。
10.一种计算机可读存储介质，所述计算机可读存储介质存储一个或多个程序，所述一个或多个程序当被包括多个应用程序的电子设备执行时，使得所述电子设备执行所述权利要求1~7之任一所述方法。
CN202211189574.9A 2022-09-28 2022-09-28 一种弓形激波外场网格的生成方法 Active CN115270378B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211189574.9A CN115270378B (zh) 2022-09-28 2022-09-28 一种弓形激波外场网格的生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211189574.9A CN115270378B (zh) 2022-09-28 2022-09-28 一种弓形激波外场网格的生成方法

Publications (2)

Publication Number Publication Date
CN115270378A true CN115270378A (zh) 2022-11-01
CN115270378B CN115270378B (zh) 2022-12-30

Family

ID=83757678

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211189574.9A Active CN115270378B (zh) 2022-09-28 2022-09-28 一种弓形激波外场网格的生成方法

Country Status (1)

Country Link
CN (1) CN115270378B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115587552A (zh) * 2022-12-12 2023-01-10 中国空气动力研究与发展中心计算空气动力研究所 网格优化方法、装置、终端设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008186133A (ja) * 2007-01-29 2008-08-14 Hokkaido Univ 凸多面体に対する重心座標を用いた自由形状変形システム
US20130328870A1 (en) * 2009-07-31 2013-12-12 Desmond T. Grenfell System and Method for Generating 2D Texture Coordinates for 3D Meshed Surfaces
CN113034694A (zh) * 2021-03-26 2021-06-25 浙江大学 自适应于流场特征的度量正交网格自动生成方法及装置
WO2021142843A1 (zh) * 2020-01-19 2021-07-22 Oppo广东移动通信有限公司 图像扫描方法及装置、设备、存储介质
CN114357906A (zh) * 2021-12-31 2022-04-15 浙江大学 激波特征曲面网格重建方法及装置、电子设备、存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008186133A (ja) * 2007-01-29 2008-08-14 Hokkaido Univ 凸多面体に対する重心座標を用いた自由形状変形システム
US20130328870A1 (en) * 2009-07-31 2013-12-12 Desmond T. Grenfell System and Method for Generating 2D Texture Coordinates for 3D Meshed Surfaces
WO2021142843A1 (zh) * 2020-01-19 2021-07-22 Oppo广东移动通信有限公司 图像扫描方法及装置、设备、存储介质
CN113034694A (zh) * 2021-03-26 2021-06-25 浙江大学 自适应于流场特征的度量正交网格自动生成方法及装置
CN114357906A (zh) * 2021-12-31 2022-04-15 浙江大学 激波特征曲面网格重建方法及装置、电子设备、存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张来平等: "矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用", 《空气动力学学报》 *
王勇达等: "法向网格增长比对双椭球气动热数值模拟的影响", 《西南科技大学学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115587552A (zh) * 2022-12-12 2023-01-10 中国空气动力研究与发展中心计算空气动力研究所 网格优化方法、装置、终端设备及存储介质

Also Published As

Publication number Publication date
CN115270378B (zh) 2022-12-30

Similar Documents

Publication Publication Date Title
CN113505443B (zh) 一种任意外形的三维绕流问题自适应笛卡尔网格生成方法
CN115270378B (zh) 一种弓形激波外场网格的生成方法
CN109376403B (zh) 一种基于笛卡尔自适应重构技术的二维结冰数值模拟方法
CN112818573B (zh) 一种用于非结构网格的获取边界层非当地变量信息的方法
Wei et al. Fgr: Frustum-aware geometric reasoning for weakly supervised 3d vehicle detection
CN111310381A (zh) 一种三维水滴收集系数计算方法
CN114996658B (zh) 一种基于投影的高超声速飞行器气动热预测方法
CN112380793A (zh) 基于gpu的湍流燃烧数值模拟并行加速实现方法
Chen et al. GPU-accelerated incremental Euclidean distance transform for online motion planning of mobile robots
CN114492250A (zh) 基于递归分解的曲面网格生成方法及系统、计算机设备
Xia et al. Highly efficient wall-distance-based parallel unstructured overset grid assembly
CN117473655A (zh) 基于边坍缩网格优化的飞行器仿真驱动设计方法和装置
Karman et al. Mesh generation challenges: A commercial software perspective
Jung et al. Reverse engineering and database of off-the-shelf propellers for middle-size multirotors
Tian et al. Efficient Parallel Dynamic Overset Unstructured Grid for Flow Simulation with Moving Bodies
CN105808930A (zh) 基于服务器集群网络的预条件共轭梯度区域网平差方法及服务器集群网络
CN114638867A (zh) 基于特征提取模块和对偶四元数的点云配准方法及系统
CN115688618A (zh) 一种模型的表面网格光顺处理方法
CN118171605B (zh) 一种面向重叠网格的网格处理方法和装置
CN104699918A (zh) 一种计算舵偏角的方法及系统
CN117274534B (zh) 一种三维结构化网格的表面流线生成方法、设备及介质
CN116628863B (zh) 车辆风阻系数确定方法、设备和介质
CN110796729B (zh) 一种基于二叉树的网格划分方法
Damigos et al. Adjoint variable‐based shape optimization with bounding surface constraints
US20230101689A1 (en) Method, apparatus, and program for matching point cloud data

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