4、热路模型基础：热阻、热容、热源的定义，与电路模型的类比关系

电机温升计算的工程难点在于：发热源（铜耗、铁耗、机械损耗）在空间上分布不均，散热路径（定子铁心→机壳→空气、转子→气隙→定子、轴承→端盖）呈现典型的多支路并联与串联结构，且绕组绝缘材料的导热系数远低于金属部件，导致电机内部温度梯度显著。若直接求解三维热传导偏微分方程，虽然精度高，但计算量大、参数获取困难，不适合工程快速评估与保护定值整定。热路模型（Thermal Network Model）正是将分布参数的热传导问题，通过集总参数（Lumped Parameter）近似，转化为与电路网络高度同构的“热阻-热容-热源”网络，从而利用成熟的电路分析方法（基尔霍夫定律、节点电压法、一阶RC暂态响应）来求解电机各关键部位的温升过程。本章将系统定义热路三要素，并深入剖析其与电路模型的类比关系，为后续章节的温升计算与热保护动作方程奠定数学与物理基础。

上图展示了热路与电路在拓扑结构和数学描述上的严格对应关系。在电机热分析中，我们通常将绕组铜损、铁心铁损和机械损耗等效为受控电流源（热源），将绝缘材料、气隙、机壳散热表面等效为热阻，将绕组铜块、铁心质量块等效为热容。这种类比的价值在于：可以直接借用电路理论中的节点分析法建立热平衡方程，利用一阶或二阶RC网络的解析解来预测电机在阶跃负载、周期性负载或堵转工况下的温升曲线，从而为热保护继电器（如基于热积累模型的过载保护）提供动作判据。

4.1 热路三要素的物理定义与工程量化

4.1.1 热源（功率损耗）

热源是热路中的激励量，对应电路中的电流源。在电机中，热源并非单一集中点，而是分布在定子绕组、定子铁心、转子导条、永磁体及轴承等区域。工程计算中，通常将总损耗按产生部位归并到对应的集总热容节点上。例如，对于一台异步电机，定子铜耗 \(P_{cus}\) 直接作用于绕组节点，铁耗 \(P_{fe}\) 作用于定子铁心节点，转子铜耗 \(P_{cur}\) 作用于转子节点，机械损耗 \(P_{fw}\) 和杂散损耗 \(P_s\) 则按经验比例分配到机壳和转子节点。热源的时变特性是热保护的核心输入——当电机过载时，铜耗与电流的平方成正比（\(P_{cu}=3I^2R\)），而铁耗近似与电压平方成正比（在变频供电时还与频率相关）。因此，热保护模型必须实时计算当前工况下的损耗值，而非采用额定损耗常数。

4.1.2 热阻（热传导与热对流的阻力）

热阻定义为温差与热流之比，单位 K/W，对应电路中的电阻。它表征热量传递路径上的阻碍能力。电机中的热阻主要分为三类：

• 导热热阻：发生在固体内部，如绕组铜线→绝缘漆→槽绝缘→铁心壁的传导路径。计算公式为 \(R_{cond} = \delta / (\lambda \cdot A)\)，其中 \(\delta\) 为材料厚度，\(\lambda\) 为导热系数，\(A\) 为垂直于热流方向的截面积。绝缘材料的导热系数仅为 0.2~0.3 W/(m·K)，远低于铜的 385 W/(m·K)，因此绝缘层是主要热阻来源。
• 对流热阻：发生在机壳外表面与冷却空气之间，以及转子端部与气隙空气之间。其热阻与表面传热系数 \(h\) 成反比：\(R_{conv} = 1/(h \cdot A)\)。\(h\) 受风速、表面粗糙度、散热片形状影响显著，强迫风冷时 \(h\) 可达 50~100 W/(m²·K)，自然冷却时仅为 5~10 W/(m²·K)。
• 辐射热阻：在电机温升不高（<200°C）时，辐射散热占比通常小于 10%，工程上常忽略或并入对流热阻。

值得注意的是，热阻并非恒定常数。例如，气隙热阻随转速变化（转子旋转增强对流），机壳散热热阻随风机转速变化。在精确热保护模型中，需要将热阻表示为转速或温度的弱函数，或者采用分段线性化处理。

4.1.3 热容（热量的储存能力）

热容定义为单位温升所吸收的热量，单位 J/K，对应电路中的电容。它表征物体储存热能的能力，计算公式为 \(C_{th} = m \cdot c_p\)，其中 \(m\) 为质量，\(c_p\) 为比热容。铜的比热容为 385 J/(kg·K)，铁为 450 J/(kg·K)，绝缘材料约为 1500 J/(kg·K)。在集总参数模型中，我们将绕组铜和绝缘视为一个复合热容节点，其等效比热容按质量加权平均。热容决定了温升的暂态响应速度——热容越大，温升变化越慢，时间常数越大。例如，一台 45 kW 电机的定子绕组热容约为 2~5 kJ/K，对应的热时间常数（绕组到机壳）通常为 5~20 分钟；而机壳到环境的热时间常数可达 30~90 分钟。这种多时间常数特性是热保护算法中必须考虑的关键因素。

物理量	电路符号	热路符号	单位	计算公式	电机中的典型值
热源（损耗）	电流源 I	热流源 P	W	\(P_{cu}=3I^2R\)	额定铜耗 1~5 kW
热阻（传导）	电阻 R	热阻 Rth	K/W	\(R_{cond}=\delta/(\lambda A)\)	绕组-铁心 0.05~0.2 K/W
热阻（对流）	电阻 R	热阻 Rth	K/W	\(R_{conv}=1/(hA)\)	机壳-空气 0.02~0.1 K/W
热容	电容 C	热容 Cth	J/K	\(C_{th}=m \cdot c_p\)	定子铁心 3~10 kJ/K
温度	电压 V	温度 T	°C 或 K	—	绕组热点 120~180°C

4.2 电路与热路的数学同构性及求解方法

从电路理论视角看，热路模型满足与电路完全相同的基尔霍夫定律：在热路节点上，流入的热流之和等于流出的热流之和（对应KCL）；沿任意闭合热路环，温差之和为零（对应KVL）。因此，我们可以直接使用节点电压法建立热平衡方程。以一个简化的单热容模型为例，设电机绕组节点温度为 \(T(t)\)，环境温度为 \(T_a\)，热源功率为 \(P(t)\)，绕组到环境的总热阻为 \(R_{th}\)，绕组热容为 \(C_{th}\)，则热平衡方程为：

Cth · dT(t)/dt = P(t) - (T(t) - Ta) / Rth

该一阶常微分方程的解为：

T(t) = Ta + P·Rth · (1 - e^(-t/τ)) + (T_initial - Ta - P·Rth) · e^(-t/τ)

其中 \(\tau = R_{th} \cdot C_{th}\) 为热时间常数。当电机从冷态启动并施加恒定负载时，\(T_{initial}=T_a\)，温升曲线呈指数上升，最终稳定在 \(T_a + P \cdot R_{th}\)。这一解析解是热保护反时限特性的物理基础——过载倍数越大，\(P\) 越大，达到跳闸温度的时间越短。

对于更精确的多节点模型（如考虑绕组、铁心、机壳三个节点），需要建立二阶或三阶微分方程组。工程上常采用状态空间法或直接数值积分（如四阶龙格-库塔法）求解。在数字式热保护继电器中，通常将微分方程离散化为差分方程，以 100 ms~1 s 的采样周期实时更新温度估计值。离散化公式为：

T[k+1] = T[k] + (Δt/Cth) · [P[k] - (T[k] - Ta)/Rth]

这种递推计算量极小，非常适合在 MCU 或 DSP 中实现。

4.3 案例拆解：基于热路模型的异步电机过载保护定值整定

背景与约束：某化工厂一台 45 kW、380 V、额定电流 85 A 的异步电机，驱动离心泵。电机安装于户外，环境温度最高可达 45°C，采用强迫风冷。用户要求热保护继电器能在 120% 额定负载下长期运行不误动，但在 150% 负载下 10 分钟内必须跳闸，且需防止冷态启动时误跳闸。

初始参数与热路建模：根据电机手册，绕组到环境的总热阻 \(R_{th} = 0.08 K/W\)（含绝缘、铁心、机壳、表面散热），绕组热容 \(C_{th} = 4.5 kJ/K\)，铁心热容 \(C_{fe} = 8 kJ/K\)。额定损耗 \(P_N = 2.8 kW\)（其中铜耗 1.8 kW，铁耗 0.7 kW，机械损耗 0.3 kW）。绕组绝缘等级 F 级，最高允许温度 155°C，设定跳闸温度 \(T_{trip} = 145°C\)（留 10°C 安全裕度）。

分析计算过程：

1. 稳态温升计算：额定负载下，绕组稳态温升 \(\Delta T_N = P_N \cdot R_{th} = 2800 \times 0.08 = 224 K\)。但此值明显过高，说明单热阻模型过于粗糙——实际电机中绕组到铁心的热阻远小于铁心到环境的热阻，绕组温度并非简单等于环境加总损耗乘总热阻。因此需采用双节点模型：绕组节点（含铜耗）与铁心节点（含铁耗），两节点间热阻 \(R_{12}=0.03 K/W\)，铁心到环境热阻 \(R_{2a}=0.06 K/W\)。列写稳态方程：

节点1（绕组）：P_cu = (T1 - T2)/R12
节点2（铁心）：P_fe + (T1 - T2)/R12 = (T2 - Ta)/R2a
代入 P_cu=1800W, P_fe=700W, Ta=45°C：
1800 = (T1 - T2)/0.03  →  T1 - T2 = 54 K
700 + 54/0.03 = (T2 - 45)/0.06  →  700 + 1800 = (T2 - 45)/0.06
T2 = 45 + 2500×0.06 = 195°C  →  T1 = 249°C

计算结果显示绕组温度高达 249°C，远超 F 级绝缘限值。这说明该电机在额定负载下实际无法连续运行——设计裕量不足或冷却条件差。经复核，实际电机额定损耗应为 2.2 kW（铜耗 1.4 kW，铁耗 0.5 kW，机械 0.3 kW），且机壳带有强制风冷，对流热阻可降低 30%。修正后：\(R_{2a}=0.042 K/W\)，重新计算得 \(T_2 = 45 + (500+1400/0.03? )\) 需重新列式：

节点1：1400 = (T1 - T2)/0.03 → T1 - T2 = 42 K
节点2：500 + 42/0.03 = (T2 - 45)/0.042 → 500 + 1400 = (T2 - 45)/0.042
T2 = 45 + 1900×0.042 = 124.8°C → T1 = 166.8°C

绕组温度 166.8°C 仍略高于 155°C 限值，但考虑到实际负载率通常为 85%，且环境温度很少达到 45°C，该电机可勉强运行。热保护整定需基于此模型。

保护定值整定：

• 长期运行区：当负载电流 ≤ 100% 额定（85 A）时，绕组稳态温度约 130°C（按 85% 负载修正），保护不动作。
• 过载跳闸曲线：设定 150% 负载（127.5 A）时，铜耗变为 \(1.4 \times (1.5)^2 = 3.15 kW\)，铁耗近似不变 0.5 kW。代入双节点模型求暂态温升，当绕组温度达到 145°C 时跳闸。利用数值求解（步长 1 s），得到跳闸时间约 420 s（7 分钟），满足用户“10 分钟内”要求。
• 冷态启动闭锁：冷态启动时，绕组初始温度等于环境温度，即使启动电流达 6 倍额定，铜耗为 \(1.4 \times 36 = 50.4 kW\)，但热容吸收热量，温升速率约为 \(P/C_{th} = 50400/4500 = 11.2 K/s\)。若按此速率，2 秒内即达 145°C 跳闸——这显然不合理，因为启动过程仅持续 3~5 秒，且启动后电流迅速下降。因此，热保护模型必须引入“启动判别”逻辑：当检测到启动电流且持续时间 < 10 s 时，将热容模型中的时间常数放大（等效为增加虚拟热容），或直接闭锁跳闸输出，仅进行热积累计算。

结果与复盘：实际投运后，电机在夏季高温时段曾出现保护误动。分析发现，环境温度从 45°C 升至 50°C 时，绕组稳态温度升高 5°C，接近跳闸阈值。解决方案是将跳闸温度设定值改为随环境温度自动调整：\(T_{trip} = 145 - (T_a - 45) \times 1.2\)。同时，将热阻 \(R_{2a}\) 修正为风速的函数（风机转速降低时热阻增大），进一步提高了保护精度。此案例表明，热路模型不仅是理论工具，其参数标定必须结合实测温升试验数据，且保护算法需考虑环境温度补偿和启动工况的特殊处理。

4.4 常见误区与工程注意事项