How to call PEEQ and stress related parameters using VUSDFLD Subroutine?

[Hossein Mamusi]

I want to model some ductile fracture criteria in ABAQUS/Explicit using VUSDFLD subroutine. For this I need to use some stress (PEEQ) and strain related parameters. I wrote a simple subroutine but it can't to recall PEEQ. Below is a part of this subrotine.

jStatus=1
call vgetvrm('SP', rData, jData, cData, jStatus)
stress = reshape(rData, (/nblock, nrData/))
jStatus=1
call vgetvrm('PEEQ', eqps, jData, cData, jStatus)
call vgetvrm( 'S', rData, jData, cData, jStatus )
jStatus=1
s11 = stress(k,1)
s22 = stress(k,2)
s33 = stress(k,3)
s12 = stress(k,4)
s13 = stress(k,5)
s23 = stress(k,6)
do i=1, nblock
pq = eqps(k)
vonMises1 = (s11-s22)**2+(s22-s33)**2+(s11-s33)**2
vonMises2 = 6*(s12**2+s23**2+s13**2)
vonMises = sqrt(0.5*(vonMises1+vonMises2))
hydrostatic = (s11+s22+s33)/3
lode=((2*s22-s11-s33)/(s11-s33))
triax = hydrostatic/vonMises
stateNew(i,1) = vonMises
stateNew(i,2) = hydrostatic
stateNew(i,3) = triax
stateNew(i,4) = lode
end do

I'm looking forward your help.
Thanks