焦耳加热的电磁-热耦合分析

焦耳加热电磁-热耦合的理论基础

焦耳加热就是电流流过时产生的热量对吧？那为什么还要专门做耦合分析呢？

问得好。电流流过确实会产生热量，这是对的。但这还不是全部故事。温度升高会导致金属电阻率增加。以铜为例，每升高1°C电阻就增加约0.4%。这样的话，同样的电流就会产生更多的热量，温度继续上升，电阻继续增加……这样就形成了正的反馈回路。要准确求出定常温度，必须对电磁场和热场进行反复迭代耦合。

那温度不会一直上升吗？

实际上，存在散热（对流、辐射、传导）。发热和散热最终达到平衡，温度就稳定在某个值。但那个值到底是多少度，只有通过"电磁场→发热→温度变化→电阻率变化→电磁场"这样的循环反复计算才能求出。这就是为什么焦耳加热的耦合分析这么重要。

具体在什么场合这会成为问题呢？

最典型的就是EV（电动汽车）的800V逆变器接线端子。400A以上的大电流流过母线，局部温度要确保不超过150°C。超过这个值，绝缘材料会提前失效，寿命缩短。如果忽略电阻率的温度依赖性，计算结果会比实际低20～30°C，设计缺陷就会被遗漏。还有保险丝熔断特性的预测、功率半导体基板设计、输电线的下垂（挠度）计算等，这些都是典型的应用场景。

为什么需要耦合分析

焦耳加热分析不能用简单的"P = I²R"计算完成，核心原因是电导率 $\sigma$ 是温度 $T$ 的函数。定电流条件下，虽然电流不变，但温度升高会改变电阻率，进而改变电流密度分布。这种双向耦合如果被忽视，分析结果会与测量值产生显著偏差：

• 定常温度的低估： 忽视温度依赖性时，铜母线的温度可能被低估20～30°C
• 电流集中的遗漏： 截面变化部、接线部会因局部电阻增加而产生热点
• 过渡响应误差： 短路电流等过渡大电流的温度上升曲线因电阻率变化而呈非线性

控制方程

那实际的方程是怎样的呢？

首先是电磁场方面。定常电流的情况下，从电荷守恒律得到连续方程。发热量由电流密度和电场的内积确定：

电场的控制方程（定常电流近似）：

单位体积的焦耳发热量：

$Q = J^2/\sigma$ 和 $Q = \sigma E^2$ 是同一个式子吗？ 看起来有点矛盾啊...

看出细节的好眼光。利用 $\mathbf{J} = \sigma \mathbf{E}$ 代入，就会得到相同的式子。具体用哪个形式看方便程度——电流密度已知时用 $J^2/\sigma$；电位已知时用 $\sigma E^2$。实务中常常以电位为未知量求解，所以 $\sigma |\nabla V|^2$ 这种形式用得比较多。

热场的控制方程（能量守恒）：

这里 $\rho$ 是密度，$c_p$ 是定压比热，$k$ 是热导率，$Q$ 是焦耳发热量（热源项）。

电阻率的温度依赖性

金属的电阻率随温度的增加几乎成线性关系（金属在室温附近）：

其中 $\rho_0$ 是基准温度 $T_0$ 处的电阻率，$\alpha$ 是温度系数。

$\alpha$ 的值对不同材料差异大吗？

差异很大。总结成表格一目了然：

镍铬的 $\alpha$ 特别小呢。加热丝用它的原因是温度升高了还能稳定发热吗？

正是这样！ 反过来，铜和铝的 $\alpha$ 大，所以在母线和电缆的设计中，温度依赖性是不能忽视的。温度上升100°C时，铜的电阻会增加约40%。设计时要是把这个忽略了，分析就没用了。

耦合的反馈结构

焦耳加热电磁-热耦合中存在明确的正反馈回路：

1. 电流密度 $\mathbf{J}$ 产生的发热 $Q = J^2/\sigma(T)$
2. 温度 $T$ 上升
3. 电阻率 $\rho(T)$ 增加（$\sigma(T)$ 减少）
4. 相同电位差下电流分布发生变化，局部发热密度改变
5. 回到第1步

定电流源和定电压源在这方面会有不同的表现吗？

眼光敏锐。定电压源情况下：温度上升→电阻增加→电流减少→发热减少，形成自限制效应。反过来，定电流源的情况：温度上升→电阻增加→发热 $Q = I^2 R$ 增加→温度更高，呈失控趋势。功率电子中，器件驱动条件由两者混合组成，所以边界条件的设置错误就会导致完全不同的结果。

焦耳加热电磁-热耦合的数值计算方法

弱耦合和强耦合

听说有"弱耦合"和"强耦合"的区分，具体是怎么用的？

焦耳加热分析的耦合方式大致分为3种：

实务中用得最多的是哪一种？

压倒性地是弱耦合（逐次反复法）。有3个原因。首先，现有的电磁场求解器和热求解器可以直接使用。其次，便于调试——各个物理场可以独立验证。再次，大规模模型时内存消耗少。单体法因为自由度翻倍，内存压力很大。

弱耦合的迭代算法（Gauss-Seidel型）：

1. 设定初始温度 $T^{(0)}$（通常取环境温度）
2. 用 $\sigma(T^{(k)})$ 求解电场方程，得 $V^{(k+1)}$
3. 计算 $Q^{(k+1)} = \sigma(T^{(k)}) |\nabla V^{(k+1)}|^2$
4. 用 $Q^{(k+1)}$ 作为热源项求解热传导方程，得 $T^{(k+1)}$
5. 收敛判定：若 $\|T^{(k+1)} - T^{(k)}\|_\infty < \varepsilon$ 则终止，否则返回第2步

收敛判定的 $\varepsilon$ 应该设成多少？

实务中按温度绝对值基准设 0.1°C～1°C 比较合适。EV母线那样容差只有10°C的设计，设成0.1°C。加热器设计那样宽松度数十°C的，1°C就足够了。通常5～15次迭代就会收敛。

有限元制定

电场的弱形式：

这里 $N_i$ 是形状函数，$J_n$ 是边界上的法向电流密度。

热场的弱形式：

右边第2项是对流边界条件（$h$：热交换系数，$T_\infty$：环境温度）。

电场和热场用的是同一个网格吗？

焦耳加热分析中，电场和热场都在同一区域存在，所以基本上用同一网格。这一点和流体-结构耦合（FSI）不同，后者需要考虑不同网格间的插值误差。COMSOL和Ansys Workbench都会自动用同一网格做耦合计算。

时间积分和收敛控制

过渡分析中，电场的响应时间（$\tau_e \sim \varepsilon/\sigma \sim 10^{-18}$ s）和热场的响应时间（$\tau_T \sim \rho c_p L^2/k \sim 1$～$100$ s）差异巨大，所以需要特殊处理。

焦耳加热电磁-热耦合的实务应用

分析流程和边界条件

开始进行焦耳加热耦合分析的时候，首先要做什么？

实务中的分析流程是这样的：

1. 造型建模： 从CAD提取导体部分。螺栓孔、倒角、圆角不要省略（影响电流集中）
2. 材料定义： 输入 $\rho(T)$、$k(T)$、$c_p(T)$ 的温度依赖表。特别是 $\rho(T)$ 是必须的
3. 电气边界条件： 电流入口面设电流密度或全电流值，出口面设GND（$V = 0$）
4. 热边界条件： 外表面设自然对流（$h = 5$～$15$ W/m²K），需要时加入辐射（$\varepsilon = 0.1$～$0.9$）
5. 耦合设定： 设置弱耦合的迭代次数和收敛基准
6. 后处理： 确认最高温度、电流密度分布、热流分布

对流的热交换系数 $h$ 怎么决定？ 教科书里5～25的范围太宽了…

实践中这样分类用：

• 密闭筐体内的自然对流： $h = 3$～$8$ W/m²K
• 开放空间自然对流（水平面朝上）： $h = 5$～$15$ W/m²K
• 开放空间自然对流（竖直面）： $h = 8$～$20$ W/m²K
• 强制风冷： $h = 25$～$250$ W/m²K（依风速）

保险做法是初次分析用 $h = 10$，再做灵敏度分析，改变 $\pm 50$%后结论是否改变。改变的话说明需要CFD耦合或实测。改变不了就可以放心。

网格划分策略

网格应该怎么划分？

焦耳加热特有的网格策略是电流密度梯度大的地方要细化：

• 截面变化部： 母线宽度或厚度变化处。电流集中，网格至少3层
• 螺栓孔周边： 电流绕过孔时密度急变。孔径的1/5以下的要素为目标
• 接触面： 评估接触电阻时，接触面至少要有2层薄要素
• 表皮效应（交流电流）： 表皮深度 $\delta = \sqrt{2/(\omega \mu \sigma)}$ 内至少3个要素

相反，导体中央部电流均匀的长直段，粗网格就OK。网格的精明配置，计算成本的80%都由此决定。

接触电阻建模

母线螺栓连接部的接触电阻，应该怎么建模？

这是实务上最难的地方。接触电阻依赖于接触面压、表面粗糙度、氧化膜的状态，数值会变动1～2个数量级。建模方法有3种：

1. 薄层要素： 接触面插入厚数μm的虚拟薄层，用其电导率表示接触电阻
2. 界面电阻条件： COMSOL的"Electric Contact"或Ansys的"Contact Resistance"，直接指定面的电阻
3. 等价电路模型： 接触部当作集中常数电阻，发热量手工加入

无论哪种方法，输入的接触电阻值必须实测。铜-铜接触通常是 $10^{-5}$～$10^{-3}$ $\Omega$·cm² 数量级，但表面处理和螺栓扭矩会大幅改变。

实务案例：EV母线

EV开发中，焦耳加热耦合分析实际怎么用？

典型例子是800V电池组母线设计。分析目标是"连续400A通电时，母线温度不超过120°C（绝缘材耐热温度-30°C裕度）"的确认。

分析条件：

• 材料：无氧铜 C1020（$\rho_0 = 1.72 \times 10^{-8}$ $\Omega$·m、$\alpha = 3.93 \times 10^{-3}$ /K）
• 截面：30 mm × 3 mm（母线）、M6螺栓4根连接
• 电流：400 A DC连续
• 环境：60°C（筐体内温度）、$h = 8$ W/m²K（自然对流）

分析结果（有无温度依赖比较）：

竟然差29°C！ 不考虑温度依赖就OK，但考虑就NG了…

正是这样。这种情况下的对策是要么把截面从3 mm加到4 mm，要么改变螺栓孔的布置来缓解电流集中。参数分析来找最优截面形状是常见做法。实际有个T1供应商用这个分析把重量削减15%、温度降低5°C同时实现的。

焦耳加热电磁-热耦合的软件对比

支持工具和功能对比

焦耳加热耦合分析用什么软件能做？

焦耳加热是相对简单的耦合，支持的软件很多：

新手首选用哪个？

COMSOL圧倒的入易。「Electric Currents」「Heat Transfer in Solids」的2選、自動的連成設定。「Busbar Heating」名前的的的的例題。COMSOL単位高額（基本+追加×2）的、会社大規模展開Ansys Workbench的方場合。

COMSOL 最容易上手。选择"Electric Currents"和"Heat Transfer in Solids"两个模块，耦合自动设置。教程里有专门的"Busbar Heating"例题。但COMSOL按模块收费（基础+追加模块），成本较高。大规模推广的话，Ansys Workbench的性价比可能更好。