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