Télécharger la plaquette Cast3M


fr en

Télécharger barksg.eso

Retour à la liste

Numérotation des lignes : 

barksg

  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