Télécharger barrig.eso

Retour à la liste

Numérotation des lignes :

  1. C BARRIG SOURCE CHAT 05/01/12 21:35:03 5004
  2. SUBROUTINE BARRIG(REL,LRE,ES,COORD,KERRE)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE
  5. C DE L ELEMENT DE BARRE
  6. C-----------------------------------------------------------------------
  7. C ENTREE
  8. C LRE DIMENSION DU TABLEAU REL
  9. C ES YOUN * SECT
  10. C COORD COORDONNEES DES 2 NOEUDS
  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. -INC CCOPTIO
  20. DIMENSION REL(LRE,*),COORD(3,*),V(3)
  21. C
  22. KERRE=0
  23. V(3)=0.
  24. DO 1 I=1,IDIM
  25. 1 V(I)=COORD(I,2)-COORD(I,1)
  26. XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  27. XL=SQRT(XL)
  28. IF(XL.EQ.0.) GO TO 999
  29. C
  30. FAC=ES/(XL**3)
  31. DO 2 I=1,IDIM
  32. II=I+IDIM
  33. DO 2 J=1,IDIM
  34. JJ=J+IDIM
  35. FAC2=FAC*V(I)*V(J)
  36. REL(I,J)=FAC2
  37. REL(II,J)=-FAC2
  38. REL(I,JJ)=-FAC2
  39. REL(II,JJ)=FAC2
  40. 2 CONTINUE
  41. RETURN
  42. 999 KERRE=1
  43. RETURN
  44. END
  45.  
  46.  

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