3、仿真基础理论：热-力耦合分析、材料本构模型、接触与摩擦模型
各位工程师朋友，大家好。这一章我们聊聊仿真的“地基”——基础理论。说实话，很多人做仿真喜欢直接上手调参数，结果算出来不对，又找不到原因。我个人习惯是，先把理论吃透，再动手。这样出了问题，你心里有数，知道该往哪个方向排查。
TCB工艺的仿真，说白了就是一场“热”与“力”的博弈。芯片被加热，材料膨胀；压头下压，产生应力。这两件事是同时发生、互相影响的。所以，我们得用热-力耦合分析来处理。另外，材料怎么变形、界面怎么摩擦，这些都需要用数学模型来描述。这就是本构模型和接触摩擦模型要干的事。
3.1 热-力耦合分析：温度与应力的“二人转”
TCB工艺中，温度变化是应力的主要来源。芯片从室温加热到300多度，再冷却下来。这期间，硅芯片、铜凸点、基板、底部填充胶，每种材料的热膨胀系数（CTE）都不一样。你想想看，它们挤在一起，谁也不让谁，应力就这么产生了。
热-力耦合分析，就是要把温度场和应力场结合起来算。主要有两种做法：
- 顺序耦合：先算温度场，再把温度结果作为载荷去算应力场。优点是计算快，适合温度变化不剧烈的情况。我早期做封装级仿真时常用这个。
- 直接耦合：温度和应力同时求解。精度高，但计算量大。TCB工艺中，压头下压速度很快，局部温度会因塑性功生热而升高，这时候就得用直接耦合。
我记得有一次做项目，客户反映芯片边缘有裂纹。我一看仿真结果，温度场分布不均匀，芯片边缘温度比中心低了20多度。这就导致边缘区域收缩更快，产生了额外的拉伸应力。后来我们调整了加热台的均温性，问题就解决了。所以，温度场的准确性直接决定了应力结果的可靠性。
3.2 材料本构模型：给材料“画像”
材料本构模型，说白了就是描述材料在受力时怎么变形的数学公式。TCB工艺涉及的材料种类多，而且都在高温下工作，所以本构模型的选择很关键。
3.2.1 弹性模型
这是最基础的模型，适用于小变形、可恢复的情况。比如硅芯片在TCB过程中，变形量很小，用弹性模型就够了。公式就是大家熟悉的胡克定律：σ = E·ε。但要注意，杨氏模量E是随温度变化的。温度越高，E值越小，材料越“软”。
3.2.2 弹塑性模型
铜凸点在TCB过程中会发生塑性变形，也就是永久变形。这时候就需要弹塑性模型。常用的有双线性模型，用屈服应力和切线模量来描述。我建议在TCB仿真中，一定要考虑铜的随动硬化，因为反复加热冷却会导致包辛格效应。
3.2.3 蠕变与粘弹性模型
焊料在高温下会“蠕变”——即使应力不变，变形也会随时间增加。TCB工艺中，焊料在熔点附近停留时间虽然短，但蠕变效应不可忽略。常用的有Norton蠕变模型：ε̇ = A·σⁿ·exp(-Q/RT)。
另外，底部填充胶是高分子材料，具有粘弹性。我习惯用广义Maxwell模型来描述它的应力松弛行为。嗯，这里要注意，粘弹性参数对温度非常敏感，一定要用DMA（动态力学分析）测试数据来拟合。
3.2.4 Anand模型（焊料专用）
对于无铅焊料（如SAC305），我个人强烈推荐使用Anand模型。它把粘塑性和蠕变统一在一个框架里，能很好地描述焊料从弹性到塑性再到蠕变的完整过程。Anand模型有9个参数，需要做不同温度和应变率下的拉伸试验来标定。
3.3 接触与摩擦模型：界面上的“博弈”
TCB工艺中，接触无处不在：压头与芯片、芯片与基板、凸点与焊盘……这些界面上的力学行为，直接影响键合质量。
3.3.1 接触类型
- 硬接触：接触面一旦闭合，就不能穿透。适合压头与芯片的接触。
- 软接触：允许一定的“弹性穿透”，通过罚函数来控制。适合凸点与焊盘的接触，因为凸点本身会变形。
3.3.2 摩擦模型
最常用的是库仑摩擦模型：τ = μ·p，其中τ是摩擦应力，μ是摩擦系数，p是接触压力。TCB工艺中，摩擦系数一般在0.1~0.3之间。但要注意，摩擦系数随温度升高会变化。我记得有一次仿真结果总是不收敛，后来发现是摩擦系数设成了常数，而实际高温下摩擦系数降低了30%。
3.3.3 接触热阻
接触界面存在热阻，会影响热量传递。TCB工艺中，压头与芯片之间的接触热阻很关键。如果热阻太大，芯片升温慢，工艺时间就会延长。我一般用Cooper-Mikic-Yovanovich (CMY) 模型来计算接触热阻，它考虑了表面粗糙度、接触压力和材料硬度。
3.4 小结：理论是“内功”，仿真才是“招式”
好了，这一章的内容就这些。热-力耦合、材料本构、接触摩擦，这三块是TCB仿真的“内功心法”。你把这些理论吃透了，后面调参数、看结果，心里就有底了。
我个人觉得，仿真不是“算出来就行”，而是要“算得对、算得准”。理论功底越扎实，你越能判断仿真结果是否合理。遇到异常结果，也能快速定位问题。
下一章，我们会把这些理论落实到具体的仿真软件操作中。到时候，我会带着大家一步步搭建TCB工艺的仿真模型。敬请期待。