粘弹塑性模型的基本概念

第七章 粘弹塑性模型的基本概念7 . 1 引言为了描述土体应力一应变关系受时间的影响，需要采用与时间有关的类模型（如粘弹胜模酬、粘塑性模型，粘弹塑隆模型）来描述土的性状。

弹性、塑性和粘性是连续介质的三种基本性质，各在定条件F 独自反映材料本构关系的一个方面的特性。

理想弹性模型、理想塑胜模型（或称刚塑性模型）和理想粘性模型是反映这三种性质的理想模型，通常称为简单模型。

实际工程材料的本构关系可以用这些简单模型的各种组合来构成。

理想弹性模型又称虎克弹性模型，通常用理想弹簧表示（图7-1( a )）。

其本构方程为虎克定律。

一维条件下，如单轴压缩和纯剪清况下，表达式分别为：E σε= （7.1.1）G τγ= （7.1.2）式中E —— 弹性模量、G ——剪切模量。

剪切模量与弹性模量和泊松比的关系如下式所示：()21E G ν=+ （7.1.3） 式中 ν ——泊松比。

三维条件下本构方程可表示为下述形式：m K νσε= （7.1.4）式中 K ——体积弹性模量。

（a ） （b ）图7-1 理想弹性模型体积弹性模量与弹性模量和泊松比的关系如下式所示：()312E K ν=- （7.1.6） 理想粘性模型又称牛顿粘滞体模型。

通常用一粘壶（或称阻尼器）表示（图7-2 ( a ) ）。

粘壶内充满粘滞液体和一个可移动的活塞。

活塞在粘滞液体中的移动速度与所受阻力成正比关系，反映了粘性介质内一点的应力与该点处应变速率成正比例关系的性质。

一维条件如单轴压缩或纯剪情况下，表达式分别为：σϕε= （7.1.7）τηγ= （7.1.8）式中 ϕ、η ——粘滞系数。

由上两式可以看出，从数学表达的形式上与理想弹性体单轴压缩和纯剪时的本构方程相类似。

与理想弹性体的方程相对应，类似式7.1.3，存在下述关系：()*21ϕην=+ （7.1.9）式中 *ν ——粘性应变速率的横向比值。

（a ） （b ）图7-2 理想粘性模型理想粘性体的体积变化与形状变化速率无关，即不具有体积粘性。

因此，*ν应等于0.5 。

于是式7.1.9成为：3ϕη= （）这与弹性不可压缩时的E=3G 相对应。

在三维条件下理想粘性体本构方程可表示为：2ij ij S e η= （）理想塑性模型又称Saint-Venant 塑性模型，或称刚塑性模型。

通常采用两块接触的粗糙面表示（图7-3 （a ））。

面上存在有一称晰脚擦阻力，与作用在面上的法向压力无关，是一常数。

若外作用力心婚此起始摩擦阻力，物体不发生变形。

一维条件如单轴压缩或此钾扮况，当轴向应力或剪应力小于某一数值时，物体不发生变形．当软祠应力或剪应力等于某数值时，物体产生流动，变形无限制增长．理想塑性模刮的体积应变等于零，即体积不发生改变。

在三维条件下理想塑性体的本构方程可表示为：（a ） （b ）图 7-3 理想塑性体模型当 ij ij S H <时，0ij e =当 ij ij S H =时，2ij ij S e λ= （）式中 ij H ——起始摩擦阻力，或称塑性条件；λ——比例常数。

式表明，理想塑性体的塑性应变偏量的变化率与应力偏量成正比。

由理想弹性模型、理想粘性模型和理想塑性模型等简单模型可以组合成许多复杂模型。

由理想弹性模型和理想塑性模型可以组合成理想弹塑性模型。

由弹性模型和粘性模型可以组合成各种粘弹性模型。

由粘性模型和塑性模型可以组合成各种粘塑性模型。

由弹性模型、粘性模型和塑性模型可以组合成各种粘弹塑性模型。

理想弹塑性模型已在第六章作了介绍。

在以下几节将对几种由简单模型组成的粘弹性模型、粘塑性模型和粘弹塑胜模型作简单介绍。

利用简单模型可以组合成各种复杂模型，从而可以建立各种材料的本构方程。

但是进一步的研究发现，许多材料的实际性状并不能满意地用简单的组合模型来描述，而目采用复杂的组合模型又常遇到数学上的困难。

因此，常常在试验的基础上，通过假设一实验一理论的方法建立材料的本构力程。

在本章的最后一节将简要介绍描述材料蠕变现象的蠕变力程。

7 . 2 粘弹性模型既具有弹性又具有粘性的性质称为粘弹性。

蠕变和应力松弛现象是人们熟悉的也是特别受重视的粘弹性胜质粘弹性性质的特点是在本构方程中除了有应力和应变项外，还包括有它们对时间导数的项。

对线性粘弹胜材料，其本构方程的一般表达式为：()()0101m n m n a a a b b b σσσεεε+++=+++ （7.2.1） 式中 ,i i a b ——与材料性质有关的参数。

下面首先介绍几种简单的粘弹性模型，然后再介绍较复杂的情况。

7.2.1Maxwell 模型Maxwell 模型又称松弛模型。

它是由线性弹簧和牛顿枯壶串联组成，如图7 -4 (a ）所示。

在串联条件下，作用在两元件上的应力相同，而总的应变应为两个元件应变的和，即εεε'''=+ （7.2.2）或εεε'''=+ （7.2.3）式中 ,εε'''——分别为线性弹簧和粘壶的应变；,εε'''——分别为线性弹簧和粘壶的应变率。

考虑到线性弹簧有/E εσ'=和牛顿粘壶有/εσϕ''=，则式7.2.3可改写成： E σσεϕ=+（7.2.4） （a ） （b ） （c ）图7-4 Maxwoll 模型写成如式7.2.1的标准形式，上式可改写为：n σσϕε+=（7.2.5）式中 n ——松驰时间，n E ϕ=，量纲为时间。

式7.2.5称为Maxwell 方程。

若物体获得初始应变0ε以后总应变保持不变（图7-4b) ，即0ε=，式7.2.5成为：0n σσ+=（7.2.6） 积分上式，得/t n Ce σ-= （7.2.7)式中 C ——积分常数。

应用初始条件，0t =，0σσ=代人式7.2.7解出C ，再代人式7.2.7 , 得 /0t n e σσ-= (7.2.8 )式7.2.8表示，Maxwell 模型在保持总应变不变的条件下，发生应力随时间衰减的松弛现象，如图7-4c 所示。

若物体获得初始应力0σ以后，保持应力不变，即0σ=0σϕε= （7.2.9 )式7.2.9表示材料应变率为常数，即应变随时间成比例地增长，因此变形随时间无限地发展。

下面讨论松弛试验的情况。

在松弛试验中，首先对试件施加应变0ε，然后保持应变为定值，进而测量作为时间函数的应力值，确定松弛规律。

松弛试验中应变可记为：()0u t εε= ( 式中 ()u t ——单位阶梯函数。

单位阶梯函数定义为：()1110,1,t t u t t t t <⎧-=⎨>⎩ （） 在松弛试验中10t =()1u t t -可表示为()u t 。

将式代人式7.2.5，得()E t n σσεδ+= （）式中 ()t δ——脉冲δ函数，()()d t u t dt δ=⎡⎤⎣⎦。

脉冲δ函数定义为：()0,0,0t t t δ≠⎧=⎨+∞=⎩ （） ()1t t dt δ-∞=⎰（） 脉冲δ函数具有下述性质，对于任何连续函数()f t ，当1t t >时，有 ()()()()111t f t d f t u t t τδττ-∞-=-⎰（）利用式 ()()/0t n t E e u t σε-= （）式表示Maxwell 模型的应力松弛规律，简记为：()()0t t σε=Φ （）式中 ()t Φ——松弛函数，其表达式为()()/t n t Ee u t -Φ= （）7.2.2 Kelvln 模型Kelvln 模型又称非松弛模型。

这种模型曾由W . Voigt 和Kelvin 提出，故又称为Voigt —Kelvin 模型。

它是由线性弹簧和牛顿粘壶并联组成，如图7-5 (a ）所示。

在并联条件下，两个元件的应变相同，而总的应力应为两个元件的应力之和，即E σσσεϕε'''=+=+ （） 若在0t =时，瞬时地加上应力0σσ=，并保持不变，则由式可得0E ϕεεσ+=积分上式，得()01t e E λσε-=- （） 式中 λ——衰减系数，1E n λϕ==； n ——滞后时间。

（a ） （b ）图7-5 Kelvln 模型由式可知，当t →∞，应变趋于个稳定值0/E σ。

若物体获得初始弹性应变0ε之后保持应变不变，即0ε=。

由式得0E σε==常量 （）上式表明在这种情况下应力不衰减。

下面讨论蠕变试验的情况。

在蠕变试验中，首先对试件施加应力0σ，然后保持应力为定值来量取作为时间函数的应变值。

若取瞬时加载的时刻为0t =，则加载过程可表示为：()0u t σσ= （） 式中 ()u t ——单位阶梯函数。

将式()0u t σελεϕ+= （） 注意到单位阶梯函数有如下性质()()()()111ttt f u t d u t f d ττττττ-∞-=-⎰⎰ 此处τ为积分变量。

积分式，得()()()01tt e u t E λσε-=- （） 式中1E n λϕ== 式表示Kelvin 模型的蠕变规律，可简记为：()()t t εσ=ψ式中 ()t ψ——蠕变函数。

蠕变函数的表达式为()()()11t t e u t Eλ-ψ=- （） 7.2.3 三元件粘弹性模型 图7-6a 表示个三元件粘弹性模型。

它是由线性弹簧和Kelvin 模型串联组成，包括二个线性弹簧和一个牛顿粘壶，共三个元件，故称三元件粘弹性模型。

用ε''表Kelvin 模型的应变，ε'表示与Kelvin 模型串联的线性弹簧的应变，σ'表示Kelvin 模型中线性弹簧中的应力，σ''表示牛顿粘壶中的应力，σ和ε分别表示总应力和总应变。

分析各元件的应力或应变相互间关系，不难得到下列各式：εεε'''=+ （）σσσ'''=+ （）E σε''= （）E σε'''''= （）σϕε''''= （）式中 E '——与Kelvin 模型串联的线性弹簧的弹性模量；E ''——Kelvin 模型中线性弹簧的弹性模量；ϕ——牛顿粘壶的粘滞系数。

结合式()E E E E E σϕσεϕε'''''''++=+ （）式还可改写为：n nH E σσεε+=+ （）式中n E E ϕ='''+ （）图7-6 三元件粘弹性模型H E '= ()E E E E E '''='''+ () 若物体作用有初始应力σ，且保持不变，即0σ=，且在0t =时，/H εσ=。