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,




  • Chris Harrold
    Chris Harrold Member, Administrator, Employee Posts: 117
    First Answer First Comment First Anniversary Ansys Employee

    @Ramdane - can you assist here?

  • Ramdane
    Ramdane Member, Employee Posts: 11
    First Comment First Anniversary Ansys Employee Name Dropper
  • Guillermo
    Guillermo Member, Employee Posts: 1
    Name Dropper First Comment

    Hi @Egbert ,
    it is my pleasure to assist you.

    I tried the code shown above for some tetrahedral mesh examples, with success.
    In order to reproduce your issue, could you send us the .rst files?

    Thank you.

    P.S. I used the following example:

    analysis = examples.download_piston_rod()
    model = dpf.Model(analysis)