Second elasticity tensor does not show minor symmetries

[dokken]

As recently discussed in Issue with dolfinx.fem.Function.interpolate - #2 by dokken, projecting discontinuous data in to continuous spaces can introduce artifacts:

Issue with dolfinx.fem.Function.interpolate

You are interpolating a discontinuous function, i.e. into a continuous space. This will cause some interpolation error which can be manifested in non-symmetry. Talking about smoothness, as the bending moment only depends on theta, which is a piecewise linear, the best you will get is a piecewise constant function. Here using paraview to illustrate this: [image] Upon mesh refinement, the stresses gets better, in the sense that the jumps between each element becomes smaller: (40x40) [im…