1 : C PO1RIG SOURCE CHAT 05/01/13 02:15:45 5004 2 : SUBROUTINE PO1RIG(REL,LRE,ES,COORD,KERRE,XDPGE,YDPGE) 3 : C----------------------------------------------------------------------- 4 : C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE 5 : C DE L ELEMENT POI1 6 : C----------------------------------------------------------------------- 7 : C ENTREE 8 : C LRE DIMENSION DU TABLEAU REL 9 : C ES YOUN * SECT 10 : C COORD COORDONNEES DU NOEUD 11 : C XDPGE COORDONNEES DU POINT AUTOUR DUQUEL SE FAIT LA 12 : C YDPGE LA ROTATION DE LA SECTION EN DEFORMATIONS PLANES 13 : C GENERALISEES 14 : C SORTIE 15 : C REL MATRICE DE RIGIDITE ELEMENTAIRE 16 : C KERRE 0 SI PAS DE PB 17 : C----------------------------------------------------------------------- 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : -INC CCREEL 21 : DIMENSION REL(LRE,*),COORD(3,*) 22 : * 23 : XX=COORD(1,1) 24 : YY=COORD(2,1) 25 : XXX=XDPGE-XX 26 : YYY=YY-YDPGE 27 : * 28 : DO 1 I=1,LRE 29 : DO 1 J=1,LRE 30 : 1 REL(I,J)=0.D0 31 : * 32 : REL(3,3)=ES 33 : REL(3,4)=ES*XXX 34 : REL(3,5)=ES*YYY 35 : REL(4,3)=REL(3,4) 36 : REL(4,4)=ES*XXX*XXX 37 : REL(4,5)=ES*XXX*YYY 38 : REL(5,3)=REL(3,5) 39 : REL(5,4)=REL(4,5) 40 : REL(5,5)=ES*YYY*YYY 41 : RETURN 42 : END 43 : 44 : 45 :
© Cast3M 2003 - All rights reserved.
Disclaimer