Consult a problem about the VUSDFLD_V1.f

[rookie-coder99]

Dear Morin,
Thank you for sharing the source code of the subroutine, I came across a question I would like to ask you, you are using vgetvrm to get the internal variables on the VUSDFLD subroutine, but I found that this array is fetched at nblock intervals, I would like to verify that this point is accurate, but I didn't find anything specific in the manual.

So I added the field variable assignment under your code to determine these values and the value of stress.

      if(nshr.gt.1)then
         do i=1,nblock
            s(i,1) = stressdata(i)
            s(i,2) = stressdata(i+nblock)
            s(i,3) = stressdata(i+nblock*2)
            s(i,4) = stressdata(i+nblock*3)
            s(i,5) = stressdata(i+nblock*4)
            s(i,6) = stressdata(i+nblock*5)
	    STATENEW(i,1) = s(i,1)
	    STATENEW(i,2) = s(i,2)
	    STATENEW(i,3) = s(i,3)
	    STATENEW(i,4) = s(i,4)
	    STATENEW(i,5) = s(i,5)
	    STATENEW(i,6) = s(i,6)
         enddo
      else
         do i=1,nblock
            s(i,1) = stressdata(i)
            s(i,2) = stressdata(i+nblock)
            s(i,3) = stressdata(i+nblock*2)
            s(i,4) = stressdata(i+nblock*3)
	    STATENEW(i,1) = s(i,1)
	    STATENEW(i,2) = s(i,2)
	    STATENEW(i,3) = s(i,3)
	    STATENEW(i,4) = s(i,4)
         enddo
      endif

But it turns out that the values of stress and these state variables are not equal in the calculation result, I am not sure if there is a problem here, so I am asking you for advice!

[rookie-coder99]

Here is some results.