求高手指教关于子程序运行的问题Problem during linking 。不知道什么毛病。急啊。 .. ...

[程序写的是一个温度粘弹性材料的本构方程，在附件里。。



abaqus运行时显示一下信息：




Problem during linking - Abaqus/Standard User Subroutines.   This error may be due to a mismatch in the Abaqus user subroutine arguments.   These arguments sometimes change from release to release, so user subroutines   used with a previous release of Abaqus may need to be adjusted.




以下是程序内容：


      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),
     2 DDSDDT(NTENS),DRPLDE(NTENS),
     3 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),
     4 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)
      DIMENSION DSTRES(6),D(3,3)
C   CONST1=αE CONST2=αμ CONST3=αλ CONST4=αC CONST5=αε CONST6=α
C   EMOD=PROPS(1)=EG EVIS=PROPS(2)=μG EREL=PROPS(3)=λG EC=PROPS(4)=CG ESTRA=PROPS(5)=TEMPG=εG PROPS(6)=TG
C
      CONST1=38.1
      CONST2=44.2
      CONST3=35.4
      CONST4=38.7
      CONST5=58.2
      CONST6=11.6**-5
      CONST7=0.3333
C
      TEMPG=PROPS(6)
      EMOD=PROSP(1)
C
      EMOD=EMOD*EXP(CONST1*(TEMPG/TEMP-1))
      EVIS=PROSP(2)
C
      EVIS=EVIS*EXP(CONST3*(TEMPG/TEMP-1))
      EREL=PROPS(3)
C
      EREL=EREL*EXP(CONST3*(TEMPG/TEMP-1))
      ESTRA=PROPS(5)
C
      ESTRA=ESTRA*EXP(-CONST5*(TEMPG/TEMP-1))
      EC=PROPS(4)
C
      EC=EC*EXP(CONST4*(TEMPG/TEMP-1))
C
      TERM1=.1*(1+CONST7)
C
      TERM2=(1-EC)/EREL
C
      TERM3=CONST7/((1+CONST7)*(1-2*CONST7))
C
      TERM4=TERM1+.5*TERM2*DTIME/(1+CONST7)
C
      TERM5=TERM3+.5*TEMP2*TERM3*DTIME
C
      TEMP6=.5*EMOD+.5*DTIME/EVIS
C
C   EVALUATE NEW STRESS TENSER
C
      EV=0.
      DEV=0.
      DO K1=1,NDI
      EV=EV+STRAN(K1)
      DEV=DEV+DSTRAN(K1)
      END DO
C
      DO K1=1,NDI
      DSTRES(K1)=EMOD*(TERM1*DSTRAN(K1)+TERM3*DEV-
    1STRESS(K1)*DTIME/EVIS+(TERM2*DSTRAN(K1)/(1+CONST7)+
    2TERM2*TERM3*EV)*DTIME-EC*ESTRA*DTIME-CONST6*DTEMP)
      STRESS(K1)=STRESS(K1)+DSTRES(K1)
      END DO
C
C
C
      I1=NDI
      DO K1=1,3
       I1=I1+1
       DSTRES(I1)=EMOD*(DSTRAN(I1)-STRESS(I1)*DTIME/EVIS+
    1TERM2*STRAN(I1)*DTIME-EC*ESTRA/EREL+CONST6*DTEMP)
    END DO
C
C
C CREATE NEW JACOBIAN
C
C
      DO K1=1,6
      DO K2=1,6
       DDSDDE(K1,K2)=0
     END DO
      END DO
C
      DO K1=1,3
      DDSDDE(K1,K1)=EMOD*TERM1
      END DO
C
      DO K1=2,3
      N2=K1-1
     DO K2=1,N2
        DDSDDE(K2,K1)=EMOD*TERM3
        DDSDDE(K1,K2)=EMOD*TERM3
     END DO
      END DO
C
      DO K1=1,3
      I1=3+1
     DDSDDE(I1,I1)=EMOD
      END DO
    end SUBROUTINE UMAT



怎么都不知道问题出在哪，关系到毕业啊，都急死了。。。求大侠帮忙！！！！！！！！！！！！  

就没人帮忙么？？？

我是有心无力啊，帮你顶吧，之前也研究过一段时间，但是没有考虑粘弹性

楼主，您现在粘弹性怎么样了，能给指导一下吗，本构方程的各变量的意义能给解释一下吗，谢啦