Esope source files

   1 : C TIFDL3    SOURCE    CHAT      05/01/13    03:38:51     5004           
   2 :        SUBROUTINE TIFDL3(RP,XL,LHOOK,DDHOOK,OMEG)
   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 :       GA = OMEG*0.5D0*(DDHOOK( 2, 2) + DDHOOK( 3, 3))
  33 :       EI = OMEG*0.5D0*(DDHOOK( 5, 5) + DDHOOK( 6, 6))
  34 : C
  35 : C  REMPLISSAGE DE LA MATRICE
  36 : C  FX1
  37 : C      RP( 1, 1)=DDHOOK( 1, 1)*RL
  38 : C      RP( 1, 7)=-RP( 1, 1)
  39 : C      RP( 7, 1)=-RP( 1, 1)
  40 : C  FY1
  41 :       RP( 2, 3)=GA*RL
  42 :       RP( 3, 2)=-GA*RL
  43 :       RP( 2, 5)=-GA/2.0D0
  44 :       RP( 5, 2)=GA/2.0D0
  45 :       RP( 2, 9)=-GA*RL
  46 :       RP( 9, 2)=GA*RL
  47 :       RP( 2,11)=-GA/2.0D0
  48 :       RP(11, 2)=GA/2.0D0
  49 : C  FZ1
  50 :       RP( 3, 6)=-GA/2.0D0
  51 :       RP( 3, 8)=GA*RL
  52 :       RP( 3, 12)=-GA/2.0D0
  53 :       RP( 6, 3)=GA/2.0D0
  54 :       RP( 8, 3)=-GA*RL
  55 :       RP( 12, 3)=GA/2.0D0
  56 : C  MX1
  57 : C      RP( 4, 4)=DDHOOK( 4, 4)*RL
  58 : C      RP( 4,10)=-RP( 4, 4)
  59 : C      RP( 10,4)=-RP( 4, 4)
  60 : C  MY1
  61 :       RP( 5, 6)=EI*RL+GA/(4.0D0*RL)
  62 :       RP( 5, 8)=-GA/2.0D0
  63 :       RP( 5,12)=-EI*RL+GA/(4.0D0*RL)
  64 :       RP( 6, 5)=-EI*RL-GA/(4.0D0*RL)
  65 :       RP( 8, 5)=GA/2.0D0
  66 :       RP(12, 5)=EI*RL-GA/(4.0D0*RL)
  67 : C  MZ1
  68 :       RP( 6, 9)=-GA/2.0D0
  69 :       RP( 6,11)=EI*RL-GA/(4.0D0*RL)
  70 :       RP( 9, 6)=GA/2.0D0
  71 :       RP( 11, 6)=-EI*RL+GA/(4.0D0*RL)
  72 : C  FX2
  73 : C      RP( 7, 7)=RP( 1, 1)
  74 : C  FY2
  75 :       RP( 8, 9)=GA*RL
  76 :       RP( 9, 8)=-GA*RL
  77 :       RP( 8, 11)=GA/2.0D0
  78 :       RP( 11, 8)=-GA/2.0D0
  79 : C  FZ2
  80 :       RP( 9, 12)=GA/2.0D0
  81 :       RP( 12, 9)=-GA/2.0D0
  82 : C  MX2
  83 : C      RP(10,10)=RP( 4, 4)
  84 : C  MY2
  85 :       RP(11,12)=EI*RL+GA/(4.0D0*RL)
  86 : C  MZ2
  87 :       RP(12,11)=-EI*RL-GA/(4.0D0*RL)
  88 : C
  89 :       RETURN
  90 :       END
  91 :  
  92 :  
  93 :