Deform二次开发步骤
DEFORM二次开发各模块介绍
材料本构模型是实现计算机数值模拟的前提条件之一，【关于计算机数值模拟技术的发展介绍】本论文所采用的有限元模拟软件DEFORM-3D进行材料的微观组织模拟介绍，DEFORM-3D 有限元软件是集成了原材料、成形、热处理和机加工为一体的软件，可用于分析各种塑性体积成形过程中金属流动以及材料的应力、应变和温度等物理场量的分布变化情况，同时提供了材料的流动、模具间的填充、成型过程的载荷量、模具所受应力、材料的纤维流向、成型过程的坯料形成、材料的韧性断裂以及金属微观组织结构等信息。
为了实现在金属成形工艺模拟过程中进行微观组织演化的定量预测，所使用的模拟软件必须有包含微观组织变化的本构模型和专门的微观组织演化分析模块。
当前国际上虽然有多个知名商业有限元软件流行，但是它们都不具备微观组织演化的预测功能;或者软件具有微观组织变化的本构模型，但仍需使用者输入材料的参数方可进行，而软件不提供材料的参数;故很多软件都淡化此微观组织演化分析模块。
庆幸的是多数商业软件都为二次开发设置了用户子程序功能，通过编制用户子程序就能实现对微观组织演化的预测功能。
用户自定义本构模型的输入方法在当今的科学研究方向中，新材料的开发占据了一个重要的角色。
不同的材料工作者开发了不同的新材料，得到了不同的本构模型，需要对这种新材料进行模拟，为了满足这种需求，DEFORM提供了两种用户自定义本构模型的输入方法：（1）以函数形式输入本构模型。
DEFORM提供了若干常用本构模型，如图4.6所示。
若用户的本构模型与系统提供的本构模型一致，则可直接输入其相关系数即可；若用户的本构模型在系统中不存在，则可通过二次开发编程的方式将用户的本构模型加入到DEFORM中，然后在图4.6中选择“User routine”并输入所调用的本构模型子程序的编号。
(2)以数据形式输入本构模型。
DEFORM还允许通过输入数据的方式来定义材料的塑性流动行为。
具体方法是根据材料的真应力一真应变曲线，取若干个数据点，逐个输入该材料在某个温度、某个变形速率和某个真应变下的真应力。
Deform详细教程
输入材料的热物理性质，如热导率、比热容、热膨胀系数 等，以便在模拟过程中考虑温度对材料性能的影响。
材料失效准则
根据实际需要，选择适当的材料失效准则，如最大主应力 准则、等效塑性应变准则等，并设置相应的失效参数。
边界条件设置
几何边界条件
定义模型的几何形状、尺寸和边 界类型，如固定边界、自由边界 、对称边界等。
04 Deform软件基本操作
用户界面介绍
主界面
包括菜单栏、工具栏、模型树、属性窗口等，提供全 面的操作功能和视图展示。
图形界面
支持多种图形显示模式，如实体、网格、轮廓等，方 便用户进行模型分析和后处理。
自定义界面
用户可根据个人习惯自定义界面布局，提高工作效率 。
基本操作命令
鼠标操作
通过鼠标左键选择、拖拽、旋转等操作，实现模型的交互操作。
未来发展趋势预测
A
随着计算机技术的不断发展，有限元分析软件 的计算能力和效率将不断提高，使得更大规模 、更复杂的仿真分析成为可能。
人工智能、机器学习等技术的引入，将为 有限元分析提供更强大的数据处理和挖掘 能力，进一步提高分析的精度和效率。
B
C
多物理场耦合分析将成为未来发展的重要方 向，Deform等软件将不断完善多物理场分 析功能，满足更广泛的应用需求。
配置环境变量和启动软件
启动软件 在完成安装和环境变量配置后，可以通过以下方式启动Deform软件 1. 点击桌面或开始菜单中的Deform图标。
配置环境变量和启动软件
2. 在命令行中输入Deform的可执行 文件名并回车。
3. 如果设置了文件关联，可以直接双 击与Deform关联的文件类型来启动 软件并打开相应文件。
基于DEFORM3D二次开发的塑性成形过程组织演化模拟
与微观耦合 ! 不具 备 微 观 组 织 演 化 的 模 拟 和 预 测 功
!! 引 ! 言
高温成形过 程 中 ! 金 属 将 发 生 动 态 和 静 态 再 结 晶 ! 产生新的晶 粒 " 这 种 微 观 组 织 的 演 变 在 很 大 程 度上决定 了 产 品 的 宏 观 力 学 性 能
’ !! "(
能 ! 或者只具有 简 单 的 预 测 能 力 ! 其 模 型 并 不 一 定 适合于所 考 察 的 问 题 " 本 文 通 过 对 ] C _ 7 > ? O ] 二次 开发 ! 将适合于材 料 的 组 织 模 型 与 成 形 的 热 力 耦 合 计算结合 ! 模拟热成形过程中的组织演化 "
"利用热加工
" -! 如果等 效 应 变 大 于 临 界 应 变 " 则 开 始 计 算 动 态 再结晶过程 # 当动态再结晶分数 F ## 9 W $ 时"则 直
接转为晶粒长大 的 有 关 计 算 ! 把 静 态 再 结 晶 模 型 与 热传导有限元耦合 " 可以分析卸载后的静态再结晶 ! ] C _ 7 > ? O ] 是由 / 0 2 公司开发的一个体积成 形有限元工艺模 拟 专 用 商 业 软 件 " 其 用 户 定 义 子 程 序的代码储存在 M C _＿ : Z > 9 _中 " 有 限 元 主 程 序 通 过 调用该文件中的 子 程 序 " 可 以 计 算 出 用 户 自 定 义 变 量的值 ! 所以关 键 问 题 在 于 该 子 程 序 的 编 写 ! 其 中 ( / )P/ E 子程序包含了 有 限 元 计 算 中 所 有 的 变 量 ! 通过这个用户 子 程 序 " 可 以 修 改 所 有 这 些 变 量 ! 该 子程序申明如下 $ / ( R ) ’(0 1 * 4 % ( / )P/ E ) F" ] ) F" () F" 0 4P." ] 0 S 万方数据 " " " " " P. ) F + ) F R 4 / 0 / 4 4 . / 0 4 . /"
Deform-3d热处理模拟操作全解
Deform-3d热处理模拟操作热处理工艺在机械制造中占有十分重要的地位。
随着机械制造现代化和热处理质量管理现代化的发展，对热处理工艺提出了更高的要求。
热处理工艺过程由于受到加热方式、冷却方式、加热温度、冷却温度、加热时间、冷却时间等影响，金属内部的组织也会发生不同的变化，因此是个十分复杂的过程，同时工艺参数的差异，也会造成热处理加工对象硬度过高过低、硬度不均匀等现象。
Deform-3d软件提供一种热处理模拟模块，可以帮助热处理工艺员，通过有限元数值模拟来获得正确的热处理参数，从而来指导热处理生产实际。
减少批量报废的质量事故发生。
热处理模拟，涉及到热应力变形、热扩散和相变等方面，因此计算很复杂，软件采用牛顿迭代法，即牛顿-拉夫逊法进行求解。
它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。
多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。
方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。
牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根等。
但由于目前Deform-3d软件的材料库只带有45钢、15NiCr13和GCr15等三种材料模型，而且受到相变模型的局限，因此只能做淬火和渗碳淬火分析，更多分析需要进行二次开发。
本例以45钢热处理淬火工艺的模拟过程为例，通过应用Deform-3d 热处理模块，让读者基本了解热处理工艺过程有限元模拟的基本方法与步骤。
1 、问题设置点击“文档”（File）或“新问题”（New problem），创建新问题。
在弹出的图框中，选择“热处理导向”(heat treatment wizard)，见图1。
sgniht图1 设置新问题2、初始化设置完成问题设置后，进入前处理设置界面。
首先修改公英制，将默认的英制t h i ng si nt he i rb ei n ga re go od fo rs （English ）修改成公制（SI ），同时选中“形变”（Deformation ）、“扩散”(Diffusion)和“相变”(Phase transformation)，见图2。
DEFORM二次开发的一点小心得!
DEFORM二次开发的一点小心得!前一段时间闲着没事,看了DEFORM的help文件.和众位兄弟交流一下.用户子程序实现过程:第一步：所需文件准备：把DEF_SIM目录下的def_usr.f ,DEF_SIM_USR_ABsoft70.gui , DEF_SIM_P4_USR_ABsoft70.gui ,DEF_SIM_P4_USR_LIB.lib , DEF_SIM_USR_LIB.lib和lib目录下的所有.lib文件拷到工作目录下。
备份DEF_SIM.exe ，DEF_SIM_P4.exe两个文件。
第二步：编写源程序。
对于新手可以使用模版自带的子程序，对于材料流动应力子程序模版里头本身就自带有一个，所以可以不需要自己编写。
(我就利用过自带的作过模拟,把前面的常数设置为10和100分别模拟)第三步：用absoft7.0或以上版本打开DEF_SIM_USR_ABsoft70.gui 文档，指定好library files，用工作目录下的lib文件来替换原来默认的lib文件，上述做完后直接点build就ok了，自动生成了DEF_SIM.exe。
重复上面的过程打开DEF_SIM_P4_USR_ABsoft70.gui 文档生成DEF_SIM_P4.exe。
第三步：用生成的DEF_SIM.exe ，DEF_SIM_P4.exe替换原来的这两个文件。
第四步，运行模型。
对于流动应力子程序，只要在material那里选择子程序就ok 了。
注:3D中必须是DEF_SIM.exe，而DEF_SIM— P4.exe不好！(我也不知道原因,师兄是这么说的!).还有一个问题,关于DEFORM的text运行模式.在安装目录下面,有如下几个需要了解的文件:<1>DEF_PRE.EXE;<2>DEF_ARM_DEF_PRE.EXE这是前处理,DOS界面输入前处理参数设置,记录每次设置时的操作,按一定格式记录在记事本上面,这个很重要,二次开发可能需要这个文件.具体参见help.DEF_ARM_ 这是运行命令.具体格式见help.下面是关键:通过fortran子程序调用DEFORM.具体如下.我用的是visual fortran6.6.result=system('E:\DEFORM3D\V5_0\DEF_PRE.EXE') (启动前处理,弹出DEF_PRE.EXE界面)result=system('E:\DEFORM3D\V5_0\DEF_PRE.EXE<INI.TXT') (进行前处理参数设置,INI.TXT是记录每次设置时的操作命令流),譬如一个简单的操作命令流:<CR>21XXXX.KEY (key文件,需要copy到你的子程序目录下)<CR>E72XXXX.DB (生成DB文件)<CR>EY<CR>result=system('E:\DEFORM3D\V5_0\DEF_ARM_ XXXX B') (运行DB文件,XXXX为文件名,B表示batch模式运行)注意运行过程中,有可能生成临时文件,如FOR003等,会影响模拟进行,可通过fortran命令open和write操作清除.以上可以实现一些特殊的操作.如运动子程序等.上述即是我的一点小体会.很多我也不是很懂,自己的课题也不是DEFORM二次开发,也不专长编程.呵呵,希望对大家有所帮助!Re:请问能不能帖个自己二次开发加本构方程的例子这是我尝试过的关于流动应力的二次开发，希望能给没有做过二次开发的朋友一点经验~~~现在要好好研究二次开发和有限元软件的核心机制了~~~共同进步~~~program USRMTRreal YS,TEPS,EFEPS,TEMread (*,*,*)TEPS,EFEPS,TEMYS=123.5237-1.238956*TEM+3.8584*LOG10(EFEPS)+6.7852*TEPS-1.29*TEM*TEPS+52.895* LOG10(EFEPS)*TEPS+0.00088442*TEM*TEMPrint *,YSEND相信看了这个程序后，大家对二次开发肯定有个初步的了解了~呵~用户子程序具体实现心得结合论坛上各位大虾的经验和自带帮助,前段时间动手试了一下用户子程序的实现过程,现将自己的一些体会与大家共享,相信对做二次开发的新手有用用户子程序实现过程:第一步：所需文件准备：把DEF_S IM目录下的def_usr.f , DEF_S IM_USR_ABsoft70.gui ,DEF_SIM_P4_USR_ABsoft70.gui , DEF_S IM_P4_USR_LIB.lib , DEF_S IM_USR_LIB.lib和lib目录下的所有.lib文件拷到工作目录下。
deform2D手册_下
11.热传导问题二次计算主要纲目：11.1打开旧问题11.2加载数据库11.3设置模拟控制和边界条件11.4写数据库11.5继续模拟11.6后处理11.6.1点轨迹捕捉11.7退出DEFORM/2D本章使用的图标阴影等值线图测量工具图标放大窗口11.热传导问题二次计算11.1打开-问题转向目录HTX并,键入DEFORM2，按回车运行DEFORM -2D,系统,单击Pre-Processor按钮进入前处理。
11.2加载数据库要对一个老问题继续进行模拟，应选择已有数据库中的最后一步作为重新计算模拟的第一步。
Pre-Processor打开如下选择提示对话框。
在上面信息窗口中，单击OK按钮加载数据库。
接着显示步数列表窗口。
由于240步是数据库中的最后一步,所以点选240，单击OK按钮,把从该步起的数据输入到前处理中。
图11.1 - Select Database Step 窗口11.3设置模拟控制参数和边界条件单击Simulation Controls按钮,打开SIMULATION CONTROLS窗口。
接着单击Stopping Step按钮,程序会自动设置Starting Step Number(NSTART)为-240,设置Number of Simulation Steps (NSTEP)为60,Step Increment to Save (STPINC)为4,改变Time per Step (DTMAX)为1秒。
完成后单击OK按钮。
单击Processing Conditions按钮,打开其窗口,设置Environment Temperature(ENVTMP) 为68,并单击OK按钮,最后单击OK按钮退出SIMULA TION CONTROLS窗口。
重新设置工件的边界条件来模拟在空气中冷确的过程,这时工件的所有表面均暴露在环境中[前例中仅部分与环境有热交换],故而需初始化边界条件来重新设置（图11.2）。
Deform 二次开发步骤
微观组织模拟：模拟步骤Deform3D二次开发步骤为了在金属成形工艺模拟过程中进行微观组织演化的定量预测，所使用的模拟软件必须有包含微观组织变化的本构模型和专门的微观组织演化分析模块。
当前国际上虽然有多个知名商业软件流行，但是它们都不具备微观组织演化的预测功能。
庆幸的是多数商业软件都为二次开发设置了用户子程序功能,通过用户子程序，用户就能根据自己的需要增加自己的微观组织预测功能。
为了使DEFORM3D软件具有微观组织演化预测功能，本研究尝试将包含动态再结晶的热刚一粘塑性材料本构模型植入到DEFORM3D中，并在模拟结果中能够显示晶粒度等用户变量在变形体内的分布。
在研究出具体开发步骤前，必须要对Defom 中的程序有所深入了解。
一、DEFORM3D二次开发基础理论1、用户子程序结构本研究的DEFORM3D二次开发涉及到的子程序有：USRMSH、USRMTR、UFLOW、USRUPD（含USR和CHAZHI）。
（1）可以改变几乎所有变量的子程序（USRMSH）子程序功能：该子程序包含了有限元计算中所有的全局变量，通过这个用户子程序，可以修改所有这些变量。
但这些全局变量的改变将直接影响有限元的计算，处理不当就会使整个程序不能正常进行。
在DEF0RM3D子程序功能中，所有的用户变量必须在USRUPD子程序中定义。
本文的用户子程序中共定义了18个用户单元变量。
各用户变量的含义如列表所示。
该子程序用于某些必要数据的获取和存储流程图如下图所示：USRMSH子程序流程图(2)流动应力子程序(USRMTR、UFLOW)SUBROUTINEUSRMTR(NPTRTN,YS,YPS,FIP,TEPS,EFEPS,TEMP) SUBROUTINEUFLOW(YS,YPS,FIP,TEPS,EFEPS,TEMP)子程序的变量含义：NPTRTN：应力模型编号；YS：流动应力；YPS：流动应力对等效应变的导数；FIP：流动应力对等效应变速率的导数；TEPS：等效应变；EFEPS：等效应变；TEMP：温度。
proe二次开发教程
proe二次开发教程——手把手教你如何二次开发本实例是针对proe4.0在VS2005环境下的开发，对于vs2003以上版本的开发本例均适用，本例从开始的打开vs2005软件开始讲解，一步不漏，并加入适当的解释和必要的说明，相信照着做完之后应该没有不成功的了！注：软件安装注意事项：proe的安装目录不能有空格，比如：“D:\Program Files\proeWildfire 4.0”是不对的，运行的时候会出错，要把路径中的空格去掉，如“D:\proeWildfire4.0”一、配置相关信息：打开VS2005，选择工具——》选项，在“项目和解决方案”中选择“VC++目录”，在“显示一下内容的目录”中选择“包含文件”，加入以下文件：D:\proeWildfire4.0\protoolkit\includesD:\proeWildfire4.0\protoolkit\protk_appls\includesD:\proeWildfire4.0\prodevelop\includes在“显示一下内容的目录”中选择“库文件”，加入以下文件：D:\proeWildfire4.0\protoolkit\i486_nt\objD:\proeWildfire4.0\prodevelop\i486_nt\obj注：这是开始前的环境配置，一次配置之后永久生效，以后再新建项目就不需要再配置了二、新建项目，进行项目设置：文件——》新建——》项目：选择MFC DLL 模板，输入项目名称：确定后点完成，项目生成后选择项目——》属性，进行属性设置：在test（本例的项目名称）属性页中选择配置属性——》连接器——》输入，在附加依赖项里加入：wsock32.libmpr.libpsapi.libprotk_dllmd.lib注：对于不同的环境附加依赖项是不同的，以上是VS2005版本的开发软件附加依赖项，如果是VC++6.0则是“protk_dll.lib protoolkit.lib mpr.lib wsock32.lib ”，同时还要加入忽略库“msvcrt.lib”，还要选择“强制输出”，否则程序会因为一些警告终止编译。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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子程序功能：计算所有的用户单元变量和用户节点变量。
在应力模型中可以使用这个子程序中定义的变量。
（4）USR子程序该子程序被USRUPD调用，它是计算材料本构模型的主体程序，具体功能是：计算流动应力、再结晶体积分数、平均晶粒度、未再结晶晶粒度、再结晶晶粒度、晶粒度差等与本构模型相关的用户自定义单元变量，并将这些变量值传到USRUPD子程序中。
该子程序的计算需要用到存储在IUSRVL数组中的热刚-粘塑性本构模型中的所有参数。
子程序中读取语句如下：READ（IUSRVL （4），*）RM1，RK1，GAMA，RMD，ROU0，ROU1，ROU2 READ（IUSRVL （5），*）RLID，RM，BEITA，A，RM2，RK2（5）CHAZHI子程序变形体网格重划后，存储在用户构造的数据文件中的建立在网格重划前单元上的变量数据必须通过插值方法转换到网格重划以后的新单元之上，不然有限元计算将不能正确进行。
因此，网格重划必须进行必要的数据转换。
为进行插值计算，需要存储网格重划前、后的单元节点信息，包括：网格重划前的单元总数、节点总数、每一时间步的节点坐标、单元节点号；网格重划后的单元总数、节点总数、节点坐标、单元节点号。
CHAZHI子程序的流程图如下图所示。
CHAZHI子程序流程图2、用户定义子程序的核心代码用户定义子程序的核心代码都储存在DEF_USR.FOR中，有限元主程序通过调用该文件中的子程序，可以计算出用户自定义变量的值。
所以关键问题在于该子程序的编写。
DEF_USR.FOR文件中的子程序有USRMSH、USRUPD 和UFLOW三个用户子程序。
下图为DEF_USR.FOR文件中子程序之间的先后调用顺序关系。
3、动态再结晶过程分析子程序的编制该子程序是在有限元模拟程序中用于计算流动应力和平均晶粒度等单元变量的，它在有限元模拟程序中的位置如流程图中加粗字体所示。
该子程序在有限元计算中的位置本研究的一个重要内容就是将公式所组成的“热刚—粘塑性本构模型”以用户子程序的方式加入到DEFORM3D中，并将该模型中的一些变量(如再结晶体积分数等)定义成用户自定义单元变量，以便从DEFORM3D的后处理中可以查看这些变量在变形体内的分布情况。
具体开发步骤如下所述。
二、开发步骤1、开发运行环境安装用虚拟光驱加载Absoft8.0.iso文件，安装Absoft 8.0编译器，安装界面如图1~图2所示。
图1 安装界面一，点击【Install Products】图2 安装界面二，顺次安装每个条目，【IMSL Libraries】可以略过。
2、安装Deform3D V6、植入子程序（1）安装Deform3D V6软件（2）植入子程序：打开Deform3D V6安装目录，如“D:\DEFORM3D”，进入子目录“V6_1\UserRoutine\DEF_SIM”中。
此时可见文件如图3所示。
图3双击打开文件“DEF_SIM_USR_Absoft75.gui”，如图4所示。
图4 DEF_SIM_USR_Absoft75在此处，我们用到的文件为：usr_mtr.f、usr_upd.f、usr_msh.f。
源文件中的各个函数在工程中文件分布情况如表1所示。
文件名包含的函数Usr_mtr.f USRMTR UFLOW1 UFLOW2 UFLOW3Usr_upd.fUSRUPD CHAZHI USR CONTOUR COLUMFOUsr_msh.f USRMSH表1将“源代码.pdf”中的代码拷贝到对应的文件中，并保存。
点击工具栏中的图标进行编译，在同一目录中生成“DEF_SIM.exe”。
将“DEF_SIM.exe”拷贝到Deform安装目录下的“V6_1”目录下，如“D:\DEFORM3D\V6_1”，此时提示该文件已经存在，覆盖即可。
当然，最好在覆盖之前，先将原有的“DEF_SIM.exe”备份一下。
经验显示：应将原“DEF_SIM.exe”文件剪切走，然后将新生成“DEF_SIM.exe”文件拷入。
此时就已经完成了准备工作，下面开始用Deform来调用自定义的模型来进行仿真模拟。
具体步骤如下：3、二次开发嵌入子程序后的仿真流程（1）新建问题打开Deform-3D，点击左上角，新建问题。
点击两次【next】，出现问题名称，如图5所示，为新建问题对话框，在新建问题名称处根据个人意愿输入，例如：输入“A”。
点击【Finish】，进入前处理环境。
图5 新建问题对话框（2）模拟控制输入选择菜单【Input】中的【Simulation Controls】，在【Main】栏更改单位制为SI，如图6所示为模拟控制输入对话框。
点击【OK】按钮。
图6模拟控制输入对话框（3）添加坯料在右边数结构中选中【[PDie] Workpiece】，并在中选择【Import Geo…】，选择Deform安装目录下“\V6_1\LABS\”中的Block_Billet.STL。
点击，划分网格，点击【Generate Mesh】生成网格。
（4）添加上模点击，添加模型，并在中选择【Import Geo…】，选择Deform安装目录下“\V6_1\LABS\”中的Block_TopDie.STL。
在中的【Constant Value】定义运行速度，如图7所示。
图7（5）添加下模点击，添加模型，并在中选择【Import Geo…】，选择Deform 安装目录下“\V6_1\LABS\”中的Block_BottomDie.STL。
（6）定义材料（关键点——应力应变本构方程）选择菜单【Input】中的【Material】，选择【load from lab】，选择合适的材料。
然后在【Plastic】中的【Flow Stress】下拉菜单中选择“User routine”.如图8所示。
图8点击，选择流动应力子函数，这里将与定义的UFLOW1、UFLOW2…对应起来，如：选择填写数字“1”，表示对应UFLOW1。
点击【Close】。
在右边数结构中选中【[PDie] Workpiece】，在【General】中的【Material】中点击，选择刚才定义的材料。
（7）定义单元数据打开【Advanced】中的Element Data，出现对话框，选择“User”,对话框中可以定义“Variable Name”和对应的“Value”。
如图9所示。
在此定义并输入18个用户单元变量和相应的初始值。
图9（8）定义步骤打开菜单【Input】中的【Simulation Controls】，选择【Step】，设置步骤信息如图10所示。
图10（9）定义接触打开菜单【Input】中的【inter-object】，弹出图11对话框，点击Yes。
弹出对话框图12。
图11图 12选中条目1，点击【Edit】，设置摩擦系数为0.12，如图13所示。
图 13点击【Close】，再点击【Generate】生成接触面。
用同样方法设置另一个接触面。
（10）定义“用户自定义变量”打开菜单【Input】中的【Simulation Controls】，选择【Advanced】中的【User Defined】，输入编号为1-6的用户变量。
其具体意义如表2所示。
表2用户变量编号对应的变量含义值1 ISE1 被调用次数02 DEFEPSE 应变率增量0.1或其他3 KSTEP 当前步数04RM1RK1GAMARMDROU0ROU1ROU2m1K1(N/m2)γ(J/m2)ρ0ρ1ρ20.170700.649760.86E31.58784E6-0.387845RLIDRMBEITAARM2RK2Ii/Idm3βiam2K2(N/m2)42.37E69.47E122.9690.19730.277466 IR 网格重新划分次数这部分用户自定义变量涉及到具体材料的相关参数。