Esope source files

   1 : C POCRLO    SOURCE    CHAT      05/01/13    02:15:57     5004           
   2 :        SUBROUTINE POCRLO(RP,XL,RHO,SD,XINX,XINY,XINZ,VROTL)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C    ROUTINE DE CALCUL DE LA MATRICE DE COUPLAGE
   6 : C    GYROSCOPIQUE ELEMENTAIRE D UNE
   7 : C    POUTRE DROITE EN AXES LOCAUX
   8 : C    (REPERE TOURNANT)
   9 : C
  10 : C
  11 : C  ENTREE
  12 : C    XL      LONGUEUR
  13 : C    SD      SECTION DROITE
  14 : C    XINX    INERTIE IX
  15 : C    XINY    INERTIE IY
  16 : C    XINZ    INERTIE IZ
  17 : C    RHO     MASSE VOLUMIQUE
  18 : C  SORTIE
  19 : C    RP      MATRICE DE COUPLAGE GYROSCOPIQUE ELEMENTAIRE 12X12
  20 : C
  21 : C    Didier COMBESCURE Mars 2003
  22 : C
  23 : C-----------------------------------------------------------------------
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8(A-H,O-Z)
  26 : -INC CCREEL     
  27 :       DIMENSION RP(12,12),VROTL(*),ROME(3,3)
  28 : C
  29 : C  INITIALISATIONS
  30 : C
  31 : C
  32 :       DO 120 I=1,144
  33 :   120 RP(I,1)=0.D00
  34 : C 
  35 : C
  36 :       RHOSD = RHO*SD
  37 :       FAC=RHOSD*XL/6.
  38 :       FAC2=FAC*2. 
  39 :       ROME(1,1) = 0.
  40 :       ROME(2,2) = 0.
  41 :       ROME(3,3) = 0.
  42 :       ROME(1,2) = (-1.)*VROTL(3)
  43 :       ROME(1,3) = VROTL(2)
  44 :       ROME(2,3) = (-1.)*VROTL(1)
  45 :       ROME(2,1) = VROTL(3)
  46 :       ROME(3,1) = (-1.)*VROTL(2)
  47 :       ROME(3,2) = VROTL(1)
  48 : C
  49 : C  REMPLISSAGE DE LA MATRICE
  50 : C  FX1
  51 :       RP(1,1)=FAC2*ROME(1,1)
  52 :       RP(1,2)=FAC2*ROME(1,2)
  53 :       RP(1,3)=FAC2*ROME(1,3)
  54 :       RP(1,7)=FAC*ROME(1,1)
  55 :       RP(1,8)=FAC*ROME(1,2)
  56 :       RP(1,9)=FAC*ROME(1,3)
  57 : C  FY1
  58 :       RP(2,1)=FAC2*ROME(2,1)
  59 :       RP(2,2)=FAC2*ROME(2,2)
  60 :       RP(2,3)=FAC2*ROME(2,3)
  61 :       RP(2,7)=FAC*ROME(2,1)
  62 :       RP(2,8)=FAC*ROME(2,2)
  63 :       RP(2,9)=FAC*ROME(2,3)
  64 : C  FZ1
  65 :       RP(3,1)=FAC2*ROME(3,1)
  66 :       RP(3,2)=FAC2*ROME(3,2)
  67 :       RP(3,3)=FAC2*ROME(3,3)
  68 :       RP(3,7)=FAC*ROME(3,1)
  69 :       RP(3,5)=FAC*ROME(3,2)
  70 :       RP(3,9)=FAC*ROME(3,3)
  71 : C  FX2
  72 :       RP(7,1)=FAC*ROME(1,1)
  73 :       RP(7,2)=FAC*ROME(1,2)
  74 :       RP(7,3)=FAC*ROME(1,3)
  75 :       RP(7,7)=FAC2*ROME(1,1)
  76 :       RP(7,8)=FAC2*ROME(1,2)
  77 :       RP(7,9)=FAC2*ROME(1,3)
  78 : C  FY2
  79 :       RP(8,1)=FAC*ROME(2,1)
  80 :       RP(8,2)=FAC*ROME(2,2)
  81 :       RP(8,3)=FAC*ROME(2,3)
  82 :       RP(8,7)=FAC2*ROME(2,1)
  83 :       RP(8,8)=FAC2*ROME(2,2)
  84 :       RP(8,9)=FAC2*ROME(2,3)
  85 : C  FZ2
  86 :       RP(9,1)=FAC*ROME(3,1)
  87 :       RP(9,2)=FAC*ROME(3,2)
  88 :       RP(9,3)=FAC*ROME(3,3)
  89 :       RP(9,7)=FAC2*ROME(3,1)
  90 :       RP(9,8)=FAC2*ROME(3,2)
  91 :       RP(9,9)=FAC2*ROME(3,3)
  92 : C
  93 :       RETURN
  94 :       END
  95 :  
  96 :  
  97 :