第4章固化反应动力学

4. 固化反应动力学：固化度概念、Kamal模型、固化放热对温度场的影响、固化收缩

各位工程师朋友，大家好。今天我们来聊聊固化反应动力学。说实话，这个主题在封装仿真里，属于那种「看着不起眼，但一错就全盘皆输」的关键环节。我刚开始做模流分析那会儿，就吃过固化参数设错的亏，后来花了整整两周才把翘曲问题调回来。嗯，咱们今天就把这块硬骨头啃下来。

固化度，简单说就是环氧树脂从液态变成固态的「完成百分比」。0% 是纯液态，100% 是完全固化。但在实际工程中，我们很少能跑到 100%，一般到 95% 以上就算固化了。

我个人习惯把固化度理解成「化学反应的进度条」。你想想看，EMC 里的环氧树脂和固化剂，在高温下发生交联反应，分子链越连越长，最终形成三维网络结构。这个网络越密，材料就越硬，模量就越高。

关键点：固化度直接影响材料的力学性能。固化度 80% 和 95% 的 EMC，模量可能差 2 倍以上。我在做 QFN 封装仿真时，发现翘曲预测不准，最后查出来是固化度模型参数没校准。

固化度的测量方法主要有两种：

在仿真中，我们通常用 DSC 数据来标定固化动力学模型。这里有个坑要注意——DSC 的升温速率会影响结果，我建议至少做 3 个不同速率（比如 5、10、20 °C/min）的测试，才能得到可靠的参数。

说到固化动力学模型，Kamal 模型绝对是应用最广的。它的形式是这样的：

dα/dt = (k1 + k2·α^m)·(1 - α)^n

其中：

这个模型的好处是能描述自催化反应——就是反应产物本身会加速反应进行的那种。EMC 的固化正好是这种类型，所以 Kamal 模型特别适合。

个人经验：我在标定 Kamal 模型参数时，发现 m 和 n 对仿真结果非常敏感。有一次我把 m 从 0.5 改成 0.6，翘曲预测值直接变了 15%。所以建议大家做参数敏感性分析，别盲目相信默认值。

下面是我常用的参数标定流程：

我曾经遇到过一个案例，客户给的 EMC 参数拟合度 R² 只有 0.85，我坚持让他们重新标定。后来发现是 DSC 样品量太大，导致热滞后严重。减少样品量后，R² 提升到了 0.98。

固化反应是放热反应，这个热量会反过来影响温度场。在厚器件或者大尺寸封装中，固化放热可能导致内部温度比设定温度高 20-30°C。你想想看，这会对固化度分布产生多大影响？

在仿真中，我们需要把固化放热作为热源项加入能量方程：

ρ·Cp·∂T/∂t = ∇·(k·∇T) + ρ·Hr·dα/dt

其中 Hr 是总反应热，dα/dt 是固化速率。这个方程把温度场和固化度场耦合起来了——温度影响固化速率，固化放热又影响温度。说白了，这是个双向耦合问题。

注意：很多初学者在做模流分析时，只开了流动模块，忘了开固化放热。结果预测的固化时间比实际短了 30% 以上。我在 Moldflow 和 Moldex3D 里都踩过这个坑，现在每次建模都会检查「Heat of Reaction」选项有没有勾上。

固化放热对温度场的影响，在以下场景尤其明显：

我记得有一次做 BGA 封装的模流分析，客户反映中心区域的芯片总是出现裂纹。我一看仿真结果，中心温度比边缘高了 18°C，固化度快了 40%。这就是典型的固化放热导致的局部应力集中。后来我们调整了模温曲线，降低了升温速率，问题就解决了。

EMC 在固化过程中会发生体积收缩，这个收缩率通常在 0.2% 到 0.8% 之间。别小看这不到 1% 的收缩，在封装结构中，它产生的应力足以导致翘曲、分层甚至芯片开裂。

固化收缩的机理主要有两个：

在仿真中，我们通常把固化收缩等效为热膨胀来处理。具体做法是：

ε_shrink = C_shrink · Δα

其中 C_shrink 是固化收缩系数，Δα 是固化度变化量。这个等效方法虽然简单，但工程上够用。

避坑指南：我曾经在仿真中直接用 EMC 供应商给的「总收缩率」数据，结果翘曲预测完全不对。后来发现，总收缩率是化学收缩和热收缩的叠加，而仿真中热收缩已经通过 CTE 考虑了。正确的做法是只输入化学收缩部分，否则会重复计算。

固化收缩对翘曲的影响，可以通过以下方式缓解：

我个人在做应力仿真时，会把固化收缩和热膨胀分开处理。先算固化过程中的化学收缩应力，再算冷却过程中的热应力，最后叠加。这样虽然麻烦一点，但能更清楚地看到每个阶段的贡献。

下面这张图总结了本章的核心内容，我建议你保存下来，做仿真时对照着看：

好了，关于固化反应动力学的内容就讲到这里。这四个知识点——固化度、Kamal 模型、固化放热、固化收缩——是环环相扣的。你只要把它们的耦合关系理清楚，模流和应力仿真就能做到八九不离十。

公众号:蓝海资料掘金营，微信deep3321