第2章有限元基础

第二章有限元基础：应力应变基础理论、热-力耦合分析原理、Ansys Workbench/APDL简介

各位同学，大家好。我是你们这门课的老朋友。今天咱们聊聊有限元分析的基础。说实话，这部分内容听起来有点枯燥，但它是后面所有仿真的根基。我个人习惯是，先把这些“硬骨头”啃下来，后面做项目才能心里有底。

咱们做2.5D封装，RDL层动不动就几微米厚，铜和介质的CTE（热膨胀系数）又不匹配。温度一变化，应力就来了。怎么算？怎么评估？这就是本章要解决的问题。

2.1 应力应变基础理论：别怕，就是拉压弯扭

先说说应力。你想想看，一根铜柱子，你用力拉它，它内部会产生抵抗变形的力。这个力除以截面积，就是应力。说白了，应力就是单位面积上的内力。单位是帕斯卡（Pa），我们常用兆帕（MPa）。

应变呢？就是变形量除以原始长度。比如你拉一根1米长的铜线，拉长了0.1毫米，那应变就是0.0001。是个无量纲的量。

在封装里，我们最关心的是热应力。为什么？因为芯片工作会发热，从室温25℃升到125℃，温差100℃。不同材料膨胀不一样，铜膨胀得快，硅膨胀得慢，中间层（比如RDL的介质）夹在中间，就会产生应力。

我记得有一次做项目，客户反馈某个封装在温度循环测试后，RDL层出现了裂纹。我一看仿真结果，最大主应力刚好超过了材料的抗拉强度。嗯，问题就出在这里。所以，理解应力应变，是判断失效的第一步。

核心公式（记住这个就够了）：

σ = E · ε

其中：σ 是应力，E 是杨氏模量（材料刚度），ε 是应变。

对于热应变：ε_thermal = α · ΔT

α 是热膨胀系数，ΔT 是温差。

这里有个坑，我提醒一下。杨氏模量E是温度相关的。温度越高，材料越软，E值会下降。如果你用室温下的E值去算高温下的应力，结果会偏大。我曾经吃过这个亏，后来老老实实查了材料在不同温度下的数据。

2.2 热-力耦合分析原理：先算温度，再算应力

热-力耦合，说白了就是“先热后力”。

第一步，算温度场。芯片发热，热量通过RDL层、基板、焊球传到外界。我们用有限元求解热传导方程，得到每个节点的温度。

第二步，把温度结果作为载荷，施加到结构上。然后算应力。这叫顺序耦合。在Ansys里，就是先做热分析，再做结构分析。

为什么不是直接耦合？因为热和力的时间尺度不一样。热传导慢，应力传播快。在大多数封装仿真中，顺序耦合已经足够精确。我个人的经验是，除非你研究的是高速冲击或者爆炸，否则顺序耦合完全够用。

你可能会问：“那热膨胀系数各向异性怎么办？”好问题。比如RDL层中的铜，电镀后晶粒取向不同，CTE可能在不同方向不一样。这时候，你需要定义正交各向异性材料。在Ansys里，用TB,ANEL命令可以定义。不过，对于大多数2.5D封装，我们假设材料是各向同性的，简化计算。

小技巧：

在设置热-力耦合时，记得把参考温度设对。参考温度通常是材料无应力时的温度，比如电镀时的温度或者固化温度。如果设错了，整个应力结果都会偏移。我一般设为25℃（室温），然后施加温度载荷到工作温度。

2.3 Ansys Workbench/APDL简介：两种武器，各有千秋

Ansys Workbench和APDL，就像自动挡和手动挡的车。Workbench方便，APDL灵活。我两个都用，看情况。

2.3.1 Ansys Workbench：图形化，上手快

Workbench适合做标准流程。你拖拽一个“Static Structural”模块，再拖一个“Thermal”模块，连起来，设置材料、网格、边界条件，求解。可视化做得好，后处理也方便。

对于2.5D封装，我通常用Workbench做参数化扫描。比如改变RDL层的厚度，看应力怎么变。在Workbench里，把厚度设成参数，然后跑几个点，自动生成曲线。这个功能很实用。

但是，Workbench也有缺点。它生成的APDL命令流是“黑盒”，你很难精细控制。比如你想在某个节点上施加一个特殊的力，或者用自定义的失效准则，Workbench就有点力不从心了。

2.3.2 Ansys APDL：命令流，全控制

APDL是Ansys的“灵魂”。所有操作都可以用命令实现。比如：

! 定义材料
MP,EX,1,1.2e5   ! 杨氏模量，单位MPa
MP,ALPX,1,1.7e-5 ! 热膨胀系数，单位1/℃
MP,NUXY,1,0.34   ! 泊松比

! 建立模型
BLOCK,0,10,0,10,0,1  ! 一个10x10x1的块

! 划分网格
ESIZE,0.5
MSHAPE,0,3D
MSHKEY,0
AMESH,ALL

! 施加热载荷
BF,ALL,TEMP,125  ! 所有节点温度设为125℃

! 求解
/SOLU
ANTYPE,0
SOLVE
FINISH

! 后处理
/POST1
PLNSOL,S,EQV,0,1.0  ! 显示等效应力云图

这段代码，我闭着眼睛都能写出来。为什么？因为做封装仿真，很多模型是重复的。我把常用的宏命令存成文件，每次改几个参数就行。效率很高。

APDL的另一个好处是批处理。你可以写一个循环，自动跑几十个工况。比如改变RDL层厚度从1μm到10μm，步长0.5μm，自动求解并提取最大应力。这在Workbench里也能做，但APDL更直接。

注意：

APDL的学习曲线比较陡。刚开始你可能觉得命令太多记不住。我的建议是：不要背，用的时候查。Ansys的帮助文档是最好的老师。我做了十年仿真，也经常查命令的语法。另外，单位制一定要统一。我习惯用mm、N、MPa、℃。如果混用了，结果会差好几个数量级。

2.4 本章知识体系：一张图看懂

下面这张SVG图，是我自己画的。它把本章的核心逻辑串起来了。从应力应变基础，到热-力耦合原理，再到Ansys工具。你保存下来，以后做仿真时对照着看。

这张图里，左边是理论，中间是原理，右边是工具。从上到下，从理论到应用。你每次做仿真前，可以看看这张图，提醒自己：我这一步是在做什么？

2.5 本章小结

好了，这一章的内容就这些。我们讲了应力应变的基本概念，热-力耦合的顺序分析方法，以及Ansys Workbench和APDL的各自特点。说白了，就是为后面的RDL层仿真打好理论基础。

我个人觉得，有限元分析不是“点一下求解”那么简单。你得理解背后的物理。应力从哪里来？温度怎么分布？材料怎么响应？想清楚了，仿真结果才有意义。

下一章，我们会进入实战，用Ansys APDL建立第一个2.5D封装模型。到时候，我会带着你一行一行写代码。别怕，跟着我来。