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