Télécharger la plaquette Cast3M


fr en

Télécharger barrig.eso

Retour à la liste

Numérotation des lignes : 

barrig

  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