一种新的扫掠方法.docxVIP

一种新的扫掠方法 有限元法是有效解决各种复杂工程问题的近乎数值解决方案的基本方法。它已成为许多领域工程计算和分析的重要手段，如结构、水流、加热、磁体、3d数字印刷等。有限网格的划分是通过有限元朝法进行分析的先决条件。工作量大，复杂的时间。其中，六部分网格在计算精度、划分数量、抗变形性和重段号的方面比其他网格具有优势。近年来，它成为金元网格分布领域的研究热点，但六部分网格的低模型适应性导致自动生成的问题没有得到真正解决。t）。 对复杂几何体进行几何分解, 在各分解块上运用映射法和扫掠法生成六面体网格是目前主流的网格划分方式, 其中扫掠法是最稳定、通用的六面体网格生成方法［３］.运用扫掠法生成六面体网格首先需要将几何体各面区分为源面、目标面与扫掠面3种.扫掠面要求必须实现结构网格 (排列有序、相邻节点位置关系和数量明确的四边形网格) 划分, 但是扫掠面通常连接其他不同形状的几何体, 这些连接体在扫掠面上形成子域边界, 使扫掠面边界结构变得复杂, 导致结构网格划分困难, 从而阻碍扫掠法的应用. 为了解决这一问题, 苏财茂等［４］运用约束引导线和节点等参光顺法实现了含有子域约束面的结构网格划分, 但是该方法对于三角形子域或多边形子域仍然会产生非结构网格;Liu等［５］针对含有子域的扫掠面进行分层, 然后在各层间填充结构网格, 但该方法应用的范围有限;Ruiz-Gironés等［６］将含有子域的几何面映射到计算空间以简化边界的复杂度, 然后对该空间内的几何面进行结构网格填充, 但该方法对于光滑边界或具有尖角的子域无法实现全结构网格划分:上述方法均致力于对含有子域约束的扫掠面进行结构网格划分, 但是受制于结构网格对节点位置和连接关系的限制, 仍不能完全实现子域约束扫掠体的六面体网格自动划分.基于前期我们对六面体网格划分研究［７－８］, 本文突破扫掠法应用中要求扫掠面必须是结构网格的限制, 将几何体的面网格进行分解并重组, 对各组内的面网格单独进行扫掠生成六面体网格, 实现子域约束扫掠体的六面体网格划分. 1 扫描和海滩法 1.1 扫路网格划分 扫掠法通常适用于形状类似于扫掠体的3D实体.扫掠体, 也称为2.5D实体, 由源面、目标面和扫掠面三部分组成.按照源面与目标面对应的数量关系, 扫掠体可分为一对多、多对一与多对多3类.一对一扫掠法是各类复杂扫掠法的基础, 诸如多对一以及多对多的扫掠体通常分解为简单的一对一扫掠体的组合.运用扫掠法对一对一扫掠体进行六面体网格划分的步骤如下: Step1.生成源面的四边形网格MＯ; Step2.将源面的四边形网格MＯ投影至目标面形成拓扑结构相同的网格MＴ; Step3.在扫掠面上生成结构网格MＬ, 形成边界节点层; Step4.计算源面与目标面各层之间的内部节点, 生成六面体网格. 1.2 结构网格划分 第1.1节所述的扫掠法能够成功应用于一对一扫掠体, 源面与扫掠面的网格关系必须满足以下2个条件: 1) 源面网格与扫掠面网格相连; 2) 源面与扫掠面角度θ近似为90°. 在满足上述要求的情况下, 运用扫掠法进行六面体网格划分的过程中, 依据第1.1节的Step3, 需要对扫掠面进行结构网格划分.设扫掠面网格MＬ的某节点V与其他节点相连的个数为D (V) .为了保证扫掠面形成的四边形网格能够形成规则的边界节点层, 依据文献, 结构网格节点必须满足以下2个条件 (如图1所示, 其中, 黑色实心节点表示扫掠面边界网格节点, 空心节点表示扫掠面非边界网格节点) : 1) 扫掠面边界上的网格节点V要求D (V) =2或D (V) =3; 2 扫抢面网格分解 当扫掠体的源面SＯ、目标面SＴ与扫掠面SＬ含有不同形状的子域约束时, 会包含子域表面 (如图2中的子域SＣ１, SＣ２和SＣ３) , SＣ２和SＣ３的存在使扫掠面SＬ不易形成规则的结构网格, 导致第1.1节Step3中无法生成规则的边界节点层, 从而影响Step4中各层内部节点的生成.为此本文将扫掠面网格进行分解, 将原扫掠面的结构网格部分单独扫掠生成六面体网格, 设置原扫掠面的非结构网格部分为源面网格, 利用先前扫掠过程中生成的部分节点作为边界层节点扫掠生成剩余的六面体网格, 整个流程如图3所示. 2.1 扫路s,s,s 按照第1.1节的方法划分扫掠面网格时, 受扫掠面上子域形状的约束, 子域内部网格通常会出现局部节点违背第1.2节结构网格节点条件的情况 (如图4b中的网格节点V０, D (V０) =3) .为此本文将扫掠面网格进行分解与重组. 以图4a所示的几何体为例 (该类几何体属于文献特别指出的不能用扫掠法处理的形状) , 源面SＯ含有三角形约束子域SＣ２, 扫掠面为SＬ１~SＬ５, 其中扫掠面SＬ１含有三角形子域SＣ１.对源面SＯ及扫掠面SＬ１,