ABAQUS模拟梁单元断裂的本构方程（VUMAT）
* TEMPNEW, STRETCHNEW, DEFGRADNEW, FIELDNEW,
C WRITE ONLY -
* STRESSNEW, STATENEW, ENERINTERNNEW, ENERINELASNEW )
C
INCLUDE 'VABA_PARAM.INC'
C
DIMENSION COORDMP(NBLOCK,*), CHARLENGTH(NBLOCK), PROPS(NPROPS),
1 DENSITY(NBLOCK), STRAININC(NBLOCK,NDIR+NSHR),
2 RELSPININC(NBLOCK,NSHR), TEMPOLD(NBLOCK),
3 STRETCHOLD(NBLOCK,NDIR+NSHR),
4 DEFGRADOLD(NBLOCK,NDIR+NSHR+NSHR),
5 FIELDOLD(NBLOCK,NFIELDV), STRESSOLD(NBLOCK,NDIR+NSHR),
6 STATEOLD(NBLOCK,NSTATEV), ENERINTERNOLD(NBLOCK),
7 ENERINELASOLD(NBLOCK), TEMPNEW(NBLOCK),
8 STRETCHNEW(NBLOCK,NDIR+NSHR),
9 DEFGRADNEW(NBLOCK,NDIR+NSHR+NSHR),
1 FIELDNEW(NBLOCK,NFIELDV),
2 STRESSNEW(NBLOCK,NDIR+NSHR), STATENEW(NBLOCK,NSTATEV),
3 ENERINTERNNEW(NBLOCK), ENERINELASNEW(NBLOCK)
C
CHARACTER*80 CMNAME
DIMENSION INTV(2)
PARAMETER ( ZERO = 0.D0, ONE = 1.D0, TWO = 2.D0, THREE = 3.D0,
* THIRD = ONE / THREE, HALF = 0.5D0, TWOTHDS = TWO / THREE,
* OP5 = 1.5D0 )
C
C THE STATE VARIABLES ARE STORED AS:
C STATE(*,1) = ELEMENT DELETE
*
*
open(100,file='E:\Temp_ABAQUS\model\00.txt')
F = PROPS(1) 【本构参数】
G = PROPS(2) 【本构参数】
H = PROPS(3) 【本构参数】
T = PROPS(4) 【判据参数】
*
*
*
DO K = 1, NBLOCK
C
STATENEW(K,2)=STATEOLD(K,2)+STRAININC(K,1)
IF(STATEOLD(K,1).EQ.ONE)THEN
STRESSNEW(K,1)=F*(exp(H*STATENEW(K,2))
1 -exp(G*STATENEW(K,2))) 【材料本构关系式，计算应力】
IF(STRESSNEW(K,1).GE.T)THEN 【应力判定】
STATENEW(K,1)=ZERO
write(100,*) 'K=',K 【写入数据】
write(100,*) 'steptime=',STEPTIME 【写入数据】
write(100,*) ' stress-crack=',STRESSNEW(K,1) 【写入数据】
ELSE
STATENEW(K,1)=ONE
STRESSPOWER = HALF *
1 ( ( STRESSOLD(K,1)+STRESSNEW(K,1) )*STRAININC(K,1) ) 【更新应力】
C
C
ENERINTERNNEW(K) = ENERINTERNOLD(K)
1 + STRESSPOWER / DENSITY(K) 【更新内能】
END IF
ELSE
STATENEW(K,1)=ZERO
END IF
END DO
RETURN
END
