Esope source files

   1 : C POGYLO    SOURCE    CHAT      05/01/13    02:16:29     5004           
   2 :       SUBROUTINE POGYLO(RP,XL,RHO,SD,XINX,XINY,XINZ)
   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    D APRES TEDEL et ROTOR
   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 : C-----------------------------------------------------------------------
  25 :       IMPLICIT INTEGER(I-N)
  26 :       IMPLICIT REAL*8(A-H,O-Z)
  27 :       DIMENSION RP(12,12)
  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 :       RHOGI = RHO*XINX
  36 :       RG1P = (RHOGI*6.D0)/(5.D0*XL)
  37 :       RG1M = (RHOGI*(-6.D0))/(5.D0*XL)
  38 :       RG2P = (RHOGI)/(10.D0)
  39 :       RG2M = ((-1.D0)*RHOGI)/(10.D0)
  40 :       RG3P = (RHOGI*XL)/(30.D0)
  41 :       RG3M = ((-1.D0)*RHOGI*XL)/(30.D0)
  42 :       RG4P = ((2.D0)*RHOGI*XL)/(15.D0)
  43 :       RG4M = (-(2.D0)*RHOGI*XL)/(15.D0)
  44 : C
  45 : C  REMPLISSAGE DE LA MATRICE
  46 : C  FY1
  47 :       RP(2,3)=RG1P
  48 :       RP(2,5)=RG2M
  49 :       RP(2,9)=RG1M
  50 :       RP(2,11)=RG2M
  51 : C  FZ1
  52 :       RP(3,2)=RG1M
  53 :       RP(3,6)=RG2M
  54 :       RP(3,8)=RG1P
  55 :       RP(3,12)=RG2M
  56 : C  MY1
  57 :       RP(5,2)=RG2P
  58 :       RP(5,6)=RG4P
  59 :       RP(5,8)=RG2M
  60 :       RP(5,12)=RG3M
  61 : C  MZ1
  62 :       RP(6,3)=RG2P
  63 :       RP(6,5)=RG4M
  64 :       RP(6,9)=RG2M
  65 :       RP(6,11)=RG3P
  66 : C  FY2
  67 :       RP(8,3)=RG1M
  68 :       RP(8,5)=RG2P
  69 :       RP(8,9)=RG1P
  70 :       RP(8,11)=RG2P
  71 : C  FZ2
  72 :       RP(9,2)=RG1P
  73 :       RP(9,6)=RG2P
  74 :       RP(9,8)=RG1M
  75 :       RP(9,12)=RG2P
  76 : C  MY2
  77 :       RP(11,2)=RG2P
  78 :       RP(11,6)=RG3M
  79 :       RP(11,8)=RG2M
  80 :       RP(11,12)=RG4P
  81 : C  MZ2
  82 :       RP(12,3)=RG2P
  83 :       RP(12,5)=RG3P
  84 :       RP(12,9)=RG2M
  85 :       RP(12,11)=RG4M
  86 : C
  87 :       RETURN
  88 :       END
  89 :  
  90 :