Sources Esope | Cast3M

tifcol
C TIFCOL    SOURCE    KK2000    05/05/25    21:15:29     5088
      SUBROUTINE TIFCOL(RP,XL,LHOOK,DDHOOK,VROTL)
C-----------------------------------------------------------------------
C
C    ROUTINE DE CALCUL DE LA MATRICE DE COUPLAGE GYROSCOPIQUE ELEMENTAIRE
C    D UNE POUTRE DE TIMOSHENKO DROITE MULTI-FIBRE EN AXES LOCAUX
C    (REPERE TOURNANT)
C
C  ENTREE
C    XL      LONGUEUR
C    DDHOOK(LHOOK,LHOOK) MATRICE DE "HOOK"
C
C  SORTIE
C    RP      MATRICE DE RAIDEUR ELEMENTAIRE 12X12
C
C-----------------------------------------------------------------------
C     Didier Combescure Mars 2003
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
-INC CCREEL
      DIMENSION RP(12,12),VROTL(3),ROME(6,6)
      DIMENSION DDHOOK(LHOOK,LHOOK)
C
C  INITIALISATIONS
C
      DO 120 I=1,12
      DO 120 J=1,12
  120 RP(I,J)=0.D00
C
      DO 121 I=1,6
      DO 121 J=1,6
  121 ROME(I,J)=0.D00
C
      ROME(1,2) = (-1.)*VROTL(3)
      ROME(1,3) = VROTL(2)
      ROME(2,3) = (-1.)*VROTL(1)
      ROME(2,1) = VROTL(3)
      ROME(3,1) = (-1.)*VROTL(2)
      ROME(3,2) = VROTL(1)
C      ROME(4,5) = (-1.)*VROTL(3)
C      ROME(4,6) = VROTL(2)
C      ROME(5,6) = (-1.)*VROTL(1)
C      ROME(5,4) = VROTL(3)
C      ROME(6,4) = (-1.)*VROTL(2)
C      ROME(6,5) = VROTL(1)
C
C
C  REMPLISSAGE DE LA MATRICE
C
      DO 201 I=1,3
        DO 201 J=1,3
        RP(I,J) =  DDHOOK( J, J)*ROME(I,J)*XL/3.
        RP(I+6,J) =DDHOOK( J, J)*ROME(I,J)*XL/6.
        RP(I,J+6) =DDHOOK( J, J)*ROME(I,J)*XL/6.
        RP(I+6,J+6) =DDHOOK( J, J)*ROME(I,J)*XL/3.
  201 CONTINUE
C
      RETURN
      END