1 : C TIMMA2 SOURCE CHAT 05/01/13 03:40:27 5004 2 : C TIMMAS SOURCE AM 93/01/12 21:36:19 781 3 : SUBROUTINE TIMMA2(REL,LRE,CARAC,COORD,WORK,KERRE) 4 : C----------------------------------------------------------------------- 5 : C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE DE L ELEMENT DE 6 : C POUTRE DE TIMOSCHENKO DANS LES AXES GLOBAUX 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 ( 10 VALEURS ) 12 : C COORD COORDONNEES DES 2 NOEUDS 13 : C 14 : C SORTIE 15 : C REL MATRICE DE MASSE ELEMENTAIRE 16 : C KERRE INDICE NON NUL SI PROBLEME , NUL SINON 17 : C TRAVAIL 18 : C WORK(300) 19 : C 20 : C D'APRES TIMMAS 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 : XINZ=CARAC(2) 29 : C 30 : C CORRECTION STCAN LE 6/7/89 : ON PREND L'INERTIE POLAIRE EN TORSION 31 : C 32 : SD=CARAC(1) 33 : RHO=CARAC(4) 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 MASSE EN AXES LOCAUX 42 : C 42 : CALL TIMAL2(WORK,XL,RHO,SD,XINZ) 43 : C 44 : C CALCUL DE LA MATRICE DE PASSAGE 45 : C 45 : CALL POUPA2(XE,YE,P,KERRE) 46 : IF(KERRE.NE.0) RETURN 47 : C 48 : C PASSAGE DE LA MATRICE DE MASSE EN AXES GLOBAUX 49 : C 49 : CALL POURO2(REL,LRE,P,WORK) 50 : RETURN 51 : 999 KERRE=1 52 : RETURN 53 : END 54 : 55 :
© Cast3M 2003 - All rights reserved.
Disclaimer