问题描述:
算一个粘弹塑性的UMAT,在进入塑性之后,计算一定增量步后,DSTRAN显示为NaN,之后应力应变计算结果均为NaN.
解决办法:
这个问题我一开始也很蒙,原先想的是DSTRAN是abaqus给的,怎么会是NaN呢?
后来看了一些资料,自己也思考了一下,认为DSTRAN并不是abaqus给的这么简单,我们这里去想一下DSTRAN是怎么样算出来的?
DSTRAN实际上是求解红框中的方程算出来的,那么这样的话,DSTRAN之所以会变成NaN,就要去两个地方找原因,要么是DDSDDE雅克比矩阵出了问题,要么是单元应力F出了问题,有了方向之后,去相应出现问题的迭代步去一步步调试看计算输出,便可以比较方便找到问题。
另外,有很多奇形怪状的问题都是由于计算步长过大引起的,建议减小步长试试。
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删