Set initial stresses from MeshNodeParameter or MeshElementParameter

[joergbuchwald]

I was not aware that MeshElements work. But sure, there the mapping to GP is straightforward.
To do this in Python you can use VTUinterface as it is done in this example:

import vtuIO
f = vtuIO.VTUIO("Decovalex_D_3D.vtu")
pts = f.cell_center_points
# voronoi is used for n.n. interpoaltion, if it is not specified interpolation is linear
sigma = f.get_set_data("sigma", pointsetarray=pts, data_type="point", interpolation_method="voronoi") 
f.add_cell_field(sigma,"sigma","Decovalex_D_3D.vtu") 

this is similar as done here:

Output medium properties (permeability, porosity, ... )

I don’t see, unfortunately, an easy/small code modification to output material properties. Best would be not to mix a Constant with MeshNode properties for, say permeability, but go with MeshNode property for both materials. Performance-wise the difference is negligible.