Esope source files

   1 : C VALINF    SOURCE    CHAT      05/01/13    04:00:36     5004           
   2 :       SUBROUTINE VALINF(X,XDEGRA,CURVE,NCURV,YOUN,Y,KERRE)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 : C============================================================
   6 : C 25/1/96  D. C. Modif de la routine de Pierre
   7 : C  Routine qui lit la valeur d'une courbe  
   8 : C  quand on lui donne l'abscisse
   9 : C  Adaptation au cas d'une famille de courbe dependant d'un parametre
  10 : C
  11 : C  ENTREES:
  12 : C    X     : Abscisse (deformation totale -E>0- ou plastique-Ereel)
  13 : C    XDEGRA: Parametre definissant la famille de courbe
  14 : C    CURVE : Courbe
  15 : C    NCURV : Nombre de points de la courbe
  16 : C    YOUN  : Module de decharge elastique E
  17 : C     
  18 : C SORTIES
  19 : C    Y     : Valeur de la fonction en X
  20 : C=============================================================
  21 :       PARAMETER (XZER = 0.D0, EPSILO=1.D-20)
  22 :       DIMENSION CURVE(2,NCURV)
  23 : C
  24 :       Y00 = XDEGRA * CURVE(2,1)
  25 :       IF (ABS(YOUN).GT.EPSILO) THEN
  26 :          X00 = (Y00/YOUN)
  27 :       ELSE
  28 :          X00 = XZER
  29 :       ENDIF
  30 :       DO IE1=2,NCURV
  31 :         YC = XDEGRA * CURVE(2,IE1)
  32 :         XC = CURVE(1,IE1) 
  33 :         IF (ABS(YOUN).GT.EPSILO) THEN
  34 :           XBIS = XC + (YC/YOUN)
  35 :         ELSE
  36 :           XBIS = XC 
  37 :         ENDIF 
  38 :         IF (X.LE.XBIS) GOTO 1
  39 :         Y00 = YC
  40 :         X00 = XBIS 
  41 :       ENDDO
  42 :       Y = Y00
  43 :       GOTO 2
  44 :  1    Y = Y00 + ((YC - Y00)/(XBIS - X00))*(X - X00) 
  45 : C
  46 :  2    CONTINUE
  47 :       RETURN
  48 :       END
  49 : 
  50 :  
  51 :