Esope source files

   1 : C TIMOLO    SOURCE    CHAT      05/01/13    03:40:40     5004           
   2 :       SUBROUTINE TIMOLO(V,CONT,XL,YOUN,CISAIL,SD,SREDY,SREDZ,
   3 :      .                  XINX,XINY,XINZ,IDEUX)
   4 : C-----------------------------------------------------------------------
   5 : C
   6 : C    ROUTINE DE CALCUL DES CONTRAINTES DE LA POUTRE DROITE
   7 : C    D APRES TEDEL
   8 : C
   9 : C  ENTREE
  10 : C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
  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    IDEUX   INDICE POUR DEFORMATIONS DU SECOND ORDRE
  21 : C  SORTIE
  22 : C    CONT    CONTRAINTES AU MILIEU 6 PAR NOEUD
  23 : C
  24 : C    PP 4/9/92
  25 : C-----------------------------------------------------------------------
  26 :       IMPLICIT INTEGER(I-N)
  27 :       IMPLICIT REAL*8(A-H,O-Z)
  28 :       DIMENSION CONT(*),V(*)
  29 :       DATA P5/0.5D0/
  30 : C
  31 :       RL=1.D00/XL
  32 :       RL2=RL*RL
  33 :       RL3= RL2*RL
  34 : C
  35 : C  SECTION "REDUITE"
  36 : C
  37 :       IF ( SREDZ. NE.0. ) THEN
  38 :            SRZ=SREDZ
  39 :       ELSE
  40 :            SRZ=SD
  41 :       ENDIF
  42 :       IF ( SREDY. NE.0. ) THEN
  43 :            SRY=SREDY
  44 :       ELSE
  45 :            SRY=SD
  46 :       ENDIF
  47 : C
  48 : C  DEFORMATIONS GENERALISEES
  49 : C
  50 :       EPS  = (V(7)-V(1))*RL
  51 :       EPTY = (V(8)-V(2))*RL - (V(6)+V(12))*P5
  52 :       EPTZ = (V(9)-V(3))*RL + (V(5)+V(11))*P5
  53 :       XSIX = (V(10)-V(4))*RL
  54 :       XSIY =-(V(5)-V(11))*RL
  55 :       XSIZ =-(V(6)-V(12))*RL
  56 : C
  57 : C  TERMES DU SECOND ORDRE pour la membrane
  58 : C
  59 :       IF(IDEUX.NE.0) THEN
  60 :       AUX1=0.5D0*(EPS**2+((V(8)-V(2))*RL)**2 +((V(9)-V(3))*RL)**2)
  61 :       EPS = EPS + AUX1
  62 :       ENDIF
  63 : C
  64 : C  CONTRAINTES GENERALISEES
  65 : C
  66 :       SIG  =YOUN*SD*    EPS
  67 :       TRANY=CISAIL*SRY* EPTY
  68 :       TRANZ=CISAIL*SRZ* EPTZ
  69 :       TOR  =CISAIL*XINX*XSIX
  70 :       XMOMY=YOUN*XINY*  XSIY
  71 :       XMOMZ=YOUN*XINZ*  XSIZ
  72 : C
  73 : C  REMPLISSAGE DU TABLEAU CONT
  74 : C
  75 :       CONT(1)=SIG
  76 :       CONT(2)=TRANY
  77 :       CONT(3)=TRANZ
  78 :       CONT(4)=TOR
  79 :       CONT(5)=XMOMY
  80 :       CONT(6)=XMOMZ
  81 :       RETURN
  82 :       END
  83 :  
  84 :