subroutine vumat(
C Read only(unmodifiable)variabies-
nblock,ndir,nshr,nstatev,nfieldv,nprops,lanneal,
step Time,totalTime,dt,cmname,coordMp,charLength,
props,density,strainInc,relSpinInc,
tempOld,stretchOld,defgradOld,fieldOld,
stressOld,stateOld,enerInternOld,enerInelasOld,
tempNew,stretchNew,defgradNew,fieldNew,
C Write only(modifiable)variables-
stressNew,stateNew,enerInternNew,enerInelasNew)
C
include ‘vaba_param.inc’
C
dimension props(nprops),density(nblock),coordMp(nblock,*),
charLength(nblock),strainInc(nblock,ndir+nshr),tempOld(nblock),
relSpinInc(nblock,nshr),stretchOld(nblock,ndir+nshr),
defgradOld(nblock,ndir+nshr+nshr),fieldOld(nblock,nfieldv),
stressOld(nblock,ndir+nshr),stateOld(nblock,nstatev),
enerInternOld(nblock),enerInelasOld(nblock),tempNew(nblock),
stretchNew(nblock,ndir+nshr),defgradNew(nblock,ndir+nshr+nshr),
fieldNew(nblock,nfieldv),
stressNew(nblock,ndir+nshr),stateNew(nblock,nstatev),
enerInternNew(nblock),enerInelasNew(nblock)
C
character*80 cmname
C
DIMENSION C(6,6),STH(6,2),strainNew(6),SN(6)
C GET THE MATERIAL PROPERTIES
E1 =PROPS(1) !E1,YOUNG\'S MODULUS IN DIRECTION 1
E2 =PROPS(2) !E2=E3,YOUNG\'S MODULUS IN DIRECTION2&3
XNU12=PROPS(3) !POISON\'S RATIO POI_12,XNU13=XNU12
XNU23=PROPS(4) !POISON\'S RATIO POI_23,XNU32=XNU23
G12 =PROPS(5) !G12=G13,SHEAR MODULUS IN 12 & 13 PLANE
G23 =PROPS(6) !G23,SHEAR MODULUS IN 23 PLANE
XNU21=XNU12*E2/E1 !POISON\'S RATIO POI_21,XNU31=XNU21
C GET THE FAILURE PROPERTIES
STH=O
STH(1,1)=PROPS(7) !FAILURE STRESS IN 1 DIRECTION IN TENSION
STH(1,2)=PROPS(8) !FAILURE STRESS IN 1 DIRECTION IN COMPRESSION
STH(2,1)=PROPS(9) !FAILURE STRESS IN 2 DIRECTION IN TENSION
STH(2,2)=PROPS(10)!FAILURE STRESS IN 2 DIRECTION IN COMPRESSION
STH(3,1)=PROPS(11)!FAILURE STRESS IN 3 DIRECTION IN TENSION
STH(3,2)=PROPS(12)!FAILURE STRESS IN 3 DIRECTION IN COMPRESSION
STH(4,1)=PROPS(13)!FAILURE STRESS IN SHEAR IN 1-2 PLANE
STH(5,1)=PROPS(14)!FAILURE STRESS IN SHEAR IN 2-3 PLANE
STH(6,1)=STH(4,1) !FAILURE STRESS IN SHEAR IN 1-3 PLANE
C ***************
C Tsai-wu Coefficients
F1=1/STH(1,1)-1/STH(1,2)
F2=1/STH(2,1)-1/STH(2,2)
F3=1/STH(3,1)-1/STH(3,2)
F11=1/(STH(1,1)*STH(1,2))
F22=1/(STH(2,1)*STH(2,2))
F33=1/(STH(3,1)*STH(3,2))
F44=1/STH(4,1)**2
F55=1/STH(5,1)**2
F66=1/STH(6,1)**2
F12=-(F11*F22)**0.5
F23=-(F22*F33)**0.5
F13=-(F11*F33)**0.5
C ***************
C STIFFNESS MATRIX C(6,6)
C=0
RNU=1/(1-2*XNU12*XNU21-XNU23**2-2*XNU12*XNU21*XNU23)
C(1,1)=E1*(1-XNU23**2)*RNU
C(2,2)=E2*(1-XNU12*XNU21)*RNU
C(1,2)=E1*(XNU21+XNU21*XNU23)*RNU
C(2,3)=E2*(XNU23+XNU12*XNU21)*RNU
C(4,4)=G12
C(5,5)=G23
C(6,6)=G12
C(2,1)=C(1,2)
C(1,3)=C(1,2)
C(3,1)=C(1,2)
C(3,3)=C(2,2)
C(3,2)=C(2,3)
C ***************
do 100k=1,nblock
do i=1,6
if(i<4)then
strainNew(i)=stateOld(k,i)+strainInc(k,i)
else
strainNew(i)=stateOld(k,i)+2*strainInc(k,i)
endif
stateNew(k,i)=strainNew(i)
enddo
do i=1,6
stressNew(k,i)=0
do j=1,6
stressNew(k,i)=stressNew(k,i)+C(i,j)*strainNew(j)
enddo
SN(i)=stressNew(k,i)
enddo
FTW=F1*SN(1)+F2*SN(2)+F3*SN(3)+F11*SN(1)**2+F22*SN(2)**2+
F33*SN(3)**2+F44*SN(4)**2+F55*SN(5)**2+F66*SN(6)**2+
F12*SN(1)*SN(2)+F23*SN(2)*SN(3)+F13*SN(1)*SN(3)
stateNew(k,7)=FTW
if(FTW<1)then
stateNew(k,8)=1
else
stateNew(k,8)=0
endif
continue
return
end
c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(11): error #5102: Cannot open include file \'vaba_param.inc\'
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(13): error #5082: Syntax error, found \'&\' when expecting one of: <IDENTIFIER>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(14): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(14): error #5082: Syntax error, found \'&\' when expecting one of: ( <IDENTIFIER> <CHAR_CON_KIND_PARAM> <CHAR_NAM_KIND_PARAM> <CHARACTER_CONSTANT> <INTEGER_CONSTANT> ...
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(15): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(15): error #5082: Syntax error, found \'&\' when expecting one of: ( <IDENTIFIER> <CHAR_CON_KIND_PARAM> <CHAR_NAM_KIND_PARAM> <CHARACTER_CONSTANT> <INTEGER_CONSTANT> ...
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(16): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(16): error #5082: Syntax error, found \'&\' when expecting one of: ( <IDENTIFIER> <CHAR_CON_KIND_PARAM> <CHAR_NAM_KIND_PARAM> <CHARACTER_CONSTANT> <INTEGER_CONSTANT> ...
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(17): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(17): error #5082: Syntax error, found \'&\' when expecting one of: ( <IDENTIFIER> <CHAR_CON_KIND_PARAM> <CHAR_NAM_KIND_PARAM> <CHARACTER_CONSTANT> <INTEGER_CONSTANT> ...
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(18): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(18): error #5082: Syntax error, found \'&\' when expecting one of: ( <IDENTIFIER> <CHAR_CON_KIND_PARAM> <CHAR_NAM_KIND_PARAM> <CHARACTER_CONSTANT> <INTEGER_CONSTANT> ...
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(19): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(19): error #5082: Syntax error, found \'&\' when expecting one of: ( <IDENTIFIER> <CHAR_CON_KIND_PARAM> <CHAR_NAM_KIND_PARAM> <CHARACTER_CONSTANT> <INTEGER_CONSTANT> ...
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(20): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(21): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(21): error #5082: Syntax error, found \'&\' when expecting one of: ( <IDENTIFIER> <CHAR_CON_KIND_PARAM> <CHAR_NAM_KIND_PARAM> <CHARACTER_CONSTANT> <INTEGER_CONSTANT> ...
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(22): error #5082: Syntax error, found \',\' when expecting one of: ( % [ . = =>
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(85): error #6410: This name has not been declared as an array or a function. [STATENEW]
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(88): error #6410: This name has not been declared as an array or a function. [STRESSNEW]
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(90): error #6410: This name has not been declared as an array or a function. [STRESSNEW]
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(97): error #6410: This name has not been declared as an array or a function. [STATENEW]
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(99): error #6410: This name has not been declared as an array or a function. [STATENEW]
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(101): error #6410: This name has not been declared as an array or a function. [STATENEW]
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(78): error #6321: An unterminated block exists.
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(78): error #6323: This label is not defined in this scoping unit. [100]
1>c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for(22): remark #7713: This statement function has not been used. [ENERINTERNNEW]
1>compilation aborted for c:\\users\\sn\\documents\\visual studio 2010\\Projects\\Console1\\Console1\\Source1.for (code 1) |