清华大学计算固体力学第五次本构模型学习教案

1、会计学1清华大学计算固体清华大学计算固体(gt)力学第五次本构力学第五次本构模型模型第一页，共94页。第第5 5章章 本构模型本构模型(mxng) (mxng) 1 1引言引言(ynyn)(ynyn)2 2应力应力- -应变曲线应变曲线 3 3一维弹性一维弹性 4 4非线性弹性（超弹性）非线性弹性（超弹性）5 5一维塑性一维塑性 6 6多轴塑性多轴塑性7 7超弹塑性模型超弹塑性模型8 8粘弹性粘弹性 9 9应力更新算法应力更新算法1010连续介质力学与本构模型连续介质力学与本构模型第1页/共94页第二页，共94页。1 1 引言引言(ynyn) (ynyn) 本构方程率形式的积分算法称为应力更新

2、算法（也称为本构更新算法）本构方程率形式的积分算法称为应力更新算法（也称为本构更新算法），包括：，包括：径向返回算法的一类图形返回算法，径向返回算法的一类图形返回算法，算法模量与基本应力更新方案一致的概念，算法模量与基本应力更新方案一致的概念，大变形问题大变形问题(wnt)的增量客观应力更新方案，的增量客观应力更新方案，基于弹性响应的应力更新方案，自动满足客观性的超弹性势能。基于弹性响应的应力更新方案，自动满足客观性的超弹性势能。 为了进行分析，选择材料模型是很重要，往往又不是很明确为了进行分析，选择材料模型是很重要，往往又不是很明确，仅有的信息可能是一般性的知识和经验，即可能是材料行为的，仅

3、有的信息可能是一般性的知识和经验，即可能是材料行为的几条应力应变曲线。几条应力应变曲线。 在有限元软件库中选择合适的本构模型，如果在有限元软件库中选择合适的本构模型，如果(rgu)(rgu)没有合没有合适的本构模型，要开发用户材料子程序。重要的是理解本构模型适的本构模型，要开发用户材料子程序。重要的是理解本构模型的关键特征，创建模型的假设，材料、荷载和变形域、以及程序的关键特征，创建模型的假设，材料、荷载和变形域、以及程序中的数值问题是否适合模型。中的数值问题是否适合模型。第2页/共94页第三页，共94页。2 应力应力(yngl)-应应变曲线变曲线 材料应力应变材料应力应变(yngbin)(y

4、ngbin)行为的许多基本特征可以从一维应力行为的许多基本特征可以从一维应力状态状态( (单轴应力或者剪切单轴应力或者剪切) )的一组应力的一组应力- -应变应变(yngbin)(yngbin)曲线中获得，曲线中获得，多轴状态的本构方程常常基于在试验中观察到的一维行为而简单生成多轴状态的本构方程常常基于在试验中观察到的一维行为而简单生成。 载荷位移载荷位移(wiy)曲线曲线 名义应力（工程应力）给出为名义应力（工程应力）给出为 定义伸长定义伸长 工程应变定义为工程应变定义为 第3页/共94页第四页，共94页。2 2 应力应力(yngl)-(yngl)-应应变曲线变曲线 CauchyCauchy

5、（或者真实（或者真实(zhnsh)(zhnsh)）应力表示为）应力表示为 以每单位当前长度以每单位当前长度(chngd)(chngd)应变的增量随长度应变的增量随长度(chngd)(chngd)的变化得到另一种应变度量的变化得到另一种应变度量 对数应变（也称为真实应变）对数应变（也称为真实应变） 对材料时间求导，表达式为对材料时间求导，表达式为一维情况，上式为变形率一维情况，上式为变形率 当前面积的表达式给出为当前面积的表达式给出为真实应力应变曲线 工程应力应变曲线第4页/共94页第五页，共94页。2 2 应力应力(yngl)-(yngl)-应变曲线应变曲线 考虑一种不可压缩材料考虑一种不可压

6、缩材料(J(J1)1)，名义应力和工程，名义应力和工程(gngchng)(gngchng)应应变的关系为变的关系为真实真实(zhnsh)(zhnsh)应力应力( (对于不可压缩材对于不可压缩材料料) )1000 xxLLLL说明了对于本构行为应用不同泛函表达式的区别，对于同样材料取决于说明了对于本构行为应用不同泛函表达式的区别，对于同样材料取决于采用何种应力和变形的度量。采用何种应力和变形的度量。 应力应变曲线的显著特征之一是非线性的度。材料线弹性行为应力应变曲线的显著特征之一是非线性的度。材料线弹性行为的范围小于应变的百分之几，就可以采用小应变理论描述。的范围小于应变的百分之几，就可以采用小

7、应变理论描述。第5页/共94页第六页，共94页。2 2 应力应力- -应变应变(yngbin)(yngbin)曲线曲线 应力应变反应与变形率无关的材料称为率无关；否则，称为率相关(xinggun)。名义应变率定义为 率无关和率相关(xinggun)材料的一维反应因为 和即名义应变率等于伸长率，例如即名义应变率等于伸长率，例如 可以看出，对于率可以看出，对于率无关材料的应力应变无关材料的应力应变曲线是应变率独立的，曲线是应变率独立的，而对于率相关材料的应而对于率相关材料的应力应变曲线，力应变曲线，当应变当应变率提高时是上升的率提高时是上升的；而；而当温度升高时是下降当温度升高时是下降的。的。第6

8、页/共94页第七页，共94页。2 2 应力应力- -应变应变(yngbin)(yngbin)曲线曲线 对于弹性材料，应力应变的对于弹性材料，应力应变的卸载曲线简单地沿加载曲线返回，卸载曲线简单地沿加载曲线返回，直到完全卸载，材料返回到了它的直到完全卸载，材料返回到了它的初始未伸长状态。然而，对于弹初始未伸长状态。然而，对于弹塑性材料，卸载曲线区别于加载曲塑性材料，卸载曲线区别于加载曲线，卸载曲线的斜率是典型的应力线，卸载曲线的斜率是典型的应力应变弹性（初始）段的斜率，卸应变弹性（初始）段的斜率，卸载后产生永久应变。其它材料的行载后产生永久应变。其它材料的行为介于这两种极端之间。由于在加为介于这

9、两种极端之间。由于在加载过程中微裂纹的形成材料已经损载过程中微裂纹的形成材料已经损伤，脆性材料的卸载行为，当荷载伤，脆性材料的卸载行为，当荷载移去后微裂纹闭合，弹性应变得到移去后微裂纹闭合，弹性应变得到(d do)(d do)恢复。卸载曲线的初始斜恢复。卸载曲线的初始斜率给出形成微裂纹损伤程度的信息率给出形成微裂纹损伤程度的信息。(a)弹性(tnxng),(b)弹-塑性,(c)弹性(tnxng)含损伤 第7页/共94页第八页，共94页。3 3 一维弹性一维弹性(tnxng) (tnxng) 弹性材料的基本性能是应力仅依赖于应变的当前水平。这意味着弹性材料的基本性能是应力仅依赖于应变的当前水平。

10、这意味着加载和卸载的应力加载和卸载的应力- -应变曲线应变曲线(qxin)(qxin)是一致的，当卸载结束时材料恢是一致的，当卸载结束时材料恢复到初始状态。称这种应变是可逆的。而且，弹性材料是率无关的复到初始状态。称这种应变是可逆的。而且，弹性材料是率无关的( (与与应变率无关应变率无关) )。弹性材料的应力和应变是一一对应的。弹性材料的应力和应变是一一对应的。小应变小应变(yngbin) 可逆和路径无关默认在变形中没有能量耗散可逆和路径无关默认在变形中没有能量耗散，在弹性材料中，储存，在弹性材料中，储存在物体中的能量全部消耗在变形中，卸载后材料恢复。在物体中的能量全部消耗在变形中，卸载后材料

11、恢复。 对于一维弹性材料，对于一维弹性材料，可逆、路径无关、无能量耗散可逆、路径无关、无能量耗散是等价的特征。是等价的特征。对于二维和三维弹性，以及超弹性材料，也类似。对于二维和三维弹性，以及超弹性材料，也类似。对于任意应变，不管如何达到应变值，上式给出唯一应力值。对于任意应变，不管如何达到应变值，上式给出唯一应力值。 第8页/共94页第九页，共94页。3 3 一维弹性一维弹性(tnxng) (tnxng) 应 变应 变 ( y n g b i n )( y n g b i n ) 能 一 般 是 应 变能 一 般 是 应 变(yngbin)(yngbin)的凸函数，例如，的凸函数，例如， (

12、a)凸应变(yngbin)能函数 (b)应力应变(yngbin)曲线 当当 公式的等号成立。公式的等号成立。 凸应变能函数的一个例子如图所示。在这种情况下，函数是凸应变能函数的一个例子如图所示。在这种情况下，函数是单调递增的，如果单调递增的，如果w 是非凸函数，则是非凸函数，则 s 先增后减，材料应变软化先增后减，材料应变软化，这是非稳定的材料反应，这是非稳定的材料反应， 如右下图。如右下图。(a)非凸应变能函数(b)相应的应力应变曲线 第9页/共94页第十页，共94页。大应变大应变(yngbin) 从弹性推广到大应变，只要选择应变度量和定义应力（功共从弹性推广到大应变，只要选择应变度量和定义

13、应力（功共轭）的弹性势能。势能的存在是默认了可逆、路径无关轭）的弹性势能。势能的存在是默认了可逆、路径无关(wgun)(wgun)和无能量耗散。如和无能量耗散。如 3 3 一维弹性一维弹性(tnxng) (tnxng) 在弹性应力在弹性应力- -应变关系中，从应变的势函数可以获得应力为超弹性。如应变关系中，从应变的势函数可以获得应力为超弹性。如一维大应变问题，以一维大应变问题，以GreenGreen应变的二次函数表示应变的二次函数表示 对于小应变问题，即为胡克定律。对于小应变问题，即为胡克定律。第10页/共94页第十一页，共94页。大应变大应变(yngbin) 一种材料的一种材料的Cauchy

14、Cauchy应力率与变形率相关，称为次弹性应力率与变形率相关，称为次弹性(tnxng)(tnxng)。这种关系一般是非线性的，给出为这种关系一般是非线性的，给出为 3 3 一维弹性一维弹性(tnxng) (tnxng) 一个特殊的线性次弹性关系给出为一个特殊的线性次弹性关系给出为 这是与路径无关的超弹性关系。对于多轴问题，一般次弹性这是与路径无关的超弹性关系。对于多轴问题，一般次弹性关系不能转换到超弹性，它仅在一维情况下是严格路径无关的。关系不能转换到超弹性，它仅在一维情况下是严格路径无关的。然而，如果是弹性小应变，其行为足以接近路径无关的弹性行为然而，如果是弹性小应变，其行为足以接近路径无关

15、的弹性行为。因为次弹性的简单性，公式。因为次弹性的简单性，公式(5.3.11)(5.3.11)的多轴一般形式常常应用在的多轴一般形式常常应用在有限元软件中，以模拟大应变弹塑性的弹性反应。有限元软件中，以模拟大应变弹塑性的弹性反应。对上式的关系积分，得到对上式的关系积分，得到第11页/共94页第十二页，共94页。4 4 非线性弹性非线性弹性(tnxng) (tnxng) 对于有限应变有许多不同的应力和变形对于有限应变有许多不同的应力和变形(bin xng)(bin xng)度量，同样的度量，同样的本构关系可以写成几种不同的形式，总是可能从一种形式的本构关系转本构关系可以写成几种不同的形式，总是可

16、能从一种形式的本构关系转换到另一种形式。换到另一种形式。 大应变弹性本构模型首先表述成大应变弹性本构模型首先表述成KirchhoffKirchhoff材料的一种特殊形式，材料的一种特殊形式，由线弹性直接生成到大变形由线弹性直接生成到大变形(bin xng)(bin xng)。满足路径无关、可逆和无能。满足路径无关、可逆和无能量耗散。因此，路径无关的程度可以视为材料模型弹性的度量。量耗散。因此，路径无关的程度可以视为材料模型弹性的度量。 次弹性材料是路径无关程度最弱的材料，遵从次弹性材料是路径无关程度最弱的材料，遵从CauchyCauchy弹性，其应力弹性，其应力是路径无关的，但是其能量不是路径

17、无关的。是路径无关的，但是其能量不是路径无关的。 超弹性材料或者超弹性材料或者GreenGreen弹性，它是路径无关和完全可逆的，应力由应弹性，它是路径无关和完全可逆的，应力由应变势能导出。变势能导出。第12页/共94页第十三页，共94页。4 4 非线性弹性非线性弹性(tnxng) (tnxng) 小应变小应变(yngbin)和大转动和大转动 式中式中 C C 为弹性模量为弹性模量( (切线模量切线模量) )的四阶张量，对的四阶张量，对KirchhoffKirchhoff材料是常数，代表了应材料是常数，代表了应力和应变的多轴状态。它可以完全力和应变的多轴状态。它可以完全(wnqun)(wnqu

18、n)反映材料的各向异性。反映材料的各向异性。 许多工程应用包括许多工程应用包括小应变和大转动小应变和大转动。在这些问题中，大变形的效果。在这些问题中，大变形的效果主要来自于大转动，如直升机旋翼、船上升降器或者钓鱼杆的弯主要来自于大转动，如直升机旋翼、船上升降器或者钓鱼杆的弯曲。由线弹性定律的简单扩展即可以模拟材料的反应，但要以曲。由线弹性定律的简单扩展即可以模拟材料的反应，但要以PK2PK2应力代替其中的应力和以应力代替其中的应力和以GreenGreen应变代替线性应变，这称为应变代替线性应变，这称为Saint-Venant- Saint-Venant- KirchhoffKirchhoff材

19、料，或者简称为材料，或者简称为KirchhoffKirchhoff材料。最一般的材料。最一般的KirchhoffKirchhoff模型为模型为第13页/共94页第十四页，共94页。4 4 非线性弹性非线性弹性(tnxng) (tnxng) E:CS klijklijECS式中式中C C为弹性模量的四阶张量，有为弹性模量的四阶张量，有8181个常数个常数(chngsh)(chngsh)。利用对。利用对称性可以显著地减少常数称性可以显著地减少常数(chngsh)(chngsh)。 一般的四阶张量有一般的四阶张量有3481个独立常数个独立常数(chngsh)，与全应力，与全应力张量的张量的9个分量和

20、全应变张量的个分量和全应变张量的9个分量有关。个分量有关。 如次弹性本构方程如次弹性本构方程第14页/共94页第十五页，共94页。这样这样C C为对称矩阵（主对称性），在为对称矩阵（主对称性），在8181个常数个常数(chngsh)(chngsh)中有中有4545个是独立的个是独立的。成为上三角或下三角矩阵。成为上三角或下三角矩阵。 4 4 非线性弹性非线性弹性(tnxng) (tnxng) 利用势能利用势能(shnng)(shnng)表示的应力应变关系和表示的应力应变关系和GreenGreen公式，公式， 故有故有 应力张量和应变张量均为对称张量（应力张量和应变张量均为对称张量（次对称性次对

21、称性），即），即 第15页/共94页第十六页，共94页。4 4 非线性弹性非线性弹性(tnxng) (tnxng) 应力应力(yngl)(yngl)张量和应变张量均为对称张量（次对称性），即张量和应变张量均为对称张量（次对称性），即 jiijjiijEESS,再利用再利用(lyng)(lyng)模量的主对称性使独立弹性常数的数目减少，由模量的主对称性使独立弹性常数的数目减少，由3636个常数减少为个常数减少为2121个，为各向异性材料。个，为各向异性材料。 应力和应变张量的对称性要求应力的应力和应变张量的对称性要求应力的6 6个独立分量仅与应变的个独立分量仅与应变的6 6个独立个独立分量有关，

22、由弹性模量的局部对称结果，独立常数的数目减少到分量有关，由弹性模量的局部对称结果，独立常数的数目减少到3636个。个。 第16页/共94页第十七页，共94页。4 4 非线性弹性非线性弹性(tnxng) (tnxng) 写成矩阵写成矩阵(j zhn)形式为（可以是上或下三角矩阵形式为（可以是上或下三角矩阵(j zhn)） 对于对于(duy)(duy)正交各向异性，具有正交的三个弹性对称面，当坐标正交各向异性，具有正交的三个弹性对称面，当坐标变号，为使应变能密度不变，有变号，为使应变能密度不变，有 这样由这样由2121个常数减少为个常数减少为1414个，为个，为正交各向异性材料正交各向异性材料。

23、若材料对称坐标平面，当沿轴平面反射时，弹性模量不变，固为若材料对称坐标平面，当沿轴平面反射时，弹性模量不变，固为正正交各向异性体交各向异性体，有，有 第17页/共94页第十八页，共94页。 对于一个由三个彼此正交的对称平面组成的正交材料对于一个由三个彼此正交的对称平面组成的正交材料( (如木材如木材或纤维或纤维(xinwi)(xinwi)增强的复合材料增强的复合材料) )，仅有，仅有9 9个独立弹性常数，个独立弹性常数，KirchhoffKirchhoff应力应变关系为材料对称坐标平面，为正交各向异性应力应变关系为材料对称坐标平面，为正交各向异性体体4 4 非线性弹性非线性弹性(tnxng)

24、(tnxng) 对于各向同性对于各向同性( xin tn xn)( xin tn xn)材料，仅有材料，仅有3 3个常数个常数 第18页/共94页第十九页，共94页。4 4 非线性弹性非线性弹性(tnxng) (tnxng) 小应变小应变(yngbin)和和大转动大转动 对于各向同性的对于各向同性的Kirchhoff材料，其应力应变材料，其应力应变(yngbin)关系可以写成为关系可以写成为式中式中Lam常数，体积模量常数，体积模量K，杨氏模量，杨氏模量 E和泊松比和泊松比的关系为的关系为 材料对称的一个重要的例子是各向同性。一个各向同性材料没有方材料对称的一个重要的例子是各向同性。一个各向同

25、性材料没有方位或者方向的选择，因此，当以任何直角坐标系表示的应力应变关系位或者方向的选择，因此，当以任何直角坐标系表示的应力应变关系是等同的。对于小应变的许多材料是等同的。对于小应变的许多材料( (如金属和陶瓷如金属和陶瓷) )可以作为各向同性进可以作为各向同性进行模拟。张量行模拟。张量C C是各向同性的。在任何坐标系统中，一个各向同性张量有是各向同性的。在任何坐标系统中，一个各向同性张量有相同的分量。相同的分量。(克罗内克)符号构成的一个线性组合： 第19页/共94页第二十页，共94页。4 4 非线性弹性非线性弹性(tnxng) (tnxng) 不可压缩性不可压缩性 在变形的过程中，不可压缩

26、材料在变形的过程中，不可压缩材料(cilio)的体积不变，密度保持常数。的体积不变，密度保持常数。不可压缩材料不可压缩材料(cilio)的运动称为等体积运动。的运动称为等体积运动。 总体总体(zngt)变形变形 等体积约束运动的率形式等体积约束运动的率形式 将应力和应变率度量写成偏量和静水（体积的）部分的和，将应力和应变率度量写成偏量和静水（体积的）部分的和，对于不可压缩材料，静水部分也称为张量的球形部分，分解式为：对于不可压缩材料，静水部分也称为张量的球形部分，分解式为： 对于不可压缩材料，压力不能从本构方程确定，而是从动量方程确定。对于不可压缩材料，压力不能从本构方程确定，而是从动量方程确

27、定。 第20页/共94页第二十一页，共94页。4 非线性弹性非线性弹性(tnxng) Kirchhoff应力应力(yngl) 由由Jacobian行列式放大行列式放大(fngd)，称它为权重，称它为权重Cauchy应力。对于等体积应力。对于等体积运动，它等同于运动，它等同于Cauchy应力。应力。 次弹性次弹性次弹性材料规律联系应力率和变形率。次弹性材料规律联系应力率和变形率。 上式是率无关、线性增加和可逆的。对于有限变形状态的微小增量，应力和应变的增量是线性关系，当卸载后可以恢复。然而，对于大变形能量不一定必须守恒，并且在闭合变形轨迹上作的功不一定必须为零。次弹性规律主要用来代表在弹-塑性规

28、律中的弹性反应，小变形弹性，且耗能效果也小。 第21页/共94页第二十二页，共94页。4 非线性弹性非线性弹性(tnxng) 切线切线(qixin)模量之间的关系模量之间的关系 对于各向同性对于各向同性( xin tn xn)材料材料Jaumann率的切线模量为率的切线模量为 某些次弹性本构关系共同应用的形式为某些次弹性本构关系共同应用的形式为对于同一种材料，切线模量不同，材料反应的率形式不同，如对于同一种材料，切线模量不同，材料反应的率形式不同，如 如果如果是常数，是常数，不是常数。不是常数。 切线模量切线模量证明见第证明见第5.4.5节，推导复杂节，推导复杂第22页/共94页第二十三页，共

29、94页。4 非线性弹性非线性弹性(tnxng) 超弹性超弹性(tnxng)材料材料 平衡方程是以物体中应力的形式建立的，应力来源于变形，如应变。如果本构行为仅是变形的当前状态的函数，为与时间无关的弹性本构。而对于接近(jijn)不可压缩的材料，仅依赖变形（应变）不一定能够得到应力。 储存在材料中的能量（功）仅取决于变形的初始和最终状态，并且是独立于变形（或荷载）路径，称这种弹性材料为超弹性（hyper-elastic）材料，或者为Green弹性，例如常用的工业橡胶。动物的肌肉也具有超弹性的力学性质。这里主要讨论橡胶材料的超弹性力学行为。第23页/共94页第二十四页，共94页。4 非线性弹性非线

30、性弹性(tnxng) 超弹性超弹性(tnxng)材料材料 对于对于(duy)功独立于荷载路径的弹性材料称之为超弹性（功独立于荷载路径的弹性材料称之为超弹性（Green弹性）材料。超弹性材料的特征是存在一个潜在（或应变弹性）材料。超弹性材料的特征是存在一个潜在（或应变）能量函数，它是应力的势能：）能量函数，它是应力的势能： 通过适当转换获得了对于不同应力度量的表达式通过适当转换获得了对于不同应力度量的表达式 由于变形梯度张量由于变形梯度张量F是不对称的，因此名义应力张量是不对称的，因此名义应力张量P的的9个分量是个分量是不对称的。不对称的。 在橡胶大变形中应用在橡胶大变形中应用多项式模型多项式模

31、型和和Ogden指数模型指数模型。第24页/共94页第二十五页，共94页。4 非线性弹性非线性弹性(tnxng) 超弹性超弹性(tnxng)材料材料 目前，世界半数以上的橡胶是合成橡胶。合成橡胶的种类很多，目前，世界半数以上的橡胶是合成橡胶。合成橡胶的种类很多，例如，制造轮胎使用的丁苯橡胶（苯乙烯和丁二烯的共聚物）或乙例如，制造轮胎使用的丁苯橡胶（苯乙烯和丁二烯的共聚物）或乙丙烯橡胶（丙烯橡胶（ERP）；用于汽车配件的有氯丁橡胶及另一种具有天然）；用于汽车配件的有氯丁橡胶及另一种具有天然橡胶各种性能的异戊橡胶。橡胶各种性能的异戊橡胶。 在众多的合成橡胶中，硅橡胶是其中的佼佼者。它具有无味在众多

32、的合成橡胶中，硅橡胶是其中的佼佼者。它具有无味(wwi)无毒，不怕高温和严寒的特点，在摄氏无毒，不怕高温和严寒的特点，在摄氏300度和零下度和零下90度时度时能够能够“泰然自若泰然自若”、“面不改色面不改色”，仍不失原有的强度和弹性。例，仍不失原有的强度和弹性。例如生物材料。如生物材料。 橡胶是提取橡胶树、橡胶草等植物的胶乳，加工后制成的具有弹性、绝缘橡胶是提取橡胶树、橡胶草等植物的胶乳，加工后制成的具有弹性、绝缘性、不透水和空气的材料。在半个世纪前，性、不透水和空气的材料。在半个世纪前，“橡胶橡胶”一词是专指生橡胶，它一词是专指生橡胶，它是从热带植物巴西三叶胶的胶乳提炼出来的。是从热带植物巴

33、西三叶胶的胶乳提炼出来的。第25页/共94页第二十六页，共94页。4 非线性弹性非线性弹性(tnxng) 超弹性超弹性(tnxng)材料材料 1839年，年，Charle Goodyear发明了橡胶的硫化方法，其姓发明了橡胶的硫化方法，其姓氏现在已经成为国际上著名橡胶轮胎的商标。氏现在已经成为国际上著名橡胶轮胎的商标。 从从19世纪中叶起橡胶就成为一种世纪中叶起橡胶就成为一种(y zhn)重要的工程材重要的工程材料。然而，橡胶材料的行为复杂，不同于金属材料仅需要几料。然而，橡胶材料的行为复杂，不同于金属材料仅需要几个参数就可以描述材料特性。橡胶材料受力以后，变形是伴个参数就可以描述材料特性。橡

34、胶材料受力以后，变形是伴随着大位移和大应变，其本构关系是非线性的，并且在变形随着大位移和大应变，其本构关系是非线性的，并且在变形过程中体积几乎保持不变。过程中体积几乎保持不变。 橡胶具有许多特殊的性能，例如电绝缘性、耐氧老化性、耐橡胶具有许多特殊的性能，例如电绝缘性、耐氧老化性、耐光老化性、防霉性、化学稳定性等。光老化性、防霉性、化学稳定性等。第26页/共94页第二十七页，共94页。4 非线性弹性非线性弹性(tnxng) 超弹性超弹性(tnxng)材料材料 由于计算机以及有限元数值分析的飞速发展，我们可以借助计算机来由于计算机以及有限元数值分析的飞速发展，我们可以借助计算机来对超弹性材料的工程

35、应用进行深入研究以及优化设计对超弹性材料的工程应用进行深入研究以及优化设计(shj)。可以用有。可以用有限元等数值方法来计算分析橡胶元件的力学性能，包括选取和拟合橡胶限元等数值方法来计算分析橡胶元件的力学性能，包括选取和拟合橡胶的本构模型，以及用有限元建模和处理计算结果等。的本构模型，以及用有限元建模和处理计算结果等。 橡胶是一种弹性聚合物，其特点是有很强的非线性粘弹性行为。它橡胶是一种弹性聚合物，其特点是有很强的非线性粘弹性行为。它的力学行为对温度、环境、应变历史、加载速率都非常敏感，这样使的力学行为对温度、环境、应变历史、加载速率都非常敏感，这样使得描述橡胶的行为变得非常复杂。橡胶的制造工

36、艺和成分也对橡胶的得描述橡胶的行为变得非常复杂。橡胶的制造工艺和成分也对橡胶的力学性能有着显著的影响。力学性能有着显著的影响。第27页/共94页第二十八页，共94页。固体橡胶材料的拉伸试验固体橡胶材料的拉伸试验(shyn)曲线与材料演化曲线与材料演化模型模型 固体橡胶是几乎不可压缩固体橡胶是几乎不可压缩(y su)的，其泊松比接近于的，其泊松比接近于0.5。可逆，大应。可逆，大应变。初始各向同性，应变增加后分子定向排列。变。初始各向同性，应变增加后分子定向排列。4 非线性弹性非线性弹性(tnxng) 超弹性材料超弹性材料 常用的橡胶性态可分为常用的橡胶性态可分为固体橡胶固体橡胶和和泡沫橡胶泡沫

37、橡胶。第28页/共94页第二十九页，共94页。4 非线性弹性非线性弹性(tnxng) 超弹性超弹性(tnxng)材料材料 一般将多孔橡胶或弹性泡沫材料统称为泡沫材料。弹性泡一般将多孔橡胶或弹性泡沫材料统称为泡沫材料。弹性泡沫材料的普通例子有多孔聚合物，如海绵、包装材料等。沫材料的普通例子有多孔聚合物，如海绵、包装材料等。 泡沫橡胶是由橡胶制成的弹性泡沫材料，能够泡沫橡胶是由橡胶制成的弹性泡沫材料，能够(nnggu)满足非常大的弹性应变要求，拉伸时的应变可以达到满足非常大的弹性应变要求，拉伸时的应变可以达到500或或更大，压缩时的应变可以达到更大，压缩时的应变可以达到90或更小。与固体橡胶的几乎

38、或更小。与固体橡胶的几乎不可压缩性相比，泡沫材料的多孔性则允许非常大的体积缩小不可压缩性相比，泡沫材料的多孔性则允许非常大的体积缩小变形，因此具有良好的能量吸收性。变形，因此具有良好的能量吸收性。泡沫橡胶材料的多面体微元模型泡沫橡胶材料的多面体微元模型 a) 开放腔室，开放腔室，b) 封闭腔室封闭腔室第29页/共94页第三十页，共94页。4 非线性弹性非线性弹性(tnxng) 超弹性超弹性(tnxng)材料材料 泡沫橡胶材料的应力泡沫橡胶材料的应力-应变应变(yngbin)曲线曲线 a)压缩压缩 b)拉伸拉伸小应变小应变 5%，线弹性，泊松比为，线弹性，泊松比为0.3 。大应变，压缩时，泊松比

39、为大应变，压缩时，泊松比为0.0； 拉伸时，泊松比大于拉伸时，泊松比大于0.0。第30页/共94页第三十一页，共94页。100%典型固体橡胶材料单轴拉伸典型固体橡胶材料单轴拉伸(l shn)应力应力-应变曲线应变曲线 橡胶橡胶(xingjio)本构模型本构模型 4 非线性弹性非线性弹性(tnxng) 第31页/共94页第三十二页，共94页。小变形小变形(bin xng) 以多项式形式以多项式形式(xngsh)本构模型为例，其应变能密度表达式为本构模型为例，其应变能密度表达式为忽略忽略(hl)二阶及二阶以上小量，变二阶及二阶以上小量，变为为弹性常数为弹性常数为 当当 橡胶本构模型橡胶本构模型 4

40、 非线性弹性非线性弹性 第32页/共94页第三十三页，共94页。0LLx定义定义(dngy)伸长伸长 工程工程(gngchng)应变定义为应变定义为 二阶张量基本二阶张量基本(jbn)不变量不变量 小变形，有小变形，有 小变形小变形 橡胶本构模型橡胶本构模型 4 非线性弹性非线性弹性 第33页/共94页第三十四页，共94页。例题例题(lt) 在超弹性计算中，橡胶使用三次减缩在超弹性计算中，橡胶使用三次减缩(jin su)多项式应变能本构模型，多项式应变能本构模型，应变能密度表达式为应变能密度表达式为若取若取（单位为（单位为MPa），求材料弹性），求材料弹性(tnxng)常数。常数。 利用公式利

41、用公式解：解：101102),(2DKCCG解出橡胶的弹性常数为解出橡胶的弹性常数为 ， E=1.384MPa，= 0.5 小变形小变形 橡胶本构模型橡胶本构模型 4 非线性弹性非线性弹性 第34页/共94页第三十五页，共94页。 常用的橡胶力学性能描述方法主要分为两类，一类是基于热常用的橡胶力学性能描述方法主要分为两类，一类是基于热力学统计力学统计(tngj)的方法，另一类是基于橡胶为连续介质的唯象学的方法，另一类是基于橡胶为连续介质的唯象学描述方法。描述方法。 热力学统计热力学统计(tngj)方法的基础为观察到橡胶中的弹性恢复力方法的基础为观察到橡胶中的弹性恢复力主要来自熵的减少。橡胶在承

42、受荷载时分子结构无序，熵的减少是主要来自熵的减少。橡胶在承受荷载时分子结构无序，熵的减少是由于橡胶伸长使得橡胶结构由高度无序变得有序。由对橡胶中分子由于橡胶伸长使得橡胶结构由高度无序变得有序。由对橡胶中分子链的长度、方向以及结构的统计链的长度、方向以及结构的统计(tngj)得到本构关系。得到本构关系。橡胶橡胶(xingjio)本构模型本构模型 唯象学描述方法假设在未变形状态下橡胶为各向同性材料，即长分子链唯象学描述方法假设在未变形状态下橡胶为各向同性材料，即长分子链方向在橡胶中是随机分布的。这种各向同性的假设是用单位体积（弹性）应方向在橡胶中是随机分布的。这种各向同性的假设是用单位体积（弹性）

43、应变能函数（变能函数（U）来描述橡胶特性）来描述橡胶特性(txng)的基础，其本构模型为多项式形式模的基础，其本构模型为多项式形式模型和型和Ogden形式模型。形式模型。第35页/共94页第三十六页，共94页。NiiiNjijiijJDIICU12121) 1(1) 3() 3(典型典型(dinxng)的本构模型为多项式形式，其应变能密度表达式为的本构模型为多项式形式，其应变能密度表达式为特殊形式可以由设定某些特殊形式可以由设定某些(mu xi)参数为参数为0来得到。如果所有来得到。如果所有 则得到则得到(d do)减缩多项式模型减缩多项式模型 对于完全多项式对于完全多项式，如果，如果, 则只

44、有线性部分的应变能量，则只有线性部分的应变能量，即即Mooney-Rivlin形式形式橡胶本构模型橡胶本构模型 第36页/共94页第三十七页，共94页。，则得到，则得到(d do)Neo-Hookean形式形式 对于减缩对于减缩(jin su)多项式，如果多项式，如果 N1100%100%Mooney-Rivlin形式和形式和Neo-Hooken形式本构模型形式本构模型(mxng)（后者是将（后者是将Hooke定律扩展至大变形）定律扩展至大变形）橡胶本构模型橡胶本构模型 第37页/共94页第三十八页，共94页。Yeoh形式形式(xngsh)本构模型是本构模型是 时减缩时减缩(jin su)多项

45、式的特殊形式多项式的特殊形式 100%典型的典型的S形橡胶形橡胶(xingjio)应力应力-应变曲线应变曲线 ，C10正值，在小变形时为切线模量；正值，在小变形时为切线模量；C20为负值，中等变形时软化；为负值，中等变形时软化；C30正值，大变形时硬化。正值，大变形时硬化。橡胶本构模型橡胶本构模型 第38页/共94页第三十九页，共94页。Ogden形式形式(xngsh)本构模型本构模型 Arruda-Boyce形式形式(xngsh)本构模型本构模型 Van der Waals模型模型(mxng) 橡胶本构模型橡胶本构模型 其他形式的本构模型有：其他形式的本构模型有：第39页/共94页第四十页，

46、共94页。试验试验(shyn)拟合本构模型系数拟合本构模型系数 橡胶类材料的本构关系除具有超弹性、大变形的特征外，其本构关系与橡胶类材料的本构关系除具有超弹性、大变形的特征外，其本构关系与生产生产(shngchn)加工过程有直接关系，如橡胶配方和硫化工艺。确定每一加工过程有直接关系，如橡胶配方和硫化工艺。确定每一批新加工出来的橡胶的本构关系，都要依赖于精确和充分的橡胶试验。批新加工出来的橡胶的本构关系，都要依赖于精确和充分的橡胶试验。第40页/共94页第四十一页，共94页。 通常在试验通常在试验(shyn)中应该测得在几种不同荷载模式下的应力中应该测得在几种不同荷载模式下的应力-应变曲线，应变

47、曲线，这样可以选择出最合适的本构模型以及描述这种模型的参数。这样可以选择出最合适的本构模型以及描述这种模型的参数。 同一种橡胶材料的三种同一种橡胶材料的三种(sn zhn)拉伸变形状态的应力拉伸变形状态的应力-应变曲线图应变曲线图，对比试验曲线，由最小二乘法拟合多项式本构模型中的系数。对比试验曲线，由最小二乘法拟合多项式本构模型中的系数。试验试验(shyn)拟合本构模型系数拟合本构模型系数第41页/共94页第四十二页，共94页。试验拟合本构模型试验拟合本构模型(mxng)系数系数 给出实验数据，应力表达式的系数通过给出实验数据，应力表达式的系数通过(tnggu)最小二乘法拟合确定最小二乘法拟合

48、确定，这样可以使得误差最小。即对于，这样可以使得误差最小。即对于n 组应力组应力-应变的试验数据，取相对误差应变的试验数据，取相对误差E 的最小值，拟合应力表达式中的系数，得到理论本构模型。的最小值，拟合应力表达式中的系数，得到理论本构模型。按照本构关系与伸长率对应的应力表达式按照本构关系与伸长率对应的应力表达式 实验数据中的应力值实验数据中的应力值 第42页/共94页第四十三页，共94页。确定确定(qudng)材料常数的经验公式材料常数的经验公式 试验拟合试验拟合(n h)本构模型系数本构模型系数 对于已经成型的橡胶元件，通常不容易通过上述试验来确定对于已经成型的橡胶元件，通常不容易通过上述

49、试验来确定其材料常数。经验公式是通过橡胶的其材料常数。经验公式是通过橡胶的IRHD硬度指标来确定材料硬度指标来确定材料的弹性模量和切变的弹性模量和切变(qi bin)模量，再由材料常数和弹性模量的模量，再由材料常数和弹性模量的关系来确定材料常数。基本公式为（小应变条件）关系来确定材料常数。基本公式为（小应变条件）将得到的材料常数代入将得到的材料常数代入Mooney-Rivlin模型进行计算。模型进行计算。 例子例子 采用氢化丁腈橡胶采用氢化丁腈橡胶H-NBR75，硬度为，硬度为75MPa，解得，解得 第43页/共94页第四十四页，共94页。 由于大型有限元软件的迅速发展，使得复杂的超弹性模型计

50、算过程由由于大型有限元软件的迅速发展，使得复杂的超弹性模型计算过程由计算机程序完成，在计算机程序完成，在ABAQUS等商用软件中给出了具体的计算。用户要等商用软件中给出了具体的计算。用户要熟悉如何输入数据文件，根据试验数据拟合和选用合适的本构模型，如何熟悉如何输入数据文件，根据试验数据拟合和选用合适的本构模型，如何处理输出结果并检验其是否正确。对于初学者来说，商用软件是一个处理输出结果并检验其是否正确。对于初学者来说，商用软件是一个“黑黑匣子匣子”，因此，掌握超弹性材料模型理论和计算方法是取得仿真成功，因此，掌握超弹性材料模型理论和计算方法是取得仿真成功(chnggng)的关键。的关键。结论结

51、论(jiln)与讨论与讨论 需要注意的是，对于不可压缩材料的平面问题，无论是解析解还是数值解，需要注意的是，对于不可压缩材料的平面问题，无论是解析解还是数值解，均不能采用平面应变解答。因为对于不可压缩材料，如果采用平面应变模型均不能采用平面应变解答。因为对于不可压缩材料，如果采用平面应变模型(mxng)，其体积不变，内力为不确定量，在有限元中的节点位移不能反映单，其体积不变，内力为不确定量，在有限元中的节点位移不能反映单元内力的变化。对于不可压缩材料或者接近于不可压缩材料的平面问题，务必应元内力的变化。对于不可压缩材料或者接近于不可压缩材料的平面问题，务必应用平面应力（或者广义平面应变）解答。

52、用平面应力（或者广义平面应变）解答。第44页/共94页第四十五页，共94页。Part3钢Part2橡胶 RsPart1钢Rr b过盈面过盈面橡胶减震轴过盈配合橡胶减震轴过盈配合(pih)的解析解和有限元的解析解和有限元解解平面应变和平面应力模型平面应变和平面应力模型过盈量过盈量1.9mm ，应力，应力(yngl)非常大，非常大，原因是平面应变模型原因是平面应变模型橡胶和钢环的解析解与橡胶和钢环的解析解与FE解的径向解的径向(jn xin)应力比较应力比较 广义平面应变平面应力问题广义平面应变平面应力问题不发生体积自锁不发生体积自锁平面应变模型平面应变模型发生体积自锁发生体积自锁第45页/共94

53、页第四十六页，共94页。5 5 一维塑性一维塑性(sxng) (sxng) 应力保持40MPa的蠕变试验数据与计算结果对比 ),( TtEE 第46页/共94页第四十七页，共94页。123= s5 5 一维塑性一维塑性 第47页/共94页第四十八页，共94页。123= s允许允许(ynx)应力应力 5 5 一维塑性一维塑性(sxng) (sxng) 第48页/共94页第四十九页，共94页。5 5 一维塑性一维塑性(sxng) (sxng) 形状形状(xngzhun)改变比能与体积改变比能与体积改变比能改变比能第49页/共94页第五十页，共94页。体积改变体积改变(gibin)(gibin)能密

54、度与形状改变能密度与形状改变(gibin)(gibin)能密度能密度2 13 3215 5 一维塑性一维塑性(sxng) (sxng) 第50页/共94页第五十一页，共94页。123= s5 5 一维塑性一维塑性 第51页/共94页第五十二页，共94页。5 5 一维塑性一维塑性 第52页/共94页第五十三页，共94页。5 一维塑性一维塑性(sxng) 对于卸载后产生永久应变对于卸载后产生永久应变(yngbin)的材料称为塑性材料。的材料称为塑性材料。 应变的每一增量分解应变的每一增量分解(fnji)成为弹性可逆部分和塑性不可逆部分成为弹性可逆部分和塑性不可逆部分 塑性理论的主要内容有：塑性理论

55、的主要内容有： 屈服函数控制塑性变形的突变和连续，是内变量和应力的函数屈服函数控制塑性变形的突变和连续，是内变量和应力的函数 流动法则控制塑性流动，即确定塑性应变增量。流动法则控制塑性流动，即确定塑性应变增量。内部变量内部变量的演化方程控制屈服函数的演化，包括应变的演化方程控制屈服函数的演化，包括应变-硬化关系。硬化关系。 弹弹-塑性定律是路径相关和耗能的，大部分的功消耗在材料塑性塑性定律是路径相关和耗能的，大部分的功消耗在材料塑性变形中，不可逆换成其它形式的能量，特别是热。应力取决于整个变形中，不可逆换成其它形式的能量，特别是热。应力取决于整个变形的历史，不能表示成为应变的单值函数；而它仅能

56、指定作为应变形的历史，不能表示成为应变的单值函数；而它仅能指定作为应力和应变的率之间的关系。力和应变的率之间的关系。第53页/共94页第五十四页，共94页。5 一维塑性一维塑性(sxng) 一维率无关一维率无关(wgun)塑性塑性 典型典型(dinxng)弹弹-塑性材料的应力塑性材料的应力-应变曲线应变曲线 应变的增量假设分解成为弹性和塑性部分的和，率形式应变的增量假设分解成为弹性和塑性部分的和，率形式 应力增量应力增量( (率率) )总是与弹性模量和弹性应变的增量总是与弹性模量和弹性应变的增量( (率率) )有关有关 非线性弹非线性弹-塑性区段，应力塑性区段，应力-应变应变切线模量切线模量

57、应力应力-应变关系的是应变关系的是率率均匀的。如均匀的。如果被任意的果被任意的时间时间因子缩放，本构关系保持因子缩放，本构关系保持不变。因此，材料反应是不变。因此，材料反应是率无关率无关的。的。 第54页/共94页第五十五页，共94页。5 一维塑性一维塑性(sxng) 一维率无关一维率无关(wgun)塑性塑性 通过流动法则给出了塑性应变通过流动法则给出了塑性应变(yngbin)率，常常表示为塑性流动势能的形式率，常常表示为塑性流动势能的形式塑性率参数塑性率参数 流动势能的一个例子是流动势能的一个例子是 等效应力等效应力 屈服条件为屈服条件为 单轴拉伸的屈服强度单轴拉伸的屈服强度 等效塑性应变等

58、效塑性应变 材料在初始屈服之后屈服强度的增加称为功硬化或者应变硬化材料在初始屈服之后屈服强度的增加称为功硬化或者应变硬化( (对应于应变软化对应于应变软化) )。硬化行为一般是塑性变形先期历史的函数。硬化行为一般是塑性变形先期历史的函数。 屈服行为是各向同性硬化；拉伸和压缩的屈服强度总是相等。屈服行为是各向同性硬化；拉伸和压缩的屈服强度总是相等。第55页/共94页第五十六页，共94页。5 一维塑性一维塑性(sxng) 一维率无关一维率无关(wgun)塑性塑性 一个一个(y )特殊的模型，特殊的模型， 塑性应变率写成为塑性应变率写成为 p)(signf 塑性模型称为关联的，否则，塑性流动是非关联

59、的。塑性模型称为关联的，否则，塑性流动是非关联的。对于关联塑性，塑性流动是沿着屈服面的法线方向。对于关联塑性，塑性流动是沿着屈服面的法线方向。 由此看出由此看出仅当满足屈服条件仅当满足屈服条件时发生塑性变形。时发生塑性变形。 0)(Yf当塑性加载时，应力必须保持在屈服面上，当塑性加载时，应力必须保持在屈服面上， 实现了实现了一致性条件一致性条件 这给出这给出塑性模量塑性模量 第56页/共94页第五十七页，共94页。5 一维塑性一维塑性(sxng) 一维率无关一维率无关(wgun)塑性塑性 ddHY)(典型的硬化曲线，典型的硬化曲线，塑性模量塑性模量对应塑性加载和纯弹性加载或卸载对应塑性加载和纯

60、弹性加载或卸载(xi zi)，切线模量为，切线模量为 塑性转换参数塑性转换参数加载卸载条件还可以写为加载卸载条件还可以写为 一致性条件的率形式一致性条件的率形式 应力状态位于塑性表面应力状态位于塑性表面 塑性率参数非负塑性率参数非负 对于塑性加载对于塑性加载 必须保持在屈服面上必须保持在屈服面上 其应力状态其应力状态对于弹性加载或者卸载对于弹性加载或者卸载 没有塑性流动没有塑性流动 因此因此 第57页/共94页第五十八页，共94页。材料硬化描述材料硬化描述(mio sh) (a) Bauschinger效果效果 (b) 屈服面的平移和扩展屈服面的平移和扩展 在循环加载中，各向同性硬化模型提供了