See for instance: How to restrict solutions or stresses to subdomains? - #5 by steve
which project the Von mises stresses into an appropriate function space, and then plots it with matplotlib.
For 3D geometries, I suggest saving the stresses to XDMF file with write_checkpoint and use Paraview to visualize them.