mapping.on_coordinates operator on tetrahedrons mesh

Hello,
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,

Egbert

Tagged:

Answers

  • Chris Harrold
    Chris Harrold Member, Administrator, Employee Posts: 183
    5 Answers 100 Comments Photogenic Name Dropper
    admin

    @Ramdane - can you assist here?

  • Ramdane
    Ramdane Member, Employee Posts: 13
    10 Comments Name Dropper Ansys Employee First Anniversary
    ✭✭✭
  • Guillermo
    Guillermo Member, Employee Posts: 1
    First Comment Name Dropper
    ✭✭✭

    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)