Esope source files

   1 : C BARMAS    SOURCE    CHAT      05/01/12    21:34:41     5004           
   2 :       SUBROUTINE BARMAS(REL,LRE,ROS,COORD)                                      
   3 : C-----------------------------------------------------------------------        
   4 : C  ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE DE L ELEMENT DE         
   5 : C  BARRE                                                                        
   6 : C                                                                               
   7 : C  ENTREE                                                                       
   8 : C          LRE     DIMENSION DU TABLEAU REL                                     
   9 : C          ROS     RHO * SECT                                                   
  10 : C          COORD   COORDONNEES DES 2 NOEUDS                                     
  11 : C                                                                               
  12 : C  SORTIE                                                                       
  13 : C          REL     MATRICE DE MASSE ELEMENTAIRE                                 
  14 : C-----------------------------------------------------------------------        
  15 :       IMPLICIT INTEGER(I-N)
  16 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  17 : -INC CCOPTIO                                                                    
  18 :       DIMENSION REL(LRE,*),COORD(3,*),V(3)                                      
  19 : C                                                                               
  20 :       V(3)=0.                                                                   
  21 :       DO 1 I=1,IDIM                                                             
  22 :   1   V(I)=COORD(I,2)-COORD(I,1)                                                
  23 :       XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))                                        
  24 :       XL=SQRT(XL)                                                               
  25 : C                                                                               
  26 :       FAC=ROS*XL/6.                                                             
  27 :       FAC2=FAC*2.                                                               
  28 :       DO 2 I=1,IDIM                                                             
  29 :       II=I+IDIM                                                                 
  30 :       REL(I,I)=FAC2                                                             
  31 :       REL(II,I)=FAC                                                             
  32 :       REL(I,II)=FAC                                                             
  33 :   2   REL(II,II)=FAC2                                                           
  34 :       RETURN                                                                    
  35 :       END                                                                       
  36 :