多物理场 — CAE术语解释

不是的，那只是"多仿真"。真正的多物理场是指多个物理场相互强烈影响的"耦合现象"，在控制方程的层面上结合它们进行求解。例如，智能手机的发热问题涉及半导体芯片的电力损耗（电学），产生热量（热传导和对流），热量导致材料膨胀（结构），反过来应力又改变电气特性。这些不是独立的，必须同时求解。

单向耦合在某些情况下可行，但双向耦合不足以应对。以制动盘的褪色现象为例：摩擦热使温度升高，热膨胀改变接触压力，而这又改变摩擦热的产生率。在这种情况下，热传导控制方程和结构运动方程中，都包含彼此的变量（温度T和位移u）。热应力项的形式如

可以用无量纲数。例如在流体-结构耦合（FSI）中，可以检查附加质量系数，或结构固有频率与流体特征时间的比。具体来说，密度比

这是求解方法的根本性差异。单体式（一体型）求解法将所有物理场的未知数（温度、位移、压力等）合并为一个大矩阵，通过一次线性求解器调用同时求解。分割式（分割型）求解法为各物理场分别准备求解器，把一个物理场的解作为另一物理场的边界条件，交替求解，重复直到收敛。COMSOL Multiphysics的"Fully Coupled"求解器是前者，Ansys的"System Coupling"是后者的代表。

强耦合问题（如前面提到的制动器或水中结构）中，单体式有优势。收敛快且稳定，但内存消耗大，矩阵条件数易恶化。分割式可以复用现有的单物理场求解器（Fluent、Mechanical），易于实现，内存效率高。但如果不用弱耦合加速方法（如Aitken加速），收敛会很慢，甚至可能发散。实务中，Ansys Workbench的热应力分析（单向）采用分割式，压电-结构（双向）采用单体式的"Piezo and MEMS"模块，根据现象灵活选择。

最近邻插值会严重损失精度。标准做法是"保守插值"。例如将流体网格的力传递到结构网格时，在流体网格的单元面上积分的力，按一定权重分配到结构网格的节点。这样做是为了保证作用-反作用定律（流体对结构的总力与结构对流体的总力相等）。MpCCI（基于网格的耦合代码接口）标准或Ansys System Coupling的"Conservative Transfer"采用这个方法。非保守传递会导致能量人为增减，引起振荡。

第一步是定性和定量估计耦合的"强度"。以这个例子为例：1) 风扇驱动的强制对流（流体）→ 2) 基板上的芯片发热被风吹散（热传递）→ 3) 基板热膨胀（结构）。这里要问：热膨胀会大幅改变流体区域的形状吗？通常答案是否，所以热-流体是双向的，对结构的影响是单向的。因此工作流可以是"用CFD求解流场和热传递"→"将热传递系数和温度传给结构分析"，两步就够了。

至少要确认以下5点：1) **单位统一**：坚持SI单位制。特别是CFD导入mm单位的CAD时，会出现kg/mm³这样的非标准密度，容易混淆。2) **物性温度依赖性**：材料的热导率或比热会随温度变化吗？变化超过10%就必须定义。3) **界面设置**：不同物理场接触的面（固体-流体界面等）是否正确应用了耦合条件（如"热耦合"或"FSI界面"）？4) **网格适配性**：数据转移的面上，网格粗细是否差异过大。5) **收敛判定**：分割式中，除了单个物理场的残差，界面数据转移的残差（如界面力变化小于0.1%）也要作为收敛条件。

分割式要怀疑**松弛系数（欠松弛系数）**。先试试大幅降低初始值，如0.75或0.5。特别是密度或刚性差异大的物理场耦合（水和结构等）时，默认值1.0易发散。其次是**时间步**。过渡分析中，除了Courant数，还要考虑耦合的更严格条件。单体式发散时，首先要调细非线性求解器的设置，比如增量加载法的步数。

对，核心理念完全不同。
**COMSOL Multiphysics**："基于方程"的方法。用户可以直接编辑偏微分方程（PDE）的"系数形式PDE"界面，在统一的框架（接近单体式）下求解各种物理场。电磁、化学反应、传热等都在同一网格上密耦合。
**Ansys**："最优品质整合"方法。流体用Fluent/CFX，结构用Mechanical，电磁用Maxwell等各领域最强专用求解器，通过Workbench平台进行耦合（分割式）。各求解器性能高，但密耦合有限制。
**Abaqus**："以结构为中心的扩展"。起源于强大的非线性结构分析，扩展支持热传导、电磁场（Abaqus/CAE+SIMULIA协同仿真）。主要通过SIMULIA协同仿真引擎进行分割式耦合。

原型设计和原理验证阶段，**COMSOL**压倒性优势。拥有"AC/DC模块""结构力学模块""微流体模块"等，可在一个模型文件中定义和耦合所有物理场。特别是微小位移改变静电容量，进而影响驱动力的那种双向耦合，相对容易设置。Ansys虽然有"Piezo and MEMS"专用工具包，但设置复杂。到了量产设计阶段，需要与现有ECAD/MCAD流程集成，**Ansys**或**Abaqus**的生态建立更成熟，优势就大了。

两种模式。一种是**专用求解器开发**。OpenFOAM包含`chtMultiRegionFoam`（传热-流体耦合）和`pimpleFluid`（流体-结构）等专用求解器，这些在C++代码内硬编码了接近单体式的方程形式。另一种是**高级描述库**。FEniCS和FreeFEM可用Python风格脚本组合弱形式表示的各物理场方程，自动生成和求解联立方程组。这类似于COMSOL的开源版本，灵活性很高，但计算效率优化需由用户负责。

这是"人工刚度不足"引起的典型不稳定。特别是薄结构（膜、薄板）受流体作用时，单靠结构求解器的刚度矩阵，面外变形的抵抗会被严重低估。对策有两个。第一，Ansys Mechanical中启用"**虚功线性项**"，考虑压力加载的几何刚性效应。第二，检查**固体单元和壳单元混合模型**。只用壳单元会只有面内刚度，丧失了现实存在的弯曲刚度。可以在厚度方向加一层固体单元，或给壳单元正确定义弯曲刚度。

几乎肯定是"**参考温度（零应力温度）**"设置错误。热应力

是时间尺度不匹配。电磁分析通常定常（AC）或非常短的时间步（微秒级），热传导时间尺度是秒～分钟级。简单耦合会导致电磁分析的1步时间内产生的发热，在热传导分析的1步中瞬间加入。若热传导1步是1秒，就相当于1秒的发热瞬间加入，无法表现连续加热。对策是先求电磁分析得到平均发热密度（W/m³），作为**体积热源用于定常热传导**，或在过渡热传导中，**正确考虑发热随时间的积分**，给出正确的负荷曲线。检查软件中的"Time Scaling Factor"。