Télécharger la plaquette Cast3M


fr en

Télécharger lump6.eso

Retour à la liste

Numérotation des lignes : 

lump6

  1. C LUMP6     SOURCE    CHAT      05/01/13    01:27:17     5004
  2.       SUBROUTINE LUMP6(REL,LRE,RE)
  3. c====================================================================
  4. c
  5. c     utilisé par l'opérateur lump pour les éléments pout timo et tuya
  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     pour les dd de translation
  14. c
  15. c     pour les ddl de rotation on prend la trace divisée par le
  16. c     nombre de termes
  17. c
  18. c     ce travail est realisé dans le repere global
  19. c
  20. c     on suppsoe que re est initialisée à 0
  21. c
  22. c====================================================================
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25.       DIMENSION REL(LRE,*),RE(LRE,*)
  26. *
  27. *
  28. *
  29. c      on traite ux uy uz
  30. c       boucle sur les noeuds
  31.       NTRA = 3
  32.       IF (LRE .NE. 12 ) THEN
  33.         WRITE(*,*) 'Erreur dans lump6.eso'
  34.         CALL ERREUR(5)
  35.       ENDIF
  36.  
  37.       DO 430 I=1,2
  38. c         boucle sur les ddl ux uy uz locaux
  39.           DO 420 J=1,NTRA
  40.               IDERIJ = (NTRA+3)*(I-1) + J
  41.               SUM = 0.D0
  42.               DO 410 K=1,2
  43.                   DO 400 L=1,NTRA
  44.                       IDERKL= (NTRA+3)*(K-1)+L
  45.                       SUM = SUM + REL(IDERIJ,IDERKL)
  46.  400              CONTINUE
  47.  410          CONTINUE
  48.               RE(IDERIJ,IDERIJ) = SUM
  49.  420      CONTINUE
  50.  430  CONTINUE
  51. c
  52. c        on traite maintenant rx ry rz
  53. c
  54.       NTRA = 3
  55.       TRAC = 0.D0
  56. c       boucle sur les noeuds
  57.       DO 530 I=1,2
  58. c         boucle sur les ddl rx ry rz
  59.           DO 520 J=1,NTRA
  60.               IDERIJ = (NTRA+3)*(I-1) + 3 + J
  61.               TRAC = TRAC + REL(IDERIJ,IDERIJ)
  62.  520      CONTINUE
  63.  530  CONTINUE
  64. c
  65.       DO 630 I=1,2
  66. c         boucle sur les ddl rx ry rz  remplir re()
  67.           DO 620 J=1,NTRA
  68.               IDERIJ = (NTRA+3)*(I-1) + 3 + J
  69.               RE(IDERIJ,IDERIJ) = TRAC / 6.D0
  70.  620      CONTINUE
  71.  630  CONTINUE
  72. *
  73.       RETURN
  74.       END
  75.  
  76.  
  77.  

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