Esope source files

   1 : C ACMALO    SOURCE    PV        07/11/23    21:15:09     5978           
   2 :       SUBROUTINE ACMALO(RP,XL,RAY,RHO,RHOR,C,CREF,LCAR)                            
   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    LCAR    LONGUEUR CARACTERISTIQUE                                                    
  16 : C  SORTIE                                                                       
  17 : C    RP      MATRICE DE MASSE ELEMENTAIRE 4X4                                 
  18 : C                                                                               
  19 : C-----------------------------------------------------------------------        
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  22 :       DIMENSION RP(4,4)
  23 :       REAL*8 LCAR
  24 : C
  25 :       DO 1 I=1,16
  26 : 1     RP(I,1)= 0.
  27 : C
  28 :       PI=2*ASIN(1.)
  29 :       SE=PI*RAY*RAY
  30 :       A=-(SE*RHOR**2*LCAR**4)/(RHO*XL)
  31 :       B=(SE*XL*RHOR**2*CREF**2*LCAR**2)/(3*RHO*C**2)
  32 : C
  33 :       RP(1,1)=A
  34 :       RP(1,2)=B
  35 :       RP(1,3)=-A
  36 :       RP(1,4)=B/2
  37 :       RP(2,3)=RP(1,4)
  38 :       RP(3,3)=A
  39 :       RP(3,4)=B
  40 : C
  41 : C SYMETRISATION
  42 : C
  43 :       DO 2 I=1,4
  44 :       DO 2 J=I,4
  45 : 2     RP(J,I)=RP(I,J)
  46 : C
  47 :       RETURN
  48 :       END
  49 : 
  50 :                                                   
  51 :  
  52 :