Deform二次开发步骤.pdfVIP

Deform 3D 二次开发步骤 为了在金属成形工艺模拟过程中进行微观组织演化的定量预测， 所使用的模 拟软件必须有包含微观组织变化的本构模型和专门的微观组织演化分析模块。 当 前国际上虽然有多个知名商业软件流行， 但是它们都不具备微观组织演化的预测 功能。庆幸的是多数商业软件都为二次开发设置了用户子程序功能， 通过用户子 程序，用户就能根据自己的需要增加自己的微观组织预测功能。 为了使 DEFORM3D 软件具有微观组织演化预测功能， 本研究尝试将包含动 态再结晶的热刚—粘塑性材料本构模型植入到 DEFORM3D 中，并在模拟结果中 能够显示晶粒度等用户变量在变形体内的分布。 在研究出具体开发步骤前， 必须 要对 Defom 中的程序有所深入了解。 一、 DEFORM3D二次开发基础理论 1、用户子程序结构 本研究的 DEFORM3D 二次开发涉及到的子程序有： USRMSH、USRMTR 、 UFLOW 、USRUPD （含USR 和 CHAZHI ）。 （1）可以改变几乎所有变量的子程序（ USRMSH ） 子程序功能： 该子程序包含了有限元计算中所有的全局变量， 通过这个用户 子程序，可以修改所有这些变量。 但这些全局变量的改变将直接影响有限元的计 算，处理不当就会使整个程序不能正常进行。 在 DEFORM3D 子程序功能中，所有的用户变量必须在 USRUPD 子程序中 定义。本文的用户子程序中共定义了 18 个用户单元变量。各用户变量的含义如 列表所示。 该子程序用于某些必要数据的获取和存储流程图如下图所示： （2 ）流动应力子程序（ USRMTR 、UFLOW ） SUBROUTINE USRMTR （NPTRTN ，YS ，YPS，FIP，TEPS，EFEPS，TEMP ） SUBROUTINE UFLOW （YS ，YPS，FIP，TEPS，EFEPS，TEMP ） 子程序的变量含义： NPTRTN ：应力模型编号； YS ：流动应力；YPS ：流动 应力对等效应变的导数； FIP：流动应力对等效应变速率的导数； TEPS ：等效应 变； EFEPS：等效应变； TEMP ：温度。 子程序 USRMTR 和 UFLOW 运行时需要输入：应力模型编号、等效应变、 等效应变速率、温度。子程序执行完后将输出：流动应力值、流动应力对等效应 变的导数， 流动应力对等效应变速率的导数。 这几个变量可以用用户定义变量来 计算。 子程序功能： USRMTR ：传入上级程序计算出的等效应变、等效应变率、温度和用户应 力模型编号，调用具有该编号的 UFLOW ，将等效应变、等效应变率、温度值传 入 UFLOW ，用于计算出的流动应力、流动应力对等效应变的导数和流动应力对 等效应变率的导数传回上级程序。 UFLOW ：用传入的等效应变、等效应变率、温度和公共块中的用户定义的 单元变量 (USRE1、USRE2)和节点变量 (USRN1 、USRN2)计算流动应力、流动应 力对等效应变的导数和流动应力对等效应变率的导数，并将这些计算值返回 USRMTR 。 （3 ）用户定义节点和单元变量 (USRUPD) 子程序格式： SUBROUTINE USRUPD 子程序功能： 计算所有的用户单元变量和用户节点变量。 在应力模型中可以 使用这个子程序中定义的变量。