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.  
  20. -INC PPARAM
  21. -INC CCOPTIO
  22. DIMENSION REL(LRE,*),COORD(3,*),V(3)
  23. C
  24. KERRE=0
  25. V(3)=0.
  26. DO 1 I=1,IDIM
  27. 1 V(I)=COORD(I,2)-COORD(I,1)
  28. XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  29. XL=SQRT(XL)
  30. IF(XL.EQ.0.) GO TO 999
  31. C
  32. FAC=ES/(XL**3)
  33. DO 2 I=1,IDIM
  34. II=I+IDIM
  35. DO 2 J=1,IDIM
  36. JJ=J+IDIM
  37. FAC2=FAC*V(I)*V(J)
  38. REL(I,J)=FAC2
  39. REL(II,J)=-FAC2
  40. REL(I,JJ)=-FAC2
  41. REL(II,JJ)=FAC2
  42. 2 CONTINUE
  43. RETURN
  44. 999 KERRE=1
  45. RETURN
  46. END
  47.  
  48.  

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