稳态导热 从串联热阻到界面温度
把概念、模型计算与实测证据分开，阅读时留意条件、单位和比较口径。
章节速览11 节 · 展开选择
热膨胀讨论均匀温升后的形变；如果两端温度不同，首先要知道温度怎样沿材料分布。扩散中的“梯度驱动输运”也可帮助理解导热，但热量与物质的变量、单位不同，不能直接混用系数。
本文只处理一维平壁、稳态、无内部热源、各层导热系数恒定、截面积相同的模型。侧面绝热，冷热端表面温度已知，层间先假定理想接触。数值全部为教学合成参数，不对应真实材料牌号或工程设计。
先画边界与热流方向，再建立热阻网络，最后用每层温降重建温度分布。
相同热流率来自串联守恒，相同温差来自共同端点；表面温度不能直接用流体温度替代。
从热流率到热阻
取x从热端指向冷端，傅里叶定律写成Q̇=−kA·dT/dx。k单位W/(m·K)，A单位m²，Q̇单位W。无热源的稳态中，各截面流入与流出相等，因此Q̇沿x不变；同一均匀层内温度梯度恒定。
对厚度L的一层积分，得到Q̇=kA(T左−T右)/L。定义热阻R=L/(kA)，单位K/W，就有温差ΔT=Q̇R。推导和电阻类比可核对MIT热阻网络讲义。类比服务于方程，不代表材料真的存在电流。
Q̇ = (T热 − T冷) / ΣR_i
第i层温降 = Q̇R_i
稳态、无内部热源、侧面绝热、各层恒定k且等面积；给定材料端面温度。R单位K/W，Q̇单位W。对流或接触热阻按后文另加，不能把流体温度直接当作表面温度。来源：文献1。
flowchart LR A[稳态且无内部热源] --> B[各层热流率相同] B --> C[每层温降等于热流率乘热阻] C --> D[各层温降相加] D --> E[总热阻为串联之和]
两层材料手算
设A=0.01m²，热端表面80°C，冷端表面20°C。第一层L₁=0.01m、k₁=1W/(m·K)；第二层L₂=0.02m、k₂=0.5W/(m·K)。
| 量 | 第一层 | 第二层 | 总体 |
|---|---|---|---|
| 热阻K/W | 1 | 4 | 5 |
| 温降K | 12 | 48 | 60 |
| 热流率W | 12 | 12 | 12，不能相加成24 |
| 两侧表面温度°C | 80→68 | 68→20 | 80→20 |
计算为Q̇=60/(1+4)=12W，界面温度80−12×1=68°C。从冷端反算20+12×4也得到68°C。温度差以摄氏度或开尔文表示时数值相同，但这不意味着绝对温度可以在任何公式中互换。
第二层热阻是第一层的四倍，承担的温降也是四倍；不是每层“平均分配30度”。在理想接触界面上温度连续，但两边梯度可以不同，因为k不同。
面积改变时，哪些量改变
把温度分布也画出来
令热端为x=0，两层界面为x=0.01m，冷端为x=0.03m。由于q″=1200W/m²，第一层dT/dx=−1200K/m，第二层dT/dx=−2400K/m。于是分段温度为：
第一层：T(x)=80−1200x，0≤x≤0.01m。
第二层：**T(x)=68−2400(x−0.01)**，0.01≤x≤0.03m；这里x以米计，T以摄氏度计。
代入界面，两段都给出68°C；代入冷端，第二段给出20°C。斜率在界面改变，并不违反温度连续：温度连续约束的是函数值，热流连续约束的是−k·dT/dx。不同导热系数需要不同梯度，才能传递同样的热流密度。
热阻的份额，决定温降的份额
将任意一层的温降Q̇R_i除以总温差Q̇ΣR_i，可得ΔT_i/ΔT总=R_i/ΣR_i。本例第一层承担1/5、第二层承担4/5的总温降。这使结果可以在不重算每个温度的情况下先做量级检查。
但这个关系描述的是热阻份额，不是厚度份额。若两层厚度相同、导热系数不同，温降仍不必相同；若厚度不同但L/k相同，等面积条件下温降反而相同。看图时也不能只凭哪层“画得更宽”判断它的隔热作用。
热流密度q″=Q̇/A，单位W/m²。本例q″=1200W/m²。单位面积热阻R″=L/k，单位m²·K/W；两层合计R″=0.05m²·K/W。因此q″=60/0.05=1200W/m²，再乘面积得到12W。
如果保持厚度、材料和端面温度不变，把面积加倍，总热阻减半，热流率加倍，热流密度不变，界面温度仍为68°C。这里不能把K/W与m²·K/W的热阻直接相加。
自检：只把第二层厚度加倍，界面温度是多少？
第二层热阻变为8K/W，总热阻9K/W，热流率60/9=20/3W。第一层温降20/3K，界面温度80−20/3=220/3°C，约73.333°C。增加冷侧热阻后，界面更接近热端温度。
接触热阻为何允许温度跳变
若界面具有集中接触热阻R_c=0.5K/W，可在两层之间串入它。此时Q̇=60/5.5=120/11W，第一层右表面温度为760/11°C≈69.091°C，第二层左表面温度为700/11°C≈63.636°C。接触处温差为60/11K≈5.455K。
这里的跳变是集中参数模型的表达；没有把界面当作温度既连续又不连续的同一个理想条件。真实接触区可能受粗糙度、压力、填充物与温度影响，R_c需要适用条件下的数据，不能把合成数值当成可直接使用的材料常数。
flowchart LR A[热端表面80度] --> B[第一层热阻1] B --> C[接触热阻0.5] C --> D[第二层热阻4] D --> E[冷端表面20度]
边界条件不一样，网络也要改变
| 给定条件 | 应怎样建模 | 不能省略什么 |
|---|---|---|
| 两端表面温度已知 | 本文串联导热网络 | 不再把它们当作流体远场温度 |
| 两侧流体温度和换热系数已知 | 在导热层外加入1/(hA)的对流热阻 | 表面温度须求解，不等于流体温度 |
| 温度随时间变化 | 加入储热与初始条件，求瞬态 | 不能只凭稳态热阻得出升温时间 |
| 有内部热源 | 使用能量方程及热源项 | 各截面热流率未必相同 |
| 热桥或显著侧向传热 | 考虑并联路径或多维模型 | 不能强行让所有路径热流率相同 |
辐射通常涉及绝对温度的四次方，也不能直接把固定导热热阻当作完整辐射模型。温度变化很大时，k的温度依赖也可能使本文的线性关系不足。
自检：热阻较大，是否就能证明升温更慢？
仅凭稳态热阻不能确定升温过程。还需要密度、比热、几何、初始条件与边界随时间的变化。稳态模型回答最终温差怎样分配，不给出到达该状态的时间。
把流体边界加入同一张网络
前面的80°C和20°C是材料表面温度。若它们改为两侧充分混合流体的温度，表面与流体之间还有换热温差，直接沿用12W就会漏掉两段热阻。
MIT导热与对流联合讲义给出了把两侧换热与固体导热串联的表达式。下面沿用这一边界模型，所有参数仍为自拟。
在指定换热系数h的集中边界模型中，先统一正方向。本例取热侧流体向冷侧流体为正，热侧Q̇=h_h A(T_h−T_s,h)，冷侧Q̇=h_c A(T_s,c−T_c)。每侧沿热流方向的温降与热流率之比均为R_conv=1/(hA)。取热侧h_h=100W/(m²·K)、冷侧h_c=50W/(m²·K)，仍使用A=0.01m²，则两侧对流热阻分别为1K/W、2K/W。
保留两层导热热阻1和4，暂不加接触热阻，总热阻为1+1+4+2=8K/W。热流率Q̇=60/8=7.5W；从热流体80°C出发，逐段减去温降即可恢复各节点：
| 位置或过程 | 热阻K/W | 本段温降K | 下一个节点温度°C |
|---|---|---|---|
| 热侧流体到热表面 | 1 | 7.5 | 72.5 |
| 第一层导热 | 1 | 7.5 | 65 |
| 第二层导热 | 4 | 30 | 35 |
| 冷表面到冷侧流体 | 2 | 15 | 20 |
这次材料的两端表面温度是72.5°C和35°C，不再是80°C和20°C；层间温度为65°C。所有温降仍相加为60K。这里的h是给定的教学参数，并没有通过流体速度、流态或实验推算它；真实换热系数未必恒定，也不能只凭材料名称确定。
练习：为什么加入对流热阻后，不能保留原来的68°C界面温度？
边界已经从材料表面改成流体远场。额外热阻改变了整个串联网络的热流率，也改变两端材料表面温度。新的热流率为7.5W，第一层左表面72.5°C，界面为72.5−7.5×1=65°C。
串联和并联：先看端点，再写公式
如果两层材料沿热流方向前后排列，热量必须依次通过它们，属于串联。若两块材料分别连接同一个热端和冷端，则属于并联：两条路径温差相同，总热流率是各支路之和。这是MIT热阻网络讲义中的串并联对照所强调的结构区别。
设并联两支路的实际支路热阻分别为1K/W和4K/W，两端温差60K。则Q̇₁=60W、Q̇₂=15W，总热流率75W，等效热阻60/75=0.8K/W，也等于1/(1/1+1/4)。
注意，这组支路热阻是另一个已给定的网络例子。若把原来0.01m²的平壁切成两条并联路径，必须根据各自面积重新计算L/(kA)，不能直接把原来的整面积热阻原封不动搬过来。
| 网络 | 相同的量 | 相加的量 | 典型误用 |
|---|---|---|---|
| 无侧向损失的串联 | 各层热流率 | 温降、热阻 | 把每层12W相加成24W |
| 共同冷热端的并联 | 支路温差 | 支路热流率、热导 | 强迫每条路径流过相同热流 |
| 串并联混合 | 逐个子网络判断 | 依连接关系化简 | 只看材料数量，不看连接端点 |
反过来设计厚度：先说明固定条件
假定仍给定两端材料表面80°C与20°C、面积0.01m²，第一层不变，第二层导热系数仍为0.5W/(m·K)。若要求热流率不超过8W，且不考虑接触热阻，必须有R总≥60/8=7.5K/W，因此R₂≥6.5K/W。
由L₂=R₂k₂A得到L₂≥6.5×0.5×0.01=0.0325m，即32.5mm。取等号时，第一层温降8K、界面72°C，第二层温降52K，合计60K。若另外加入确定的接触或表面换热热阻，应先写出新网络再求厚度，不能随意把另一种边界下算出的32.5mm当成普适值。
这只是合成模型的逆向计算，不是保温层选型建议。允许热流、材料参数、温度范围、接触状态和其他失效模式都需在真实任务中另行确定。
练习：如果把两层的导热系数都加倍，界面温度会改变吗？
在原先给定表面温度、无接触热阻的模型中，两层热阻都减半，热流率加倍为24W；每层温降Q̇R_i不变，因此界面仍为68°C。若还存在没有同步缩小的接触或对流热阻，这个结论就未必成立。
怎样判断哪一项最值得继续核查
在固定总温差的串联模型中，Q̇=ΔT/R总。如果只有第i项热阻产生小变化，则相对变化近似为dQ̇/Q̇=−dR_i/R总；若写成该项自身的相对变化，就是−(R_i/R总)·dR_i/R_i。因此同样的相对误差，发生在占总热阻份额较大的项目上时，对热流率的影响通常更大。
本例R₂占总热阻4/5。如果它增加10%，总热阻从5变为5.4，热流率从12降为100/9≈11.111W，实际降低约7.407%；一阶近似给出降低8%。有限变化不等于微分，不能把近似百分比当成精确答案。
这也提醒读者不要只报告许多小数位：若k、接触热阻或边界条件并不精确，算出十位小数也不代表材料温度同样可靠。更有用的是报告输入、假设、主要热阻份额，以及结论对这些条件的敏感性。
一张从题目到结果的核算单
端点是流体还是表面？侧面能否散热？有无热源和储热变化？面积是否相同？哪些界面具有接触热阻？
串联热流率一致吗？温降总和正确吗？界面两侧的温度是否符合接触模型？极端参数下结果是否合理？
先检查单位：米、平方米、导热系数与热阻是否一致。再检查能量：无热源串联模型中的热流率应相同。最后检查温降：所有层和接触处的温降之和，应等于给定的端面温差。
本例依次得到12W、68°C和60K的总温降，三项相互约束。若只得到一个看起来合理的温度，却没检查热流与总温差，仍可能把面积或接触条件放错位置。继续阅读扩散的时间尺度，再对比热膨胀的约束条件。