Télécharger rigbex.eso

Retour à la liste

Numérotation des lignes :

  1. C RIGBEX SOURCE CHAT 05/01/13 02:56:35 5004
  2. C RIGBEX SOURCE CHAT 89/05/26 22:52:13
  3. SUBROUTINE RIGBEX(REL,LRE,ES,AL,TGL1)
  4. C-----------------------------------------------------------------------
  5. C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE
  6. C DE L ELEMENT DE BARRE
  7. C-----------------------------------------------------------------------
  8. C ENTREE
  9. C LRE DIMENSION DU TABLEAU REL
  10. C ES YOUN * SECT
  11. C
  12. C SORTIE
  13. C REL MATRICE DE RIGIDITE ELEMENTAIRE
  14. C KERRE 0 SI PAS DE PB
  15. C 1 SI LONGUEUR NULLE
  16. C-----------------------------------------------------------------------
  17. IMPLICIT INTEGER(I-N)
  18. IMPLICIT REAL*8(A-H,O-Z)
  19. PARAMETER(XZER=0.D0)
  20. C
  21. DIMENSION REL(LRE,LRE),TGL1(2,LRE),C(2),RE(2,2),C1(2)
  22. C
  23. FAC2=ES/AL
  24. RE(1,1)=FAC2
  25. RE(1,2)=-FAC2
  26. RE(2,1)=-FAC2
  27. RE(2,2)=FAC2
  28. C
  29. DO 5 I =1,LRE
  30. C
  31. DO 6 I1=1,2
  32. C1(I1)=TGL1(I1,I)
  33. 6 CONTINUE
  34.  
  35. DO 1 J=1,2
  36. S=XZER
  37. DO 2 K=1,2
  38. S=S+C1(K)*RE(K,J)
  39. 2 CONTINUE
  40. C(J)=S
  41. 1 CONTINUE
  42. C
  43. DO 3 J=1,LRE
  44. S=XZER
  45. DO 4 K=1,2
  46. S=S+C(K)*TGL1(K,J)
  47. 4 CONTINUE
  48. REL(I,J)=S
  49. 3 CONTINUE
  50. C
  51. 5 CONTINUE
  52. * WRITE(*,*) ' matrice K elem ='
  53. * WRITE(*,1190) ((REL(IA,IC),IA=1,LRE),IC=1,LRE)
  54. *1190 FORMAT(12(12(1X,E8.2),/))
  55. RETURN
  56. END
  57.  
  58.  
  59.  

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