Télécharger j2.eso

Retour à la liste

Numérotation des lignes :

  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