模型公式
幂律模型:
$$\eta = m \dot{\gamma}^{n-1}$$
Carreau模型:
$$\eta = \eta_\infty + (\eta_0 - \eta_\infty)[1+(\lambda\dot{\gamma})^2]^{\frac{n-1}{2}}$$
Cross模型:
$$\eta = \eta_\infty + \frac{\eta_0 - \eta_\infty}{1+(m\dot{\gamma})^n}$$
在同一图表上实时比较幂律、Carreau和Cross三种粘度模型。调整流变参数,即时查看剪切速率0.001~1000 s⁻¹范围内的对数-对数粘度曲线。
幂律模型:
$$\eta = m \dot{\gamma}^{n-1}$$
Carreau模型:
$$\eta = \eta_\infty + (\eta_0 - \eta_\infty)[1+(\lambda\dot{\gamma})^2]^{\frac{n-1}{2}}$$
Cross模型:
$$\eta = \eta_\infty + \frac{\eta_0 - \eta_\infty}{1+(m\dot{\gamma})^n}$$
幂律模型 (Power Law):最简单的经验模型,描述粘度与剪切速率之间的幂次关系,适用于中等剪切速率范围。
$$\eta = m \dot{\gamma}^{n-1}$$$\eta$:表观粘度 (Pa·s);$m$:稠度系数 (Pa·s),代表粘度大小;$n$:幂律指数,$n<1$为剪切稀化;$\dot{\gamma}$:剪切速率 (s⁻¹)。
Carreau模型:基于分子网络理论的物理模型,能完整描述零剪切牛顿区、剪切稀化区和第二牛顿区。
$$\eta = \eta_\infty + (\eta_0 - \eta_\infty)[1+(\lambda\dot{\gamma})^2]^{\frac{n-1}{2}}$$$\eta_0$:零剪切粘度;$\eta_\infty$:无限剪切粘度;$\lambda$:松弛时间 (s),决定过渡的快慢;$n$:幂律指数,控制剪切稀化区的斜率。
Cross模型:另一种广泛使用的经验模型,形式略有不同,同样能描述三个区域,在涂料、粘合剂领域应用广泛。
$$\eta = \eta_\infty + \frac{\eta_0 - \eta_\infty}{1 + (m \dot{\gamma})^n}$$$m$:Cross常数 (s),具有时间量纲,与Carreau模型的$\lambda$物理意义类似;其他变量含义与Carreau模型相同。
注塑成型CAE仿真:在Moldflow、C-Mold等软件中,Carreau-WLF模型是模拟塑料熔体充填过程的核心。准确的粘度模型能预测熔体前沿的流动形态、注射压力、以及可能出现的短射或飞边等缺陷。
挤出机与流道设计:设计螺杆和模具流道时,工程师需要知道熔体在不同剪切速率(对应螺杆不同位置)下的粘度变化,以优化混合效果、降低能耗并确保出料均匀。幂律模型常用于初步估算。
涂料与粘合剂配方:对于油漆、胶水等材料,其“刷涂性”和“抗流挂性”分别对应中剪切和低剪切下的粘度。Cross模型常被用来拟合实验数据,以调整配方达到理想的施工性能。
纤维增强复合材料成型:在树脂传递模塑(RTM)中,树脂浸润纤维毡的过程涉及极低的剪切速率。此时零剪切粘度$\eta_0$是关键参数,它直接影响充填时间和最终制品质量。
开始使用此工具时,CAE初学者常会陷入几个误区。首先是“预设材料参数并非绝对常数”。例如即便是“PC”材料,其粘度特性也会因生产商和牌号(高流动型、阻燃型等)产生显著差异。工具的预设值仅为代表性数据。关键是要养成将实际成型材料的 datasheet 中的粘度曲线与工具结果进行比对,必要时对λ和n进行微调的习惯。
第二点是忽视浇口压力损失计算的前提条件。计算公式基于“完全发展的圆形流道内层流”这一假设。这意味着实际浇口入口处的收缩损失、以及模温导致的壁面固化层(冻结皮层)的影响均未包含在内。计算值应视为理论上的最小压力损失,实际应用中通常需预留1.5至2倍的安全余量,这是现场经验的智慧。
第三点是切勿过度依赖冷却时间计算。工具的冷却时间仅是基于恒定模壁温度热传导的理论计算。但实际成型中涉及冷却水路布局、流量、树脂结晶放热等复杂因素。该计算值应作为“把握冷却时间量级”的参考。例如若计算值为10秒,实际周期时间通常在12~15秒左右,建议以这种关联性来理解。
本工具涉及的非牛顿流体粘度模型与流动分析,是除注塑成型外众多工程领域的基础。例如涂料与化妆品的流变学:从管中缓慢挤出的发胶与取到手上快速抹开时触感的差异,正是剪切速率依赖性的体现。在质量控制中,会使用与本工具Carreau模型类似的模型来评估产品特性。
另一领域是食品工程。番茄酱和蛋黄酱缓慢倾倒时不易流出(高零剪切粘度),但摇晃或施加较大外力时就会流动(剪切稀化)。这虽伴随触变性等时间依赖性,但对基本剪切稀化行为的理解是相通的。此外在血液流变学(生物流变学)中,血管内剪切速率导致血液粘度变化的现象,在数学上与聚合物熔体分析具有相似性。
由此可见,通过一个工具学到的“非牛顿流体思维方式”,是即使材料改变仍可通用的普适性工程语言。
熟悉本工具计算后若产生“为什么?”的疑问,建议进入下一阶段学习。首先推荐概览“流变学”这门学科。特别要注意,粘度曲线在低剪切速率区出现的第一牛顿区,与分子链缠结状态下的松弛现象密切相关。要理解参数“λ(松弛时间)”的物理意义,学习麦克斯韦模型等基础粘弹性模型会令人豁然开朗。
数学背景方面,建议探究非牛顿流体力学中幂律公式的推导过程。从假设圆管内流速分布出发,由剪切应力平衡导出的下式是起点: $$ \Delta P \cdot \pi r^2 = \tau \cdot 2\pi r L $$ 由此,通过幂律关系($\tau = m \dot{\gamma}^n$)连接剪切应力τ与剪切速率$\dot{\gamma}$,即可推导出工具中使用的压力损失公式。理解此推导过程后,公式中“3n+1”等系数的含义也会变得清晰。
接下来的实践性课题,可尝试“拓展至模具内多样流道形状(矩形、薄平板)的解析”。存在利用水力半径概念近似计算圆形流道公式的方法。掌握此法后,就能更接近实际模具型腔或流道内更真实的压力损失预测。