Télécharger pomal2.eso

Retour à la liste

Numérotation des lignes :

  1. C POMAL2 SOURCE CHAT 05/01/13 02:17:24 5004
  2. C POMAL2 SOURCE CASTEM 85/06/13 12:33:23
  3. SUBROUTINE POMAL2(RP,XL,RHO,SD,XINZ)
  4. C-----------------------------------------------------------------------
  5. C
  6. C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE
  7. C POUTRE DROITE EN AXES LOCAUX
  8. C D APRES TEDEL
  9. C DANS CETTE VERSION,ON PREND EN COMPTE LES INERTIES DE ROTATION
  10. C MAIS PAS LA DEFORMABILITE A L EFFORT TRANCHANT
  11. C
  12. C ENTREE
  13. C XL LONGUEUR
  14. C SD SECTION DROITE
  15. C XINZ INERTIE IZ
  16. C RHO MASSE VOLUMIQUE
  17. C SORTIE
  18. C RP MATRICE DE MASSE ELEMENTAIRE 12X12
  19. C
  20. C VIENT DE TEDEL AOUT 84
  21. C
  22. C-----------------------------------------------------------------------
  23. IMPLICIT INTEGER(I-N)
  24. IMPLICIT REAL*8(A-H,O-Z)
  25. DIMENSION RP(6,6)
  26. C
  27. C INITIALISATIONS
  28. C
  29. DO 120 I=1,36
  30. 120 RP(I,1)=0.D00
  31. RL=1.D00/XL
  32. XL2=XL*XL
  33. RHOSL=RHO*SD*XL
  34. RSP1= RHOSL*0.333333333333333333
  35. RSP2= RHOSL*0.166666666666666666
  36. U 1= 0.3714285714285714D+00 * RHOSL
  37. U 2= 0.5238095238095238D-01 * RHOSL *XL
  38. U 3= 0.1285714285714286D+00 * RHOSL
  39. U 4= -0.3095238095238095D-01 * RHOSL *XL
  40. U 5= 0.9523809523809524D-02 * RHOSL *XL2
  41. U 6= 0.3095238095238095D-01 * RHOSL *XL
  42. U 7= -0.7142857142857143D-02 * RHOSL *XL2
  43. U 8= 0.3714285714285714D+00 * RHOSL
  44. U 9= -0.5238095238095238D-01 * RHOSL *XL
  45. U10= 0.9523809523809524D-02 * RHOSL *XL2
  46. U11=1.2*RHO*RL*XINZ
  47. U12=0.1*RHO*XINZ
  48. U18=0.13333333333333333*XINZ*RHO*XL
  49. U19=0.1*XINZ*RHO
  50. U20=0.03333333333333333*XINZ*XL*RHO
  51. C
  52. RP(1,1)=RSP1
  53. RP(1,4)=RSP2
  54. RP(2,2)=U1+U11
  55. RP(2,3)=U2+U12
  56. RP(2,5)=U3-U11
  57. RP(2,6)=U4+U12
  58. RP(3,3)=U5+U18
  59. RP(3,5)=U6-U19
  60. RP(3,6)=U7-U20
  61. RP(4,4)=RSP1
  62. RP(5,5)=U8+U11
  63. RP(5,6)=U9-U12
  64. RP(6,6)=U10+U18
  65. DO 9 I=1,6
  66. DO 9 J=I,6
  67. 9 RP(J,I)= RP(I,J)
  68. RETURN
  69. END
  70.  
  71.  
  72.  

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