首先你要知道分析的问题NDI,NSHR和NTENS的值,然后根据相应的值来编写UMAT,所以很多例子中都有一句判断: IF (NDI.NE.3) THEN
WRITE(6,1)
1 FORMAT(//,30X,'***ERROR - THIS UMAT MAY ONLY BE USED FOR ',
1 'ELEMENTS WITH THREE DIRECT STRESS COMPONENTS')
ENDIF
如果能统一,应该就可以去掉判断了,这取决于你的本构模型和相应的算法.
第一个问题就是1994的那本《塑性与建模》的“加载状态的确定》,照斑竹这么说那个r的确定是隐式本构积分算法了?那么大牛Tedbelyscle的那本非线性有限元书上的算法也是隐式算法?而且他的算法就是应力驱动算法?
第二个问题,不知道版主的本构是三维的,还是二维的?
ABAQUS对于平面应变和轴对称的处理,好像不是一般有限元书上讲的,对于平面应变,ABAQUS 还是按照3个正应力1个剪应力来做的,当然轴对称是4个应力分量(不过z方向变成了2方向,而不是三维问题中的3方向),个人认为对于轴对称不需作特殊处理,只需记住2方向的不同即可;但是平面应变问题,按照ABAQUS ,3方向是z方向,在这个方向上应变始终为0,那么在umat中怎么应力更新,得到非0的应力?在一般有限元书上,E,V都要作一下转换的(对于平面应变)。
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删