2. 热仿真物理基础

高热流密度电子冷却仿真的可信度，首先取决于是否把三条基本传热路径写对：热传导负责芯片内部与封装体中的热量扩散，对流负责把热量从固体表面带入冷却介质，辐射则在高温差或真空环境中不可忽略。仿真软件并不会自动替我们判断物理机制，它只是求解我们给出的方程与边界条件。因此，本节的重点不是背公式，而是理解每个方程中哪些量是几何与材料决定的、哪些量是工况决定的，以及它们如何被组织成可计算的热阻网络。

一、三条传热路径的控制方程

热传导由傅里叶定律描述。一维稳态形式为 q = -k·dT/dx，其中 q 是热流密度（W/m²），k 是材料导热系数（W/(m·K)），dT/dx 是温度梯度。负号表示热量沿温度降低方向传递。对三维各向同性材料，稳态导热方程为 ∇·(k∇T) + Q = 0，Q 为体积热源（W/m³）。在芯片有源区，Q 就是功耗密度；在铜层和焊料层，Q 通常为零。瞬态时需加入热容项：ρc_p·∂T/∂t = ∇·(k∇T) + Q，ρ 为密度，c_p 为比热容。这解释了为什么仿真中必须同时输入 k、ρ、c_p，而不仅仅是导热系数。

对流由牛顿冷却定律描述：q = h·(T_s - T_f)，h 为对流换热系数（W/(m²·K)），T_s 为固体表面温度，T_f 为流体参考温度。h 不是材料常数，而是流场、流速、特征尺寸和流体物性的综合结果。自然对流下 h 常在 5～25 W/(m²·K)，强制风冷可达 25～100 W/(m²·K)，水冷微通道可超过 5000 W/(m²·K)。仿真中若把 h 设为常数，等于假设流场均匀；若使用共轭传热求解，则 h 由流场自动算出，但计算量显著增加。

辐射由斯特藩-玻尔兹曼定律描述：q = εσ(T_s⁴ - T_sur⁴)，ε 为发射率，σ = 5.67×10⁻⁸ W/(m²·K⁴)。辐射与绝对温度的四次方成正比，因此在 25℃ 附近、温差 20℃ 时，辐射热流通常只有几 W/m²，远小于对流；但当表面温度升到 150℃ 以上，或处于真空、低气压环境时，辐射可能成为主要散热路径。仿真中需要设置发射率、视角因子和周围环境温度，忽略辐射会低估高温器件的散热能力。

二、热阻网络与热路建模

热阻的定义是 R = ΔT / Q，单位为 K/W 或 ℃/W。它把复杂的传热问题转化为类似电路的串并联网络：热量 Q 相当于电流，温差 ΔT 相当于电压，热阻 R 相当于电阻。对于一维导热，R_cond = L/(k·A)，L 为传热路径长度，A 为截面积。对于对流，R_conv = 1/(h·A)。对于辐射，在温差不大时可线性化为 R_rad = 1/(h_r·A)，其中 h_r = εσ(T_s² + T_sur²)(T_s + T_sur)。

以典型 FCBGA 封装为例，从结到环境的总热阻可写成串联形式：R_ja = R_jc + R_cs + R_sa。若已知芯片功耗 Q = 15 W，环境温度 T_a = 45℃，R_jc = 0.35 K/W，R_cs = 0.25 K/W，R_sa = 1.8 K/W，则 R_ja = 2.4 K/W，结温 T_j = 45 + 15×2.4 = 81℃。这个估算在方案阶段非常有用，但它假设热量全部沿单一路径传递。实际封装中，热量还会通过 PCB 铜层、焊球和侧面路径分流，形成并联热阻。仿真建模时，若只建串联网络，会高估结温；若把 PCB 作为并联路径显式建模，结果更接近真实。

热路建模的关键判断是：哪些路径必须显式建模，哪些可以合并为等效热阻。一般原则是，热流密度高、温差大的路径必须显式；热流密度低、温差小的路径可合并。例如，芯片到散热器的主路径必须保留 TIM 层和均热板；而芯片到 PCB 的次要路径，若 PCB 铜层面积大、导热好，可用一个等效热阻代替。仿真中常见的错误是把 TIM 层厚度设为 0 或忽略接触热阻，这会导致结温被严重低估。

三、电子器件热特性参数与温度限值

器件数据手册中与热仿真直接相关的参数包括：R_jc（结到壳热阻）、R_ja（结到环境热阻）、R_jb（结到板热阻）、Ψ_jt（结到顶面特征参数）和 Ψ_jb（结到板特征参数）。R_jc 由封装材料和结构决定，是器件级参数；R_ja 依赖测试板条件，不能直接用于不同 PCB 设计；Ψ_jt 和 Ψ_jb 是 JEDEC 标准下的特征参数，用于估算结温，但不等同于热阻。仿真中应优先使用 R_jc 和 Ψ_jt，并结合实际散热路径计算。

温度限值方面，硅结温通常不超过 125℃（工业级）或 150℃（汽车级），部分功率器件可到 175℃。焊料层温度受熔点限制，SnAgCu 焊料回流温度约 217℃，长期工作温度建议低于 125℃。PCB 材料 Tg（玻璃化转变温度）常见为 130～180℃，超过 Tg 后机械性能下降。TIM 材料也有工作温度上限，硅脂类通常不超过 150℃。仿真结果判断时，不能只看结温平均值，还要检查热点温度、温度梯度和界面温差。若结温接近限值，应优先降低 R_cs 或提高 R_sa 的散热能力，而不是简单提高风速。

四、仿真中的参数设置与结果判读

在仿真软件中设置材料属性时，导热系数 k 随温度变化不可忽略。硅的 k 从 25℃ 的 148 W/(m·K) 降到 125℃ 的约 100 W/(m·K)；铜的 k 从 401 降到约 380 W/(m·K)。若仿真温度跨度大，应使用温度相关材料曲线。对流边界条件若用常数 h，需根据流速和特征尺寸估算：平板强制对流 h ≈ 10.45 - v + 10√v，v 为风速（m/s）。例如 v = 2 m/s 时，h ≈ 10.45 - 2 + 14.1 = 22.6 W/(m²·K)。辐射设置中，抛光金属 ε 约 0.05，氧化或涂漆表面 ε 约 0.8～0.9。

结果判读时，先检查能量平衡：总发热量是否等于总散热量，误差通常应小于 1%。再检查温度场是否合理：热点是否出现在预期位置，界面温差是否与热阻估算一致。若仿真结温比热阻估算低很多，可能是漏掉了接触热阻或 TIM 层；若高很多，可能是网格过粗或对流系数设置过低。最后，把仿真结果与热阻网络估算、实测数据（如有）交叉验证，才能确认模型可信。