Showing the integration point values is an interesting feature.
Yet, I guess that the integration point values displayed as described above aren’t the true IP values, because they are back-interpolated from averaged nodal values if you don’t suppress the averaging.
You will most likely have different values in ParaView and in the dat file. In particular, in ParaView you still can have IPs with the equivalent stress exceeding the yield stress, whereas in the dat file this should not happen.
To avoid this, you need to output the nodal values for each element separately. How to do this has been discussed in other threads here. I made a python script for this purpose.