Télécharger coumal.eso

Retour à la liste

Numérotation des lignes :

  1. C COUMAL SOURCE CB215821 18/06/19 21:15:04 9862
  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.  
  22. -INC CCREEL
  23.  
  24. DIMENSION RP(10,10)
  25. REAL*8 LCAR
  26. C
  27. DO J=1,10
  28. DO I=J,10
  29. RP(I,J)=0.D0
  30. ENDDO
  31. ENDDO
  32. C
  33. SE=XPI*RAY*RAY
  34. A=RHO*SE*XL*0.3714285714285714D0
  35. B=-RHOR*LCAR*LCAR*SE
  36. C=RHO*SE*XL*0.1285714285714286D0
  37. C
  38. RP(4,1) =-B
  39. RP(2,2) = A
  40. RP(3,3) = A
  41. RP(7,2) = C
  42. RP(8,3) = C
  43. RP(9,6) = B
  44. RP(7,7) = A
  45. RP(8,8) = A
  46. C
  47. IF (RC.NE.0) THEN
  48. C =(RHOR*LCAR*LCAR*SE*XL)/(3.D0*RC)
  49. CS2=C/2.D0
  50. RP(4,3) = C
  51. RP(9,3) = CS2
  52. RP(8,4) = CS2
  53. RP(9,8) = C
  54. ENDIF
  55. C
  56. C SYMETRISATION
  57. C
  58. DO I=1,10-1
  59. DO J=I+1,10
  60. RP(I,J)=RP(J,I)
  61. ENDDO
  62. ENDDO
  63. RETURN
  64. END
  65.  
  66.  

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