A1-01 一维/多维、稳态与瞬态导热的工程建模

课程代码 A1-01 · 板块 A 基础理论与共性技术 / A1 传热与热力学基础 时长 约 3 小时（5 讲 + 1 次母排热阻网络实操） 适合对象 热管理系统 / 仿真工程师入门；想把教科书传热落到工程建模的人 前置 高等数学（偏微分方程基础）、大学物理传热概念 本讲义定位 讲师授课蓝本 / 学员自学讲义，是 A1-01 大纲的完整展开版

引言：教科书里的傅里叶定律，和你要交的那份仿真报告，差在哪

翻开任何一本传热学教材，第一章都是傅里叶定律 q = −k·∇T，干净、优雅、只有三个量。但你真正坐到工位上、要为一段母排或一块冷板出温升结论时，会立刻发现教科书没告诉你最要命的几件事：这个件，我该按一维算、二维算，还是必须建三维？稳态结论够不够，还是得跑瞬态？那些教材里默认忽略、工程里却能让结果偏差数倍的隐藏项——接触热阻、各向异性——到底进不进模型、怎么赋值？

这门课不是重讲一遍传热学。它要解决的是一个建模决策问题：同一个物理定律，落到不同的部件、不同的工况，该用多复杂的模型才算「刚刚好」。模型建复杂了，网格上百万、算一晚上、还未必收敛；建简单了，一维热阻网络三分钟出结果，却在接头拐角处偏差可观，把设计带沟里。真正的功夫，是用几个无量纲判据（Biot 数、傅里叶数、时间常数）在动手之前就判断出：这个问题的「复杂度下限」在哪里。

贯穿全课有一条主线——用最简单、但不至于失真的模型。我们会从傅里叶定律和导热微分方程出发（第 1 讲），先建立「什么时候一维够用」的判据（第 2 讲），再讲清「什么时候必须上多维」以及数值离散怎么做（第 3 讲），然后用时间常数和傅里叶数解决稳态与瞬态的取舍（第 4 讲），最后落到整车部件的真实建模，把那些藏在角落里、却能翻盘的隐藏项一条条揪出来（第 5 讲）。

第 1 讲　傅里叶定律与导热微分方程

一切导热建模都建立在傅里叶定律和由它推出的导热微分方程之上。这一讲不是背公式，而是把每个符号背后的物理讲清楚——因为后面所有判据（一维够不够、稳态还是瞬态）都是从这里推出来的。

1.1 热流密度、导热系数、温度梯度：一个定律，三个量的物理关系

傅里叶定律 q = −k·∇T 说的是一件很朴素的事：热量往温度低的方向流，流得多快，取决于温度变化有多陡、材料导热有多好。拆开看三个量：

• 温度梯度 ∇T 是「驱动力」。没有温差就没有导热，梯度越陡（同样距离内温差越大），热流越强。负号的物理含义很关键：热流方向与温度梯度相反——梯度指向温度升高的方向，而热量自发地从高温流向低温，所以要加负号让 q 指对方向。
• 导热系数 k 是材料的「导热能力」，单位 W/(m·K)。金属高、塑料低、空气更低，跨度能有几个数量级。它决定了同样温差下，铜排和塑料件的传热能力天差地别——这也是为什么母排用铜/铝、隔热件用塑料。
• 热流密度 q 是单位面积的热流，是我们最终要算的结果，也是和边界换热、焦耳产热对接的物理量。

理解这三者的比例关系，比记公式重要：热流正比于梯度、正比于导热系数。这意味着——想减小某处温升，要么增大导热系数（换材料），要么缩短传热路径（在缩短距离的同时压低总温差），要么增大传热面积。整个热管理设计，本质上都在这三个旋钮上做文章。

1.2 直角/柱/球坐标：件形决定用哪套方程

傅里叶定律是矢量式，落到具体坐标系就展开成导热微分方程。为什么要分直角、柱、球三套？因为方程的坐标系应当贴合部件的几何形状，选对了坐标系，问题往往能降维、能拿到解析解：

• 直角坐标：适合板、壁、方块类件形——冷板基板、电池壳壁、母排的矩形截面。绝大多数整车热管理件默认用它。
• 柱坐标：适合管、圆棒、圆柱类——冷却液管路、圆柱电芯、圆形连接器。用柱坐标，径向导热问题能写成一维（只沿半径变化），比硬套直角坐标简洁得多。
• 球坐标：适合球形或近球形——工程中较少，但某些颗粒、传感器探头模型会用到。

选坐标系的判断动作很直接：看这个件在哪个方向上具有对称性。圆管沿周向对称，用柱坐标就能把二维问题降成一维径向问题；这一步降维，往往就是「能不能手算」和「必须上有限元」的分界。

1.3 三类边界条件：模型和真实世界的接口

导热微分方程本身只描述物体内部，它要有解，必须在边界上给定条件。边界条件是模型与外部世界的接口，选错了，内部算得再准也是错的。 三类边界条件对应三种物理情形：

• 第一类（定温 / Dirichlet）：边界温度已知、被固定。适合贴着大热容恒温体、或有强控温的界面。用它要小心——真实工程中「温度被完全钉死」的边界很少，滥用会人为地给模型「喂」进一个它本没有的散热能力。
• 第二类（定热流 / Neumann）：边界热流密度已知。典型是已知焦耳产热注入面、或绝热面（热流为零）。母排通电产热、隔热包覆的外表面，常用这一类。
• 第三类（对流 / Robin）：边界与流体对流换热，热流 = h·(T_壁 − T_流体)。这是整车热管理最常见的边界——冷板对冷却液、壳体对空气、换热器对气流，几乎都是对流边界。这里的 h（对流换热系数）本身就是一门课（见 A1-02 对流换热：自然对流、强制对流与换热系数估算），它的取值精度直接决定第三类边界的可信度。

一个实操提醒：同一个面，用第一类还是第三类，结论可能差很多。把一个实际是对流散热的面，图省事写成定温边界，等于假设它散热能力无限强，温升必然被低估。选边界条件本身就是建模的一部分。

这三类管的都是「模型与外部世界」的接口，物体内部两块固体贴合的那个面不在其中。 母排接头的两段导体、电芯与冷板之间的 TIM 层，都是两个求解域之间的界面，不是外部边界：它既不是定温、也不是定热流，更不是对流——两侧都是固体，没有流体温度可对。在场模型里它由界面条件表达（软件里常见的名字是接触热导 thermal contact conductance、接触热阻 contact resistance 或热接触对；很多求解器的默认设置是「绑定 / 理想接触」，等价于把界面热阻设成零）。它怎么设、取值从哪来，见 5.4。这里先把边界记清楚：把接头面当成第三类对流边界去凑一个虚构的流体温度，等于给它编了一条现实中不存在的散热路径；而直接用默认绑定，则是把这一段热阻整个删掉——两种偷懒都让接头温升朝乐观侧偏，而接头恰恰是母排最热的地方。

1.4 热扩散率 α：温度「扩散」得多快

导热微分方程的通式是 ρc·∂T/∂t = ∇·(k∇T) + q̇：左边是升温存进去的热，右边第一项是导进来的热，第二项 q̇ 是内热源（单位体积产热，W/m³）。整车热管理里的导热问题几乎都带 q̇——母排通电的焦耳热、电芯充放电的产热、电机绕组铜损，都是分布在体内的产热，不是加在表面的边界条件（第 5 讲的母排 I²R 正是这样一个体热源）。同理，k 写在散度算子里面是有讲究的：k 可以是张量，叠片电芯面内与厚度方向差 10–30 倍量级，只有 ∇·(k∇T) 装得下（第 5 讲的各向异性隐藏项）。

只有当无内热源、常物性、各向同性三个条件同时成立，通式才化简成教科书上那个漂亮的形式 ∂T/∂t = α·∇²T；若只是常物性、各向同性而有内热源，对应的化简式是 ∂T/∂t = α·∇²T + q̇/(ρc)。用化简式之前先确认这三条在你的件上成立——拿它去算一段通着几百安培的母排，等于把产热整个丢掉了。

化简式里出现了一个新的组合量——热扩散率 α = k/(ρc)。它和导热系数 k 是两回事，务必分清：

• k 管的是「稳态传多少热」——温度场定下来后，热流多大。
• α 管的是「温度场变化多快」——一个温度扰动在材料里「扩散」、铺开的速度。

α = k/(ρc) 的物理含义可以这样理解：分子 k 是导热能力（传得快），分母 ρc 是体积热容（存得多）。导热能力越强、热容越小，温度变化传得越快。这解释了很多直觉现象：金属 α 大，加热一端另一端很快就感觉到；而蓄热材料 ρc 大、α 小，温度扰动传得慢、被「吸收」掉。

为什么现在就要建立 α 的概念？因为它是第 4 讲所有瞬态判据的基石——傅里叶数 Fo = αt/L² 里的 α 就是它。一个部件该用稳态还是瞬态、瞬态发展到哪一步了，本质上都由 α 和特征尺寸 L 一起决定。把 α 理解透，后面的瞬态判据就是水到渠成。

本讲小结：傅里叶定律 q = −k·∇T 的三个量（梯度是驱动力、k 是材料能力、q 是结果）决定了热管理设计的三个旋钮；坐标系要贴合件形以便降维；三类边界条件是模型与外界的接口，选错则全盘错；热扩散率 α=k/(ρc) 管的是「温度变化多快」，是后面瞬态判据的基石。