abaqus1用户材料子程序(8)

STATEV(1+2*NTENS)=EQPLAS C

RETURN END C C C

SUBROUTINE USERHARD(SYIELD,HARD,EQPLAS,TABLE) C

INCLUDE 'ABA_PARAM.INC' C

DIMENSION TABLE(3) C

C GET PARAMETERS, SET HARDENING TO ZERO C

A=TABLE(1) B=TABLE(2) EN=TABLE(3) HARD=0.0 C

C CALSULATE CURRENT YIELD STRESS AND HARDENING RATE C

IF(EQPLAS.EQ.0.0) THEN SYIELD=A ELSE

HARD=EN*B*EQPLAS**(EN-1) SYIELD=A+B*EQPLAS**EN END IF RETURN END

18.5.2 UMATHT（包含材料的热行为）

SUBROUTINE UMATHT(U,DUDT,DUDG,FLUX,DFDT,DFDG,STATEV,TEMP,

$ DTEMP,DTEMDX,TIME,DTIME,PREDEF,DPRED,CMNAME,NTGRD,NSTATV,

$ PROPS,NPROPS,COORDS,PNEWDT,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) C

INCLUDE 'ABA_PARAM.INC'

19-36

C

CHARACTER*80 CMNAME C

DIMENSION DUDG(NTGRD),FLUX(NTGRD),DFDT(NTGRD),

$ DFDG(NTGRD,NTGRD),STATEV(NSTATV),DTEMDX(NTGRD),TIME(2), $ PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3) C C

COND=PROPS(1) SPECHT=PROPS(2) C

C INPUT SPECIFIC HEAT DUDT=SPECHT DU=DUDT*DTEMP U=U+DU C

C INPUT FLUX = -[K]*{DTEMDX} DO I=1, NTGRD

FLUX(I)=-COND*DTEMDX(I) END DO C

C INPUT ISOTROPIC CONDUCTIVITY C

DO I=1, NTGRD

DFDG(I,I)=-COND END DO C

RETURN END

本章主要内容引自：

卢剑锋，冲击载荷作用下材料和结构力学行为有限元模拟，清华大学硕士论文，2003。

19-37