第3章：模流分析理论基础

各位工程师朋友，大家好。今天我们来聊聊模流分析的理论基础。说实话，很多做封装的人一听到「流体力学」就头大，觉得那是搞CFD的人该操心的事。但我要说，不懂这些基本方程，你连仿真结果对不对都判断不了。

我自己刚入行那会儿，就吃过这个亏。有一次做QFN的模流分析，结果死活跟实际填充对不上。折腾了两周，最后发现是粘度模型参数设错了。嗯，从那以后，我老老实实把基础补了一遍。

3.1 流体力学基本方程

模流分析，说白了就是求解三个方程：连续性方程、动量方程、能量方程。这三个方程一起，构成了流体运动的基本描述。你想想看，任何流动问题，本质上都是在解这三组方程。

3.4.1 连续性方程（质量守恒）

这个方程讲的是：流入控制体的质量，等于流出控制体的质量加上内部质量的变化。用数学表达就是：

∂ρ/∂t + ∇·(ρv) = 0

其中ρ是密度，v是速度矢量。对于环氧塑封料，我们通常假设为不可压缩流体，所以密度ρ是常数，方程简化为：

∇·v = 0

我在项目中遇到过一个问题：有人把EMC当成完全不可压缩来处理，结果仿真出来的填充时间总是偏短。后来发现，EMC在高压下其实有微弱的可压缩性，特别是对于大尺寸封装，这个效应不能忽略。

3.4.2 动量方程（牛顿第二定律）

动量方程，其实就是F=ma在流体中的表现形式。对于粘性流体，它的完整形式是Navier-Stokes方程：

ρ(∂v/∂t + v·∇v) = -∇p + ∇·τ + ρg

左边是惯性项，右边依次是压力项、粘性项、重力项。对于EMC模流，我们通常忽略重力，因为熔体粘度很大，重力影响微乎其微。

这里有个关键点：EMC的流动是典型的低雷诺数流动（Re < 1），所以惯性项也可以忽略。这就变成了Stokes流动，求解起来简单很多。

3.4.3 能量方程（能量守恒）

能量方程描述的是温度场的变化。对于EMC模流，温度直接影响粘度，粘度又影响流动，所以能量方程和动量方程是耦合的。

ρCp(∂T/∂t + v·∇T) = ∇·(k∇T) + τ:∇v + Q

左边是温度变化和对流项，右边依次是热传导、粘性耗散、反应热。其中粘性耗散项在高速流动中很重要，但EMC模流速度慢，通常可以忽略。

反应热Q来自环氧树脂的固化反应。这个项很关键，因为固化反应是放热的，会导致局部温度升高，反过来加速固化。我曾经见过一个案例，就是因为忽略了反应热，仿真出来的固化度比实际低了20%。

3.2 非牛顿流体模型

EMC不是水，也不是油。它是典型的非牛顿流体，而且是非牛顿流体中的剪切变稀型。什么意思呢？就是流速越快，它越稀。你想想看，这跟番茄酱有点像——使劲挤的时候，它反而更容易流出来。

3.2.1 为什么EMC是非牛顿流体？

EMC由环氧树脂、固化剂、填料、添加剂组成。其中填料颗粒（二氧化硅）占60-90%的重量。这些颗粒在流动中会相互作用，形成复杂的微观结构。当剪切速率低时，颗粒之间形成网络，粘度很高；当剪切速率高时，网络被破坏，粘度下降。

这就是剪切变稀的物理本质。我刚开始做仿真时，直接用牛顿流体模型算，结果填充模式完全不对。后来换成非牛顿模型，才跟实际对得上。

3.2.2 Cross-WLF粘度模型

描述非牛顿流体的粘度模型有很多，比如Power Law、Carreau、Cross等。但在EMC模流分析中，最常用的是Cross-WLF模型。为什么？因为它能同时描述剪切变稀和温度依赖两个效应。

Cross-WLF模型的数学形式：

η(T, γ̇) = η₀(T) / [1 + (η₀γ̇/τ*) ^ (1-n)]

其中η₀是零剪切粘度，由WLF方程描述：

η₀(T) = D₁ exp[-A₁(T - T*)/ (A₂ + T - T*)]

参数说明：

参数	含义	典型值范围
n	幂律指数（n<1为剪切变稀）	0.3 - 0.6
τ*	临界剪切应力	10⁴ - 10⁵ Pa
D₁	参考粘度	10⁶ - 10¹⁰ Pa·s
A₁	WLF常数	20 - 40
A₂	WLF常数	50 - 100 K
T*	参考温度（通常取Tg）	100 - 150 °C

3.2.3 固化度对粘度的影响

EMC在流动过程中，固化反应一直在进行。随着固化度α增加，粘度会急剧上升。这个效应在Cross-WLF模型中通常通过一个固化因子来考虑：

η(α) = η₀ × (α_gel / (α_gel - α)) ^ (C₁ + C₂α)

其中α_gel是凝胶点固化度，C₁、C₂是材料常数。当α接近α_gel时，粘度趋于无穷大，流动停止。

我曾经犯过一个错误：忽略了固化度对粘度的影响，结果仿真显示填充很顺利，但实际生产时出现了短射。后来发现，是因为流动路径太长，固化度已经接近凝胶点，粘度太大导致填充不完全。

3.3 知识体系总览

为了让大家更直观地理解本章的知识结构，我画了一张图：

这张图把本章的核心内容串起来了。左边是三大基本方程，中间是非牛顿特性，右边是具体的Cross-WLF模型参数。你想想看，这三个部分缺一不可。

好了，这一章的内容就到这里。记住，理论是基础，但最终要落到实际参数标定上。下一章我们会讲如何用流变仪测量这些参数，以及如何把数据导入仿真软件。