Télécharger la plaquette Cast3M


fr en

Télécharger ottrot.eso

Retour à la liste

Numérotation des lignes : 

ottrot

  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