📚 课程目录第4章热传导与热辐射基础

4. 热传导与热辐射基础

在电磁制动器的电磁热耦合仿真中，电磁场计算给出的是焦耳损耗、涡流损耗与磁滞损耗的分布，而温度场计算要回答的是“这些损耗把哪里的温度抬升了多少”。两者之间的桥梁，就是热传导、对流与辐射三类换热机制，以及描述它们的偏微分方程与边界条件。本节把这条桥梁拆开讲清楚：先建立热传导方程，再逐项说明对流与辐射的处理方式，最后落到热物性参数的温度依赖性上——这是电磁热耦合区别于普通传热问题的关键。

一、核心原理：能量守恒与热传导方程

对制动器中的任一微元体做能量平衡：单位时间进入微元体的净热量，加上内热源生成的热量，等于微元体内能的增量。由此得到瞬态热传导方程（傅里叶导热定律 + 能量守恒）：

ρ·c_p·(∂T/∂t) = ∇·(k∇T) + q_v

其中：
  ρ      密度            kg/m³
  c_p    定压比热容      J/(kg·K)
  k      导热系数        W/(m·K)
  T      温度            K 或 ℃
  q_v    体积热源密度    W/m³   ← 电磁损耗耦合项
  t      时间            s

稳态时左端为零，方程退化为 ∇·(k∇T) + q_v = 0。对电磁制动器而言，q_v 主要来自三部分：线圈铜损（I²R）、制动盘或衔铁中的涡流损耗、以及铁磁材料的磁滞损耗。在仿真中，q_v 由电磁场求解器逐单元映射到热网格上，因此热网格与电磁网格的匹配精度会直接影响温升结果。

热传导方程本身只描述固体内部导热。要让方程有唯一解，必须补充初始条件（t=0 时的温度分布）和边界条件（表面如何与外界换热）。

二、对流换热与热辐射定律

制动器外表面与空气之间的换热通常用牛顿冷却定律描述：

q_conv = h·(T_s − T_∞)

  h      对流换热系数   W/(m²·K)
  T_s    表面温度       K
  T_∞    环境温度       K

h 不是物性参数，而是流场、风速、几何形状与温差的综合结果。自然对流下 h 常在 5～15 W/(m²·K)，强制风冷可到 20～100 W/(m²·K)。仿真中若把 h 设为常数，往往在高温段低估散热，需要按经验关联式修正。

辐射换热遵循斯特藩–玻尔兹曼定律，实际工程表面用灰体修正：

q_rad = ε·σ·(T_s⁴ − T_sur⁴)

  ε      发射率（0～1）
  σ      斯特藩–玻尔兹曼常数 5.67×10⁻⁸ W/(m²·K⁴)
  T_sur  环境辐射温度  K

注意温度必须用开尔文。辐射项与温度的四次方成正比，在 200 ℃ 以上时其贡献迅速上升；对制动盘这类高温部件，忽略辐射会明显高估温升。工程上常把辐射线性化为等效辐射换热系数 h_rad = εσ(T_s²+T_sur²)(T_s+T_sur)，与对流合并成总换热系数 h_total = h_conv + h_rad，便于在仿真中统一施加。

三、边界条件类型与设置方法

热仿真中常见的边界条件有三类，选择哪一类取决于你对外界环境的已知程度：

类型	数学形式	物理含义	制动器典型场景
第一类（Dirichlet）	T = T_w	表面温度已知	与恒温水冷板贴合的面
第二类（Neumann）	−k∂T/∂n = q_w	表面热流已知	摩擦面输入热流、加热膜
第三类（Robin）	−k∂T/∂n = h(T_s−T_∞)	对流/辐射换热	外壳、端盖、线圈外露表面

实际建模中，同一部件不同面往往用不同边界条件：线圈外圆面用 Robin，铁芯与安装座接触面用等效接触热阻，对称面用绝热（q=0）。接触热阻是常被忽略的一项，它由表面粗糙度、接触压力与界面材料决定，量级可达 10⁻⁴～10⁻³ m²·K/W，对局部温升影响显著。

四、热物性参数与温度依赖性

电磁热耦合仿真的非线性主要来自两方面：一是损耗随温度变化（铜电阻率随温度上升），二是热物性随温度变化。若把 k、c_p、ρ 都当作常数，计算会快很多，但在宽温域下误差不可接受。

导热系数 k：金属随温度升高通常缓慢下降；硅钢片叠压方向与叠片方向差异大，需用各向异性 k。
比热容 c_p：随温度升高而增大，决定瞬态温升速率，对短时制动工况尤为关键。
密度 ρ：变化较小，一般可视为常数。
发射率 ε：随温度与表面氧化状态变化，抛光金属低（0.05～0.2），氧化后显著升高（0.6～0.9）。
电阻率 ρ_e：铜约 0.0039 /K 的温度系数，直接改变焦耳损耗，是电磁热双向耦合的核心。

建议在仿真中至少把铜线圈的电阻率、制动盘材料的 k 与 c_p 设为温度函数，其余参数可先取常数做敏感性分析，再决定是否升级为全温度相关。

五、案例：制动盘温升的边界条件配置

以一个电磁制动器制动盘为例，仿真中可这样组织边界与参数：

材料：铸铁制动盘
  k(T)   = 52 − 0.02·(T−20)      W/(m·K)
  c_p(T) = 460 + 0.30·(T−20)     J/(kg·K)
  ρ      = 7200                  kg/m³
  ε      = 0.75（氧化表面）

边界：
  摩擦环面      q_w = 摩擦生热热流（第二类）
  外圆面/端面   h_conv = 12 W/(m²·K)，T_∞ = 25 ℃
  辐射          ε = 0.75，T_sur = 25 ℃
  安装配合面    接触热阻 R_c = 5×10⁻⁴ m²·K/W
  对称面        q = 0（绝热）

求解后应重点检查：最高温度是否出现在摩擦环而非几何中心；辐射项在峰值温度时的占比是否超过 15%；接触热阻是否造成安装侧温度梯度异常。若发现温升远高于实测，优先排查 h 是否偏小、辐射是否被忽略、接触热阻是否漏设。

六、常见误区与复盘

误区一：温度单位混用。辐射计算必须用开尔文，若把 200 ℃ 直接代入 T⁴，结果会严重偏小。

误区二：把 h 当物性常数。h 随流场与温差变化，高温段应改用经验关联式或与 CFD 耦合。

误区三：忽略接触热阻。装配面直接设为理想接触，会导致热量被“顺畅”导出，低估线圈与铁芯温度。

误区四：热物性全取常温值。铜电阻率随温度上升约 40%（100 K 温升），损耗与温升形成正反馈，忽略后瞬态温升曲线会明显偏低。

复盘时可用一张“参数—边界—结果”对照清单：逐项确认 q_v 来源、h 取值依据、ε 与 T_sur、接触热阻、k 与 c_p 的温度函数是否与工况匹配。只有这些输入可信，电磁热耦合的温度场结果才具备工程参考价值。