Esope source files

   1 : C TIRILO    SOURCE    CHAT      05/01/13    03:41:32     5004           
   2 :       SUBROUTINE TIRILO(RP,XL,YOUN,CISAIL,SD,SREDY,SREDZ,
   3 :      .                  XINX,XINY,XINZ)
   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,144
  33 :   120 RP(I,1)=0.D00
  34 :       RL=1.D00/XL
  35 :       RL2=RL*RL
  36 : C
  37 : C  SECTION "REDUITE"
  38 : C
  39 :       IF ( SREDZ. NE.0. ) THEN
  40 :            SRZ=SREDZ
  41 :       ELSE
  42 :            SRZ=SD
  43 :       ENDIF
  44 :       IF ( SREDY. NE.0. ) THEN
  45 :            SRY=SREDY
  46 :       ELSE
  47 :            SRY=SD
  48 :       ENDIF
  49 : C
  50 :       ESL=YOUN*SD*RL
  51 :       EIZ=YOUN*XINZ
  52 :       GAZ=CISAIL*SRZ
  53 :       EIY=YOUN*XINY
  54 :       GAY=CISAIL*SRY
  55 :       GJL=CISAIL*XINX*RL
  56 : C
  57 : C  REMPLISSAGE DE LA MATRICE
  58 : C  FX1
  59 :        RP(1,1)=ESL
  60 :       RP(1,7)=-ESL
  61 : C  FY1
  62 :        RP(2,2)= GAY*RL
  63 :       RP(2,6)=  P5*GAY
  64 :        RP(2,8)= -GAY*RL
  65 :       RP(2,12)= P5*GAY
  66 : C  FZ1
  67 :        RP(3,3)= GAZ*RL
  68 :       RP(3,5)=  -P5*GAZ
  69 :        RP(3,9)= -GAZ*RL
  70 :       RP(3,11)= -P5*GAZ
  71 : C  MX1
  72 :        RP(4,4)= GJL
  73 :       RP(4,10)=-GJL
  74 : C  MY1
  75 :       RP(5,5)=  P25*GAZ*XL + EIY*RL
  76 :        RP(5,9)= P5*GAZ
  77 :       RP(5,11)= P25*GAZ*XL - EIY*RL
  78 : C  MZ1
  79 :       RP(6,6)=  P25*GAY*XL + EIZ*RL
  80 :        RP(6,8)= -P5*GAY
  81 :       RP(6,12)= P25*GAY*XL - EIZ*RL
  82 : C  FX2
  83 :       RP(7,7)=  ESL
  84 : C  FY2
  85 :        RP(8,8)= GAY*RL
  86 :       RP(8,12)= -P5*GAY
  87 : C  FZ2
  88 :        RP(9,9)= GAZ*RL
  89 :       RP(9,11)= P5*GAZ
  90 : C  MX2
  91 :       RP(10,10)= GJL
  92 : C  MY2
  93 :       RP(11,11)= P25*GAZ*XL + EIY*RL
  94 : C  MZ2
  95 :       RP(12,12)= P25*GAY*XL + EIZ*RL
  96 : C
  97 :       DO 9 I=1,12
  98 :       DO 9 J=I,12
  99 :    9  RP(J,I)=RP(I,J)
 100 :       RETURN
 101 :       END
 102 :