Numérotation des lignes :

acomas
C ACOMAS    SOURCE    PV        18/06/18    21:15:01     9860                 SUBROUTINE ACOMAS(REL,LRE,CARAC,COORD,KERRE)C-----------------------------------------------------------------------C  ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE DE L ELEMENT D'C  ACOUSTIQUE PURE LSE2CC  ENTREEC          LRE     DIMENSION DU TABLEAU RELC          CARAC   TABLEAU DE CARACTERISTIQUES GEOMETRIQUES ETC                  MATERIELLES DU TUYAU FLUIDE  ( 7 VALEURS )C          COORD   COORDONNEES DES 2 NOEUDSCC  SORTIEC          REL     MATRICE DE MASSE ELEMENTAIREC          KERRE   INDICE NON NUL SI PROBLEME , NUL SINONC-----------------------------------------------------------------------      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)      DIMENSION REL(LRE,*),CARAC(*),COORD(3,*),     &          XE(2),YE(2),ZE(2)      REAL*8 LCARC      KERRE=0C      RAY  = CARAC(1)      RC   = CARAC(2)      RHO  = CARAC(10)      C    = CARAC(11)      RHOR = CARAC(12)      CREF = CARAC(13)      LCAR = CARAC(14)C      DO 2 I=1,2      XE(I)=COORD(1,I)      YE(I)=COORD(2,I)  2   ZE(I)=COORD(3,I)      XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2      XL=SQRT(XL)      IF(XL.EQ.0.) GO TO 999CC MODIFICATION DE LA LONGUEUR EFFECTIVE DU TUYAU DANS LE CAS D'UN COUDEC      IF (RC.NE.0) THEN      XL = 2*RC*ASIN(XL/(2*RC))      ENDIFCC CALCUL DE LA MATRICE MASSE ELEMENTAIRE EN AXES LOCAUXC      CALL ACMALO(REL,XL,RAY,RHO,RHOR,C,CREF,LCAR,LRE)      RETURN999   KERRE=1      RETURN      END     

© Cast3M 2003 - Tous droits réservés.
Mentions légales