Esope source files

   1 : C DPLAS3D   SOURCE    AM        12/02/27    23:39:38     7289           
   2 :       SUBROUTINE dplas3D(x1,x2,P2,pi2,ag,fc,fb,
   3 :      .                  parahot3,idimpara3,dplas,lcp) 
   4 : 
   5 : c     This subroutine calculates the INCREMENT OF PLASTIC STRAIN corresponding to
   6 : c      the plastic multipliers Dlambdat = x1 and Dlambdac = x2
   7 : 
   8 :       IMPLICIT REAL*8 (A-B,D-H,O-Z)
   9 :       implicit integer (I-K,M,N)
  10 :       implicit logical (L)
  11 :       implicit character*10 (C)
  12 : 
  13 :       dimension sigmai(6)
  14 :       dimension pi2(6)
  15 :       dimension dgc(6),dgt(6),dplas(6),P2(6,6)
  16 :       dimension parahot3(idimpara3)
  17 : 
  18 :       i0 = 0
  19 :       i1 = 1
  20 :       i2 = 2
  21 :       i3 = 3
  22 :       i4 = 4
  23 :       i5 = 5
  24 :       i6 = 6
  25 :       i7 = 7
  26 : 
  27 :       r2 = 2.
  28 : 
  29 :       Dlambdat = x1
  30 :       Dlambdac = x2
  31 :       
  32 : c     take the effective stress at the previous converged step
  33 :       do iloc=i1,i6
  34 :         sigmai(iloc) = parahot3(idimpara3-48+iloc)
  35 :       end do
  36 : 
  36 :       call dsiggc(sigmai,dgc,ag,P2,pi2,i1,i6,fc,fb,lcp)
  36 :       call dsigRank(sigmai,dgt,i3,i6,parahot3,idimpara3,rkappait,lcp)
  37 :            
  38 :       do iloc=i1,i6
  39 :         dplas(iloc) = Dlambdat*dgt(iloc) + Dlambdac*dgc(iloc)
  40 :       end do
  41 :       
  42 :       RETURN
  43 :       END
  44 :       
  45 : 
  46 : 
  47 : 
  48 :       
  49 :