Télécharger ottrot.eso

Retour à la liste

Numérotation des lignes :

  1. C OTTROT SOURCE CHAT 05/01/13 02:08:08 5004
  2. SUBROUTINE OTTROT(SIGMA0,DDEF0,SIGMA,DDEF,
  3. & VCA1,VCA2,VCA3,ICAS)
  4. C
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. C
  8. DIMENSION SIGMA0(6),DDEF0(6)
  9. DIMENSION SIGMA(6), DDEF(6)
  10. DIMENSION VCA1(3),VCA2(3),VCA3(3)
  11. C
  12. C
  13. C========================================================================
  14. C CALCUL DE DDEF0 ET SIGMA0 DANS LE REPERE DE CALCUL
  15. C RESULTAT NOTE : DDEF ET SIGMA
  16. C========================================================================
  17. C
  18. IF(ICAS.EQ.1) THEN
  19. C
  20. SIGMA(1)=SIGMA0(1)*VCA1(1)*VCA1(1)+SIGMA0(4)*VCA1(2)*VCA1(1)
  21. &+SIGMA0(5)*VCA1(3)*VCA1(1)+SIGMA0(4)*VCA1(1)*VCA1(2)
  22. &+SIGMA0(2)*VCA1(2)*VCA1(2)+SIGMA0(6)*VCA1(3)*VCA1(2)
  23. &+SIGMA0(5)*VCA1(1)*VCA1(3)+SIGMA0(6)*VCA1(2)*VCA1(3)
  24. &+SIGMA0(3)*VCA1(3)*VCA1(3)
  25. C
  26. SIGMA(2)=SIGMA0(1)*VCA2(1)*VCA2(1)+SIGMA0(4)*VCA2(2)*VCA2(1)
  27. &+SIGMA0(5)*VCA2(3)*VCA2(1)+SIGMA0(4)*VCA2(1)*VCA2(2)
  28. &+SIGMA0(2)*VCA2(2)*VCA2(2)+SIGMA0(6)*VCA2(3)*VCA2(2)
  29. &+SIGMA0(5)*VCA2(1)*VCA2(3)+SIGMA0(6)*VCA2(2)*VCA2(3)
  30. &+SIGMA0(3)*VCA2(3)*VCA2(3)
  31. C
  32. SIGMA(3)=SIGMA0(1)*VCA3(1)*VCA3(1)+SIGMA0(4)*VCA3(2)*VCA3(1)
  33. &+SIGMA0(5)*VCA3(3)*VCA3(1)+SIGMA0(4)*VCA3(1)*VCA3(2)
  34. &+SIGMA0(2)*VCA3(2)*VCA3(2)+SIGMA0(6)*VCA3(3)*VCA3(2)
  35. &+SIGMA0(5)*VCA3(1)*VCA3(3)+SIGMA0(6)*VCA3(2)*VCA3(3)
  36. &+SIGMA0(3)*VCA3(3)*VCA3(3)
  37. C
  38. SIGMA(4)=SIGMA0(1)*VCA1(1)*VCA2(1)+SIGMA0(4)*VCA1(2)*VCA2(1)
  39. &+SIGMA0(5)*VCA1(3)*VCA2(1)+SIGMA0(4)*VCA1(1)*VCA2(2)
  40. &+SIGMA0(2)*VCA1(2)*VCA2(2)+SIGMA0(6)*VCA1(3)*VCA2(2)
  41. &+SIGMA0(5)*VCA1(1)*VCA2(3)+SIGMA0(6)*VCA1(2)*VCA2(3)
  42. &+SIGMA0(3)*VCA1(3)*VCA2(3)
  43. C
  44. SIGMA(5)=SIGMA0(1)*VCA1(1)*VCA3(1)+SIGMA0(4)*VCA1(2)*VCA3(1)
  45. &+SIGMA0(5)*VCA1(3)*VCA3(1)+SIGMA0(4)*VCA1(1)*VCA3(2)
  46. &+SIGMA0(2)*VCA1(2)*VCA3(2)+SIGMA0(6)*VCA1(3)*VCA3(2)
  47. &+SIGMA0(5)*VCA1(1)*VCA3(3)+SIGMA0(6)*VCA1(2)*VCA3(3)
  48. &+SIGMA0(3)*VCA1(3)*VCA3(3)
  49. C
  50. SIGMA(6)=SIGMA0(1)*VCA2(1)*VCA3(1)+SIGMA0(4)*VCA2(2)*VCA3(1)
  51. &+SIGMA0(5)*VCA2(3)*VCA3(1)+SIGMA0(4)*VCA2(1)*VCA3(2)
  52. &+SIGMA0(2)*VCA2(2)*VCA3(2)+SIGMA0(6)*VCA2(3)*VCA3(2)
  53. &+SIGMA0(5)*VCA2(1)*VCA3(3)+SIGMA0(6)*VCA2(2)*VCA3(3)
  54. &+SIGMA0(3)*VCA2(3)*VCA3(3)
  55. C
  56. C
  57. DDEF(1)=DDEF0(1)*VCA1(1)*VCA1(1)+DDEF0(4)*VCA1(2)*VCA1(1)
  58. &+DDEF0(5)*VCA1(3)*VCA1(1)+DDEF0(4)*VCA1(1)*VCA1(2)
  59. &+DDEF0(2)*VCA1(2)*VCA1(2)+DDEF0(6)*VCA1(3)*VCA1(2)
  60. &+DDEF0(5)*VCA1(1)*VCA1(3)+DDEF0(6)*VCA1(2)*VCA1(3)
  61. &+DDEF0(3)*VCA1(3)*VCA1(3)
  62. C
  63. DDEF(2)=DDEF0(1)*VCA2(1)*VCA2(1)+DDEF0(4)*VCA2(2)*VCA2(1)
  64. &+DDEF0(5)*VCA2(3)*VCA2(1)+DDEF0(4)*VCA2(1)*VCA2(2)
  65. &+DDEF0(2)*VCA2(2)*VCA2(2)+DDEF0(6)*VCA2(3)*VCA2(2)
  66. &+DDEF0(5)*VCA2(1)*VCA2(3)+DDEF0(6)*VCA2(2)*VCA2(3)
  67. &+DDEF0(3)*VCA2(3)*VCA2(3)
  68. C
  69. DDEF(3)=DDEF0(1)*VCA3(1)*VCA3(1)+DDEF0(4)*VCA3(2)*VCA3(1)
  70. &+DDEF0(5)*VCA3(3)*VCA3(1)+DDEF0(4)*VCA3(1)*VCA3(2)
  71. &+DDEF0(2)*VCA3(2)*VCA3(2)+DDEF0(6)*VCA3(3)*VCA3(2)
  72. &+DDEF0(5)*VCA3(1)*VCA3(3)+DDEF0(6)*VCA3(2)*VCA3(3)
  73. &+DDEF0(3)*VCA3(3)*VCA3(3)
  74. C
  75. DDEF(4)=DDEF0(1)*VCA1(1)*VCA2(1)+DDEF0(4)*VCA1(2)*VCA2(1)
  76. &+DDEF0(5)*VCA1(3)*VCA2(1)+DDEF0(4)*VCA1(1)*VCA2(2)
  77. &+DDEF0(2)*VCA1(2)*VCA2(2)+DDEF0(6)*VCA1(3)*VCA2(2)
  78. &+DDEF0(5)*VCA1(1)*VCA2(3)+DDEF0(6)*VCA1(2)*VCA2(3)
  79. &+DDEF0(3)*VCA1(3)*VCA2(3)
  80. C
  81. DDEF(5)=DDEF0(1)*VCA1(1)*VCA3(1)+DDEF0(4)*VCA1(2)*VCA3(1)
  82. &+DDEF0(5)*VCA1(3)*VCA3(1)+DDEF0(4)*VCA1(1)*VCA3(2)
  83. &+DDEF0(2)*VCA1(2)*VCA3(2)+DDEF0(6)*VCA1(3)*VCA3(2)
  84. &+DDEF0(5)*VCA1(1)*VCA3(3)+DDEF0(6)*VCA1(2)*VCA3(3)
  85. &+DDEF0(3)*VCA1(3)*VCA3(3)
  86. C
  87. DDEF(6)=DDEF0(1)*VCA2(1)*VCA3(1)+DDEF0(4)*VCA2(2)*VCA3(1)
  88. &+DDEF0(5)*VCA2(3)*VCA3(1)+DDEF0(4)*VCA2(1)*VCA3(2)
  89. &+DDEF0(2)*VCA2(2)*VCA3(2)+DDEF0(6)*VCA2(3)*VCA3(2)
  90. &+DDEF0(5)*VCA2(1)*VCA3(3)+DDEF0(6)*VCA2(2)*VCA3(3)
  91. &+DDEF0(3)*VCA2(3)*VCA3(3)
  92. C
  93. ENDIF
  94. C
  95. C
  96. C=======================================================================
  97. C CALCUL DE SIGMA0 ET DDEF0 DANS LE REPERE GLOBAL
  98. C RESULTAT NOTE : SIGMA ET DDEF
  99. C=======================================================================
  100. C
  101. IF(ICAS.EQ.2) THEN
  102. C
  103. SIGMA(1)=SIGMA0(1)*VCA1(1)*VCA1(1)+SIGMA0(4)*VCA2(1)*VCA1(1)
  104. &+SIGMA0(5)*VCA3(1)*VCA1(1)+SIGMA0(4)*VCA1(1)*VCA2(1)
  105. &+SIGMA0(2)*VCA2(1)*VCA2(1)+SIGMA0(6)*VCA3(1)*VCA2(1)
  106. &+SIGMA0(5)*VCA1(1)*VCA3(1)+SIGMA0(6)*VCA2(1)*VCA3(1)
  107. &+SIGMA0(3)*VCA3(1)*VCA3(1)
  108. C
  109. SIGMA(2)=SIGMA0(1)*VCA1(2)*VCA1(2)+SIGMA0(4)*VCA2(2)*VCA1(2)
  110. &+SIGMA0(5)*VCA3(2)*VCA1(2)+SIGMA0(4)*VCA1(2)*VCA2(2)
  111. &+SIGMA0(2)*VCA2(2)*VCA2(2)+SIGMA0(6)*VCA3(2)*VCA2(2)
  112. &+SIGMA0(5)*VCA1(2)*VCA3(2)+SIGMA0(6)*VCA2(2)*VCA3(2)
  113. &+SIGMA0(3)*VCA3(2)*VCA3(2)
  114. C
  115. SIGMA(3)=SIGMA0(1)*VCA1(3)*VCA1(3)+SIGMA0(4)*VCA2(3)*VCA1(3)
  116. &+SIGMA0(5)*VCA3(3)*VCA1(3)+SIGMA0(4)*VCA1(3)*VCA2(3)
  117. &+SIGMA0(2)*VCA2(3)*VCA2(3)+SIGMA0(6)*VCA3(3)*VCA2(3)
  118. &+SIGMA0(5)*VCA1(3)*VCA3(3)+SIGMA0(6)*VCA2(3)*VCA3(3)
  119. &+SIGMA0(3)*VCA3(3)*VCA3(3)
  120. C
  121. SIGMA(4)=SIGMA0(1)*VCA1(1)*VCA1(2)+SIGMA0(4)*VCA2(1)*VCA1(2)
  122. &+SIGMA0(5)*VCA3(1)*VCA1(2)+SIGMA0(4)*VCA1(1)*VCA2(2)
  123. &+SIGMA0(2)*VCA2(1)*VCA2(2)+SIGMA0(6)*VCA3(1)*VCA2(2)
  124. &+SIGMA0(5)*VCA1(1)*VCA3(2)+SIGMA0(6)*VCA2(1)*VCA3(2)
  125. &+SIGMA0(3)*VCA3(1)*VCA3(2)
  126. C
  127. SIGMA(5)=SIGMA0(1)*VCA1(1)*VCA1(3)+SIGMA0(4)*VCA2(1)*VCA1(3)
  128. &+SIGMA0(5)*VCA3(1)*VCA1(3)+SIGMA0(4)*VCA1(1)*VCA2(3)
  129. &+SIGMA0(2)*VCA2(1)*VCA2(3)+SIGMA0(6)*VCA3(1)*VCA2(3)
  130. &+SIGMA0(5)*VCA1(1)*VCA3(3)+SIGMA0(6)*VCA2(1)*VCA3(3)
  131. &+SIGMA0(3)*VCA3(1)*VCA3(3)
  132. C
  133. SIGMA(6)=SIGMA0(1)*VCA1(2)*VCA1(3)+SIGMA0(4)*VCA2(2)*VCA1(3)
  134. &+SIGMA0(5)*VCA3(2)*VCA1(3)+SIGMA0(4)*VCA1(2)*VCA2(3)
  135. &+SIGMA0(2)*VCA2(2)*VCA2(3)+SIGMA0(6)*VCA3(2)*VCA2(3)
  136. &+SIGMA0(5)*VCA1(2)*VCA3(3)+SIGMA0(6)*VCA2(2)*VCA3(3)
  137. &+SIGMA0(3)*VCA3(2)*VCA3(3)
  138. C
  139. C
  140. DDEF(1)=DDEF0(1)*VCA1(1)*VCA1(1)+DDEF0(4)*VCA2(1)*VCA1(1)
  141. &+DDEF0(5)*VCA3(1)*VCA1(1)+DDEF0(4)*VCA1(1)*VCA2(1)
  142. &+DDEF0(2)*VCA2(1)*VCA2(1)+DDEF0(6)*VCA3(1)*VCA2(1)
  143. &+DDEF0(5)*VCA1(1)*VCA3(1)+DDEF0(6)*VCA2(1)*VCA3(1)
  144. &+DDEF0(3)*VCA3(1)*VCA3(1)
  145. C
  146. DDEF(2)=DDEF0(1)*VCA1(2)*VCA1(2)+DDEF0(4)*VCA2(2)*VCA1(2)
  147. &+DDEF0(5)*VCA3(2)*VCA1(2)+DDEF0(4)*VCA1(2)*VCA2(2)
  148. &+DDEF0(2)*VCA2(2)*VCA2(2)+DDEF0(6)*VCA3(2)*VCA2(2)
  149. &+DDEF0(5)*VCA1(2)*VCA3(2)+DDEF0(6)*VCA2(2)*VCA3(2)
  150. &+DDEF0(3)*VCA3(2)*VCA3(2)
  151. C
  152. DDEF(3)=DDEF0(1)*VCA1(3)*VCA1(3)+DDEF0(4)*VCA2(3)*VCA1(3)
  153. &+DDEF0(5)*VCA3(3)*VCA1(3)+DDEF0(4)*VCA1(3)*VCA2(3)
  154. &+DDEF0(2)*VCA2(3)*VCA2(3)+DDEF0(6)*VCA3(3)*VCA2(3)
  155. &+DDEF0(5)*VCA1(3)*VCA3(3)+DDEF0(6)*VCA2(3)*VCA3(3)
  156. &+DDEF0(3)*VCA3(3)*VCA3(3)
  157. C
  158. DDEF(4)=DDEF0(1)*VCA1(1)*VCA1(2)+DDEF0(4)*VCA2(1)*VCA1(2)
  159. &+DDEF0(5)*VCA3(1)*VCA1(2)+DDEF0(4)*VCA1(1)*VCA2(2)
  160. &+DDEF0(2)*VCA2(1)*VCA2(2)+DDEF0(6)*VCA3(1)*VCA2(2)
  161. &+DDEF0(5)*VCA1(1)*VCA3(2)+DDEF0(6)*VCA2(1)*VCA3(2)
  162. &+DDEF0(3)*VCA3(1)*VCA3(2)
  163. C
  164. DDEF(5)=DDEF0(1)*VCA1(1)*VCA1(3)+DDEF0(4)*VCA2(1)*VCA1(3)
  165. &+DDEF0(5)*VCA3(1)*VCA1(3)+DDEF0(4)*VCA1(1)*VCA2(3)
  166. &+DDEF0(2)*VCA2(1)*VCA2(3)+DDEF0(6)*VCA3(1)*VCA2(3)
  167. &+DDEF0(5)*VCA1(1)*VCA3(3)+DDEF0(6)*VCA2(1)*VCA3(3)
  168. &+DDEF0(3)*VCA3(1)*VCA3(3)
  169. C
  170. DDEF(6)=DDEF0(1)*VCA1(2)*VCA1(3)+DDEF0(4)*VCA2(2)*VCA1(3)
  171. &+DDEF0(5)*VCA3(2)*VCA1(3)+DDEF0(4)*VCA1(2)*VCA2(3)
  172. &+DDEF0(2)*VCA2(2)*VCA2(3)+DDEF0(6)*VCA3(2)*VCA2(3)
  173. &+DDEF0(5)*VCA1(2)*VCA3(3)+DDEF0(6)*VCA2(2)*VCA3(3)
  174. &+DDEF0(3)*VCA3(2)*VCA3(3)
  175. C
  176. ENDIF
  177. C
  178. RETURN
  179. END
  180.  
  181.  

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