首先你要知道分析的問題NDI,NSHR和NTENS的值，然后根據(jù)相應(yīng)的值來編寫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

如果能統(tǒng)一，應(yīng)該就可以去掉判斷了，這取決于你的本構(gòu)模型和相應(yīng)的算法.

第一個(gè)問題就是1994的那本《塑性與建模》的“加載狀態(tài)的確定》，照斑竹這么說那個(gè)r的確定是隱式本構(gòu)積分算法了?那么大牛Tedbelyscle的那本非線性有限元書上的算法也是隱式算法？而且他的算法就是應(yīng)力驅(qū)動算法？

第二個(gè)問題，不知道版主的本構(gòu)是三維的，還是二維的?

ABAQUS對于平面應(yīng)變和軸對稱的處理，好像不是一般有限元書上講的，對于平面應(yīng)變，ABAQUS 還是按照3個(gè)正應(yīng)力1個(gè)剪應(yīng)力來做的，當(dāng)然軸對稱是4個(gè)應(yīng)力分量（不過z方向變成了2方向，而不是三維問題中的3方向），個(gè)人認(rèn)為對于軸對稱不需作特殊處理，只需記住2方向的不同即可；但是平面應(yīng)變問題，按照ABAQUS ，3方向是z方向，在這個(gè)方向上應(yīng)變始終為0，那么在umat中怎么應(yīng)力更新，得到非0的應(yīng)力？在一般有限元書上，E,V都要作一下轉(zhuǎn)換的（對于平面應(yīng)變）。