C     混凝土(水灰比为0.3)MPS徐变模型UMAT子程序
C      
C
      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
     1 RPL,DDSDDT,DRPLDE,DRPLDT,
     2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
     3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
     4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME
      DIMENSION 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),DFGRD0(3,3),DFGRD1(3,3)   
C 
C    模型材料参数 
      EQ1=6.7E-12                     !q1--弹性模量的倒数
      EQ2=34.23E-12                   !q2
      EQ3=49.2E-12                    !q3
      EQ4=4.5E-12                     !q4
      TILO=1.0                        !加载时刻的时间
      ENU=0.2                         !泊松比
      TAU1=0.00001                    !第1个开尔文单元延迟时间
      TAU2=0.0001                     !第2个开尔文单元延迟时间
      TAU3=0.001                      !第3个开尔文单元延迟时间
      TAU4=0.01                       !第4个开尔文单元延迟时间
      TAU5=0.1                        !第5个开尔文单元延迟时间
      TAU6=1                          !第6个开尔文单元延迟时间
      TAU7=10                         !第7个开尔文单元延迟时间
      TAU8=100                        !第8个开尔文单元延迟时间
      TAU9=1000                       !第9个开尔文单元延迟时间
      TAU10=10000                     !第10个开尔文单元延迟时间
      A0=0.2794*EQ2                   !第0个开尔文单元参数A
      A1=0.0535*EQ2                   !第1个开尔文单元参数A       
      A2=0.0636*EQ2                   !第2个开尔文单元参数A      
      A3=0.0747*EQ2                   !第3个开尔文单元参数A       
      A4=0.0868*EQ2                   !第4个开尔文单元参数A       
      A5=0.0996*EQ2                   !第5个开尔文单元参数A      
      A6=0.1128*EQ2                   !第6个开尔文单元参数A   
      A7=0.1261*EQ2                   !第7个开尔文单元参数A       
      A8=0.1391*EQ2                   !第8个开尔文单元参数A       
      A9=0.1514*EQ2                   !第9个开尔文单元参数A      
      A10=0.1629*EQ2                  !第10个开尔文单元参数A
      EC=0.7E-19
      EK1=2.5E6
      EC0=2*EC/EQ4      
C 
      call dll_UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
     1 RPL,DDSDDT,DRPLDE,DRPLDT,
     2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
     3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
     4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC, 
     5 EQ1,EQ2,EQ3,EQ4,TILO,ENU,TAU1,TAU2,TAU3,TAU4,TAU5,,TAU6,
     6 TAU7,TAU8,TAU9,TAU10,A0,A1,A2,A3,A4,A5,A6,A7,A8,A9,A10,
     7 EC,EK1,EC0)
                     
      RETURN
      END  
C      
