2. 封装结构中的材料行为

封装翘曲与残余应力仿真能否给出可信结果，首先取决于材料模型是否选对、参数是否给对。封装体是一个多材料叠层结构：硅芯片、塑封料（EMC）、有机基板、铜布线/铜柱、焊料凸点/焊球，各自在温度循环、固化过程和机械约束下表现出完全不同的力学行为。若把塑封料当成线弹性材料、把焊料当成弹性体、忽略固化收缩，仿真得到的翘曲方向和焊点应力往往与实测相差数倍。本节按“材料—主导行为—关键参数—在仿真中的处理方式”这条线索，逐一梳理五类关键材料。

一、硅：近似线弹性的热弹性基准

硅在封装温度区间（约 −55 ℃～260 ℃）内不发生明显屈服和蠕变，是最“干净”的材料，通常按各向异性或各向同性线弹性处理。其关键参数为：弹性模量约 130～170 GPa（各向异性时 E₁₀₀≈130 GPa、E₁₁₁≈188 GPa），泊松比约 0.22～0.28，热膨胀系数 CTE 约 2.6～3.3 ppm/℃（随温度略升）。硅的 CTE 远小于塑封料和基板，这正是翘曲与界面剪应力的主要来源。

仿真中硅一般设为线弹性 + 温度相关 CTE，不需要塑性或粘弹性。需要注意的是：硅片减薄到 50～100 μm 后，其弯曲刚度显著下降，芯片自身变形会明显增大，因此厚度必须按实际值建模，不能沿用裸片标准厚度。此外，硅的断裂属于脆性断裂，若需评估开裂风险，应补充最大主应力或断裂力学判据，而不是看塑性应变。

二、塑封料：粘弹性 + 固化收缩 + 温度相关

塑封料是封装中行为最复杂的材料，也是翘曲仿真的核心难点。它同时具备四种特征：

• 温度相关弹性：玻璃化转变温度 T_g 以下模量约 15～25 GPa，T_g 以上骤降到 0.1～1 GPa。因此必须给出模量—温度曲线，而不是单一模量。
• 粘弹性：在 T_g 附近及以上的应力松弛和蠕变不可忽略，常用 Prony 级数（广义 Maxwell 模型）描述剪切松弛模量。
• 固化收缩：环氧树脂交联固化时体积收缩约 0.2%～0.5%，是脱模后室温翘曲的重要来源，必须在固化阶段以化学收缩应变或等效热应变施加。
• 温度相关 CTE：T_g 以下约 8～12 ppm/℃，T_g 以上约 30～45 ppm/℃，呈明显双段折线。

在仿真中，塑封料通常采用“热-化学-力学”耦合：先用固化动力学方程（如 Kamal 模型）计算固化度 α(t,T)，再由 α 得到化学收缩应变和瞬时模量，最后叠加粘弹性松弛。若只做温度循环分析，可简化为温度相关线弹性 + 粘弹性，但固化收缩必须单独作为初始应变施加，否则室温翘曲会被低估。

三、有机基板：各向异性热弹性 + 吸湿

基板（如 BT、FR-4）是纤维增强复合材料，面内（x-y）与厚度方向（z）性能差异显著：面内 CTE 约 12～18 ppm/℃，厚度方向可达 40～70 ppm/℃；面内模量约 20～30 GPa，厚度方向约 8～12 GPa。仿真中应至少按正交各向异性输入 E_x、E_y、E_z、ν 和三个方向的 CTE。

此外，基板在回流焊高温下会吸湿膨胀，并存在 T_g（约 170～190 ℃）附近的模量下降。对于翘曲分析，基板厚度方向 CTE 大、面内 CTE 小，会与铜层形成明显的“双金属效应”，是封装四角翘曲的重要贡献项。若忽略各向异性，翘曲形态可能从“马鞍形”被错误预测为“圆柱形”。

四、铜：弹塑性 + 温度相关

铜布线、铜柱、铜箔在封装中承担导电与部分承载功能。其弹性模量约 110～130 GPa，泊松比 0.34，CTE 约 16～17 ppm/℃。铜在回流焊和温度循环中会进入塑性：屈服强度随温度升高而下降，室温约 200～300 MPa，200 ℃ 时可能降至 100 MPa 以下，且存在加工硬化。

仿真中铜通常采用双线性或非线性等向硬化塑性模型，输入弹性模量、泊松比、屈服强度和切线模量（或真应力—塑性应变曲线）。铜的塑性变形会吸收部分热失配应变，降低界面应力峰值，但会留下残余塑性应变，导致封装体在回到室温后出现不可恢复的翘曲偏移。因此铜的塑性参数对“残余应力”结果影响很大，不能只给弹性模量。

五、焊料：粘塑性 + 蠕变 + 温度相关

焊料（SAC305、SnPb 等）在室温下即处于同系温度 0.5～0.7 T_m 以上，蠕变和应力松弛非常显著。其行为通常用 Anand 粘塑性模型或 Garofalo 双曲正弦蠕变模型描述，关键参数包括：弹性模量（约 30～50 GPa，随温度下降）、屈服强度（约 30～45 MPa）、蠕变参数 A、n、Q/R 等。

焊点在温度循环中承受剪切与拉伸组合载荷，塑性应变和蠕变应变在角点处集中，是疲劳裂纹的起源。仿真中若把焊料简化为线弹性，会高估应力峰值、低估累积非弹性应变，导致寿命预测偏差。正确做法是采用粘塑性本构，并输出等效塑性应变、蠕变应变和应变能密度作为疲劳判据输入。

六、材料参数获取与仿真落地要点

上述参数并非全部来自手册。实际项目中，硅和铜的弹性参数可直接查表，但塑封料的粘弹性 Prony 级数、固化收缩率、焊料的 Anand 参数通常需要 DMA、TMA、动态力学分析和拉伸/蠕变试验拟合。若没有试验数据，可先用供应商 datasheet 中的模量—温度曲线和 CTE 双段值做温度相关线弹性分析，但必须清楚：这种简化会低估塑封料在 T_g 以上的松弛，从而高估翘曲幅值。

在有限元实现层面，建议按以下顺序检查材料设置：

1. 确认每个部件的材料模型类型（线弹性 / 正交各向异性 / 粘弹性 / 塑性 / 粘塑性）与物理行为匹配。
2. 确认所有温度相关参数以表格或曲线形式输入，而不是取室温单值。
3. 确认固化收缩以初始应变或等效热应变方式施加，并检查参考温度设置是否与固化结束温度一致。
4. 确认焊料和铜的塑性/蠕变参数单位一致（MPa、s、K），避免量纲错误导致非弹性应变异常。
5. 用简单叠层梁解析解或文献算例校核一次翘曲量级，再进入全模型分析。

总结来看，硅提供热失配的“刚性基准”，塑封料贡献粘弹松弛与固化收缩，基板引入各向异性和吸湿，铜通过塑性吸收并锁定部分应变，焊料则以蠕变主导疲劳损伤。五类材料的本构选择共同决定了翘曲形态和残余应力分布，任何一项被简化成线弹性，都会使仿真结果偏离真实封装行为。