Télécharger coumal.eso

Retour à la liste

Numérotation des lignes :

  1. C COUMAL SOURCE PV 11/03/08 21:15:17 6888
  2. SUBROUTINE COUMAL(RP,XL,RAY,RC,RHO,LCAR,RHOR)
  3. C-----------------------------------------------------------------------
  4. C
  5. C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D'UN
  6. C ELEMENT DE RACCORD LITU EN AXES LOCAUX
  7. C
  8. C ENTREE
  9. C XL LONGUEUR
  10. C RAY RAYON DU TUYAU FLUIDE
  11. C RC RAYON DE COURBURE
  12. C RHO MASSE VOLUMIQUE DU FLUIDE
  13. C LCAR LONGUEUR CARACTERISTIQUE
  14. C RHOR MASSE VOLUMIQUE DE REFERENCE
  15. C SORTIE
  16. C RP MATRICE DE MASSE ELEMENTAIRE 10X10
  17. C
  18. C-----------------------------------------------------------------------
  19. IMPLICIT INTEGER(I-N)
  20. IMPLICIT REAL*8(A-H,O-Z)
  21. DIMENSION RP(10,10)
  22. REAL*8 LCAR
  23. C
  24. DO 1 I=1,100
  25. 1 RP(I,1)=0.
  26. C
  27. PI=2*ASIN(1.)
  28. SE=PI*RAY*RAY
  29. A=RHO*SE*XL*0.3714285714285714D+00
  30. B=-RHOR*LCAR*LCAR*SE
  31. C=RHO*SE*XL*0.1285714285714286D+00
  32. C
  33. RP(2,2)=A
  34. RP(3,3)=A
  35. RP(2,7)=C
  36. RP(3,8)=C
  37. RP(7,7)=A
  38. RP(8,8)=A
  39. RP(1,4)=-B
  40. RP(6,9)=B
  41. C
  42. IF (RC.NE.0) THEN
  43. C=(RHOR*LCAR*LCAR*SE*XL)/(3*RC)
  44. RP(3,4)=C
  45. RP(3,9)=C/2
  46. RP(4,8)=RP(3,9)
  47. RP(8,9)=C
  48. ENDIF
  49. C
  50. C SYMETRISATION
  51. C
  52. DO 2 I=1,10
  53. DO 2 J=I,10
  54. 2 RP(J,I)=RP(I,J)
  55. C
  56. RETURN
  57. END
  58.  
  59.  
  60.  

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