Télécharger pomal2.eso

Retour à la liste

Numérotation des lignes :

pomal2
  1. C POMAL2 SOURCE PV 18/06/18 21:15:24 9860
  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 J=1,6
  30. DO 120 I=1,6
  31. 120 RP(I,J)=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. U 1= 0.3714285714285714D+00 * RHOSL
  38. U 2= 0.5238095238095238D-01 * RHOSL *XL
  39. U 3= 0.1285714285714286D+00 * RHOSL
  40. U 4= -0.3095238095238095D-01 * RHOSL *XL
  41. U 5= 0.9523809523809524D-02 * RHOSL *XL2
  42. U 6= 0.3095238095238095D-01 * RHOSL *XL
  43. U 7= -0.7142857142857143D-02 * RHOSL *XL2
  44. U 8= 0.3714285714285714D+00 * RHOSL
  45. U 9= -0.5238095238095238D-01 * RHOSL *XL
  46. U10= 0.9523809523809524D-02 * RHOSL *XL2
  47. U11=1.2*RHO*RL*XINZ
  48. U12=0.1*RHO*XINZ
  49. U18=0.13333333333333333*XINZ*RHO*XL
  50. U19=0.1*XINZ*RHO
  51. U20=0.03333333333333333*XINZ*XL*RHO
  52. C
  53. RP(1,1)=RSP1
  54. RP(1,4)=RSP2
  55. RP(2,2)=U1+U11
  56. RP(2,3)=U2+U12
  57. RP(2,5)=U3-U11
  58. RP(2,6)=U4+U12
  59. RP(3,3)=U5+U18
  60. RP(3,5)=U6-U19
  61. RP(3,6)=U7-U20
  62. RP(4,4)=RSP1
  63. RP(5,5)=U8+U11
  64. RP(5,6)=U9-U12
  65. RP(6,6)=U10+U18
  66. DO 9 I=1,6
  67. DO 9 J=I,6
  68. 9 RP(J,I)= RP(I,J)
  69. RETURN
  70. END
  71.  
  72.  
  73.  
  74.  

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