Esope source files

   1 : C TIFILO    SOURCE    CHAT      05/01/13    03:39:02     5004           
   2 :       SUBROUTINE TIFILO(RP,XL,LHOOK,DDHOOK)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
   6 : C    POUTRE DROITE DE TIMOSHENKO MULTI-FIBRE EN AXES LOCAUX
   7 : C    "D APRES TEDEL"
   8 : C
   9 : C  ENTREE
  10 : C    XL      LONGUEUR
  11 : C    DDHOOK(LHOOK,LHOOK) MATRICE DE HOOK
  12 : C
  13 : C  SORTIE
  14 : C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 12X12
  15 : C
  16 : C-----------------------------------------------------------------------
  17 : C     Pierre Pegon (ISPRA) Juillet/Aout 1993
  18 : C-----------------------------------------------------------------------
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)
  21 :       DIMENSION RP(12,12)
  22 :       DIMENSION DDHOOK(LHOOK,LHOOK)
  23 : C
  24 : C  INITIALISATIONS
  25 : C
  26 :       DO 120 I=1,12
  27 :       DO 120 J=1,12
  28 :   120 RP(I,J)=0.D00
  29 : C
  30 :       RL=1.D00/XL
  31 : C
  32 : C  REMPLISSAGE DE LA MATRICE
  33 : C  FX1
  34 :       RP( 1, 1)=DDHOOK( 1, 1)*RL
  35 :       RP( 1, 5)=DDHOOK( 1, 5)*RL
  36 :       RP( 1, 6)=DDHOOK( 1, 6)*RL
  37 :       RP( 1, 7)=-RP( 1, 1)
  38 :       RP( 1,11)=-RP( 1, 5)
  39 :       RP( 1,12)=-RP( 1, 6)
  40 : C  FY1
  41 :       RP( 2, 2)=DDHOOK( 2, 2)*RL
  42 :       RP( 2, 4)=DDHOOK( 2, 4)*RL
  43 :       RP( 2, 6)=DDHOOK( 2, 2)/2.0D0
  44 :       RP( 2, 8)=-RP( 2, 2)
  45 :       RP( 2,10)=-RP( 2, 4)
  46 :       RP( 2,12)= RP( 2, 6)
  47 : C  FZ1
  48 :       RP( 3, 3)=DDHOOK( 3, 3)*RL
  49 :       RP( 3, 4)=DDHOOK( 3, 4)*RL
  50 :       RP( 3, 5)=-DDHOOK( 3, 3)/2.0D0
  51 :       RP( 3, 9)=-RP( 3, 3)
  52 :       RP( 3,10)=-RP( 3, 4)
  53 :       RP( 3,11)= RP( 3, 5)
  54 : C  MX1
  55 :       RP( 4, 4)=DDHOOK( 4, 4)*RL
  56 :       RP( 4, 5)=-DDHOOK( 4, 3)/2.0D0
  57 :       RP( 4, 6)=DDHOOK( 4, 2)/2.0D0
  58 :       RP( 4, 8)=-DDHOOK( 4, 2)*RL
  59 :       RP( 4, 9)=-DDHOOK( 4, 3)*RL
  60 :       RP( 4,10)=-RP( 4, 4)
  61 :       RP( 4,11)=RP( 4, 5)
  62 :       RP( 4,12)=RP( 4, 6)
  63 : C  MY1
  64 :       RP( 5, 5)=DDHOOK( 5, 5)*RL+DDHOOK( 3, 3)/(4.0D0*RL)
  65 :       RP( 5, 6)=DDHOOK( 5, 6)*RL
  66 :       RP( 5, 7)=-DDHOOK( 5, 1)*RL
  67 :       RP( 5, 9)=DDHOOK( 3, 3)/2.0D0
  68 :       RP( 5,10)=DDHOOK( 3, 4)/2.0D0
  69 :       RP( 5,11)=-DDHOOK( 5, 5)*RL+DDHOOK( 3, 3)/(4.0D0*RL)
  70 :       RP( 5,12)=-RP( 5, 6)
  71 : C  MZ1
  72 :       RP( 6, 6)=DDHOOK( 6, 6)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  73 :       RP( 6, 7)=-DDHOOK( 6, 1)*RL
  74 :       RP( 6, 8)=-DDHOOK( 2, 2)/2.0D0
  75 :       RP( 6,10)=-DDHOOK( 2, 4)/2.0D0
  76 :       RP( 6,11)=-DDHOOK( 6, 5)*RL
  77 :       RP( 6,12)=-DDHOOK( 6, 6)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  78 : C  FX2
  79 :       RP( 7, 7)=RP( 1, 1)
  80 :       RP( 7,11)=RP( 1, 5)
  81 :       RP( 7,12)=RP( 1, 6)
  82 : C  FY2
  83 :       RP( 8, 8)=RP( 2, 2)
  84 :       RP( 8,10)=RP( 2, 4)
  85 :       RP( 8,12)=-RP( 2, 6)
  86 : C  FZ2
  87 :       RP( 9, 9)=RP( 3, 3)
  88 :       RP( 9,10)=RP( 3, 4)
  89 :       RP( 9,11)=-RP( 3, 5)
  90 : C  MX2
  91 :       RP(10,10)=RP( 4, 4)
  92 :       RP(10,11)=-RP( 4, 5)
  93 :       RP(10,12)=-RP( 4, 6)
  94 : C  MY2
  95 :       RP(11,11)=RP( 5, 5)
  96 :       RP(11,12)=RP( 5, 6)
  97 : C  MZ2
  98 :       RP(12,12)=RP( 6, 6)
  99 : C
 100 :       DO 9 I=1,11
 101 :       DO 9 J=I+1,12
 102 :    9  RP(J,I)=RP(I,J)
 103 :       RETURN
 104 :       END
 105 :