心得 对VUMAT与UMAT相互转换中应注意的应变含义的探讨

sr楚风

cuba20052001 发表于 2014-2-16 01:57
本构调用VUMAT子程序，odb文件中的剪切应变应该还是工程剪切应变吧

请问您这的句活能肯定吗？

princewhb

谢谢楼主的解答。

taishanbuzuo

借宝地补充一点：
除了2倍关系之外，UMAT和VUMAT中应力  应变存储顺序是不一样的，
以应力为例，
UMAT中存储顺序：S11，S22，S33，S12，S13，S23
VUMAT中存储顺序：S11，S22，S33，S12, S23, S13

因此，在分别编写两种子程序时，需要注意两者的区别。

112233a

关于abaqus接触参数的设置abaqus本身就是库伦模型么 还是要重新设置编程呢

chuguangxiang

学习中

chuguangxiang

amani 发表于 2013-10-1 11:26
这个问题我来说两句。在力学领域，理论推导一般是用张量，其非常优雅而简洁，但是矩阵数值计算只认可数字 ...

您好，我最近访问论坛，看到您回复过很多信息，我看庄茁那篇文章做的效果，我没做出来，总会出现不收敛，不知道为什么？由于本人等级较低，不能加您为好友，所以只能通过此种方式和您交流

snowwave02

在UMAT中abaqus的剪切应变用的是工程应变gamm12=eps12+eps21,Jacobian的理论公式是正常的剪切应变eps12，那么到UMAT中转换为工程剪切应变gamm12后，Jacobian矩阵在剪切项的元素都应该除以2才正常。
但abaqus自带的塑性UMAT的for文件是在求前面部分时正常除了2，但后面部分没有除以2，为什么？
        do k1=1, ndi
          do k2=1, ndi
            ddsdde(k2, k1)=efflam
          end do
          ddsdde(k1, k1)=effg2+efflam
        end do
        do k1=ndi+1, ntens
          ddsdde(k1, k1)=effg  （说明：和直接分量ndi相比除以2，正常）
        end do
        do k1=1, ntens
          do k2=1, ntens
            ddsdde(k2, k1)=ddsdde(k2, k1)+effhrd*flow(k2)*flow(k1)（说明：没有将剪切分量除以2,为什么？）
          end do
        end do

慌张的小铁匠

版主太给力了

枯藤老树昏鸦

学习一下