USDFLD subroutine, friction angle , plastic strain

[Majid_F93]

Hi everyone
i wrote a USDFLD subroutine in abaqus and i want to see is it correct or not?

my friction angle and dilation angle varying with plastic strain(PEEQ)


subroutine usdfld(field,statev,pnewdt,direct,t,celent,
1 time,dtime,cmname,orname,nfield,nstatv,noel,npt,layer,
2 kspt,kstep,kinc,ndi,nshr,coord,jmac,jmatyp,matlayo,
3 laccfla)
c
include 'aba_param.inc'
c
character*80 cmname,orname
character*3 flgray(15)
dimension field(nfield),statev(nstatv),direct(3,3),
1 t(3,3),time(2)
dimension array(15),jarray(15),jmac(*),jmatyp(*),
1 coord(*)

c --------------------------------------------------------------------------
c USER CODE START

c it does only for specific material than mention below
if (cmname .eq. 'Material-1') then

c get plastic strain (value PE) and save it into 'array'
call getvrm('PE',array,jarray,flgray,jrcd,jmac,jmatyp,
& matlayo,laccfla)

c PEEQ output = array(7)
Pstrain = abs(array(7))

c assign strain value to field variable
c Friction angle and Dilation angle depend on field variable
field(1) = Pstrain


end if

c USER CODE END
c --------------------------------------------------------------------------

return
end


pictures of material property and a part of .inp file is available in below link:

thanks for helping