1 : C TIKCLO SOURCE CHAT 05/01/13 03:39:47 5004 2 : SUBROUTINE TIKCLO(RP,XL,RHO,SD,XINX,XINY,XINZ,VROTL) 3 : C----------------------------------------------------------------------- 4 : C 5 : C ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR 6 : C CENTRIFUGE ELEMENTAIRE D UNE 7 : C POUTRE DE TIMOSHENKO EN AXES LOCAUX 8 : C 9 : C 10 : C ENTREE 11 : C XL LONGUEUR 12 : C SD SECTION DROITE 13 : C XINX INERTIE IX 14 : C XINY INERTIE IY 15 : C XINZ INERTIE IZ 16 : C RHO MASSE VOLUMIQUE 17 : C SORTIE 18 : C RP MATRICE DE RAIDEUR CENTRIFUGE ELEMENTAIRE 12X12 19 : C 20 : C Didier COMBESCURE mars 2003 21 : C 22 : C----------------------------------------------------------------------- 23 : IMPLICIT INTEGER(I-N) 24 : IMPLICIT REAL*8(A-H,O-Z) 25 : -INC CCREEL 26 : DIMENSION RP(12,12),VROTL(3),ROME(3,3) 27 : C 28 : C INITIALISATIONS 29 : C 30 : C 31 : DO 120 I=1,144 32 : 120 RP(I,1)=0.D00 33 : C 34 : RHOSD = RHO*SD 35 : FAC=RHOSD*XL/6.D0 36 : FAC2=FAC*2.D0 37 : TOR1= RHO*XINX*XL/3.D0 38 : TOR2= RHO*XINX*XL/6.D0 39 : FLZ1= RHO*XINY*XL/3.D0 40 : FLZ2= RHO*XINY*XL/6.D0 41 : FLY1= RHO*XINZ*XL/3.D0 42 : FLY2= RHO*XINZ*XL/6.D0 43 : C 44 : C 45 : ROME(1,1) = (-1.)*((VROTL(2)**2) + (VROTL(3)**2)) 46 : ROME(2,2) = (-1.)*((VROTL(1)**2) + (VROTL(3)**2)) 47 : ROME(3,3) = (-1.)*((VROTL(1)**2) + (VROTL(2)**2)) 48 : ROME(1,2) = VROTL(1)*VROTL(2) 49 : ROME(1,3) = VROTL(1)*VROTL(3) 50 : ROME(2,3) = VROTL(2)*VROTL(3) 51 : ROME(2,1) = ROME(1,2) 52 : ROME(3,1) = ROME(1,3) 53 : ROME(3,2) = ROME(2,3) 54 : C 55 : C REMPLISSAGE DE LA MATRICE 56 : C FX1 57 : RP(1,1)=FAC2*ROME(1,1) 58 : RP(1,2)=FAC2*ROME(1,2) 59 : RP(1,3)=FAC2*ROME(1,3) 60 : RP(1,7)=FAC*ROME(1,1) 61 : RP(1,8)=FAC*ROME(1,2) 62 : RP(1,9)=FAC*ROME(1,3) 63 : C FY1 64 : RP(2,1)=FAC2*ROME(2,1) 65 : RP(2,2)=FAC2*ROME(2,2) 66 : RP(2,3)=FAC2*ROME(2,3) 67 : RP(2,7)=FAC*ROME(2,1) 68 : RP(2,8)=FAC*ROME(2,2) 69 : RP(2,9)=FAC*ROME(2,3) 70 : C FZ1 71 : RP(3,1)=FAC2*ROME(3,1) 72 : RP(3,2)=FAC2*ROME(3,2) 73 : RP(3,3)=FAC2*ROME(3,3) 74 : RP(3,7)=FAC*ROME(3,1) 75 : RP(3,5)=FAC*ROME(3,2) 76 : RP(3,9)=FAC*ROME(3,3) 77 : C FX2 78 : RP(7,1)=FAC*ROME(1,1) 79 : RP(7,2)=FAC*ROME(1,2) 80 : RP(7,3)=FAC*ROME(1,3) 81 : RP(7,7)=FAC2*ROME(1,1) 82 : RP(7,8)=FAC2*ROME(1,2) 83 : RP(7,9)=FAC2*ROME(1,3) 84 : C FY2 85 : RP(8,1)=FAC*ROME(2,1) 86 : RP(8,2)=FAC*ROME(2,2) 87 : RP(8,3)=FAC*ROME(2,3) 88 : RP(8,7)=FAC2*ROME(2,1) 89 : RP(8,8)=FAC2*ROME(2,2) 90 : RP(8,9)=FAC2*ROME(2,3) 91 : C FZ2 92 : RP(9,1)=FAC2*ROME(3,1) 93 : RP(9,2)=FAC2*ROME(3,2) 94 : RP(9,3)=FAC2*ROME(3,3) 95 : RP(9,7)=FAC*ROME(3,1) 96 : RP(9,8)=FAC*ROME(3,2) 97 : RP(9,9)=FAC*ROME(3,3) 98 : C 99 : C On rajoute le moment 100 : C 101 : C MX1 102 : C RP(4,4)=TOR1*ROME(1,1) 103 : C RP(4,5)=TOR1*ROME(1,2) 104 : C RP(4,6)=TOR1*ROME(1,3) 105 : C RP(4,10)=TOR2*ROME(1,1) 106 : C RP(4,11)=TOR2*ROME(1,2) 107 : C RP(4,12)=TOR2*ROME(1,3) 108 : C MY1 109 : C RP(5,4)=FLY1*ROME(2,1) 110 : C RP(5,5)=FLY1*ROME(2,2) 111 : C RP(5,6)=FLY1*ROME(2,3) 112 : C RP(5,10)=FLY2*ROME(2,1) 113 : C RP(5,11)=FLY2*ROME(2,2) 114 : C RP(5,12)=FLY2*ROME(2,3) 115 : C MZ1 116 : C RP(6,4)=FLZ1*ROME(3,1) 117 : C RP(6,5)=FLZ1*ROME(3,2) 118 : C RP(6,6)=FLZ1*ROME(3,3) 119 : C RP(6,10)=FLZ2*ROME(3,1) 120 : C RP(6,11)=FLZ2*ROME(3,2) 121 : C RP(6,12)=FLZ2*ROME(3,3) 122 : C MX2 123 : C RP(10,10)=TOR1*ROME(1,1) 124 : C RP(10,11)=TOR1*ROME(1,2) 125 : C RP(10,12)=TOR1*ROME(1,3) 126 : C RP(10,4)=TOR2*ROME(1,1) 127 : C RP(10,5)=TOR2*ROME(1,2) 128 : C RP(10,6)=TOR2*ROME(1,3) 129 : C MY2 130 : C RP(11,4)=FLY2*ROME(2,1) 131 : C RP(11,5)=FLY2*ROME(2,2) 132 : C RP(11,6)=FLY2*ROME(2,3) 133 : C RP(11,10)=FLY1*ROME(2,1) 134 : C RP(11,11)=FLY1*ROME(2,2) 135 : C RP(11,12)=FLY1*ROME(2,3) 136 : C MZ2 137 : C RP(12,4)=FLZ2*ROME(3,1) 138 : C RP(12,5)=FLZ2*ROME(3,2) 139 : C RP(12,6)=FLZ2*ROME(3,3) 140 : C RP(12,10)=FLZ1*ROME(3,1) 141 : C RP(12,11)=FLZ1*ROME(3,2) 142 : C RP(12,12)=FLZ1*ROME(3,3) 143 : C 144 : C 145 : RETURN 146 : END 147 : 148 : 149 :
© Cast3M 2003 - All rights reserved.
Disclaimer