4、集总参数热网络法：LPTN基本原理、节点划分原则、热网络方程建立
各位工程师朋友，咱们今天聊点实在的。电机热分析，说白了就是搞清楚热量从哪来、往哪去、在哪堵。我做了这么多年电机控制，发现一个有意思的现象：很多搞控制的同事，对电磁算得门清，一提到热就头大。其实没那么玄乎，集总参数热网络法（LPTN）就是帮你把复杂的热问题，简化成一张看得懂的“热路图”。
4.1 LPTN的基本原理：把热问题变成电路问题
先说说基本原理。你想想看，电流在电路里流动，有电阻、有电压、有电流。热量在电机里流动呢？也有类似的东西——热阻、温差、热流。LPTN的核心思想，就是把电机内部连续的温度场，离散成若干个“节点”，每个节点代表一个区域的平均温度。节点之间用热阻连接，热源（比如铜耗、铁耗）就相当于电流源。
我习惯这么理解：热路 = 电路。
- 温度差（ΔT） ↔ 电压差（ΔV）
- 热流（Q） ↔ 电流（I）
- 热阻（Rth） ↔ 电阻（R）
- 热容（Cth） ↔ 电容（C）
这样一来，你就能用基尔霍夫定律来解热问题了。是不是很熟悉？
核心公式：
稳态：ΔT = Q × Rth
瞬态：Cth × dT/dt = Q_in - Q_out
嗯，这里要注意：LPTN是“集总参数”的，意味着我们假设每个节点内部温度是均匀的。这当然是个近似，但对于工程应用来说，精度完全够用。我记得有一次做项目，用LPTN算出来的温升和实测只差了3度，当时心里还挺得意。
4.2 节点划分原则：怎么切才合理？
节点怎么划分？这是LPTN里最考验经验的地方。我刚开始做的时候，恨不得把电机切成一百块，结果模型复杂得要命，算出来还不准。后来慢慢悟出一个道理：节点不是越多越好，关键是要抓住主要矛盾。
我个人习惯遵循以下几个原则：
- 几何对称性原则：电机结构通常是对称的，比如定子齿、定子轭、转子磁钢。利用对称性，可以大幅减少节点数。我一般只建1/4或1/2模型。
- 材料均一性原则：同一种材料、温度梯度不大的区域，可以合并成一个节点。比如整个定子铁芯，如果轴向温度变化不大，就只用一个节点。
- 热源集中原则：热源所在的位置，必须单独设节点。比如绕组铜耗集中在槽内，那绕组必须是一个独立节点。
- 边界条件清晰原则：与外界换热的面（如机壳表面、端盖），要单独设节点，方便施加对流或辐射边界。
举个例子，一个典型的水冷永磁同步电机，我通常会划分这么几个节点：
| 节点编号 | 名称 | 说明 |
|---|---|---|
| 1 | 绕组（铜） | 主要热源，铜耗集中在此 |
| 2 | 定子齿 | 铁耗热源，与绕组有导热 |
| 3 | 定子轭 | 铁耗热源，与机壳接触 |
| 4 | 机壳 | 与水冷道换热 |
| 5 | 转子磁钢 | 涡流损耗热源 |
| 6 | 转轴 | 与轴承、空气换热 |
| 7 | 端盖（前后） | 与空气自然对流 |
你看，7个节点就把一台电机的主要热特性抓住了。我曾经见过有人把绕组分成上中下三层，每个槽单独一个节点，结果模型跑起来慢得要死，精度也没提高多少。说白了，工程问题，够用就好。
4.3 热网络方程建立：从图到数学
节点划好了，接下来就是建立方程。这一步其实很机械，但容易出错。我建议你画一张热网络图，把节点、热阻、热源、热容都标清楚。然后对着图写方程，不容易漏。
咱们还是以上面那个7节点模型为例。先画一张热网络图：
有了这张图，写方程就简单了。对每个节点，应用基尔霍夫电流定律：流入节点的热流之和 = 流出节点的热流之和。对于瞬态问题，还要加上热容的储能项。
以“绕组节点”为例，它的热平衡方程是这样的：
// 绕组节点的瞬态热平衡方程
C_w * dT_w/dt = Q_cu - (T_w - T_t)/R_wt - (T_w - T_r)/R_gap
其中：
C_w = 绕组热容
T_w = 绕组温度
T_t = 定子齿温度
T_r = 转子温度
Q_cu = 铜耗（热源）
R_wt = 绕组到定子齿的热阻
R_gap = 绕组到转子的气隙热阻把所有节点的方程写出来，就得到一个微分方程组。写成矩阵形式：
[C] * d{T}/dt = [G] * {T} + {Q}
其中：
[C] = 热容矩阵（对角阵）
[G] = 热导矩阵（由热阻的倒数构成）
{T} = 温度向量
{Q} = 热源向量我的小技巧：写方程的时候，先写稳态的（去掉热容项），再扩展成瞬态的。这样不容易乱。另外，热阻的取值很关键，我一般会先查材料手册，再用实验数据校准一次。尤其是气隙热阻，它和转速、气隙长度都有关系，千万别随便估一个值。
避坑指南：我曾经犯过一个错误——把定子轭和机壳之间的接触热阻忽略了。结果算出来的机壳温度比实测低了8度。后来才发现，定子铁芯和机壳之间是过盈配合，但微观上还是有空气间隙，这个接触热阻不能省。建议你查一下“Cooper-Mikic-Yovanovich”模型，专门算接触热阻的。
最后，方程建好了，怎么解？稳态问题直接求逆矩阵就行。瞬态问题，我习惯用四阶龙格-库塔法（RK4），步长取0.1秒，精度和速度都还行。当然，你要是用Matlab或Python，直接调用ode求解器也行。
好了，LPTN的基本原理、节点划分、方程建立，咱们就聊到这儿。说白了，这就是个“化整为零、以简驭繁”的思路。你只要把电机拆成几块，搞清楚每块之间的热关系，剩下的就是解方程了。下次你遇到电机过热的问题，不妨先画一张热网络图，说不定问题就一目了然了。