ANSYS结构非线性分析指南 - 第四章(8)

ANSYS结构非线性分析指南

*dim,xval,table,100

strn(1,1)=-0.45,-0.4,-0.35,-0.3,-0.25,-0.2,-0.15,-0.1,-0.05 strss(1,1)=-256,-128,-64,-32,-16,-8,-4,-2,-1

strn(10,1)=0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9 strss(10,1)=0.0,1,1.5,2.0,2.9,3.6,5,7.5,9.7,17

tb,mooney,,,,1

*mooney,strn(1,1),strss(1,1),,const(1),calc(1,1),sortsn(1,1),sortss(1,1)

*eval,1,1,const(1),-0.2,0,xval(1),ecalc(1) *vplot,xval(1),ecalc(1) fini

/solu alls

nsel,s,loc,y,0 d,all,all,0 alls

nsel,s,loc,x,0 d,all,ux,0 d,all,uz,0 alls

nsel,s,loc,y,h3 d,all,ux,0 d,all,uz,0 d,all,uy,-d alls

antype,static nlgeom,on nropt,,,on outpr,all,all outres,all,all autots,on time,1

deltim,0.03,0.01,0.3 cnvtol,f,,0.02,2 lnsrch,on pred,on alls solve fini

4-36

ANSYS结构非线性分析指南

4.4 蠕变分析

4.4.1 蠕变理论 4.4.1.1 定义

蠕变是率相关材料非线性，即在常荷载作用下，材料连续变形的特性。相反如果位移固定，反力或应力将随时间而变小，这种特性有时也称为应力松驰，见图4-18a。

图4-18 应力松弛和蠕变

蠕变的三个阶段如图4-18b所示。在初始蠕变阶段，应变率随时间而减小，这个阶段一般发生在一个相当短的时期。在第二期蠕变阶段，有一个常应变率，所以应变以常速率发展，在第三期蠕变阶段，应变率迅速增加直到材料失效。

由于第三期蠕变阶段所经历的时间很短，材料将失效，所以通常情况下，我们感兴趣的是初始蠕变和第二期蠕变。ANSYS程序中的蠕变行为用来模拟初始蠕变和第二期蠕变。蠕变系数可以是应力、应变、温度、时间或其它变量的函数。

在高温应力分析中(如核反应堆等)，蠕变分析非常重要。例如，假设在核反应堆中施加了预荷载，以保证与相邻部件保持接触而不松开。在高温下过了一段时间后，预荷载将降低(应力松驰)，可能使接触部件松开。对于一些材料如预应力砼，蠕变也可能十分重要。最重要的是要记住，蠕变是永久变形。 4.4.1.2 理论介绍

蠕变方程：我们通过一个方程来模拟蠕变行为，此方程描述了在实验中观测到的主要特征（特别是在一维的拉伸实验中）。这个方程以蠕应变率的方式表示出来，其形式如下：

??cr?A?B?ctD

上式中，A、B、C、D是从实验中得到的材料常数，常数本身也可能是应力，应变，时间或温度的函数，这种形式的方程被称为状态方程。

上式中，当常数D为负值时，蠕应变率随时间下降，材料处于初始蠕变阶段，当D为0时，蠕应变率为常值，材料处于第二期蠕变阶段。

对于2－D或3－D应力状态，使用VON Mises方程计算蠕应变率方程中所使用的标量等效应力和等效应变。

对蠕变方程积分时，我们使用经过修改的总应变，其表达式为:

????????????????????

plthcrnnnnn?1 经过修改的等效总应变为:

4-37

ANSYS结构非线性分析指南

?et?1???2??2??23?23?23?2?

(???)?(???)?(???)?(?)?(?)?(?)?xyyzzxxyyzzx?2222??其等效应力由下式算出：

?e?2G?et

其中：G=剪切模量＝E2(1??) 等效蠕应变增量??cr由程序给出的某一种公式进行计算，一般为正值，如果在数据表中

C11?1，则使用的是衰减的蠕应变率而不是常蠕变率，但这个选项一般不被推荐，因为在初始蠕变

所产生的应力为主的情况下，它可能会严重的低估蠕变值。如果C11?1，程序使用修正的等效蠕应变增量来代替蠕应变增量。

??mcr??et(1? 其中：e=2.718281828（自然对数的底数）

1) AeA???cr?et

下面是计算积分点的蠕应变率与弹性应变比率的公式：

CS???cr?et

将本次迭代的所有单元的所有积分点的的最大值记为CS，并且作为“CREEPRATIO”输出。 计算出等效蠕应变增量后，可将它转换成分量的形式，假设Nc是某个特定单元类型的应变分量的个数。

如果Nc?1则有：

??x 如果Nc?4，则有： ??xcrcr?x????()

?etcr?????cr(2?x??y??z)? ?et2(1??)?????cr(2?y??z??x)? ?et2(1??) ??y ??zcrcr????xcr???ycr

4-38

ANSYS结构非线性分析指南

??xycr???cr?et3?xy?

2(1??) 上式中，前三个为正应变分量，第四个是剪应变分量。 如果Nc?6，前四个分量与上式相同，另两个剪应变分量为： ??yzcr???cr?et?et?3?yz?

2(1??)3?xz?

2(1??) ??xzcr???cr 接下来，可以按下式来计算弹性应变和总的蠕应变（以X方向的分量为例）： (?x)n?(?x)n???x (?x)n?(?x)n?1???x 为了从标量??crelcrcrcrcr来计算分量?x?cr,?y?cr,?z?cr,程序使用相关流动准则：Prandtl-Reuss方程，与

塑性应变相同，蠕应变只有偏差分量（剪分量），没有由于蠕变引起的体积应变。

为了考虑应力随时间的变化，使用两种强化准则，时间强化和应变强化。我们以一简单拉伸试验来说明：刚开始时，杆被加载到应力为?1，在时间t它被卸载到应力为?2。

(a) 时间强化 (b) 应变强化

图4-19 典型的单轴蠕变曲线

时间强化假定蠕应变率仅仅依赖于蠕应变过程开始的时间。当应力从?1变到?2时，材料的蠕变率由点A表示（相当于曲线向上移动）。

应变强化假定蠕应变率仅仅依赖于材料中的应变，当应力从?1变到?2时，材料的蠕变率由点B表示（相当于曲线左移）。大多数实验数据与应变强化准则吻合得更好。

4-39

ANSYS结构非线性分析指南

4.4.2 求解算法

ANSYS使用隐式和显式积分二种方法来进行蠕变分析，均可应用于静态和瞬态分析。隐式蠕变分析方法更强大、更快、更精确，一般推荐使用隐式蠕变分析。它可以处理温度相关蠕变常数，同时模拟蠕变与等向强化塑性模型。

对于需要很小时间步的情况，显式蠕变分析就非常有用。蠕变常数不能有温度相关性，而与其他塑性材料模型的耦合只能应用迭加法。

注意--蠕变分析中的“隐式”和“显式”，与“显式动力分析”或“显式单元”没有任何关系。

隐式蠕变分析方法支持下列单元：PLANE42,SOLID45,PLANE82,SOLID92, SOLID95,LINK180，SHELL181，PLANE182，PLANE183，SOLID185，SOLID186，SOLID187，BEAM188 和 BEAM189。

显式蠕变分析方法支持下列单元：LINK1,PLANE2,LINK8,PIPE20,BEAM23, BEAM24,PLANE42，SHELL43，SOLID45，SHELL51，PIPE60，SOLID62，SOLID65，PLANE82，SOLID92 和 SOLID95。 蠕变应变率可以是应力、应变、温度、电子流水平的函数。蠕变应变率方程已按初始蠕变、第二期蠕变和辐射引起的蠕变在ANSYS中建立。参见《ANSYS Elements Reference》中关于这些蠕变方程的讨论和输入方法。有一些方程需要特殊的单位。特别是，对于显式蠕变选项，蠕变方程中的温度应当基于绝对温度。 4.4.2.1 隐式蠕变方法

隐式蠕变方 …… 此处隐藏：2210字，全部文档内容请下载后查看。喜欢就下载吧 ……