C FIBETA    SOURCE    CB215821  17/11/30    21:16:12     9639           
C FIBETA      SOURCE    AM        96/12/24    21:19:32     2448
C======================================================================
C       FIBETA  - D. COMBESCURE et P. PEGON - ELSA- 1996
C======================================================================
C
C                    MODELE GLOBAL D'ETAGE
C           (Sur des elements de poutre TIMO - Fibre
C                 (Contrainte de cisaillement/Cisail.)
C
      SUBROUTINE FIBETA(XMAT,XCAR,SIG0,VAR0,SIGF,VARF,DEPS,
     . WRK2,NCURVP,NCURVN,KERRE)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
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 DEPS     = 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
      PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16)
       DIMENSION XMAT(*),XCAR(*),SIG0(*),VAR0(*),SIGF(*),VARF(*),DEPS(*)
C
C=======================================================================
C      VARIABLES ET SEGMENTS NECESSAIRES
C=========================================================================

-INC PPARAM
-INC CCOPTIO
C  Segment de la courbe
       SEGMENT WRK2
         REAL*8 TRAC(LTRAC)
       ENDSEGMENT
C
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
        SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(2)
C
        XDDEP = DEPS(3)
        XFOR0 = SIG0(3)
C
       ELSE
C
        SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(1)
C
        XDDEP = DEPS(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) = XZER
        SIGF(2) = XZER
        SIGF(3) = XFORF
        SIGF(4) = XZER
        SIGF(5) = XZER
        SIGF(6) = XZER
       ELSE
        SIGF(1) = XZER
        SIGF(2) = XFORF
        SIGF(3) = XZER
       ENDIF
C
C=====================================================
C            FIN DE LA ROUTINE DU MODELE
C======================================================
      RETURN

      END