Download uo2rot.eso

Back to the list

   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 :  

© Cast3M 2003 - All rights reserved.
Disclaimer