Numérotation des lignes :

C TIFCOR    SOURCE    CHAT      05/01/13    03:38:45     5004      SUBROUTINE TIFCOR(REL,LRE,CARA,COORD,VROT,     >           WORK,LHOOK,DDHOOK,KERRE)C-----------------------------------------------------------------------C  ROUTINE DE CALCUL DE LA MATRICE DE COUPLAGE GYROSCOPIQUE DE L ELEMENT DEC  POUTRE DE TIMOSCHENKO EN BIFLEXION DANS LES AXES GLOBAUXCC  ENTREEC          LRE     DIMENSION DU TABLEAU RELC          VECT    DIRECTION DE L'OY LOCALC          OMEG    VITESSE DE ROTATIONC          COORD   COORDONNEES DES 2 NOEUDSC          DDHOOK(LHOOK,LHOOK) MATRICE DE "HOOK"CC  SORTIEC          REL     MATRICE DE RAIDEUR ELEMENTAIREC          KERRE   INDICE NON NUL SI PROBLEME , NUL SINONC  TRAVAILC          WORK(300)CC-----------------------------------------------------------------------C     Didier Combescure Mars 2003C-----------------------------------------------------------------------      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)      DIMENSION REL(LRE,*),COORD(3,*),CARA(*),VECT(3),P(3,3),     .          XE(2),YE(2),ZE(2),WORK(*),VROT(*),VROTL(3)C      VECT(1) = CARA(2)      VECT(2) = CARA(3)      VECT(3) = CARA(4)C      KERRE=0      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 999CCC  CALCUL DE LA MATRICE DE PASSAGEC      CALL POUPAS(XE,YE,ZE,VECT,P,KERRE)      IF(KERRE.NE.0) RETURNCC   CALCUL DU VECTEUR ROTATION DANS LE REPERE LOCALC      CALL VECROT(VROT,VROTL,P)CC  CALCUL DE LA MATRICE DE COUPLAGE GYROSCOPIQUE EN AXES LOCAUXC      CALL TIFCOL(WORK,XL,LHOOK,DDHOOK,VROTL)CC  PASSAGE DE LA MATRICE DE MASSE EN AXES GLOBAUXC      CALL POUROG(REL,LRE,P,WORK)CCC      RETURN999   KERRE=1      RETURN      END

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