1. 课程导入与边界一致性检查总览
在多物理场仿真中,我们很少只面对一个孤立的物理场。热-结构耦合、流-固耦合、电磁-热耦合、声-结构耦合,几乎每一个真实工程问题都要求两个甚至更多物理场在同一套几何、同一套网格、同一套时间推进逻辑下协同求解。而多物理场模型边界一致性检查,正是针对这种"协同"所必须做的第一道质量把关:它要确认不同物理场之间共享的边界在几何位置、网格离散、变量传递和求解时序上是否真正对齐,从而保证耦合界面上的物理量能够被正确、守恒、稳定地交换。
需要强调的是,边界一致性并不等于"几何看起来重合"。一个热-结构模型里,热分析用的实体单元表面和结构分析用的壳单元中面,在 CAD 里可能完全贴合,但离散后一个是面网格、一个是线网格,节点根本对不上;一个流-固耦合模型里,流体域出口和固体域入口在几何上共面,但两侧网格疏密相差数倍,插值误差会直接污染界面力。因此,边界一致性检查的对象是"几何—离散—求解"三个层面的对齐状态,而不是单一的几何重合度。
为什么边界不一致会带来严重后果
边界不一致引发的工程问题,通常以三种形式暴露出来,且往往在仿真后期才被发现,返工代价很高。
- 误差放大:耦合界面上的插值误差会随迭代累积。例如流-固耦合中,若固体侧界面网格比流体侧粗,流体压力向固体节点的映射会丢失高频分量,导致界面力分布失真,结构响应出现系统性偏差。
- 求解发散:强耦合问题(如热-结构双向耦合)对界面能量守恒极为敏感。当两侧界面网格不匹配且传递算法未做守恒修正时,界面总热流可能不守恒,迭代残差无法收敛,表现为耦合迭代振荡甚至发散。
- 物理失真:即使数值上收敛,错误的边界条件也可能给出"看起来合理但物理错误"的结果。典型如电磁-热耦合中,若损耗映射的边界区域与热分析加热区域错位,热点位置会整体偏移,热应力评估随之失效。
下面这张图概括了从几何到离散再到求解的检查框架,以及各层不一致所对应的典型后果。它把"检查什么"和"不一致会怎样"放在同一条链路上,便于建立整体判断。
三层检查框架的具体含义
几何层关注的是"边界在哪里"。检查内容包括:共享面或共享线是否在同一位置、法向是否一致(法向反向会导致界面力符号错误)、两侧模型是否使用同一坐标系与同一单位制。工程上常见的坑是:热模型用毫米、结构模型用米,导入后几何尺寸相差千倍,表面看是"网格不匹配",实质是单位不一致。因此几何层检查应先于任何网格操作。
离散层关注的是"边界被怎样切开"。核心问题是界面两侧网格是否匹配。匹配分两种策略:一是共形网格(conformal),即界面节点一一对应,传递无插值误差;二是非共形网格(non-conformal),通过映射算法(如最近邻、形函数插值、守恒型映射)传递。前者精度高但网格生成约束强,后者灵活但必须验证映射精度与守恒性。检查时应记录界面节点数、单元类型(面单元/线单元)、网格尺寸比,并确认映射算法与物理量类型匹配——例如热流密度这类通量量适合守恒型映射,而温度这类状态量适合形函数插值。
求解层关注的是"边界在求解过程中如何被使用"。包括耦合变量的传递方向与顺序(单向还是双向)、时间步长是否一致或存在子循环、强耦合迭代的收敛判据是否包含界面残差。一个典型检查项是:双向耦合中,界面残差是否被纳入收敛判据;若只检查各物理场自身残差,可能出现"各自收敛、耦合不收敛"的假收敛。
一个可操作的检查顺序示例
以热-结构双向耦合为例,假设热分析采用四面体实体单元、结构分析采用壳单元,界面为一块平板的上表面。合理的检查顺序是:先确认两者共享面的几何位置与法向一致、单位统一;再检查界面处热侧面网格与结构侧壳网格的节点分布,若采用非共形映射,记录映射算法并做一次"单位热流输入—界面总热流输出"的守恒性抽检;最后确认耦合迭代中温度与热流是否双向传递、时间步是否同步、界面残差是否进入收敛判据。每一步都应有明确的通过/不通过判据,而不是凭经验"看着差不多"。