第七章粘弹塑性模型的基本概念

1、第七章粘弹塑性模型的基本概念7 . 1 引言为了描述土体应力一应变关系受时间的影响，需要采用与时间有关的类模型（如粘弹胜模酬、粘塑性模型，粘弹塑隆模型）来描述土的性状。弹性、塑性和粘性是连续介质的三种基本性质，各在定条件 F独自 反映材料本构关系的一个方面的特性。理想弹性模型、理想塑胜模型（或称刚塑 性模型）和理想粘性模型是反映这三种性质的理想模型，通常称为简单模型。实 际工程材料的本构关系可以用这些简单模型的各种组合来构成。理想弹性模型又称虎克弹性模型，通常用理想弹簧表示（图7-1（ a ）。其本构方程为虎克定律。一维条件下，如单轴压缩和纯剪清况下， 表达式分别为：E（）G（ ）式中E弹性模

2、量、G剪切模量。剪切模量与弹性模量和泊松比的关系如下式所示：GE（ ）2 1式中泊松比。)三维条件下本构方程可表示为下述形式：m K式中K 体积弹性模量。(a)(b)图7-1理想弹性模型体积弹性模量与弹性模量和泊松比的关系如下式所示:)理想粘性模型又称牛顿粘滞体模型。通常用一粘壶(或称阻尼器)表示(图 7-2 ( a )。粘壶内充满粘滞液体和一个可移动的活塞。活塞在粘滞液体中的 移动速度与所受阻力成正比关系，反映了粘性介质内一点的应力与该点处应变速 率成正比例关系的性质。一维条件如单轴压缩或纯剪情况下，表达式分别为：&()&()式中、粘滞系数 由上两式可以看出，从数学表达的形式

3、上与理想弹性体单轴压缩和纯剪时的本构 方程相类似。与理想弹性体的方程相对应，类似式 ，存在下述关系：)式中粘性应变速率的横向比值。(a)图7-2理想粘性模型理想粘性体的体积变化与形状变化速率无关,即不具有体积粘性。因此， 应等于。于是式成为:)这与弹性不可压缩时的E=3G相对应。在三维条件下理想粘性体本构方程可表示为:Sj2 &)理想塑性模型又称Saint-Venant塑性模型，或称刚塑性模型。通常采用 两块接触的粗糙面表示（图7-3 （a）。面上存在有一称晰脚擦阻力，与作用 在面上的法向压力无关，是一常数。若外作用力心婚此起始摩擦阻力，物体不发 生变形。一维条件如单轴压缩或此钾扮况，

4、当轴向应力或剪应力小于某一数值时， 物体不发生变形.当软祠应力或剪应力等于某数值时， 物体产生流动，变形无限 制增长.理想塑性模刮的体积应变等于零， 即体积不发生改变。在三维条件下理 想塑性体的本构方程可表示为：(a)(b)图7-3理想塑性体模型当qm时，列o当Sj可时,Sj 2 & J(7.1.12 )式中Hj 起始摩擦阻力，或称塑性条件;比例常数。式表明，理想塑性体的塑性应变偏量的变化率与应力偏量成正比由理想弹性模型、理想粘性模型和理想塑性模型等简单模型可以组合成许多复杂模型。由理想弹性模型和理想塑性模型可以组合成理想弹塑性模型。由弹性 模型和粘性模型可以组合成各种粘弹性模型。 由

5、粘性模型和塑性模型可以组合成 各种粘塑性模型。 由弹性模型、 粘性模型和塑性模型可以组合成各种粘弹塑性模 型。理想弹塑性模型已在第六章作了介绍。 在以下几节将对几种由简单模型组成 的粘弹性模型、粘塑性模型和粘弹塑胜模型作简单介绍。利用简单模型可以组合成各种复杂模型， 从而可以建立各种材料的本构 方程。但是进一步的研究发现， 许多材料的实际性状并不能满意地用简单的组合 模型来描述， 而目采用复杂的组合模型又常遇到数学上的困难。 因此，常常在试 验的基础上， 通过假设一实验一理论的方法建立材料的本构力程。 在本章的最后 一节将简要介绍描述材料蠕变现象的蠕变力程。7 . 2 粘弹性模型既具有弹性又具

6、有粘性的性质称为粘弹性。蠕变和应力松弛现象是人们熟 悉的也是特别受重视的粘弹性胜质粘弹性性质的特点是在本构方程中除了有应 力和应变项外， 还包括有它们对时间导数的项。 对线性粘弹胜材料， 其本构方程 的一般表达式为：mna0 a1 & L amb0 b1& L bn（7.2.1 ）式中 ai ,bi 与材料性质有关的参数。下面首先介绍几种简单的粘弹性模型，然后再介绍较复杂的情况。7.2.1Maxwell 模型Maxwell 模型又称松弛模型。它是由线性弹簧和牛顿枯壶串联组成，如图 7 -4 （a ）所示。在串联条件下，作用在两元件上的应力相同，而总的应变应为两 个元件应变的和，

7、即（ 7.2.2 ） 或& & &（ 7.2.3 ）式中 , 分别为线性弹簧和粘壶的应变；&, & 分别为线性弹簧和粘壶的应变率。考虑到线性弹簧有& &E和牛顿粘壶有&/ ，则式723可改写成:(7.2.4 )OEIT-(a)(b)图 7-4 Maxwoll 模型(c)写成如式的标准形式，上式可改写为:n& & () 式中n 松驰时间，n己，量纲为时间。式7.2.5 称为Maxwell方程。若物体获得初始应变 °以后总应变保持不变（图7-4b），即& 0,式成为：n& 0（）积分上式，得C

8、e t/n（ 7.2.7）式中C积分常数。应用初始条件，t 0，0代人式解出C，再代人式,得t/n0e式表示，Maxwell模型在保持总应变不变的条件下，发生应力随时间衰减 的松弛现象，如图7-4c所示若物体获得初始应力 °以后，保持应力不变，即 & 0,贝U式成为:0&（7.2.9 ）式表示材料应变率为常数,即应变随时间成比例地增长，因此变形随时间无限地发展。下面讨论松弛试验的情况。在松弛试验中，首先对试件施加应变0 ,然后保持应变为定值，进而测量作为时间函数的应力值,确定松弛规律。松弛试验中应变可记为：out式中 u t单位阶梯函数。单位阶梯函数定义为：0,tt1

9、（）ut t11,tt1在松弛试验中t1 0 u t t1可表示为u t。将式代人式，得& 一 nE t（）式中t脉冲函数，d tu t 。dt脉冲函数定义为：t0,t0（）,t0ttdt 1（）脉冲 函数具有下述性质，对于任何连续函数ft,当t tl时，有tfti d f ti u t ti（7.2.15 ）利用式，积分式，可得t E 0e t/nu t（ ）式表示Maxwell模型的应力松弛规律，简记为：t t o（ ）式中 t 松弛函数，其表达式为(7218 )t/ nt Ee u t722 Kelvl n模型Kelvln模型又称非松弛模型。这种模型曾由W . Voigt和Kel

10、vin提出，故又称为Voigt Kelvin模型。它是由线性弹簧和牛顿粘壶并联组成，如图7-5(a )所示。在并联条件下，两个元件的应变相同，而总的应力应为两个元件的应 力之和，即(7219)若在t 0时，瞬时地加上应力0，并保持不变，则由式可得()积分上式，得式中(a)滞后时间。rLaaf衰减系数,，应变趋于个稳定值0/E。由式可知，当t若物体获得初始弹性应变0之后保持应变不变，即& 0。由式得上式表明在这种情况下应力不衰减。下面讨论蠕变试验的情况。在蠕变试验中，首先对试件施加应力0，然后保持应力为定值来量取作为时间函数的应变值。若取瞬时加载的时刻为t 0 ,则加载过程可表示为：oU

11、 t（7223）式中 u t单位阶梯函数。将式代人式，得&u t（）注意到单位阶梯函数有如下性质ttfut| d ut1ft1d（）此处为积分变量。积分式，得t0 1Et丄e u t（）式中1 E_n式 7.2.26 表示 Kelvin模型的蠕变规律，可简记为：tt（）式中t蠕变函数。蠕变函数的表达式为1t - 1Et.e u t（）三元件粘弹性模型图7-6a表示个三元件粘弹性模型。它是由线性弹簧和Kelvin模型串联组成，包括二个线性弹簧和一个牛顿粘壶，共三个元件，故称三元件粘弹性模型。用表Kelvin模型的应变，表示与Kelvin模型串联的线性弹簧的应变，表示Kelvin模型中线性

12、弹簧中的应力，表示牛顿粘壶中的应力，和分别表示总应力和总应变。分析各元件的应力或应变相互间关系，不难得到下列各式：（）(7230 )(7231 )E(7.2.32 )(7.2.33 )式中E 与Kelvin模型串联的线性弹簧的弹性模量；E Kelvin模型中线性弹簧的弹性模量；牛顿粘壶的粘滞系数。结合式至式各式，消去组成元件中的应力和应变，得E E& EE E &()式还可改写为：n& nH & E( )式中n()E ErLiia厂壬图7-6三元件粘弹性模型H EE EE E E若物体作用有初始应力，且保持不变，即&于是，由式可求得应变的变化规律为：H

13、E 彳 1 eH HE(7237)(7238)0 ,且在t 0时,/H。Et/Hn(7239 )上式表示的应变随时间的变化规律如图 7-6 (b )所示。图中应变起始值为 /H , 最终值为/ E，其应变速率由起始时的最大值逐渐趋于零。0后总应变保持不变，即0，& o且在t 0时,若物体获得初始弹性应变H 0。于是，由式可求得应力随时间的变化规律为:Et/n0 H E °e(7240)上式表示的应力变化规律如图7-6(b)所示。由图可以看到，物体中的应力从最初的H 0衰减到最终值E 0若物体初始时作用有应力以后随时间变化作用有应力t。根据叠加原理，由式可以得到在时刻t时物体的

14、变形,H E 1 0 1 H HEe Et/HnH E E t /HnHE(7241)对上式右端进行分部积分,则式可改写为H2nH EH2n/Hn(7242)(7243)(7244)式通常称为线性遗传方程。式中H称为瞬时弹性模量,称为遗传函数，它表示在 时刻作用的应力对时刻t的变形的影响。三元件粘弹性模型除了上述介绍的基本形式外，还有其它组成方式的三元件 粘弹性模型。如由Maxwell模型与一个粘壶并联组成，或由一个粘壶与Kelvin模 型串联组成。这些形式的本构方程读者自己不妨加以推导。广义Maxwell模型和广义Kelvin 模型增加组成模型的元件数，可以得到更为复杂的模型应用得较多的是广

15、义Maxwell模型和广义Kelvin模型。rEi图7-7广义Maxwell模型广义Maxwell模型是由一个线性弹簧和一系列Maxwell 模型 并联而成，如图7-7所示。若t 0时模型获得单位弹性应变0 1后，保持总应变不变，模型中的应力随时间的变化应等于各简单模型之和，即(7245 )G tE0Ee t/nii式中G t松弛弹性模量，等于单位总应变所对应的应力;ni 松弛时间，nii/Ei若模型的应变可用表示.其本构力程可由叠加原理得到,利用分部积分法,上式可改写为:oG(7246 )(7247 )dG td t上式又可简写为:(7248 )式中R tdG tr4LifHI图7-8广义K

16、elvln模型广义Kelvin模型是由一个Maxwell模型和一系列Kelvin 模型串联而成,如图7-8所示。若t 0时模型受到单位应力1后保持不变，它的总应变等于各个简单模型的应变之和，即EoiL"(7249)式中 J t 蠕变柔度，等于单位应力引起的应变;衰减系数，i Ei/ i，其倒数为延迟时间若模型的应力用 t表示，其本构力程可由叠加原理得到,t0J t I t d d00d(7250)利用分部积分法，上式可改写为:dJ ttttK tdH0(7252 )(7251 )记J 0右,皆K t这样就得到了与式7244相同的线性遗传力程，粘塑性模型既具有粘性又具有塑性性质称为粘塑

17、性。粘塑性体在荷载作用下，当应力达 到某临界值时，屈服和流动现象发生，其变形速率与物体的粘性有关。材料的粘 塑性可由粘性元件(粘壶)和塑性元件(摩擦件)组合来描述。Bingham 模型是由理想刚塑性模型和牛顿粘壶并联而成，如图7-9 (a)所示。显然，Bingham模型只有当应力达到屈服极限时，才开始变形。在此以前 表现为刚性，屈服以后，呈现出粘塑性性质。其本构关系为：s &()当 s时，& 0,物体不发生变形。当s时，由式，得(b)(a)图 7-9 Bingham 模型对Bingham模型，应力s时，应变为零。如应力s时，应力可由式确定，而应变无限地增大。7.4 粘弹塑性模型

18、粘弹塑性是包含了弹性、粘性和塑性三力面的性质。粘弹塑性可以由弹簧、 粘壶和摩擦元件的各种组合来描述。下面简略介绍一个三元件粘弹塑性模型。图7-10表示一个三元件粘弹塑性模型，由线性弹簧、牛顿粘壶和一 个摩擦件组成。首先考虑线性强化情况，然后再分析理想粘弹塑性情况。 对这一 模型，总的应变为：e vp(741)式中 e 弹性应变;vp粘塑性应变。弹簧中应力与总的应力相等，即(742)摩擦件中应力 p取决于是否已经达到屈服应力s，可表示为:vp(743)(744)式中B 强化参数，定义为：B滸Et/(745)式中Et 切线模量。s时，还有vp)图7-10三元件粘弹塑性模型结合式744和式，得vpB

19、 vp(747)结合式、式和式，得BE(748)1记 一，称为介质流动参数，则式可改写为:s B vp因此，粘塑性应变率为:&ps b vp(7410)（称为“过式表明粘塑性应变率是由超过稳态屈服应力的那部分应力值 应力”）所决定的。若作用于模型的应力为常值A时，即/to，式可改写为:(7411 )式的解为:exp(7412)_AA sE B(7413)对于理想粘塑性材料，B 0,利用罗比达法则，式可改写为W A s t对于更复杂的弹粘塑性模型读者可参阅有关专着，这里不作进步介绍了。7.5 蠕变物体的蠕变现象可以采用由一定数量的弹性、 粘性和塑性元件组成的模型来 描述。但元件多了，计算

20、相当复杂，且其关系不容易由试验确定。在实际应用中， 常常直接由试验来确定应力、应变和时间之间的关系。下面简单介绍几种主要的 蠕变方程的形式。1 老化理论老化理论假设蠕变应变与应力、时间之间具有某种函数关系，即c f ,t(7.5.1 )式中c蠕变应变物体的总应变将由三部分组成，即(7.5.2 )式中 c, p 分别为弹性应变和塑性应变当应力未超过屈服应力时,0，则式可改写为:,t(7.5.3 )0时间r »nun30604011120008060o O4 27-11 )，即hap clogt(7.5.4 )图7-11应力-应变-时间实测曲线Buisman（ 1936）根据大量的试验资

21、料，认为饱和粘土的时间与 沉降关系在半对数坐标上呈线性关系（图式中沉降;ha 土层或试样起始高度;竖向应力增量;p主固结系数；时间效应系数。式可改写成:(7.5.5 )a log t2 流动理论流动理论认为蠕变应变速率与应力、时间之间存在某种简单的 函数关系，即,t(7.5.6 )式中蠕变应变速率。物体的总应变速率可表示为:(7.5.7 )& & &p &式中 &e &p分别表示弹性应变速率和塑性应变速率。试验资料表明，正常固结粘土和超固结粘土，在排水或不排水条件下，其应变速率和蠕变时间的关系在全对数坐标上呈线性关系。应变速率与应力的关系也是线性关系（图7-12 ）。应变速率与时问的关系为：In & In &t1,q mln t1(7.5.8 )式中 &应变速率;& bq 单位时间的应变速率，为应力 q的函数;t1 单位时间（即1分钟）；m关系曲线的坡度，一般为1。应变速率与应力的关系可表示为:ln & ln &1，q。(7.5.9 )式中&U时应变速率，为蠕变时间t的函数;q偏应力，q i关系曲线的坡度。式也可