5.2 耦合策略取舍：从框图到求解器

5.2 耦合策略取舍：从框图到求解器

本节摘要：多物理场的工程决策集中在两处：耦合是单向还是双向，求解用全耦合还是分离式。本节给出耦合框图画法与策略选择表，用软件里的"多物理场"分支演示预定义耦合的挂接，为下一节的铜排实战备好决策依据。

先画框图，再谈策略

拿到多场问题，第一动作不是开软件，是画耦合框图：把涉及的场各画一个方框，箭头表示作用方向，然后盯住三个问题——有哪些箭头？箭头围不围成圈？圈上的增益有多大？以铜排为例：电流场生热（箭头一），温度改变电阻率（箭头二），温度引发热膨胀（箭头三）。前两个箭头围成圈，是个带反馈的闭环；第三个箭头只是出口。增益也好估：铜电阻率温度系数 0.0039/K，温升几十开尔文，电阻变化一成上下——弱反馈，但没弱到可以装看不见。框图画完，策略已经定了大半，这五分钟是全章性价比最高的时间。

图：铜排问题耦合框图与求解策略

预定义耦合：软件替你写弱形式

5.1 节说过耦合即分块矩阵，本节看它在界面上的落点。新建模型时在模型向导的物理场选择页按住 Ctrl 同时勾选"电流""固体传热""固体力学"，完成向导后，模型树里会多出一个"多物理场"分支。右键它 → 添加"电磁热"耦合，选择范围盖住铜排域——这个节点自动在传热方程里生成焦耳热源项 Q = Q_j，在电流方程里保持温度相关的材料参数；再添加"热膨胀"耦合，选好结构域与温度场，热应变项自动进入固体力学方程。两个节点，五分钟，三场联立的方程组就位，这就是"软件替你写弱形式"的兑现。

没有预定义耦合的场景怎么办？两种出路。单向作用（A 影响 B，B 不回头）最简单：先解 A，在 B 的源项里直接引用 A 的解变量，比如热源填一个含温度的解析表达式。双向但无现成节点：用非局部耦合算子（平均、最大值）与变量表定义公共中间量，把 A 的输出转成 B 认识的表达式。不过动手手搓之前先翻一遍多物理场分支的耦合清单——电磁热、热膨胀、非等温流动、流固耦合、声固耦合这些高频组合都有成品，手搓是最后的手段，不是炫技的舞台。

全耦合与分离式：拿什么换什么

同一套方程组，求解器菜单里有两条路。全耦合把三场的未知数排成一列联立求解，一次牛顿迭代同时更新温度、电位、位移。好处是收敛行为整体可控——场间的反馈不必靠外层迭代来来回回，强耦合问题（反馈增益大、非线性猛，如流固耦合大变形、相变）几乎只能走它；代价是内存：分块矩阵的非对角块要整体参与分解，内存开销随未知数总数增长。分离式把方程组按场拆开轮流解：先解电流，把焦耳热递给传热；传热解完把温度递回电流；如此轮转直到场间传递的量不再变化。好处是每一步只解一个小方程组，内存友好、单步快，还能给不同场配不同的求解策略；代价是外层迭代轮数不可控，反馈强时轮转十来次才收敛，甚至干脆发散。

选择逻辑因此归一句话：看反馈增益与内存预算。铜排这类弱反馈问题两条路都通，分离式更快省内存；风扇叶片的流固大变形、结冰相变这类强反馈，直接全耦合别犹豫。判断不准时有个便宜的试错顺序：先分离式跑，观察收敛历史里场间量的变化幅度，轮数多且不稳就升级全耦合——升比降容易，别一上来就摆满汉全席。

本节要点回顾

• 框图五分钟定策略：数箭头、看圈、估增益，围成圈且增益大才考虑全耦合。
• 预定义耦合节点自动生成场间弱项，电磁热与热膨胀是本案例的两块拼图；手搓表达式只用于单向或无成品节点时。
• 全耦合换收敛稳、付内存；分离式换内存省、付迭代轮数；铜排属弱反馈，两可，先试分离式。
• 升级容易降级难：先分离式观察收敛历史，不稳再全耦合。
• 策略想清楚了，下一节把铜排案例从头算到尾，三笔手算对账一个不少。

耦合策略的取舍本质上是仿真精度与计算成本的经济学。单向耦合成本低、计算快、结果稳定，但它忽略了反向影响——比如温度场影响材料属性从而影响应力分布，再反过来影响变形和接触状态，这个闭环在单向耦合中被人为切断。双向耦合把这个闭环打开，代价是每次迭代都要交替求解两个物理场，收敛时间和内存消耗成倍增长。判断用哪个的决策标准：如果反馈效应的量级小于正向效应的十分之一，单向就够了；如果反馈效应与正向效应相当（如流固耦合中流体压力显著改变结构变形），必须用双向。还有一个中间态叫"伪双向"——用上一轮单向的结果作为下一轮的输入做外部迭代，不做内部同时求解，这种方式计算成本介于两者之间，收敛性也更好控制。选择的原则是"用最低成本回答问题"——不是精度越高越好，而是精度与问题需求匹配就好。给初学者的建议：先做单向看结果，如果发现反馈效应不可忽略（比如结构变形改变了流道截面），再升级为双向——逐级升级比一步到位更稳健。

补一个耦合策略选择的实际案例：压电传感器的仿真。压电材料同时具有电学和力学响应——施加电压产生变形（逆压电效应），施加力产生电压（正压电效应）。这个双向耦合的本质是电-力强耦合：电场直接驱动变形，变形直接改变电荷分布，两个方向都不能忽略。在 COMSOL 中使用压电接口（已内置电-力双向耦合），不需要手动设置耦合项——但你需要理解耦合强度由压电系数矩阵（d 矩阵）决定，不同压电材料的 d 系数差异很大（石英约 2 pC/N、PZT-5H 约 590 pC/N）。仿真时的注意点：压电耦合的频率响应是仿真的核心输出——共振频率处变形急剧放大，这个频率由材料的弹性、密度和几何尺寸决定。压电传感器的设计优化就是在"共振频率"与"灵敏度"之间找平衡——灵敏度正比于 d 系数，共振频率反比于尺寸和密度。掌握了压电的双向耦合，你就打开了 MEMS 仿真（微镜、加速度计、超声换能器）的大门。

一个耦合选择的量化判断方法值得记录：定义一个"耦合强度比"指标——反馈效应产生的物理量变化与正向效应引起的物理量变化之比。以热-流耦合为例：温度升高导致流体密度下降（正向），密度下降引起浮升力驱动的自然对流（反馈），自然对流又改变温度分布。耦合强度比就是自然对流带来的温度变化除以不加自然对流时的温度变化——如果这个比值小于 10%，单向耦合（不含自然对流）的误差就在工程可接受范围内。这个比值的计算不需要双向耦合仿真——你可以先做单向得到一个温度分布，再手动估算浮升力大小，从而得到耦合强度比的量级估计。这种"先估算再决定"的方法比直接上双向耦合省时省力，而且估算过程中的物理分析本身就会加深你对耦合机制的理解。把耦合强度比的概念和估算方法写进你的仿真方法论笔记，它会在你面对每一个多物理场问题时帮你快速做出"要不要双向"的判断。