C COUPLA    SOURCE    CHAT      05/01/12    22:24:58     5004           
C COUPLA     SOURCE    INSL       24/10/96    
      SUBROUTINE COUPLA(EPS1,COURBE,IPCOUR,YOUN,SIGX,PENTE)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION COURBE(IPCOUR)
C
C         **************************************
C         *  CAS ELASTOPLASTIQUE ECROUISSABLE  *
C         *  COURBE INTRODUITE POINT PAR POINT *
C         *  CALCUL DE LA PENTE DE LA COURBE   *
C         **************************************
C
      A1=1.D0
      IF(EPS1.LT.0.D0) A1=-1.D0
      EPSX=DABS(EPS1)
      SIGX1=DABS(SIGX)
C
      IF(EPSX.LE.COURBE(4).AND.SIGX1.LT.COURBE(3)) THEN
C          --------------------------------------------
C          * PREMIER TRONCON DE LA COURBE DE TRACTION *
C          --------------------------------------------
         PENTE=YOUN
         SIGX=A1*YOUN*EPSX
         EPPX=0.D0
      ELSE
C
       DO 1 I=6,IPCOUR,2
C          --------------------------------------------------------
C          * RECHERCHE DU POINT DE LA COURBE CORRESPONDANT A SET *
C          --------------------------------------------------------
         IF(EPSX.LE.COURBE(I)) THEN
            A=COURBE(I-3)
            B=COURBE(I-2)
            C=COURBE(I-1)
            D=COURBE(I)
            GOTO 3
         ENDIF
    1  CONTINUE
C
       PENTE=0.D0
       SIGX=0.D0
       EPPX=EPS1-SIGX/YOUN
       GOTO 4
C     ------------------------------------------------------------
C     * PENTEC = PENTE DE LA COURBE DE TRACTION (BP= A-B*PENTEC) *
C     ------------------------------------------------------------
    3  CONTINUE
       PENTE=(C-A)/(D-B)
       SIGX=A1*(A+(PENTE*(EPSX-B)))
       EPPX=EPS1-SIGX/YOUN
      ENDIF
   4  CONTINUE
C
      RETURN
      END