Is there a systematic manner in dolfinx in which I can solve a linear elasticity problem where I assign the effective stress and tangent stiffness at the individual quadrature points? The idea is to solve an auxilliary boundary value problem (given a set of boundary conditions) for each quadrature point and use the homogenized stress and corresponding stiffness for the macroscale problem. This approach is also known as “finite element squared”.

You can use the `Quadrature`

-element, as used in dolfiny: dolfiny/solid_plasticity_monolithic.py at b23bdeb5f202e5c240a20e23b93daed0b2294d61 · michalhabera/dolfiny · GitHub