Télécharger la plaquette Cast3M


fr en

Télécharger barbsg.eso

Retour à la liste

Numérotation des lignes : 

barbsg

  1. C BARBSG    SOURCE    CHAT      05/01/12    21:33:33     5004
  2.       SUBROUTINE BARBSG(FORCE,COORD,EFFORT,KERRE)
  3. C-----------------------------------------------------------------------
  4. C  ROUTINE DE CALCUL DES FORCES INTERNES B*SIGMA
  5. C             POUR LA BARRE
  6. C
  7. C  ENTREE
  8. C          EFFORT  L'EFFORT DANS LA BARRE
  9. C          COORD   COORDONNEES DES 2 NOEUDS
  10. C
  11. C  SORTIE
  12. C          FORCE   FORCES INTERNES
  13. C          KERRE   INDICE D'ERREUR ( 0 SI PAS DE PB )
  14. C
  15. C-----------------------------------------------------------------------
  16.       IMPLICIT INTEGER(I-N)
  17.       IMPLICIT REAL*8(A-H,O-Z)
  18.  
  19. -INC PPARAM
  20. -INC CCOPTIO
  21.       DIMENSION COORD(3,*),V(3),FORCE(*)
  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.       DO 2 I=1,IDIM
  31.   2   V(I)=V(I)*EFFORT/XL
  32. C
  33. C  REMPLISSAGE DU VECTEUR FORCE
  34. C
  35.       DO 3 I=1,IDIM
  36.       FORCE(I)=-V(I)
  37.    3  FORCE(I+IDIM)=V(I)
  38.       RETURN
  39. 999   KERRE=1
  40.       RETURN
  41.       END
  42.  
  43.  

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