Esope source files

   1 : C BAMAEX    SOURCE    CHAT      05/01/12    21:31:50     5004           
   2 : C BAMAEX    SOURCE    DJER      93/03/02    22:53:00
   3 :       SUBROUTINE BAMAEX(REL,LRE,ROS,AL,TGL)
   4 : C-----------------------------------------------------------------------        
   5 : C  ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE DE L ELEMENT DE         
   6 : C  BARRE                                                                        
   7 : C                                                                               
   8 : C  ENTREE                                                                       
   9 : C          LRE     DIMENSION DU TABLEAU REL                                     
  10 : C          ROS     RHO * SECT                                                   
  11 : C          COORD   COORDONNEES DES 2 NOEUDS                                     
  12 : C                                                                               
  13 : C  SORTIE                                                                       
  14 : C          REL     MATRICE DE MASSE ELEMENTAIRE                                 
  15 : C-----------------------------------------------------------------------        
  16 :       IMPLICIT INTEGER(I-N)
  17 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  18 :       DIMENSION REL(LRE,LRE),TGL(2,LRE),RE(2,2),C(2),C1(2)
  19 : C                                                                               
  20 :       FAC=ROS*AL/6.D0
  21 :       FAC2=FAC*2.D0
  22 :       RE(1,1)=FAC2
  23 :       RE(1,2)=FAC
  24 :       RE(2,1)=FAC
  25 :       RE(2,2)=FAC2
  26 : C     
  27 :       DO 5 I=1,LRE
  28 :       DO 6 J1=1,2
  29 :         C1(J1)=TGL(J1,I)
  30 :     6 CONTINUE
  31 :       DO 1 J=1,2
  32 :       S=0.D0
  33 :       DO 2 K=1,2
  34 :     2 S=S+C1(K)*RE(K,J)
  35 :     1 C(J)=S
  36 :       DO 3 J=1,LRE
  37 :       S=0.D0
  38 :       DO 4 K=1,2
  39 :     4 S=S+C(K)*TGL(K,J)
  40 :     3 REL(I,J)=S
  41 :     5 CONTINUE
  42 : C
  43 :       RETURN
  44 :       END                                                                       
  45 :  
  46 :