Numérotation des lignes :

C POURHG    SOURCE    AM        08/05/19    21:15:13     6108      SUBROUTINE POURHG(REL,LRE,CARAC,COORD,WORK,IMAT,     &                    LHOOK, DDHOOK,KERRE)C-----------------------------------------------------------------------C  ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR ELEMENTAIRE DE L ELEMENT DC  POUTRE DROITE DANS LES AXES GLOBAUXCC  ENTREEC          LRE     DIMENSION DU TABLEAU RELC          CARAC   TABLEAU DE CARACTERISTIQUES GEOMETRIQUES ETC                  MATERIELLES DE LA POUTRE  ( 11 VALEURS )C          COORD   COORDONNEES DES 2 NOEUDSCC  SORTIEC          REL     MATRICE DE RIGIDITE ELEMENTAIREC          KERRE   CODE D ERREUR  NON NUL SI PROBLEME , NUL SINONC  TRAVAILC          WORK(300)CC     VIENT DE TEDEL AOUT 84C-----------------------------------------------------------------------      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)      DIMENSION REL(LRE,*),CARAC(*),COORD(3,*),VECT(3),P(3,3),     .          XE(2),YE(2),ZE(2),WORK(*)      DIMENSION DDHOOK(LHOOK,*)C      KERRE=0      XINX=CARAC(1)      XINY=CARAC(2)      XINZ=CARAC(3)      SD=CARAC(4)      SREDY=CARAC(5)      SREDZ=CARAC(6)      DO 1 I=1,31     VECT(I)=CARAC(6+I)      YOUN=CARAC(10)      CISAIL=CARAC(11)      DO 2 I=1,2      XE(I)=COORD(1,I)      YE(I)=COORD(2,I)  2   ZE(I)=COORD(3,I)      XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2      XL=SQRT(XL)      IF(XL.EQ.0.) GO TO 999CC  CALCUL DE LA RIGIDITE EN AXES LOCAUXC       CALL PORIMO(WORK,XL,YOUN,CISAIL,SD,SREDY,SREDZ,XINX,XINY,XINZ,     &            IMAT,LHOOK,DDHOOK) CC  CALCUL DE LA MATRICE DE PASSAGEC      CALL POUPAS(XE,YE,ZE,VECT,P,KERRE)      IF(KERRE.NE.0) RETURNCC  PASSAGE DE LA MATRICE DE RIGIDITE EN AXES GLOBAUXC      CALL POUROT(REL,LRE,P,WORK)      RETURN999   KERRE=1      RETURN      END

© Cast3M 2003 - Tous droits réservés.
Mentions légales