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

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