1 : C TIMRI2 SOURCE CHAT 05/01/13 03:40:46 5004 2 : C TIMRIG SOURCE AM 93/01/12 21:36:38 781 3 : SUBROUTINE TIMRI2(REL,LRE,CARAC,COORD,WORK,KERRE) 4 : C----------------------------------------------------------------------- 5 : C ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR ELEMENTAIRE DE L ELEMENT D 6 : C POUTRE DE TIMOSCHENKO 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 ( 5 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 D'APRES TIMRIG DC 98 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 : C 27 : KERRE=0 28 : SD=CARAC(1) 29 : XINZ=CARAC(2) 30 : SREDY=CARAC(3) 31 : YOUN=CARAC(4) 32 : CISAIL=CARAC(5) 33 : DO 2 I=1,2 34 : XE(I)=COORD(1,I) 35 : 2 YE(I)=COORD(2,I) 36 : XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2 37 : XL=SQRT(XL) 38 : IF(XL.EQ.0.) GO TO 999 39 : C 40 : C CALCUL DE LA RIGIDITE EN AXES LOCAUX 41 : C 41 : CALL TIRIL2(WORK,XL,YOUN,CISAIL,SD,SREDY,XINZ) 42 : C 43 : C CALCUL DE LA MATRICE DE PASSAGE 44 : C 44 : CALL POUPA2(XE,YE,P,KERRE) 45 : IF(KERRE.NE.0) RETURN 46 : C 47 : C PASSAGE DE LA MATRICE DE RIGIDITE EN AXES GLOBAUX 48 : C 48 : CALL POURO2(REL,LRE,P,WORK) 49 : RETURN 50 : 999 KERRE=1 51 : RETURN 52 : END 53 : 54 :
© Cast3M 2003 - All rights reserved.
Disclaimer