Deform二次开发步骤.docVIP

  1. 1、本文档共7页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多

Deform3D二次开发步骤

为了在金属成形工艺模拟过程中进行微观组织演化的定量预测,所使用的模拟软件必须有包含微观组织变化的本构模型和专门的微观组织演化分析模块。当前国际上虽然有多个知名商业软件流行,但是它们都不具备微观组织演化的预测功能。庆幸的是多数商业软件都为二次开发设置了用户子程序功能,通过用户子程序,用户就能根据自己的需要增加自己的微观组织预测功能。

为了使DEFORM3D软件具有微观组织演化预测功能,本研究尝试将包含动态再结晶的热刚—粘塑性材料本构模型植入到DEFORM3D中,并在模拟结果中能够显示晶粒度等用户变量在变形体内的分布。在研究出具体开发步骤前,必须要对Defom中的程序有所深入了解。

一、DEFORM3D二次开发基础理论

1、用户子程序结构

本研究的DEFORM3D二次开发涉及到的子程序有:USRMSH、USRMTR、UFLOW、USRUPD(含USR和CHAZHI)。

(1)可以改变几乎所有变量的子程序(USRMSH)

子程序功能:该子程序包含了有限元计算中所有的全局变量,通过这个用户子程序,可以修改所有这些变量。但这些全局变量的改变将直接影响有限元的计算,处理不当就会使整个程序不能正常进行。

在DEFORM3D子程序功能中,所有的用户变量必须在USRUPD子程序中定义。本文的用户子程序中共定义了18个用户单元变量。各用户变量的含义如列表所示。

该子程序用于某些必要数据的获取和存储流程图如下图所示:

(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:温度。

子程序USRMTR和UFLOW运行时需要输入:应力模型编号、等效应变、等效应变速率、温度。子程序执行完后将输出:流动应力值、流动应力对等效应变的导数,流动应力对等效应变速率的导数。这几个变量可以用用户定义变量来计算。

子程序功能:

USRMTR:传入上级程序计算出的等效应变、等效应变率、温度和用户应力模型编号,调用具有该编号的UFLOW,将等效应变、等效应变率、温度值传入UFLOW,用于计算出的流动应力、流动应力对等效应变的导数和流动应力对等效应变率的导数传回上级程序。

UFLOW:用传入的等效应变、等效应变率、温度和公共块中的用户定义的单元变量(USRE1、USRE2)和节点变量(USRN1、USRN2)计算流动应力、流动应力对等效应变的导数和流动应力对等效应变率的导数,并将这些计算值返回USRMTR。

(3)用户定义节点和单元变量(USRUPD)

子程序格式:SUBROUTINEUSRUPD

子程序功能:计算所有的用户单元变量和用户节点变量。在应力模型中可以使用这个子程序中定义的变量。

(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、用户定义子程序的核心代码

3、动态再结晶过程分析子程序的编制

该子程序是在有限元模拟程序中用于计算流动应力和平均晶粒度等单元变量的,它在有限元模拟程序中的位置如流程图中加粗字体所示。

该子程序在有限元计算中的位置

本研究的一个重要内容就是将公式所组成的“热刚—粘塑性本构模型”以用户子程序的方式加入到DEFORM3D中,并将该模型中的一些变量(如再结晶体积分数等)定义成用户自定义单元变量,

文档评论(0)

趁早学习 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档