---------------------------------------------------------------
C 瞬态温度场下修正 Burgers 模型 UMAT 子程序源代码
C 给状态变量数组赋初值为零,调用 ABAQUS 子程序 SDVINI
C GIVE STATEV THE INITIAL VALUE OF ZERO
SUBROUTINE SDVINI(STATEV,COORDS,NSTATV,NCRDS,NOEL,NPT,LAYER,KSPT)
INCLUDE 'ABA_PARAM.INC'
DIMENSION STATEV(NSTATV),COORDS(NCRDS)
DO K=1,NSTATV
STATEV(K)=0.0
END DO
RETURN
END
C 瞬态温度场下修正 Burgers 模型 UMAT 子程序
C UMAT FOR MODIFIED BURGERS MODEL
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,
1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,
2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,
3 PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
INCLUDE 'ABA_PARAM.INC'
CHARACTER*80 CMNAME
DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),
1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),
2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),
3 DFGRD0(3,3),DFGRD1(3,3)
C
C LOCAL ARRAYS(定义局部数组,即用户自己定义的数组)
C ----------------------------------------------------------------------------------------------------------------
C EELAS - ELASTIC STRAINS
C DEELA - ELASTIC STRAINS INCREMENT(dt)
C ECREE - CREEP STRAINS
C DECRE - CREEP STRAINS INCREMENT(dt)
C STREST - TEMPORARY ARRAY FOR SAVED STRESS(t+dt)
C DECRT - TEMPORARY ARRAY FOR SAVED CREEP STRAINS INCREMENT(dt)
C DSTREST- STRESS INCREMENT(dt)
C STREST2- TEMPORARY ARRAY FOR SAVED STRESS(t+θdt)
C DVSTRESS-DEVIATORIC STRESS(t+θdt)
C PARAM - ARRAY FOR SAVED MATERIAL PARAMETERS AT TEMPERATURE
C OF THE KSTEP
C ----------------------------------------------------------------------------------------------------------------
DIMENSION EELAS(6),ECREE(6),DECRE(6),DVSTRESS(6),STREST(6),
1 DECRT(6),DEELA(6),DSTREST(6),STREST2(6),PARAM(7),
2 TABLE(NPROPS/7,7),TABLE0(7),TABLE1(7)
C
评论0