Esope source files

   1 : C TIFAL2    SOURCE    CHAT      05/01/13    03:38:29     5004           
   2 : C TIFAL2    SOURCE    AM        94/02/25    22:55:09     1061           
   3 :       SUBROUTINE TIFAL2(RP,XL,LHOOK,DDHOOK)
   4 : C-----------------------------------------------------------------------
   5 : C
   6 : C    ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE
   7 : C    POUTRE DE TIMOSHENKO DROITE MULTI-FIBRE EN AXES LOCAUX
   8 : C    D APRES TEDEL
   9 : C
  10 : C  ENTREE
  11 : C    XL      LONGUEUR
  12 : C    DDHOOK(LHOOK,LHOOK) MATRICE DE "HOOK"
  13 : C
  14 : C  SORTIE
  15 : C    RP      MATRICE DE MASSE ELEMENTAIRE 12X12
  16 : C
  17 : C D'APRES TIFALO DC 98
  18 : C-----------------------------------------------------------------------
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)
  21 :       DIMENSION RP(6,6)
  22 :       DIMENSION DDHOOK(LHOOK,LHOOK)
  23 : C
  24 : C  INITIALISATIONS
  25 : C
  26 :       DO 120 I=1,6
  27 :       DO 120 J=1,6
  28 :   120 RP(I,J)=0.D00
  29 : C
  30 : C  REMPLISSAGE DE LA MATRICE
  31 : C  FX1
  32 :       RP( 1, 1)=DDHOOK( 1, 1)*XL/3.0D0
  33 :       RP( 1, 3)=DDHOOK( 1, 3)*XL/3.0D0
  34 :       RP( 1, 4)=RP( 1, 1)/2.0D0
  35 :       RP( 1, 6)=RP( 1, 3)/2.0D0
  36 : C  FY1
  37 :       RP( 2, 2)=DDHOOK( 2, 2)*XL/3.0D0
  38 :       RP( 2, 5)=RP( 2, 2)/2.0D0
  39 : C  MZ1
  40 :       RP( 3, 3)=DDHOOK( 3, 3)*XL/3.0D0
  41 :       RP( 3, 4)=DDHOOK( 3, 1)*XL/6.0D0
  42 :       RP( 3, 6)=RP( 3, 3)/2.0D0
  43 : C  FX2
  44 :       RP( 4, 4)=RP( 1, 1)
  45 :       RP( 4, 6)=RP( 1, 3)
  46 : C  FY2
  47 :       RP( 5, 5)=RP( 2, 2)
  48 : C  MZ2
  49 :       RP( 6, 6)=RP( 3, 3)
  50 : C
  51 :       DO 9 I=1,5
  52 :       DO 9 J=I+1,6
  53 :   9   RP(J,I)= RP(I,J)
  54 :       RETURN
  55 :       END
  56 :  
  57 :