Page 1 of 1

Negative Jacobian in some nodes issue

Posted: Fri Mar 24, 2023 2:58 am
by zaherkurdi
Dear,

I am reaching out to you regarding a problem I encountered while generating mesh for my model. The model is about a cube of concrete with a steel bar and an interface between them. I received an error message saying "There are elements with Negative Jacobian in some nodes issue".

I have previously created the same model without encountering this issue, and everything was fine. Could you please suggest some solutions to solve this problem? Additionally, I am curious if there are alternative ways to mesh the model, such as using a brick type mesh for the cube. What do you think?

For your reference, the problem type used is Atena, and I am using GID v 16.1.3d. I have attached the model in a ZIP format for your convenience.

Thank you for your time and assistance. I look forward to your response.
Screenshot 2023-03-24 021640.jpg
Screenshot 2023-03-24 021640.jpg (96.02 KiB) Viewed 2609 times

Re: Negative Jacobian in some nodes issue

Posted: Tue Mar 28, 2023 3:29 pm
by escolano
It is raising the message 'There are elements with negative Jacobian in some node'
because of the flat prism elements of the contact volumes.
It seems that the mesh of tetrahedra doesn't have any negative Jacobian.

You can see it with the tool: Mesh->Mesh quality
and set Quality criteria=Minimum Jacobian

the tab of Prism show Minimum Jacobian range from -5.39e-24 to 1.52e-24
this really is numerical noise, can be considered as 0 (this is normal, because the contact elements have volume zero)

About meshing it with cubes, in particular this case is topologically a prismatic shape (an extrusion), and it is possible to generate a semi-structured mesh of hexahedra: Mesh_>SemiStructured->Volumes (and a number of divisions in the extrusion direction),
by default the mesh will be of tetrahedra, can set
Mesh->Elemtype->Hexahedra
semistructured_mesh_hexas.png
semistructured_mesh_hexas.png (39.13 KiB) Viewed 2597 times
I attach your GiD
pullouttest_semistructured.gid.zip
(402.73 KiB) Downloaded 173 times
model meshed semi-structured:

Re: Negative Jacobian in some nodes issue

Posted: Wed Mar 29, 2023 11:31 am
by abelcs
Hello,
It seems there are some tetrahedra with jacobian negative (or almost zero). You can detect them using the Mesh-Mesh Quality tool.
They are in volume 2, and it seems that is due to streched elements in the contour surfaces. Checking it, if you assign a finner mesh to the surfaces (the long ones), the final mesh is ok.

Note that as these surfaces are part of contact volumes (the volumes number from 2 to 63), you should assign the size to the surfaces of that contact volumes. Assigning a size of 0.04 to these surfaces, for instance, the final mesh is ok.
Regards,
ABEL