ansys级非线性分析蠕变课件

蠕变四章蠕变四章

隐式和显式蠕变

本章综述前一章探讨了率无关塑性,本章将讨论ANSYS中蠕变行为的分析方法。尽管从材料的观点看蠕变和粘塑性是相同的,但本构模型的使用不同。因此,率相关的塑性主题被分为两部分,这部分属于蠕变。本章将介绍ANSYS中可用的大量的隐式和显式蠕变法则。主要讨论金属的蠕变。然而,讨论的各种观点也适用于塑料或混凝土等其它材料的蠕变。ANSYS有隐式和显式两种蠕变分析过程。首先讨论蠕变的一般知识，然后是进行隐式或显式蠕变分析的细节。

隐式和显式蠕变

…本章综述本章包括下列主题：A.蠕变现象背景B.术语定义C.一般蠕变方程D.隐式蠕变过程E.显式蠕变过程F.ANSYS蠕变模型的求解过程G.隐式蠕变与显式蠕变的比较

隐式和显式蠕变

A.蠕变背景晶体材料中,如金属,蠕变机理与空隙的扩散流动和位错运动有关。空隙是点缺陷,倾向于形成与所施加应力方向垂直(而不是平行)的晶界。空隙由高集中区向低集中区运动。在低应力状态下发生扩散流动，但通常需要高温条件。晶粒的位错是线缺陷.位错运动(攀升、滑动、偏移)在高应力状态下被激活,尽管在中温时也可能发生位错运动。有时晶界滑移被认为是一种独立的导致蠕变变形的机理。

隐式和显式蠕变

...蠕变背景尽管对材料科学的详细论述超越了本书的范围，也足以说明前面提及的物理机理导致蠕变。蠕变变形与应力、应变、时间和温度的相关性一般用与下式相似的形式来模拟：函数f1-4

与选择的蠕变法则有关。通常由不同应变速率和温度条件下的各种拉伸实验确定相关蠕变常数。假设各向同性行为,vonMises方程用于计算有效应力,在蠕变应变率方程中使用等效应变(与率无关塑性相似)。

隐式和显式蠕变

...蠕变背景当计算弹性、塑性和蠕变应变时,ANSYS使用附加的应变分量:

计算塑性应变(流动法则)类似于前一章的描述。根据蠕变应变率方程计算蠕变应变，它的详细形式将在后面讨论。弹性、蠕变和塑性应变都基于(当前的)应力状态，但它们是独立计算的(彼此互不依赖)。注意在利用隐式蠕变与显式蠕变进行计算时有差别。附加分量应力-应变

隐式和显式蠕变

...蠕变背景与塑性类似,蠕变是一种基于偏量行为的不可逆(非弹性)应变。在蠕变流动条件下,假设材料是不可压缩的。另一方面,与率无关塑性不同，蠕变在发生非弹性应变时没有屈服面。因此,蠕变不需要高应力值来发生更多的蠕变应变。假设在所有非零应力值时都会发生蠕变应变。前面提到,从材料的角度看,蠕变和粘塑性是相同的。工程应用中,通常蠕变用于描述低应变速率的热激活过程。率无关塑性和隐式蠕变应变以弱耦合方式处理。相反,ANSYS中粘塑性本构模型用于描述高应变速率的应用(例如,冲击载荷)。非弹性应变以强耦合方式处理。

隐式和显式蠕变

...蠕变分析实例 焊球蠕变分析的例子(热循环)。BretZahn,ChipPAC授权的ANSYS模型(/)单元185(B-Bar),双曲正弦隐式蠕变模型

隐式和显式蠕变

B.术语的定义蠕变的三个阶段:在常载荷条件下,蠕变的单轴应变与时间的关系如下图所示。第一阶段,应变率随时间减小，该阶段在较短时间内完成。第二阶段具有常应变率。第三阶段,应变率迅速增加直到破坏(断裂)。t第二阶段第三阶段第一阶段断裂

隐式和显式蠕变

...术语的定义蠕变的三个阶段(续)蠕变应变率可能是应力、应变、温度、和/或时间的函数。对于工程分析,通常最关心蠕变的第一和第二阶段，第三阶段蠕变通常与开始破坏有关(颈缩、损坏),且时间较短，因此,ANSYS中不模拟第三阶段。第一阶段蠕变的应变率通常远大于第二阶段蠕变。然而,应变率在第一阶段逐渐降低而在第二阶段几乎为常值(对于前面提及的常应力、恒温下单轴试验情况而言)，而且第一阶段蠕变时间比第二阶段短。

隐式和显式蠕变

...术语的定义蠕变在施加恒应力条件下,蠕变应变增加。应力松弛在施加恒应变条件下,应力降低。tets

隐式和显式蠕变

...术语的定义时间强化假设蠕变应变率仅与从蠕变过程开始的时间有关，也就是说,该曲线上/下移动，当应力从s1

到s2变化时,计算A到B点的不同蠕变速率。应变强化假设蠕变速率仅与材料现有的应变有关，也就是说,曲线左/右移动.当应力从s1

到s2变化时,计算A到B点的不同蠕变应变率。tes1s2ABtes1s2AB

隐式和显式蠕变

...术语的定义显式蠕变显式蠕变是指应用Euler向前法进行蠕变应变演化的计算。在每个时间步使用的蠕变应变率与该时间步开始时的速率一致，并假设在整个时间步Dt内为常量，因此，需要很小的时间步以使减小误差。对于有塑性的显式蠕变,首先进行塑性修正,然后进行蠕变修正。两种修正在不同应力值处进行，因此精度较差。

隐式和显式蠕变

...术语的定义隐式蠕变隐式蠕变应用了Euler向后积分法求解蠕变应变。该方法在数值上无条件稳定，这意味着不必象显式蠕变方法那样,使用小的时间步,所以总体上会更快。对于隐式蠕变加上率无关塑性,塑性修正和蠕变修正同时进行,而不是分别进行。因此,隐式蠕变一般比显式蠕变更精确,但它仍与时间步大小有关，必须使用足够小的时间步来精确捕捉路径相关行为。基于上述原因,隐式蠕变是ANSYS推荐的方法(高效、精确)。后面将详细讨论两种蠕变过程。

隐式和显式蠕变

C.一般蠕变方程如前所述,蠕变方程通常是一种速率形式,类似于下式:然而,被分析的材料类型决定了具体蠕变方程的选择。现在讨论一些一般的特性,具体模型包括在隐式蠕变和显式蠕变两部分中。单元手册第2.5章节也包括了隐式和显式蠕变方程。

隐式和显式蠕变

...一般蠕变方程温度相关性蠕变效应被热激活,它的温度相关性通常由Arrhenius法则来表示:

式中Q为激活能,R为普适气体常数,T为绝对温度。应力相关性蠕变应变通常也与应力有关,尤其是位错蠕变。Norton法则为:

对上述幂定律的常见修正如下:

隐式和显式蠕变

...一般蠕变方程通常第一阶段蠕变显示时间或应变强化。时间强化包含在一个时间相关项中:应变强化包含在一个应变相关项中:由可用的材料数据来决定使用哪一项(应变或时间强化)。第二阶段不具有时间或应变强化，第二阶段的蠕变应变率通常是常数。

隐式和显式蠕变

...一般蠕变方程ANSYS中可用的蠕变法则汇总如下:

如前所述,只要可能,推荐采用隐式蠕变。

隐式蠕变D节隐式蠕变D节

隐式和显式蠕变

D.隐式蠕变过程本节讨论进行隐式蠕变分析的过程。如前所述,因为隐式蠕变比显式蠕变更稳定、精确和高效(更快),所以是首选的方法。这是由于该公式中使用了隐式积分法(Euler向后法)。有些情况下采用的蠕变法则或单元类型需要使用显式蠕变方法。

隐式和显式蠕变

...支持的单元类型隐式蠕变材料支持的单元类型:“核心”单元:PLANE42,SOLID45,PLANE82,SOLID92和SOLID95“18x”系列单元:LINK180,SHELL181,PLANE182,PLANE183,SOLID185,SOLID186,SOLID187,BEAM188和BEAM189隐式蠕变分析中,推荐选择“18x”系列单元。因为在“18x”系列单元中有大量的单元技术可用(参考第二章),这些单元提供了更强的灵活性和能力。这些公式包括B-Bar,URI,增强应变及混合U-P。

“18x”系列单元比核心单元支持更多的本构模型,包括超弹性。

隐式和显式蠕变

...支持的塑性模型前已述及,蠕变与率无关塑性解耦。隐式蠕变允许与下述率无关塑性模型相结合(见第三章):有蠕变的BISO,MISO和NLISO(等向强化蠕变)有蠕变的BKIN(随动强化蠕变)蠕变与HILL(各向异性蠕变)蠕变与BISO,MISO和有HILL的NLISO(等向强化的各向异性蠕变)蠕变与有HILL的BKIN(随动强化的各向异性蠕变)

隐式和显式蠕变

...定义隐式蠕变模型可以用命令或通过GUI定义隐式蠕变模型(在下面的幻灯片中讨论)。对于隐式蠕变,以两种方式定义温度相关性:由TBDATA(或材料GUI)定义温度相关常数很多蠕变方程包括前面提及的Arrhenius方程由用户确定使用任一个或两个方法来包含温度相关性。

隐式和显式蠕变

...定义隐式蠕变模型通过命令定义隐式蠕变时,使用TB,CREEP,具体蠕变模型由TBOPT定义。TB,CREEP,mat,ntemp,npts,TBOPTTBTEMP定义与温度相关的常数TBDATA定义实际常数

对下面的例子,TBOPT=2是指采用时间-强化蠕变方程.使用TBTEMP命令指定温度相关常数,与该方程相关的四个常数由TBDATA命令指定为自变量。

TB,CREEP,1,1,4,2

TBTEMP,100

TBDATA,1,C1,C2,C3,C4

隐式和显式蠕变

...定义隐式蠕变模型所有隐式蠕变模型可以在如下的材料GUI中选择:Structural

>Nonlinear

>Inelastic

>RateDependent

>Creep确保首先定义必需的线弹性材料属性(EX和PRXY)。

隐式和显式蠕变

...定义隐式蠕变模型当选择合适的隐式蠕变模型后,出现一个独立对话框显示所需的输入项。下面例子中,已经定义了第一阶段蠕变方程,提示用户输入四个蠕变常数。也可以输入温度相关的常数。

隐式与显式蠕变

...可用的隐式蠕变模型下表显示可用的隐式蠕变法则，方程在下面幻灯片中说明。

隐式和显式蠕变

...可用的隐式蠕变模型1)

应变强化

TBOPT=1

第一阶段蠕变2)

时间强化

TBOPT=2

第一阶段蠕变3) 广义指数

TBOPT=3

第一阶段蠕变4) 广义Graham

TBOPT=4

第一阶段蠕变

隐式和显式蠕变

...可用的隐式蠕变模型5) 广义Blackburn

TBOPT=5

第一阶段蠕变6) 修正的时间强化

TBOPT=6

第一阶段蠕变7) 修正的应变强化

TBOPT=7

第一阶段蠕变

隐式和显式蠕变

...可用的隐式蠕变模型8) 广义Garofalo

TBOPT=8

第二阶段蠕变9) 指数形式

TBOPT=9

第二阶段蠕变10) Norton

TBOPT=10

第二阶段蠕变

隐式和显式蠕变

...可用的隐式蠕变模型11) 时间强化

TBOPT=11

第一阶段+第二阶段12) 有理多项式

TBOPT=12

第一阶段+第二阶段13)

广义时间强化

TBOPT=13

第一阶段蠕变

隐式和显式蠕变

...练习请参考附加练习题:练习6:应力松弛

蠕变四章蠕变四章

隐式和显式蠕变

本章综述前一章探讨了率无关塑性,本章将讨论ANSYS中蠕变行为的分析方法。尽管从材料的观点看蠕变和粘塑性是相同的,但本构模型的使用不同。因此,率相关的塑性主题被分为两部分,这部分属于蠕变。本章将介绍ANSYS中可用的大量的隐式和显式蠕变法则。主要讨论金属的蠕变。然而,讨论的各种观点也适用于塑料或混凝土等其它材料的蠕变。ANSYS有隐式和显式两种蠕变分析过程。首先讨论蠕变的一般知识，然后是进行隐式或显式蠕变分析的细节。

隐式和显式蠕变

…本章综述本章包括下列主题：A.蠕变现象背景B.术语定义C.一般蠕变方程D.隐式蠕变过程E.显式蠕变过程F.ANSYS蠕变模型的求解过程G.隐式蠕变与显式蠕变的比较

隐式和显式蠕变

A.蠕变背景晶体材料中,如金属,蠕变机理与空隙的扩散流动和位错运动有关。空隙是点缺陷,倾向于形成与所施加应力方向垂直(而不是平行)的晶界。空隙由高集中区向低集中区运动。在低应力状态下发生扩散流动，但通常需要高温条件。晶粒的位错是线缺陷.位错运动(攀升、滑动、偏移)在高应力状态下被激活,尽管在中温时也可能发生位错运动。有时晶界滑移被认为是一种独立的导致蠕变变形的机理。

隐式和显式蠕变

...蠕变背景尽管对材料科学的详细论述超越了本书的范围，也足以说明前面提及的物理机理导致蠕变。蠕变变形与应力、应变、时间和温度的相关性一般用与下式相似的形式来模拟：函数f1-4

与选择的蠕变法则有关。通常由不同应变速率和温度条件下的各种拉伸实验确定相关蠕变常数。假设各向同性行为,vonMises方程用于计算有效应力,在蠕变应变率方程中使用等效应变(与率无关塑性相似)。

隐式和显式蠕变

...蠕变背景当计算弹性、塑性和蠕变应变时,ANSYS使用附加的应变分量:

计算塑性应变(流动法则)类似于前一章的描述。根据蠕变应变率方程计算蠕变应变，它的详细形式将在后面讨论。弹性、蠕变和塑性应变都基于(当前的)应力状态，但它们是独立计算的(彼此互不依赖)。注意在利用隐式蠕变与显式蠕变进行计算时有差别。附加分量应力-应变

隐式和显式蠕变

...蠕变背景与塑性类似,蠕变是一种基于偏量行为的不可逆(非弹性)应变。在蠕变流动条件下,假设材料是不可压缩的。另一方面,与率无关塑性不同，蠕变在发生非弹性应变时没有屈服面。因此,蠕变不需要高应力值来发生更多的蠕变应变。假设在所有非零应力值时都会发生蠕变应变。前面提到,从材料的角度看,蠕变和粘塑性是相同的。工程应用中,通常蠕变用于描述低应变速率的热激活过程。率无关塑性和隐式蠕变应变以弱耦合方式处理。相反,ANSYS中粘塑性本构模型用于描述高应变速率的应用(例如,冲击载荷)。非弹性应变以强耦合方式处理。

隐式和显式蠕变

...蠕变分析实例 焊球蠕变分析的例子(热循环)。BretZahn,ChipPAC授权的ANSYS模型(/)单元185(B-Bar),双曲正弦隐式蠕变模型

隐式和显式蠕变

B.术语的定义蠕变的三个阶段:在常载荷条件下,蠕变的单轴应变与时间的关系如下图所示。第一阶段,应变率随时间减小，该阶段在较短时间内完成。第二阶段具有常应变率。第三阶段,应变率迅速增加直到破坏(断裂)。t第二阶段第三阶段第一阶段断裂

隐式和显式蠕变

...术语的定义蠕变的三个阶段(续)蠕变应变率可能是应力、应变、温度、和/或时间的函数。对于工程分析,通常最关心蠕变的第一和第二阶段，第三阶段蠕变通常与开始破坏有关(颈缩、损坏),且时间较短，因此,ANSYS中不模拟第三阶段。第一阶段蠕变的应变率通常远大于第二阶段蠕变。然而,应变率在第一阶段逐渐降低而在第二阶段几乎为常值(对于前面提及的常应力、恒温下单轴试验情况而言)，而且第一阶段蠕变时间比第二阶段短。

隐式和显式蠕变

...术语的定义蠕变在施加恒应力条件下,蠕变应变增加。应力松弛在施加恒应变条件下,应力降低。tets

隐式和显式蠕变

...术语的定义时间强化假设蠕变应变率仅与从蠕变过程开始的时间有关，也就是说,该曲线上/下移动，当应力从s1

到s2变化时,计算A到B点的不同蠕变速率。应变强化假设蠕变速率仅与材料现有的应变有关，也就是说,曲线左/右移动.当应力从s1

到s2变化时,计算A到B点的不同蠕变应变率。tes1s2ABtes1s2AB

隐式和显式蠕变

...术语的定义显式蠕变显式蠕变是指应用Euler向前法进行蠕变应变演化的计算。在每个时间步使用的蠕变应变率与该时间步开始时的速率一致，并假设在整个时间步Dt内为常量，因此，需要很小的时间步以使减小误差。对于有塑性的显式蠕变,首先进行塑性修正,然后进行蠕变修正。两种修正在不同应力值处进行，因此精度较差。

隐式和显式蠕变

...术语的定义隐式蠕变隐式蠕变应用了Euler向后积分法求解蠕变应变。该方法在数值上无条件稳定，这意味着不必象显式蠕变方法那样,使用小的时间步,所以总体上会更快。对于隐式蠕变加上率无关塑性,塑性修正和蠕变修正同时进行,而不是分别进行。因此,隐式蠕变一般比显式蠕变更精确,但它仍与时间步大小有关，必须使用足够小的时间步来精确捕捉路径相关行为。基于上述原因,隐式蠕变是ANSYS推荐的方法(高效、精确)。后面将详细讨论两种蠕变过程。

隐式和显式蠕变

C.一般蠕变方程如前所述,蠕变方程通常是一种速率形式,类似于下式:然而,被分析的材料类型决定了具体蠕变方程的选择。现在讨论一些一般的特性,具体模型包括在隐式蠕变和显式蠕变两部分中。单元手册第2.5章节也包括了隐式和显式蠕变方程。

隐式和显式蠕变

...一般蠕变方程温度相关性蠕变效应被热激活,它的温度相关性通常由Arrhenius法则来表示:

式中Q为激活能,R为普适气体常数,T为绝对温度。应力相关性蠕变应变通常也与应力有关,尤其是位错蠕变。Norton法则为:

对上述幂定律的常见修正如下:

隐式和显式蠕变

...一般蠕变方程通常第一阶段蠕变显示时间或应变强化。时间强化包含在一个时间相关项中:应变强化包含在一个应变相关项中:由可用的材料数据来决定使用哪一项(应变或时间强化)。第二阶段不具有时间或应变强化，第二阶段的蠕变应变率通常是常数。

隐式和显式蠕变

...一般蠕变方程ANSYS中可用的蠕变法则汇总如下:

如前所述,只要可能,推荐采用隐式蠕变。

隐式蠕变D节隐式蠕变D节

隐式和显式蠕变

D.隐式蠕变过程本节讨论进行隐式蠕变分析的过程。如前所述,因为隐式蠕变比显式蠕变更稳定、精确和高效(更快),所以是首选的方法。这是由于该公式中使用了隐式积分法(Euler向后法)。有些情况下采用的蠕变法则或单元类型需要使用显式蠕变方法。

隐式和显式蠕变

...支持的单元类型隐式蠕变材料支持的单元类型:“核心”单元:PLANE42,SOLID45,PLANE82,SOLID92和SOLID95“18x”系列单元:LINK180,SHELL181,PLANE182,PLANE183,SOLID185,SOLID186,SOLID187,BEAM188和BEAM189隐式蠕变分析中,推荐选择“18x”系列单元。因为在“18x”系列单元中有大量的单元技术可用(参考第二章),这些单元提供了更强的灵活性和能力。这些公式包括B-Bar,URI,增强应变及混合U-P。

“18x”系列单元比核心单元支持更多的本构模型,包括超弹性。

隐式和显式蠕变

...支持的塑性模型前已述及,蠕变与率无关塑性解耦。隐式蠕变允许与下述率无关塑性模型相结合(见第三章):有蠕变的BISO,MISO和NLISO(等向强化蠕变)有蠕变的BKIN(随动强化蠕变)蠕变与HILL(各向异性蠕变)蠕变与BISO,MISO和有HILL的NLISO(等向强化的各向异性蠕变)蠕变与有HILL的BKIN(随动强化的各向异性蠕变)

隐式和显式蠕变

...定义隐式蠕变模型可以用命令或通过GUI定义隐式蠕变模型(在下面的幻灯片中讨论)。对于隐式蠕变,以两种方式定义温度相关性:由TBDATA(或材料GUI)定义温度相关常数很多蠕变方程包括前面提及的Arrhenius方程由用户确定使用任一个或两个方法来包含温度相关性。

隐式和显式蠕变

...定义隐式蠕变模型通过命令定义隐式蠕变时,使用TB,CREEP,具体蠕变模型由TBOPT定义。TB,CREEP,mat,