Esope source files

   1 : C TIDHL3    SOURCE    CHAT      05/01/13    03:38:25     5004           
   2 :       SUBROUTINE TIDHL3(RP,XL,YOUN,CISAIL,SD,SREDY,SREDZ,
   3 :      .                  XINX,XINY,XINZ,KERRE)
   4 : C-----------------------------------------------------------------------
   5 : C
   6 : C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
   7 : C    POUTRE DROITE DE TIMOSHENKO EN AXES LOCAUX
   8 : C    "D APRES TEDEL"
   9 : C
  10 : C  ENTREE
  11 : C    XL      LONGUEUR
  12 : C    SD      SECTION DROITE
  13 : C    SREDY   SECTION REDUITE A L EFFORT TRANCHANT TY
  14 : C    SREDZ   SECTION REDUITE A L EFFORT TRANCHANT TZ
  15 : C    XINX    INERTIE IX
  16 : C    XINY    INERTIE IY
  17 : C    XINZ    INERTIE IZ
  18 : C    YOUN    MODULE D YOUNG
  19 : C    CISAIL  MODULE DE CISAILLEMENT
  20 : C  SORTIE
  21 : C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 12X12
  22 : C
  23 : C     PP 3/9/92
  24 : C-----------------------------------------------------------------------
  25 :       IMPLICIT INTEGER(I-N)
  26 :       IMPLICIT REAL*8(A-H,O-Z)
  27 :       DIMENSION RP(12,12)
  28 :       DATA P5,P25/.5D0,.25D0/
  29 : C
  30 : C  INITIALISATIONS
  31 : C
  32 :       DO 120 I=1,12
  33 :       DO 120 J=1,12
  34 :   120 RP(I,J)=0.D00
  35 :       RL=1.D00/XL
  36 :       RL2=RL*RL
  37 : C
  38 : C  SECTION "REDUITE"
  39 : C
  40 :       IF ( SREDZ. NE.0. ) THEN
  41 :            SRZ=SREDZ
  42 :       ELSE
  43 :            SRZ=SD
  44 :       ENDIF
  45 :       IF ( SREDY. NE.0. ) THEN
  46 :            SRY=SREDY
  47 :       ELSE
  48 :            SRY=SD
  49 :       ENDIF
  50 :       IF (ABS (SRY - SRZ).GT.(MAX(SRY,SRZ)*1.D-8)) GOTO 999 
  51 :       IF (ABS (XINY - XINZ).GT.(MAX(XINY,XINZ)*1.D-8)) GOTO 999 
  52 : C
  53 :       ESL=YOUN*SD*RL
  54 :       EIZ=YOUN*XINZ
  55 :       GAZ=CISAIL*SRZ
  56 :       EIY=YOUN*XINY
  57 :       GAY=CISAIL*SRY
  58 :       GJL=CISAIL*XINX*RL
  59 : C
  60 : C
  61 : C
  62 : C
  63 : C  REMPLISSAGE DE LA MATRICE
  64 : C  FX1
  65 : C      RP(1,1)=ESL
  66 : C      RP(1,7)=-ESL
  67 : C      RP(7,1)=-ESL
  68 : C  FY1
  69 :       RP(2,3)= GAY*RL
  70 :       RP(3,2)= -GAY*RL
  71 :       RP(2,5)=  -P5*GAY
  72 :       RP(5,2)=  P5*GAY
  73 :       RP(2,9)= -GAY*RL
  74 :       RP(9,2)= GAY*RL
  75 :       RP(2,11)= -P5*GAY
  76 :       RP(11,2)= P5*GAY
  77 : C  FZ1
  78 :       RP(3,6)=  -P5*GAZ
  79 :       RP(6,3)=  P5*GAZ
  80 :       RP(3,8)= GAZ*RL
  81 :       RP(8,3)= -GAZ*RL
  82 :       RP(3,12)= -P5*GAZ
  83 :       RP(12,3)= P5*GAZ
  84 : C  MX1
  85 : C      RP(4,4)= -GJL
  86 : C      RP(4,10)=GJL
  87 : C      RP(10,4)=-GJL
  88 : C  MY1
  89 :       RP(5,6)=  P25*GAZ*XL + EIY*RL
  90 :       RP(6,5)=  -P25*GAZ*XL - EIY*RL
  91 :       RP(5,8)= -P5*GAZ
  92 :       RP(8,5)= P5*GAZ
  93 :       RP(5,12)= P25*GAZ*XL - EIY*RL
  94 :       RP(12,5)=  -P25*GAZ*XL + EIY*RL
  95 : C  MZ1
  96 :       RP(6,9)= -P5*GAY
  97 :       RP(9,6)= P5*GAY
  98 :       RP(6,11)= -P25*GAY*XL + EIZ*RL
  99 :       RP(11,6)= P25*GAY*XL - EIZ*RL
 100 : C  FX2
 101 : C      RP(7,7)=  ESL
 102 : C  FY2
 103 :       RP(8,9)= GAY*RL
 104 :       RP(9,8)= -GAY*RL
 105 :       RP(8,11)= P5*GAY
 106 :       RP(11,8)= -P5*GAY
 107 : C  FZ2
 108 :       RP(9,12)= P5*GAZ
 109 :       RP(12,9)= -P5*GAZ
 110 : C  MX2
 111 : C      RP(10,10)= -GJL
 112 : C  MY2
 113 :       RP(11,12)= P25*GAZ*XL + EIY*RL
 114 : C  MZ2
 115 :       RP(12,11)= -P25*GAY*XL - EIZ*RL
 116 : C
 117 :       RETURN
 118 : 999   KERRE=1
 119 :       RETURN
 120 :       END
 121 :  
 122 :  
 123 :