Télécharger pokclo.eso

Retour à la liste

Numérotation des lignes :

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

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