mapping.on_coordinates operator on tetrahedrons mesh
I have a question regarding the "mapping.on_coordinates" operator. To test my code I use a simple ANSYS model consisting of a beam, meshed with volume elements, fixed on one side and loaded with a lateral force on the other side.
My aim is to determine the stresses along the outside-surface node normals. I have the following piece of code to read in the model and plot the mesh.
and to determine the mapped stresses on the normals:
When I use a hexahedron mesh (sweep mesh) in the ANSYS model, the mapping routine works fine, and I can extract the mapped normal stresses for each node:
However when I re-mesh the ANSYS model with tetrahedron elements:
the same script fails:
I played around with the "use_quadratic_elements" option in the operator but did not find a solution for this behaviour. Any help is appreciated.
With kind regards,