Sources Esope | Cast3M

cetag
C CETAG     SOURCE    PV        17/12/08    21:15:35     9660           
C======================================================================
C       ETAGE  - D. COMBESCURE et P. PEGON - ELSA- 1996
C======================================================================
C
C                    MODELE GLOBAL D'ETAGE
C           (Sur des elements de poutre TIMO - Effort tranchant/Cisail.)
C
      SUBROUTINE CETAG(WRK52,WRK53,WRK54,WRK2,NCURVP,NCURVN)
c      SUBROUTINE ETAGE(WRK0,WRK1,WRK2,NCURVP,NCURVN,KERRE)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C ETAGE     SOURCE    DC        98/11/10    21:19:23     3351
C
C=======================================================================
C             CETTE ROUTINE EST APPELE DANS ECO100
C
C
C WRK0      = Segment materiaux
C WRK1      = Segment contraintes
C WRK2      = Segment courbe
C NCURVP    = Longueur courbe sens positif
C NCURVN    = Longueur courbe sens negatif
C
C======================================================================
C XCAR      = Caracteristique de la section
C DEPST     = Increment de deformation axiale
C SIG0      = Contrainte initiale
C VAR0      = Variables internes initiales
C SIGF      = Contrainte finale
C VARF      = Variables internes finales
C DEFP      = Deformation plastique
C
C=======================================================================
C
      REAL*8 INRX,INRY,INRZ
      PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16)
C
C=======================================================================
C      VARIABLES ET SEGMENTS NECESSAIRES
C=========================================================================
 
-INC PPARAM
-INC CCOPTIO
-INC DECHE
C  Segment de la courbe
       SEGMENT WRK2
         REAL*8 TRAC(LTRAC)
       ENDSEGMENT
C
C  Lecture materiau
C
       YOUN = XMAT(1)
       XNU = XMAT(2)
       XDELAP = XMAT(5)
       XDELAN = (-1.D0)*XMAT(6)
       XDMAXP = XMAT(7)
       XDMAXN = XMAT(8)
       XBETA  = XMAT(9)
       XALPH  = XMAT(10)
       XTETA  = XMAT(11)
C
C
       IF (IDIM.EQ.3) THEN
        INRX = YOUN*xcarb(1)
        INRY = YOUN*xcarb(2)
        INRZ = YOUN*xcarb(3)
        SECT = YOUN*xcarb(4)
        SRDY = (YOUN/(2.D0*(1.D0+XNU)))*xcarb(5)
        SECZ = (YOUN/(2.D0*(1.D0+XNU)))*xcarb(6)
C
        XDDEP = DEPST(3)
        XFOR0 = SIG0(3)
C
       ELSE
C
        SECT = YOUN*xcarb(1)
        INRZ = YOUN*xcarb(2)
        SECZ = (YOUN/(2.D0*(1.D0+XNU)))*xcarb(3)
C
        XDDEP = DEPST(2)
        XFOR0 = SIG0(2)
       ENDIF
C
C Appel de la loi
C
        CALL ETAGEL(XDDEP,XFOR0,XFORF,VAR0,VARF,SECZ,
     & XDELAP,XDELAN,XDMAXP,XDMAXN,XBETA,XALPH,XTETA,
     & WRK2,NCURVP,NCURVN,KERRE)
C
C==================================================================
C
C  On remplit les tableaux avant de sortir
C
C==================================================================
       IF (IDIM.EQ.3) THEN
        SIGF(1) = SIG0(1) + SECT*DEPST(1)
        SIGF(2) = SIG0(2) + SRDY*DEPST(2)
        SIGF(3) = XFORF
        SIGF(4) = SIG0(4) + INRX*DEPST(4)
        SIGF(5) = SIG0(5) + INRY*DEPST(5)
        SIGF(6) = SIG0(6) + INRZ*DEPST(6)
       ELSE
        SIGF(1) = SIG0(1) + SECT*DEPST(1)
        SIGF(2) = XFORF
        SIGF(3) = SIG0(3) + INRZ*DEPST(3)
       ENDIF
C
C=====================================================
C            FIN DE LA ROUTINE DU MODELE
C======================================================
      RETURN
 
      END