I am working on an automation python script in Abaqus 2020HF9.
I need to get the PEEQ value at nodes for further post-processing. I am using below logic:
 

1. fieldOutput = odb.steps[step_name].frames[-1].fieldOutputs[result_name].getSubset(position=ELEMENT_NODAL, readOnly=True)
2. Simple average the values from different elements to a single value at a node.
    

However, I noticed that, I am getting negative values after averaging. After investigating, I came to know that values after getting subset ELEMENT_NODAL are negative - which is strange.
To confirm, I output the values ELEMENT_NODAL using Abaqus-Viewer using Report -> Field Output. And in the report file, there were no negative values.

I do not understand the reason for this deviation. I have also confirmed the same in Abaqus 2024HF7 (though with the same ODB).