C PO1RIG SOURCE CHAT 05/01/13 02:15:45 5004 SUBROUTINE PO1RIG(REL,LRE,ES,COORD,KERRE,XDPGE,YDPGE) C----------------------------------------------------------------------- C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE C DE L ELEMENT POI1 C----------------------------------------------------------------------- C ENTREE C LRE DIMENSION DU TABLEAU REL C ES YOUN * SECT C COORD COORDONNEES DU NOEUD C XDPGE COORDONNEES DU POINT AUTOUR DUQUEL SE FAIT LA C YDPGE LA ROTATION DE LA SECTION EN DEFORMATIONS PLANES C GENERALISEES C SORTIE C REL MATRICE DE RIGIDITE ELEMENTAIRE C KERRE 0 SI PAS DE PB C----------------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) -INC CCREEL DIMENSION REL(LRE,*),COORD(3,*) * XX=COORD(1,1) YY=COORD(2,1) XXX=XDPGE-XX YYY=YY-YDPGE * DO 1 I=1,LRE DO 1 J=1,LRE 1 REL(I,J)=0.D0 * REL(3,3)=ES REL(3,4)=ES*XXX REL(3,5)=ES*YYY REL(4,3)=REL(3,4) REL(4,4)=ES*XXX*XXX REL(4,5)=ES*XXX*YYY REL(5,3)=REL(3,5) REL(5,4)=REL(4,5) REL(5,5)=ES*YYY*YYY RETURN END