Esope source files

   1 : C DPFUNC    SOURCE    AM        12/02/27    23:39:37     7289           
   2 :       FUNCTION DPfunc(sigma,alpha)
   3 : 
   4 : c     Drucker Prager function to describe the compressive part of a concrete material.
   5 : 
   6 : c     sigmaeq = dsqrt (3 J2) + alpha I1  with  J2 : 2nd invariant of the deviatoric stress                                           
   7 : c                                        and   I1 : 1st invariant of the stress 
   8 : 
   9 :       IMPLICIT REAL*8 (A-B,D-H,O-Z)
  10 :       implicit integer (I-K,M,N)
  11 :       implicit logical (L)
  12 :       implicit character*10 (C)
  13 : 
  14 :       dimension sigma(6)
  15 : 
  16 :       sigmaxx = sigma(1)
  17 :       sigmayy = sigma(2)
  18 :       sigmazz = sigma(3)
  19 :       sigmaxy = sigma(4) 
  20 :       sigmaxz = sigma(5) 
  21 :       sigmayz = sigma(6)
  22 :        
  23 :       rloc = (sigmaxx-sigmayy)**2 + (sigmayy-sigmazz)**2 +
  24 :      .            (sigmazz-sigmaxx)**2
  25 :       rloc2 = ((1d0)/(6d0))*rloc
  26 :       rj2 = rloc2 + sigmaxy**2 + sigmayz**2 + sigmaxz**2
  27 :       ri1 = sigmaxx + sigmayy + sigmazz
  28 :       DPfunc = dsqrt (3d0*rj2) + alpha*ri1
  29 : 
  30 :       return
  31 :       end
  32 : 
  33 :       
  34 : 
  35 : 
  36 :