VUMAT错误：zero or negative initial dilatational modulus

wzhswjtu

用VUMAT时，程序报错： zero or negative initial dilatational modulus caused by bad material data. Please check your material input and any initial conditions if necessary.
搜索半天，发现这个问题很多人遇到，但无解决方法。

首先，什么是dilatational modulus ？膨胀模量？如何计算？
查了abaqus帮助，这个词未出现过，只在2.4.5 Explicit dynamic analysis中有dilatational wave，说明，此模量与材料膨胀波速有关，波速与材料弹性模量及密度有关，但我的程序中这个参数是正常的。

此错误是否与稳定时间步长有关呢？稳定时间步长如何输出呢？
加大了密度1e7倍，增大稳定时间步长，无效！

这个错误在check的时候出现，然后输出了过程中各个变量，程序给定的check的strainInc为-0.001，最终计算的stressNew是正常的。

请大家关注，谢谢！

wzhswjtu

1# wzhswjtu
查到了一点信息：Dilatational modulus is the Bulk modulus or the ) modulus of incompressibility. If you get an error of 7 \"Bad Material definition: zero or negative ..." that 6 w0 dsimply means that your material data is problematic. 0 If you get it for a VUMAT your material response as a ; bconsequence of bad material data is unphysical. 5 {. B There are mathematical restrictions on how " C7 |. the values of E, G, and nu relate to each other.  For isotropic
materials, the two equations are:  G=E/(2*(1+nu)) and -1 < nu < 1/2.  " G1 _4 o( F1 Y6 d* u3 N) @For orthotropic/anisotropic material they get much more complex.
若Dilational modulus是体积模量，体积模量公式为：

程序中输入弹性模量为32173000000Pa，nu=0.18，计算出K=16756770833，不会为负。

wzhswjtu

1# wzhswjtu
部分人说：ABAQUS在检查VUMAT时，试探地施加一个很小的变形，然后检查相应的应力。如果根据你的VUMAT计算得到的应力为0，就得到上面的错误消息。所以，可以检查一下你的应以更新算法，或者稍微更改一下，应该就可以消除这个错误。


在check过程中应变增量为：（-0.001，-0.001，-0.001，0.001，0.001，0.001）


应力计算结果为：（-103214603.568，-103214603.568，-103214604.093，25154255.573，25154255.628，25154255.628）

wzhswjtu

3# wzhswjtu
在“VUMAT Initial calculations and checks”中说：
In the data check phase of the analysis Abaqus/Explicit calls user subroutine VUMAT with a set of fictitious strains and a totalTime and stepTime both equal to 0.0. This is done as a check on your constitutive relation and to calculate the equivalent initial material properties, based upon which the initial elastic wave speeds are computed.


翻译一下：在explicit的模型信息检查阶段内，调用子程序VUMAT时，施加一虚应变，totalTime和stepTime均为0。这是为了检查用户本构关系；计算材料的初始等效属性，基于此，计算初始弹性波速。


是不是可认为此错误的产生仍然是与弹性波速有关，庄老师的书中给出了
对于泊松比为0的线弹性材料，波速cd=sqrt（E/density）
那对于泊松比不为0的材料，波速如何计算呢？

wzhswjtu

4# wzhswjtu
问题解决了，但是原因还是不明。

我在子程序的最后添加了
      do 70 j=1,6
        stressOld(i,j) = stressNew(i,j)
  70  continue
便会出现这样的error

我看了其他的例子，stress是不需要像上面的语句那样来赋值的，那么下一迭代步或荷载步的stressOld是如何得到的呢？在主程序中实现的吗？

采用了上面的赋值方式，为什么会牵涉到材料属性的问题呢？

wzhswjtu

5# wzhswjtu

试验证明，stressOld不需要在子程序最后进行赋值更新。

lxzzwn1004

wzhswjtu 发表于 2011-2-20 12:07
5# wzhswjtu  

试验证明，stressOld不需要在子程序最后进行赋值更新。

那楼主，最后不应该这么对old赋值，那你的问题最终是这么解决的？
我现在也遇到了同样的问题！

猪公快跑

wzhswjtu 发表于 2011-2-20 12:07
5# wzhswjtu  

试验证明，stressOld不需要在子程序最后进行赋值更新。

这一步的stressnew在下一步的再传进来的时候就变成了stressold了，所以不需要赋值

猪公快跑

wzhswjtu 发表于 2011-2-20 10:54
4# wzhswjtu  
问题解决了，但是原因还是不明。

楼主的这个问题是怎么改正通过的，我用的是大变形的本构关系，需要有什么特殊考虑吗？

MENGZHEN

wzhswjtu 发表于 2011-2-20 10:54
4# wzhswjtu  
问题解决了，但是原因还是不明。

同问，楼主问题是怎么解决的？

陌丶上

没有那么麻烦，只是子程序开头需要个判断steptime=0时的应力更新，补上就可以了

zyhnwpu

陌丶上 发表于 2015-2-6 21:01
没有那么麻烦，只是子程序开头需要个判断steptime=0时的应力更新，补上就可以了 ...

step=0时的应力更新是怎么定义的？

nuaa2013

wzhswjtu 发表于 2011-2-20 10:54
4# wzhswjtu  
问题解决了，但是原因还是不明。

请问楼主，这个报错是如何解决的呢？

nuaa2013

nuaa2013 发表于 2022-6-17 22:29
请问楼主，这个报错是如何解决的呢？

问题已解决，是子程序中材料参数的变量名字的问题，应该采用大写字母