2. 3D IC封装典型结构与材料体系
3D IC 封装的热结构耦合仿真,第一步不是打开软件画网格,而是把封装体“拆解”成一组在几何、材料和热源上都能独立描述的对象。只有先建立正确的结构—材料映射,后续的热边界条件、材料非线性、应力参考温度才有物理依据。本节围绕 TSV、微凸点、RDL、硅中介层、堆叠芯片、基板与散热盖这条典型链路,梳理各部分的几何角色、材料选择及其热力学参数的温度依赖性。
一、典型结构链路与各层功能
以常见的 2.5D/3D 混合封装为例,从下往上依次是:封装基板(有机层压板或陶瓷)→ 硅中介层(Interposer)→ 微凸点与 RDL 布线层 → 堆叠芯片(Die)→ 热界面材料(TIM)→ 散热盖(Lid/Heat Spreader)。TSV 则垂直贯穿中介层或芯片,把上下表面的信号与电源引出。这条链路中,每一层的厚度跨度极大:RDL 与微凸点在微米量级,硅中介层约 50~200 μm,基板可达数百微米至毫米级,散热盖则在毫米量级。这种跨尺度特征直接决定了仿真中必须采用分层建模或子模型策略。
从热流路径看,芯片产生的焦耳热主要有两条通道:一是经 TIM 与散热盖向上传导至外部散热器;二是经微凸点、TSV、中介层向下传导至基板。两条路径的热阻并联,决定了结温分布。因此,各层材料的导热系数、界面热阻以及温度依赖性,是影响仿真精度的核心参数。
二、关键材料的热力学参数
硅是芯片与中介层的主体材料,室温下导热系数约 148 W/(m·K),热膨胀系数(CTE)约 2.6 ppm/K。铜用于 TSV、RDL 和微凸点,导热系数约 400 W/(m·K),CTE 约 17 ppm/K。两者 CTE 相差近 6 倍,这是 TSV 周围产生热应力的根本原因。焊料(如 SAC305)导热系数约 58 W/(m·K),熔点约 217 ℃,在接近熔点前力学行为显著软化。有机基板(如 FR4 或 BT 树脂)导热系数仅 0.3~0.5 W/(m·K),CTE 在 15~18 ppm/K,且存在玻璃化转变温度 Tg(约 130~180 ℃),超过 Tg 后 CTE 会突增。TIM 的导热系数跨度很大,从普通硅脂的 1~5 W/(m·K) 到金属基 TIM 的 30~80 W/(m·K)。散热盖常用铜或铝,铜的导热系数约 400 W/(m·K),铝约 237 W/(m·K)。
这些参数并非恒定值。硅的导热系数在 300~500 K 范围内会下降约 20%~30%;铜的导热系数随温度升高略有下降,但幅度较小;焊料在接近熔点时导热系数和弹性模量都会明显降低;有机基板超过 Tg 后,CTE 从约 16 ppm/K 跃升至 40 ppm/K 以上,同时弹性模量下降一个数量级。若仿真中把这些参数当作常数,在高温工况下会高估或低估热应力,误差可能达到 30% 以上。
三、温度依赖性的处理策略
在仿真中处理温度依赖性,通常有三种做法。第一种是常数近似,适用于温度变化范围小(如 ΔT < 50 K)且对精度要求不高的场景。第二种是分段线性,把 k(T) 或 CTE(T) 曲线离散成若干温度点,在求解器中以表格形式输入,适用于大多数工程分析。第三种是用户子程序,把材料本构写成温度的函数,适用于焊料蠕变、基板粘弹性等强非线性问题。
以硅中介层为例,若工作温度从 25 ℃ 升到 125 ℃,硅的导热系数从约 148 W/(m·K) 降到约 120 W/(m·K),降幅约 19%。如果仿真中固定用 148,中介层的热阻会被低估,导致结温预测偏低。对于有机基板,若结温超过 Tg,CTE 从 16 ppm/K 跳到 45 ppm/K,热应力会急剧增大,此时必须使用分段线性或子程序描述。
建模提示:在定义材料库时,建议为每种材料至少准备两个温度点的数据(如 25 ℃ 和 125 ℃),并在求解器中开启温度依赖选项。对于焊料和基板,建议准备 25 ℃、Tg、熔点附近三个关键温度点。
四、界面热阻与等效处理
除了体材料参数,界面热阻(TBR)往往被忽略却影响显著。TIM 与芯片、TIM 与散热盖之间的接触热阻,典型值在 1~10 mm²·K/W 量级。在仿真中,界面热阻可以用薄层单元模拟,也可以直接作为接触对的热阻边界条件输入。对于微凸点阵列,由于焊料与 Cu 焊盘之间存在界面,等效导热系数通常低于纯焊料,需要根据凸点直径、间距和高度折算。
一个实用的等效方法是:把微凸点层视为各向异性层,垂直方向导热系数由焊料和空隙面积加权,水平方向则受 Cu 布线影响。这样既保留了主要热路径,又避免了逐个建凸点带来的网格量爆炸。
五、材料参数表(示例)
下表汇总了本节涉及的主要材料在室温附近的典型参数,供建模时参考。表中数值为文献常见范围,实际项目应以供应商数据表或实测为准。
| 材料/部件 | 导热系数 k [W/(m·K)] | CTE [ppm/K] | 关键温度点 |
|---|---|---|---|
| Si(芯片/中介层) | 148(25℃)→120(125℃) | 2.6 | — |
| Cu(TSV/RDL) | 400(25℃)→385(125℃) | 17 | — |
| SAC305 焊料 | 58 | 21 | 熔点 217℃ |
| 有机基板 | 0.3~0.5 | 16(<Tg)→45(>Tg) | Tg 130~180℃ |
| TIM(硅脂) | 1~5 | — | — |
| TIM(金属基) | 30~80 | — | — |
| 散热盖 Cu | 400 | 17 | — |
| 散热盖 Al | 237 | 23 | — |
注意:表中 CTE 为线膨胀系数,仿真中若使用体膨胀系数需乘以 3。焊料的力学参数(弹性模量、泊松比、蠕变参数)随温度变化剧烈,热结构耦合分析中不能只用导热系数和 CTE 两个参数。
建立正确的结构—材料映射后,下一步才是几何清理、网格划分和边界条件设置。本节给出的参数范围和温度依赖处理策略,将直接决定后续热仿真与应力仿真的可信度。建议在项目初期就建立一份材料参数清单,标注数据来源、温度范围和不确定度,避免在结果分析阶段才发现参数缺失或矛盾。