Télécharger lump6.eso

Retour à la liste

Numérotation des lignes :

  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