第三章粘弹性流体的本构方程

第三章非线性粘弹流体的本构方程

1．本构方程概念

本构方程（constitutive equation），又称状态方程——描述一大类材料所遵循的与材料结构属性相关的力学响应规律的方程。

不同材料以不同本构方程表现其最基本的物性，对高分子材料流变学来讲，寻求能够正确描述高分子液体非线性粘弹响应规律的本构方程无疑为其最重要的中心任务，这也是建立高分子材料流变学理论的基础。

两种。

唯象性方法，一般不追求材料的微观结构，而是强调实验事实，现象性地推广流体力学、弹性力学、高分子物理学中关于线性粘弹性本构方程的研究结果，直接给出描写非线性粘弹流体应力、应变、应变率间的关系。以本构方程中的参数，如粘度、模量、松弛时间等，表征材料的特性。

分子论方法，重在建立能够描述高分子材料大分子链流动的正确模型，研究微观结构对材料流动性的影响。采用热力学和统计力学方法，将宏观流变性质与分子结构参数（如分子量，分子量分布，链段结构参数等）联系起来。为此首先提出能够描述大分子链运动的正确模型是问题关键。

根据研究对象不同，

象性方法和分子论方法虽然出发点不同，逻辑推理的思路不尽相同，而最终的结论却十分接近，表明这是一个正确的科学的研究基础。

目前关于高分子材料，特别浓厚体系本构方程的研究仍十分活跃。

同时，大量的实验积累着越来越多的数据，它们是检验本构方程优劣的最重要标志。

从形式上分，

速率型本构方程，方程中包含应力张量或形变速率张量的时间微商，或同时包含这两个微商。

积分型本构方程，利用迭加原理，把应力表示成应变历史上的积分，或者用一系列松弛时间连续分布的模型的迭加来描述材料的非线性粘弹性。积分又分为单重积分或多重积分。

判断一个本构方程的优劣主要考察：

1）方程的立论是否科学合理，论据是否充分，结论是否简单明了。

2）一个好的理论，不仅能正确描写已知的实验事实，还应能预言至今未知，但可能发生的事实。

3）有承前启后的功能。例如我们提出一个描写非线性粘弹流体的本构方程，当条件简化时，它应能还原为描写线性粘弹流体的本构关系。

4）最后也是最重要的一条，即实验事实（实验数据）是判断一个本构方程优劣的出发点和归宿。实践是检验真理的唯一标准。

本章重点介绍用唯象论方法对一般非线性粘弹流体建立的本构方程。分子论方法在第四章介绍。

2． 速率型本构方程

2．1 经典的线性粘弹性模型——Maxwell 模型

已知高分子本体的线性粘弹行为可以用一些力学模型，如Maxwell 模型、Voigt 模型及它们的恰当组合进行描述。其中Maxwell 模型由一个虎克型弹簧和一个牛顿型粘壶串联而成（图3-1）。由于形变时粘壶不受弹簧约束，可产生大形变。原则上Maxwell 模型可用于描述液体流动的性质。

图3-1 Maxwell 模型

设液体在剪切力作用下发生流动，弹簧、粘壶同时发生形变。

对弹簧有 11γσG =

对粘壶有 202γησ&=

因为串联，总应力 21σσσ==

总应变 σησγγγ0

2111+=+=&&&&G 所以有 γησλσ&&01=+ （3-1）

式中 G /01ηλ= 称松弛时间 ，单位为秒； （3-2） t

∂∂=σσ&

（3-3） 将（3-1）式推广写成三维形式，以张量表示，则有

d 012ηλ=+σσ& （3-4）

式中：σ 为偏应力张量； d 为形变率张量

()

2/T L L d += （3-5）

L 为速度梯度张量。注意这儿的推广是将方程简单地从一维形式推广到三维形式，并无深刻物理意义。公式中系数2的出现是由于采用了张量描述的缘故。

例1 Maxwell 模型用于描述稳态简单剪切流场。

简单剪切流场形式见图2-3，其中速度场方程见公式（2-46）。

我们在固定坐标系中考察流场中某一确定点上材料流过时的应力状态。由于流场是稳定的，因此该点的应力状态不随时间变化，故有 0=σ

& 对于稳态简单剪切流场，其形变率张量为

⎪⎪⎪⎭

⎫ ⎝⎛=000002/02/0γγ&&d （3-6）

代入（3-4）式，得到

⎪⎪⎪⎭⎫ ⎝⎛333231232221131211σσσσσσσσσ=20η⎪⎪⎪⎭

⎫ ⎝⎛000002/02/0γγ&&

这是一个由九个方程组成的方程组。由此解得：

⎪⎪⎭

⎪⎪⎬⎫

=-=-======000332222113113322302112σσσσσσσσγησσ& （3-7）

结果表明，采用Maxwell模型确实能描述材料在稳态简单剪切流场中的流动，但是模型的描述能力很有限。实际上它只能描述具有常数粘度η0的牛顿型流体的粘性行为，高分子液体在剪切速率极低情况下（γ&→0）的流动状态（具有常数粘度）也可用该模型近似描述。

对于非牛顿型流体在一般流场中的非线性粘弹行为，Maxwell模型无能为力。既不能描述高分子液体典型的剪切变稀（即结构粘性）行为，也不能描述流动中存在法向应力差（即具有弹性）的事实。（3-7）式中给出的两个法向应力差值均等于零。

分析可知，Maxwell模型有限的描述能力与方程的推广方式有关，特别与方程中应力张量的导数形式有关。（3-1）式中描述的应力变化的导数形式σ&是应力对时间的一般偏微商，这种偏微商通常只能描述无穷小形变行为，或流动中体系性质无变化的形变行为。对于描述高分子液体在大形变下的非线性粘弹行为，必须对应力张量的导数形式审慎定义和推广。

另外，在考察流场中流体流动时，紧盯着固定坐标系的一点考察（注意在不同时刻流经该点的流体元不同）和紧跟着一个流体元考察（该流体元在不同时刻占据空间不同位置）是大不相同的。为此我们首先介绍流体力学中描写材料元流动的空间描述法和物质描述法，然后再讨论经典Maxwell模型的推广。

2．2空间描述法和物质描述法

流体力学中，在固定的空间坐标系描写一个材料元的流动有两种不同方法：

一是物质描述法，观察者的视点集中于一个具体的流体元及其邻域所发生的事件，研究它在不同时刻所处的位置，以及它的速度，加速度等，