基于体素有限元法、共形分解有限元法与虚拟域方法的复杂几何处理

《Finite Elements in Analysis and Design》：Treating complex geometries using voxel-based FEM, CDFEM, and fictitious domain method

　　

研究人员针对计算机辅助设计(CAD)与计算机辅助工程(CAE)集成中存在的网格生成瓶颈问题，系统评估了三种基于水平集(level-set)表征的非拟合网格方法在图像驱动仿真中的性能。研究表明，传统有限元法(FEM)依赖贴合边界的协调网格，在处理增材制造复杂晶格结构、生物医学影像衍生几何时面临几何保真度低与预处理耗时长的挑战，网格划分可占复杂仿真总耗时的80%。体素有限元法(voxel-based FEM)虽能通过体素到六面体单元的映射实现影像数据到分析模型的高效转换，但其阶梯状几何近似易引发应力振荡。共形分解有限元法(Conformal Decomposition FEM, CDFEM)通过将背景网格单元沿水平集界面剖分为协调子单元，显著提升了几何精度，但存在病态条件数风险且高阶扩展困难。虚拟域方法(fictitious domain method)采用非贴合背景网格，结合子网格(sub-grid)积分、非负矩拟合(non-negative moment fitting)及ALGOIM高度函数积分规则，实现了几何表征与分析解耦。在相同代码环境(mfem库)下的二维与三维热学、力学算例表明：对于图像基问题，CDFEM与虚拟域方法的精度均优于体素有限元法；光滑解场景下，高阶虚拟域方法的收敛效率超越CDFEM；针对规则几何，ALGOIM积分规则在计算效率上表现突出。该研究为自动化网格生成与复杂几何高精度仿真提供了方法选择依据。

《Finite Elements in Analysis and Design》刊发的这项研究聚焦于复杂几何数值模拟的效率与精度瓶颈。在现代工程工作流中，CAD与CAE的集成是航空航天、汽车、生物医学及能源领域数值仿真的基石，但从CAD模型到仿真可用网格的转换仍是持久瓶颈，常伴随数据丢失、几何误差及繁重的预处理工作。研究显示，复杂仿真中网格划分可占总分析时间的80%，间隙、干涉与重叠等问题进一步降低关注量估算质量。这一问题在增材制造中尤为突出，建筑材料和晶格结构的几何复杂性对传统网格技术构成严峻挑战，尖角、未倒圆角边缘及随机微结构若未被精确表征，将导致应力集中与早期失效；而解析此类细节所需的细密网格会显著提升计算成本，限制大规模或高分辨率分析的可扩展性。除增材制造外，基于影像的计算在生物医学工程与材料科学中日趋重要，CT或MRI衍生的仿真需精确表征多孔骨结构、软组织界面及微结构材料等复杂内部几何，无缝转换影像数据为保真网格对患者特异性医疗器械设计、组织生物力学及非均质复合材料分析至关重要。传统有限元法(FEM)虽稳健通用，但在应用于几何复杂或影像衍生域时存在固有局限——其依赖与域边界对齐的协调网格，在涉及复杂几何、形状优化、大变形或流固耦合的场景中难以平衡精度与计算效率，手动网格细化与几何简化虽常用，却可能引入额外误差与延迟。为应对这些挑战，研究人员系统评估了三种基于水平集(level-set)方法（该方法通过标量函数ψ的零等值线隐式表征界面，可自然适应拓扑变化）的非拟合网格方法：体素有限元法(voxel-based FEM)、共形分解有限元法(Conformal Decomposition FEM, CDFEM)与虚拟域方法(fictitious domain method)。研究在同一代码环境（基于开源高性能有限元库mfem实现）下，通过二维与三维热学、力学算例，对比了三类方法在精度、收敛性及计算成本上的差异，明确了各方法的适用场景。

研究采用的关键技术方法包括：1. 水平集(level-set)几何表征技术，通过符号距离函数隐式描述界面，支持线性至高阶插值；2. 体素有限元法，通过标记体素分类构建计算网格，支持结构化网格与四叉树/八叉树自适应；3. CDFEM，将背景网格中与界面相交的单元剖分为协调子单元（保留四边形/六面体子单元以改善方程条件数）；4. 虚拟域方法，采用非贴合背景网格，集成三种积分策略——递归细分相交单元的子网格(sub-grid)积分、保证权重非负的矩拟合(moment fitting)积分，以及基于隐式高度函数的ALGOIM积分规则；5. 数值验证体系，涵盖二维（带孔方板、花瓣形自由面、渗流域均匀化、多材料问题、泡沫影像力学计算）与三维（含球孔Laplace问题、陀螺微结构、工业齿轮箱扫描模型）算例，通过h收敛与p收敛分析量化性能。

研究结果如下：

3.1 二维算例

3.1.1 带圆孔的方板：稳态热学问题验证显示，体素法因自由表面正则性低，无论多项式阶数如何均呈亚优化收敛；CDFEM在线性近似下最优收敛，二次近似收敛阶提升但仍亚优化，三次近似收敛性下降；虚拟域方法在线性至三次近似下均最优收敛，子网格积分在几何精度足够时可实现指数收敛，ALGOIM积分规则在全p范围内保持指数收敛。

3.1.2 花瓣形自由面方板：该几何包含大曲率变化，体素法仍无最优收敛；CDFEM在线性与二次近似下最优收敛，三次近似仍未达最优；虚拟域方法在p=1,2,3时收敛阶优于CDFEM，且随p呈指数收敛；ALGOIM积分规则虽降低 quadrature点数量，但对光滑几何的效率优势在此未显现，子网格积分经压缩后效率可与之相当。

3.1.3 渗流域数值均匀化：针对含40个随机孔洞的二维域，以有效热导率张量Λ为指标，CDFEM作为参考收敛最快，体素法收敛缓慢，虚拟域方法收敛速度介于两者之间；对角项Λ_ii因能量单调性呈单调收敛，非对角项收敛非单调。

3.1.4 多材料问题：针对同心双材料圆环，仅线性CDFEM达到预期收敛率；体素法因几何误差、虚拟域法因无法表征材料界面的弱间断（需单位分解富集或分层基函数修正）均呈亚优化收敛；二次CDFEM仅在预渐近阶段改善收敛率。

3.1.5 力学影像计算：基于836×613像素泡沫影像的线弹性分析显示，三类方法均收敛至相同有效模量，虚拟域方法收敛最快，体素法最慢；低阶方法严重低估应力水平，体素法出现应力振荡。

3.2 三维算例

3.2.1 三维泊松问题：含球形孔洞的单位立方体热学分析中，所有方法在p=1时均达最优收敛；二次近似时体素法亚优化，CDFEM与虚拟域法最优收敛（后者误差更低）；三次近似时仅虚拟域法最优收敛；ALGOIM积分规则对光滑几何效率显著。

3.2.2 三维陀螺微结构：基于STL三角网格重构水平集的力学分析显示，三类方法收敛至同一有效模量Λ₃₃，体素法因几何解析困难收敛最慢，CDFEM在几何捕获后收敛较快，虚拟域法p收敛效率最高，自由度节省一个量级。

3.2.3 工业三维案例：针对Artec 3D重型卡车变速箱扫描模型（1,145,406个面片），水平集网格化耗时不足1分钟；采用20×20×20四阶四面体网格与3层积分细化的虚拟域法求解热学问题，温度场与热通量场无需后处理即呈现光滑性，能量在四次近似时已基本实现收敛。

讨论与结论部分，研究人员指出：体素有限元法实现最简单，但CDFEM与虚拟域技术在效率上显著更优——低阶近似下CDFEM表现突出，高阶近似下虚拟域方法性能更优，数值均匀化中亦呈现相同趋势。虚拟域方法的积分策略选择至关重要：子网格积分成本高于ALGOIM，但现有压缩技术可缓解开销；ALGOIM规则依赖几何的高阶水平集描述，而影像采集数据固有的非光滑性使其难以满足这一要求，高阶近似易引发几何振荡或需过度平滑处理，反而降低几何保真度。研究未涉及共形界面结构自适应网格加密(CISAMR)与X-Mesh等新兴方法，建议未来可在本文基准问题上开展对比评估。该研究为CAD/CAE集成中的复杂几何处理提供了系统的性能对比数据，尤其为图像驱动仿真的方法选择提供了实证依据，推动了非拟合网格方法在工业级复杂几何分析中的应用。