卢剑锋umat子程序的疑问，对JOhn-cook公式进行求导？

woaigechang

而在卢剑锋的子程序中，
SYIELD=SYIEL0
          DEQPL=(SMISES-SYIELD)/EG3
          DSTRES=TOLER*SYIEL0/EG3
          DEQMIN=HALF*DTIME*EXP(1.0D-4/C)                           
          DO 130 KEWTON=1,NEWTON
               DEQPL=MAX(DEQPL,DEQMIN)
               CALL USERHARD(SYIELD,HARD,EQPLAS+DEQPL,PROPS(4))
               TVP=C*LOG(DEQPL/DTIME)
               TVP1=TVP+ONE
               HARD1=HARD*TVP1+SYIELD*C/DEQPL 这个就应该是h了吧，求导好像不太对啊
               SYIELD=SYIELD*TVP1
               RHS=SMISES-EG3*DEQPL-SYIELD
               DEQPL=DEQPL+RHS/(EG3+HARD1)
          IF(ABS(RHS/EG3) .LE. DSTRES ) GOTO 140


根据后面userhead的子程序，

woaigechang

求斑竹来回答

woaigechang

自己顶啊，这个问题这么重要，不能没人回答啊

songhengxu

你看求导是上面那个式子求导么？userhard里面的HARD和SYIELD是这个表达式么？

woaigechang

songhengxu 发表于 2011-11-7 13:48
你看求导是上面那个式子求导么？userhard里面的HARD和SYIELD是这个表达式么？

是啊，难道我看错了？

蓝牙

我见过这个umat 但是没有深究过，我提一个问题也许对你有帮助
hard应该是等效应力的增量和等效塑性应变的增量的比值，不是等效应力的增量和等效应变的增量的比值，这个要用牛顿迭代的方式实现！如果楼主注意到这个问题，请略过！

billowriver

把应变率按照差分格式展开以后再求导。应变率是应变的函数，这一项的影响不能忽略。还有，弱问一下，卢剑锋是谁

woaigechang

billowriver 发表于 2011-11-8 13:24
把应变率按照差分格式展开以后再求导。应变率是应变的函数，这一项的影响不能忽略。还有，弱问一下，卢剑锋 ...

现在看的是卢剑锋的UMAT程序。论坛里能下到。
怎样按照差分格式展开？能说的具体一点吗

billowriver

向前差分，E^dot=(E^k-E^k-1)/dt，屈服函数也是E^k的函数，整个表达式对E^k求导，

bhx_gs

不知问题是否解决？
各人认为求导没有问题！
因为应变速率等于应变增量（deqpl）除以时间增量（dtime）

只是很难收敛总是提示too many increments needed to complete the step,不知为何？

tjuabraham

JC模型其实有两种形式：（1）直接将应变率效应视为一个屈服应力的放大scale，此时应变率项是等效应变（2）使用VISO-PLASTIC即粘塑性形式，此时应变率项是等效塑性应变。这时候求导应该就能看懂了。
上面的结论是从DYNA keyword手册里面看到的。再填写card的时候dyna中有vp项进行选择。不知道abaqus和dyna是不是一样

wangcongkang

改程序没有考虑温度软化部分 ，我将temp变量根据J-C公式写入程序，发现仿真结果例如应力和力都小于软件自带程序计算所得结果，请问大家有无遇到这个问题。

小小书虫

wangcongkang 发表于 2012-3-3 09:57
改程序没有考虑温度软化部分 ，我将temp变量根据J-C公式写入程序，发现仿真结果例如应力和力都小于软件自带 ...

抱着试一试的态度，不知道楼主还登陆论坛不？如果看到了，希望能交流一下怎么把温度项加进去。。。。谢谢楼主？

good1990

感觉只是在硬化那里用了吧

littlehuster

感觉论坛里没什么人了，楼主还在吗，我写了对这个问题的看法，文章链接是https://blog.csdn.net/weixin_42887395/article/details/107882614