Esope source files

   1 : C PRANDT    SOURCE    CHAT      05/01/13    02:21:23     5004           
   2 : c*****************************************************************************
   3 : c       sous rouine utilisee dans l'ecoulement plastique de la loi
   4 : c       de GURSON 
   5 : c*********************************************************************
   6 :         subroutine prandt (wrkgur)
   7 : c       
   8 : c       this subroutine is used whenever phin < phi0 and Pn+1(phin)>0
   9 : c       it performs a classical radial return algorithm for the Prandtl Reuss model
   9 : c       rho and sqrtj2 must have been updated for the trial state before the call 

  10 : c
  11 : c---    variables
  12 :       IMPLICIT INTEGER(I-N)
  13 :         IMPLICIT REAL*8 (A-H,O-Z)
  14 :         segment wrkgur
  15 :          real*8 sigbar, sy0,phi0,rho0,g,b,h
  16 :          real*8 epn,phin,sqrtj2,rho,sig(6)
  17 :          real*8 e(7),dt
  18 :          real*8 conv,tol1,tol2
  19 :         endsegment
  20 : 
  21 : *        common/prop/sigbar,sy0,phi0,rho0,g,b,h
  22 : *        common/state/epn,phin,sqrtj2,rho,sig(6)
  23 : c               e is the strain rate:
  24 : c                       e(1-6) is the deviatoric part
  25 : c                       e(7) is the trace/3
  26 : c               dt ist the time increment
  27 : *        common /delta/ e(7),dt
  28 : 
  29 : c
  30 : c---    equivalent strain increment
  31 :         dep=(sqrtj2-sy0-h*epn)/(h+3.d0*g)
  32 :         epn=epn+dep
  33 :         sqrtj2=sqrtj2-3.d0*g*dep
  34 : c
  35 :         fac=1.d0+3.d0*g*dep/sqrtj2
  36 :         do 100 i=1,6
  37 :         sig(i)=sig(i)/fac
  38 :  100    continue
  39 : c
  40 :         return
  41 :         end
  42 : c*****************************************************************************
  43 : 
  44 : 
  45 :  
  46 :