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 : -INC CCOPTIO 19 : DIMENSION COORD(3,*),V(3),FORCE(*) 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 : DO 2 I=1,IDIM 29 : 2 V(I)=V(I)*EFFORT/XL 30 : C 31 : C REMPLISSAGE DU VECTEUR FORCE 32 : C 33 : DO 3 I=1,IDIM 34 : FORCE(I)=-V(I) 35 : 3 FORCE(I+IDIM)=V(I) 36 : RETURN 37 : 999 KERRE=1 38 : RETURN 39 : END 40 :
© Cast3M 2003 - All rights reserved.
Disclaimer