Esope source files

   1 : C J2        SOURCE    CHAT      12/08/06    21:15:27     7479           
   2 :       SUBROUTINE J2(A,B)
   3 : C====================================================================== 
   4 : C Routine permettant de calculer le deuxiÃ¨me invariant du dÃ©viateur 
   5 : C de A 
   6 : C     OBJET                                                                     
   7 : C     -----                                                                     
   8 : C     CALCUL DE B = J2(A)                            
   9 : C                                                                               
  10 : C     ENTREES                                                                   
  11 : C     -------                                                                   
  12 : C     A(3,3) = MATRICE SYMETRIQUE          
  13 : C
  14 : C     SORTIES                                                                   
  15 : C     -------                                                                   
  16 : C     B      = J2(A) 
  17 : C                                                                               
  18 : C===============================================================                
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8 (A-H,O-Z)
  21 :       
  22 :       INTEGER NCAS
  23 :       REAL*8 A(3,3),TRAV(3,3),UNIT(3,3)
  24 : 
  25 :         DO I=1,3
  26 :          DO J=1,3
  27 :             IF (I.EQ.J) THEN
  28 :                UNIT(I,J) = 1.0D0
  29 :             ELSE
  30 :                UNIT(I,J) = 0.0D0
  31 :             ENDIF
  32 :          ENDDO
  33 :       ENDDO
  34 : 
  35 :       B = 0.0D0
  36 : 
  37 :       TRA_SE = A(1,1) + A(2,2) + A(3,3)
  38 : 
  39 :       DO I=1,3
  40 :          DO J=1,3
  41 :             TRAV(I,J) = 0.0D0
  42 :             TRAV(I,J) = A(I,J)-(1.0D0/3.0D0)*TRA_SE*UNIT(I,J)
  43 :          ENDDO
  44 :       ENDDO
  45 : 
  45 :       CALL DBLECONT(TRAV,TRAV,B)
  46 :       B = (3.0D0/2.0D0)*B
  47 :       B = DSQRT(B)
  48 : 
  49 :       RETURN
  50 :       END
  51 : 
  52 :  
  53 :  
  54 :  
  55 :  
  56 :  
  57 :