Télécharger barksg.eso

Retour à la liste

Numérotation des lignes :

  1. C BARKSG SOURCE CHAT 05/01/12 21:34:27 5004
  2. SUBROUTINE BARKSG(REL,LRE,EFFORT,COORD,KERRE)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR GEOMETRIQUE DE L ELEMENT D
  5. C POUTRE DROITE DANS LES AXES GLOBAUX
  6. C
  7. C ENTREE
  8. C LRE DIMENSION DU TABLEAU REL
  9. C EFFORT EFFORT
  10. C COORD COORDONNEES DES 2 NOEUDS
  11. C
  12. C SORTIE
  13. C REL MATRICE DE RIGIDITE GEOMETRIQUE ELEMENTAIRE
  14. C KERRE CODE D ERREUR NON NUL SI PROBLEME , NUL SINON
  15. C-----------------------------------------------------------------------
  16. IMPLICIT INTEGER(I-N)
  17. IMPLICIT REAL*8(A-H,O-Z)
  18. -INC CCOPTIO
  19. DIMENSION REL(LRE,*),COORD(3,*),V(3)
  20. C
  21. KERRE=0
  22. V(3)=0.
  23. DO 1 I=1,IDIM
  24. 1 V(I)=COORD(I,2)-COORD(I,1)
  25. XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  26. XL=SQRT(XL)
  27. IF(XL.EQ.0.) GO TO 999
  28. C
  29. FAC=EFFORT/XL
  30. DO 2 I=1,IDIM
  31. II=I+IDIM
  32. REL(I,I)=FAC
  33. REL(II,I)=-FAC
  34. REL(I,II)=-FAC
  35. 2 REL(II,II)=FAC
  36. RETURN
  37. 999 KERRE=1
  38. RETURN
  39. END
  40.  
  41.  

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