C TIFDL3    SOURCE    CB215821  26/08/24    21:18:34     12622          
       SUBROUTINE TIFDL3(RP,XL,LHOOK,DDHOOK,OMEG)
C-----------------------------------------------------------------------
C
C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
C    POUTRE DROITE DE TIMOSHENKO MULTI-FIBRE EN AXES LOCAUX
C    "D APRES TEDEL"
C
C  ENTREE
C    XL      LONGUEUR
C    DDHOOK(LHOOK,LHOOK) MATRICE DE HOOK
C
C  SORTIE
C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 12X12
C
C-----------------------------------------------------------------------
C     Pierre Pegon (ISPRA) Juillet/Aout 1993
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION RP(12,12)
      DIMENSION DDHOOK(LHOOK,LHOOK)
C
C  INITIALISATIONS
C
      DO 121 I=1,12
      DO 120 J=1,12
      RP(I,J)=0.D00
  120 CONTINUE
  121 CONTINUE
C
      RL=1.D00/XL
C
      GA = OMEG*0.5D0*(DDHOOK( 2, 2) + DDHOOK( 3, 3))
      EI = OMEG*0.5D0*(DDHOOK( 5, 5) + DDHOOK( 6, 6))
C
C  REMPLISSAGE DE LA MATRICE
C  FX1
C      RP( 1, 1)=DDHOOK( 1, 1)*RL
C      RP( 1, 7)=-RP( 1, 1)
C      RP( 7, 1)=-RP( 1, 1)
C  FY1
      RP( 2, 3)=GA*RL
      RP( 3, 2)=-GA*RL
      RP( 2, 5)=-GA/2.0D0
      RP( 5, 2)=GA/2.0D0
      RP( 2, 9)=-GA*RL
      RP( 9, 2)=GA*RL
      RP( 2,11)=-GA/2.0D0
      RP(11, 2)=GA/2.0D0
C  FZ1
      RP( 3, 6)=-GA/2.0D0
      RP( 3, 8)=GA*RL
      RP( 3, 12)=-GA/2.0D0
      RP( 6, 3)=GA/2.0D0
      RP( 8, 3)=-GA*RL
      RP( 12, 3)=GA/2.0D0
C  MX1
C      RP( 4, 4)=DDHOOK( 4, 4)*RL
C      RP( 4,10)=-RP( 4, 4)
C      RP( 10,4)=-RP( 4, 4)
C  MY1
      RP( 5, 6)=EI*RL+GA/(4.0D0*RL)
      RP( 5, 8)=-GA/2.0D0
      RP( 5,12)=-EI*RL+GA/(4.0D0*RL)
      RP( 6, 5)=-EI*RL-GA/(4.0D0*RL)
      RP( 8, 5)=GA/2.0D0
      RP(12, 5)=EI*RL-GA/(4.0D0*RL)
C  MZ1
      RP( 6, 9)=-GA/2.0D0
      RP( 6,11)=EI*RL-GA/(4.0D0*RL)
      RP( 9, 6)=GA/2.0D0
      RP( 11, 6)=-EI*RL+GA/(4.0D0*RL)
C  FX2
C      RP( 7, 7)=RP( 1, 1)
C  FY2
      RP( 8, 9)=GA*RL
      RP( 9, 8)=-GA*RL
      RP( 8, 11)=GA/2.0D0
      RP( 11, 8)=-GA/2.0D0
C  FZ2
      RP( 9, 12)=GA/2.0D0
      RP( 12, 9)=-GA/2.0D0
C  MX2
C      RP(10,10)=RP( 4, 4)
C  MY2
      RP(11,12)=EI*RL+GA/(4.0D0*RL)
C  MZ2
      RP(12,11)=-EI*RL-GA/(4.0D0*RL)
C
      RETURN
      END