1 : C POUKS2 SOURCE CHAT 05/01/13 02:19:07 5004 2 : C POUKS2 SOURCE CHAT 89/05/26 21:32:11 3 : SUBROUTINE POUKS2(REL,LRE,CONT,CARAC,COORD,WORK,KERRE) 4 : C----------------------------------------------------------------------- 5 : C ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR GEOMETRIQUE DE L ELEMENT D 6 : C POUTRE DROITE DANS LES AXES GLOBAUX 7 : C 8 : C ENTREE 9 : C LRE DIMENSION DU TABLEAU REL 10 : C CONT CONTRAINTES 11 : C CARAC TABLEAU DE CARACTERISTIQUES GEOMETRIQUES 12 : C DE LA POUTRE ( 9 VALEURS ) 13 : C COORD COORDONNEES DES 2 NOEUDS 14 : C 15 : C SORTIE 16 : C REL MATRICE DE RIGIDITE GEOMETRIQUE ELEMENTAIRE 17 : C KERRE CODE D ERREUR NON NUL SI PROBLEME , NUL SINON 18 : C TRAVAIL 19 : C WORK(300) 20 : C 21 : C MAI 85 MILLARD 22 : C----------------------------------------------------------------------- 23 : IMPLICIT INTEGER(I-N) 24 : IMPLICIT REAL*8(A-H,O-Z) 25 : DIMENSION REL(LRE,*),CARAC(*),COORD(3,*),P(2,2), 26 : . XE(2),YE(2),WORK(*),CONT(*) 27 : C 28 : KERRE=0 29 : DO 2 I=1,2 30 : XE(I)=COORD(1,I) 31 : 2 YE(I)=COORD(2,I) 32 : XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2 33 : XL=SQRT(XL) 34 : IF(XL.EQ.0.) GO TO 999 35 : C 36 : C CALCUL DE LA RIGIDITE EN AXES LOCAUX 37 : C 38 : SIGM=0.5D00*(CONT(1)+CONT(4)) 38 : CALL POKSL2(WORK,XL,SIGM) 39 : C 40 : C CALCUL DE LA MATRICE DE PASSAGE 41 : C 41 : CALL POUPA2(XE,YE,P,KERRE) 42 : IF(KERRE.NE.0) RETURN 43 : C 44 : C PASSAGE DE LA MATRICE DE RIGIDITE EN AXES GLOBAUX 45 : C 45 : CALL POURO2(REL,LRE,P,WORK) 46 : RETURN 47 : 999 KERRE=1 48 : RETURN 49 : END 50 : 51 :
© Cast3M 2003 - All rights reserved.
Disclaimer