刚哥经典UMAT求释疑

volcanozhang

刚哥的UMAT里面有这么一段代码：
         DO J=1,NSLPTL

            DDEMSD(4,J)=DDEMSD(4,J)-SLPSPN(1,J)*STRESS(1)
            DDEMSD(5,J)=DDEMSD(5,J)+SLPSPN(2,J)*STRESS(1)

            IF (NDI.GT.1) THEN
               DDEMSD(4,J)=DDEMSD(4,J)+SLPSPN(1,J)*STRESS(2)
               DDEMSD(6,J)=DDEMSD(6,J)-SLPSPN(3,J)*STRESS(2)
            END IF

            IF (NDI.GT.2) THEN
               DDEMSD(5,J)=DDEMSD(5,J)-SLPSPN(2,J)*STRESS(3)
               DDEMSD(6,J)=DDEMSD(6,J)+SLPSPN(3,J)*STRESS(3)
            END IF

            IF (NSHR.GE.1) THEN
               DDEMSD(1,J)=DDEMSD(1,J)+SLPSPN(1,J)*STRESS(NDI+1)
               DDEMSD(2,J)=DDEMSD(2,J)-SLPSPN(1,J)*STRESS(NDI+1)
               DDEMSD(5,J)=DDEMSD(5,J)-SLPSPN(3,J)*STRESS(NDI+1)
               DDEMSD(6,J)=DDEMSD(6,J)+SLPSPN(2,J)*STRESS(NDI+1)
            END IF

            IF (NSHR.GE.2) THEN
               DDEMSD(1,J)=DDEMSD(1,J)-SLPSPN(2,J)*STRESS(NDI+2)
               DDEMSD(3,J)=DDEMSD(3,J)+SLPSPN(2,J)*STRESS(NDI+2)
               DDEMSD(4,J)=DDEMSD(4,J)+SLPSPN(3,J)*STRESS(NDI+2)
               DDEMSD(6,J)=DDEMSD(6,J)-SLPSPN(1,J)*STRESS(NDI+2)
            END IF

            IF (NSHR.EQ.3) THEN
               DDEMSD(2,J)=DDEMSD(2,J)+SLPSPN(3,J)*STRESS(NDI+3)
               DDEMSD(3,J)=DDEMSD(3,J)-SLPSPN(3,J)*STRESS(NDI+3)
               DDEMSD(4,J)=DDEMSD(4,J)-SLPSPN(2,J)*STRESS(NDI+3)
               DDEMSD(5,J)=DDEMSD(5,J)+SLPSPN(1,J)*STRESS(NDI+3)
            END IF

         END DO
这段代码的意思应该是求L+σ*W+W*σ（1）。
如果我没看错的话，DDEMSD(1,J)=DDEMSD(1,J)+SLPSPN(1,J)*STRESS(NDI+1)-SLPSPN(2,J)*STRESS(NDI+2)
DDEMSD(2,J)=DDEMSD(2,J)-SLPSPN(1,J)*STRESS(NDI+1)+SLPSPN(3,J)*STRESS(NDI+3)
DDEMSD(3,J)=DDEMSD(3,J)+SLPSPN(2,J)*STRESS(NDI+2)-SLPSPN(3,J)*STRESS(NDI+3)
但是从式（1）来看，后面的值应该乘一个2才对，也就是说例如DDEMSD(1,J)，应该是
DDEMSD(1,J)=DDEMSD(1,J)+2*SLPSPN(1,J)*STRESS(NDI+1)-2*SLPSPN(2,J)*STRESS(NDI+2)
但不知道为什么刚哥没有乘2，谢谢。呵呵。

volcanozhang

自己顶一下

陌陌

刚哥是谁？

fr_zengtao

近来我也在看，没有完全看懂，如果可以 能否交流交流

仔细看了一些你说的那个问题，似乎应该要乘以2.，不明白为什么不乘

volcanozhang

再顶三次，如果没人回答就让它沉底吧。
哪位哥们知道如果不方便写，就留个电话吧，电话请教。在家靠父母在外靠朋友，多谢了。

volcanozhang

现在很多人在讨论刚哥的程序。搭这个便车，再顶两次。

monolithlee

σ*W+W*σ   ?=    2W*σ
黄勇刚的程序我没仔细研究过，我考虑是不是这里的问题

volcanozhang

你说的等式成立是有前提的，前提是σ*W是对称张量
7# monolithlee

watsonminer

来学习的，楼主的问题解决了么~~~求解答~~

ccy7270801

谁能提供以下程序？没有找到啊

volcanozhang

http://forum.simwe.com/viewthread.php?tid=765756&highlight=%BE%A7%CC%E5%CB%DC%D0%D4
10# ccy7270801

billowriver

L+σ*W+W*σ（1）那个笑脸是什么东西，如果是求Cauchy 应力的Jaumann Rate的话为什么要乘以2，楼主的出发点是什么，DDEMSD又是什么

yuhan0172009

有没有人解答啊，我也想知道

pangshaobing

有没有人解答啊

liliu328

具体内容没有看，不过ABAQUS里面确实有一个量是两倍的关系：失量写法中存储的是工程剪应变=2*张量中存储的是剪应变，不知道是不是这个问题。