umat到vumat转化程序_VUMAT与UMAT相互转换注意事项 -

VUMAT与UMAT相互转换注意事项

问题描述：最近在学习和编写子程序(V)UMAT的过程中，发现VUMAT和UMAT程序中应变增量的含意有不同。这个不同点对于需要在VUMAT和UMAT之间做转换的人来说容易搞错，虽然很基础，不清楚也容易搞死人，所以我在本帖写个心得，提醒朋友们注意后处理输出的应变、UMAT的应变和VUMAT的应变有如下差别

理论准备：根据弹性理论，应力应变关系如下：

公式.jpg (20.08 KB, 下载次数: 1) 下载附件 保存到相册 2013-4-5 14:03 上传

一、在Abaqus/Standard 和Abaqus/Explicit中，输出的数据中的应变分量均是工程应变：

公式2.jpg (6.37 KB, 下载次数: 0) 下载附件 保存到相册 2013-4-5 14:06 上传

1111.jpg (10.83 KB, 下载次数: 0) 下载附件 保存到相册 2013-4-5 14:03 上传

二、在VUAMT中传入的应变增量strainInc的各个分量为：

公式3.jpg (3.12 KB, 下载次数: 0) 下载附件 保存到相册 2013-4-5 14:11 上传

，而不是工程应变

公式4.jpg (3.11 KB, 下载次数: 2) 下载附件 保存到相册 2013-4-5 14:11 上传

，所以在进行应力更新的时候应为：

EMU = E/(2.0*(1.0+ENU)) TWOMU = E/ (1.0+ENU) BULK3 = E/(1.0-TWO*ENU) LAM = TWOMU*BULK3*ENU/E

TRACE = STRAININC(NP,1) + S