Télécharger pocrlo.eso

Retour à la liste

Numérotation des lignes :

  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.  
  98.  

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