2. 聚合物材料流变与热性能基础
注塑成形的流动与翘曲仿真，本质上是在求解“熔体在型腔中如何流动、热量如何传递、体积如何随温压变化”这三类耦合问题。而决定这三类问题求解精度的，不是网格有多密，而是材料模型是否真实。材料数据是仿真的物理边界条件，材料参数错了，后续所有计算都只是在精确地算错。本节围绕粘度模型、剪切速率与温度依赖性、PVT数据、比热与导热系数，以及材料数据库的使用与测试方法展开。
一、粘度模型：把“流动阻力”写成可计算的函数
聚合物熔体是典型的非牛顿流体，其粘度不是常数，而是剪切速率、温度和压力的函数。仿真中最常用的两个模型是幂律模型和Cross-WLF 模型。
幂律模型形式简单：η = m·γ̇^(n−1)，其中 m 为稠度系数，n 为幂律指数。对大多数聚合物熔体，n 在 0.2～0.7 之间，n 越小剪切变稀越明显。幂律模型在中等剪切速率区间拟合良好，但在极低剪切速率下会高估粘度、在极高剪切速率下缺乏上平台，因此更适合做快速估算。
Cross-WLF 模型则同时描述低剪切牛顿平台、剪切变稀区和温度依赖：
η(T, γ̇, P) = η0(T,P) / [1 + (η0·γ̇ / τ*)^(1−n)]
η0(T,P) = D1 · exp[ −A1·(T − T*) / (A2 + (T − T*)) ]
T* = D2 + D3·P
A2 = Ã2 + D3·P其中 η0 是零剪切粘度，τ* 是剪切应力转折点，n 为幂律指数，D1、D2、D3、A1、Ã2 为材料常数。Cross-WLF 能覆盖注塑中从浇口附近的高剪切到型腔末端的低剪切全范围，是 Moldflow、Moldex3D 等软件的默认首选。选择模型时，若仿真关注填充末端滞流、保压补缩，应优先用 Cross-WLF；若只做初步流道平衡估算，幂律模型足够。
二、剪切速率与温度依赖性
注塑过程中的剪切速率跨度极大：浇口处可达 10⁴～10⁵ s⁻¹，型腔末端可能低至 10⁰～10¹ s⁻¹。粘度对剪切速率的敏感性决定了“剪切变稀”带来的流动性改善，也决定了剪切生热与分子取向。温度方面，聚合物粘度随温度升高而下降，可用 Arrhenius 或 WLF 方程描述。WLF 方程在接近玻璃化转变温度时更准确，Arrhenius 在远离 Tg 的高温区更简便。
实际仿真中要特别注意：温度依赖性参数必须与所用模型匹配。Cross-WLF 自带 WLF 形式的温度项，若再叠加一个 Arrhenius 修正，会造成重复计算。此外，压力对粘度的影响在高压保压阶段不可忽略，D3 项正是为此设置。
三、PVT 数据：体积收缩与翘曲的根源
PVT（压力-比容-温度）关系直接决定制品的体积收缩率，是翘曲预测的核心输入。常用 Tait 方程描述：
v(T,P) = v0(T)·[1 − C·ln(1 + P/B(T))] + vT(T,P)
v0(T) = b1 + b2·(T − b5) （T > Tg）
v0(T) = b1 + b2·(T − b5) + b7·... （T < Tg）其中 b1～b9 为材料常数，Tg 为玻璃化转变温度。PVT 曲线在 Tg 处出现斜率突变，对应比容随温度变化的“拐点”。仿真中若忽略这一拐点，脱模后的收缩预测会明显偏差。判断 PVT 数据是否可用，可检查三点：曲线是否覆盖实际注塑温度与压力范围；Tg 拐点是否清晰；冷却段数据是否完整。
四、比热与导热系数
比热 Cp 决定熔体降温速率和凝固层生长速度，导热系数 k 决定热量从熔体向模具钢传递的效率。两者都随温度变化：Cp 在 Tg 附近出现阶跃，k 通常随温度升高略有下降。仿真中若把 Cp、k 当作常数，冷却时间预测可能偏差 15%～30%。对于结晶型聚合物（如 PP、PA），还需注意结晶潜热的释放，通常以等效比热峰的形式并入 Cp(T) 曲线。
五、材料数据库使用与测试方法
主流仿真软件自带材料数据库，但数据库中的牌号数据来源不一，必须核对测试标准与测试条件。使用材料数据库时建议按以下顺序核查：
- 牌号匹配：确认供应商、牌号、颜色、增强比例完全一致，玻纤含量差 5% 就可能显著改变粘度与收缩。
- 测试方法：粘度用毛细管流变仪或旋转流变仪测定，PVT 用高压膨胀计（PVT 测试仪），Cp 用 DSC，k 用激光闪射法。不同方法得到的数据适用范围不同。
- 数据范围：检查测试温度、压力、剪切速率是否覆盖实际工艺窗口。若数据库只测到 240°C，而实际注塑 280°C，则需外推并标注不确定性。
- 拟合质量：查看软件拟合曲线与原始数据点的重合度，Cross-WLF 拟合在低剪切区偏差大时，应手动调整 D1 或 τ*。
六、参数敏感性：先抓主要矛盾
材料参数众多，但并非同等重要。对填充阶段，粘度模型中的 n 和 τ* 影响最大；对保压与收缩，PVT 的 b1、b2、b4 和 Tg 最关键；对冷却时间，Cp 与 k 共同决定。建议在正式优化前做一次单因素敏感性扫描：将每个关键参数分别上下浮动 10%，观察填充压力、体积收缩率和翘曲量的变化幅度，优先修正敏感度最高的参数。这样可以用最少的测试成本，把仿真精度提升到可用于判断浇口位置、流道平衡和冷却方案的水平。