Importing Tetrahedron mesh from Abaqus into dealii

[Masoud Ahmadi]

Hi all,

I import Hexahedron mesh from Abaqus *inp file into dealii by the following lines of code:

Triangulation<dim> mesh;

GridIn<dim> grid_in;

grid_in.attach_triangulation(mesh);

std::ifstream input_mesh("cube.inp");

grid_in.read_abaqus(input_mesh);

and it works fine:

But, with the same lines of code, for importing Tetrahedron mesh I get:

which obviously is not correct. What should I do here?

Best regards,

Masoud

[Peter Munch]

Hi Masoud,

it seems that read_abaqus() calls read_ucd() and the the latter has been not generalized for simplex meshes yet. See: https://github.com/dealii/dealii/blob/ee0a4995de1009875a4a5572ca14bde8a1979ad5/source/grid/grid_in.cc#L948-L950. I am a bit surprised that you did not get an error message. Would you be willing to generalize the function? You could take a look at read_vtk() or read_msh(), which have been generalized for arbitrary mesh types. Alternatively, you could use one of the latter two function.

Hope this helps,
Peter