C ACMALO SOURCE PV 18/06/19 21:15:01 9861 SUBROUTINE ACMALO(RP,XL,RAY,RHO,RHOR,C,CREF,LCAR,LREL) C----------------------------------------------------------------------- C C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D'UN C TUYAU ACOUSTIQUE PUR EN AXES LOCAUX C C ENTREE C XL LONGUEUR C RAY RAYON DU TUYAU FLUIDE C RHO MASSE VOLUMIQUE DU FLUIDE C RHOR MASSE VOLUMIQUE DE REFERENCE C C CELERITE DU FLUIDE C CREF CELERITE DE REFERENCE C LCAR LONGUEUR CARACTERISTIQUE C SORTIE C RP MATRICE DE MASSE ELEMENTAIRE 4X4 C C----------------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) DIMENSION RP(LREL,*) REAL*8 LCAR C DO 1 J=1,LREL DO 1 I=1,LREL 1 RP(I,J)= 0. C PI=2*ASIN(1.) SE=PI*RAY*RAY A=-(SE*RHOR**2*LCAR**4)/(RHO*XL) B=(SE*XL*RHOR**2*CREF**2*LCAR**2)/(3*RHO*C**2) C RP(1,1)=A RP(1,2)=B RP(1,3)=-A RP(1,4)=B/2 RP(2,3)=RP(1,4) RP(3,3)=A RP(3,4)=B C C SYMETRISATION C DO 2 I=1,4 DO 2 J=I,4 2 RP(J,I)=RP(I,J) C RETURN END