ABAQUS反复加载卸载应力应变计算说明1、塑性模型采用ArmstrongandFrederickmodel(AF模型)[1]屈服准则(Mises屈服)3~~~~屈服函数为f=2(~_<~)：(~_(~)-Y2式中：S为应力偏量，为总背应力，Y为屈服极限流动准则(theassociatedflowrule)~p=迟d<~式中：~P为塑性应变对时间的微分，血为待定量[2],&为应力张量(3)硬化准则<1>除)=3h(i)~p-匚(i)(~(i)p式中：pi)为背应力分量对时间的微分；h(i)，:(i)为材料常数为已知量，P为等效塑性应变对时间的微分。<2>Y=Y+兰r(i)(1一e-Q⑴p)0i=1式中：Y为对应于等效塑性应变p的屈服极限，Y为初始屈服极限为已知量，r(i)为材料常0数为已知量。2、反复加载卸载应力应变计算过程说明假设在受载前，物体的初始应力、应变以及背应力均为零加载过程计算<1>外力不足以使得物体中的任何一点的Mises应力值大于屈服极限此时：塑性应变、背应力均保持为零，屈服极限保持不变。应力由&=De:~e计算，总应变值等于弹性应变。<2>外力使得物体中的任何一点的Mises应力值大于屈服极限为了说明ABAQUS是如何确定应变增量A~，有必要对ABAQUS求解材料非线性问n+1题进行简单介绍[3]。ABAQUS首先将载荷分为若干个微小增量，如图1所示。当结构收到一个微小增量AP时，ABAQUS用与初始结构位移相对应的初始刚度矩阵%和载荷增量AP计算出结构在这一增量后的位移修正c、修正后的位移值U和相应的新的刚度矩阵K。ABAQUS用新的aaa刚度矩阵计算结构的内力I，载荷P和I的差值为迭代的残余力R。如果R在模型内的aaaa每一个自由度上的值都为零或小于一个给定的容差，如图1所示的a点，则结构处于平衡状态即ABAQUS计算到的内外力是平衡的。若假设整个物体只有一个单元，则位移修正c就为应变增量A~，刚度矩阵就是an+1UMAT程序中的雅可比矩阵。根据塑性变形时的雅可比矩阵计算公式：Dp=De—4G2L-1:/其中De为弹性状态下的雅可比矩阵，可知Dp较De小。UMAT计d算出的应力便是结构的内力I。a在UMAT程序中通过计算：3~~~~Ftri=(~tri—(~)：(~tri—(~)—Y2TOC\o"1-5"\h\zn+12n+1nn+1n式中：C~tri=C~+De:Af~n+1nn+1并判断若Ft”是否大于零来决定节点是否达到屈服。若Ftri小于零则节点未达到屈服,n+1n+1若Ftri大于等于零节点达到屈服，需要进行塑性迭代与重新计算雅可比矩阵。n+1图1增量法迭代原理ABAQUS处理塑性加载过程问题可以描述为：首先读入上一状态的应力、应变、背应力及屈服强度，使用下标n表示（对应时间为t）。再根据上一步的刚度矩阵K和载荷增nn量AP计算出给定的应变增量A~，然后在UMAT程序中求出t（t二t+At）时n+1n+1n+1nn+1刻满足所用塑性模型给出的屈服准则、流动准则、硬化准则下的应力、应变、背应力及屈服强度。UMAT程序利用已知的分，分⑴，~，~p，P，Y和给定的A~首先求出可满足nnnnnnn+1所用塑性模型给出的屈服准则、流动准则、硬化准则的Ap，再使用公式依次求出其余量。n+1等效塑性应变：P=P+Apn+1nn+1屈服极限：Y二Y+n+10i=1n+1塑性流动：nn+1(~stri—艺e(i)o(i))n+1n+1nIi1.(Stri—艺e(i)((i)):(S艺eIn+1n+1nn+1n+1n'i=1i=1式中：e(i)n+11—(i)Apn+1〜厅塑性应变增量：A〜p=Apnn+1*2n+1n+1~~2~背应力：O(i)=e(i)(O(i)+h(i)AP)n+1n+1n3n+1+ICJtrimC2应力：&=(~(i)+：Ynn+1n+1\l3n+1n+1式中：Ctri=](Ctri+Ctri+Ctri)m3112233再计算出：E~，E~Pnn加载过程按照此方式进行计算同样已知：C~，o~(i)，E~，E~P，nnnnnn当外力卸载时P,Y并且给定：丈。这里考虑卸载过程的第n+1一步迭代。由于外力开始卸载，计算得到的(~tri—()：(Stri—()<(S—()：(S—()。n+1nn+1nnnnn由于上一步迭代材料处于屈服状态，所以在UMAT程序运行结束时有下面关系成立：3~~~~3~~~~F=—(s—oc):(s—()—Y2=0所以此时的F=—(stri—OL)：(stri—oc)—Y2<02nnnnnn+12n+1nn+1nn根据程序判断不进行塑性计算，只计算弹性雅可比矩阵即按照弹性方式进行卸载。各个参数有关系：o(i)=o(i)，EP=EP，P=P，Y=Yn+1nn+1nn+1nn+1n弹性应变：Ee=Ee+Aen+1nn+1总应变：E=Ee+EPnn+1n+1按照此方式，直到卸载完成。这里不讨论反向加载。再次加载时由于（~⑴，~p,Y在卸载过程中不变化，所以在第二次加载过程中，读入的&,&⑴，nnnnn~~p，p，Y中（~⑴，~p一定是不为零的，并且Y>Y。对于&是否为零，与残余nnnnnnn0n应变是否引起了残余应力有关。对于一个单元残余塑性应变是无法引起残余应力的；对于多个单元构成的物体，当卸载完成后，一些单元的残余塑性应变使得它们与周围的单元之间产生作用力从而产生残余应力。此时读入的是不为零的。这是与第一次加载初始应力状态n有区别的。根据UMAT弹性试应力的计算公式&衍=<~+De:A~，若物体中含有残余TOC\o"1-5"\h\zn+1nn+1应力则在计算过程中直接代入程序计算。~3~~~~由给定的应变增量As计算出下一时刻的Ftrl=（stri-（~）：（s加）—Y2。若n+1n+12n+1nn+1nFtri小于零则节点未达到屈服按照中的弹性加载过程进行计算；若Ftri大于等于零按照中n+1n+1的塑性加载过程进行计算。分析若要此模型出现图2所示的应力应变曲线必须进行双向加载。在单向反复加载时，若所受最大值是不变的则，在后续的加载过程中均为弹性变形。图2应力应变曲线问题：我这样的分析合理吗？[1]2002-implementationofcyclicplasticitymodelsbasedonageneralformofkinematichardening有限单元法王勋成清华大学2005-ABAQUS用户材料子程序开发及应用（硕士论文-华中科技大学）
本文档为【ABAQUS在反复加载卸载应力应变计算说明】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483