Télécharger pomalo.eso

Retour à la liste

Numérotation des lignes :

  1. C POMALO SOURCE PV 18/06/18 21:15:25 9860
  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 J=1,12
  31. DO 120 I=1,12
  32. 120 RP(I,J)=0.D00
  33. RL=1.D00/XL
  34. XL2=XL*XL
  35. RHOSL=RHO*SD*XL
  36. RSP1= RHOSL*0.333333333333333333
  37. RSP2= RHOSL*0.166666666666666666
  38. TOR1= RHO*XINX*XL*0.333333333333333333
  39. TOR2= RHO*XINX*XL*0.166666666666666666
  40. U 1= 0.3714285714285714D+00 * RHOSL
  41. U 2= 0.5238095238095238D-01 * RHOSL *XL
  42. U 3= 0.1285714285714286D+00 * RHOSL
  43. U 4= -0.3095238095238095D-01 * RHOSL *XL
  44. U 5= 0.9523809523809524D-02 * RHOSL *XL2
  45. U 6= 0.3095238095238095D-01 * RHOSL *XL
  46. U 7= -0.7142857142857143D-02 * RHOSL *XL2
  47. U 8= 0.3714285714285714D+00 * RHOSL
  48. U 9= -0.5238095238095238D-01 * RHOSL *XL
  49. U10= 0.9523809523809524D-02 * RHOSL *XL2
  50. U11=1.2*RHO*RL*XINZ
  51. U12=0.1*RHO*XINZ
  52. U13=1.2*RHO*RL*XINY
  53. U14=0.1*RHO*XINY
  54. U15=0.13333333333333333*XINY*RHO*XL
  55. U16=0.1*XINY*RHO
  56. U17=0.03333333333333333*XINY*XL*RHO
  57. U18=0.13333333333333333*XINZ*RHO*XL
  58. U19=0.1*XINZ*RHO
  59. U20=0.03333333333333333*XINZ*XL*RHO
  60. C
  61. RP(1,1)=RSP1
  62. RP(1,7)=RSP2
  63. RP(2,2)=U1+U11
  64. RP(2,6)=U2+U12
  65. RP(2,8)=U3-U11
  66. RP(2,12)=U4+U12
  67. RP(3,3)=U1+U13
  68. RP(3,5)=-U2-U14
  69. RP(3,9)=U3-U13
  70. RP(3,11)=-U4-U14
  71. RP(4,4)= TOR1
  72. RP(4,10)= TOR2
  73. RP(5,5)=U5+U15
  74. RP(5,9)=-U6+U16
  75. RP(5,11)=U7-U17
  76. RP(6,6)=U5+U18
  77. RP(6,8)=U6-U19
  78. RP(6,12)=U7-U20
  79. RP(7,7)=RSP1
  80. RP(8,8)=U8+U11
  81. RP(8,12)=U9-U12
  82. RP(9,9)=U8+U13
  83. RP(9,11)=-U9+U14
  84. RP(10,10)=TOR1
  85. RP(11,11)=U10+U15
  86. RP(12,12)=U10+U18
  87. DO 9 I=1,12
  88. DO 9 J=I,12
  89. 9 RP(J,I)= RP(I,J)
  90. RETURN
  91. END
  92.  
  93.  
  94.  

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