3.1 广义牛顿流体模型

本节摘要：广义牛顿流体模型把黏度写成剪切速率的标量函数，用一条流动曲线直接换算黏度规律，是工程上用得最多的快速破案工具。幂律、Carreau、Bingham 与 Herschel-Bulkley 各有适用边界，参数辨识则是决定拟合是否可信的关键一环。

学习目标

阅读完本节，你应当能够：

说明广义牛顿模型与牛顿模型之间的推广逻辑；
为给定曲线挑选合适的广义牛顿模型，并说出各参数的物理含义；
判断哪些场景该用带屈服应力的模型，哪些不该；
检查自己的拟合是否存在参数不可辨识的问题。

一个函数的诞生

牛顿模型的全部家当就是一个比例常数：剪应力 τ 等于黏度 η 乘以剪切速率 γ̇。可实测曲线一出来，这个"常数"立刻露馅——番茄酱静置时几乎不流，一拍瓶底倾泻而出，黏度显然在变。把常数升级成函数，问题就解决了：让黏度成为局部剪切速率的函数 η(γ̇)，本构关系随之推广为：

\boldsymbol{\tau} = \eta(\dot{\gamma})\,\mathbf{D}

这一步看着只是把 η 后面的下标去掉，实际上承认了两件事：一是材料响应是非线性的，黏度随流动强弱而变；二是这个模型仍然"没记性"——应力只由当前应变速率决定，不记得过去，也没有法向应力差。代价要记清：广义牛顿模型只适合粘性主导、弹性可忽略的过程。挤出口模的胀大、搅拌棒的爬杆，它一律答不上来。

幂律：简单，但别外推

最朴素的非线性扩展是幂律模型：

\eta(\dot{\gamma}) = K \dot{\gamma}^{\,n-1}

指数 n	行为	典型材料
n 等于 1	退化为牛顿流体	水、溶剂
n 小于 1	剪切变稀	聚合物溶液、熔体
n 大于 1	剪切增稠	浓颗粒悬浮液

优点就是简单：两个参数，解析解好算，管道流动的压降能直接积分出来。致命伤也在简单——γ̇ 趋于零时黏度发散到无穷，趋于无穷时黏度压到零，两头都不像话。所以幂律只配在中段用，把它外推到零切区或无穷剪切区，是新手最容易犯的错。我的建议：幂律用来快速估工艺窗口可以，但报告里务必标注适用范围。

双平台家族：Carreau 与 Cross

要描述真实的完整曲线，得让黏度在两头都有平台。Carreau 模型提供了从零切黏度到无穷切黏度的平滑过渡：

\eta = \eta_\infty + (\eta_0 - \eta_\infty)\left[1 + (\lambda\dot{\gamma})^2\right]^{(n-1)/2}

Cross 模型则是另一套写法，分母从幂变为 1 加某次方：

\eta = \eta_\infty + \frac{\eta_0 - \eta_\infty}{1 + (\lambda\dot{\gamma})^m}

参数	含义	对应曲线特征
η0	零切黏度	低速平台高度
η∞	无穷切黏度	高速平台高度
λ	特征时间	拐点位置，约等于 1 除以拐点速率
n 或 m	变稀指数	中段下坡的陡峭度

两个模型都能画双平台曲线，差别在过渡区的写法。Carreau 的过渡更圆滑，数值求解时高阶次幂可能带来稳定性麻烦；Cross 结构更省事，商业软件里用得更多。选择标准很实际：哪个拟合残差小、哪个在后续求解里稳，就用哪个，不必迷信哪家。

有门槛的：Bingham 与 Herschel-Bulkley

上面几位的共同设定是"只要有一点应力就开始流动"。可发胶、水泥浆这类材料有门槛——应力低于屈服值时纹丝不动，超过才流。于是有了带屈服应力的家族：

\tau > \tau_y \text{ 时，}\ \tau = \tau_y + K\dot{\gamma}^n

n 取 1 就是 Bingham 模型，n 小于 1 就是 Herschel-Bulkley。这里有个争议值得点破：屈服应力到底是不是材料的固有属性，学术界吵了很多年——有人论证它只是"极低剪切速率下还没测到流动"的测量假象。工程上别纠结这个哲学问题，看场合：涂布停挂、灌装保形这类"静置不能流"的设计，带屈服项才靠谱；但如果只是描述稳态剪切曲线，先试不带屈服项的模型，数据不够再加。

拟合不是玄学：参数辨识

模型挑好，参数怎么定？四条初值经验：零切黏度用最低速率点的读数，无穷切黏度用最高速率点的读数，特征时间取曲线拐点附近，变稀指数量中段斜率。然后做非线性最小二乘，收敛后必须检查每个参数的标准差——这一步最容易被跳过，也最容易出事。


# 演示：数据只覆盖中段时，零切黏度不可辨识
# 两组悬殊的参数都能把中段拟合得很像，外推却天差地别
def eta_carreau(rate, eta0, lam, n):
    return eta0 * (1 + (lam * rate) ** 2) ** ((n - 1) / 2)

rates = [0.1, 0.5, 1.0, 2.0, 5.0, 10.0]   # 只有拐点附近的点
data =  [8.2, 6.4, 5.5, 4.8, 4.0, 3.6]    # 假想的实测黏度

# 组一：eta0 取 10，lam 取 0.6，n 取 0.4
# 组二：eta0 取 30，lam 取 0.15，n 取 0.4
# 两组在中段的拟合线几乎重合，外推到低速段却差出好几倍

判读原则只有一句：若某个参数的标准差与自身同量级，说明数据里没有支撑它的信息——补测对应段落的实验，而不是硬收这个拟合值。零切黏度标准差大，就去补低速率点；无穷切黏度标准差大，就去补高速率点。这个检查做完，拟合才算数。

要点小结

推广逻辑：牛顿常数变函数，代价是丢掉记忆与法向应力差。
按范围选型：只有中段用幂律，有双平台用 Carreau 或 Cross，有门槛加屈服项。
参数对特征：η0 对低速平台、λ 对拐点、指数对坡度，一一对应。
辨识优先：参数标准差大就先补数据，别硬拟合。

广义牛顿模型审不了"记性"。下一节让弹簧与粘壶登场，把记忆写进方程。

作者与出处

原作者：灏天文库

来源：灏天文库

整理：灏天文库整理

由灏天文库平台收录，内容或由平台用户上传，仅供学习交流