Hello all, I am in a desperate need of some expert help
For a newbie at Abaqus with subroutines,
After months of trying, finally installed Abaqus, verified results all passed.
(Abaqus 6.13, Microsoft Visual Studio 2010, Intel Parallel Studio XE 2013)
Error: Problem during compilation: What is wrong with the code below subroutine uhyper? (A former student had this subroutine working, may be in earlier versions but had it working)
SUBROUTINE UHYPER(BI1,BI2,AJ,U,UI1,UI2,UI3,TEMP,NOEL,CMNAME,
\$ INCMPFLAG,NUMSTATEV,STATEV,NUMFIELDV,
\$ FIELDV,FIELDVINC,NUMPROPS,PROPS)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME
DIMENSION U(2),UI1(3),UI2(6),UI3(6),STATEV(*),FIELDV(*),
\$ FIELDVINC(*),PROPS(*)
C
PARAMETER (ZERO=0.0D0,ONE=1.0D0, TWO=2.0D0, THREE=3.0D0,SIX=6.0D0)
C
IF (CMNAME(1:6) .EQ. 'ABC') THEN
C10 = 0.08423
C01 = 0.00505
C11 = 0.76506
C20 = 1.5
C30 = 0.04238
D1 = 0.000008928
U=C10*(BI1-THREE)+C01*(BI2-THREE)+(C20*(BI1-THREE)**2)+C11*
\$ (BI1-THREE)*(BI2-THREE)+(C30*(BI1-THREE)**3)+(((AJ-ONE)**2)/D1)
UI1(1)=C10+C11*(BI2-THREE)+C20*(TWO*BI1-SIX)+
\$ THREE*C30*(BI1-THREE)**2
UI1(2)=C01+C11*(BI1-THREE)
UI1(3)=(TWO*AJ-TWO)/D1
UI2(1)=TWO*C20+THREE*C30*(TWO*BI1-SIX)
UI2(2)= ZERO
UI2(3)=TWO/D1
UI2(4)=C11
UI2(5)=ZERO
UI2(6)=ZERO
UI3(1)=ZERO
UI3(2)=ZERO
UI3(3)=ZERO
UI3(4)=ZERO
UI3(5)=ZERO
UI3(6)=ZERO
ELSE IF (CMNAME(1:5) .EQ. 'DEF') THEN
C10 = 0.00355
C01 = 0.00066
C11 = 0.01868
C20 = 0.02154
C30 = 0.01977
D1=0.000000421
U=C10*(BI1-THREE)+C01*(BI2-THREE)+(C20*(BI1-THREE)**2)+C11*
\$ (BI1-THREE)*(BI2-THREE)+(C30*(BI1-THREE)**3)+(((AJ-ONE)**2)/D1)
UI1(1)=C10+C11*(BI2-THREE)+C20*(TWO*BI1-SIX)+
\$ THREE*C30*(BI1-THREE)**2
UI1(2)=C01+C11*(BI1-THREE)
UI1(3)=(TWO*AJ-TWO)/D1
UI2(1)=TWO*C20+THREE*C30*(TWO*BI1-SIX)
UI2(2)= ZERO
UI2(3)=TWO/D1
UI2(4)=C11
UI2(5)=ZERO
UI2(6)=ZERO
UI3(1)=ZERO
UI3(2)=ZERO
UI3(3)=ZERO
UI3(4)=ZERO
UI3(5)=ZERO
UI3(6)=ZERO
ELSE IF (CMNAME(1:10) .EQ. 'GHI') THEN
C10 = 0.00714
C01 = 0.00063
C11 = 0.09579
C20 = 0.00803
C30 = 0.09931
D1=0.00000077
U=C10*(BI1-THREE)+C01*(BI2-THREE)+(C20*(BI1-THREE)**2)+C11*
\$ (BI1-THREE)*(BI2-THREE)+(C30*(BI1-THREE)**3)+(((AJ-ONE)**2)/D1)
UI1(1)=C10+C11*(BI2-THREE)+C20*(TWO*BI1-SIX)+
\$ THREE*C30*(BI1-THREE)**2
UI1(2)=C01+C11*(BI1-THREE)
UI1(3)=(TWO*AJ-TWO)/D1
UI2(1)=TWO*C20+THREE*C30*(TWO*BI1-SIX)
UI2(2)= ZERO
UI2(3)=TWO/D1
UI2(4)=C11
UI2(5)=ZERO
UI2(6)=ZERO
UI3(1)=ZERO
UI3(2)=ZERO
UI3(3)=ZERO
UI3(4)=ZERO
UI3(5)=ZERO
UI3(6)=ZERO
ELSE IF (CMNAME(1:7) .EQ. 'JKL') THEN
C10=0.00753
C01=0.00999
C20=0.02063
C11=0.00039
C30=0.00078
D1=0.000001752
U =C10*(BI1-THREE)+C01*(BI2-THREE)+C20*(BI1-THREE)**2+C11*
\$ (BI1-THREE)*(BI2-THREE)+C30*(BI2-THREE)**2+(((AJ-ONE)**2)/D1)
UI1(1)=C10+C11*(BI2-THREE)+C20*(2*BI1-SIX)
UI1(2)=C01+C11*(BI1-THREE)+C30*(TWO*BI2-SIX)
UI1(3)=(TWO*AJ-TWO)/D1
UI2(1)=TWO*C20
UI2(2)= TWO*C30
UI2(3)=TWO/D1
UI2(4)=C11
UI2(5)=ZERO
UI2(6)=ZERO
UI3(1)=ZERO
UI3(2)=ZERO
UI3(3)=ZERO
UI3(4)=ZERO
UI3(5)=ZERO
UI3(6)=ZERO
ELSE IF (CMNAME(1:5) .EQ. 'MNO') THEN
C10=0.00116
C01=0.01168
C20=0.00858
C11=0.0006
C30=0.00022
D1=0.000001284
U =C10*(BI1-THREE)+C01*(BI2-THREE)+C20*(BI1-THREE)**2+C11*
\$ (BI1-THREE)*(BI2-THREE)+C30*(BI2-THREE)**2+(((AJ-ONE)**2)/D1)
UI1(1)=C10+C11*(BI2-THREE)+C20*(2*BI1-SIX)
UI1(2)=C01+C11*(BI1-THREE)+C30*(TWO*BI2-SIX)
UI1(3)=(TWO*AJ-TWO)/D1
UI2(1)=TWO*C20
UI2(2)= TWO*C30
UI2(3)=TWO/D1
UI2(4)=C11
UI2(5)=ZERO
UI2(6)=ZERO
UI3(1)=ZERO
UI3(2)=ZERO
UI3(3)=ZERO
UI3(4)=ZERO
UI3(5)=ZERO
UI3(6)=ZERO
ELSE
C10=0.001144
C01=0.01357
C20=0.06238
C11=0.02626
C30=0.00761
D1=0.0000014714
U =C10*(BI1-THREE)+C01*(BI2-THREE)+C20*(BI1-THREE)**2+C11*
\$ (BI1-THREE)*(BI2-THREE)+C30*(BI2-THREE)**2+(((AJ-ONE)**2)/D1)
UI1(1)=C10+C11*(BI2-THREE)+C20*(2*BI1-SIX)
UI1(2)=C01+C11*(BI1-THREE)+C30*(TWO*BI2-SIX)
UI1(3)=(TWO*AJ-TWO)/D1
UI2(1)=TWO*C20
UI2(2)= TWO*C30
UI2(3)=TWO/D1
UI2(4)=C11
UI2(5)=ZERO
UI2(6)=ZERO
UI3(1)=ZERO
UI3(2)=ZERO
UI3(3)=ZERO
UI3(4)=ZERO
UI3(5)=ZERO
UI3(6)=ZERO
C123456
END IF
RETURN
END
