Download tifil2.eso

Back to the list

   1 : C TIFIL2    SOURCE    CHAT      05/01/13    03:38:58     5004           
   2 : C TIFIL2    SOURCE    AM        94/02/25    22:55:47     1061           
   3 :       SUBROUTINE TIFIL2(RP,XL,LHOOK,DDHOOK)
   4 : C-----------------------------------------------------------------------
   5 : C
   6 : C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
   7 : C    POUTRE DROITE DE TIMOSHENKO 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 RIGIDITE ELEMENTAIRE 12X12
  16 : C
  17 : C D'APRES TIFILO 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 :       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, 3)=DDHOOK( 1, 3)*RL
  36 :       RP( 1, 4)=-RP( 1, 1)
  37 :       RP( 1, 6)=-RP( 1, 3)
  38 : C  FY1
  39 :       RP( 2, 2)=DDHOOK( 2, 2)*RL
  40 :       RP( 2, 3)=DDHOOK( 2, 2)/2.0D0
  41 :       RP( 2, 5)=-RP( 2, 2)
  42 :       RP( 2, 6)= RP( 2, 3)
  43 : C  MZ1
  44 :       RP( 3, 3)=DDHOOK( 3, 3)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  45 :       RP( 3, 4)=-DDHOOK( 3, 1)*RL
  46 :       RP( 3, 5)=-DDHOOK( 2, 2)/2.0D0
  47 :       RP( 3, 6)=-DDHOOK( 3, 3)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  48 : C  FX2
  49 :       RP( 4, 4)=RP( 1, 1)
  50 :       RP( 4, 6)=-RP( 1, 6)
  51 : C  FY2
  52 :       RP( 5, 5)=RP( 2, 2)
  53 :       RP( 5, 6)=-RP( 2, 3)
  54 : C  MZ2
  55 :       RP( 6, 6)=RP( 3, 3)
  56 : C
  57 :       DO 9 I=1,5
  58 :       DO 9 J=I+1,6
  59 :    9  RP(J,I)=RP(I,J)
  60 :       RETURN
  61 : 
  62 :       END
  63 :  
  64 :  

© Cast3M 2003 - All rights reserved.
Disclaimer