Esope source files

   1 : C TIMOL2    SOURCE    CHAT      05/01/13    03:40:35     5004           
   2 : C TIMOL2    SOURCE    IANIS     95/09/18    21:20:16     1823           
   3 :       SUBROUTINE TIMOL2(V,CONT,XL,YOUN,CISAIL,SD,SREDY,XINZ,IDEUX)
   4 : C-----------------------------------------------------------------------
   5 : C
   6 : C    ROUTINE DE CALCUL DES CONTRAINTES DE LA POUTRE DROITE EN 2D
   7 : C    D APRES TEDEL
   8 : C
   9 : C  ENTREE
  10 : C    XL      LONGUEUR
  11 : C    SD      SECTION DROITE
  12 : C    SREDY   SECTION REDUITE A L EFFORT TRANCHANT TY
  13 : C    XINZ    INERTIE IZ
  14 : C    YOUN    MODULE D YOUNG
  15 : C    CISAIL  MODULE DE CISAILLEMENT
  16 : C    IDEUX   INDICE POUR DEFORMATIONS DU SECOND ORDRE
  17 : C  SORTIE
  18 : C    CONT    CONTRAINTES AU MILIEU 6 PAR NOEUD
  19 : C
  20 : C    DC 14/4/98 d'apres TIMOL2
  21 : C-----------------------------------------------------------------------
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8(A-H,O-Z)
  24 :       DIMENSION CONT(*),V(*)
  25 :       DATA P5/0.5D0/
  26 : C
  27 :       RL=1.D00/XL
  28 :       RL2=RL*RL
  29 :       RL3= RL2*RL
  30 : C
  31 : C  SECTION "REDUITE"
  32 : C
  33 :       IF ( SREDY. NE.0. ) THEN
  34 :            SRY=SREDY
  35 :       ELSE
  36 :            SRY=SD
  37 :       ENDIF
  38 : C
  39 : C  DEFORMATIONS GENERALISEES
  40 : C
  41 :       EPS  = (V(4)-V(1))*RL
  42 :       EPTY = (V(5)-V(2))*RL - (V(3)+V(6))*P5
  43 :       XSIZ =-(V(3)-V(6))*RL
  44 : C
  45 : C  TERMES DU SECOND ORDRE pour la membrane
  46 : C
  47 :       IF(IDEUX.NE.0) THEN
  48 :       AUX1=0.5D0*(EPS**2+((V(5)-V(2))*RL)**2)
  49 :       EPS = EPS + AUX1
  50 :       ENDIF
  51 : C
  52 : C  CONTRAINTES GENERALISEES
  53 : C
  54 :       SIG  =YOUN*SD*    EPS
  55 :       TRANY=CISAIL*SRY* EPTY
  56 :       XMOMZ=YOUN*XINZ*  XSIZ
  57 : C
  58 : C  REMPLISSAGE DU TABLEAU CONT
  59 : C
  60 :       CONT(1)=SIG
  61 :       CONT(2)=TRANY
  62 :       CONT(3)=XMOMZ
  63 :       RETURN
  64 :       END
  65 :  
  66 :  
  67 :