VUMAT与UMAT相互转换注意事项.doc

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/ETRACE = STRAININC(NP,1) + STRAININC(NP,2) + STRAININC(NP,3)STRESSNEW(NP,1) = STRESSOLD(NP,1) + LAM*TRACE + TWOMU*STRAININC(NP,1)STRESSNEW(NP,2) = STRESSOLD(NP,2) + LAM*TRACE + TWOMU*STRAININC(NP,2)STRESSNEW(NP,3) = STRESSOLD(NP,3) + LAM*TRACE + TWOMU*STRAININC(NP,3)STRESSNEW(NP,4) = STRESSOLD(NP,4)? ?? ?? ?? ? + TWOMU*STRAININC(NP,4)STRESSNEW(NP,5) = STRESSOLD(NP,5)? ?? ?? ?? ? + TWOMU*STRAININC(NP,5)STRESSNEW(NP,6) = STRESSOLD(NP,6)? ?? ?? ?? ? + TWOMU*STRAININC(NP,6) 验证：我做了一个简单的弹性剪切的例子，输出了E.E23和VUAMT的23方向的应变，正好呈现两倍的关系： 2222.jpg (18.55 KB, 下载次数: 0) 下载附件 ?保存到相册 2013-4-5 14:03 上传 三、在UAMT中传入的应变增量DSTRAN的各个分量却为工程应变 ，与VUAMT不同，所以在进行应力更新的时候，应为：EMU? ? = E/(2.0*(1.0+ENU))TWOMU??= E/ (1.0+ENU)BULK3??= E/(1.0-TWO*ENU)LAM? ? = TWOMU*BULK3*ENU/ETRACE = DSTRAN(1) + DSTRAN (2) + DSTRAN (3)STRESS (1) = STRESS (1) + LAM*TRACE + TWOMU * DSTRAN(1)STRESS (2) = STRESS (2) + LAM*TRACE + TWOMU * DSTRAN(2)STRESS (3) = STRESS (3) + LAM*TRACE + TWOMU * DSTRAN(3)STRESS (4) = STRESS (4)? ?? ?? ?? ? + EMU * DSTRAN(4)STRESS (5) = STRESS (5)? ?? ?? ?? ? + EMU * DSTRAN(5)STRESS (6) = STRESS (6)? ?? ?? ?? ? + EMU * DSTRAN(6) 或写成如下提供的UMAT程序验证如下：输出了E.E23和UAMT的23方向的应变，正好重合：  3333.jpg (18.44 KB, 下载次数: 0) 下载附件 ?保存到相册 2013-4-5 14:03 上传 2013-4-5 14:28 上传 点击文件名下载附件