第三章粘弹性流体的本构方程.docx
- 文档编号:26227374
- 上传时间:2023-06-17
- 格式:DOCX
- 页数:21
- 大小:76.48KB
第三章粘弹性流体的本构方程.docx
《第三章粘弹性流体的本构方程.docx》由会员分享,可在线阅读,更多相关《第三章粘弹性流体的本构方程.docx(21页珍藏版)》请在冰豆网上搜索。
第三章粘弹性流体的本构方程
第三章非线性粘弹流体的本构方程
1.本构方程概念
本构方程(constitutiveequation),又称状态方程描述一大类
材料所遵循的与材料结构属性相关的力学响应规律的方程。
不同材料以不同本构方程表现其最基本的物性,对高分子材料流变学来讲,寻求能够正确描述高分子液体非线性粘弹响应规律的本构方程无疑为其最重要的中心任务,这也是建立高分子材料流变学理论的基础。
寻求流变本构方程的基本方法大致可分为唯象性方法和分子论方法两种。
唯象性方法,一般不追求材料的微观结构,而是强调实验事实,现象性地推广流体力学、弹性力学、高分子物理学中关于线性粘弹性本构方程的研究结果,直接给出描写非线性粘弹流体应力、应变、应变率间的关系。
以本构方程中的参数,如粘度、模量、松弛时间等,表征材料的特性。
分子论方法,重在建立能够描述高分子材料大分子链流动的正确模型,研究微观结构对材料流动性的影响。
采用热力学和统计力学方法,将宏观流变性质与分子结构参数(如分子量,分子量分布,链段结构参数等)联系起来。
为此首先提出能够描述大分子链运动的正确模型是问题关键。
根据研究对象不同,高分子流变本构方程又分为稀溶液理论及浓—
系理论,两部分的理论和实验研究工作都取得巨大的进展。
有趣的是唯象性方法和分子论方法虽然出发点不同,逻辑推理的思路不尽相同,而最终的结论却十分接近,表明这是一个正确的科学的研究基础。
目前关于高分子材料,特别浓厚体系本构方程的研究仍十分活跃。
同时,大量的实验积累着越来越多的数据,它们是检验本构方程优劣的最重要标志。
从形式上分,非线性粘弹流体的本构方程主要分为两大类:
I速率型1(亦称微商型)本构方程和积分型本构方程。
速率型本构方程,方程中包含应力张量或形变速率张量的时间微商,或同时包含这两个微商。
积分型本构方程,利用迭加原理,把应力表示成应变历史上的积分,或者用一系列松弛时间连续分布的模型的迭加来描述材料的非线性粘弹性。
积分又分为单重积分或多重积分。
判断一个本构方程的优劣主要考察:
1)方程的立论是否科学合理,论据是否充分,结论是否简单明了。
2)一个好的理论,不仅能正确描写已知的实验事实,还应能预言至今未知,但可能发生的事实。
3)有承前启后的功能。
例如我们提出一个描写非线性粘弹流体的本构方程,当条件简化时,它应能还原为描写线性粘弹流体的本构关系。
4)最后也是最重要的一条,即实验事实(实验数据)是判断一个本构方程优劣的出发点和归宿。
实践是检验真理的唯一标准。
本章重点介绍用唯象论方法对一般非线性粘弹流体建立的本构方程。
分子论方法在第四章介绍。
2.速率型本构方程
2.1经典的线性粘弹性模型——Maxwell模型
已知高分子本体的线性粘弹行为可以用一些力学模型,如Maxwell模
型、Voigt模型及它们的恰当组合进行描述。
其中Maxwell模型由一个虎克型弹簧和一个牛顿型粘壶串联而成(图3-1)。
由于形变时粘壶不受弹
簧约束,可产生大形变。
原则上Maxwell模型可用于描述液体流动的性
质。
图3-1Maxwell模型
设液体在剪切力作用下发生流动,
弹簧、
粘壶同时发生形变。
对弹簧有
1
G1
对粘壶有
2
02
因为串联,
总应力
1
2
总应变
1
2
丄
1
G
0
所以有
1
0
(3-1)
式中1
0/G
称松弛时间
单位为秒;(3-2)
——为应力。
对时间的一般偏微商。
(3-3)
t
将(3-1)式推广写成三维形式,以张量表示,则有
(X12od(3-4)
式中:
为偏应力张量;d为形变率张量
dLLT/2(3-5)
L为速度梯度张量。
注意这儿的推广是将方程简单地从一维形式推广到三
维形式,并无深刻物理意义。
公式中系数2的出现是由于采用了张量描述的缘故。
例1Maxwell模型用于描述稳态简单剪切流场。
简单剪切流场形式见图2-3,其中速度场方程见公式(2-46)。
我们在固定坐标系中考察流场中某一确定点上材料流过时的应力状
态。
由于流场是稳定的,因此该点的应力状态不随时间变化,故有
0
对于稳态简单剪切流场,其形变率张量为
0
/2
0
d/2
0
0
(3-6)
0
0
0
代入(3-4)式,得到
1112
13
0
/2
0
2122
23
=20
/2
0
0
3132
33
0
0
0
这是一个由九个方程组成的方程组。
由此解得:
12210
23
32
13
310
(3-7)
11
22
0
22
33
0
结果表明,采用Maxwell模型确实能描述材料在稳态简单剪切流场中的流动,但是模型的描述能力很有限。
实际上它只能描述具有常数粘度n0的牛顿型流体的粘性行为,高分子液体在剪切速率极低情况下(-0)
的流动状态(具有常数粘度)也可用该模型近似描述。
对于非牛顿型流体在一般流场中的非线性粘弹行为,Maxwell模型无
能为力。
既不能描述高分子液体典型的剪切变稀(即结构粘性)行为,也不能描述流动中存在法向应力差(即具有弹性)的事实。
(3-7)式中给
出的两个法向应力差值均等于零。
分析可知,Maxwell模型有限的描述能力与方程的推广方式有关,特别与方程中应力张量的导数形式有关。
(3-1)式中描述的应力变化的导数形式是应力对时间的一般偏微商,这种偏微商通常只能描述无穷小形变行为,或流动中体系性质无变化的形变行为。
对于描述高分子液体在大形变下的非线性粘弹行为,必须对应力张量的导数形式审慎定义和推广。
另外,在考察流场中流体流动时,紧盯着固定坐标系的一点考察(注意在不同时刻流经该点的流体元不同)和紧跟着一个流体元考察(该流体元在不同时刻占据空间不同位置)是大不相同的。
为此我们首先介绍流体力学中描写材料元流动的空间描述法和物质描述法,然后再讨论经典Maxwell模型的推广。
2.2空间描述法和物质描述法
流体力学中,在固定的空间坐标系描写一个材料元的流动有两种不同方法:
一是物质描述法,观察者的视点集中于一个具体的流体元及其邻域所发生的事件,研究它在不同时刻所处的位置,以及它的速度,加速度等,
与通常力学中集中于一个质点的方法相同。
这种方法又称拉格朗日描述
法。
在该方法中一般以流体元在参考构型中的物质坐标Xr(R=1,2,3)为
自变量,以便区别不同的材料元。
另一个方法称空间描述法,观察者的视点集中于坐标空间某一特殊点及其邻域所发生的事件,不针对一个具体的流体元。
这种方法又称欧拉描述法。
在该方法中,往往以固定坐标系的空间坐标Xi(i=1,2,3)为自变量。
一般偏导数和物质导数:
流场中的任一物理量u都是时间t和空间坐标Xi(i=1,2,3)的函数,记成ut,x1,x2,x3。
当求u的时间导数时,应当区分两种情况。
展开来写,有
(3-16)
是固定空间坐标Xi
数,称一般偏导数。
(i=1,2,3)不变(空间描述法),只对时间t求偏导
ut,X1,X2,X3
ut,X1,X2,X3
t
一是米用物质描述法,
紧盯着个材料兀求时间导数。
由于材料兀的
坐标也在变化(为时间t的函数),因此求导时不仅要对t求,也要对
也称u对时间求全导数。
(3-16)还可记成以下矢量形式:
(3-17)
3丄V?
u
Dtt
2.3广义Maxwell模型
考虑将经典的Maxwell模型进行推广。
推广的方法是唯象的。
在唯象方法中,强调建立描述应力分量与形变分量或形变率分量间正确关系的方程,而对材料的物质结构和其他性质不作深究。
下面介绍几种广义Maxwell模型。
2.3.1White-Metzner模型
该模型的主要特点是在Maxwell模型方程(3-4)中,采用对应力张量求Oldroyd随流微商代替一般偏微商。
随流坐标系(convectedframeofreferenee)。
对于纯粘性流体,由于无记忆特性,应力只依赖于形变速率的瞬时值,因此采用固定空间坐标系计算是方便的。
对于粘弹性流体,其应力不仅依
赖于即时形变,还依赖于形变历史,流体元有“记忆”能力,因此采用固定空间坐标系描述就很麻烦。
另外在固定坐标系中考察流动时,材料元的形变往往总与平动、转动牵扯在一起,讨论也不方便。
为此,人们采用一种镶嵌在所考察的材料元上,随材料元一起运动的坐标系作为参照系,称为随流坐标系,在此参照系中考察流体元的形变。
这种参照系最初是由Oldroyd提出的。
由于在随流坐标系中定义的任何形变的度量总是针对同一个材料元的,可摆脱平动和转动速率的影响,故讨论流体元的形变问题有明显的优越性。
重要的是,我们必须建立随流坐标系和固定空间坐标系中各种物理量之间的转换关系。
因为所有的实验仪器都安装在固定坐标系中,所有对流
体性质的测量也都在固定坐标系中进行,只有建立起随流坐标系和固定坐标系中各物理量之间的转换关系,才能将随流坐标系中讨论的结果转换到实验室系中加以验证,以确定本构方程的优劣。
随流坐标系中,质点的随流坐标不变,为常数,故采用随流坐标对流体元的描述为物质描述。
同样在随流坐标系中,对物理量求时间导数时保持随流坐标不变,因此对任何物理量所求的时间导数均为物质导数。
Oldroyd随流微商即其中一种,记作「。
按照上面的讨论,这种随流
微商需要转换到固定的空间坐标系中。
二阶应力张量
Tij的Oldroyd随流
微商转换到固定坐标系后的形式为:
Dvi
TjTjiT
tDtxk
kjjTik
Xk
(3-20)
式中等号右边第一项为
D3
TjTjVk
DtJtJk1
Tj
xk
(3-21)
即二阶应力张量在固定坐标系的物质微商,可以理解为在固定坐标系中观
察者见到的某一材料元的应力张量对时间的变化率。
第二、三项中含有速
度梯度丄的影响,速度梯度中含有形变率张量d和旋转速率张量3两xk
部分,它描述了材料元对于固定坐标系的有限形变和旋转运动。
White-Metzner推广经典的Maxwell模型,其方法就是在方程(3-4)中采用对应力张量求Oldroyd随流微商代替一般偏微商。
White-Metzner模型的方程形式为:
a
c12°d(3-22)
此公式在形式上虽然与方程(3-4)相仿,但物理意义不同。
在这儿应力张量的时间变化率是在随流坐标系中计算的,它与在固定的空间坐标系中所求的一般偏微商以及物质微商都不相同。
2.3.2DeWitt模型
另一种广义Maxwell模型一DeWitt模型,是在Maxwell方程中对应
力张量求时间微商这一项,用共旋随流微商(Jaumann微商)代替一般偏微商。
二阶应力张量的共旋随流微商变换到固定坐标系后的形式为:
式中等号右边第一项和的意义同(3-21)式’为二阶张量Tj在固定坐
1
(3-29)
(ViVk)
ik()
2xkxi
它代表了材料元对于固定坐标系的有限旋转。
DeWitt推广Maxwell模型,在Maxwell方程中用对偏应力张量求共
旋随流微商(Jaumann微商)代替般偏微商,得到的
DeWitt模型方程
形式为:
a
1a20d
(3-30)
Dt
式中
D_Dt
D
ijDtijikjkjkik
(3-31)
我们用
DeWitt模型处理一下稳态简单剪切流场中应力应变关系,以
检验模型的说明能力。
已知在(3-23)式描绘的稳态简单剪切流场中,旋转速度张量为
0/20
(3-32)
3/200
000
计算(3-31)式中的Jaumann微商:
与在计算Oldroyd微商中同样的
原因,首先确定(3-31)式中第一项等于零,即
Dt=°;第二、三项分
别为
0
/2
0
11
2131
ik
jk3?
彷
T
/2
0
0?
12
2232
0
0
0
13
2333
12
22
——
32
(3-33)
22
2
11
21
——
31
22
2
01
0
0
和
11
12
13
0
/20
jk
T
:
ik彷?
3
21
22
23
?
/2
00
31
32
33
0
00
—12—
110
(3-34)
22
T22:
210
22
32
310
22
这样,方程(3-30)可展开写成:
0—0
2
20㊁°°
000
(3-35)
这同样是
9个方程联立的方程组,解此方程组得到粘度和法向应力差
系数为:
12/0/112
2
1
1122/2201/122(3-36)
2
2
/2/122
2233/01/I1
这一结果使人惊奇,首先(
3-36)式描述了剪切粘度与法向应力差系
数的剪切速率依赖性:
当剪切速率0时,材料表现出常数粘度n0(牛
顿性)和常数法向应力差系数;当剪切速率升高,粘度和法向应力差系
数均趋于下降,呈现出剪切变稀行为。
其次,公式表明,这里描述的液体是弹性液体,10,20。
与大多数高分子液体的实验事实一致的还有,公式给出的第一法向应力差系数为正(>0),第二法向应力差系数为负(v0),第一法向应力差系数的值大于第二法向应力差系数的绝对值(!
|2)(对比图2-11)。
这些结果,都只是因为我们将Maxwell方程
由实验室系推广到共旋随流坐标系,并且对时间的微商采用了共旋随流微商而得到的。
基于上述唯象性推广的结果,为了使数学模型更精确,更全面的描述非牛顿型流体复杂流变性质,人们又陆续提出许多类似的更复杂的模型。
2.3.3其他类型的微分模型
Jeffreys模型,方程形式为:
(3-37)
a
a1()20d
模型的特点是在原始Maxwell模型基础上,引入对形变率张量的偏微
分,同时引入第二时间参数2,使材料常数成为三个:
1,2,0
Oldroyd模型,方程形式为:
模型的特点是将Jeffreys模型中的时间微商由一般偏微商推广为
Oldroyd随流微商。
材料常数保留为三个:
广义Jeffreys模型,方程形式为:
(3-39)
Dd2Dt
模型的特点是将Jeffreys模型中的时间微商由一般偏微商推广为
Jaumann微商(共旋随流微商)。
材料常数保留为三个:
Oldroyd八常数模型,方程形式为:
Da
1(_dT)
1(adda)0tra?
dv1tradI
Dd
Dt
22
22d2v2trd2I
(3-40)
该模型的时间微商采用
Jaumann微商,并设置了八个常数:
入
1,入2,
卩0,卩1,卩2,v1,v2,n0,用于考虑偏应力张量和应变张量的各种关系。
其中前七个常数的量纲均为时间(s)。
运用此模型确实可以描写非牛顿型
流体的可变粘度和法向应力差效应,方程适用的范围也比较宽广。
Oldroyd
三常数模型:
Williams
和Bird
对Oldroyd
条件:
1
2
13v1
2
2
23V2
从而得到如下方程:
(X
1
Da
(Dt)
CTdda
八常数模型中的常数提出了如下限制性
00
(3-41)
(3-42)
itradI
20d2理2d2-trd2I
02Dt3
此方程中只剩三个材料常数:
1,2,0,相对来说比较容易确定。
可以
看出,DeWitt模型只是上述模型的一个特例
(3-43)
2.3.4Maxwell模型的叠加
前面讨论的Maxwell模型是单一的“弹簧和粘壶串联”的力学模型,因而只给出一个松弛时间1。
事实上由于高分子液体结构和运动单元的多重性和复杂性,其松弛时间往往有多个,形成一个松弛时间谱。
由此可以提出另一种推广Maxwell模型的方法,即从材料结构特点的观点出发进行推广。
虽然考虑了材料的结构特点,但推广的思路仍是唯象的。
推广方法为用多个Maxwell模型并联代替单一Maxwell模型来描写材料的非线性流变性,见图3-3。
其中每一个Maxwell模型都可用方程(3-4)表示,有独自的松弛时间,相应于大分子链的一种运动模式。
由于各个Maxwell
模型并联,因此总模型上承受的总应力。
应当等于各个Maxwell模型上
的应力之和。
模型方程可写成:
(7ap
p1
式中入p,np均为材料常数。
在此基础上,后人还有把方程中偏应力张量对时间的一般偏微商推广为Oldroyd随流微商或Jaumann微商(共旋随
流微商),此外并无其他新观点
图3-3并联Maxwell模型
3.积分型本构方程
与微商型本构方程相对应,本构方程也可写成积分型式,两者等价。
这种积分可以对应变历史进行,也可以对连续变化的运动模式求积分。
3.1Bolzman叠加原理
原理表述:
对于时间序列中一系列阶跃应变(或应力)的输入,体系
在时刻t的应力(或应变)响应,可以表示为不同时刻t'(t 列个别响应的线性叠加。 设t为现在时刻,t'为过去的时刻序列(t从非常遥远的-=演进到现在时刻t)。 对一系列的过去应变e(t)体系在现在时刻的应力总响应,按照上述原理记为: a(t)2t(tt)e(t)dt(3-50) 公式中(tt')dt'dG(tt')(3-51) 为弹性模量的微分,而G(tt')称材料的松弛(模量)函数,它是一个随 时间间隔(tt')变化的弹性模量。 注意t。 dt 需要指出的是,Bolzmann叠加原理属于线性叠加理论,原则上只适用于小形变过程。 如果要对大形变过程也适用,必须加以推广或再加说明。 或者a、把问题变换到恰当的坐标系下去讨论(如选择在随流坐标系中讨论);或者b、假定对应于大应变过程,其分割的每一个子应变过程的应力响应足够小,小到还是可以进行线性叠加。 如把材料在一段历史中受 到的应变近似分割为若干小阶跃应变,每一个子应变过程对体系现在时刻的应力响应可以进行线性叠加,然后求和取极限(图3-4)。 图3-4应力-应变响应的Bolzmann线性叠加 3.2Maxwell模型的积分形式 这儿我们不加推导地给出Maxwell模型本构方程的积分形式,并证明它与微分型本构方程完全等价。 Maxwell模型本构方程的积分形式以Bolzmann叠加原理的形式写出: t c(t)2G(tt)d(t)dt(3-52) 虫 式中G(tt')』e(3-53) 为与时间过程相关的松弛弹性模量,其中°为常数粘度,为松弛时间,均为材料参数。 d(t')为过去时刻t'勺体系所受的形变率张量,(Tt为体系在现在时刻t的应力响应。 从(3-53)式可知,G(tt')是一个随时间过程衰减的弹性模量,表征着材料的弹性记忆能力随时间衰减。 可以证明,公式(3-52)与公式(3-4)描写的Maxwell方程的微分形式完全等价。 证明如下: 对公式(3-52)利用定积分的微分公式: —f(x)dxf(t)(3-54) dta 仃1 得到: 2G(0)d(t)(-)a(3-55) 将以上结果代入微分模型(3-4)式中, 1 证毕。 方程左方二a[2G(0)d(t)-a]=2G(O)d(t) =20dt方程右方 考虑到高分子材料的分子链具有多种运动模式,所以其流变性质也可 以采用一系列松弛时间不同的Maxwell模型并联而成的复杂模型来加以 描述(见图3-3)。 对于积分型本构方程,设应力张量的各分量与各模型应变张量的各分量线性相关,则总应力响应可以表示为若干个分响应的线性叠加,方程形式为: tt/ tN (3-56) at2^eid(t/)dt' i1i t 2m(tt')d(t')dt' 弛时间。 N (3-57) m(tt')=Gi i1 (t-F)称为材料的记忆函数,表示材料在遥远过去时刻(t'所承 受的应变,至今仍保留着的对现在时刻应力响应的贡献。 除Maxwell模型外,可以证明,对其他微分型本构模型,如 White-Metzner模型,DeWitt模型,Jefreys模型等,对应的都有其积分型本构模型,且两者等价。 限于篇幅及数学工具的缺陷,此处删略,有兴趣的读者可参阅有关专著。 4.流变模型对高分子科学和高分子工程问题的意义 对高分子液体流变本构方程理论和实验规律的研究对于促进高分子材料科学,尤其高分子物理的发展和解决聚合物工程中(包括聚合反应工程和聚合物加工工程)若干重要理论和技术问题都具有十分重要的意义。 一则由于高分子材料复杂的流变性质需要精确地加以描述,二则由于高新技术对聚合物制品的精密加工和完美设计提出越来越高的要求,因此 以往那些对材料流动性质的经验的定性的粗糙认识已远远不够。 众所周知,高分子结构研究(包括链结构、聚集态结构研究)以及这种结构与高分子材料作为材料使用时所体现出来的性能、功能间的关系研 究始终是高分子物理研究的主要线索。 与“静态”的结构研究相比,高分子“动态”结构的研究,诸如分子链运动及动力学行为、聚集态变化的动力学规律、高分子流体的非线性粘弹行为等,更是近年来引人注目的前沿领域。 按现代凝聚态物理学的概念,高分子体系被称为软物质(soft matter)或复杂流体(complexfluids)。 所谓软物质,即材料在很小的应变下就会出现强烈的非线性响应,表现出独特的形态选择特征。 这正是高分子流体的本征特点。 如果能精确描述出高分子液体的复杂应力-应变关 系,找出这种关系与材料的各级结构间的联系,无疑对高分子凝聚态理论的发展具有
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第三 粘弹性 流体 方程