《第七章-粘弹塑性模型的基本概念(20页).doc》由会员分享，可在线阅读，更多相关《第七章-粘弹塑性模型的基本概念(20页).doc（17页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-第 1 页第七章-粘弹塑性模型的基本概念-第 2 页第七章 粘弹塑性模型的基本概念7.1 引言为了描述土体应力一应变关系受时间的影响，需要采用与时间有关的类模型（如粘弹胜模酬、粘塑性模型，粘弹塑隆模型）来描述土的性状。弹性、塑性和粘性是连续介质的三种基本性质，各在定条件 F 独自反映材料本构关系的一个方面的特性。理想弹性模型、理想塑胜模型（或称刚塑性模型）和理想粘性模型是反映这三种性质的理想模型，通常称为简单模型。实际工程材料的本构关系可以用这些简单模型的各种组合来构成。理想弹性模型又称虎克弹性模型，通常用理想弹簧表示（图 7-1(a)）。其本构方程为虎克定律。一维条件下，如单轴压缩和纯剪清
2、况下，表达式分别为：E（7.1.1）G（7.1.2）式中 E 弹性模量、G剪切模量。剪切模量与弹性模量和泊松比的关系如下式所示：2 1EG（7.1.3）式中泊松比。三维条件下本构方程可表示为下述形式：mK（7.1.4）式中K体积弹性模量。（a）（b）图 7-1 理想弹性模型体积弹性模量与弹性模量和泊松比的关系如下式所示：3 1 2EK（7.1.6）-第 3 页理想粘性模型又称牛顿粘滞体模型。通常用一粘壶（或称阻尼器）表示（图7-2(a)）。粘壶内充满粘滞液体和一个可移动的活塞。活塞在粘滞液体中的移动速度与所受阻力成正比关系，反映了粘性介质内一点的应力与该点处应变速率成正比例关系的性质。一维条件
3、如单轴压缩或纯剪情况下，表达式分别为：（7.1.7）（7.1.8）式中、粘滞系数。由上两式可以看出，从数学表达的形式上与理想弹性体单轴压缩和纯剪时的本构方程相类似。与理想弹性体的方程相对应，类似式 7.1.3，存在下述关系：*2 1（7.1.9）式中*粘性应变速率的横向比值。（a）（b）图 7-2 理想粘性模型理想粘性体的体积变化与形状变化速率无关，即不具有体积粘性。因此，*应等于 0.5。于是式 7.1.9 成为：3（7.1.10）这与弹性不可压缩时的 E=3G 相对应。在三维条件下理想粘性体本构方程可表示为：2ijijSe（7.1.11）理想塑性模型又称 Saint-Venant 塑性模型
4、，或称刚塑性模型。通常采用两块接触的粗糙面表示（图 7-3（a）。面上存在有一称晰脚擦阻力，与作用在面上的法向压力无关，是一常数。若外作用力心婚此起始摩擦阻力，物体不发生变形。一维条件如单轴压缩或此钾扮况，当轴向应力或剪应力小于某一数值时，物体不发生变形当软祠应力或剪应力等于某数值时，物体产生流动，变形无限-第 4 页制增长理想塑性模刮的体积应变等于零，即体积不发生改变。在三维条件下理想塑性体的本构方程可表示为：（a）（b）图 7-3 理想塑性体模型当ijijSH时，0ije 当ijijSH时，2ijijSe（7.1.12）式中ijH起始摩擦阻力，或称塑性条件；比例常数。式 7.1.12 表明
5、，理想塑性体的塑性应变偏量的变化率与应力偏量成正比。由理想弹性模型、理想粘性模型和理想塑性模型等简单模型可以组合成许多复杂模型。由理想弹性模型和理想塑性模型可以组合成理想弹塑性模型。由弹性模型和粘性模型可以组合成各种粘弹性模型。由粘性模型和塑性模型可以组合成各种粘塑性模型。由弹性模型、粘性模型和塑性模型可以组合成各种粘弹塑性模型。理想弹塑性模型已在第六章作了介绍。在以下几节将对几种由简单模型组成的粘弹性模型、粘塑性模型和粘弹塑胜模型作简单介绍。利用简单模型可以组合成各种复杂模型，从而可以建立各种材料的本构方程。但是进一步的研究发现，许多材料的实际性状并不能满意地用简单的组合模型来描述，而目采用
6、复杂的组合模型又常遇到数学上的困难。因此，常常在试验-第 5 页的基础上，通过假设一实验一理论的方法建立材料的本构力程。在本章的最后一节将简要介绍描述材料蠕变现象的蠕变力程。7.2 粘弹性模型既具有弹性又具有粘性的性质称为粘弹性。蠕变和应力松弛现象是人们熟悉的也是特别受重视的粘弹性胜质粘弹性性质的特点是在本构方程中除了有应力和应变项外，还包括有它们对时间导数的项。对线性粘弹胜材料，其本构方程的一般表达式为：0101mnmnaaabbb（7.2.1）式中,iia b与材料性质有关的参数。下面首先介绍几种简单的粘弹性模型，然后再介绍较复杂的情况。7.2.1Maxwell 模型Maxwell 模型又
7、称松弛模型。它是由线性弹簧和牛顿枯壶串联组成，如图 7-4(a）所示。在串联条件下，作用在两元件上的应力相同，而总的应变应为两个元件应变的和，即（7.2.2）或（7.2.3）式中,分别为线性弹簧和粘壶的应变；,分别为线性弹簧和粘壶的应变率。考虑到线性弹簧有/E 和牛顿粘壶有/，则式 7.2.3 可改写成：E（7.2.4）（a）（b）（c）图 7-4 Maxwoll 模型-第 6 页写成如式 7.2.1 的标准形式，上式可改写为：n（7.2.5）式中n松驰时间，nE，量纲为时间。式 7.2.5 称为 Maxwell 方程。若物体获得初始应变0以后总应变保持不变（图 7-4b)，即0，式 7.2.
8、5成为：0n（7.2.6）积分上式，得/t nCe（7.2.7)式中C积分常数。应用初始条件，0t，0代人式 7.2.7 解出C，再代人式 7.2.7,得/0t ne(7.2.8)式 7.2.8 表示，Maxwell 模型在保持总应变不变的条件下，发生应力随时间衰减的松弛现象，如图 7-4c 所示。若物体获得初始应力0以后，保持应力不变，即0，则式 7.2.5 成为：0（7.2.9)式 7.2.9 表示材料应变率为常数，即应变随时间成比例地增长，因此变形随时间无限地发展。下面讨论松弛试验的情况。在松弛试验中，首先对试件施加应变0，然后保持应变为定值，进而测量作为时间函数的应力值，确定松弛规律。
9、松弛试验中应变可记为：0u t(7.2.10)式中 u t单位阶梯函数。单位阶梯函数定义为：1110,1,ttu tttt（7.2.11）在松弛试验中10t 1u tt可表示为 u t。-第 7 页将式 7.2.10 代人式 7.2.5，得 Etn（7.2.12）式中 t脉冲函数，dtu tdt。脉冲函数定义为：0,0,0ttt（7.2.13）1tt dt（7.2.14）脉冲函数具有下述性质，对于任何连续函数 f t，当1tt时，有 111tftdf t u tt （7.2.15）利用式 7.2.15，积分式 7.2.12，可得 /0t ntEeu t（7.2.16）式 7.2.16 表示 M
10、axwell 模型的应力松弛规律，简记为：0tt（7.2.17）式中 t松弛函数，其表达式为 /t ntEeu t（7.2.18）7.2.2 Kelvln 模型Kelvln 模型又称非松弛模型。这种模型曾由 W.Voigt 和 Kelvin 提出，故又称为 VoigtKelvin 模型。它是由线性弹簧和牛顿粘壶并联组成，如图 7-5(a）所示。在并联条件下，两个元件的应变相同，而总的应力应为两个元件的应力之和，即E（7.2.19）若在0t 时，瞬时地加上应力0，并保持不变，则由式 7.2.19 可得0E（7.2.20）积分上式，得01teE（7.2.21）式中衰减系数，1En；-第 8 页n滞
11、后时间。（a）（b）图 7-5 Kelvln 模型由式 7.2.21 可知，当t ，应变趋于个稳定值0/E。若物体获得初始弹性应变0之后保持应变不变，即0。由式 7.2.19 得0E常量（7.2.22）上式表明在这种情况下应力不衰减。下面讨论蠕变试验的情况。在蠕变试验中，首先对试件施加应力0，然后保持应力为定值来量取作为时间函数的应变值。若取瞬时加载的时刻为0t，则加载过程可表示为：0u t（7.2.23）式中 u t单位阶梯函数。将式 7.2.23 代人式 7.2.19，得 0u t（7.2.24）注意到单位阶梯函数有如下性质 111tttfutdutfd（7.2.25）此处为积分变量。积分
12、式 7.2.24，得 01tteu tE（7.2.26）式中1En式 7.2.26 表示 Kelvin 模型的蠕变规律，可简记为：tt（7.2.27）式中 t蠕变函数。蠕变函数的表达式为 11tteu tE（7.2.28）7.2.3 三元件粘弹性模型-第 9 页图 7-6a 表示个三元件粘弹性模型。它是由线性弹簧和 Kelvin 模型串联组成，包括二个线性弹簧和一个牛顿粘壶，共三个元件，故称三元件粘弹性模型。用表 Kelvin 模型的应变，表示与 Kelvin 模型串联的线性弹簧的应变，表示Kelvin 模型中线性弹簧中的应力，表示牛顿粘壶中的应力，和分别表示总应力和总应变。分析各元件的应力或
13、应变相互间关系，不难得到下列各式：（7.2.29）（7.2.30）E（7.2.31）E（7.2.32）（7.2.33）式中E与 Kelvin 模型串联的线性弹簧的弹性模量；EKelvin 模型中线性弹簧的弹性模量；牛顿粘壶的粘滞系数。结合式 7.2.29 至式 7.2.33 各式，消去组成元件中的应力和应变，得EEE EE（7.2.34）式 7.2.34 还可改写为：nnHE（7.2.35）式中nEE（7.2.36）图 7-6 三元件粘弹性模型HE(7.2.37)E EEEE(7.2.38)若物体作用有初始应力，且保持不变，即0，且在0t 时，/H。于是，由式 7.2.35 可求得应变的变化规
14、律为：/1Et HnHEeHHE（7.2.39）上式表示的应变随时间的变化规律如图 7-6(b）所示。图中应变起始值为/H，最终值为/E，其应变速率由起始时的最大值逐渐趋于零。若物体获得初始弹性应变0后总应变保持不变，即0，0且在0t 时，0H。于是，由式 7.2.35 可求得应力随时间的变化规律为：/00t nEHEe（7.2.40）-第 10 页上式表示的应力变化规律如图 7-6（b）所示。由图可以看到，物体中的应力从最初的0H衰减到最终值0E。若物体初始时作用有应力0，以后随时间变化作用有应力 t。根据叠加原理，由式 7.2.39 可以得到在时刻t时物体的变形，/000111tE tHn
15、Et HnHEHEdeedHHEHHEd（7.2.41）对上式右端进行分部积分，得/020tE tHnHEedHH n（7.2.42）记/2E tHnHEeK tH n（7.2.43）则式 7.2.42 可改写为 0ttK tdH（7.2.44）式 7.2.44 通常称为线性遗传方程。式中 H 称为瞬时弹性模量，K t称为遗传函数，它表示在时刻作用的应力对时刻t的变形的影响。三元件粘弹性模型除了上述介绍的基本形式外，还有其它组成方式的三元件粘弹性模型。如由Maxwell 模型与一个粘壶并联组成，或由一个粘壶与 Kelvin 模型串联组成。这些形式的本构方程读者自己不妨加以推导。7.2.4 广义
16、 Maxwell 模型和广义 Kelvin 模型增加组成模型的元件数，可以得到更为复杂的模型应用得较多的是广义Maxwell 模型和广义 Kelvin 模型。图 7-7 广义 Maxwell 模型广义 Maxwell 模型是由一个线性弹簧和一系列 Maxwell 模型并联而成，如图 7-7 所示。若0t 时模型获得单位弹性应变01后，保持总应变不变，模型中的应力随时间的变化应等于各简单模型之和，即/0it niiG tEEe（7.2.45）式中 G t松弛弹性模量，等于单位总应变所对应的应力；-第 11 页in松弛时间，/iiinE。若模型的应变可用 t表示其本构力程可由叠加原理得到，00td
17、tG tG tdd（7.2.46）利用分部积分法，上式可改写为：00tdG ttt Gdd t（7.2.47）上式又可简写为：0ttt ER td（7.2.48）式中dG tR td t；图 7-8 广义 Kelvln 模型广义 Kelvin 模型是由一个 Maxwell 模型和一系列 Kelvin 模型串联而成，如图 7-8 所示。若0t 时模型受到单位应力1后保持不变，它的总应变等于各个简单模型的应变之和，即 001111itiiJ tteEE（7.2.49）式中 J t蠕变柔度，等于单位应力引起的应变；i衰减系数，/iiiE，其倒数为延迟时间。若模型的应力用 t表示，其本构力程可由叠加原
18、理得到，00tdtJ tJ tdd（7.2.50）利用分部积分法，上式可改写为：00tdJ ttt Jdd t（7.2.51）记 10JH，dJ tK td t这样就得到了与式 7.2.44 相同的线性遗传力程，0tttK tdH（7.2.52）-第 12 页7.3 粘塑性模型既具有粘性又具有塑性性质称为粘塑性。粘塑性体在荷载作用下，当应力达到某临界值时，屈服和流动现象发生，其变形速率与物体的粘性有关。材料的粘塑性可由粘性元件（粘壶）和塑性元件（摩擦件）组合来描述。Bingham 模型是由理想刚塑性模型和牛顿粘壶并联而成，如图 7-9(a）所示。显然，Bingham 模型只有当应力达到屈服极限
19、时，才开始变形。在此以前表现为刚性，屈服以后，呈现出粘塑性性质。其本构关系为：s（7.3.1）当s时，0，物体不发生变形。当s时，由式 7.3.1，得s（7.3.2）（a）（b）图 7-9 Bingham 模型对 Bingham 模型，应力s时，应变为零。如应力s时，应力可由式 7.3.1确定，而应变无限地增大。7.4 粘弹塑性模型粘弹塑性是包含了弹性、粘性和塑性三力面的性质。粘弹塑性可以由弹簧、粘壶和摩擦元件的各种组合来描述。下面简略介绍一个三元件粘弹塑性模型。图 7-10 表示一个三元件粘弹塑性模型，由线性弹簧、牛顿粘壶和一个摩擦件组成。首先考虑线性强化情况，然后再分析理想粘弹塑性情况。对
20、这一模型，总的应变为：evp（7.4.1）式中e弹性应变；vp粘塑性应变。弹簧中应力与总的应力相等，即eeE（7.4.2）摩擦件中应力p取决于是否已经达到屈服应力s，可表示为：peps（7.4.3）-第 13 页pvpspsB（7.4.4）式中B强化参数，定义为：/1ttppEddBEdEdd（7.4.5）式中tE切线模量。当ps时，还有vppvt（7.4.6）图 7-l0 三元件粘弹塑性模型结合式 7.4.4 和式 7.4.6，得vpvpsBt（7.4.7）结合式 7.4.1、式 7.4.2 和式 7.4.7，得sBEEBEtt（7.4.8）记1，称为介质流动参数，则式 7.4.8 可改写为
21、：vpsBE（7.4.9）因此，粘塑性应变率为：vpvpsB（7.4.10）式 7.4.10 表明粘塑性应变率是由超过稳态屈服应力的那部分应力值（称为“过应力”）所决定的。若作用于模型的应力为常值A时，即/0t，式 7.4.8 可改写为：AAsBBtE （7.4.11）式 7.4.11 的解为：1 expAsAB tEB（7.4.12）对于理想粘塑性材料，0B，利用罗比达法则，式 7.4.12 可改写为-第 14 页AAstE（7.4.13）对于更复杂的弹粘塑性模型读者可参阅有关专著，这里不作进步介绍了。7.5 蠕变物体的蠕变现象可以采用由一定数量的弹性、粘性和塑性元件组成的模型来描述。但元件
22、多了，计算相当复杂，且其关系不容易由试验确定。在实际应用中，常常直接由试验来确定应力、应变和时间之间的关系。下面简单介绍几种主要的蠕变方程的形式。1 老化理论老化理论假设蠕变应变与应力、时间之间具有某种函数关系，即,cft（7.5.1）式中c蠕变应变。物体的总应变将由三部分组成，即epc（7.5.2）式中,cp分别为弹性应变和塑性应变。当应力未超过屈服应力时，0p，则式 7.5.2 可改写为：,ftE（7.5.3）图 7-11 应力-应变-时间实测曲线Buisman（1936）根据大量的试验资料，认为饱和粘土的时间与沉降关系在半对数坐标上呈线性关系（图 7-11)，即logapcht（7.5.
23、4）式中沉降；ah土层或试样起始高度；竖向应力增量；p主固结系数；c时间效应系数。式 7.5.4 可改写成：logtat（7.5.5）-第 15 页式中/tah；2 流动理论流动理论认为蠕变应变速率与应力、时间之间存在某种简单的函数关系，即,cft（7.5.6）式中蠕变应变速率。物体的总应变速率可表示为：epc（7.5.7）式中,ep分别表示弹性应变速率和塑性应变速率。试验资料表明，正常固结粘土和超固结粘土，在排水或不排水条件下，其应变速率和蠕变时间的关系在全对数坐标上呈线性关系。应变速率与应力的关系也是线性关系（图 7-12）。应变速率与时问的关系为：11lnln,lntt qmt（7.5.
24、8）式中应变速率；1,t q单位时间的应变速率，为应力 q 的函数；1t单位时间（即 1 分钟）；m关系曲线的坡度，一般为 0.751。应变速率与应力的关系可表示为：10,lnlnt qq（7.5.9）式中10,t q时应变速率，为蠕变时间 t 的函数；q偏应力，13q；关系曲线的坡度。式 7.5.9 也可改写为：1expmtAqt（7.5.10）式中A1tt时的曲线延长至0q 时，在坐标上的截距。-第 16 页式 7.5.10 表示应变速率，时间 t 和应力 q 三者之间的关系，反映土的蠕变特性。积分式 7.5.10，得 1exp1mAtqm1m（7.5.11）图 7-12 应变速率与时间、
25、应力的关系（引自 Singh and Mitchell,1968)1ln/expAq1,1mt（7.5.12）式中1exp1Aqm；11t 时的应变值。3 遗传理论1874 年 Boltzmann 首先提出了线性遗传理论的本构力程，即 0ttK tdE （7.5.13）式中时间 t 时的总应变；t时间 t 时的应力；K t在瞬时作用的应力对时间的变形的影响函数（遗传函数），当t增加时，函数值单调减小。为了进一步了解遗传函数的意义，现考虑在0t 时施加常应力的情况。当常数时，式 7.5.13 可改写为：01tK tdE（7.5.14）进行变量置换，令t，则式 7.5.14 可改写为：01tKdE（7.5.15）上式对时间 t 求导，得蠕变速率为：K t（7.5.16）于是 K t（7.5.17）上式表明函数 K t是试件在常应力1作用下的应变速率。由此可得到在常荷载作用下用试验确定遗传函数的简单方法。-第 17 页思考题与习题1.简述理想弹性模型、理想塑性模型和理想粘性模型的主要内并说明其物理意义。2.什么是 Maxwell 模型、Kelvin 模型和三元件粘弹性模型？并能用简图表示。3.试简要介绍 Bingham 模型。4.举例说明粘弹塑性模型的主要特点。5.试说明粘性流动与塑胜流动有什么区别？