ABAQUS Subroutine Error

[ehsansharei]

Hello everybody.
I am totally new to implement ABAQUS user subroutine in my project and I have started to implement a simple linear elastic UMAT-subroutine within a one-element model to learn the procedure. I have written the subroutine which is attached to this message with the input file. I am simply unable to run the job or debug the code getting an error: "Problem during compilation - D:\Temp\Elasticity.for". I dont know if wrong Fortran compiler (32bit instead of 64) is running, or is there a code mistake in the subroutine. Could anyone check my subroutine within the givel input file please? Any help would work me a lot and is really appreciated.
Regards
Ehsan

SUBROUTINE UMAT(STRESS, STATEV, DDSDDE, SSE, SPD, SCD, RPL,
1 DDSDDT, DRPLDE, DRPLDT, STRAN, DSTRAN, TIME, DTIME, TEMP, DTEMP,
2 PREDEF, DPRED, CMNAME, NDI, NSHR, NTENS, NSTATV, PROPS, NPROPS,
3 COORDS, DROT, PNEWDT, CELENT, DFGRD0, DFGRD1, NOEL, NPT, LAYER,
4 KSPT, KSTEP, KINC)
C
INCLUDE ’ABA_PARAM.INC’
C
CHARACTER*8 CMNAME
C
DIMENSION STRESS(NTENS), STATEV(NSTATV), DDSDDE(NTENS, NTENS),
1 DDSDDT(NTENS), DRPLDE(NTENS), STRAN(NTENS), DSTRAN(NTENS),
2 PREDEF(1), DPRED(1), PROPS(NPROPS), COORDS(3), DROT(3, 3),
3 DFGRD0(3, 3), DFGRD1(3, 3)
C
DIMENSION EELAS(6), EPLAS(6), FLOW(6)
C
PARAMETER(ZERO=0.D0, ONE=1.D0, TWO=2.D0, THREE=3.D0, SIX=6.D0,
1 ENUMAX=.4999D0, NEWTON=10, TOLER=1.0D-6)
C ----------------------------------------------------------------
C UMAT FOR ISOTROPIC ELASTICITY
C CANNOT BE USED FOR PLANE STRESS
C ----------------------------------------------------------------
C PROPS(1) - E
C PROPS(2) - NU
C ----------------------------------------------------------------
C
IF (NDI.NE.3) THEN
WRITE (7, *) ’THIS UMAT MAY ONLY BE USED FOR ELEMENTS
1 WITH THREE DIRECT STRESS COMPONENTS’
CALL XIT
ENDIF
C
C ELASTIC PROPERTIES
EMOD=PROPS(1)
ENU=PROPS(2)
EBULK3=EMOD/(ONE-TWO*ENU)
EG2=EMOD/(ONE+ENU)
EG=EG2/TWO
EG3=THREE*EG
ELAM=(EBULK3-EG2)/THREE
C
C ELASTIC STIFFNESS
C
DO K1=1, NDI
DO K2=1, NDI
DDSDDE(K2, K1)=ELAM
END DO
DDSDDE(K1, K1)=EG2+ELAM
END DO
DO K1=NDI+1, NTENS
DDSDDE(K1 ,K1)=EG
END DO
C
C CALCULATE STRESS
C
DO K1=1, NTENS
DO K2=1, NTENS
STRESS(K2)=STRESS(K2)+DDSDDE(K2, K1)*DSTRAN(K1)
END DO
END DO
C
RETURN
END

 

[IceBreakerSours]

You must get comfortable with simpler subroutines before working with UMATs. I recommend playing with disp, dload, uhyper etc. Also, search the forum on how to set up your machine to run subroutines.

http://www.eng-tips.com/faqs.cfm?fid=376

http://www.eng-tips.com/faqs.cfm?fid=1083

 

[ehsansharei]

@ PostIceBreakerSours: Thanks for your advice. I will take a look on them to see if it is a better starting point. I want to thankfully ask for help from somebody who has the write establishment of his ABAQUS and visual fortran compiler to check my subroutine. I want to undrestand if my abaqus and visual fortran are correctly working, scince I am doubtfull if there is an error within the subroutine and model.
Regards
Ehsan