资源描述
第第5 5章章 本构模型本构模型 1 1引言引言2 2应力应力-应变曲线应变曲线 3 3一维弹性一维弹性 4 4非线性弹性（超弹性）非线性弹性（超弹性）5 5一维塑性一维塑性 6 6多轴塑性多轴塑性7 7超弹塑性模型超弹塑性模型8 8粘弹性粘弹性 9 9应力更新算法应力更新算法1010连续介质力学与本构模型连续介质力学与本构模型1 1 引言引言 本本构构方方程程率率形形式式的的积积分分算算法法称称为为应应力力更更新新算算法法（也也称称为为本本构构更更新新算法），包括：算法），包括：径向返回算法的一类图形返回算法，径向返回算法的一类图形返回算法，算法模量与基本应力更新方案一致的概念，算法模量与基本应力更新方案一致的概念，大变形问题的增量客观应力更新方案，大变形问题的增量客观应力更新方案，基于弹性响应的应力更新方案，自动满足客观性的超弹性势能。基于弹性响应的应力更新方案，自动满足客观性的超弹性势能。为为了了进进行行分分析析，选选择择材材料料模模型型是是很很重重要要，往往往往又又不不是是很很明明确确，仅仅有有的的信信息息可可能能是是一一般般性性的的知知识识和和经经验验，即即可可能能是是材材料料行行为为的的几几条条应力应变曲线。应力应变曲线。在在有有限限元元软软件件库库中中选选择择合合适适的的本本构构模模型型，如如果果没没有有合合适适的的本本构构模模型型，要要开开发发用用户户材材料料子子程程序序。重重要要的的是是理理解解本本构构模模型型的的关关键键特特征征，创创建建模模型型的的假假设设，材材料料、荷荷载载和和变变形形域域、以以及及程程序序中中的的数数值值问问题题是是否适合模型。否适合模型。2 应力应力-应变曲线应变曲线 材料应力应变行为的许多基本特征可以从一维应力状态(单轴应力或者剪切)的一组应力-应变曲线中获得，多轴状态的本构方程常常基于在试验中观察到的一维行为而简单生成。载荷位移曲线载荷位移曲线 名义应力（工程应力）给出为名义应力（工程应力）给出为 定义伸长定义伸长 工程应变定义为工程应变定义为 2 2 应力应力-应变曲线应变曲线 CauchyCauchy（或者真实）应力（或者真实）应力表示为表示为 以以每单位当前长度应变的增量每单位当前长度应变的增量随长度随长度的变化得到另一种应变度量的变化得到另一种应变度量 对数应变（也称为真实应变）对数应变（也称为真实应变）对材料时间求导，表达式为对材料时间求导，表达式为一维情况，上式为变形率一维情况，上式为变形率 当前面积的表达式给出为当前面积的表达式给出为真实应力应变曲线 工程应力应变曲线2 2 应力应力-应变曲线应变曲线 考虑一种不可压缩材料考虑一种不可压缩材料(J J1)1)，名义应力和工程应变的，名义应力和工程应变的关系为关系为真实应力真实应力(对于不可压缩材料对于不可压缩材料)说说明明了了对对于于本本构构行行为为应应用用不不同同泛泛函函表表达达式式的的区区别别，对对于于同同样样材材料取决于采用何种应力和变形的度量。料取决于采用何种应力和变形的度量。应应力力应应变变曲曲线线的的显显著著特特征征之之一一是是非非线线性性的的度度。材材料料线线弹弹性性行行为的范围小于应变的百分之几，就可以采用小应变理论描述。为的范围小于应变的百分之几，就可以采用小应变理论描述。2 2 应力应力-应变曲线应变曲线 应应力力应应变变反反应应与与变变形形率率无无关关的的材材料料称称为为率率无无关关；否否则则，称为称为率相关率相关。名义应变率定义为。名义应变率定义为 率无关和率相关材料的一维反应因为 和即名义应变率等于伸长率，例如即名义应变率等于伸长率，例如 可可以以看看出出，对对于于率率无无关关材材料料的的应应力力应应变变曲曲线线是是应应变变率率独独立立的的，而而对对于于率率相相关关材材料料的的应应力力应应变变曲曲线线，当当应应变变率率提提高高时时是是上上升升的的；而而当温度升高时是下降的。当温度升高时是下降的。2 2 应力应力-应变曲线应变曲线 对对于于弹弹性性材材料料，应应力力应应变变的的卸卸载载曲曲线线简简单单地地沿沿加加载载曲曲线线返返回回，直直到到完完全全卸卸载载，材材料料返返回回到到了了它它的的初初始始未未伸伸长长状状态态。然然而而，对对于于弹弹塑塑性性材材料料，卸卸载载曲曲线线区区别别于于加加载载曲曲线线，卸卸载载曲曲线线的的斜斜率率是是典典型型的的应应力力应应变变弹弹性性（初初始始）段段的的斜斜率率，卸卸载载后后产产生生永永久久应应变变。其其它它材材料料的的行行为为介介于于这这两两种种极极端端之之间间。由由于于在在加加载载过过程程中中微微裂裂纹纹的的形形成成材材料料已已经经损损伤伤，脆脆性性材材料料的的卸卸载载行行为为，当当荷荷载载移移去去后后微微裂裂纹纹闭闭合合，弹弹性性应应变变得得到到恢恢复复。卸卸载载曲曲线线的的初初始始斜斜率率给给出出形形成微裂纹损伤程度的信息。成微裂纹损伤程度的信息。(a)弹性,(b)弹-塑性,(c)弹性含损伤 3 3 一维弹性一维弹性 弹弹性性材材料料的的基基本本性性能能是是应应力力仅仅依依赖赖于于应应变变的的当当前前水水平平。这这意意味味着着加加载载和和卸卸载载的的应应力力-应应变变曲曲线线是是一一致致的的，当当卸卸载载结结束束时时材材料料恢恢复复到到初初始始状状态态。称称这这种种应应变变是是可可逆逆的的。而而且且，弹弹性性材材料料是是率率无无关关的的(与应变率无关与应变率无关)。弹性材料的应力和应变是一一对应的。弹性材料的应力和应变是一一对应的。小应变小应变 可可逆逆和和路路径径无无关关默默认认在在变变形形中中没没有有能能量量耗耗散散，在在弹弹性性材材料料中中，储存在物体中的能量全部消耗在变形中，卸载后材料恢复。储存在物体中的能量全部消耗在变形中，卸载后材料恢复。对于一维弹性材料，对于一维弹性材料，可逆、路径无关、无能量耗散可逆、路径无关、无能量耗散是等价的特征。是等价的特征。对于二维和三维弹性，以及超弹性材料，也类似。对于二维和三维弹性，以及超弹性材料，也类似。对于任意应变，不管如何达到应变值，上式给出唯一应力值。对于任意应变，不管如何达到应变值，上式给出唯一应力值。3 3 一维弹性一维弹性 应变能一般是应变的凸函数，例如，应变能一般是应变的凸函数，例如，(a)凸应变能函数(b)应力应变曲线 当当 公式的等号成立。公式的等号成立。凸应变能函数的一个例子如图所示。在这种情况下，函数凸应变能函数的一个例子如图所示。在这种情况下，函数是单调递增的，如果是单调递增的，如果w 是非凸函数，则是非凸函数，则 s 先增后减，材料应变先增后减，材料应变软化，这是非稳定的材料反应，软化，这是非稳定的材料反应，如右下图。如右下图。(a)非凸应变能函数(b)相应的应力应变曲线 大应变大应变 从从弹弹性性推推广广到到大大应应变变，只只要要选选择择应应变变度度量量和和定定义义应应力力（功功共共轭轭）的的弹弹性性势势能能。势势能能的的存存在在是是默默认认了了可可逆逆、路路径径无无关关和和无无能能量量耗散耗散。如。如 3 3 一维弹性一维弹性 在在弹弹性性应应力力-应应变变关关系系中中，从从应应变变的的势势函函数数可可以以获获得得应应力力为为超弹性。如一维大应变问题，以超弹性。如一维大应变问题，以GreenGreen应变的二次函数表示应变的二次函数表示 对于小应变问题，即为胡克定律。对于小应变问题，即为胡克定律。大应变大应变 一一种种材材料料的的CauchyCauchy应应力力率率与与变变形形率率相相关关，称称为为次次弹弹性性。这这种关系一般是非线性的，给出为种关系一般是非线性的，给出为 3 3 一维弹性一维弹性 一个特殊的线性次弹性关系给出为一个特殊的线性次弹性关系给出为 这这是是与与路路径径无无关关的的超超弹弹性性关关系系。对对于于多多轴轴问问题题，一一般般次次弹弹性性关关系系不不能能转转换换到到超超弹弹性性，它它仅仅在在一一维维情情况况下下是是严严格格路路径径无无关关的的。然然而而，如如果果是是弹弹性性小小应应变变，其其行行为为足足以以接接近近路路径径无无关关的的弹弹性性行行为为。因因为为次次弹弹性性的的简简单单性性，公公式式(5.3.11)(5.3.11)的的多多轴轴一一般般形形式式常常常常应应用用在在有限元软件中，以模拟大应变弹塑性的弹性反应。有限元软件中，以模拟大应变弹塑性的弹性反应。对上式的关系积分，得到对上式的关系积分，得到4 4 非线性弹性非线性弹性 对对于于有有限限应应变变有有许许多多不不同同的的应应力力和和变变形形度度量量，同同样样的的本本构构关关系系可可以以写写成成几几种种不不同同的的形形式式，总总是是可可能能从从一一种种形形式式的的本本构构关关系转换到另一种形式。系转换到另一种形式。大大应应变变弹弹性性本本构构模模型型首首先先表表述述成成KirchhoffKirchhoff材材料料的的一一种种特特殊殊形形式式，由由线线弹弹性性直直接接生生成成到到大大变变形形。满满足足路路径径无无关关、可可逆逆和和无无能能量量耗耗散散。因因此此，路路径径无无关关的的程程度度可可以以视视为为材材料料模模型型弹弹性性的的度度量。量。次次弹弹性性材材料料是是路路径径无无关关程程度度最最弱弱的的材材料料，遵遵从从CauchyCauchy弹弹性性，其应力是路径无关的，但是其能量不是路径无关的。其应力是路径无关的，但是其能量不是路径无关的。超超弹弹性性材材料料或或者者GreenGreen弹弹性性，它它是是路路径径无无关关和和完完全全可可逆逆的的，应力由应变势能导出。应力由应变势能导出。4 4 非线性弹性非线性弹性 小应变和大转动小应变和大转动 式式中中 C C 为为弹弹性性模模量量(切切线线模模量量)的的四四阶阶张张量量，对对KirchhoffKirchhoff材材料料是是常常数数，代代表表了了应应力力和和应应变变的的多多轴轴状状态态。它它可可以以完完全全反反映映材材料料的的各各向异性。向异性。许许多多工工程程应应用用包包括括小小应应变变和和大大转转动动。在在这这些些问问题题中中，大大变变形形的的效效果果主主要要来来自自于于大大转转动动，如如直直升升机机旋旋翼翼、船船上上升升降降器器或或者者钓钓鱼鱼杆杆的的弯弯曲曲。由由线线弹弹性性定定律律的的简简单单扩扩展展即即可可以以模模拟拟材材料料的的反反应应，但但要要以以PK2PK2应应力力代代替替其其中中的的应应力力和和以以GreenGreen应应变变代代替替线线性性应应变变，这这称称为为Saint-Venant-Saint-Venant-KirchhoffKirchhoff材材料料，或或者者简简称称为为KirchhoffKirchhoff材材料料。最一般的最一般的KirchhoffKirchhoff模型为模型为4 4 非线性弹性非线性弹性 式中式中C C为弹性模量的四阶张量，有为弹性模量的四阶张量，有8181个常数。利用对称性可以显个常数。利用对称性可以显著地减少常数。著地减少常数。一般的四阶张量有一般的四阶张量有3 34 48181个独立常数，与全应力张量的个独立常数，与全应力张量的9 9个个分量和全应变张量的分量和全应变张量的9 9个分量有关。个分量有关。如次弹性本构方程如次弹性本构方程这这样样C C为为对对称称矩矩阵阵（主主对对称称性性），在在8181个个常常数数中中有有4545个个是是独独立立的的。成为上三角或下三角矩阵。成为上三角或下三角矩阵。4 4 非线性弹性非线性弹性 利用势能表示的应力应变关系和利用势能表示的应力应变关系和GreenGreen公式，公式，故有故有 应力张量和应变张量均为对称张量（应力张量和应变张量均为对称张量（次对称性次对称性），即），即 4 4 非线性弹性非线性弹性 应力张量和应变张量均为对称张量（应力张量和应变张量均为对称张量（次对称性次对称性），即），即 再再利利用用模模量量的的主主对对称称性性使使独独立立弹弹性性常常数数的的数数目目减减少少，由由3636个个常常数数减减少少为为2121个，为个，为各向异性材料各向异性材料。应应力力和和应应变变张张量量的的对对称称性性要要求求应应力力的的6 6个个独独立立分分量量仅仅与与应应变变的的6 6个个独独立立分分量量有有关关，由由弹弹性性模模量量的的局局部部对对称称结结果果，独独立立常常数数的的数数目目减少到减少到3636个。个。4 4 非线性弹性非线性弹性 写成矩阵形式为（可以是上或下三角矩阵）写成矩阵形式为（可以是上或下三角矩阵）对对于于正正交交各各向向异异性性，具具有有正正交交的的三三个个弹弹性性对对称称面面，当当坐坐标标变变号号，为使应变能密度不变，有为使应变能密度不变，有 这样由这样由2121个常数减少为个常数减少为1414个，为个，为正交各向异性材料正交各向异性材料。若材料对称坐标平面，当沿轴平面反射时，弹性模量不变，若材料对称坐标平面，当沿轴平面反射时，弹性模量不变，固为固为正交各向异性体正交各向异性体，有，有 对对于于一一个个由由三三个个彼彼此此正正交交的的对对称称平平面面组组成成的的正正交交材材料料(如如木木材材或或纤纤维维增增强强的的复复合合材材料料)，仅仅有有9 9个个独独立立弹弹性性常常数数，KirchhoffKirchhoff应力应变关系为材料对称坐标平面，为应力应变关系为材料对称坐标平面，为正交各向异性体正交各向异性体4 4 非线性弹性非线性弹性 对于对于各向同性材料各向同性材料，仅有，仅有3 3个常数个常数 4 4 非线性弹性非线性弹性 小应变和大转动小应变和大转动 对于各向同性的于各向同性的Kirchhoff材料，其材料，其应力力应变关系可以写成关系可以写成为式中式中Lam常数，体常数，体积模量模量K，杨氏模量氏模量 E和泊松比和泊松比的关系为的关系为 材材料料对对称称的的一一个个重重要要的的例例子子是是各各向向同同性性。一一个个各各向向同同性性材材料料没没有有方方位位或或者者方方向向的的选选择择，因因此此，当当以以任任何何直直角角坐坐标标系系表表示示的的应应力力应应变变关关系系是是等等同同的的。对对于于小小应应变变的的许许多多材材料料(如如金金属属和和陶陶瓷瓷)可可以以作作为为各各向向同同性性进进行行模模拟拟。张张量量C C是是各各向向同同性性的的。在在任任何何坐坐标标系统中，一个各向同性张量有相同的分量。系统中，一个各向同性张量有相同的分量。(克罗内克)符号构成的一个线性组合：4 4 非线性弹性非线性弹性 不可压缩性不可压缩性 在变形的过程中，不可压缩材料的体积不变，密度保持常数。在变形的过程中，不可压缩材料的体积不变，密度保持常数。不可压缩材料的运动称为等体积运动。不可压缩材料的运动称为等体积运动。总体变形总体变形 等体积约束运动的率形式等体积约束运动的率形式 将应力和应变率度量写成偏量和静水（体积的）部分的和，将应力和应变率度量写成偏量和静水（体积的）部分的和，对于不可压缩材料，静水部分也称为张量的球形部分，分解式为：对于不可压缩材料，静水部分也称为张量的球形部分，分解式为：对于不可压缩材料，压力不能从本构方程确定，而是从动量对于不可压缩材料，压力不能从本构方程确定，而是从动量方程确定。方程确定。4 非线性弹性非线性弹性 Kirchhoff应力应力 由由Jacobian行列式放大，称它为权重行列式放大，称它为权重Cauchy应力。对于应力。对于等体积等体积运动运动，它等同于，它等同于Cauchy应力。应力。次弹性次弹性次弹性材料规律联系应力率和变形率。次弹性材料规律联系应力率和变形率。上式是率无关、线性增加和可逆的。对于有限变形状态的微小增量，应力和应变的增量是线性关系，当卸载后可以恢复。然而，对于大变形能量不一定必须守恒，并且在闭合变形轨迹上作的功不一定必须为零。次弹性规律主要用来代表在弹-塑性规律中的弹性反应，小变形弹性，且耗能效果也小。4 非线性弹性非线性弹性 切线模量之间的关系切线模量之间的关系 对于各向同性材料对于各向同性材料Jaumann率的切线模量为率的切线模量为 某些次弹性本构关系共同应用的形式为某些次弹性本构关系共同应用的形式为对于同一种材料，切线模量不同，材料反应的率形式不同，如对于同一种材料，切线模量不同，材料反应的率形式不同，如 如果如果是常数，是常数，不是常数。不是常数。切切线模量模量证明明见第第5.4.5节，推，推导复复杂4 非线性弹性非线性弹性 超弹性材料超弹性材料 平衡方程是以物体中应力的形式建立的，应力来源于变形，如应变。如果本构行为仅是变形的当前状态的函数，为与时间无关的弹性本构。而对于接近不可压缩的材料，仅依赖变形（应变）不一定能够得到应力。储存在材料中的能量（功）仅取决于变形的初始和最终状态，并且是独立于变形（或荷载）路径，称这种弹性材料为超弹性（hyper-elastic）材料，或者为Green弹性，例如常用的工业橡胶。动物的肌肉也具有超弹性的力学性质。这里主要讨论橡胶材料的超弹性力学行为。4 非线性弹性非线性弹性 超弹性材料超弹性材料 对于功独立于荷载路径的弹性材料称之为超弹性（对于功独立于荷载路径的弹性材料称之为超弹性（Green弹弹性）材料。超弹性材料的特征是存在一个潜在（或应变）能量性）材料。超弹性材料的特征是存在一个潜在（或应变）能量函数，它是应力的势能：函数，它是应力的势能：通过适当转换获得了对于不同应力度量的表达式通过适当转换获得了对于不同应力度量的表达式 由于变形梯度张量由于变形梯度张量F是不对称的，因此名义应力张量是不对称的，因此名义应力张量P的的9个个分量是不对称的。分量是不对称的。在橡胶大变形中应用在橡胶大变形中应用多项式模型多项式模型和和Ogden指数模型指数模型。4 非线性弹性非线性弹性 超弹性材料超弹性材料 目目前前，世世界界半半数数以以上上的的橡橡胶胶是是合合成成橡橡胶胶。合合成成橡橡胶胶的的种种类类很很多多，例例如如，制制造造轮轮胎胎使使用用的的丁丁苯苯橡橡胶胶（苯苯乙乙烯烯和和丁丁二二烯烯的的共共聚聚物物）或或乙乙丙丙烯烯橡橡胶胶（ERP）；用用于于汽汽车车配配件件的的有有氯氯丁丁橡橡胶胶及及另另一种具有天然橡胶各种性能的异戊橡胶。一种具有天然橡胶各种性能的异戊橡胶。在在众众多多的的合合成成橡橡胶胶中中，硅硅橡橡胶胶是是其其中中的的佼佼佼佼者者。它它具具有有无无味味无无毒毒，不不怕怕高高温温和和严严寒寒的的特特点点，在在摄摄氏氏300度度和和零零下下90度度时时能能够够“泰泰然然自自若若”、“面面不不改改色色”，仍仍不不失失原原有有的的强强度度和和弹弹性性。例例如生物材料。如生物材料。橡橡胶胶是是提提取取橡橡胶胶树树、橡橡胶胶草草等等植植物物的的胶胶乳乳，加加工工后后制制成成的的具具有有弹弹性性、绝绝缘缘性性、不不透透水水和和空空气气的的材材料料。在在半半个个世世纪纪前前，“橡橡胶胶”一一词词是是专专指指生生橡橡胶胶，它它是是从从热热带带植植物物巴巴西西三三叶叶胶胶的的胶胶乳乳提提炼出来的。炼出来的。4 非线性弹性非线性弹性 超弹性材料超弹性材料 1839年年，Charle Goodyear发发明明了了橡橡胶胶的的硫硫化化方方法法，其其姓姓氏现在已经成为国际上著名橡胶轮胎的商标。氏现在已经成为国际上著名橡胶轮胎的商标。从从19世世纪纪中中叶叶起起橡橡胶胶就就成成为为一一种种重重要要的的工工程程材材料料。然然而而，橡橡胶胶材材料料的的行行为为复复杂杂，不不同同于于金金属属材材料料仅仅需需要要几几个个参参数数就就可可以以描描述述材材料料特特性性。橡橡胶胶材材料料受受力力以以后后，变变形形是是伴伴随随着着大大位位移移和和大大应应变变，其其本本构构关关系系是是非非线线性性的的，并并且且在在变变形形过过程程中中体体积积几几乎乎保保持不变。持不变。橡胶具有许多特殊的性能，例如电绝缘性、耐氧老化性、耐橡胶具有许多特殊的性能，例如电绝缘性、耐氧老化性、耐光老化性、防霉性、化学稳定性等。光老化性、防霉性、化学稳定性等。4 非线性弹性非线性弹性 超弹性材料超弹性材料 由由于于计计算算机机以以及及有有限限元元数数值值分分析析的的飞飞速速发发展展，我我们们可可以以借借助助计计算算机机来来对对超超弹弹性性材材料料的的工工程程应应用用进进行行深深入入研研究究以以及及优优化化设设计计。可可以以用用有有限限元元等等数数值值方方法法来来计计算算分分析析橡橡胶胶元元件件的的力力学学性性能能，包包括括选选取取和和拟拟合合橡橡胶胶的的本本构构模模型型，以以及及用用有有限限元元建建模模和和处处理理计计算算结果等。结果等。橡橡胶胶是是一一种种弹弹性性聚聚合合物物，其其特特点点是是有有很很强强的的非非线线性性粘粘弹弹性性行行为为。它它的的力力学学行行为为对对温温度度、环环境境、应应变变历历史史、加加载载速速率率都都非非常常敏敏感感，这这样样使使得得描描述述橡橡胶胶的的行行为为变变得得非非常常复复杂杂。橡橡胶胶的的制制造造工工艺和成分也对橡胶的力学性能有着显著的影响。艺和成分也对橡胶的力学性能有着显著的影响。固体橡胶材料的拉伸试验曲线与材料演化模型固体橡胶材料的拉伸试验曲线与材料演化模型 固体橡胶是几乎不可压缩的，其泊松比接近于固体橡胶是几乎不可压缩的，其泊松比接近于0.5。可逆，可逆，大应变。大应变。初始各向同性，应变增加后分子定向排列。初始各向同性，应变增加后分子定向排列。4 非线性弹性非线性弹性 超弹性材料超弹性材料 常用的橡胶性态可分为常用的橡胶性态可分为固体橡胶固体橡胶和和泡沫橡胶泡沫橡胶。4 非线性弹性非线性弹性 超弹性材料超弹性材料 一一般般将将多多孔孔橡橡胶胶或或弹弹性性泡泡沫沫材材料料统统称称为为泡泡沫沫材材料料。弹弹性性泡泡沫材料的普通例子有多孔聚合物，如海绵、包装材料等。沫材料的普通例子有多孔聚合物，如海绵、包装材料等。泡泡沫沫橡橡胶胶是是由由橡橡胶胶制制成成的的弹弹性性泡泡沫沫材材料料，能能够够满满足足非非常常大大的的弹弹性性应应变变要要求求，拉拉伸伸时时的的应应变变可可以以达达到到500或或更更大大，压压缩缩时时的的应应变变可可以以达达到到90或或更更小小。与与固固体体橡橡胶胶的的几几乎乎不不可可压压缩缩性性相相比比，泡泡沫沫材材料料的的多多孔孔性性则则允允许许非非常常大大的的体体积积缩缩小小变变形形，因因此此具有良好的能量吸收性。具有良好的能量吸收性。泡沫橡胶材料的多面体微元模型泡沫橡胶材料的多面体微元模型 a)开放腔室，开放腔室，b)封闭腔室封闭腔室4 非线性弹性非线性弹性 超弹性材料超弹性材料 泡沫橡胶材料的应力泡沫橡胶材料的应力-应变曲线应变曲线 a)压缩压缩 b)拉伸拉伸小应变小应变 5%，线弹性，泊松比为，线弹性，泊松比为0.3。大应变，压缩时，泊松比为大应变，压缩时，泊松比为0.0；拉伸时，泊松比大于拉伸时，泊松比大于0.0。典型固体橡胶材料单轴拉伸应力典型固体橡胶材料单轴拉伸应力-应变曲线应变曲线 橡胶本构模型橡胶本构模型 4 非线性弹性非线性弹性 小变形小变形 以多项式形式本构模型为例，其应变能密度表达式为以多项式形式本构模型为例，其应变能密度表达式为忽略二阶及二阶以上小量，变为忽略二阶及二阶以上小量，变为弹性常数为弹性常数为 当当 橡胶本构模型橡胶本构模型 4 非线性弹性非线性弹性 定义伸长定义伸长 工程应变定义为工程应变定义为 二阶张量基本不变量二阶张量基本不变量 小变形，有小变形，有 小变形小变形 橡胶本构模型橡胶本构模型 4 非线性弹性非线性弹性 例题例题 在超弹性计算中，橡胶使用三次减缩多项式应变能本构模型，在超弹性计算中，橡胶使用三次减缩多项式应变能本构模型，应变能密度表达式为应变能密度表达式为若取若取（单位为（单位为MPa），求材料弹性常数。），求材料弹性常数。利用公式利用公式解：解：解出橡胶的弹性常数为解出橡胶的弹性常数为 ，E=1.384MPa，=0.5 小变形小变形 橡胶本构模型橡胶本构模型 4 非线性弹性非线性弹性 常常用用的的橡橡胶胶力力学学性性能能描描述述方方法法主主要要分分为为两两类类，一一类类是是基基于于热热力力学学统统计计的的方方法法，另另一一类类是是基基于于橡橡胶胶为为连连续续介介质质的的唯唯象象学学描描述方法。述方法。热热力力学学统统计计方方法法的的基基础础为为观观察察到到橡橡胶胶中中的的弹弹性性恢恢复复力力主主要要来来自自熵熵的的减减少少。橡橡胶胶在在承承受受荷荷载载时时分分子子结结构构无无序序，熵熵的的减减少少是是由由于于橡橡胶胶伸伸长长使使得得橡橡胶胶结结构构由由高高度度无无序序变变得得有有序序。由由对对橡橡胶胶中中分子链的长度、方向以及结构的统计得到本构关系。分子链的长度、方向以及结构的统计得到本构关系。橡胶本构模型橡胶本构模型 唯唯象象学学描描述述方方法法假假设设在在未未变变形形状状态态下下橡橡胶胶为为各各向向同同性性材材料料，即即长长分分子子链链方方向向在在橡橡胶胶中中是是随随机机分分布布的的。这这种种各各向向同同性性的的假假设设是是用用单单位位体体积积（弹弹性性）应应变变能能函函数数（U）来来描描述述橡橡胶胶特特性性的的基基础，其本构模型为多项式形式模型和础，其本构模型为多项式形式模型和Ogden形式模型。形式模型。典型的典型的本构模型本构模型为多项式形式，其应变能密度表达式为为多项式形式，其应变能密度表达式为特殊形式可以由设定某些参数为特殊形式可以由设定某些参数为0来得到。如果所有来得到。如果所有 则得到减缩多项式模型则得到减缩多项式模型 对于完全多于完全多项式式，如果，如果,则只有只有线性部分的性部分的应变能量，能量，即即Mooney-Rivlin形式形式橡胶本构模型橡胶本构模型，则得到，则得到Neo-Hookean形式形式 对于减缩多项式对于减缩多项式，如果，如果 Mooney-Rivlin形式和形式和Neo-Hooken形式本构模型形式本构模型（后者是将（后者是将Hooke定律扩展至大变形）定律扩展至大变形）橡胶本构模型橡胶本构模型 Yeoh形式本构模型是形式本构模型是 时减缩多项式的特殊形式时减缩多项式的特殊形式 典型的典型的S形橡胶应力形橡胶应力-应变曲线应变曲线，C10正值，在小变形时为切线模量；正值，在小变形时为切线模量；C20为负值，中等变形时软化；为负值，中等变形时软化；C30正值，大变形时硬化。正值，大变形时硬化。橡胶本构模型橡胶本构模型 Ogden形式本构模型形式本构模型 Arruda-Boyce形式本构模型形式本构模型 Van der Waals模型模型 橡胶本构模型橡胶本构模型 其他形式的本构模型有：其他形式的本构模型有：试验拟合本构模型系数试验拟合本构模型系数 橡橡胶胶类类材材料料的的本本构构关关系系除除具具有有超超弹弹性性、大大变变形形的的特特征征外外，其其本本构构关关系系与与生生产产加加工工过过程程有有直直接接关关系系，如如橡橡胶胶配配方方和和硫硫化化工工艺艺。确确定定每每一一批批新新加加工工出出来来的的橡橡胶胶的的本本构构关关系系，都都要要依依赖赖于于精精确和充分的橡胶试验。确和充分的橡胶试验。通常在试验中应该测得在几种不同荷载模式下的应力通常在试验中应该测得在几种不同荷载模式下的应力-应变曲线，应变曲线，这样可以选择出最合适的本构模型以及描述这种模型的参数。这样可以选择出最合适的本构模型以及描述这种模型的参数。同一种橡胶材料的三种拉伸变形状态的应力同一种橡胶材料的三种拉伸变形状态的应力-应变曲线图，应变曲线图，对比试验曲线，由最小二乘法拟合多项式本构模型中的系数。对比试验曲线，由最小二乘法拟合多项式本构模型中的系数。试验拟合本构模型系数试验拟合本构模型系数试验拟合本构模型系数试验拟合本构模型系数 给给出出实实验验数数据据，应应力力表表达达式式的的系系数数通通过过最最小小二二乘乘法法拟拟合合确确定定，这这样样可可以以使使得得误误差差最最小小。即即对对于于n 组组应应力力-应应变变的的试试验验数数据据，取取相相对对误误差差E 的的最最小小值值，拟拟合合应应力力表表达达式式中中的的系系数数，得得到到理论本构模型。理论本构模型。按照本构关系与伸长率对应的应力表达式按照本构关系与伸长率对应的应力表达式 实验数据中的应力值实验数据中的应力值 确定材料常数的经验公式确定材料常数的经验公式 试验拟合本构模型系数试验拟合本构模型系数 对对于于已已经经成成型型的的橡橡胶胶元元件件，通通常常不不容容易易通通过过上上述述试试验验来来确确定定其其材材料料常常数数。经经验验公公式式是是通通过过橡橡胶胶的的IRHD硬硬度度指指标标来来确确定定材材料料的的弹弹性性模模量量和和切切变变模模量量，再再由由材材料料常常数数和和弹弹性性模模量量的的关关系系来来确确定材料常数。基本公式为（小应变条件）定材料常数。基本公式为（小应变条件）将得到的材料常数代入将得到的材料常数代入Mooney-Rivlin模型进行计算。模型进行计算。例子例子 采用氢化丁腈橡胶采用氢化丁腈橡胶H-NBR75，硬度为，硬度为75MPa，解得，解得 由由于于大大型型有有限限元元软软件件的的迅迅速速发发展展，使使得得复复杂杂的的超超弹弹性性模模型型计计算算过过程程由由计计算算机机程程序序完完成成，在在ABAQUS等等商商用用软软件件中中给给出出了了具具体体的的计计算算。用用户户要要熟熟悉悉如如何何输输入入数数据据文文件件，根根据据试试验验数数据据拟拟合合和和选选用用合合适适的的本本构构模模型型，如如何何处处理理输输出出结结果果并并检检验验其其是是否否正正确确。对对于于初初学学者者来来说说，商商用用软软件件是是一一个个“黑黑匣匣子子”，因因此此，掌掌握握超超弹弹性性材材料模型理论和计算方法是取得仿真成功的关键。料模型理论和计算方法是取得仿真成功的关键。结论与讨论结论与讨论 需需要要注注意意的的是是，对对于于不不可可压压缩缩材材料料的的平平面面问问题题，无无论论是是解解析析解解还还是是数数值值解解，均均不不能能采采用用平平面面应应变变解解答答。因因为为对对于于不不可可压压缩缩材材料料，如如果果采采用用平平面面应应变变模模型型，其其体体积积不不变变，内内力力为为不不确确定定量量，在在有有限限元元中中的的节节点点位位移移不不能能反反映映单单元元内内力力的的变变化化。对对于于不不可可压压缩缩材材料料或或者者接接近近于于不不可可压压缩缩材材料料的的平平面面问问题题，务务必必应应用用平平面面应应力力（或或者广义平面应变）解答者广义平面应变）解答。Part3钢Part2橡胶 RsPart1钢Rr b过盈面过盈面橡胶减震轴过盈配合的解析解和有限元解橡胶减震轴过盈配合的解析解和有限元解平面应变和平面应力模型平面应变和平面应力模型过盈量过盈量1.9mm，应力非常大，应力非常大，原因是平面应变模型原因是平面应变模型橡胶和钢环的解析解与橡胶和钢环的解析解与FE解的径向应力比较解的径向应力比较 广义平面应变平面应力问题广义平面应变平面应力问题不发生体积自锁不发生体积自锁平面应变模型平面应变模型发生体积自锁发生体积自锁问题：问题：问题：问题：在在在在有有有有限限限限元元元元力力力力学学学学模模模模型型型型中中中中，加加加加载载载载是是是是任任任任意意意意的的的的（如如如如三三三三维维维维），材材材材料料料料实实实实验验验验数数数数据据据据是是是是单单单单轴轴轴轴拉拉拉拉伸伸伸伸（如如如如一一一一维维维维），如如如如何何何何在在在在有有有有限限限限元元元元计计计计算算算算中中中中建建建建立立立立联联联联系系系系，实实实实现现现现对对对对应应应应的的的的应应应应力力力力状状状状态态态态，直直直直到到到到发发发发生生生生屈屈屈屈服和破坏？服和破坏？服和破坏？服和破坏？5 5 一维塑性一维塑性 从屈服准则的建立来从屈服准则的建立来从屈服准则的建立来从屈服准则的建立来回答这样的问题。回答这样的问题。回答这样的问题。回答这样的问题。应力保持40MPa的蠕变试验数据与计算结果对比 最大切应力屈服准则最大切应力屈服准则(Trescas Criterion)Trescas Criterion)无无无无论论论论材材材材料料料料处处处处于于于于什什什什么么么么应应应应力力力力状状状状态态态态，只只只只要要要要发发发发生生生生屈屈屈屈服服服服，都都都都是是是是由由由由于于于于微微微微元元元元内内内内的的的的最最最最大大大大切切切切应应应应力力力力达达达达到到到到了了了了某某某某一一一一共同的极限值。共同的极限值。共同的极限值。共同的极限值。1 1 2 2 3 3=s s拉伸屈服拉伸屈服拉伸屈服拉伸屈服试验确定试验确定试验确定试验确定任意状任意状任意状任意状态应力态应力态应力态应力5 5 一维塑性一维塑性 1 1 2 2 3 3=s s失效判据失效判据设计准则允许应力允许应力 5 5 一维塑性一维塑性 在在在在有有有有限限限限元元元元计计计计算算算算中中中中，材材材材料料料料的的的的应应应应力力力力和和和和应应应应变变变变状状状状态态态态等等等等价价价价于于于于单单单单轴轴轴轴拉拉拉拉伸伸伸伸实实实实验验验验数数数数据据据据的的的的对对对对应应应应值值值值，与与与与加加加加载载载载历历历历史史史史相相相相关关关关，只只只只要要要要发发发发生生生生屈屈屈屈服服服服，都都都都是是是是由由由由于于于于单单单单元元元元内内内内的的的的最最最最大大大大切切切切应应应应力力力力达达达达到到到到了了了了某某某某一一一一共共共共同同同同的的的的极限值。极限值。极限值。极限值。形状改变比能准则形状改变比能准则(Misess Criterion)Misess Criterion)无论材料处于什么应力状态，只要发生屈服，无论材料处于什么应力状态，只要发生屈服，无论材料处于什么应力状态，只要发生屈服，无论材料处于什么应力状态，只要发生屈服，都是由于微元的形状改变比能达到了一个共同的都是由于微元的形状改变比能达到了一个共同的都是由于微元的形状改变比能达到了一个共同的都是由于微元的形状改变比能达到了一个共同的极限值。极限值。极限值。极限值。5 5 一维塑性一维塑性 形状形状改变比能与改变比能与体积体积改变比能改变比能体积改变能密度与形状改变能密度体积改变能密度与形状改变能密度+5 5 一维塑性一维塑性 形状改变比能准则形状改变比能准则 1 1 2 2 3 3=s s单向应力单向应力三向应力三向应力5 5 一维塑性一维塑性 形状改变比能准则失效判据失效判据失效判据失效判据设计准则设计准则设计准则设计准则5 5 一维塑性一维塑性 5 一维塑性一维塑性 对于卸载后产生永久应变的材料称为塑性材料。对于卸载后产生永久应变的材料称为塑性材料。应变的每一增量分解成为弹性可逆部分和塑性不可逆部分应变的每一增量分解成为弹性可逆部分和塑性不可逆部分 塑性理论的主要内容有：塑性理论的主要内容有：屈服函数控制塑性屈服函数控制塑性变形的突形的突变和和连续，是内，是内变量和量和应力的函数力的函数 流动法则控制塑性流动，即确定塑性应变增量。流动法则控制塑性流动，即确定塑性应变增量。内部变量内部变量的演化方程控制屈服函数的演化，包括应变的演化方程控制屈服函数的演化，包括应变-硬化关系。硬化关系。弹弹-塑塑性性定定律律是是路路径径相相关关和和耗耗能能的的，大大部部分分的的功功消消耗耗在在材材料料塑塑性性变变形形中中，不不可可逆逆换换成成其其它它形形式式的的能能量量，特特别别是是热热。应应力力取取决决于于整整个个变变形形的的历历史史，不不能能表表示示成成为为应应变变的的单单值值函函数数；而而它它仅仅能能指指定定作作为为应应力力和应变的率之间的关系。和应变的率之间的关系。5 一维塑性一维塑性 一维率无关塑性一维率无关塑性 典型弹典型弹-塑性材料的应力塑性材料的应力-应变曲线应变曲线 应变的增量假设分解成为弹性和塑性部分的和，率形式应变的增量假设分解成为弹性和塑性部分的和，率形式 应力增量应力增量(率率)总是与弹性模量和弹性应变的增量总是与弹性模量和弹性应变的增量(率率)有关有关 非线性弹非线性弹-塑性区段，应力塑性区段，应力-应变应变切线模量切线模量 应应力力-应应变变关关系系的的是是率率均均匀匀的的。如如果果被被任任意意的的时时间间因因子子缩缩放放，本本构构关关系系保保持持不不变变。因因此此，材材料料反反应是应是率无关率无关的。的。5 一维塑性一维塑性 一维率无关塑性一维率无关塑性 通通过流流动法法则给出了塑性出了塑性应变率，常常表示率，常常表示为塑性流塑性流动势能的形式能的形式塑性率参数塑性率参数 流动势能的一个例子是流动势能的一个例子是 等效应力等效应力 屈服条件为屈服条件为 单轴拉伸的屈服强度单轴拉伸的屈服强度 等效塑性应变等效塑性应变 材料在初始屈服之后屈服强度的增加称为功硬化或者应变硬化材料在初始屈服之后屈服强度的增加称为功硬化或者应变硬化(对应于应变软化对应于应变软化)。硬化行为一般是塑性变形先期历史的函数。硬化行为一般是塑性变形先期历史的函数。屈服行为是各向同性硬化；拉伸和压缩的屈服强度总是相等。屈服行为是各向同性硬化；拉伸和压缩的屈服强度总是相等。5 一维塑性一维塑性 一维率无关塑性一维率无关塑性 一个特殊的模型，一个特殊的模型，塑性应变率写成为塑性应变率写成为 塑性模型称为关联的，否则，塑性流动是非关联的。塑性模型称为关联的，否则，塑性流动是非关联的。对于关联塑性，塑性流动是沿着屈服面的法线方向。对于关联塑性，塑性流动是沿着屈服面的法线方向。由此看出由此看出仅当当满足屈服条件足屈服条件时发生塑性生塑性变形。形。当塑性加载时，应力必须保持在屈服面上，当塑性加载时，应力必须保持在屈服面上，实现了实现了一致性条件一致性条件 这给出出塑性模量塑性模量 5 一维塑性一维塑性 一维率无关塑性一维率无关塑性 典型的硬化曲典型的硬化曲线，塑性模量塑性模量对应塑性加载和纯弹性加载或卸载，切线模量为对应塑性加载和纯弹性加载或卸载，切线模量为 塑性转换参数塑性转换参数加载卸载条件还可以写为加载卸载条件还可以写为 一致性条件的率形式一致性条件的率形式 应力状态位于塑性表面应力状态位于塑性表面 塑性率参数非负塑性率参数非负 对于塑性加载对于塑性加载 必须保持在屈服面上必须保持在屈服面上 其应力状态其应力状态对于弹性加载或者卸载对于弹性加载或者卸载 没有塑性流动没有塑性流动 因此因此 材料硬化描述材料硬化描述(a)Bauschinger效果效果(b)屈服面的平移和扩展屈服面的平移和扩展 在在循循环环加加载载中中，各各向向同同性性硬硬化化模模型型提提供供了了金金属属应应力力应应变变反反应应的的粗粗糙糙模模型型。图图a为为Bauschinger效效果果，在在拉拉伸伸初初始始屈屈服服之之后后的的压压缩缩屈屈服服强强度度降降低低。认认识识这这种种行行为为的的方方法法之之一一是是观观察察屈屈服服表表面面的的中中心心沿沿着着塑塑性性流流动动方方向向移移动动。图图b为为多多轴轴应应力力状状态态圆圆环环屈屈服服表表面面扩扩张张对对应应于于各各向向同性硬化（幂硬化）同性硬化（幂硬化），它的中心平移对应于，它的中心平移对应于运动硬化运动硬化。5 一维塑性一维塑性 混合硬化混合硬化 屈服面积改变，屈服中心不变，各向同性硬化；屈服面积改变，屈服中心不变，各向同性硬化；屈服面积不变，屈服中心平移，运动硬化。屈服面积不变，屈服中心平移，运动硬化。背应力的内部变量背应力的内部变量 Stress-strain curve under cyclic loadsCombined hardening model 混合硬化混合硬化 5 一维塑性一维塑性 屈服面积改变，屈服中心不变，各向同性硬化；屈服面积改变，屈服中心不变，各向同性硬化；屈服面积不变，屈服中心平移，运动硬化。屈服面积不变，屈服中心平移，运动硬化。5 一维塑性一维塑性 运动硬化运动硬化 塑性流动关系塑性流动关系 背应力背应力的内部变量的内部变量 屈服条件屈服条件一维率相关塑性一维率相关塑性 在率相关塑性中，材料的塑性反应取决于加载率在率相关塑性中，材料的塑性反应取决于加载率，一种方法是过应力模型，等效塑性应变率取决于超过多少屈服应力一种方法是过应力模型，等效塑性应变率取决于超过多少屈服应力 等效塑性应变率的一种交换形式等效塑性应变率的一种交换形式 粘度粘度 过应力过应力5 一维塑性一维塑性 应变软化应变软化 单调凸本构曲线不再成立。应变软化如何加载？单调凸本构曲线不再成立。应变软化如何加载？位移加载位移加载6 多轴塑性多轴塑性 Tresca屈服准则屈服准则Mises屈服准则屈服准则在有限元程序中一般应用哪种屈服准则？为什么？在有限元程序中一般应用哪种屈服准则？为什么？摩擦滑移屈服表面摩擦滑移屈服表面 6 多轴塑性多轴塑性 Mohr-Coulomb本构模型本构模型 滑移方向（塑性流动）是水平的（沿Q 的方向）而不是垂直屈服面。这是非关联塑性流动的例子。对于连续体和多轴应力应变状态的行为，M-C准则具有普适性。它应用于模拟土壤和岩石。M-C准则是基于这样的概念，即当任意面上的切应力和平均法准则是基于这样的概念，即当任意面上的切应力和平均法向应力达到临界组合时在材料中发生屈服向应力达到临界组合时在材料中发生屈服 c是内聚力，通过是内聚力，通过 定义内摩擦角定义内摩擦角 6 多轴塑性多轴塑性 Mohr-Coulomb屈服行为屈服行为Mohr-Coulomb屈服表面屈服表面Drucker-Prager屈服表面屈服表面 在在Mohr平面上的两条直线代表了方程式，它们是平面上的两条直线代表了方程式，它们是Mohr圆的包络圆的包络并称为并称为Mohr破坏或者失效包络。假设主应力破坏或者失效包络。假设主应力 应力状态应力状态屈服准则屈服准则6 多轴塑性多轴塑性 考虑考虑 的特殊情况并让的特殊情况并让，代表剪切屈服强度，代表剪切屈服强度，上式成为上式成为 即为即为Tresca准则。准则。在在Tresca和和M-C屈屈服服表表面面上上的的直直线线线线段段便便于于塑塑性性问问题题的的解解析析处处理理。然然而而，从从计计算算的的观观点点看看，夹夹角角使使得得本本构构方方程程难难以以建建立立(例例如如，计计算算屈屈服服面面的的法法线线)。通通过过改改进进von Mises屈屈服服准准则则结结合合压压力力的的影影响，响，Drucker-Prager屈服准则避免了与夹角有关的问题：屈服准则避免了与夹角有关的问题：这是一个光滑圆锥的方程，这是一个光滑圆锥的方程，为等效为等效Cauchy应力，选择常数有应力，选择常数有 D-P屈服表面通过了屈服表面通过了M-C屈服表面上的内部或者外部顶点（取加屈服表面上的内部或者外部顶点（取加号对应于内部顶点，而取减号对应于外部顶点）。号对应于内部顶点，而取减号对应于外部顶点）。9 应力更新算法应力更新算法本本构构方方程程率率形形式式的的积积分分算算法法称称为为应应力力更更新新算算法法（也也称称为为本本构构更更新新算算法），包括：法），包括：径向返回算法的一类图形返回算法，径向返回算法的一类图形返回算法，算法模量与基本应力更新方案一致的概念，算法模量与基本应力更新方案一致的概念，大变形问题的增量客观应力更新方案，大变形问题的增量客观应力更新方案，基于弹性响应的应力更新方案，即自动满足客观性的超弹性势能。基于弹性响应的应力更新方案，即自动满足客观性的超弹性势能。给出描述本构模型的某些其它连续介质力学观点，给出描述本构模型的某些其它连续介质力学观点，展示展示Eulerian，Lagrangian和两点拉伸的概念，和两点拉伸的概念，描述后拉、前推和描述后拉、前推和Lie导数的运算，导数的运算，材材料料框框架架客客观观性性，材材料料的的对对称称性性，以以本本构构行行为为的的张张量量表表示示讨讨论论了不变性的某些方面，了不变性的某些方面，讨讨论论由由于于热热力力学学第第二二定定律律和和某某些些附附加加的的稳稳定定性性必必要要条条件件对对材材料料行为的约束。行为的约束。9 应力更新算法应力更新算法 对对于于积积分分率率本本构构方方程程的的数数值值算算法法称称为为本本构构积积分分算算法法或或者者应应力力更更新新算法算法。对于率无关和率相关材料提供。对于率无关和率相关材料提供了了本构积分算法。本构积分算法。讨讨论论简简单单的的小小应应变变塑塑性性，将将小小应应变变算算法法扩扩展展至至大大变变形形，将将大大变变形形分析的积分算法保持在基于本构方程客观性的基础上。分析的积分算法保持在基于本构方程客观性的基础上。展示了关于大变形塑性的逐步客观积分算法。展示了关于大变形塑性的逐步客观积分算法。讨讨论论关关于于大大变变形形超超弹弹塑塑性性材材料料的的应应力力更更新新算算法法，回回避避对对应应力力率率方程的积分。方程的积分。描描述述了了与与本本构构积积分分算算法法相相关关的的计计算算模模量量，采采用用隐隐式式求求解解算算法法发发展展材料的切线刚度矩阵。材料的切线刚度矩阵。率无关塑性的图形返回算法率无关塑性的图形返回算法 9 应力更新算法应力更新算法小应变、率无关弹塑性的本构方程小应变、率无关弹塑性的本构方程 应力应变反应与变形率无关的一种材料称为率无关；否则为率相关。应力应变反应与变形率无关的一种材料称为率无关；否则为率相关。，Kuhn-Tucker条件，上面第一个条件表明塑性率参数是非条件，上面第一个条件表明塑性率参数是非负的，的，第二个条件表明第二个条件表明当塑性加载时，当塑性加载时，应力状力状态必必须位于或限制在塑性表面上，位于或限制在塑性表面上，最后条件也可以作最后条件也可以作为由已知一致性条件由已知一致性条件的率形式。的率形式。塑性流塑性流动方向方向经常特指常特指为，这里里称称为塑性流塑性流动势 屈服条件屈服条件 是是标量塑性流量塑性流动率，率，是塑性流动方向是塑性流动方向h 塑性模量塑性模量 q 内变量内变量）应力状力状态必必须保持在屈服面保持在屈服面因此因此。对于于弹性加性加载或者卸或者卸载，没有塑性流，没有塑性流动。对于塑性加于塑性加载（率无关塑性的图形返回算法率无关塑性的图形返回算法 9 应力更新算法应力更新算法上，上，在时刻在时刻n 给出一组给出一组 和和应变增量增量 本构本构积分算法的目的是分算法的目的是计算算并并满足加卸足加卸载条件条件 在在时刻的刻的应力力给出出为 求解的一致性条件给出 设想能够应用这个塑性参数值以提供更新的应力率、塑性应变设想能够应用这个塑性参数值以提供更新的应力率、塑性应变率和内变量率，并且写出简单的率和内变量率，并且写出简单的向前向前Euler积分公式算法积分公式算法率无关塑性的图形返回算法率无关塑性的图形返回算法 9 应力更新算法应力更新算法但在下一步，但在下一步，这些些应力和内力和内变量的更新量的更新值并不并不满足屈服条件，所以足屈服条件，所以 由由于于解解答答从从屈屈服服表表面面漂漂移移，常常常常导导致致不不精精确确的的结结果果，因因此此不不受受人人青青睐睐。公公式式也也称称为切切线模模量量更更新新算算法法，形形成成了了计算算率率无无关关塑塑性早期工作的基性早期工作的基础。率无关塑性的图形返回算法率无关塑性的图形返回算法 9 应力更新算法应力更新算法 这导致考致考虑另外一些方法另外一些方法进行率本构方程的行率本构方程的积分，目的之一是分，目的之一是强强化在化在时间步步结束束时的的一致性一致性，例如，例如，为避避免免离离开开屈屈服服面面的的漂漂移移。有有许多多不不同同的的积分分本本构构算算法法，这里里主主要要关关注注一一类方方法法返返回回图形形算算法法，它它是是强强健健和和精精确确的的，被被广广泛泛应用。著名的用。著名的von Mises塑性塑性径向返回方法径向返回方法是返回是返回图形算法的特例。形算法的特例。返回图形算法包括：返回图形算法包括：一个初始的弹性预测步，包含（在应力空间）对屈服表面的偏离，一个初始的弹性预测步，包含（在应力空间）对屈服表面的偏离，以及塑性调整步使应力返回到更新后的屈服表面。以及塑性调整步使应力返回到更新后的屈服表面。方法的两个组成部分是：方法的两个组成部分是：一个积分算法，它将一组本构方程转
展开阅读全文