非线性粘弹流体本构方程.docVIP

非线性粘弹流体的本构方程 本构方程概念 本构方程（constitutive equation），又称状态方程——描述一大类材料所遵循的与材料结构属性相关的力学响应规律的方程。 不同材料以不同本构方程表现其最基本的物性，对高分子材料流变学来讲，寻求能够正确描述高分子液体非线性粘弹响应规律的本构方程无疑为其最重要的中心任务，这也是建立高分子材料流变学理论的基础。 寻求流变本构方程的基本方法大致可分为唯象性方法和分子论方法两种。 唯象性方法，一般不追求材料的微观结构，而是强调实验事实，现象性地推广流体力学、弹性力学、高分子物理学中关于线性粘弹性本构方程的研究结果，直接给出描写非线性粘弹流体应力、应变、应变率间的关系。以本构方程中的参数，如粘度、模量、松弛时间等，表征材料的特性。 分子论方法，重在建立能够描述高分子材料大分子链流动的正确模型，研究微观结构对材料流动性的影响。采用热力学和统计力学方法，将宏观流变性质与分子结构参数（如分子量，分子量分布，链段结构参数等）联系起来。为此首先提出能够描述大分子链运动的正确模型是问题关键。 根据研究对象不同，高分子流变本构方程又分为稀溶液理论及浓厚体系理论，两部分的理论和实验研究工作都取得巨大的进展。有趣的是唯象性方法和分子论方法虽然出发点不同，逻辑推理的思路不尽相同，而最终的结论却十分接近，表明这是一个正确的科学的研究基础。 目前关于高分子材料，特别浓厚体系本构方程的研究仍十分活跃。 同时，大量的实验积累着越来越多的数据，它们是检验本构方程优劣的最重要标志。 从形式上分，非线性粘弹流体的本构方程主要分为两大类：速率型（亦称微商型）本构方程和积分型本构方程。 速率型本构方程，方程中包含应力张量或形变速率张量的时间微商，或同时包含这两个微商。 积分型本构方程，利用迭加原理，把应力表示成应变历史上的积分，或者用一系列松弛时间连续分布的模型的迭加来描述材料的非线性粘弹性。积分又分为单重积分或多重积分。 判断一个本构方程的优劣主要考察： 1）方程的立论是否科学合理，论据是否充分，结论是否简单明了。 2）一个好的理论，不仅能正确描写已知的实验事实，还应能预言至今未知，但可能发生的事实。 3）有承前启后的功能。例如我们提出一个描写非线性粘弹流体的本构方程，当条件简化时，它应能还原为描写线性粘弹流体的本构关系。 4）最后也是最重要的一条，即实验事实（实验数据）是判断一个本构方程优劣的出发点和归宿。实践是检验真理的唯一标准。 本章重点介绍用唯象论方法对一般非线性粘弹流体建立的本构方程。分子论方法在第四章介绍。 速率型本构方程 经典的线性粘弹性模型——Maxwell模型 已知高分子本体的线性粘弹行为可以用一些力学模型，如Maxwell模型、Voigt模型及它们的恰当组合进行描述。其中Maxwell模型由一个虎克型弹簧和一个牛顿型粘壶串联而成（图3-1）。由于形变时粘壶不受弹簧约束，可产生大形变。原则上Maxwell模型可用于描述液体流动的性质。 图3-1 Maxwell模型 设液体在剪切力作用下发生流动，弹簧、粘壶同时发生形变。 对弹簧有 对粘壶有 因为串联，总应力 总应变 所以有 （3-1） 式中 称松弛时间 ，单位为秒； （3-2） 为应力σ对时间的一般偏微商。 （3-3） 将（3-1）式推广写成三维形式，以张量表示，则有 （3-4） 式中： 为偏应力张量； d 为形变率张量 （3-5） L为速度梯度张量。注意这儿的推广是将方程简单地从一维形式推广到三维形式，并无深刻物理意义。公式中系数2的出现是由于采用了张量描述的缘故。 Maxwell模型用于描述稳态简单剪切流场。 简单剪切流场形式见图2-3，其中速度场方程见公式（2-46）。 我们在固定坐标系中考察流场中某一确定点上材料流过时的应力状态。由于流场是稳定的，因此该点的应力状态不随时间变化，故有 对于稳态简单剪切流场，其形变率张量为 （3-6） 代入（3-4）式，得到 =2 这是一个由九个方程组成的方程组。由此解得： （3-7） 结果表明，采用Maxwell模型确实能描述材料在稳态简单剪切流场中的流动，但是模型的描述能力很有限。实际上它只能描述具有常数粘度η0的牛顿型流体的