| 查看: 2508 | 回复: 1 | ||
[求助]
子程序调试错误求助
|
|
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) |
回复此楼» 猜你喜欢
实验室接单子
已经有6人回复
-
假如你的研究生提出不合理要求
已经有11人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
shuishengy
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 798.4
- 散金: 215
- 红花: 2
- 帖子: 208
- 在线: 159.8小时
- 虫号: 3293850
- 注册: 2014-06-26
- 专业: 交通工程
2楼2019-11-25 19:41:50