一种基于对称方向场的三维翼面四边形网格自动生成方法技术

本发明专利技术公开了一种基于对称方向场的三维翼面四边形网格自动生成方法，其特征在于，方法包括：导入三维翼型的截面形状参数，设定空间几何参数，计算生成控制点，并基于控制点获取B样条三维几何模型；针对获取的三维翼型B样条几何模型，通过采样点离散化生成三维翼型三角网格，并提取三角网格中的前缘、后缘和翼尖的离散化曲线特征；根据生成的三维翼型三角网格，计算曲线特征约束下的对称方向场，通过优化参数化梯度和方向场的误差得到三角网格的参数化映射，利用libQEx软件提取参数化坐标得到四边形网格。本申请实现了几何生成、网格剖分的自动化改进，大大减少了气动外形设计中的人工干预，起到了加速机翼气动外形设计流程的作用。作用。作用。

【技术实现步骤摘要】
一种基于对称方向场的三维翼面四边形网格自动生成方法


[0001]本专利技术涉及计算机辅助工程几何前处理
，具体涉及一种基于对称方向场的三维翼面四边形网格自动生成方法。

技术介绍

[0002]面向飞行器机翼气动外形设计的几何前处理技术，是一个涉及到计算机图形学和计算流体力学等多学科领域的典型多学科交叉问题。针对飞行器机翼气动外形设计，高效实现计算机辅助设计(Computer Aided Design，CAD)模型到离散化网格模型的转换，是后续开展计算流体力学(Computer Fluid Dynamics，CFD)偏微分方程组高效解算，以得到翼面气动特性的重要一环。
[0003]目前飞行器机翼气动外形设计的一大主流范式为：首先针对截面翼型进行气动设计，接着通过参数化拉伸获得三维翼型。在传统的计算机辅助工程(Computer Aided Engineering，CAE)流程中，针对该三维翼型的几何前处理需要借助网格划分软件来手动划分网格，整个几何前处理流程占整个CAE计算任务全部人力时间的60％
‑
80％，而其中网格划分又可以占到几何前处理流程耗时的80％。由于需要大量人工干预，因而自动化程度较低是当前CAE领域的一大技术瓶颈。

技术实现思路

[0004]本申请实施例针对上述情况，提出了一种基基于对称方向场的三维翼面四边形网格自动生成方法，以克服或部分克服现有技术的不足之处。本方法通过三维翼型的截面形状参数和空间几何参数生成B样条曲线模型，根据B样条曲线模型采样点离散化生成三角网格，并提取三角网格中的前缘、后缘和翼尖曲线特征，接着通过优化曲线特征约束下的方向场和参数化，提取四边形网格，实现从二维截面翼型到三维翼面再到高质量四边形网格的自动生成。
[0005]为了实现上述目的，本申请公开了一种基于对称方向场的三维翼面四边形网格自动生成方法，所述方法包括如下步骤：
[0006]步骤S1：导入三维翼型的截面形状参数，设定空间几何参数，计算生成控制点，并基于控制点获取B样条三维几何模型；
[0007]步骤S2：针对获取的三维翼型B样条几何模型，通过采样点离散化生成三维翼型三角网格，并提取所述三角网格中的前缘、后缘和翼尖的离散化曲线特征；
[0008]步骤S3：根据生成的三维翼型三角网格，计算曲线特征约束下的对称方向场，通过优化参数化梯度和方向场的误差得到三角网格的参数化映射，利用libQEx软件提取参数化坐标得到四边形网格。
[0009]可选地，所述步骤S1中，设定空间几何参数具体包括：设定包括翼根弦长、翼尖弦长、前缘后掠角、半翼展和后缘厚度在内的空间集合参数。
[0010]可选地，所述步骤S1中，基于控制点获取B样条三维几何模型的算法流程包括：
[0011]步骤S1
‑
1：导入NACA翼型描述格式的三维翼型截面形状参数文本文件，根据后缘厚度设置截面翼型上沿和下沿中的后缘点，接着对上下沿靠近后缘点线性插值，使得前缘到后缘钝面之间平滑过渡，生成具有后缘钝面的截面翼型；
[0012]步骤S1
‑
2：将三维翼型截面分为上下沿两部分，将其包含的翼型二维点作为控制点分别生成3阶B样条上沿曲线和下沿曲线；
[0013]步骤S1
‑
3：分别对上沿曲线和下沿曲线进行等数目均匀采样，使得上下沿采样点一一对应；
[0014]步骤S1
‑
4：根据翼根弦长、翼尖弦长、前缘后掠角和半翼展，计算翼尖相对翼根偏移；使用上一步骤S1
‑
3中的二维上下沿采样点创建三维翼型翼根上沿控制点、翼根下沿控制点、翼尖上沿控制点和翼尖下沿控制点，其中上下沿、翼根与翼尖的控制点具有一一对应关系；组合上沿控制点为上翼面控制点，组合下沿控制点为下翼面控制点，选择所有控制点中位于后缘钝面的四个点组成后缘钝面控制点；
[0015]步骤S1
‑
5：创建节点向量，设置上下翼面B样条曲面u方向为机翼弦长方向，v方向为翼展方向，阶数为3
×
1，创建上翼面和下翼面B样条曲面；设置后缘钝面B样条曲面u方向为机翼厚度方向，v方向为翼展方向，阶数为1
×
1，创建后缘钝面B样条曲面；
[0016]步骤S1
‑
6：根据上下翼面控制点计算平均水平方向向量和竖直方向向量，进而计算上下翼尖面控制点；
[0017]步骤S1
‑
7：设置上下翼尖面B样条曲面u方向为机翼弦长方向，v方向为机翼厚度方向，阶数为3
×
3，创建上下翼尖面B样条曲面。
[0018]可选地，所述步骤S2中，针对获取的三维翼型B样条几何模型，通过采样点离散化生成三维翼型三角网格，并提取所述三角网格中的前缘、后缘和翼尖的离散化曲线特征，其具体步骤包括：
[0019]步骤S2
‑
1：设置三角网格目标边长度；
[0020]步骤S2
‑
2：根据三角网格目标边长度计算累计长度函数和边缘累计面积函数，分别根据计算所得的累计长度函数和边缘累计面积函数在表面边界曲线上和曲面内部生成均匀采样点；
[0021]步骤S2
‑
3：根据生成的采样点，利用Delaunay剖分获得三维翼型三角网格，并提取前缘、后缘和翼尖的离散化曲线特征。
[0022]可选地，所述步骤S2
‑
2中，根据三角网格目标边长度计算累计长度函数和边缘累计面积函数，分别根据计算所得的累计长度函数和边缘累计面积函数在表面边界曲线上和曲面内部生成均匀采样点具体包括：
[0023]步骤S2
‑
21：对步骤S1生成的每个B样条曲面边界曲线，计算累计长度函数并据此估算样本点个数；使用等差累计长度对应的u为节点向量；根据节点向量对曲线进行采样得到边界采样点，并生成采样点的曲面参数坐标；筛选其中位于前缘曲线、后缘钝面边界曲线和翼尖上下沿曲线上的采样点为边界特征采样点；
[0024]步骤S2
‑
22：对步骤S1生成的每个B样条曲面，计算边缘累计面积函数并据此估算样本点个数；使用等差边缘累计面积对应的u为节点向量K
u
，使用等差边缘累计面积对应的v为节点向量K
v
；使用K
u
×
K
v
在每个B样条曲面采样；根据B样条节点参数，创建采样点的曲面参数坐标，得到曲面内部特征采样点；
[0025]步骤S2
‑
23：合并所述边界特征采样点和曲面内部特征采样点，得到所有特征采样点的集合。
[0026]可选地，步骤S2
‑
3中，根据生成的采样点，利用Delaunay剖分获得三维翼型三角网格，并提取前缘、后缘和翼尖的离散化曲线特征，具体包括：
[0027]步骤S2
‑
31：根据上翼面、下翼面、上翼尖面、下翼尖面本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于对称方向场的三维翼面四边形网格自动生成方法，其特征在于，所述方法包括如下步骤：步骤S1：导入三维翼型的截面形状参数，设定空间几何参数，计算生成控制点，并基于控制点获取B样条三维几何模型；步骤S2：针对获取的三维翼型B样条几何模型，通过采样点离散化生成三维翼型三角网格，并提取所述三角网格中的前缘、后缘和翼尖的离散化曲线特征；步骤S3：根据生成的三维翼型三角网格，计算曲线特征约束下的对称方向场，通过优化参数化梯度和方向场的误差得到三角网格的参数化映射，利用libQEx软件提取参数化坐标得到四边形网格。2.根据权利要求1所述的方法，其特征在于，所述步骤S1中，设定空间几何参数具体包括：设定包括翼根弦长、翼尖弦长、前缘后掠角、半翼展和后缘厚度在内的空间集合参数。3.根据权利要求2所述的方法，其特征在于，所述步骤S1中，基于控制点获取B样条三维几何模型的算法流程包括：步骤S1
‑
1：导入NACA翼型描述格式的三维翼型截面形状参数文本文件，根据后缘厚度设置截面翼型上沿和下沿中的后缘点，接着对上下沿靠近后缘点线性插值，使得前缘到后缘钝面之间平滑过渡，生成具有后缘钝面的截面翼型；步骤S1
‑
2：将三维翼型截面分为上下沿两部分，将其包含的翼型二维点作为控制点分别生成3阶B样条上沿曲线和下沿曲线；步骤S1
‑
3：分别对上沿曲线和下沿曲线进行等数目均匀采样，使得上下沿采样点一一对应；步骤S1
‑
4：根据翼根弦长、翼尖弦长、前缘后掠角和半翼展，计算翼尖相对翼根偏移；使用上一步骤S1
‑
3中的二维上下沿采样点创建三维翼型翼根上沿控制点、翼根下沿控制点、翼尖上沿控制点和翼尖下沿控制点，其中上下沿、翼根与翼尖的控制点具有一一对应关系；组合上沿控制点为上翼面控制点，组合下沿控制点为下翼面控制点，选择所有控制点中位于后缘钝面的四个点组成后缘钝面控制点；步骤S1
‑
5：创建节点向量，设置上下翼面B样条曲面u方向为机翼弦长方向，v方向为翼展方向，阶数为3
×
1，创建上翼面和下翼面B样条曲面；设置后缘钝面B样条曲面u方向为机翼厚度方向，v方向为翼展方向，阶数为1
×
1，创建后缘钝面B样条曲面；步骤S1
‑
6：根据上下翼面控制点计算平均水平方向向量和竖直方向向量，进而计算上下翼尖面控制点；步骤S1
‑
7：设置上下翼尖面B样条曲面u方向为机翼弦长方向，v方向为机翼厚度方向，阶数为3
×
3，创建上下翼尖面B样条曲面。4.根据权利要求1所述的方法，其特征在于，所述步骤S2中，针对获取的三维翼型B样条几何模型，通过采样点离散化生成三维翼型三角网格，并提取所述三角网格中的前缘、后缘和翼尖的离散化曲线特征，其具体步骤包括：步骤S2
‑
1：设置三角网格目标边长度；步骤S2
‑
2：根据三角网格目标边长度计算累计长度函数和边缘累计面积函数，分别根据计算所得的累计长度函数和边缘累计面积函数在表面边界曲线上和曲面内部生成均匀采样点；
步骤S2
‑
3：根据生成的采样点，利用Delaunay剖分获得三维翼型三角网格，并提取前缘、后缘和翼尖的离散化曲线特征。5.根据权利要求4所述的方法，其特征在于，所述步骤S2
‑
2中，根据三角网格目标边长度计算累计长度函数和边缘累计面积函数，分别根据计算所得的累计长度函数和边缘累计面积函数在表面边界曲线上和曲面内部生成均匀采样点具体包括：步骤S2

【专利技术属性】
技术研发人员：李妮，余延超，龚光红，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：