billowriver 发表于 2006-8-30 22:48

求教VUMAT如何定义材料失效

最近用VUMAT处理断裂问题，本构程序用的就是ABAQUS例子自带的2D线性硬化模型，只是
加了处理断裂的程序，按照document里的说法，failure state variable- 1==not
failed, 0==failed。所以.f里加的程序：
         if ( stateNEW(i,5).LT.0.5)then
         statenew(i,6)=1
         else
         statenew(i,6)=0
         endif         
*其中StateNEW(i,5)表示的是等效塑性应变(这个从不加断裂准则输出结果来看是没问题
的), .INP文件里加的程序
   *MATERIAL,NAME=Alu
   *USER MATERIAL,CONSTANTS=4
   68.66E9,.3,200E6,0E6
   *DENSITY
   2700,
   *DEPVAR,DELETE=6
   6,
但是计算结果显示材料点并没有被删除，虽然当等效应变大于0.5的时候单元的应力都变
为0但是当应变小于0.5的时候，单元又重新承载,与直接用ABAQUS自带的失效模型给出的
结果相差很大，请教这样是什么原因造成的，是不是VUMAT用的不对。

[ 本帖最后由 billowriver 于 2006-8-31 16:14 编辑 ]

mafia00 发表于 2006-9-1 10:44

应变小于0.5的时候程序里面又重新进行判断，结果删除标记置为1，就是不删除，所以单元又恢复了。如果不想让它恢复就检查一下stateold的状态，如果是0的话statenew就强制为0

billowriver 发表于 2006-9-1 15:57

原帖由 mafia00 于 2006-9-1 10:44 发表
应变小于0.5的时候程序里面又重新进行判断，结果删除标记置为1，就是不删除，所以单元又恢复了。如果不想让它恢复就检查一下stateold的状态，如果是0的话statenew就强制为0


谢谢你的答复，把vumat里的程序改了一下，但是问题依然存在，主要是单元无法断开，一直连在一起，是不是vumat的用法上的问题呢？
                     if(stateNew(i,5).LT.0.5)then
             if(stateOld(i,6).EQ.0)then
             stateNew(i,6)=0
             else
             stateNew(i,6)=1
             endif
          else
          stateNew(i,6)=0
          endif

billowriver 发表于 2007-3-7 00:44

由于前段时间忙其他的事情，这个问题现在还没有解决，希望各位用过VUMAT的高手能给解答一下

liupenghuild 发表于 2011-10-30 22:17

您好！我现在也在做关于材料的实效问题，在vumat子程序中导入实效准则。但现在我不知道改怎么编译实效准则，您能把你关于这个算例的vumat子程序发我一份让我参考下嘛？我的邮箱是312512518@163.com。谢谢！

niweiyaa 发表于 2012-12-5 16:24

liupenghuild 发表于 2011-10-30 22:17 static/image/common/back.gif
您好！我现在也在做关于材料的实效问题，在vumat子程序中导入实效准则。但现在我不知道改怎么编译实效准则， ...

您好，问下你做出来了吗？上面的程序对不？

liupenghuild 发表于 2012-12-6 19:18

niweiyaa 发表于 2012-12-5 16:24 static/image/common/back.gif
您好，问下你做出来了吗？上面的程序对不？

您好，我们可以QQ上交流。我的QQ号是312512518.

Lorraine 发表于 2012-12-10 11:19

一份相关的资料材料破坏的算法，这里采用了B.P. Gearing, L. Anand / International Journal of Solids and Structures 41 (2004) 3125–3150中推荐的方法。原理很简单:
（1）根据应力状态，计算主应力。如果hydrostatic stress为正（就是受拉了），最大主应力为正且大于一个值（材料参数），破坏开始形成（即initiation），并在最大主应力方向上产生一个塑性变形，文章里给出了一个公式计算塑性变形的大小（记为P1），当然自己也可以提出合适的公式。否则剪切变形占主导，采用一般的本构方程计算应力应变(等效塑性应变记为P2)。
（2）当p1_cr或p2_cr时（p1_cr,p2_cr均为材料参数），积分点失效，在VUMAT中通知ABAQUS删除就行了（这是原方法），当一个单元内所有的积分点失效后，该单元自动失效。
（3）为了计算稳定，并且更符合实际破坏的过程，我做了一些修改，不是在p1_cr或p2_cr时立即将单元删除，而是定义了一个我自己称之为质量因子（q_factor）的自变量，当p1_cr或p2_cr后，q_factor很快衰减，材料的模量等抵抗变形的参数也随之快速衰减，直至没有能力承受载荷，然后删除积分点
（4）几点说明
(a)破坏的initiation,growth直至材料breakdown，基本上没有统一的标准，所以可以根据实际材料的性能自己提出标准；
(b)这种failure model只能在abaqus/explicit中使用，所以计算效率可能不令人满意；

techique 发表于 2014-5-21 11:09

谢谢，下载看看，我正在学Abaqus，可能用得着。

vten1 发表于 2014-7-29 21:21

first log,interesting!

vten1 发表于 2014-7-29 21:21

伤脑筋啊

bhulele 发表于 2014-8-18 16:57

{:{10}:}{:{10}:}哈哈哈，我也做相关的研究

bhulele 发表于 2014-8-18 17:34

看看看，研究一下

kyletao 发表于 2014-9-1 09:29

谢谢，下载看看，我正在学Abaqus
页: [1]
查看完整版本: 求教VUMAT如何定义材料失效