Esope source files

   1 : C CETAG     SOURCE    KICH      09/09/11    21:15:45     6505           
   2 : C======================================================================
   3 : C       ETAGE  - D. COMBESCURE et P. PEGON - ELSA- 1996
   4 : C======================================================================
   5 : C
   6 : C                    MODELE GLOBAL D'ETAGE 
   7 : C           (Sur des elements de poutre TIMO - Effort tranchant/Cisail.)
   8 : C
   9 :       SUBROUTINE CETAG(WRK52,WRK53,WRK54,WRK2,NCURVP,NCURVN)
  10 : c      SUBROUTINE ETAGE(WRK0,WRK1,WRK2,NCURVP,NCURVN,KERRE)
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)
  13 : C ETAGE     SOURCE    DC        98/11/10    21:19:23     3351           
  14 : C
  15 : C=======================================================================
  16 : C             CETTE ROUTINE EST APPELE DANS ECO100
  17 : C
  18 : C
  19 : C WRK0      = Segment materiaux
  20 : C WRK1      = Segment contraintes
  21 : C WRK2      = Segment courbe
  22 : C NCURVP    = Longueur courbe sens positif
  23 : C NCURVN    = Longueur courbe sens negatif
  24 : C
  25 : C======================================================================
  26 : C XCAR      = Caracteristique de la section 
  27 : C DEPST     = Increment de deformation axiale
  28 : C SIG0      = Contrainte initiale
  29 : C VAR0      = Variables internes initiales
  30 : C SIGF      = Contrainte finale
  31 : C VARF      = Variables internes finales
  32 : C DEFP      = Deformation plastique
  33 : C
  34 : C=======================================================================
  35 : C
  36 :       REAL*8 INRX,INRY,INRZ
  37 :       PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16)
  38 : C
  39 : C=======================================================================
  40 : C      VARIABLES ET SEGMENTS NECESSAIRES
  41 : C=========================================================================
  42 : -INC CCOPTIO
  43 : -INC DECHE
  44 : C  Segment de la courbe
  45 :        SEGMENT WRK2
  46 :          REAL*8 TRAC(LTRAC)
  47 :        ENDSEGMENT
  48 : C
  49 : C  Lecture materiau
  50 : C
  51 :        YOUN = XMAT(1)
  52 :        XNU = XMAT(2)
  53 :        XDELAP = XMAT(5)
  54 :        XDELAN = (-1.D0)*XMAT(6)
  55 :        XDMAXP = XMAT(7)
  56 :        XDMAXN = XMAT(8)
  57 :        XBETA  = XMAT(9)
  58 :        XALPH  = XMAT(10)
  59 :        XTETA  = XMAT(11)
  60 : C
  61 : C
  62 :        IF (IDIM.EQ.3) THEN
  63 :         INRX = YOUN*xcarb(1)
  64 :         INRY = YOUN*xcarb(2)
  65 :         INRZ = YOUN*xcarb(3)
  66 :         SECT = YOUN*xcarb(4)
  67 :         SRDY = (YOUN/(2.D0*(1.D0+XNU)))*xcarb(5)
  68 :         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*xcarb(6)
  69 : C
  70 :         XDDEP = DEPST(3)
  71 :         XFOR0 = SIG0(3)
  72 : C
  73 :        ELSE
  74 : C
  75 :         SECT = YOUN*xcarb(1)
  76 :         INRZ = YOUN*xcarb(2)
  77 :         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*xcarb(3)
  78 : C
  79 :         XDDEP = DEPST(2)
  80 :         XFOR0 = SIG0(2)
  81 :        ENDIF
  82 : C
  83 : C Appel de la loi 
  84 : C
  84 :         CALL ETAGEL(XDDEP,XFOR0,XFORF,VAR0,VARF,SECZ,
  85 :      & XDELAP,XDELAN,XDMAXP,XDMAXN,XBETA,XALPH,XTETA,
  86 :      & WRK2,NCURVP,NCURVN,KERRE)
  87 : C
  88 : C==================================================================
  89 : C
  90 : C  On remplit les tableaux avant de sortir
  91 : C
  92 : C==================================================================
  93 :        IF (IDIM.EQ.3) THEN
  94 :         SIGF(1) = SIG0(1) + SECT*DEPST(1) 
  95 :         SIGF(2) = SIG0(2) + SRDY*DEPST(2)
  96 :         SIGF(3) = XFORF
  97 :         SIGF(4) = SIG0(4) + INRX*DEPST(4)
  98 :         SIGF(5) = SIG0(5) + INRY*DEPST(5)
  99 :         SIGF(6) = SIG0(6) + INRZ*DEPST(6)
 100 :        ELSE
 101 :         SIGF(1) = SIG0(1) + SECT*DEPST(1) 
 102 :         SIGF(2) = XFORF
 103 :         SIGF(3) = SIG0(3) + INRZ*DEPST(3)
 104 :        ENDIF
 105 : C
 106 : C=====================================================
 107 : C            FIN DE LA ROUTINE DU MODELE
 108 : C======================================================
 109 :       RETURN
 110 : 
 111 :       END
 112 :  
 113 :  
 114 :  
 115 :  
 116 :  
 117 :  
 118 :  
 119 :