1 : C POURH2 SOURCE AM 08/05/19 21:15:11 6108 2 : SUBROUTINE POURH2(REL,LRE,CARAC,COORD,WORK, 3 : & IMAT,LHOOK,DDHOOK,KERRE) 4 : C----------------------------------------------------------------------- 5 : C ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR ELEMENTAIRE DE L ELEMENT D 6 : C POUTRE DROITE DANS LES AXES GLOBAUX EN 2D 7 : C 8 : C ENTREE 9 : C LRE DIMENSION DU TABLEAU REL 10 : C CARAC TABLEAU DE CARACTERISTIQUES GEOMETRIQUES ET 11 : C MATERIELLES DE LA POUTRE ( 11 VALEURS ) 12 : C COORD COORDONNEES DES 2 NOEUDS 13 : C 14 : C SORTIE 15 : C REL MATRICE DE RIGIDITE ELEMENTAIRE 16 : C KERRE CODE D ERREUR NON NUL SI PROBLEME , NUL SINON 17 : C TRAVAIL 18 : C WORK(300) 19 : C 20 : C VIENT DE TEDEL AOUT 84 21 : C----------------------------------------------------------------------- 22 : IMPLICIT INTEGER(I-N) 23 : IMPLICIT REAL*8(A-H,O-Z) 24 : DIMENSION REL(LRE,*),CARAC(*),COORD(3,*),P(2,2), 25 : . XE(2),YE(2),WORK(*) 26 : DIMENSION DDHOOK(LHOOK,*) 27 : C 28 : KERRE=0 29 : XINZ=CARAC(2) 30 : SD=CARAC(1) 31 : SREDY=CARAC(3) 32 : YOUN=CARAC(4) 33 : CISAIL=CARAC(5) 34 : DO 2 I=1,2 35 : XE(I)=COORD(1,I) 36 : 2 YE(I)=COORD(2,I) 37 : XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2 38 : XL=SQRT(XL) 39 : IF(XL.EQ.0.) GO TO 999 40 : C 41 : C CALCUL DE LA RIGIDITE EN AXES LOCAUX 42 : C 42 : CALL PORIM2(WORK,XL,YOUN,CISAIL,SD,SREDY,XINZ, 43 : & IMAT,LHOOK,DDHOOK) 44 : C 45 : C CALCUL DE LA MATRICE DE PASSAGE 46 : C 46 : CALL POUPA2(XE,YE,P,KERRE) 47 : IF(KERRE.NE.0) RETURN 48 : C 49 : C PASSAGE DE LA MATRICE DE RIGIDITE EN AXES GLOBAUX 50 : C 50 : CALL POURO2(REL,LRE,P,WORK) 51 : RETURN 52 : 999 KERRE=1 53 : RETURN 54 : END 55 : 56 : 57 : 58 :
© Cast3M 2003 - All rights reserved.
Disclaimer