Esope source files

   1 : C EPSEQU    SOURCE    CHAT      05/01/12    23:41:24     5004           
   2 : C EPSEQU     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE EPSEQU(EPEQ,SEQ,IPLA,ICAL,IPOS,EX,RB,EPO1,EPO,EMAX)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : C
   8 : C                  *******************************
   9 : C                  *** COURBE DE COMPRESSION   ***
  10 : C                  ***     UNIAXIALE           ***
  11 : C                  *******************************
  12 : C-------------------------------------------------------------------
  13 :       TU=RB
  14 :       EPO2=EPO
  15 :       IF(IPOS.EQ.1) EPO2=EPO1
  16 :       IF(SEQ.GE.TU) THEN
  17 :          SEQ =TU
  18 :          EPEQ=EPO2
  19 :          IPLA=2
  20 :          GOTO 100
  21 :       ENDIF
  22 : C-------------------------------------------------------------------
  23 : C               PRE-PIC  ET POST-PIC LINEAIRES
  24 : C-------------------------------------------------------------------
  25 :       IF(IPLA.NE.0) THEN
  26 :         EPEQ=EMAX-SEQ*(EMAX-EPO2)/TU
  27 :         GOTO 100
  28 :       ENDIF
  29 :       IF(ICAL.NE.0) THEN
  30 :         EPEQ=SEQ/EX
  31 :         GOTO 100
  32 :       ENDIF
  33 : C-------------------------------------------------------------------
  34 : C                INSA_LBS1  ( PRE-PIC HYPERBOLE)
  35 : C-------------------------------------------------------------------
  36 :         EPEQ=0.D0
  37 :         IF(SEQ.EQ.0.D0.OR.EPO2.EQ.0.D0.OR.TU.EQ.0.D0) GOTO 100
  38 : C
  39 :         CC3=EPO2*EPO2
  40 :         BB3=CC3*(EX/SEQ-EX/TU+2.D0/EPO2)
  41 :         DELT=BB3*BB3-4.D0*CC3
  42 : C
  43 :         IF(DELT.EQ.0.D0) THEN
  44 :           EPEQ=EPO2
  45 :           GOTO 100
  46 :         ENDIF
  47 :         IF(DELT.GT.0.D0) EPEQ=(BB3-DSQRT(DELT))/2.D0
  48 : C-------------------------------------------------------------------
  49 :   100 CONTINUE
  50 :       RETURN
  51 :       END
  52 :  
  53 :