ABAQUS在反复加载卸载应力应变计算说明.docVIP

ABAQUS反复加载卸载应力应变计算说明 塑性模型 采用Armstrong and Frederick model (AF模型)[1] 屈服准则（Mises屈服） 屈服函数为 式中：为应力偏量，为总背应力，为屈服极限 流动准则（the associated flow rule） 式中：为塑性应变对时间的微分，为待定量[2]，为应力张量 硬化准则 1 式中：为背应力分量对时间的微分；，为材料常数为已知量，为等效塑性应变对时间的微分。 2 式中：为对应于等效塑性应变的屈服极限，为初始屈服极限为已知量，为材料常数为已知量。 反复加载卸载应力应变计算过程说明 假设在受载前，物体的初始应力、应变以及背应力均为零 2.1加载过程计算 1 外力不足以使得物体中的任何一点的Mises应力值大于屈服极限 此时：塑性应变、背应力均保持为零，屈服极限保持不变。应力由计算，总应变值等于弹性应变。 2 外力使得物体中的任何一点的Mises应力值大于屈服极限 为了说明ABAQUS是如何确定应变增量，有必要对ABAQUS求解材料非线性问题进行简单介绍[3]。 ABAQUS首先将载荷分为若干个微小增量，如图1所示。当结构收到一个微小增量时，ABAQUS用与初始结构位移相对应的初始刚度矩阵和载荷增量计算出结构在这一增量后的位移修正、修正后的位移值和相应的新的刚度矩阵。ABAQUS用新的刚度矩阵计算结构的内力，载荷和的差值为迭代的残余力。如果在模型内的每一个自由度上的值都为零或小于一个给定的容差，如图1所示的点，则结构处于平衡状态即ABAQUS计算到的内外力是平衡的。 若假设整个物体只有一个单元，则位移修正就为应变增量，刚度矩阵就是UMAT程序中的雅可比矩阵。根据塑性变形时的雅可比矩阵计算公式：其中为弹性状态下的雅可比矩阵，可知较小。UMAT计算出的应力便是结构的内力。 在UMAT程序中通过计算： 式中： 并判断若是否大于零来决定节点是否达到屈服。若小于零则节点未达到屈服，若大于等于零节点达到屈服，需要进行塑性迭代与重新计算雅可比矩阵。 图1 增量法迭代原理 ABAQUS处理塑性加载过程问题可以描述为：首先读入上一状态的应力、应变、背应力及屈服强度，使用下标表示（对应时间为）。再根据上一步的刚度矩阵和载荷增量计算出给定的应变增量，然后在UMAT程序中求出（）时刻满足所用塑性模型给出的屈服准则、流动准则、硬化准则下的应力、应变、背应力及屈服强度。 UMAT程序利用已知的，，，，，和给定的首先求出可满足所用塑性模型给出的屈服准则、流动准则、硬化准则的，再使用公式依次求出其余量。 等效塑性应变： 屈服极限： 塑性流动： 式中： 塑性应变增量： 背应力： 应力： 式中： 再计算出：， 加载过程按照此方式进行计算。 2.2当外力卸载时 同样已知：，，，，，并且给定：。这里考虑卸载过程的第一步迭代。由于外力开始卸载，计算得到的。由于上一步迭代材料处于屈服状态，所以在UMAT程序运行结束时有下面关系成立： 所以此时的根据程序判断不进行塑性计算，只计算弹性雅可比矩阵即按照弹性方式进行卸载。 各个参数有关系： ， ，， 弹性应变： 总应变： 按照此方式，直到卸载完成。这里不讨论反向加载。 2.3 再次加载时 由于，，在卸载过程中不变化，所以在第二次加载过程中，读入的，，，，，中，一定是不为零的，并且。对于是否为零，与残余应变是否引起了残余应力有关。对于一个单元残余塑性应变是无法引起残余应力的；对于多个单元构成的物体，当卸载完成后，一些单元的残余塑性应变使得它们与周围的单元之间产生作用力从而产生残余应力。此时读入的是不为零的。这是与第一次加载初始应力状态有区别的。根据UMAT弹性试应力的计算公式，若物体中含有残余应力则在计算过程中直接代入程序计算。 由给定的应变增量计算出下一时刻的。若小于零则节点未达到屈服按照2.1中的弹性加载过程进行计算；若大于等于零按照2.1中的塑性加载过程进行计算。 分析：若要此模型出现图2所示的应力应变曲线必须进行双向加载。在单向反复加载时，若所受最大值是不变的则，在后续的加载过程中均为弹性变形。 图2 应力应变曲线 问题：我这样的分析合理吗？ [1] 2002-implementation of cyclic plasticity models based on a general form of kinematic hardening [2] 有限单元法 王勛成 清华大学 [3] 2005-ABAQUS用户材料子程序开发及应用（硕士论文-华中科技大学）