C UNICO1    SOURCE    CHAT      05/01/13    03:58:30     5004
C     UNICO1    SOURCE   NECH   99/05/26
      SUBROUTINE UNICO1(DK,PAEC,IC,SEQ,BETJEF,NECH0,NECH1)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*      COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,
*     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
*     &                ITER,ISIM,IBB,IGAU1,IZON
*      COMMON /DBETNECH0/DT,DC,ALFG,S0
*      COMMON /DBETNECH1/ENDL
      SEGMENT BETJEF
       REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,TDEF,
     &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
       INTEGER ICT,ICC,IMOD,IVISS,ITER,
     &                ISIM,IBB,IGAU1,IZON
      ENDSEGMENT
      SEGMENT NECH0
       REAL*8 DT,DC,ALFG,S0,ENDO
      ENDSEGMENT
      SEGMENT NECH1
       REAL*8 ENDL
      ENDSEGMENT
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
*                                                                    *
*    CALCUL DU PARAMETRE D'ECROUISSAGE                               *
*    _________________________________                               *
*                                                                    *
*    ENTREES:                                                        *
*                                                                    *
*        DK        VARIABLE D'ECROUISSAGE                            *
*        IC        =1 : COURBE de TRACTION                           *
*                  =2 : COURBE de COMPRESSION                        *
*                                                                    *
*    SORTIES:                                                        *
*                                                                    *
*        SEQ       CONTRAINTE EQUIVALENTE                            *
*        PAEC      PENTE                                             *
*                                                                    *
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
C
        BB=1.D0/AA
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C                      Direction de traction
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C
           IF (IC.EQ.1) THEN
       F=ALFA*RB
            A= -0.5D0
            B=F*(HC/GFT)*(1.D0+A/2.D0)
            C=B*LOG(1.D0-DT)/(LOG(((1+A)-SQRT(1+A*A))/(2.D0*A)))
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C                      Direction de compression
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C
           ELSE
       F=AA*RB
            A=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)
            B=F*(HC/GFC)*(1.D0+A/2.D0)
            C=B*LOG(1.D0-DC)/(LOG((1+A)/(2.D0*A)))
           ENDIF
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       CALCUL DE LA DEFORMATION PLASTIQUE LIMITE
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C
        DKL=LOG(1.D0-ENDL)/(-C)
C
           IF(DK.GE.DKL)THEN
            DK=DKL
           ENDIF
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       CALCUL LA CONTRAINTE UNIAXIALE ET DE LA PENTE
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       1: CONTRAINTE UNIAXIALE
C
            EXP1=EXP(-B*DK)
            DD=1.D0-C/B
            SEQ=F*((1.D0+A)*EXP1**(DD)-A*EXP1**(1.D0+DD))
C
C       2: CALCUL DE LA (+)*PENTE
C
            PAEC=-B*F*(DD*(1.D0+A)*EXP1**(DD)-
     &                  (1.D0+DD)*A*EXP1**(1.D0+DD))
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C
           IF(ABS(SEQ).LT.1.D-8)THEN
            SEQ=1.D-8
            PAEC=0.D0
            GOTO 10
           ENDIF
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C
   10   CONTINUE
C
        RETURN
        END