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