​​有限元网格划分技术：原理、方法与前沿综述​​

有限元网格划分是进行有限元数值模拟分析的关键预处理步骤，其结果（单元形状、拓扑类型、单元类型选择、网格密度、节点编号以及几何建模方式）直接决定了后续数值计算的精度和效率。在有限元求解过程中，单元的等效节点力、刚度矩阵、质量矩阵等通常通过数值积分生成：连续体单元以及壳、板、梁单元的面内方向多采用高斯(Gauss)积分；壳、板、梁单元的厚度方向则常采用辛普森(Simpson)积分。因此，网格质量是保证计算精度的基石。

有限元方法的核心思想在于结构离散化——用简化的几何单元逼近连续体，并在满足变形协调条件下求解。构建合理的有限元模型需兼顾几何描述准确性与变形梯度描述的精确性。以下为网格划分需遵循的关键原则：

1. 其曲线/曲面边界能更精确地贴合实际几何轮廓。

2. 高阶插值函数能更准确地逼近复杂场函数。
   因此，增加单元阶次可有效提升计算精度。然而，高阶单元单个单元节点数增多，同等网格数量下模型总规模（总自由度）更大，计算耗时增加。选择时需综合权衡精度要求与计算成本。

1. 纵横比（Aspect Ratio）:

   单元边长比（1D）、面积比（2D）或体积比（3D），以正三角形、正四面体、正六面体为理想标准。比值越接近1越好。

   
   

2. 扭曲度（Skewness / Distortion）:

   表征单元面内角度或面外翘曲的畸变程度。

   
   

3. 单元内部角:

   (特别是三角形单元)避免出现过尖或过钝的内角。
   此外，节点编号策略对求解过程中总体刚度矩阵的带宽或波前规模影响显著，进而影响求解耗时和存储需求，也属网格质量优化范畴。

1. 单元的节点必须同时也是相邻单元的节点，不能是其内部点或孤立的边界点。

2. 相邻单元在共享节点处具有相同类型（性质）的自由度。
   需注意，节点自由度数一致是协调的必要条件但非充分条件，单元形函数在共享边界（线或面）上的连续性也需保证。

有限元网格生成方法可按多种方式分类（如生成单元类型、几何维度、自动化程度）。下文结合典型方法分析其机理、特点与适用范围：

• 机理：

  基于实际几何域与参数空间标准模板（如正方形、立方体）的双向映射。

  
  

1. 利用映射函数将复杂物理域映射到规则参数域。

2. 在规则参数域内生成结构化网格。

3. 将参数域网格映射回物理空间得到最终网格。

• 特点：

  可生成高质量、方向一致的结构化网格（全四边形/六面体）。

  
  

• 局限：

  对几何适应性差（复杂形体无法直接映射），需繁琐的人工子区域分解（将目标域拆分为若干可映射的子域），局部网格控制能力弱，各子域网格密度设定相互掣肘，难以实现自动化（尤其3D问题）。

  
  

• 适用范围：

  形状规则或可被方便分解为规则子域的结构（如管道、简单的拉伸体）。

• 机理：

• 保留完全在目标域内的栅格单元。

• 删除完全在目标域外的栅格单元。

• 裁剪与边界相交的栅格单元以精确逼近边界。

1. 生成一个覆盖目标区域的背景栅格（Cartesian Grid）。

   识别栅格单元：

   对内部栅格单元直接网格化，对裁剪后的边界栅格单元进行网格化或进一步细分。

   特点：

   自动化程度高，局部网格控制能力较强（通过局部栅格细化），易于处理非流形边界。

   
   

   局限：

   边界逼近精度依赖栅格单元大小（需权衡精度与效率），边界处单元质量可能较差（生成过渡棱锥、棱柱等单元或需细化），生成的网格单元大小通常依赖于背景栅格。

   
   

   适用范围：

   自动生成四面体网格，或作为六面体网格生成（如结构化模板）的基础。通用性好。

   3.3 节点连元法

   
   

• 机理：

1. 在目标区域（边界和内部）按预设密度分布节点（布点算法是关键）。

2. 基于几何拓扑关系和特定连接准则（如 Delaunay 三角化空球/空圆准则、推进波前法）连接节点形成单元（如三角形、四面体）。

• 特点：

  概念直接，适用于任意形状区域，自动化程度高（尤其 Delaunay 法和 Advancing Front 法）。

  
  

• 局限：

  2D 三角形网格技术成熟，3D 四面体网格质量控制和生成速度仍是挑战；不易生成高质量的全四边形/六面体网格。

  
  

• 适用范围：

  复杂几何的非结构化三角形/四面体网格自动生成。常用作基础网格生成方法。

• 机理：

  从形体的拓扑特征（如顶点、边、环、孔洞）入手，忽略具体几何形状细节，将目标域分解为若干简单的拓扑单元（如三角形、四边形面片），然后对各拓扑子域网格化。核心目标是使用最少的基本单元（如三角形）完全覆盖目标域。

  
  

• 特点：

  强调拓扑特征识别，处理含有孔洞等复杂拓扑结构的几何具有一定优势。

  
  

• 局限：

  实际应用常需结合几何信息进行细化；分解策略较抽象，分解质量影响后续网格质量。

  
  

• 适用范围：

  特别适用于拓扑结构复杂的 2D 区域（如带孔面域）的分解引导网格生成。

• 机理：

  区别于拓扑分解，此法更侧重于利用几何特征（如尖角、窄带、曲面曲率变化）进行分解。节点和单元通常在分解过程中同步生成。

  
  

• 特点：

  分解策略紧密联系几何特征，通常能直接生成质量较高的单元，边界拟合度好。

  
  

• 局限：

  算法复杂度较高，通用化自动生成难度大，对于极其复杂的几何仍需人工干预引导。

  
  

• 适用范围：

  尤其适用于对单元质量要求高、几何特征显著的区域的四边形/六面体网格生成（如特征驱动方法）。

• 机理：

  将二维截面（通常已划分好四边形网格）通过平移、旋转、缩放等扫描操作生成三维实体网格（主要为六面体或棱柱体）。

  
  

• 特点：

  操作直观，在CAD/CAE软件中普遍集成，容易理解和操作，能高效生成结构化或部分结构化的六面体网格。

  
  

• 局限：

  应用范围受限，仅适用于具有明显拉伸、旋转或扫掠路径特征的体（即扫掠体），不支持任意三维实体，自动化程度不高（常需人工定义源/目标面、路径线等）。

  
  

• 适用范围：

  典型扫掠实体（如型材、旋转对称体、沿路径拉伸体）的六面体网格划分。

  
  

有限元分析在工程领域的深入应用推动了网格生成技术的发展。当前研究重心已从二维平面转向三维实体，并聚焦于四边形(Quadrilateral)和六面体(Hexahedral)网格的全自动生成技术、网格自适应技术、高质量（边界和内部）网格技术等。

• 映射法：

  先将复杂体人工或自动分解为可映射的区块（类似于大六面体），然后在每个区块内映射细小的六面体单元。挑战在于自动化分解复杂拓扑。

  
  

• 单元转换法：

  将初步生成的其他类型单元（如四面体）转换为六面体单元，或将低阶六面体转换为高阶单元。

  
  

• 基于栅格法（如Octree/Block-based）：

  利用背景网格或多块结构化模板覆盖目标体并处理边界相交，常需辅助单元（棱锥、棱柱）或进一步分解生成全六面体网格。

  
  

• 多子区域法（Multi-block Method）：

  将复杂区域分解为简单的拓扑子块，对各子块独立划分六面体网格（常用映射法），最后组合全局网格。关键在子域分解与交界处网格协调。

  
  

• 扫描法：

  (前文3.6所述) 适用于具有扫掠特征的几何。

  
  

• 投影法：

  (一种模板法) 利用模板六面体网格或高质量四面体网格作为基础，通过投影或插值到目标体并调整节点位置。需处理投影失真和边界贴合问题。

  
  

• 现状：

  高质量、全自动、适应任意复杂几何的六面体网格生成仍是公认难题，是当前研究的前沿热点。

• 直接法：

  在物理空间曲面直接操作网格划分算法（如曲面Delaunay、曲面推进波前法）。该方法直接依据曲面的局部几何性质（法向、曲率）确定节点位置和单元连接策略，无需参数化映射。

  
  

• 优点：

  网格质量较高，能更好地保持局部曲率特征。

  
  

• 缺点：

  算法实现复杂，对曲面自交、高度扭曲等情况处理难度大。

  
  

• 映射法（或称参数化法）：

• 优点：

  可利用成熟的平面网格生成器。

  
  

• 缺点：

  网格质量高度依赖于参数化的均匀性和保角性（参数域形状和间距失真会传递到物理空间曲面网格）；参数化本身对复杂拓扑曲面（如多洞环状曲面）可能是NP难问题。

  
  

• 现状：

  参数化方法更成熟、应用更广，尤其适合拓扑简单的曲面；直接法因更少受参数化扭曲影响而在高质量曲面网格生成中日益受到重视。

  
  

1. 建立曲面到二维参数域的可接受映射（如UV参数化）。

2. 在规则参数域中生成高质量平面网格（通常为三角形或四边形）。

3. 将参数域网格节点反向映射回物理空间得到曲面网格。

尽管有限元网格技术日趋成熟（尤其在三维四面体网格的全自动化方面），以下领域仍存在挑战和广阔发展空间：

1. 高度自动化全六面体网格技术：

   克服复杂拓扑和多尺度特征障碍，实现高质量全六面体网格的完全自动生成是终极目标之一。需要更智能的几何分解、单元转换及模板生成方法。

   
   

2. 高效高质量曲面网格生成：

   提升直接法的鲁棒性和对高复杂度曲面的处理能力；改进参数化方法使其更适应复杂拓扑和曲率变化；发展自动生成全四边形曲面网格的稳健方法。

   
   

3. 网格自适应技术：

   发展更精准可靠的误差估计方法，并基于此实现高效、高鲁棒性、无缝的网格自适应细化/粗化策略（特别是对六面体和非线性分析）。

   
   

4. 几何复杂性和鲁棒性：

   增强网格算法处理工程现实中存在的“脏”几何（如小特征、缝隙、重叠面、不精确模型）的健壮性（Robustness），减少修复预处理需求。

   
   

5. 人工智能/机器学习辅助：

   探索应用深度学习和强化学习等AI技术优化网格生成（如智能布点、预测高质量网格参数、优化网格密度分布、生成六面体分解策略、快速判断网格质量等）。

   
   

6. 多尺度与混合网格：

   高效处理多尺度问题，将局部高分辨率网格与全局粗网格无缝耦合；发展更优化的混合网格（如四面体与六面体混合）接口技术。

   
   

7. 并行化与大规模网格处理：

   优化并行算法以适应超大规模模型的高效网格生成和自适应。

   
   

8. 无网格法结合：

   探索网格生成与无网格法(如SPH, Peridynamics)的前后处理流程结合点，或在多物理场、大变形分析中提供互补解决方案。

   
   

9. 网格技术的可访问性与易用性：

   在商业和开源CAE软件中进一步集成和优化先进的网格技术，简化用户操作，提升工程师的使用体验和效率。