I think you need to identify the nodes on your region of interest and you can use this code to extract data
session.xyDataListFromField(odb=odb, outputPosition=NODAL, variable=(('S', INTEGRATION_POINT, ((INVARIANT, 'Max. Principal'), )), ), nodeLabels=(('Part-1-1', ('5', '10', '11', '21', )), ))...