基于DEFORM3D二次开发的塑性成形过程组织演化模拟【文献】

摘 要:

金属热成形方法可以有效改善产品的综合机械性能,利用有限元模拟可以为控制锻造和产品质量提供理 论依据。在Deform3D的热力耦合刚粘塑性有限元模拟技术的基础上,进行了微观组织演化的二次开发,可以扩展 有限元软件的组织模拟能力,并利用该方法对20CrMnTi钢镦锻热成形过程进行了计算机模拟,得到了热力参数的 分布状况和内部晶粒度变化的规律。通过摇臂轴的镦锻成形模拟证明了组织模拟能够为工艺改进提供理论依据。 关键词:刚粘塑性;有限元;晶粒尺寸;显微组织演化;热镦锻 

中图分类号:TG111.7 文献标识码:A 文章编号:1007-2012 (2005)Sup-0040-04

1 引言

高温成形过程中，金属将发生动态和静态再结晶，产生新的晶粒。这种微观组织的演变在很大程度上决定了产品的宏观力学性能\(^{[1-2]}\)。利用热加工过程控制晶粒大小，细化微观组织，是提高产品力学性能的重要手段。因此，研究材料在热成形过程中的宏观力学行为和微观组织的变化，揭示其相互之间的关系，并依据优化工艺参数、设计塑性成形工艺和锻后冷却方案，对解决目前的工艺问题，提高产品质量是很有意义的，同时也是变形过程全面模拟的前沿课题\(^{[3]}\)。有限元数值模拟技术是随着物理模拟设备的完善以及计算机技术的发展而发展起来的。鉴于有限元法是目前唯一能对塑性加工过程给出全面且较为精确数值解的分析方法，本文对材料组织性能所进行的数值模拟均采用该分析方法。数值模拟软件是求解塑性加工问题的一个基本工具。现在市场上已有许多成熟的用于金属塑性加工的商业软件。如 DEFORM，MSC.MARC，MSC.SUPERFORM，Dynaform 等，但这些软件都只进行宏观变形和温度的分析计算，没有考虑宏观与微观耦合，不具备微观组织演化的模拟和预测功能，或者只具有简单的预测能力，其模型并不一定适合于所考察的问题。本文通过对 Deform3D 二次开发，将适合于材料的组织模型与成形的热力耦合计算结合，模拟热成形过程中的组织演化。

2 模型建立

材料在热塑性变形中除了应力和应变的变化外，还要发生复杂的微观组织变化，即动态再结晶、静态再结晶、晶粒生长等。研究表明，再结晶的晶粒尺寸和再结晶晶粒百分比，除了和原始晶粒尺寸与微量元素含量有关外，主要取决于变形及冷却过程的温度、应变和应变速率。

由于热成形过程的复杂性，很难进行现场测试，而在实验室进行研究，只能得到简单条件下的显微组织与宏观热力参数间的定量关系。通过有限元模拟，可以十分方便经济地获得场量信息，因此可以将有限元用于组织演化的预报和控制。图 1 给出了用有限元模拟材料热变形过程中微观组织演化的计算流程。

在实验室通过热力模拟试验，可以获得材料的准经验数学模型，其经验模型大都是通过试验数据回归获得 Avrami 形式的数学模型。在对有限元软件开发过程中，将组织模型离散，在每一迭代步中计算，其计算过程如图 2 所示。

把再结晶模型与刚粘塑性有限元程序耦合，即在子程序计算时，首先分析变形过程中的动态再结晶，在实际计算中首先判断再结晶的体积百分数，当 X＜0.95 时，判断等效应变 ε 是否大于临界应变ε_c。如果等效应变大于临界应变，则开始计算动态再结晶过程；当动态再结晶分数 \(X>0.95\) 时，则直接转为晶粒长大的有关计算。把静态再结晶模型与热传导有限元耦合，可以分析卸载后的静态再结晶。

Deform3D 是由 SFTC 公司开发的一个体积成形有限元工艺模拟专用商业软件，其用户定义子程序的代码储存在 def_usr.f 中，有限元主程序通过调用该文件中的子程序，可以计算出用户自定义变量的值。所以关键问题在于该子程序的编写。其中 USRMSH 子程序包含了有限元计算中所有的变量。通过这个用户子程序，可以修改所有这些变量。该子程序申明如下：

SUBROUTINE

USRMSH (RZ, DRZ, URZ, TEMP, DTMP, FRZA, FRZB, EFSTS, EFEPS, TEPS,RDTY, STS, EPS, DCRP, TSRS, DAMG, USRVE,USRVN, ATOM, HEATND,EPRE, VOLT, WEAR, DUM1, PRZA, DUM2, DUM3,HDNS,VF,DVF,VFN,TICF,GRN- SZ, CURTIM, DTMAbC, BCD, NBCDT, NOD, MATR, NBDRY, KOBJ, NUMEL, NUMNP, NDSTART, NDEND, NUMFAC, NUSRVE, NUSRND, NPHASE, ISTATUS, NROUTINE,NONP) 

其中 TEMP 为温度,DTMP 为该迭代步的温 度增量,EFEPS (NUMEL)是等效应变率,TEPS (NUMEL)是塑性应变增量。在子程序中需读入组 织模拟需要的以上4个基本参量。对于组织模拟来 

说,通过 USRMSH 可以修改以下参量: VF (NPHASE,*):再结晶分数 DFV (NPHASE,*):再结晶分数增量 DFN (NPHASE,*):静态再结晶分数 GRNSZ (*):晶粒尺寸 

3 实例分析

扇形齿轮摇臂轴广泛应用于汽车等各种机动车 的发动机、转向系统、循环球方向系统。由于在工 作时受力复杂,工作环境恶劣,要求摇臂轴有较高 的力学性能。因此对其热成形过程中的组织演化过 程进行模拟,指导工艺有着重要的意义[4,5]。 

摇臂轴采用一步成形,材料采用 20CrMnTi, 模具材料为5CrNiMo。因为工件为平面对称,所 以,取其1/2进行模拟,坯料划分为52111个六面体单元,共7085个节点,上模1918个单元,下模 5823个单元。坯料的材料为20CrMnTi,其组织模型如下所示: 

将组织数学模型离散化,通过修改 USRMSH 子程序获得了模拟结果。

图3是模拟使用的模型图。反映材料热成形过程的主要特征参量为温度场 (见图4)、应变与应变速率分布 (见图5、图6)等。

而热成形过程中材料内部显微组织的变化,可 由再结晶程度和晶粒度大小来说明 (图7和图8)。 图中坯料内的温度场与图中应变与应变率的等值线 分布符合塑性加工过程中热力因素的变化规律。由图7可以看出,在这种等效应变分布条件下坯料内各区域都发生了不同程度的再结晶。

4 结 论

采用热力耦合刚粘性有限元微观组织模拟技术, 对20CrMnTi结构钢摇臂轴镦锻成形进行了模拟研 究,所得到的有关场量信息符合变形规律。数值模 拟结果表明,有限元组织模拟可以为摇臂轴热成形 工艺的实施及产品质量的控制,生产合格产品提供 理论依据。 

原资料见附件

免责声明：

本页面/内容部分素材来源于互联网公 开 信 息，旨在传递更多信息，不代表本平台立场。

版权归原作者或机构所有，如涉及侵权，请通过平台联系我们，我们将在核实后第一时间处理。

本平台对转载内容的真实性、准确性不作任何保证，用户需自行判断并承担使用风险。