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 : -INC CCOPTIO 20 : DIMENSION REL(LRE,*),COORD(3,*),V(3) 21 : C 22 : KERRE=0 23 : V(3)=0. 24 : DO 1 I=1,IDIM 25 : 1 V(I)=COORD(I,2)-COORD(I,1) 26 : XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3)) 27 : XL=SQRT(XL) 28 : IF(XL.EQ.0.) GO TO 999 29 : C 30 : FAC=ES/(XL**3) 31 : DO 2 I=1,IDIM 32 : II=I+IDIM 33 : DO 2 J=1,IDIM 34 : JJ=J+IDIM 35 : FAC2=FAC*V(I)*V(J) 36 : REL(I,J)=FAC2 37 : REL(II,J)=-FAC2 38 : REL(I,JJ)=-FAC2 39 : REL(II,JJ)=FAC2 40 : 2 CONTINUE 41 : RETURN 42 : 999 KERRE=1 43 : RETURN 44 : END 45 :
© Cast3M 2003 - All rights reserved.
Disclaimer