zwt umat子程序 一个非线性弹簧并联4个maxwell体

SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,

     1 RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN,

     2 TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,MATERL,NDI,NSHR,NTENS,

     3 NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT,

     4 DFGRD0,DFGRD1,NOEL,NPT,KSLAY,KSPT,KSTEP,KINC)

      INCLUDE 'ABA_PARAM.INC'

      CHARACTER*80 MATERL

       real*8 STRESS(NTENS),STATEV(NSTATV),

     1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),

     2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),

     3 PROPS(NPROPS),COORDS(3),DROT(3,3),

     4 DFGRD0(3,3),DFGRD1(3,3),PIAN(6),PIAN1(6),PIAN2(6),PIAN3(6),

     5 SHUZU(NTENS)

      !props(1)=E0,props(2)=v,泊松比，props(3)=a,props(4)=b,!非线性弹簧参数与泊松比

         ! 4个Maxwell体参数。props(5)=E1,props(6)=o1, props(7)=E2,props(8)=o2,

      ! props(9)=E3,props(10)=o3, props(11)=E4,props(12)=o4

      REAL*8 TONGSHI