VUMAT与UMAT相互转换注意事项

1、VUMAT与UMAT相互转换注意事项问题描述：最近在学习和编写子程序(V)UMAT的过程中，发现 VUMAT 和UMAT程序中应变增量的含意有不同。这个不同点对于需要在 VUMAT 和UMAT之间做转换的人来说容易搞错，虽然很基础，不清楚也容易搞 死人，所以我在本帖写个心得，提醒朋友们注意后处理输出的应变、UMAT 的应变和VUMA的应变有如下差别理论准备：根据弹性理论，应力应变关系如下:51 = 十 2Gsn ; ”耳= 2G&22 ; b# =+ 2Gs33=2G甩=Grn ； cr；3其中G是剪切模量,FROV www.Siri'bWecom公式.jpg (20.08 K

2、B,下载次数:1)下载附件保存到相册2013-4-5 14:03 上传一、在Abaqus/Standard禾口 Abaqus/Explicit中，输出的数据中的应变分量均是工程应变：匚成缶迟12, £1 玄= 11? 22? 33Z135/23公式2.jpg （6.37 KB,下载次数：0）下载附件保存到相册2013-4-5 14:06 上传 E: Strain componentsMax. PrincipalMid, Principal n Min. PrincipalEll E22E33/' E12 -,BEURcImjVeuaifM nerrts1111.jpg （10

3、.83 KB,下载次数：0）下载附件保存到相册2013-4-5 14:03 上传的各个分量为:二、在VUAM中传入的应变增量 strainine巧“ £黑養®空戸电? 23 公式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

4、= 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) =

5、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和VUAM的23方向的应变，正好呈现两倍的关系：LkI.E-8 D.VUMT中石方向应变www. si：nW丁即已2222.jpg（18.55 KB, 下载次数：0）、E23D1Xsow下载附件保存到相册2013-4-5 14:03 上传三、在UAM

6、中传入的应变增量DSTRAI的各个分量却为工程应变，与VUAM不同，所以在进行应力更新的时候，应为：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)+ EMU * DSTRAN(4)+ EMU * DSTRAN(5)+ EMU * DSTRAN(6)STRESS (4) = STRESS (4)STRESS (5) = STRESS (5)STRESS (6) = STRESS (6)或写成如下提供的UMAT程序验证如下：输出了 E.E23和UAMT勺23方向的应变，正好重合:D. 4 :MDEMBD4D