网格生成文献卷 04:六面体网格、体参数化与块结构

为什么这条线三十多年仍然难

一个六面体不是六个四边形面拼起来就够了。全局上,奇异边构成嵌入体内的图;局部上,体参数映射必须保持正 Jacobian;边界上,六面体面要与 CAD 特征和边界 quad 相容;数值上,还要控制 scaled Jacobian、尺寸渐变与边界层方向。

2022 年系统综述的核心判断仍然成立:已有家族分别擅长鲁棒性、规则结构、几何质量、粗网格或自动化,但没有一种方法同时在任意复杂输入上满足所有目标。

先认识五个家族

家族核心对象优势典型瓶颈
octree/template空间层次与局部模板自动、稳健、易并行表面阶梯、奇异结构复杂、单元偏多
sweeping/mapping源面到目标面的层叠块内结构和边界层优秀需要可扫掠拓扑和可靠块分解
plastering/whisker边界 quad 与组合对偶直接控制边界、强调拓扑闭合、自交和几何实现困难
polycube/体参数化到轴对齐域或整数格的映射可采样结构网格、与参数化统一polycube 构造、双射与奇异图
frame-field三维正交方向标架自动推断局部方向场可积性、奇异线可实现性和网格提取

1996–2010:组合拓扑与 polycube

级别工作历史作用主要限制
AWhisker Weaving, 1996在对偶空间逐层连接 whisker sheet,代表组合式 all-hex 思路对偶自交、闭合和几何嵌入十分困难
APolycube-Maps, 2004用轴对齐 polycube 域组织全局参数化原工作面向纹理;用于 hex 时必须增加体双射和有效性
BHexahedral Mesh Generation Constraints, 2008系统讨论边界 quad、奇异结构和可填充性的约束必要条件与可构造算法之间仍有距离
BFeature-Preserving Octree Hex Meshing, 2009以 refinement templates 改善尖锐特征处理规则性和单元数量仍受背景树控制

2011–2016:体参数化和 frame field 快速发展

级别论文解决的问题还留下什么
ACubeCover, 2011将分支覆盖和体参数化用于三维整数格结构场、奇异线和无翻转体映射的联合难题
AAll-Hex Mesh Generation via Volumetric Polycube Deformation, 2011从 polycube 体变形得到 all-hexpolycube 分区和变形双射依赖输入与优化
BBoundary Aligned Smooth 3D Cross-Frame Field, 2011构造边界对齐的三维 frame field平滑场不保证奇异结构可网格化
BPolyCut, 2013以可重用模板/切割把复杂体映射到 polycube 结构分割与模板匹配决定自动化上限
AFrame Fields: Anisotropic and Non-Orthogonal Cross Fields in 3D, 2014扩展三维 frame 表示与优化空间表示连续不等于全局可积或可提取
AAll-Hex Meshing Using Closed-Form Induced Polycube, 2016用解析诱导场构造 polycube 参数化复杂拓扑和边界有效性仍依赖后续处理
APractical 3D Frame Field Generation, 2016改善 frame field 优化的实用性与边界对齐实用场仍可能产生不可接受奇异图
CHexEx, 2016从有缺陷的体参数化中鲁棒提取 hex提取器不能替代上游映射有效性

2017–2023:奇异图、量化与基准意识

级别工作研究推进关键阅读问题
BRobust Hex-Dominant Meshing, 2017以四面体聚合获得实用 hex-dominant 输出混合单元比例、奇异结构与 FEM 表现
ASingularity-Constrained Octahedral Fields, 2018把用户/拓扑期望的奇异结构纳入场优化给定奇异线是否可被有效体参数实现
AIntrinsic Mixed-Integer Polycubes, 2022用内蕴度量和混合整数变量构造 polycube整数求解规模与复杂几何鲁棒性
AVolume Parametrization Quantization for Hex Meshing, 2022研究体参数周期/尺度的离散量化量化和局部单射需要共同维护
AHex-Mesh Generation and Processing: A Survey, 2022以输入、输出、方法家族和后处理统一三十多年工作用其分类回读原论文,不要只停在表格
BHex Me If You Can, 2022用公共模型和指标暴露不同方法的成功/失败边界数据集仍不能覆盖工业 CAD 全分布
ALocally Meshable Frame Fields, 2023从局部可网格化约束重新审视 frame field 设计局部可实现仍需拼成全局有效网格

2024–2026:拓扑优化与新的理论构造

级别工作新意当前证据边界
AHex Singularity Structure Simplification via ILP, 2024将 sheet 操作、价数约束和特征保持写成整数优化优化已有 hex,不是任意 CAD 的从零生成
BFeature-Aware Singularity Structure Optimization, 2025联合 sheet collapse/inflate,使奇异结构靠近几何特征局部贪心与复杂拓扑下的全局质量
BVolume Parametrization Quantization with Flexible Singularities, 2026放宽奇异结构以改善量化可行性新方法需在更广泛 CAD 和求解任务上验证
BHex Meshing Based on Meromorphic Quadratic Differentials, 2026用叶状结构、闭轨线和体分解连接复分析与 hex固定/可得骨架、体映射正性与通用性
BCubeFlipNet, 2026用强化学习选择 cube-flip 与关键点需和确定性 ILP/启发式在失败率及泛化上比较

奇异边为什么是核心对象

规则 hex 内部边通常有 4 个六面体相邻。价数偏离 4 的边形成奇异边,并在体内连接成图。这个图同时约束:

  • 块结构如何连接;
  • 参数坐标绕闭环的 holonomy;
  • sheet/chord 能否穿过体域;
  • 局部 Jacobian 是否有优化空间;
  • 最终求解中的方向连续性。

只优化每个 hex 的几何形状,无法自由改变这张全局图;只简化图,也不保证存在几何上无翻转的实现。

评价 all-hex 方法的最低报告表

类别至少报告
输入CAD/三角汤/四面体、流形性、亏格、特征角和容差
有效性负 Jacobian 数、最小与低分位 scaled Jacobian
结构奇异边长度/分支数、base complex 块数、边界对齐
几何Hausdorff/法向误差、尖角保持、薄层分辨率
规模单元数、内存、运行时间、人工交互分钟数
数值同自由度下误差、条件数、求解迭代和总时间

可复现资源

推荐阅读次序

  1. 先读 2022 survey,建立五大家族;
  2. 用 Whisker Weaving 和 Polycube Maps 对比“组合对偶”与“参数域”思路;
  3. 读 CubeCover → frame field → HexEx,跟踪场到网格的断点;
  4. 读奇异约束、量化和 locally meshable fields,理解可实现性;
  5. 读 ILP/sheet 优化与 feature-aware 工作,理解已有网格的后处理;
  6. 最后把学习方法和共形方法放到同一基准中,而不是按论文图示判断。

返回文献地图 · 对应教学正文