Esope source files

   1 : C GURSON    SOURCE    CHAT      05/01/13    00:22:25     5004           
   2 : c*****************************************************************************
   3 : c       sous rouine utilisee dans l'ecoulement plastique de la loi
   4 : c       de gurson
   5 : c**********************************************************************
   6 :         subroutine gurson (s,p,depsp,gamma,poro,densit,epsi,dtime,
   7 :      &  wrkgur)
   8 : c
   9 : c       this subroutine is the interface between the finite element program
  10 : c       and the rest of the material update program
  11 : c       given the history variables and the starin rate it updates the
  12 : c       state
  13 : c               by arnaud de gayffier 06/13/93 in san diego california
  14 : c
  15 : c---    arguments
  16 : c               s........deviatoric stress  in gpa
  17 : c                               s(1)=sxx|
  18 : c                               s(2)=syy| s(1)+s(2)+s(3)=0
  19 : c                               s(3)=szz|
  20 : c                               s(4)=sxy
  21 : c                               s(5)=sxz
  22 : c                               s(6)=syz
  23 : c               p........pressure (if >0 compression)  in gpa
  24 : c               gamma....cumulative plastic strain
  25 : c               poro.....porosity or void fraction
  26 : c               densit...material density  in gr/cm3
  27 : c               epsi.....total strain rate same conventions as "s"
  28 : c               depsp....total plasti strain rate same convention
  29 : c               dtime....time increment
  30 : c
  31 : c         warning!! the argument epsi might be modified
  32 : c
  33 : c---    variables
  34 : c
  35 :       IMPLICIT INTEGER(I-N)
  36 :         IMPLICIT REAL*8 (A-H,O-Z)
  37 :         segment wrkgur
  38 :          real*8 sigbar, sy0,phi0,rho0,g,b,h
  39 :          real*8 epn,phin,sqrtj2,rho,sig(6)
  40 :          real*8 e(7),dt
  41 :          real*8 conv,tol1,tol2
  42 :         endsegment
  43 :           
  44 : *        COMMON/PROP/SIGBAR,SY0,PHI0,RHO0,G,B,H
  45 : *        COMMON/STATE/EPN,PHIN,SQRTJ2,RHO,SIG(6)
  46 : c               e is the strain rate:
  47 : c                       e(1-6) is the deviatoric part
  48 : c                       e(7) is the trace/3
  49 : c               dt ist the time increment
  50 : *        COMMON /DELTA/ E(7),DT
  51 :         dimension s(6),epsi(6),depsp(6)
  52 : c
  53 : c---    identification
  54 :         DPHI=0.D0
  55 :         SQRTJ2=0.D0
  56 :         epn=gamma
  57 :         phin=poro
  58 :         rho=densit
  59 :         do 100 i=1,6
  60 :           e(i)=epsi(i)
  61 :           SIG(I)=S(I)
  62 :   100   CONTINUE 
  63 :         dt=dtime
  64 : c
  65 : c---    partition in deviatoric and hydrostatic
  66 :         E(7)=(E(1)+E(2)+E(3))/3.D0
  67 :         E(1)=E(1)-E(7)
  68 :         E(2)=E(2)-E(7)
  69 :         E(3)=E(3)-E(7)
  70 :         E(4)=E(4)/2.D0
  71 :         E(5)=E(5)/2.D0
  72 :         E(6)=E(6)/2.D0
  73 : c
  74 :   101   FORMAT (A,7(1PE12.5,1X))
  75 :   102   FORMAT (5(A,1PE12.5))
  76 : c
  77 : c---    material update
  77 :         call newsta(wrkgur)
  78 : c
  79 : c---    arguments update
  80 :         DPHI = PHIN-PORO
  81 :         DEPN = EPN-GAMMA
  82 : c
  83 :         P=(1.D0-PHIN)*B*(1.D0-(1.D0-PHIN)*RHO0/((1.D0-PHI0)*RHO))
  84 :         gamma=epn
  85 :         poro=phin
  86 :         densit=rho
  87 :         do 200 i=1,6
  88 :   200   S(I)=SIG(I)
  89 : c                plastic strain
  90 :         DO 300 I=1,6
  91 :         IF (SQRTJ2 .GE. 1.D-6 ) THEN
  92 :           DEPSP(I)=1.5D0*SIG(I)/SQRTJ2*DEPN
  93 :         ELSE
  94 :           DEPSP(I) =0.D0
  95 :         ENDIF
  96 :   300   CONTINUE
  97 : c
  98 :         DEPSP(4)=2.D0*DEPSP(4)
  99 :         DEPSP(5)=2.D0*DEPSP(5)
 100 :         DEPSP(6)=2.D0*DEPSP(6)
 101 : c
 102 :         DO 400 I=1,3
 103 :          DEPSP(I)=DEPSP(I)+DPHI/(1.D0-PHIN)/3.D0
 104 :   400   CONTINUE
 105 : c
 106 :         RETURN
 107 :         end     
 108 : c------------------------------------------------------------------------------
 109 : 
 110 : 
 111 : 
 112 :  
 113 : 
 114 : 
 115 : 
 116 :  
 117 :