Télécharger la plaquette Cast3M


fr en

Télécharger lump1.eso

Retour à la liste

Numérotation des lignes : 

lump1

  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.  
  50.  

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