Télécharger acmalo.eso

Retour à la liste

Numérotation des lignes :

  1. C ACMALO SOURCE PV 18/06/19 21:15:01 9861
  2. SUBROUTINE ACMALO(RP,XL,RAY,RHO,RHOR,C,CREF,LCAR,LREL)
  3. C-----------------------------------------------------------------------
  4. C
  5. C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D'UN
  6. C TUYAU ACOUSTIQUE PUR EN AXES LOCAUX
  7. C
  8. C ENTREE
  9. C XL LONGUEUR
  10. C RAY RAYON DU TUYAU FLUIDE
  11. C RHO MASSE VOLUMIQUE DU FLUIDE
  12. C RHOR MASSE VOLUMIQUE DE REFERENCE
  13. C C CELERITE DU FLUIDE
  14. C CREF CELERITE DE REFERENCE
  15. C LCAR LONGUEUR CARACTERISTIQUE
  16. C SORTIE
  17. C RP MATRICE DE MASSE ELEMENTAIRE 4X4
  18. C
  19. C-----------------------------------------------------------------------
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22. DIMENSION RP(LREL,*)
  23. REAL*8 LCAR
  24. C
  25. DO 1 J=1,LREL
  26. DO 1 I=1,LREL
  27. 1 RP(I,J)= 0.
  28. C
  29. PI=2*ASIN(1.)
  30. SE=PI*RAY*RAY
  31. A=-(SE*RHOR**2*LCAR**4)/(RHO*XL)
  32. B=(SE*XL*RHOR**2*CREF**2*LCAR**2)/(3*RHO*C**2)
  33. C
  34. RP(1,1)=A
  35. RP(1,2)=B
  36. RP(1,3)=-A
  37. RP(1,4)=B/2
  38. RP(2,3)=RP(1,4)
  39. RP(3,3)=A
  40. RP(3,4)=B
  41. C
  42. C SYMETRISATION
  43. C
  44. DO 2 I=1,4
  45. DO 2 J=I,4
  46. 2 RP(J,I)=RP(I,J)
  47. C
  48. RETURN
  49. END
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  

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