第四章：材料属性——硅、铜、介电材料的热力学参数与温度相关属性设置
做RDL层应力仿真，说白了就是在跟材料打交道。你想想看，一个2.5D封装里，硅、铜、介电材料，这三兄弟热胀冷缩的程度完全不一样。温度一变，应力就来了。我刚开始做封装仿真那会儿，总觉得材料参数随便填填就行，结果算出来的应力值跟实测差了30%以上。嗯，从那以后我再也不敢马虎了。
4.1 硅材料：看似简单，实则暗藏玄机
硅作为衬底材料，它的力学行为其实挺复杂的。很多人以为硅是各向同性的，其实不然。单晶硅是各向异性材料，不同晶向上的杨氏模量能差出将近50%。
我个人习惯用以下参数作为基准：
| 参数 | 数值 | 单位 | 备注 |
|---|---|---|---|
| 杨氏模量（<100>方向） | 130 | GPa | 各向异性，需注意方向 |
| 杨氏模量（<110>方向） | 169 | GPa | 比<100>方向高约30% |
| 泊松比 | 0.28 | - | 变化不大，取常数 |
| 热膨胀系数（CTE） | 2.6 | ppm/K | 随温度略有变化 |
| 热导率 | 148 | W/(m·K) | 温度升高会下降 |
我的经验： 如果做全局应力分析，用<100>方向的各向同性近似就够了。但如果你要分析RDL拐角处的局部应力集中，我建议还是用各向异性模型。我在一个项目中偷懒用了各向同性，结果拐角应力低估了15%，后来被工艺工程师追着问了好久。
4.2 铜：温度依赖性是关键
铜在RDL中扮演着导电路径的角色，但它也是最容易出问题的材料。为什么？因为铜的CTE高达17 ppm/K，是硅的6倍多。温度一变化，铜和硅之间就会产生巨大的热失配应力。
更关键的是，铜的力学性能对温度非常敏感。我整理了一份温度相关的参数表，这是我在多个项目中验证过的：
| 温度（°C） | 杨氏模量（GPa） | 屈服强度（MPa） | CTE（ppm/K） |
|---|---|---|---|
| -55 | 128 | 280 | 16.2 |
| 25 | 120 | 250 | 17.0 |
| 125 | 110 | 200 | 18.1 |
| 200 | 102 | 160 | 19.0 |
| 300 | 92 | 120 | 20.2 |
注意： 铜在高温下会发生明显的塑性变形。如果你只做线弹性分析，算出来的应力会偏大。我曾经在一个项目中只用了弹性模型，结果预测的寿命比实际测试短了3倍。后来改用弹塑性模型，才跟实测数据对得上。
在Ansys APDL中设置温度相关的铜材料属性，我一般这样写：
! 定义铜的温度相关杨氏模量
MPTEMP,1,-55,25,125,200,300
MPDATA,EX,1,1,128e9,120e9,110e9,102e9,92e9
! 定义铜的温度相关CTE
MPDATA,ALPX,1,1,16.2e-6,17.0e-6,18.1e-6,19.0e-6,20.2e-6
! 定义双线性各向同性硬化参数
TB,BISO,1,5
TBTEMP,-55
TBDATA,1,280e6,10e9
TBTEMP,25
TBDATA,1,250e6,8e9
TBTEMP,125
TBDATA,1,200e6,6e9
TBTEMP,200
TBDATA,1,160e6,4e9
TBTEMP,300
TBDATA,1,120e6,2e94.3 介电材料：最容易被低估的角色
RDL层中的介电材料，比如聚酰亚胺（PI）或苯并环丁烯（BCB），它们的模量比硅和铜低好几个数量级。很多人觉得它们不重要，随便给个常数就完事了。但我想说，这种想法很危险。
为什么？因为介电材料的CTE通常很高，而且模量随温度变化剧烈。在回流焊过程中（260°C），PI的模量可能降到室温时的十分之一。这意味着什么？意味着高温下介电层几乎失去了承载能力，应力全部由铜承担。
我常用的PI材料参数如下：
| 温度（°C） | 杨氏模量（GPa） | CTE（ppm/K） | 玻璃化转变温度Tg |
|---|---|---|---|
| 25 | 3.5 | 35 | 约280°C |
| 150 | 2.8 | 40 | |
| 260 | 0.4 | 55 | |
| 300 | 0.2 | 60 |
关键点： 介电材料的玻璃化转变温度（Tg）是一个分水岭。低于Tg时，材料处于玻璃态，模量较高；高于Tg时，进入橡胶态，模量骤降。这个转变对RDL层的应力分布影响巨大。我建议你在仿真中一定要考虑Tg效应，否则高温下的应力结果会完全失真。
4.4 温度相关属性的设置策略
好了，现在三种材料都说完了。你可能会问：这么多温度点，我该怎么选？
我的建议是：
- 至少设置5个温度点：-55°C（低温存储）、25°C（室温）、125°C（高温工作）、260°C（回流焊峰值）、300°C（极限温度）
- 关注Tg附近的突变：在Tg前后各加密一个温度点，比如250°C和280°C
- 使用插值函数：Ansys会自动在温度点之间线性插值，但如果你发现线性插值不够准确，可以用MP命令的TBOPT选项指定对数插值
下面这张图是我总结的材料属性设置流程，你可以参考一下：
4.5 避坑指南：我踩过的那些坑
最后，分享几个我亲身经历过的教训：
- 不要忽略蠕变效应：铜在高温下会发生蠕变，尤其是在回流焊过程中。如果你只做静态应力分析，可能会低估RDL的变形量。我建议在高温段（>200°C）考虑蠕变模型。
- 注意介电材料的吸湿性：PI等材料会吸收水分，导致CTE和模量发生变化。如果你做的是湿度-热-力耦合分析，一定要考虑这个因素。
- 验证参数的一致性：不同批次的材料参数可能有差异。我习惯在仿真前先做一组DMA（动态力学分析）测试，确认实际材料的模量和CTE与数据库中的一致。
一个小技巧： 如果你不确定某个温度点的参数，可以用Ansys的MP命令配合TBOPT=LINEAR做线性插值。但记住，在Tg附近，线性插值可能不够准确。我一般会在Tg前后各加一个实测点，确保曲线平滑过渡。
好了，材料属性这块就讲到这里。记住一句话：仿真结果的准确性，90%取决于材料参数的准确性。花时间把材料属性搞扎实了，后面的应力分析才能让人信服。