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    CB215821  26/08/24    21:15:12     12622          
  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.       V(I)=COORD(I,2)-COORD(I,1)
  28.     1 CONTINUE
  29.       XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  30.       XL=SQRT(XL)
  31.       IF(XL.EQ.0.) GO TO 999
  32. C
  33.       FAC=ES/(XL**3)
  34.       DO 1000 I=1,IDIM
  35.       II=I+IDIM
  36.       DO 2 J=1,IDIM
  37.       JJ=J+IDIM
  38.       FAC2=FAC*V(I)*V(J)
  39.       REL(I,J)=FAC2
  40.       REL(II,J)=-FAC2
  41.       REL(I,JJ)=-FAC2
  42.       REL(II,JJ)=FAC2
  43.   2   CONTINUE
  44.  1000 CONTINUE
  45.       RETURN
  46.  999  KERRE=1
  47.       RETURN
  48.       END
  49.  
  50.  
  51.  

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