1 : C POUBSG SOURCE CHAT 05/01/13 02:18:27 5004 2 : SUBROUTINE POUBSG(FORCE,CARAC,COORD,CONT,WORK,KERRE) 3 : C----------------------------------------------------------------------- 4 : C ROUTINE DE CALCUL DES FORCES INTERNES B*SIGMA 5 : C POUR LA POUTRE DROITE 6 : C 7 : C ENTREE 8 : C CONT CONTRAINTES 9 : C CARAC TABLEAU DE CARACTERISTIQUES GEOMETRIQUES ET 10 : C MATERIELLES DE LA POUTRE ( 11 VALEURS ) 11 : C COORD COORDONNEES DES 2 NOEUDS 12 : C 13 : C SORTIE 14 : C FORCE FORCES INTERNES 15 : C KERRE INDICE D'ERREUR ( 0 SI PAS DE PB ) 16 : C TRAVAIL 17 : C WORK(300) 18 : C 19 : C MAI 85 MILLARD 20 : C----------------------------------------------------------------------- 21 : IMPLICIT INTEGER(I-N) 22 : IMPLICIT REAL*8(A-H,O-Z) 23 : DIMENSION CONT(*),CARAC(*),COORD(3,*),VECT(3),P(3,3), 24 : . XE(2),YE(2),ZE(2),FORCE(*),WORK(*) 25 : C 26 : KERRE=0 27 : DO 1 I=1,3 28 : 1 VECT(I)=CARAC(6+I) 29 : DO 2 I=1,2 30 : XE(I)=COORD(1,I) 31 : YE(I)=COORD(2,I) 32 : 2 ZE(I)=COORD(3,I) 33 : XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2 34 : XL=SQRT(XL) 35 : IF(XL.EQ.0.) GO TO 999 36 : C 37 : C VALEURS MOYENNES DANS L ELEMENT 38 : C 39 : SIGM=0.5D00*(CONT(1)+CONT(7)) 40 : TORM=0.5D00*(CONT(4)+CONT(10)) 41 : WORK(1)=-SIGM 42 : WORK(7)= SIGM 43 : WORK(4)=-TORM 44 : WORK(10)=TORM 45 : WORK(2)=-CONT(2) 46 : WORK(8)= CONT(8) 47 : WORK(3)=-CONT(3) 48 : WORK(9)= CONT(9) 49 : WORK(5)= -CONT(5) 50 : WORK(11)= CONT(11) 51 : WORK(6)= -CONT(6) 52 : WORK(12)= CONT(12) 53 : C 54 : C CALCUL DE LA MATRICE DE PASSAGE 55 : C 55 : CALL POUPAS(XE,YE,ZE,VECT,P,KERRE) 56 : IF(KERRE.NE.0) RETURN 57 : C 58 : C PASSAGE DES FORCES EN AXES GLOBAUX 59 : C 59 : CALL POUVEC(WORK,FORCE,P,2) 60 : RETURN 61 : 999 KERRE=1 62 : RETURN 63 : END 64 :
© Cast3M 2003 - All rights reserved.
Disclaimer