4. 面向数字孪生的几何建模与网格策略

数字孪生体要承载结构、流体、热、电磁等多物理场计算,几何模型与网格就是它的“空间骨架”。几何决定物理域边界在哪里、耦合面如何对应;网格决定方程在何处离散、精度与算力如何平衡。本节围绕参数化几何、轻量化处理、多域网格划分与网格映射四条主线,说明如何为多物理场集成准备可复用、可对齐、可计算的空间离散基础。

一、参数化几何:让孪生体“可调、可追溯”

面向数字孪生的几何不应是一次性画死的实体,而应是由参数驱动的模型。以一台工业换热器为例,可把壳体长度 L=1200 mm、管束外径 d=16 mm、管间距 p=22 mm、折流板间距 B=200 mm 定义为顶层参数,管束阵列数量、折流板缺口位置由这些参数通过表达式派生。这样当工况或设计变更时,只需修改参数即可重建几何,孪生体与物理实体保持同步。

参数化建模的关键是建立“参数—特征—约束”的层级关系:顶层参数控制草图尺寸,草图约束控制拉伸、阵列、抽壳等特征,特征再生成最终实体。建议遵循三条实践原则:其一,命名规范统一,如 P_Shell_L、P_Tube_d,便于脚本批量读写;其二,关键尺寸全部参数化,避免“硬编码”导致改一处崩一片;其三,保留特征树与版本记录,使几何变更可回溯。参数化几何的价值不仅在于快速改型,更在于为后续网格与物理场设置提供稳定的命名接口——边界名、面名、体名一旦固定,网格与求解器配置就能脚本化复用。

二、轻量化处理:在保真与算力之间取舍

CAD 模型常包含螺栓孔、倒角、小圆角、铭牌、装饰面等对物理场影响很小的细节。若直接划分网格,会产生大量畸变单元,显著增加计算量甚至导致求解失败。轻量化处理的目标是“去掉不影响物理结论的细节,保留影响流动、传热、应力路径的关键特征”。

常用手段包括:小特征抑制(小于网格尺寸 1/3 的圆角、倒角可移除)、孔洞填充(非关键安装孔可封堵)、面合并(共面碎片面合并为单一面)、中面抽取(薄壁结构用壳单元替代实体)、以及包围盒或凸包简化(用于外部流场域)。以换热器壳体为例,若最小网格尺寸设为 3 mm,则半径小于 1 mm 的倒角可安全抑制;而折流板缺口、管板开孔等直接影响流动与传热的特征必须保留。轻量化不是“越简越好”,而是以物理场敏感度为判据:先做一次粗网格试算,观察哪些细节处梯度变化剧烈,再决定保留或简化。

轻量化前后应记录几何体积、表面积、特征数量等指标,并对比关键物理量(如压降、最高温度)的变化。若简化后压降偏差小于 3%、温度偏差小于 2%,通常可认为简化可接受。

三、多域网格划分:不同物理场用不同策略

多物理场集成的难点在于:结构场需要贴合壁面的高阶单元,流体场需要边界层加密,热场需要与流体共轭传热,电磁场可能需要规则六面体。若强行用同一套网格,往往顾此失彼。合理做法是“分域划分、界面共享”。

以换热器管束区域为例:流体域采用多面体或六面体核心网格,壁面处布置 5~8 层边界层,首层高度按 y+≈1 估算;固体管壁采用扫掠六面体,沿壁厚方向 3~5 层;管板与壳体连接处采用四面体过渡。各域在耦合界面(管外壁面、管板端面)上必须保证节点或面片一一对应,否则插值误差会破坏能量守恒。

网格质量指标需重点关注:正交性(建议大于 0.2)、长宽比(边界层外不宜超过 100)、雅可比(大于 0.6)、以及最大偏斜度(小于 0.85)。对于瞬态问题,还需控制全局网格量,使单步求解时间与孪生体更新频率匹配。例如孪生体要求 1 s 更新一次,而单步求解需 0.8 s,则网格量已接近上限,需通过局部加密与自适应细化来平衡。

四、网格映射:让多域数据“对得上、传得准”

网格映射解决的是“不同网格之间如何传递物理量”。常见场景包括:流体网格计算出的壁面温度映射到结构网格作为热载荷;结构网格计算出的位移映射到流体网格作为动网格边界;电磁网格的损耗密度映射到热网格作为热源。映射方法主要有三类:节点插值(最近邻、线性、径向基函数)、面片守恒映射(面积加权)、以及体映射(体积加权)。

工程上优先选择守恒型映射,确保总热量、总力等积分量不因插值而凭空增减。映射前需检查两套网格的空间重叠度:若耦合面法向偏差超过 5°,或面片覆盖率低于 95%,应重新划分或调整几何。映射后应做“守恒性校验”:分别对源面和目标面积分,比较总量偏差,一般要求小于 1%。

几何—网格—映射:多物理场空间离散链路 参数化几何 L=1200 mm d=16 mm p=22 mm B=200 mm 参数→特征→约束 轻量化处理 抑制 r<1 mm 倒角 封堵非关键孔 合并共面碎片 保留:折流板缺口 多域网格划分 流体:边界层 y+≈1 固体:扫掠六面体 过渡:四面体 正交性>0.2 偏斜<0.85 网格映射 温度→结构热载荷 位移→流体动网格 损耗→热源 守恒偏差<1% 耦合界面网格映射示意(管壁—流体—结构) 流体域网格 多面体核心 + 边界层 首层高度 0.02 mm 5~8 层边界层 输出:壁面温度 T_w 壁面热流 q_w 固体管壁网格 扫掠六面体 壁厚方向 3~5 层 输入:T_w、q_w 输出:位移 u、应力 σ 温度场 T_s 结构域网格 高阶四面体/六面体 管板连接处加密 输入:T_s、压力载荷 输出:应力、变形 反馈:位移→动网格 T_w, q_w T_s, u 位移反馈 → 流体动网格更新(双向耦合)

五、从几何到网格的工程检查清单

把上述方法落到实际项目,建议在提交计算前完成以下检查:几何参数是否全部命名并可从脚本访问;轻量化后关键物理量偏差是否在允许范围;各域网格质量指标是否达标;耦合界面覆盖率与法向偏差是否合格;映射守恒性校验是否通过。只有这些环节都闭环,多物理场集成才有可靠的空间离散基础。

检查项典型判据不通过时的处理
几何参数化关键尺寸均有命名参数补充参数与表达式,重建特征树
轻量化偏差压降偏差<3%,温度偏差<2%恢复敏感细节或减小网格尺寸
网格正交性>0.2调整节点分布或改用多面体
边界层 y+≈1(壁面解析)调整首层高度与层数
耦合面覆盖率>95%重新划分或修补几何缝隙
映射守恒偏差<1%改用守恒型映射或加密目标网格

几何与网格策略不是一次性的预处理,而是贯穿孪生体全生命周期的空间离散基础。参数化让模型可调,轻量化让计算可行,多域划分让各物理场各得其所,网格映射让数据在域间准确传递。四者协同,才能支撑后续结构、流体、热、电磁等多物理场的稳定集成与实时更新。