一种面向结构化四边网格生成的二维区域自动分解方法技术

本发明专利技术公开了一种面向结构化四边网格生成的二维区域自动分解方法。当前，针对复杂几何外形的结构化网格生成，主要依赖人工将复杂区域手动分解为规则的子区域，增加了出错的可能性。本发明专利技术步骤：建立分析实体对象的二维几何模型，并离散成三角形网格模型；根据几何边界约束建立二维几何模型内部区域光滑矢量场的控制偏微分方程，并采用边界元数值方法进行求解；建立与光滑矢量场对应的光滑标架场，并通过分析光滑标架场的奇异结构实现对二维几何模型的自动分解。本发明专利技术可以自动快速地获得任意复杂二维几何模型的高质量分区结果，为后续多块结构化四边形网格生成提供输入。

A two-dimensional region automatic decomposition method for structured quadrangle mesh generation

The invention discloses a two-dimensional region automatic decomposition method for generating structured quadrangle grids. At present, for the structured mesh generation of complex geometric shapes, it mainly relies on manual decomposition of complex regions into regular sub-regions, which increases the possibility of errors. The steps of the invention are as follows: establishing a two-dimensional geometric model for analyzing a solid object and discretizing it into a triangular mesh model; establishing a control partial differential equation of the smooth vector field in the inner region of the two-dimensional geometric model according to the geometric boundary constraints, and solving it by the boundary element numerical method; and establishing a smooth frame field corresponding to the smooth vector field. By analyzing the singular structure of the smooth frame field, we can automatically decompose the two-dimensional geometric model. The invention can automatically and quickly obtain high quality partitioning results of any complex two-dimensional geometric model, and provides input for subsequent multi-block structured quadrilateral mesh generation.

【技术实现步骤摘要】
一种面向结构化四边网格生成的二维区域自动分解方法
本专利技术涉及数值模拟领域前处理中的分块结构化四边形网格生成的区域分解过程，具体涉及一种面向结构化四边网格生成的二维区域自动分解方法。
技术介绍
网格生成是有限元法、有限体积法和有限差分法等数值模拟技术中的前处理过程，该过程将连续的几何区域剖分成有限个基本几何形体的组合，这些基本几何形体被称为网格单元。常用的单元类型有三角形单元、四边形单元、四面体单元和六面体单元等；根据网格单元的拓扑结构是否规则可将网格分为结构化网格和非结构化网格。在数值模拟中，结构化网格由于具有比非结构化网格更高的求解精度和更少的求解自由度而备受青睐。相对其它类型的网格，结构网格具有很强的拓扑约束。因而，针对复杂几何模型，生成一套高质量的结构网格是一项非常费时且困难的工作。通过将模型区域分解为多个子域，并在各子域中生成结构化网格，而各子域之间的拓扑不要求保持结构化特性，将得到分块结构化网格，这将在一定程度上降低网格生成难度，也能避免产生完全结构化网格中可能出现的扭曲单元。正因如此，在实际工程中，分块结构化网格比完全结构化网格应用范围更广。对于复杂几何模型，对其进行区域分解具有和生成其结构化网格一样的难度，自动、可控、高质量的区域分解至今仍是一个开放性难题。当前，针对复杂几何外形的结构化网格生成，主要依赖人工将复杂区域手动分解为规则的子区域，再在各子区域里生成规则的四边形网格。人工环节的引入一方面要求工程人员具有较强的工程经验，另一方面非常耗时且增加了出错的可能性。这些缺点制约了自动化分块结构化网格生成工具的诞生。
技术实现思路
为避免手工区域分解存在的缺点，本专利技术提供一种面向结构化四边网格生成的二维区域自动分解方法，能够针对任意复杂二维几何模型快速自动获得高质量的区域分解结果，从而提升分块结构化网格生成效率和分块结构化网格质量。本专利技术采用的技术方案是：建立分析实体对象的二维几何模型，并离散成三角形网格模型；根据几何边界约束建立二维几何模型内部区域光滑矢量场的控制偏微分方程，并采用边界元数值方法进行求解；建立与光滑矢量场对应的光滑标架场，并通过分析光滑标架场的奇异结构实现对二维几何模型的自动分解，具体包括以下步骤：步骤1、建立分析实体对象的二维几何模型，在二维几何模型区域上定义矢量场与标架场之间的映射关系，具体如下：建立分析实体对象的二维几何模型，将二维几何模型区域Ω离散成三角形网格Th，并在Th上定义矢量场和标架场，分别由u(p1)和vk(p1)表示，0≤k≤3，v0(p1)垂直于v1(p1)且与v2(p1)反向，v3(p1)与v1(p1)反向，p1为网格顶点；记θd(p1)为x轴沿逆时针旋转到与u(p1)方向相同时所转过的最小角度；θc,k(p1)为x轴沿逆时针旋转到达与vk(p1)方向相同时所转过的最小角度，其中θc,k(p1)，0≤k≤3为四个互不相等的角度值，矢量场和标架场的映射关系用式(1)表示：θd(p1)＝4θ,θ＝min{θc,0(p1),θc,1(p1),θc,2(p1),θc,3(p1)}(1)步骤2、在二维几何模型区域上建立光滑矢量场最优化模型，具体如下：以二维几何模型区域边界上矢量向区域内推进时满足最小变化为目标，以二维几何模型区域边界的几何特征约束为约束条件建立光滑矢量场最优化模型，即目标函数为约束条件为：其中，为ui(p)的拉普拉斯算子，p为二维几何模型区域内的点，i为点p的坐标分量x或y，ui(p)为u(p)在x或y方向的分量，和分别为狄利克雷边界条件和诺依曼边界条件，n为边界上的点的外法向量；ui,0(p)为二维几何模型区域边界上的矢量沿i方向的分量，为ui,0(p)的法向导数。光滑矢量场最优化模型的最优解等价为如下拉普拉斯边值问题的解：式(3)中，Δui(p)为ui(p)的拉普拉斯算子的另一种表达形式。为使二维几何模型区域分解后的区域分割线垂直于区域分解后的边界，三角形网格边界点的标架满足如下条件：式中，p0为三角形网格边界上的网格顶点，θh(p0)和n(p0)分别为三角形网格边界上以p0点为交点的两条边朝向区域内的夹角值和该夹角平分线的单位向量，t(p0)为垂直于n(p0)的单位向量，vk(p0)为p0处的标架，0≤k≤3。通过公式(4)求得三角形网格边界上的标架信息后，通过步骤1定义的矢量场与标架场之间的映射关系求解三角形网格边界上的矢量边界条件，其中，三角形网格边界矢量的模为1。步骤3、求解矢量拉普拉斯方程根据公式(3)，二维几何模型区域内部点p2处的矢量沿i方向的分量通过如下边界积分公式求解：其中，二维几何模型区域边界上的点dsp表示二维几何模型区域边界上p点邻域的微分，G(p3,p2)为Laplace方程的基本解，即n3为p3的外法向量。由公式(3)，在二维几何模型区域边界上每个点的ui(p3)和中有一个值是已知的，并且另一个值通过式(6)得到：其中，二维几何模型区域边界上的点且p4与p3为不同点，由于在点p3处光滑，则n4为p4的外法向量。采用边界元数值方法求解公式(6)，将划分成N段，记为k1在1～N中取值，N＞2，得到如下离散形式：其中，qk1为位于上的点，qj为位于上的点，nqj为qj的外法向量。公式(7)表示一个由N个未知数和N个方程组成的方程组，通过求解该方程组，得到位于二维几何模型区域边界点上的边值ui(p)或随后，通过公式(5)求解位于二维几何模型区域内的内部点处的矢量，最终获得整个二维几何模型区域内的矢量场。步骤4、建立光滑标架场并利用其实现区域分解通过步骤1定义的矢量场与标架场之间的映射关系，求得步骤3中获得的矢量场对应的标架场。通过分析该标架场获得奇异点的位置，并以该标架场为输入数据，利用龙格-库塔法寻找以奇异点为起始点和终止点的流线。最后，以这些流线将二维几何模型区域分解为多个四边结构的子区域。本专利技术的有益效果是，本专利技术通过建立覆盖几何模型区域的矢量场物理模型，并采用边界元法进行数值求解，随后建立与之一一映射的标架场，有效地增强了获得标架场的精度和速度，最终通过分析标架场的奇异结构，将二维几何模型区域精确且自动分解成多个四边形结构的子区域。该方法实现了任意复杂二维几何模型区域的自动快速分解，提高了分块结构化网格生成效率。本专利技术面向分块结构化网格生成需求，提出一种针对任意复杂二维几何模型的自动区域分解方法，该专利技术可以分析二维飞机翼形外流场区域结构并将其自动分解为多个具有四边结构的子区域，并且该专利技术也适合其他复杂二维几何模型的区域分解。附图说明图1a和图1b分别为位于p1点处的方向向量和标架；图2a为本实施例中二维几何模型区域边界示意图；图2b图为图2a中二维几何模型区域离散成三角形网格的示意图；图2c为本专利技术方法针对图2a中二维几何模型区域形成的标架场边界图；图2d为本专利技术方法针对图2a中二维几何模型区域形成的矢量场边界图；图3a为本专利技术方法针对图2a中二维几何模型区域形成的矢量场图；图3b为本专利技术方法针对图2a中二维几何模型区域形成的标架场图；图3c为本专利技术方法针对图2a中二维几何模型区域形成的多个四边结构子区域分解图；图4为本实施例中的二维飞机翼形外流场区域示意图；图5为本专利技术方法针对图4中二维飞机翼形外流本文档来自技高网...

【技术保护点】
1.一种面向结构化四边网格生成的二维区域自动分解方法，其特征在于：该方法具体包括以下步骤：步骤1、建立分析实体对象的二维几何模型，在二维几何模型区域上定义矢量场与标架场之间的映射关系，具体如下：建立分析实体对象的二维几何模型，将二维几何模型区域Ω离散成三角形网格Th，并在Th上定义矢量场和标架场，分别由u(p1)和vk(p1)表示，0≤k≤3，v0(p1)垂直于v1(p1)且与v2(p1)反向，v3(p1)与v1(p1)反向，p1为网格顶点；记θd(p1)为x轴沿逆时针旋转到与u(p1)方向相同时所转过的最小角度；θc，k(p1)为x轴沿逆时针旋转到达与vk(p1)方向相同时所转过的最小角度，其中θc，k(p1)，0≤k≤3为四个互不相等的角度值，矢量场和标架场的映射关系用式(1)表示：θd(p1)＝4θ，θ＝min{θc，0(p1)，θc，1(p1)，θc，2(p1)，θc，3(p1)}     (1)步骤2、在二维几何模型区域上建立光滑矢量场最优化模型，具体如下：以二维几何模型区域边界

【技术特征摘要】
1.一种面向结构化四边网格生成的二维区域自动分解方法，其特征在于：该方法具体包括以下步骤：步骤1、建立分析实体对象的二维几何模型，在二维几何模型区域上定义矢量场与标架场之间的映射关系，具体如下：建立分析实体对象的二维几何模型，将二维几何模型区域Ω离散成三角形网格Th，并在Th上定义矢量场和标架场，分别由u(p1)和vk(p1)表示，0≤k≤3，v0(p1)垂直于v1(p1)且与v2(p1)反向，v3(p1)与v1(p1)反向，p1为网格顶点；记θd(p1)为x轴沿逆时针旋转到与u(p1)方向相同时所转过的最小角度；θc，k(p1)为x轴沿逆时针旋转到达与vk(p1)方向相同时所转过的最小角度，其中θc，k(p1)，0≤k≤3为四个互不相等的角度值，矢量场和标架场的映射关系用式(1)表示：θd(p1)＝4θ，θ＝min{θc，0(p1)，θc，1(p1)，θc，2(p1)，θc，3(p1)}(1)步骤2、在二维几何模型区域上建立光滑矢量场最优化模型，具体如下：以二维几何模型区域边界上矢量向区域内推进时满足最小变化为目标，以二维几何模型区域边界的几何特征约束为约束条件建立光滑矢量场最优化模型，即目标函数为约束条件为：其中，为ui(p)的拉普拉斯算子，p为二维几何模型区域内的点，i为点p的坐标分量x或y，ui(p)为u(p)在x或y方向的分量，和分别为狄利克雷边界条件和诺依曼边界条件，n为边界上的点的外法向量；ui，0(p)为二维几何模型区域边界上的矢量沿i方向的分量，为ui，0(p)的法向导数；光滑矢量场最优化模型的最优解等价为如下拉普拉斯边值问题的解：式(3)中，Δui(p)为ui(p)的拉普拉斯算子的另一种表达形式；为使二维几何模型区域分解后的区域分割线垂直于区域分解后的边界，三角形网格边界点的标架满足如下条件：式中，p0...

【专利技术属性】
技术研发人员：徐岗，何寿平，肖周芳，许金兰，吴卿，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江,33