Esope source files

   1 : C DFHARDCC  SOURCE    AM        12/02/27    23:39:36     7289           
   2 :       FUNCTION dfhardCC(rkappac,fc,fc0,rkappa1,ac,bc)
   3 : 
   4 : c     This subroutine calculates the derivative of the compressive 
   5 : c     hardening function with respect to rkappac
   6 : 
   7 :       IMPLICIT REAL*8 (A-B,D-H,O-Z)
   8 :       implicit integer (I-K,M,N)
   9 :       implicit logical (L)
  10 :       implicit character*10 (C)
  11 : 
  12 :       r0 = 0.
  13 :       hc = bc
  14 :    
  15 :       if (rkappac.le.rkappa1) then
  16 :         rn = (rkappac/rkappa1)
  17 :         rloc = 2.*(fc-fc0)*rn/(1.0d0+(rn*rn))
  18 :         ry = (dexp(ac*rkappac))*(fc0+rloc)
  19 :         rloc2 = 2.*(fc-fc0)*(1.0d0-(rn*rn))/
  20 :      .          (rkappa1*(1.0d0+(rn*rn))*(1.0d0+(rn*rn)))
  21 :         dry = ac*ry + (dexp(ac*rkappac))*rloc2
  22 :       else
  23 :         rloc = (1.0d0 + hc*(rkappac-rkappa1)) * 
  24 :      .    (dexp(-bc*(rkappac-rkappa1)))
  25 :         ry = (dexp(ac*rkappac))*fc*rloc
  26 :         rloc2 = fc * (hc-bc-hc*bc*(rkappac-rkappa1))
  27 :         dry = ac*ry + (dexp((ac-bc)*rkappac + bc*rkappa1)) * rloc2
  28 :       endif
  29 : 
  30 :       dfhardcc = dry
  31 : 
  32 :       RETURN
  33 :       END
  34 : 
  35 : 
  36 :