Esope source files

   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 :