@nate ,
Thank you for the response. I believe that I’ve overlooked the element family I chose for the previous stress state \sigma{}_{n}. Thanks to @bleyerj 's page, specifically Elasto-plastic analysis of a 2D von Mises material and the post Quadrature function space - I realized that for my problem declaring:
is the issue and instead:
Vs_CG_3 = fem.VectorFunctionSpace(domain, (“DG”, 0), dim = 3)
is sufficient given my choice of function space order for the displacement u. I also am using fem.function.interpolate() for going between two different function spaces rather than project() such as from Dolfinx Discontinuous Expressions SIGSEGV.
Thanks,
Sven