第2章热机械仿真基础

2. 热机械仿真基础：热传导与热应力基础理论、有限元方法（FEM）在封装中的应用、常用仿真软件对比

各位工程师朋友，欢迎来到第二章。这一章咱们聊聊热机械仿真的基础。说白了，就是搞清楚“热量怎么跑”、“应力怎么来”，以及“我们用什么工具去算”。

我刚开始做封装仿真那会儿，也踩过不少坑。比如有一次，一个硅桥互联的模型，算出来的热应力总是偏大，后来发现是边界条件设错了。嗯，这些经验教训，今天一并分享给大家。

热传导，说白了就是热量在固体内部“搬家”的过程。在2.5D封装里，芯片产生的热量通过硅桥、基板、凸点一层层传出去。这个过程的数学描述，就是大家熟悉的傅里叶定律。

核心公式：

q = -k · ∇T

其中 q 是热流密度，k 是导热系数，∇T 是温度梯度。负号表示热量从高温向低温传递。

我个人习惯把热传导问题分成两类：

在硅桥互联中，我特别关注的是界面热阻。不同材料之间（比如硅和铜凸点）的接触面，其实存在额外的热阻。这个值如果忽略，仿真结果可能差20%以上。我曾经在一个项目中，就是因为没考虑界面热阻，导致温度预测偏低，后来实测才发现问题。

热应力怎么来的？很简单——材料受热会膨胀，但不同材料膨胀系数不一样，互相“拉扯”就产生了应力。

举个例子：硅的CTE（热膨胀系数）大约是2.6 ppm/℃，而铜的CTE是17 ppm/℃。温度变化时，铜凸点膨胀得比硅快，于是界面处就产生了剪切应力。时间长了，可能引发裂纹或疲劳失效。

避坑指南：我曾经在仿真中忽略了CTE随温度的变化，结果算出来的应力值在高温区严重失真。后来改用分段线性CTE模型，才把结果校准回来。所以，CTE的温度依赖性一定要考虑。

热应力的数学描述，核心是热弹性力学方程：

σ = E · (ε - α · ΔT)

其中 σ 是应力，E 是弹性模量，ε 是总应变，α 是热膨胀系数，ΔT 是温度变化量。

你想想看，如果ΔT很大，或者α差异很大，应力就会很大。这就是为什么封装中要尽量选用CTE匹配的材料。

有限元方法，说白了就是把一个连续体切成很多小块（单元），然后对每个小块求解，最后组装起来。为什么封装仿真离不开FEM？因为封装结构太复杂了——有芯片、基板、凸点、硅桥、底部填充胶……解析解根本算不出来。

在2.5D封装中，FEM主要用来做三件事：

关键步骤：

我记得有一次，一个同事用粗网格算硅桥应力，结果最大值只有50 MPa。我让他把网格加密一倍，结果应力飙升到120 MPa。这就是网格敏感性的典型例子。所以，网格收敛性检查是必须的。

市面上主流的封装仿真软件有三款：Ansys、Comsol、Abaqus。我用过它们各自一段时间，说说我的感受。

软件	优势	劣势	适用场景
Ansys	热-结构耦合成熟；Workbench界面友好；网格工具强大	非线性分析稍弱；学习曲线较陡	封装级热应力、焊点疲劳、系统级散热
Comsol	多物理场耦合方便；自定义方程灵活；后处理可视化好	大规模模型计算慢；网格能力一般	多物理场耦合（热电、流固）、研究型分析
Abaqus	非线性分析最强；接触算法优秀；子模型功能强大	热分析模块较弱；前处理繁琐	非线性热应力、裂纹扩展、大变形问题

我的建议：

我个人习惯是：先用Ansys做快速评估，再用Abaqus做精细的非线性分析。这样既保证了效率，又兼顾了精度。

注意：不同软件的材料库和单元类型定义有差异。同一个模型，在Ansys和Abaqus中算出来的结果可能差5-10%。所以，跨软件对比时，务必统一材料属性和边界条件。

下面这张图，是我自己总结的本章知识框架。你可以把它当作一个思维导图来用。

这张图把本章内容串起来了。从热传导和热应力的理论出发，到FEM的应用，再到软件选型，每一步都环环相扣。你可以在后续章节中，随时回来对照这张图，看看自己学到哪一步了。