Esope source files

   1 : C COUPLA    SOURCE    CHAT      05/01/12    22:24:58     5004           
   2 : C COUPLA     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE COUPLA(EPS1,COURBE,IPCOUR,YOUN,SIGX,PENTE)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 :       DIMENSION COURBE(IPCOUR)
   8 : C
   9 : C         **************************************
  10 : C         *  CAS ELASTOPLASTIQUE ECROUISSABLE  *
  11 : C         *  COURBE INTRODUITE POINT PAR POINT *
  12 : C         *  CALCUL DE LA PENTE DE LA COURBE   *
  13 : C         **************************************
  14 : C
  15 :       A1=1.D0
  16 :       IF(EPS1.LT.0.D0) A1=-1.D0
  17 :       EPSX=DABS(EPS1)
  18 :       SIGX1=DABS(SIGX)
  19 : C
  20 :       IF(EPSX.LE.COURBE(4).AND.SIGX1.LT.COURBE(3)) THEN
  21 : C          --------------------------------------------
  22 : C          * PREMIER TRONCON DE LA COURBE DE TRACTION *
  23 : C          --------------------------------------------
  24 :          PENTE=YOUN
  25 :          SIGX=A1*YOUN*EPSX
  26 :          EPPX=0.D0
  27 :       ELSE
  28 : C
  29 :        DO 1 I=6,IPCOUR,2
  30 : C          --------------------------------------------------------
  31 : C          * RECHERCHE DU POINT DE LA COURBE CORRESPONDANT A SET *
  32 : C          --------------------------------------------------------
  33 :          IF(EPSX.LE.COURBE(I)) THEN
  34 :             A=COURBE(I-3)
  35 :             B=COURBE(I-2)
  36 :             C=COURBE(I-1)
  37 :             D=COURBE(I)
  38 :             GOTO 3
  39 :          ENDIF
  40 :     1  CONTINUE
  41 : C
  42 :        PENTE=0.D0
  43 :        SIGX=0.D0
  44 :        EPPX=EPS1-SIGX/YOUN
  45 :        GOTO 4
  46 : C     ------------------------------------------------------------
  47 : C     * PENTEC = PENTE DE LA COURBE DE TRACTION (BP= A-B*PENTEC) *
  48 : C     ------------------------------------------------------------
  49 :     3  CONTINUE
  50 :        PENTE=(C-A)/(D-B)
  51 :        SIGX=A1*(A+(PENTE*(EPSX-B)))
  52 :        EPPX=EPS1-SIGX/YOUN
  53 :       ENDIF
  54 :    4  CONTINUE
  55 : C
  56 :       RETURN
  57 :       END
  58 :  
  59 :  
  60 :