基于abaqus的Huang晶体塑性UMAT改VUMAT

您所在的位置：网站首页 › abaqus有限滑移 › 基于abaqus的Huang晶体塑性UMAT改VUMAT

基于abaqus的Huang晶体塑性UMAT改VUMAT

2024-06-17 11:52| 来源: 网络整理| 查看: 265

黄永刚院士编写的单晶晶体塑性UMAT，主要用于在Abaqus有限元仿真中进行单晶及多晶晶体塑性变形的计算，是许多科研工作者学习晶体塑性模拟的教学资源。可以在其基础上对硬化模型进行修改，甚至引入损伤。

UMAT主要应用于隐式分析，而对于大变形接触问题，隐式分析往往计算效率较低。对于接触、碰撞、冲击等问题采用VUMAT往往具有更高的计算效率和收敛速度。本文旨在将Huang编写的UMAT改写为VUMAT，并进行对比验证。

将UMAT改写为VUMAT需要从以下方面考虑：（1）UMAT是在积分点上调用的，而VUMAT一次调用会计算很多个积分点上的变量，需要对子程序接口形式进行修改；（2）Huang本构中的转动张量DROT是用于对滑移面和滑移方向进行旋转的，在UMAT中，Abaqus会提供转动张量DROT，在VUMAT中，子程序接口没有提供DROT，需要通过VUMAT传入的变形梯度更新滑移面和滑移方向；（3）显示分析采用了Green-Naghdi率，而隐式分析采用Jaumann率，需要对应变率进行修改。

采用了两个多晶模型进行一致性的验证，第一个模型是125个网格的单位长度代表体积单元，每25个网格设置1个取向。第二个模型是采用Voronoi方法获得的15个不同取向晶粒的多晶模型。

（1）15个不同取向晶粒的多晶模型

15个不同取向晶粒的多晶模型，采用狗骨单轴拉伸试件进行数值试验，有限元模型如下图所示。开展单轴拉伸，UMAT采用隐式分析，VUMAT采用显式分析。