各向同性硬化参数对板料成形极限曲线计算结果影响.docVIP

各向同性硬化参数对板料成形极限曲线计算结果影响 　　摘 要：板料的成形极限曲线可以通过材料的本构方程结合塑性失稳理论进行计算得到。本文使用了金属的弹塑性本构模型结合极限点分叉理论进行计算，研究了各向同性硬化指数n及强度系数k对金属材料的成形极限的影响。 　　关键词：各向同性硬化参数;成形极限曲线;塑性失稳理论 　　1 导论 　　 板料的成形极限反映了它在加工达到塑性失稳前能获得的最大的变形程度，是板材成形的重要性能指标。1965年Keeler[1]及Goodwin[2]提出的以极限应变构成的成形极限曲线（Forming Limit Curve，FLC）为定量衡量板材拉伸失稳极限提供了一种简便实用的工具，至今仍广泛应用于板材成形领域。成形极限曲线可以通过试验测定和理论计算两种方法获得，其中理论计算主要是通过采用不同的屈服准则和塑性本构关系，利用拉伸失稳准则作为颈缩或断裂发生的判据来进行的[3]。随着有限元法及计算机技术的发展，板材成形极限图的理论计算更容易实现，其应用也更加广泛。 　　 材料的应变硬化能力决定了材料在颈缩前依靠硬化使材料均匀变形的能力的大小，显然，材料的应变硬化参数对它的成形极限有着很大的影响。本文用弹塑性本构模型表征材料的力学行为，以极限点分叉理论作为颈缩发生的判据，使用ABAQUS进行了金属材料的成形极限曲线的预测计算，通过控制变量法研究了各向同性硬化指数n及强度参数k对成形极限曲线计算结果的影响。 　　2 材料的力学行为 　　 计算使用的材料为弹塑性材料。当金属材料经历大变形时，其弹性变形远小于塑性变形，其力学行为一般采用次弹塑性模型来表征[4]，次弹性响应定律为： 　　 （1） 　　 其中，为柯西应力张量的Jaumann客观率;C为四阶弹性切线模量张量，它具有主对称性与次对称性;D表示总变形率张量，即速度梯度张量的对称部分;De与Dp则分别表示变形率张量的弹性部分及塑性部分。该材料满足塑性流动法则： 　　 （2） 　　 其中表示标量塑性流动率，它确定了塑性变形率的大小;表示塑性流动势，对于金属材料，取塑性流动势与屈服函数f相同。材料服从Mises各向同性屈服准则： 　　 （3） 　　 其中表示偏应力张量，Y为用来确定屈服面的大小的各向同性强化准则。本文计算时使用如下式所示的Ludwig各向同性硬化式： 　　 （4） 　　 其中为初始屈服极限，n为应变硬化指数，k为应变硬化强度系数，为等效塑性应变。 　　 以上各式结合一致性条件可以推导出弹塑性材料的应力率与总变形率之间的关系用矩阵形式描述为： 　　 （5） 　　 其中HY为一个标量模量，它与材料各向同性强化常数有关，Cep为四阶弹塑性切线模量的Voigt矩阵形式，a为表征加载方式的指数，当a=1时表示塑性加载，a=0表示弹性加载或卸载。 　　 一般情况下，材料在断裂前内部组织会发生恶化，因此本文引入Lemaitre损伤理论提出的损伤标量d来表征金属材料中的微裂纹、孔隙等微损伤，变量d的定义为： 　　 （6） 　　 式中S、Sdef 、Seff分别表示单位体积单元上某截面的面积、该截面的损伤面积及该截面的有效承载面积。在Lemaitre模型中，如果考虑微裂纹和微空洞的方向性，那么d就是法向量的函数，本文中假设损伤是各向同性的，即金属在各方向上微缺陷都是相同的，因此将损伤变量d作为标量考虑且取值区间为[0，1]。d等于0时表示单元完全无损伤，d等于1时表示单元完全受损分裂成两部分。 　　 那么在此模型下材料某截面上的有效应力与名义应力之间则具有如下关系： 　　 （7） 　　 二者关系的物理解释为损伤系数为d的单位体积元某截面名义应力对应的应变等效于无损伤的单位体积元该截面上有效应力对应的应变。损伤参数d的演化方程为： 　　 （8） 　　 其中Hd为一个标量模量，与材料常数β、Sa、sb，弹性应变能释放率Ye及其阈值Yei有关。 　　 使用式（6）新定义的有效应力取代公式（1）（2）（3）中的应力张量，经过推导可以得到新的应力率与总变形率之间的关系为： 　　 （9） 　　 式（5）与式（9）即为我们计算时使用的材料本构方程，计算时将以上本构关系通过UMAT用户子程序实现，程序中采用4阶龙格库塔算法来进行应力迭代更新。 　　3 Valanis极限点分叉理论 　　 极限点分叉理论是一种预测材料的分散性失稳的理论，最早由Valanis在1989年提出，所以又称为Valanis极限点分叉理论，它认为颈缩发生时名义应力是不变的，即有： 　　 （10） 　　 其中，N为名义应力张量，T为第一类P-K应力张量，F为变形梯度，Lt为连接变形梯度材料时间导数与第一类P-K应力率的切线模量，在使用