Esope source files

   1 : C ACRILO    SOURCE    CHAT      05/01/12    21:16:10     5004           
   2 :       SUBROUTINE ACRILO(RP,XL,RAY,RHO,RHOR,C,CREF)                            
   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  SORTIE                                                                       
  16 : C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 4X4                                 
  17 : C                                                                               
  18 : C-----------------------------------------------------------------------        
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  21 :       DIMENSION RP(4,4)
  22 : C
  23 :       DO 1 I=1,16
  24 : 1     RP(I,1)= 0.
  25 : C
  26 :       PI=2*ASIN(1.)
  27 :       SE=PI*RAY*RAY
  28 :       A=(SE*XL*RHOR**2*CREF**4)/(3*RHO*C**2)
  29 : C
  30 :       RP(2,2)=A
  31 :       RP(2,4)=A/2
  32 :       RP(4,4)=A
  33 : C
  34 : C SYMETRISATION
  35 : C
  36 :       DO 2 I=1,4
  37 :       DO 2 J=I,4
  38 : 2     RP(J,I)=RP(I,J)
  39 : C
  40 :       RETURN
  41 :       END
  42 : 
  43 :                                                   
  44 :