Télécharger pomalo.eso

Retour à la liste

Numérotation des lignes :

  1. C POMALO SOURCE CHAT 05/01/13 02:17:28 5004
  2. SUBROUTINE POMALO(RP,XL,RHO,SD,XINX,XINY,XINZ)
  3. C-----------------------------------------------------------------------
  4. C
  5. C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE
  6. C POUTRE DROITE EN AXES LOCAUX
  7. C D APRES TEDEL
  8. C DANS CETTE VERSION,ON PREND EN COMPTE LES INERTIES DE ROTATION
  9. C MAIS PAS LA DEFORMABILITE A L EFFORT TRANCHANT
  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 MASSE ELEMENTAIRE 12X12
  20. C
  21. C VIENT DE TEDEL AOUT 84
  22. C
  23. C-----------------------------------------------------------------------
  24. IMPLICIT INTEGER(I-N)
  25. IMPLICIT REAL*8(A-H,O-Z)
  26. DIMENSION RP(12,12)
  27. C
  28. C INITIALISATIONS
  29. C
  30. DO 120 I=1,144
  31. 120 RP(I,1)=0.D00
  32. RL=1.D00/XL
  33. XL2=XL*XL
  34. RHOSL=RHO*SD*XL
  35. RSP1= RHOSL*0.333333333333333333
  36. RSP2= RHOSL*0.166666666666666666
  37. TOR1= RHO*XINX*XL*0.333333333333333333
  38. TOR2= RHO*XINX*XL*0.166666666666666666
  39. U 1= 0.3714285714285714D+00 * RHOSL
  40. U 2= 0.5238095238095238D-01 * RHOSL *XL
  41. U 3= 0.1285714285714286D+00 * RHOSL
  42. U 4= -0.3095238095238095D-01 * RHOSL *XL
  43. U 5= 0.9523809523809524D-02 * RHOSL *XL2
  44. U 6= 0.3095238095238095D-01 * RHOSL *XL
  45. U 7= -0.7142857142857143D-02 * RHOSL *XL2
  46. U 8= 0.3714285714285714D+00 * RHOSL
  47. U 9= -0.5238095238095238D-01 * RHOSL *XL
  48. U10= 0.9523809523809524D-02 * RHOSL *XL2
  49. U11=1.2*RHO*RL*XINZ
  50. U12=0.1*RHO*XINZ
  51. U13=1.2*RHO*RL*XINY
  52. U14=0.1*RHO*XINY
  53. U15=0.13333333333333333*XINY*RHO*XL
  54. U16=0.1*XINY*RHO
  55. U17=0.03333333333333333*XINY*XL*RHO
  56. U18=0.13333333333333333*XINZ*RHO*XL
  57. U19=0.1*XINZ*RHO
  58. U20=0.03333333333333333*XINZ*XL*RHO
  59. C
  60. RP(1,1)=RSP1
  61. RP(1,7)=RSP2
  62. RP(2,2)=U1+U11
  63. RP(2,6)=U2+U12
  64. RP(2,8)=U3-U11
  65. RP(2,12)=U4+U12
  66. RP(3,3)=U1+U13
  67. RP(3,5)=-U2-U14
  68. RP(3,9)=U3-U13
  69. RP(3,11)=-U4-U14
  70. RP(4,4)= TOR1
  71. RP(4,10)= TOR2
  72. RP(5,5)=U5+U15
  73. RP(5,9)=-U6+U16
  74. RP(5,11)=U7-U17
  75. RP(6,6)=U5+U18
  76. RP(6,8)=U6-U19
  77. RP(6,12)=U7-U20
  78. RP(7,7)=RSP1
  79. RP(8,8)=U8+U11
  80. RP(8,12)=U9-U12
  81. RP(9,9)=U8+U13
  82. RP(9,11)=-U9+U14
  83. RP(10,10)=TOR1
  84. RP(11,11)=U10+U15
  85. RP(12,12)=U10+U18
  86. DO 9 I=1,12
  87. DO 9 J=I,12
  88. 9 RP(J,I)= RP(I,J)
  89. RETURN
  90. END
  91.  
  92.  

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