一种复杂三维曲面共型网格生成方法技术

本发明专利技术公开了一种复杂三维曲面共型网格生成方法，包括以下步骤：S1.给定规则曲面，并在XOY平面上进行投影得到矩形区域，并确定横向与纵向的剖分长度，进行XOY平面的网格剖分；S2.将剖分得到的每一个网格的中心点坐标从Z平面映射到W平面；S3.确定空间曲边四边形到平面三角形单元的映射方案，并生成共形平面三角形网格。生成的复杂曲面网格在二维平面上的投影网格具有保角特性，能够实现曲面结构的平面三角形网格剖分问题，且完成了四边形网格到三角形网格的转换。角形网格的转换。角形网格的转换。

【技术实现步骤摘要】
一种复杂三维曲面共型网格生成方法


[0001]本专利技术涉及网格生成，特别是涉及一种复杂三维曲面共型网格生成方法。

技术介绍

[0002]随着计算机软硬件技术的飞速发展,工业界普遍使用先进的CAD和CAE技术构建目标电磁几何模型和网格模型以完成产品的模拟仿真以预测与验证产品的电磁性能,不仅提高了产品设计的灵活性,也显著降低了产品的研发周期与经济成本.使用CAD技术构建产品模型,再采用CAE技术对产品进行模拟仿真以评估与验证产品性能,已成为工业界研发产品的一个标准流程。然而,CAD通常采用连续的边界表示法(boundary representation,B
‑
Rep)表示实体模型,CAE通常以离散的网格模型作为数据输入；同时,由于高质量的网格模型对电磁仿真环节至关重要,因此如何精确有效地构造高质量的且有一定物理和几何意义的离散网格模型,成为电磁计算过程中一个至关重要的问题。

技术实现思路

[0003]本专利技术的目的在于克服现有技术的不足，提供一种复杂三维曲面共型网格生成方法，生成的复杂曲面网格在二维平面上的投影网格具有保角特性，能够实现曲面结构的平面三角形网格剖分问题，且完成了四边形网格到三角形网格的转换。
[0004]本专利技术的目的是通过以下技术方案来实现的：一种复杂三维曲面共型网格生成方法包括以下步骤：
[0005]S1.给定规则曲面，并在XOY平面上进行投影得到矩形区域，并确定横向与纵向的剖分长度，进行XOY平面的网格剖分；
[0006]S2.将剖分得到的每一个网格的中心点坐标从Z平面映射到W平面；
[0007]S3.确定空间曲边四边形到平面三角形单元的映射方案，并生成共形平面三角形网格
[0008]本专利技术的有益效果是：本专利技术由于共形网格具有保形性，保留了几何网格变形前后的夹角属性，如果变形之前的网格线之间相互正交，则变形后的网格线也具有正交性，这样的正交特性能够方便电磁计算数值算法(比如矩量法和有限元法)在网格上的离散化处理，简化数值计算过程，而且网格线的正交性能够极大地降低系统矩阵的规模，降低计算耗时和资源，进一步提升计算精度。
附图说明
[0009]图1为本专利技术的方法流程图；
[0010]图2为XOY平面矩形网格示意图；
[0011]图3为Z平面到W平面的映射原理示意图；
[0012]图4为XOY平面共形网格示意图；
[0013]图5为四边形网格与三角形网格映射关系示意图；
[0014]图6为实施例中采用方法一得到共形网格示意图；
[0015]图7为实施例中采用方法二得到共形网格示意图。
具体实施方式
[0016]下面结合附图进一步详细描述本专利技术的技术方案，但本专利技术的保护范围不局限于以下所述。
[0017]如图1所示，一种复杂三维曲面共型网格生成方法，包括以下步骤：
[0018]S1.给定规则曲面，并在XOY平面上进行投影得到矩形区域，并确定横向与纵向的剖分长度，进行XOY平面的网格剖分；
[0019]S101.已知规则曲面及边界显示表示为：z＝f(x,y)x1≤x≤x
2 y1≤y≤y2；曲面在世界坐标系XYOZ的XOY平面上的投影为一长L＝x2‑
x1且宽W＝y2‑
y1的矩形区域，则有该区域中的x和y坐标映射为Z平面，即Z＝x+Jy，其中其中，Z平面就是XOY平面；X轴就是Z平面的实轴，Y轴就是Z平面的虚轴；这是复分析的一种惯用的坐标建立方式；
[0020]S102.已知XOY平面上的剖分最大长度为Δl，则有横向的剖分数nL＝MAX(10,ceil(L/Δl))，纵向剖分数nW＝MAX(10,ceil(W/Δl))，其中MAX为两个数取大值函数以保证横纵轴的最小网格数不低于10，ceil为向上取整函数；得到横向的剖分长度Δx＝L/(nL+1)，纵向的剖分长度Δy＝W/(nW+1)；
[0021]S103.由横向的剖分长度Δx建立x1≤x≤x2范围内的包括端点在内的(nL+1)个顶点和nL个区间范围，每个顶点组成了纵向网格线，第j个纵向网格线表示为：
[0022][0023]由纵向的剖分长度Δy建立y1≤y≤y2范围内的包括端点在内的(nW+1)个顶点和nW个区间范围，每个顶点组成了纵向网格线，第i个纵向网格线表示为：
[0024][0025]由(nL+1)个纵向网格线与(nW+1)个横向网格线交错形成(nL+1)
×
(nW+1)个顶点和nL
×
nW个矩形区间范围的矩形网格单元；如图2所示；
[0026]图中，设G
ij
为第i个纵向网格线和第i+1个纵向网格线，第j个横向网格线以及第j+1个横向网格线围成的矩形区域；P
ij
(x
ij
,y
ij
)为G
ij
的中心点二维坐标，即：
[0027][0028][0029]S2.将剖分得到的每一个网格的中心点坐标从Z平面映射到W平面；
[0030]如图3所示，已知Z平面到W平面的共形映射函数为g，该映射表示为W＝g(Z)＝ξ+Jη；
[0031]则有ξ＝Re(g(x+Jy))和η＝Im(g(x+Jy))，其中Re为取实部函数，Im为取虚部函数；映射得到的W面曲面网格如图4所示：
[0032]对于任一个网格有G
ij
，中心点坐标P
ij
(x
ij
,y
ij
)从XOY平面映射到ξOη平面的C
ij
(ξ
ij
,η
ij
)，则有：
[0033]ξ
ij
＝Re(g(x
ij
+Jy
ij
))和η
ij
＝Im(g(x
ij
+Jy
ij
))；
[0034]映射方法采用如下任意一种：
[0035]映射方法一：先利用x
ij
和y
ij
计算z
ij
＝f(x
ij
,y
ij
)，再利用x
ij
和y
ij
计算ξ
ij
和η
ij
[0036]则有得到三维映射(x
ij
,y
ij
,z
ij
)
→
(ξ
ij
,η
ij
,f(x
ij
,y
ij
))
[0037]映射方法二：先利用x
ij
和y
ij
计算ξ
ij
和η
ij
，再利用z
ij
＝f(ξ
ij
,η
ij
)计算本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种复杂三维曲面共型网格生成方法，其特征在于：包括以下步骤：S1.给定规则曲面，并在XOY平面上进行投影得到矩形区域，并确定横向与纵向的剖分长度，进行XOY平面的网格剖分；S2.将剖分得到的每一个网格的中心点坐标从Z平面映射到W平面；S3.确定空间曲边四边形到平面三角形单元的映射方案，并生成共形平面三角形网格。2.根据权利要求1所述的一种复杂三维曲面共型网格生成方法，其特征在于：所述步骤S1包括：S101.已知规则曲面及边界显示表示为：z＝f(x,y)x1≤x≤x2y1≤y≤y2；曲面在世界坐标系XYOZ的XOY平面上的投影为一长L＝x2‑
x1且宽W＝y2‑
y1的矩形区域，则有该区域中的x和y坐标映射为Z平面，即Z＝x+Jy，其中其中，Z平面就是XOY平面；S102.已知XOY平面上的剖分最大长度为Δl，则有横向的剖分数nL＝MAX(10,ceil(L/Δl))，纵向剖分数nW＝MAX(10,ceil(W/Δl))，其中MAX为两个数取大值函数以保证横纵轴的最小网格数不低于10，ceil为向上取整函数；得到横向的剖分长度Δx＝L/(nL+1)，纵向的剖分长度Δy＝W/(nW+1)；S103.由横向的剖分长度Δx建立x1≤x≤x2范围内的包括端点在内的(nL+1)个顶点和nL个区间范围，每个顶点组成了纵向网格线，第j个纵向网格线表示为：由纵向的剖分长度Δy建立y1≤y≤y2范围内的包括端点在内的(nW+1)个顶点和nW个区间范围，每个顶点组成了纵向网格线，第i个纵向网格线表示为：由(nL+1)个纵向网格线与(nW+1)个横向网格线交错形成(nL+1)
×
(nW+1)个顶点和nL
×
nW个矩形区间范围的矩形网格单元；设G
ij
为第i个纵向网格线和第i+1个纵向网格线，第j个横向网格线以及第j+1个横向网格线围成的矩形区域；P
ij
(x
ij
,y
ij
)为G
ij
的中心点二维坐标，即：即：3.根据权利要求1所述的一种复杂三维曲面共型网格生成方法，其特征在于：所述步骤S2包括：已知Z平面到W平面的共形映射函数为g，该映射表示为W＝g(Z)＝ξ+Jη；则有ξ＝Re(g(x+Jy))和η＝Im(g(x+Jy))，其中Re为取实部函数，Im为取虚部函数；对于任一个网格有G
ij
，中心点坐标P
ij
(x
ij
,y
ij
)从XOY平面映射到ξOη平面的C
ij
(ξ
ij
,η
ij
)，则有：和映射方法采用如下任意一种：映射方法一：先利用x
ij
和y
ij
计算z
ij
＝f(x
ij
,y
ij
)，再利用x
ij
和y
ij
计算ξ
ij
和η
ij
则有得到三维映射(x
ij
,y
ij
,z
ij

【专利技术属性】
技术研发人员：李尧尧，曹成，
申请(专利权)人：东莘电磁科技成都有限公司，
类型：发明
国别省市：