Esope source files

   1 : C COUMAS    SOURCE    PV        11/03/07    21:16:09     6885           
   2 :       SUBROUTINE COUMAS(REL,LRE,CARAC,COORD,KERRE)                         
   3 : C-----------------------------------------------------------------------        
   4 : C  ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE DE L ELEMENT DE         
   5 : C  COUPLAGE LITU                                         
   6 : C                                                                               
   7 : C  ENTREE                                                                       
   8 : C          LRE     DIMENSION DU TABLEAU REL                                     
   9 : C          CARAC   TABLEAU DE CARACTERISTIQUES GEOMETRIQUES ET                  
  10 : C                  MATERIELLES DU TUYAU FLUIDE  ( 8 VALEURS )                     
  11 : C          COORD   COORDONNEES DES 2 NOEUDS                                     
  12 : C                                                                               
  13 : C  SORTIE                                                                       
  14 : C          REL     MATRICE DE MASSE ELEMENTAIRE                                 
  15 : C          KERRE   INDICE NON NUL SI PROBLEME , NUL SINON                       
  16 : C-----------------------------------------------------------------------        
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  19 :       DIMENSION REL(LRE,LRE),CARAC(200),COORD(3,2),P(10,10),                 
  20 :      &          XE(2),YE(2),ZE(2),XMATLO(10,10),VECT(3)
  21 :       REAL*8 LCAR
  22 : C                                                                               
  23 :       KERRE=0                                                                   
  24 : C
  25 :       RAY  = CARAC(1)
  26 :       RC   = CARAC(2)
  27 :       RHO  = CARAC(10)
  28 :       LCAR = CARAC(11)
  29 :       RHOR = CARAC(12)
  30 : C
  31 :       DO 1 I=1,3
  32 :       VECT(I)=CARAC(2+I)
  33 :    1  CONTINUE
  34 : C
  35 :       DO 2 I=1,2                                                                
  36 :       XE(I)=COORD(1,I)                                                          
  37 :       YE(I)=COORD(2,I)                                                          
  38 :       ZE(I)=COORD(3,I)
  39 :    2  CONTINUE                                                        
  40 :       XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2                     
  41 :       XL=SQRT(XL)                                                               
  42 :       IF (XL.EQ.0.) GO TO 999
  43 : C
  44 : C MODIFICATION DE LA LONGUEUR EFFECTIVE DU TUYAU DANS LE CAS D'UN COUDE
  45 : C                              
  46 :       IF (RC.NE.0) THEN
  47 :       XL = 2*RC*ASIN(XL/(2*RC))
  48 :       ENDIF
  49 : C
  50 : C CALCUL DE LA MATRICE MASSE ELEMENTAIRE EN AXES LOCAUX
  51 : C
  51 :       CALL COUMAL(XMATLO,XL,RAY,RC,RHO,LCAR,RHOR)
  52 : C
  53 : C CALCUL DE LA MATRICE DE PASSAGE
  54 : C
  54 :       CALL COUPAS(XE,YE,ZE,VECT,P,KERRE)
  55 : C
  56 : C PASSAGE DE LA MATRICE DE MASSE EN AXES GLOBAUX
  57 : C
  57 :       CALL COUROT(REL,LRE,P,XMATLO)
  58 :       RETURN
  59 : 999   KERRE=1
  60 :       RETURN
  61 :       END
  62 :  
  63 :