Télécharger la plaquette Cast3M


fr en

Télécharger coumal.eso

Retour à la liste

Numérotation des lignes : 

coumal

  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