You are projecting your stress into a CG 1 space, when the stressed are contained in a DG0 space. You should project them into an appropriate function space and use XDMFFile.write_checkpoint to write it to a suitable file format for visualization in paraview. See for instance: When using CG degree 1, the result corresponds to degree 2 - #10 by dokken
2 Likes