Esope source files

   1 : C LUMP1     SOURCE    CHAT      05/01/13    01:26:54     5004           
   2 :       SUBROUTINE LUMP1(REL,LRE,RE,IFOUR)
   3 : c====================================================================
   4 : c
   5 : c     utilisé par l'opérateur LUMP pour les éléments massifs
   6 : c
   7 : c     remplissage de RE à partir de 
   8 : c     rel (lre,lre) matrice de rigidite triangulaire
   9 : c     inferieure.
  10 : c
  11 : c     on diagonalise la matrice en somme les termes de la
  12 : c     meme lignes, somme que l'on met sur la diagonale
  13 : c
  14 : c     on suppsoe que RE est initialisée à 0    
  15 : c
  16 : c====================================================================
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)
  19 :       DIMENSION REL(LRE,*),RE(LRE,*)
  20 : *     
  21 : *
  22 :       IF ( IFOUR .EQ. -3) LRE=LRE-3 
  23 : *
  24 :       DO 200 IA=1,LRE
  25 :           SUM = 0.D0
  26 :           DO 100 IB=1,LRE
  27 :               I = MAX(IA,IB)
  28 :               J = MIN(IA,IB)
  29 :               SUM = SUM + REL(I,J)   
  30 :  100      CONTINUE
  31 :           RE(IA,IA) = SUM 
  32 :  200  CONTINUE
  33 : * 
  34 : *  cas des déformations planes généralisées 
  35 : *  la diagonalisation est différente pour les ddl transverses
  36 :       IF ( IFOUR .EQ. -3) THEN
  37 : *       inconnue RX
  38 :         RE(LRE+2,LRE+2)=(REL(LRE+2,LRE+2)+REL(LRE+3,LRE+3))/2.D0
  39 : *       inconue RY
  40 :         RE(LRE+3,LRE+3)=RE(LRE+2,LRE+2)
  41 : *       inconnue UZ
  42 :         RE(LRE+1,LRE+1)=REL(LRE+1,LRE+1)     
  43 :         LRE = LRE + 3
  44 :       ENDIF
  45 : *
  46 :       RETURN
  47 :       END
  48 :  
  49 :