Sources Esope | Cast3M

tifal2
C TIFAL2    SOURCE    CHAT      05/01/13    03:38:29     5004
C TIFAL2    SOURCE    AM        94/02/25    22:55:09     1061
      SUBROUTINE TIFAL2(RP,XL,LHOOK,DDHOOK)
C-----------------------------------------------------------------------
C
C    ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE
C    POUTRE DE TIMOSHENKO DROITE 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 MASSE ELEMENTAIRE 12X12
C
C D'APRES TIFALO DC 98
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION RP(6,6)
      DIMENSION DDHOOK(LHOOK,LHOOK)
C
C  INITIALISATIONS
C
      DO 120 I=1,6
      DO 120 J=1,6
  120 RP(I,J)=0.D00
C
C  REMPLISSAGE DE LA MATRICE
C  FX1
      RP( 1, 1)=DDHOOK( 1, 1)*XL/3.0D0
      RP( 1, 3)=DDHOOK( 1, 3)*XL/3.0D0
      RP( 1, 4)=RP( 1, 1)/2.0D0
      RP( 1, 6)=RP( 1, 3)/2.0D0
C  FY1
      RP( 2, 2)=DDHOOK( 2, 2)*XL/3.0D0
      RP( 2, 5)=RP( 2, 2)/2.0D0
C  MZ1
      RP( 3, 3)=DDHOOK( 3, 3)*XL/3.0D0
      RP( 3, 4)=DDHOOK( 3, 1)*XL/6.0D0
      RP( 3, 6)=RP( 3, 3)/2.0D0
C  FX2
      RP( 4, 4)=RP( 1, 1)
      RP( 4, 6)=RP( 1, 3)
C  FY2
      RP( 5, 5)=RP( 2, 2)
C  MZ2
      RP( 6, 6)=RP( 3, 3)
C
      DO 9 I=1,5
      DO 9 J=I+1,6
  9   RP(J,I)= RP(I,J)
      RETURN
      END