If you instead of using subdomains, use a piecewise constant function to define G, it is easy to obtain the stresses. See for instance How to define different materials / importet 3D geometry (gmsh).
If you instead of using subdomains, use a piecewise constant function to define G, it is easy to obtain the stresses. See for instance How to define different materials / importet 3D geometry (gmsh).