Télécharger la plaquette Cast3M


fr en

Télécharger uo2rot.eso

Retour à la liste

Numérotation des lignes : 

uo2rot

  1. C UO2ROT    SOURCE    CHAT      05/01/13    03:59:26     5004
  2.       SUBROUTINE UO2ROT(SIGMA0,SIGMA,VCA1,VCA2,VCA3,ICAS)
  3. C
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. C
  7.       DIMENSION SIGMA0(6),SIGMA(6)
  8.       DIMENSION VCA1(3),VCA2(3),VCA3(3)
  9. C
  10. C
  11. C========================================================================
  12. C       CALCUL DE SIGMA0 DANS LE REPERE DE CALCUL VCA1 VCA2 VCA3
  13. C                    RESULTAT NOTE : SIGMA
  14. C========================================================================
  15. C
  16.       IF(ICAS.EQ.1) THEN
  17. C
  18.       SIGMA(1)=SIGMA0(1)*VCA1(1)*VCA1(1)+SIGMA0(4)*VCA1(2)*VCA1(1)
  19.      &+SIGMA0(5)*VCA1(3)*VCA1(1)+SIGMA0(4)*VCA1(1)*VCA1(2)
  20.      &+SIGMA0(2)*VCA1(2)*VCA1(2)+SIGMA0(6)*VCA1(3)*VCA1(2)
  21.      &+SIGMA0(5)*VCA1(1)*VCA1(3)+SIGMA0(6)*VCA1(2)*VCA1(3)
  22.      &+SIGMA0(3)*VCA1(3)*VCA1(3)
  23. C
  24.       SIGMA(2)=SIGMA0(1)*VCA2(1)*VCA2(1)+SIGMA0(4)*VCA2(2)*VCA2(1)
  25.      &+SIGMA0(5)*VCA2(3)*VCA2(1)+SIGMA0(4)*VCA2(1)*VCA2(2)
  26.      &+SIGMA0(2)*VCA2(2)*VCA2(2)+SIGMA0(6)*VCA2(3)*VCA2(2)
  27.      &+SIGMA0(5)*VCA2(1)*VCA2(3)+SIGMA0(6)*VCA2(2)*VCA2(3)
  28.      &+SIGMA0(3)*VCA2(3)*VCA2(3)
  29. C
  30.       SIGMA(3)=SIGMA0(1)*VCA3(1)*VCA3(1)+SIGMA0(4)*VCA3(2)*VCA3(1)
  31.      &+SIGMA0(5)*VCA3(3)*VCA3(1)+SIGMA0(4)*VCA3(1)*VCA3(2)
  32.      &+SIGMA0(2)*VCA3(2)*VCA3(2)+SIGMA0(6)*VCA3(3)*VCA3(2)
  33.      &+SIGMA0(5)*VCA3(1)*VCA3(3)+SIGMA0(6)*VCA3(2)*VCA3(3)
  34.      &+SIGMA0(3)*VCA3(3)*VCA3(3)
  35. C
  36.       SIGMA(4)=SIGMA0(1)*VCA1(1)*VCA2(1)+SIGMA0(4)*VCA1(2)*VCA2(1)
  37.      &+SIGMA0(5)*VCA1(3)*VCA2(1)+SIGMA0(4)*VCA1(1)*VCA2(2)
  38.      &+SIGMA0(2)*VCA1(2)*VCA2(2)+SIGMA0(6)*VCA1(3)*VCA2(2)
  39.      &+SIGMA0(5)*VCA1(1)*VCA2(3)+SIGMA0(6)*VCA1(2)*VCA2(3)
  40.      &+SIGMA0(3)*VCA1(3)*VCA2(3)
  41. C
  42.       SIGMA(5)=SIGMA0(1)*VCA1(1)*VCA3(1)+SIGMA0(4)*VCA1(2)*VCA3(1)
  43.      &+SIGMA0(5)*VCA1(3)*VCA3(1)+SIGMA0(4)*VCA1(1)*VCA3(2)
  44.      &+SIGMA0(2)*VCA1(2)*VCA3(2)+SIGMA0(6)*VCA1(3)*VCA3(2)
  45.      &+SIGMA0(5)*VCA1(1)*VCA3(3)+SIGMA0(6)*VCA1(2)*VCA3(3)
  46.      &+SIGMA0(3)*VCA1(3)*VCA3(3)
  47. C
  48.       SIGMA(6)=SIGMA0(1)*VCA2(1)*VCA3(1)+SIGMA0(4)*VCA2(2)*VCA3(1)
  49.      &+SIGMA0(5)*VCA2(3)*VCA3(1)+SIGMA0(4)*VCA2(1)*VCA3(2)
  50.      &+SIGMA0(2)*VCA2(2)*VCA3(2)+SIGMA0(6)*VCA2(3)*VCA3(2)
  51.      &+SIGMA0(5)*VCA2(1)*VCA3(3)+SIGMA0(6)*VCA2(2)*VCA3(3)
  52.      &+SIGMA0(3)*VCA2(3)*VCA3(3)
  53. C
  54.       ENDIF
  55. C
  56. C
  57. C=======================================================================
  58. C            CALCUL DE SIGMA0 DANS LE REPERE GLOBAL
  59. C                    RESULTAT NOTE : SIGMA
  60. C=======================================================================
  61. C
  62.       IF(ICAS.EQ.2) THEN
  63. C
  64.       SIGMA(1)=SIGMA0(1)*VCA1(1)*VCA1(1)+SIGMA0(4)*VCA2(1)*VCA1(1)
  65.      &+SIGMA0(5)*VCA3(1)*VCA1(1)+SIGMA0(4)*VCA1(1)*VCA2(1)
  66.      &+SIGMA0(2)*VCA2(1)*VCA2(1)+SIGMA0(6)*VCA3(1)*VCA2(1)
  67.      &+SIGMA0(5)*VCA1(1)*VCA3(1)+SIGMA0(6)*VCA2(1)*VCA3(1)
  68.      &+SIGMA0(3)*VCA3(1)*VCA3(1)
  69. C
  70.       SIGMA(2)=SIGMA0(1)*VCA1(2)*VCA1(2)+SIGMA0(4)*VCA2(2)*VCA1(2)
  71.      &+SIGMA0(5)*VCA3(2)*VCA1(2)+SIGMA0(4)*VCA1(2)*VCA2(2)
  72.      &+SIGMA0(2)*VCA2(2)*VCA2(2)+SIGMA0(6)*VCA3(2)*VCA2(2)
  73.      &+SIGMA0(5)*VCA1(2)*VCA3(2)+SIGMA0(6)*VCA2(2)*VCA3(2)
  74.      &+SIGMA0(3)*VCA3(2)*VCA3(2)
  75. C
  76.       SIGMA(3)=SIGMA0(1)*VCA1(3)*VCA1(3)+SIGMA0(4)*VCA2(3)*VCA1(3)
  77.      &+SIGMA0(5)*VCA3(3)*VCA1(3)+SIGMA0(4)*VCA1(3)*VCA2(3)
  78.      &+SIGMA0(2)*VCA2(3)*VCA2(3)+SIGMA0(6)*VCA3(3)*VCA2(3)
  79.      &+SIGMA0(5)*VCA1(3)*VCA3(3)+SIGMA0(6)*VCA2(3)*VCA3(3)
  80.      &+SIGMA0(3)*VCA3(3)*VCA3(3)
  81. C
  82.       SIGMA(4)=SIGMA0(1)*VCA1(1)*VCA1(2)+SIGMA0(4)*VCA2(1)*VCA1(2)
  83.      &+SIGMA0(5)*VCA3(1)*VCA1(2)+SIGMA0(4)*VCA1(1)*VCA2(2)
  84.      &+SIGMA0(2)*VCA2(1)*VCA2(2)+SIGMA0(6)*VCA3(1)*VCA2(2)
  85.      &+SIGMA0(5)*VCA1(1)*VCA3(2)+SIGMA0(6)*VCA2(1)*VCA3(2)
  86.      &+SIGMA0(3)*VCA3(1)*VCA3(2)
  87. C
  88.       SIGMA(5)=SIGMA0(1)*VCA1(1)*VCA1(3)+SIGMA0(4)*VCA2(1)*VCA1(3)
  89.      &+SIGMA0(5)*VCA3(1)*VCA1(3)+SIGMA0(4)*VCA1(1)*VCA2(3)
  90.      &+SIGMA0(2)*VCA2(1)*VCA2(3)+SIGMA0(6)*VCA3(1)*VCA2(3)
  91.      &+SIGMA0(5)*VCA1(1)*VCA3(3)+SIGMA0(6)*VCA2(1)*VCA3(3)
  92.      &+SIGMA0(3)*VCA3(1)*VCA3(3)
  93. C
  94.       SIGMA(6)=SIGMA0(1)*VCA1(2)*VCA1(3)+SIGMA0(4)*VCA2(2)*VCA1(3)
  95.      &+SIGMA0(5)*VCA3(2)*VCA1(3)+SIGMA0(4)*VCA1(2)*VCA2(3)
  96.      &+SIGMA0(2)*VCA2(2)*VCA2(3)+SIGMA0(6)*VCA3(2)*VCA2(3)
  97.      &+SIGMA0(5)*VCA1(2)*VCA3(3)+SIGMA0(6)*VCA2(2)*VCA3(3)
  98.      &+SIGMA0(3)*VCA3(2)*VCA3(3)
  99. C
  100.       ENDIF
  101. C
  102.       RETURN
  103.       END
  104.  
  105.  
  106.  

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