2、热传导基础：傅里叶定律、热阻网络模型、热容与瞬态热响应
在 SiP 封装中，热管理仿真的核心是准确描述热量从芯片结到外壳、再到环境的传递过程。本章从热传导的基本定律出发，建立热阻网络模型，并引入热容概念以分析瞬态热响应，为后续芯片布局与热协同优化奠定理论基础。
2.1 傅里叶定律：热传导的基石
傅里叶定律描述了稳态热传导中热流密度与温度梯度之间的线性关系。对于各向同性材料，其矢量形式为：
q = -k · ∇T其中：
- q — 热流密度 (W/m²)，单位面积上的热流量
- k — 导热系数 (W/(m·K))，材料固有属性
- ∇T — 温度梯度 (K/m)，负号表示热量从高温向低温传递
在一维稳态条件下，通过厚度为 L、截面积为 A 的平板，热流量 Q 可简化为：
Q = k · A · (T₁ - T₂) / L
工程要点： 导热系数 k 是温度的函数。在 SiP 封装中，硅（≈150 W/(m·K)）、铜（≈400 W/(m·K)）、模塑化合物（≈0.7 W/(m·K)）的 k 值差异巨大，导致热流路径严重不均。
2.2 热阻网络模型：从物理场到电路类比
热阻网络模型将热传导问题类比为电路问题，极大简化了 SiP 封装的热分析。热阻 R 定义为温差与热流量的比值：
R = ΔT / Q在封装中，常见的热阻类型包括：
| 热阻类型 | 物理含义 | 典型值范围 (SiP) | 电路类比 |
|---|---|---|---|
| Rjc | 结到外壳热阻 | 0.5 – 5 °C/W | 电阻 |
| Rjb | 结到基板热阻 | 1 – 10 °C/W | 电阻 |
| Rja | 结到环境热阻 | 10 – 50 °C/W | 电阻 |
| Rspread | 扩散热阻（非一维效应） | 0.1 – 2 °C/W | 非线性电阻 |
建模技巧： 在 EDA 工具（如 Ansys Icepak、FloTHERM）中，热阻网络模型通常使用 RC 网络 进行瞬态仿真。对于多芯片 SiP，每个芯片对应一个热节点，通过热阻矩阵耦合。
2.3 热容与瞬态热响应
热容 C 描述了材料储存热量的能力，定义为：
C = m · cp其中 m 为质量 (kg)，cp 为比热容 (J/(kg·K))。在瞬态热分析中，热容导致温度变化滞后于功率变化，形成 RC 时间常数：
τ = R · C对于 SiP 封装，瞬态热响应通常表现为多阶指数衰减，其温度响应函数为：
T(t) = Tamb + P · Rja · (1 - e-t/τ)
常见误区： 在芯片布局优化中，仅关注稳态热阻而忽略热容效应，会导致对功率脉冲（如 CPU 瞬时睿频）的热响应评估不准确。瞬态仿真必须包含热容网络。
2.4 热阻网络模型的工程应用
在 SiP 封装设计中，热阻网络模型用于：
- 结温估算： Tj = Tamb + P · Rja
- 芯片布局评估： 通过热阻矩阵分析芯片间热耦合
- 散热方案选型： 比较不同散热器、TIM 材料的热阻贡献
- 降额设计： 确定最大允许功耗以保证结温低于限值
协同优化实践： 在芯片布局阶段，将高功耗芯片置于靠近散热路径的位置（如基板中心或散热通孔上方），可显著降低 Rjc。热阻网络模型为这种布局决策提供了量化依据。
| 参数 | 符号 | 单位 | 典型 SiP 值 |
|---|---|---|---|
| 导热系数（硅） | kSi | W/(m·K) | 150 |
| 导热系数（模塑化合物） | kM | W/(m·K) | 0.7 |
| 比热容（硅） | cp,Si | J/(kg·K) | 700 |
| 密度（硅） | ρSi | kg/m³ | 2330 |
| 热容（典型芯片） | Cchip | J/K | 0.01 – 0.1 |