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

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