Esope source files

   1 : C FHARDCC   SOURCE    AM        12/02/28    00:28:17     7289           
   2 :       FUNCTION fhardCC(rkappac,fc,fc0,rkappa1,ac,bc,lerror)
   3 : 
   4 : c     hardening in effective stress!
   5 : c     strain hardening inspired from concrete ETC curve in uniaxial compression
   6 : c     for the ascending branch, i.e. until peak stress and peak stress strain
   7 : c     then the descending branch is an exponential  
   8 :         
   9 : c     Elastic from (0;0) to (0;fc0)
  10 : c     The curve begins on (0;fc0),
  11 : c       then ascending branch based on ETC until peak
  12 : c       then descending branch exponential 
  13 : c       and at the junction, the tangent is horizontal
  14 : 
  15 :       IMPLICIT REAL*8 (A-B,D-H,O-Z)
  16 :       implicit integer (I-K,M,N)
  17 :       implicit logical (L)
  18 :       implicit character*10 (C)
  19 : 
  20 :       hc = bc
  21 : 
  22 :       if (rkappac.le.rkappa1) then
  23 :         rn = (rkappac/rkappa1)
  24 :         rloc = 2.*(fc-fc0)*rn/(1.+(rn*rn))
  25 :         ry = (dexp(ac*rkappac))*(fc0+rloc)
  26 :       else
  27 :         rloc = (1. + hc*(rkappac-rkappa1)) * 
  28 :      .   (dexp(-bc*(rkappac-rkappa1)))
  29 :         ry = (dexp(ac*rkappac))*fc*rloc
  30 :       endif
  31 : 
  32 :       if (dabs(ry).gt.10.E+25) then
  33 :         lerror = .true.
  34 :         return
  35 :       endif
  36 :       
  37 :       fhardCC= ry
  38 :    
  39 :       RETURN
  40 :       END
  41 : 
  42 :