A1-06 热阻网络法与集总参数建模入门

课程代码 A1-06 · 板块 A 基础理论与共性技术 / A1 传热与热力学基础 时长 约 3.5 小时（5 讲 + 1 次建模作业） 适合对象 仿真 / 系统工程师；做实时热模型、BMS/MCU 温度估计的人 前置 A1-01 导热的工程建模；A1-02 对流换热 h 估算；A1-03 辐射换热与发射率（第 1 讲的 h_r 线性化直接用它的结论） 本讲义定位 讲师授课蓝本 / 学员自学讲义，是 A1-06 大纲的完整展开版

引言：CFD 太慢的时候，你还有一张热阻网络

方案阶段的热管理工程师常常处在一个尴尬的位置：架构还在天天变，可总师会上就要一个"电机大概能热到多少度""电池包这么排会不会超温"的答案。这时候你去搭 CFD、划网格、跑瞬态，一轮下来两三天，等结果出来架构又改了。热阻网络（RC 模型）就是为这个场景而生的——把系统抽象成几个节点、几根热阻、几个热容，几分钟就能给出系统级的温度估算。它换来的速度不是免费的：它抹掉了空间细节。但只要你知道它在什么条件下可信、什么时候会骗你，它就是概念阶段最锋利的一把快刀。

更重要的是，热阻网络不只是"草稿工具"。因为它足够轻，它能被降阶后塞进 BMS 或电机控制器里实时运行，去估算那些你根本装不了传感器的地方的温度——比如电机绕组最深处的热点。嵌入式的测温点永远是滞后的、是"事后诸葛亮"；而一个标定好的 RC 模型，能在温度还没到之前就预测到它，让控制器提前降额。

这门课的主线，就是把这件事从头讲透：怎么把导热、对流、辐射三种传热统统换成热阻；怎么把热容、热源、边界拼成一张网络；怎么判断"单节点集总"这个偷懒假设什么时候成立、什么时候会翻车；怎么让网络跑瞬态；最后怎么用试验或 CFD 数据把这张"快模型"反标定得可信。贯穿全课的一个心法是：热阻网络的价值不在于它多精确，而在于它用最少的自由度抓住了主导热流——抓对了，几个参数就够；抓错了，再多节点也是错的。

第 1 讲　三种传热的热阻表达

要搭热阻网络，第一步是把三种传热机理——导热、对流、辐射——全都翻译成同一种语言：热阻。翻译成功的标志是，每一段传热都能写成"温差 = 热流 × 热阻"的欧姆定律形式，这样它们才能像电阻一样串起来、并起来。

1.1 导热热阻 R = L/(kA)：最"干净"的一段

导热热阻直接从傅里叶定律重排而来。一维稳态导热下，热流 Q = kA·ΔT/L，把它写成 ΔT = Q·(L/kA)，括号里那一坨就是导热热阻 R_cond = L/(kA)。它和欧姆定律 V = I·R 一模一样：温差对应电压，热流对应电流，L/(kA) 对应电阻。

理解每个量的角色很重要。L 是热流走过的路径长度，A 是垂直于热流方向的截面积，k 是材料导热系数。路径越长、材料越不导热、截面越窄，热阻越大。这三者里 L 和 A 是几何、由结构决定，k 是材料属性、相对确定，所以导热热阻是三种热阻里最"干净"、最可算的一段——参数不确定性最小。

但它有一个隐藏前提：热流真的是一维的、截面真的是均匀的。当热量会横向铺开（比如一小块功率模块贴在一大片散热板、一片翅片基板上），热流就不再是一根等截面的"管道"，这时 A 到底取哪个面，直接决定你算出来的热阻是偏大还是偏小——这个口径不交代清楚，连误差方向都是不确定的。

按本节的定义（A 是垂直于热流方向的截面积，对一块散热板就是整块板的截面），这一段真实热阻其实是两段串联：

R_真实 = L/(k·A_板截面) + R_扩散

其中 R_扩散 是扩散热阻（spreading resistance；热流反过来从大截面收拢到一个小接触面时，同一个东西叫收缩热阻）。它恒为正，只有当热源铺满整个截面、真的没有横向铺开时才为零。于是两种取法各自给出一个方向确定的界：

• A 取整块板的截面 ⇒ 偏小。它把 R_扩散 整段丢掉了，给出的是真实值的下界；
• A 取热源投影面积 ⇒ 偏大。它等于假设热量被关在热源正下方那根柱子里、一点也不横向铺开；而横向铺开出去的每一条路径都是并联支路，只会让热阻更小，所以这个值是真实值的上界。

⇒ 落地动作：两个都算，当上下界夹逼用。两者接近（热源面积与板截面同量级），说明扩散不重要，取哪个都行；两者差出几倍（小热源 + 大散热板，功率模块贴散热器就是典型），说明 R_扩散 才是这一段的主导项，必须把它显式写成网络里的一段——查形状因子 / 扩散热阻关联式，或干脆把这块板沿横向离散成几个节点（多维导热什么时候必须做、怎么离散，见前置课 A1-01 一维/多维、稳态与瞬态导热的工程建模）。⛔ 尤其不要只算 L/(k·A_板截面) 就以为自己"算保守了"：它是下界，方向是低估热阻、进而低估热源温升——照它放行的散热方案，实车上可能是超温的。所以用这个公式前先问两句：这段热流是被约束在一根"管道"里的，还是会自由扩散的？如果会扩散，我这个 A 取的是哪一个面、它把结果推向哪一边？

1.2 对流热阻 R = 1/(hA)：网络里最不确定的一段

对流热阻从牛顿冷却定律来：Q = hA·ΔT，重排即 R_conv = 1/(hA)。h 是对流换热系数，A 是换热的润湿面积。形式上它比导热热阻还简单，但工程上它恰恰是整张网络里最不确定的一段——因为 h 本身不是常数。

h 取决于流态（层流还是湍流）、流速、流道几何、流体物性，甚至壁面温度。同一个换热面，流速翻一倍，h 可能变化几十个百分点；从自然对流切到强制对流，h 可能差一个量级。这正是前置课 A1-02 对流换热：自然对流、强制对流与换热系数估算 花整门课去估 h 的原因。放到热阻网络里，这意味着：你的模型精度往往被 h 这一段卡住。如果标定后模型和实测对不上，第一个该怀疑的就是对流热阻——它是最该投入标定精力、也最值得去查手册关联式或做试验回归的地方。而导热那一段，通常不用你操心。

1.3 辐射线性化热阻 R = 1/(h_r A)：把 T⁴ 掰直

辐射是三种传热里唯一"天生非线性"的：斯特藩-玻尔兹曼定律里辐射热流正比于 T⁴ 的差，Q ∝ εσ(T₁⁴ − T₂⁴)。这个 T⁴ 让网络不再是线性的——欧姆定律要求热流正比于温差的一次方，塞进一个四次方项，手算的串并联、叠加原理就全失效了，只解线性方程组的求解器也解不动。⚠ 但注意，这说的是求解手段的限制，不是"辐射进不了热网络"：支持非线性迭代的求解器可以直接带着 T⁴ 解（见第 4 讲 4.3）。

解决办法是线性化：把 (T₁⁴ − T₂⁴) 做因式分解成 (T₁ − T₂)(T₁ + T₂)(T₁² + T₂²)，于是可以定义一个"等效辐射换热系数" h_r = εσ(T₁+T₂)(T₁²+T₂²)，让辐射热流写成 Q = h_r·A·(T₁ − T₂)——形式上和对流一模一样，就能塞进网络了。代价是：h_r 里含温度，它只在某个工作点附近成立。温度变化不大时（围绕一个稳态波动），线性化误差可以忽略；温度大范围扫动时，必须重新计算 h_r，否则辐射那一段就算错了。在多数液冷 TMS 里辐射相对对流是小项，但在机舱、电机端部这类高温、低风速、大温差的场合，辐射不容忽视——该建就得建，只是记得先掰直。

除了"只在工作点附近成立"，这个式子还有两条前提必须跟着它一起写在纸面上，漏掉任何一条，辐射那一段都不可信：

• T 一律取绝对温标（开尔文 K），绝不能代摄氏度。斯特藩-玻尔兹曼定律本身就定义在绝对温标上，T⁴ 对温标平移不是线性的，代完摄氏度没有任何事后补救的系数。差多少？以 T₁=100 ℃、T₂=40 ℃、ε≈0.9 为例：取 K（373 K / 313 K）算出的 h_r 在 8 W/(m²·K) 量级；若把 100 和 40 直接代进去，(T₁+T₂)(T₁²+T₂²) 这一坨小了约两个数量级，h_r 只剩 0.1 W/(m²·K) 不到。危险之处在于这个错不会报错，只会让辐射悄悄消失——你会得到一个"辐射果然可以忽略"的结论，并且理直气壮地把它从网络里删掉。
• ε 在这里是"有效发射率"，它带着几何前提。h_r = εσ(T₁+T₂)(T₁²+T₂²) 这个最简形式对应的几何是小物体被大空腔包围（物体发出去的辐射几乎不返回它自己），此时有效发射率恰好等于物体自身的发射率 ε，式子里才只出现一个 ε。换一种几何就要换 ε：两个大平行面之间要用 ε_eff = 1/(1/ε₁ + 1/ε₂ − 1)；一般的两个表面之间还要再乘一个角系数 F₁₂（不是发出去的辐射都打到对面）。这个差别不是小数点后的事：两平行面发射率同为 ε 时 ε_eff = ε/(2−ε)，ε=0.9 的黑漆面只差约一成，而 ε=0.05 的抛光铝箔差了近一倍——低发射率箔恰恰是隔热屏最常用的东西。所以落笔 h_r 之前先问一句：这两个面到底是什么几何关系？（发射率、角系数与灰体假设的完整讲法见前置课 A1-03 辐射换热与整车热害中的辐射屏蔽）

1.4 接触热阻与界面材料：最容易被漏掉的大电阻

前三种热阻讲的是"热量在一段介质里怎么走"，但真实系统里热量还要跨过两个固体的接触面。理想中两个贴合的固体应该无缝传热，实际上微观看两个面只是零星几个凸点在接触，中间是导热极差的空气间隙。这道"接触热阻 R_c"经常是整条热流路径上最大、又最容易被忽略的一段——比如功率模块贴到散热器上，芯片到冷板的一路热阻里，接触界面可能占了相当比重。

关键认知有两条。其一，接触热阻算不出来——它不由材料体导热系数决定，而由表面粗糙度、平面度、装配压力、有没有涂界面材料（TIM，thermal interface material）共同决定，只能实测或用供应商数据。其二，它是可设计的：涂上导热硅脂、导热垫、相变材料填满那些空气间隙，接触热阻能大幅下降；加大装配压力也有效。所以在网络里，接触面必须作为一段独立热阻显式画出来，而不能想当然并进相邻的导热段——漏掉它，你的模型会系统性地低估温升。

本讲小结：三种传热都能翻译成热阻——导热 R=L/(kA) 最干净、对流 R=1/(hA) 最不确定、辐射要先用 h_r 线性化才能进网络；此外别忘了接触面那段隐藏的大电阻。翻译成功的标志是每段都写成"温差=热流×热阻"，能像电阻一样串并联。