Negative Jacobian in some nodes issue

Moderator: GiD Team

Post Reply
zaherkurdi
Posts: 4
Joined: Tue Mar 21, 2023 6:00 pm

Negative Jacobian in some nodes issue

Post 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 2432 times
Attachments
pullouttest.gid.zip
(1.02 MiB) Downloaded 167 times
User avatar
escolano
Posts: 1918
Joined: Sun Sep 05, 1982 10:51 pm

Re: Negative Jacobian in some nodes issue

Post 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 2420 times
I attach your GiD
pullouttest_semistructured.gid.zip
(402.73 KiB) Downloaded 165 times
model meshed semi-structured:
User avatar
abelcs
Posts: 146
Joined: Mon Oct 13, 2014 9:14 am

Re: Negative Jacobian in some nodes issue

Post 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
Attachments
contact_volumes.png
contact_volumes.png (32.53 KiB) Viewed 2411 times
streched_elems.png
streched_elems.png (62.85 KiB) Viewed 2411 times
Post Reply