Esope source files

   1 : C TIRIL2    SOURCE    CHAT      05/01/13    03:41:29     5004           
   2 : C TIRILO    SOURCE    AM        93/01/12    21:36:57     781            
   3 :       SUBROUTINE TIRIL2(RP,XL,YOUN,CISAIL,SD,SREDY,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    XINZ    INERTIE IZ
  15 : C    YOUN    MODULE D YOUNG
  16 : C    CISAIL  MODULE DE CISAILLEMENT
  17 : C  SORTIE
  18 : C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 6X6
  19 : C
  20 : C     D'PRES TIRILO DC 98
  21 : C-----------------------------------------------------------------------
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8(A-H,O-Z)
  24 :       DIMENSION RP(6,6)
  25 :       DATA P5,P25/.5D0,.25D0/
  26 : C
  27 : C  INITIALISATIONS
  28 : C
  29 :       DO 120 I=1,36
  30 :   120 RP(I,1)=0.D00
  31 :       RL=1.D00/XL
  32 :       RL2=RL*RL
  33 : C
  34 : C  SECTION "REDUITE"
  35 : C
  36 :       IF ( SREDY. NE.0. ) THEN
  37 :            SRY=SREDY
  38 :       ELSE
  39 :            SRY=SD
  40 :       ENDIF
  41 : C
  42 :       ESL=YOUN*SD*RL
  43 :       EIZ=YOUN*XINZ
  44 :       GAY=CISAIL*SRY
  45 : C
  46 : C  REMPLISSAGE DE LA MATRICE
  47 : C  FX1
  48 :        RP(1,1)=ESL
  49 :        RP(1,4)=-ESL
  50 : C  FY1
  51 :        RP(2,2)= GAY*RL
  52 :        RP(2,3)=  P5*GAY
  53 :        RP(2,5)= -GAY*RL
  54 :        RP(2,6)= P5*GAY
  55 : C  MZ1
  56 :       RP(3,3)=  P25*GAY*XL + EIZ*RL
  57 :        RP(3,5)= -P5*GAY
  58 :       RP(3,6)= P25*GAY*XL - EIZ*RL
  59 : C  FX2
  60 :       RP(4,4)=  ESL
  61 : C  FY2
  62 :        RP(5,5)= GAY*RL
  63 :       RP(5,6)= -P5*GAY
  64 : C  MZ2
  65 :       RP(6,6)= P25*GAY*XL + EIZ*RL
  66 : C
  67 :       DO 9 I=1,6
  68 :       DO 9 J=I,6
  69 :    9  RP(J,I)=RP(I,J)
  70 :       RETURN
  71 :       END
  72 :  
  73 :