Télécharger la plaquette Cast3M


fr en

Télécharger pocrlo.eso

Retour à la liste

Numérotation des lignes : 

pocrlo

  1. C POCRLO    SOURCE    PV        18/06/18    21:15:19     9860           
  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 J=1,12
  33.       DO 120 I=1,12
  34.   120 RP(I,J)=0.D00
  35. C
  36. C
  37.       RHOSD = RHO*SD
  38.       FAC=RHOSD*XL/6.
  39.       FAC2=FAC*2.
  40.       ROME(1,1) = 0.
  41.       ROME(2,2) = 0.
  42.       ROME(3,3) = 0.
  43.       ROME(1,2) = (-1.)*VROTL(3)
  44.       ROME(1,3) = VROTL(2)
  45.       ROME(2,3) = (-1.)*VROTL(1)
  46.       ROME(2,1) = VROTL(3)
  47.       ROME(3,1) = (-1.)*VROTL(2)
  48.       ROME(3,2) = VROTL(1)
  49. C
  50. C  REMPLISSAGE DE LA MATRICE
  51. C  FX1
  52.       RP(1,1)=FAC2*ROME(1,1)
  53.       RP(1,2)=FAC2*ROME(1,2)
  54.       RP(1,3)=FAC2*ROME(1,3)
  55.       RP(1,7)=FAC*ROME(1,1)
  56.       RP(1,8)=FAC*ROME(1,2)
  57.       RP(1,9)=FAC*ROME(1,3)
  58. C  FY1
  59.       RP(2,1)=FAC2*ROME(2,1)
  60.       RP(2,2)=FAC2*ROME(2,2)
  61.       RP(2,3)=FAC2*ROME(2,3)
  62.       RP(2,7)=FAC*ROME(2,1)
  63.       RP(2,8)=FAC*ROME(2,2)
  64.       RP(2,9)=FAC*ROME(2,3)
  65. C  FZ1
  66.       RP(3,1)=FAC2*ROME(3,1)
  67.       RP(3,2)=FAC2*ROME(3,2)
  68.       RP(3,3)=FAC2*ROME(3,3)
  69.       RP(3,7)=FAC*ROME(3,1)
  70.       RP(3,5)=FAC*ROME(3,2)
  71.       RP(3,9)=FAC*ROME(3,3)
  72. C  FX2
  73.       RP(7,1)=FAC*ROME(1,1)
  74.       RP(7,2)=FAC*ROME(1,2)
  75.       RP(7,3)=FAC*ROME(1,3)
  76.       RP(7,7)=FAC2*ROME(1,1)
  77.       RP(7,8)=FAC2*ROME(1,2)
  78.       RP(7,9)=FAC2*ROME(1,3)
  79. C  FY2
  80.       RP(8,1)=FAC*ROME(2,1)
  81.       RP(8,2)=FAC*ROME(2,2)
  82.       RP(8,3)=FAC*ROME(2,3)
  83.       RP(8,7)=FAC2*ROME(2,1)
  84.       RP(8,8)=FAC2*ROME(2,2)
  85.       RP(8,9)=FAC2*ROME(2,3)
  86. C  FZ2
  87.       RP(9,1)=FAC*ROME(3,1)
  88.       RP(9,2)=FAC*ROME(3,2)
  89.       RP(9,3)=FAC*ROME(3,3)
  90.       RP(9,7)=FAC2*ROME(3,1)
  91.       RP(9,8)=FAC2*ROME(3,2)
  92.       RP(9,9)=FAC2*ROME(3,3)
  93. C
  94.       RETURN
  95.       END
  96.  
  97.  
  98.  
  99.  
  100.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales