Télécharger lump1.eso

Retour à la liste

Numérotation des lignes :

  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