How to restrict solutions or stresses to subdomains?

[dokken]

@steve

Your example, dokker, worked perfectly well as long as the bodies didn’t have different material properties

Note that you can use a DG0 function for the Lame parameters, such that you only need one definition og sigma, see: How to define different materials / importet 3D geometry (gmsh) - #2 by dokken

@steve

By the way, the element order of the FunctionSpace for von Mises stresses are reduced by one compared to that of the displacments, because of the derivative in sigma and epsilon. Did I get that right?

That is correct.