| ²é¿´: 1366 | »Ø¸´: 1 | |||
[½»Á÷]
ÎÒÓá°dimension"¶¨ÒåÁËÊý×飬¿ÉÊǵ÷ÊÔÌáʾÎÒûÓÐÉùÃ÷´ËÊý×é¡£
|
|
DIMENSION EELAS(6), EELASP(3),BBAR(6), BBARP(3), BBARN(3, 3), 1 DISTGR(3,3), GRF(3,3) PARAMETER(ZERO=0.D0, ONE=1.D0, TWO=2.D0, THREE=3.D0, FOUR=4.D0, 1 SIX=6.D0) C---------------------------------------------------------------------------------------------- C------------------------------------------------------------------------------------------------ EMOD=PROPS(1) ENU=PROPS(2) LAME1=EMOD*ENU/((1+ENU)*(1-2*ENU)) LAME2=0.5*EMOD/(1+ENU) C---------------------------------------------------------------------------------------------- C JACOBIAN OR DET DET=DFGRD1(1, 1)*DFGRD1(2, 2)*DFGRD1(3, 3) 1 -DFGRD1(1, 2)*DFGRD1(2, 1)*DFGRD1(3, 3) IF(NSHR.EQ.3) THEN DET=DET+DFGRD1(1, 2)*DFGRD1(2, 3)*DFGRD1(3, 1) 1 +DFGRD1(1, 3)*DFGRD1(3, 2)*DFGRD1(2, 1) 2 -DFGRD1(1, 3)*DFGRD1(3,1)*DFGRD1(2, 2) 3 -DFGRD1(2, 3)*DFGRD1(3, 2)*DFGRD1(1, 1) END IF C------------------------------------------------------------------------------------------------- C CALCULATE THE GROWTH FACTOR IF TIME(1)=1 THEN STATEV(NTENS+1)=1.0005 GR0=STATEV(NTENS+1) GO TO 110 GR0=STATEV(NTENS+1) AB=DET/GR0**3 IF (AB.LE.1.0) THEN STATEV(NTENS+1)=GR0 GO TO 110 END IF C---------------------------------------------------------------------- EPS=1.0E-08 CALL DNEWT(GR0,EPS,RES,L) IF (L.NE.0) THEN WRITE(*,10) GR0 STATEV(NTENS+1) =GR0 END IF 10 FORMAT(001X,'GR0=',E15.6) END C----------------------------------------------------------------------- C CALCULATE GROTH TENSOR 110 DO K1=1,3 DO K2=1,3 GRF(K1,K2)=0 END DO END DO DO K3=1,3 GRF(K3,K3)=GR0 END DO C--------------------------------------------------------------------------- C CALCULATE ELASTIC TENSOR SCALE=ONE/GR0 DO K1=1, 3 DO K2=1, 3 DISTGR(K2,K1)=SCALE*DFGRD1(K2,K1) END DO END DO C---------------------------------------------------------------------------- C------------------------------------------------------------------------------ BBAR(1)=DISTGR(1, 1)**2+DISTGR(1, 2)**2+DISTGR(1, 3)**2 BBAR(2)=DISTGR(2, 1)**2+DISTGR(2, 2)**2+DISTGR(2, 3)**2 BBAR(3)=DISTGR(3, 3)**2+DISTGR(3, 1)**2+DISTGR(3, 2)**2 BBAR(4)=DISTGR(1, 1)*DISTGR(2, 1)+DISTGR(1, 2)*DISTGR(2, 2) 1 +DISTGR(1, 3)*DISTGR(2, 3) IF(NSHR.EQ.3) THEN BBAR(5)=DISTGR(1, 1)*DISTGR(3, 1)+DISTGR(1, 2)*DISTGR(3, 2) 1 +DISTGR(1, 3)*DISTGR(3, 3) BBAR(6)=DISTGR(2, 1)*DISTGR(3, 1)+DISTGR(2, 2)*DISTGR(3, 2) 1 +DISTGR(2, 3)*DISTGR(3, 3) ¿ÉÊǵ÷ÊÔµÄʱºòÀÏÌáʾ 1>G:\Fortran\Console22\Console22\Source1.for(70): error #6410: This name has not been declared as an array or a function. [GRF] 1>G:\Fortran\Console22\Console22\Source1.for(74): error #6410: This name has not been declared as an array or a function. [GRF] 1>G:\Fortran\Console22\Console22\Source1.for(81): error #6410: This name has not been declared as an array or a function. [DISTGR] 1>G:\Fortran\Console22\Console22\Source1.for(86): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(87): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(88): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(89): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(92): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(94): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(102): error #6410: This name has not been declared as an array or a function. [STRESS] 1>G:\Fortran\Console22\Console22\Source1.for(105): error #6410: This name has not been declared as an array or a function. [STRESS] |
» ²ÂÄãϲ»¶
²ÄÁÏ¿¼Ñе÷¼Á
ÒѾÓÐ3È˻ظ´
²ÄÁϵ÷¼Á
ÒѾÓÐ12È˻ظ´
Ó¢Ò»ÊýÒ»408£¬×Ü·Ö284£¬¶þÕ½Õæ³ÏÇóµ÷¼Á
ÒѾÓÐ14È˻ظ´
085410 Ò»Ö¾Ô¸211 22408·ÖÊý359Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
271Çóµ÷¼Á
ÒѾÓÐ19È˻ظ´
385·Ö ÉúÎïѧ£¨071000£©Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
Ò»Ö¾Ô¸°²»Õ´óѧ¼ÆËã»ú¿ÆÑ§Óë¼¼Êõѧ˶£¬331·ÖÇóµ÷¼Á
ÒѾÓÐ3È˻ظ´
318Çóµ÷¼Á£¬¼ÆËã²ÄÁÏ·½Ïò
ÒѾÓÐ8È˻ظ´
291Çóµ÷¼Á
ÒѾÓÐ25È˻ظ´
Ò»Ö¾Ô¸±±¾©¿Æ¼¼´óѧ085601²ÄÁϹ¤³ÌÓ¢Ò»Êý¶þ³õÊÔ×Ü·Ö335Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´


2Â¥2015-05-03 10:03:09














»Ø¸´´ËÂ¥