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 :
© Cast3M 2003 - All rights reserved.
Disclaimer