1 : C POUSTR SOURCE CHAT 05/01/13 02:20:42 5004 2 : SUBROUTINE POUSTR(COORD,DEPLA,CARAC,CONT,WORK,IDEUX) 3 : C----------------------------------------------------------------------- 4 : C ROUTINE DE CALCUL DES CONTRAINTES POUR LA POUTRE DROITE 5 : C 6 : C ENTREE 7 : C DEPLA DEPLACEMENTS 8 : C CARAC TABLEAU DE CARACTERISTIQUES GEOMETRIQUES ET 9 : C MATERIELLES DE LA POUTRE ( 11 VALEURS ) 10 : C COORD COORDONNEES DES 2 NOEUDS 11 : C IDEUX INDICE POUR DEFORMATIONS DU SECOND ORDRE 12 : C 13 : C SORTIE 14 : C CONT CONTRAINTES AUX 2 NOEUDS 6 PAR NOEUD 15 : C TRAVAIL 16 : C WORK(300) 17 : C VIENT DE TEDEL AOUT 84 18 : C----------------------------------------------------------------------- 19 : IMPLICIT INTEGER(I-N) 20 : IMPLICIT REAL*8(A-H,O-Z) 21 : DIMENSION CONT(*),CARAC(*),COORD(3,*),VECT(3),P(3,3), 22 : . XE(2),YE(2),ZE(2),DEPLA(*) 23 : C 24 : XINX=CARAC(1) 25 : XINY=CARAC(2) 26 : XINZ=CARAC(3) 27 : SD=CARAC(4) 28 : SREDY=CARAC(5) 29 : SREDZ=CARAC(6) 30 : DO 1 I=1,3 31 : 1 VECT(I)=CARAC(6+I) 32 : YOUN=CARAC(10) 33 : CISAIL=CARAC(11) 34 : DO 2 I=1,2 35 : XE(I)=COORD(1,I) 36 : YE(I)=COORD(2,I) 37 : 2 ZE(I)=COORD(3,I) 38 : XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2 39 : XL=SQRT(XL) 40 : C 41 : C CALCUL DE LA MATRICE DE PASSAGE 42 : C 42 : CALL POUPAS(XE,YE,ZE,VECT,P,KERRE) 43 : C 44 : C PASSAGE DES DEPLACEMENTS EN AXES LOCAUX 45 : C 45 : CALL POUVEC(DEPLA,WORK,P,1) 46 : C 47 : C CALCUL DES CONTRAINTES LOCALES 48 : C 48 : CALL POCOLO(WORK,CONT,XL,YOUN,CISAIL,SD,SREDY,SREDZ, 49 : . XINX,XINY,XINZ,IDEUX) 50 : RETURN 51 : END 52 :
© Cast3M 2003 - All rights reserved.
Disclaimer