Télécharger la plaquette Cast3M


fr en

Télécharger j2.eso

Retour à la liste

Numérotation des lignes : 

j2

  1. C J2        SOURCE    CB215821  16/04/21    21:17:19     8920
  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.  
  46.       CALL DBLECONT(TRAV,TRAV,B)
  47.       B = (3.0D0/2.0D0)*B
  48.       B = SQRT(B)
  49.  
  50.       RETURN
  51.       END
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales