Issue converting INP to VTK file format

[prince2110]

I have attached the INP file for your reference. I had to compress it to zip as I was not allowed to upload INP file in this issue form.

The INP file consists of nodes and only volume elements along with element sets and solid sections. I am aiming to convert it to a file format such that I can view element sets in Paraview.

I am trying to convert this INP file to VTK format, but it throws me 'IndexError'. Following is the error details:

$ meshio-convert mesh_TET.inp mesh_TET.vtk

Number of points: 94
Number of cells:
tetra: 221
Traceback (most recent call last):
File "/home/prince/.virtualenvs/python101/bin/meshio-convert", line 8, in
sys.exit(convert())
File "/home/prince/.virtualenvs/python101/lib/python3.5/site-packages/meshio/_cli.py", line 50, in convert
write(args.outfile, mesh, file_format=args.output_format)
File "/home/prince/.virtualenvs/python101/lib/python3.5/site-packages/meshio/_helpers.py", line 125, in write
if value.shape[1] != num_nodes_per_cell[key]:
IndexError: tuple index out of range

mesh_file.zip

[nschloe]

It works in the latest revision of meshio; just upgrade.

[prince2110]

Thanks work your quick reply. It worked but VTK file has only surface information. Can you suggest me to which format I can convert this INP file such that I can view volume elements along with element sets in Paraview ?

[nschloe]

but VTK file has only surface information

Not true, all tets are contained:

[prince2110]

Sorry Sir, for my lack of information. Just a last question, I had created element sets in INP file, but now in VTK how can I view those group of elements again ?

[nschloe]

Aha. Well, VTK doesn't support element sets as abaqus does. I suppose one could map those to integer cell functions. (It's really the reverse of what #716 is about.) I'll take a look.

[nschloe]

Alright, -s should help you out now.

[prince2110]

Respected Sir, I am still struggling to figure out, how to display element sets using VTK file in paraview. Regards

…

On Thu, 20 Feb 2020 at 17:15, Nico Schlömer ***@***.***> wrote: Alright, -s should help you out now. — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#719?email_source=notifications&email_token=AL5K3EPTFZTWGBCGZIS5V3TRD2UAPA5CNFSM4KYN7OG2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEMO6E5Q#issuecomment-589161078>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AL5K3EO3STGT72BHS65CS7LRD2UAPANCNFSM4KYN7OGQ> .

[nschloe]

Respected Sir,

We call each other by our first names here. I'm Nico.

I am still struggling to figure out, how to display element sets using VTK file in paraview.

It should all work now. What have you tried?

[prince2110]

I imported VTK file to paraview. Then created group dataset and tried if vtkBlockColor displays element set. After my failure, I tried to do Color Map. But still I am unable to figure out the way to display set of elements with different colours as it is in Abaqus. Regards

…

On Thu 20 Feb, 2020, 11:18 PM Nico Schlömer, ***@***.***> wrote: Respected Sir, We call each other by our first names here. I'm Nico. I am still struggling to figure out, how to display element sets using VTK file in paraview. It should all work now. What have you tried? — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#719?email_source=notifications&email_token=AL5K3ENICTUMW6HM4KENFBDRD36TNA5CNFSM4KYN7OG2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEMQTSJY#issuecomment-589379879>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AL5K3EIYFKBU2YNUYD7SOM3RD36TNANCNFSM4KYN7OGQ> .

[nschloe]

It works for me.

meshio-convert mesh_TET.inp out.vtk -s

and then open the vtk file.

[keurfonluu]

Looks like a conversion issue. The scale goes completely through the roof in the VTK one.

[nschloe]

Works fine here?

[prince2110]

I do not know if I am the most dumb person on this planet and I think I am. I convert INP to VTK using

$ meshio-convert mesh_TET.inp out.vtk

I have attached the VTK file as well.

When I open this directly on Paraview, I can see only this:

I don't understand if you guys are using some filters or some settings related to Color map in paraview.

Meshio Version

$ meshio-info --version
meshio 4.0.2 [Python 3.5.3]
Copyright (c) 2015-2020 meshio developers

Paraview version
Paraview 5.7.0

VTK.zip

[keurfonluu]

I do not know if I am the most dumb person on this planet and I think I am. I convert INP to VTK using

$ meshio-convert mesh_TET.inp out.vtk

Your command is missing the option -s (see the answers above).

[keurfonluu]

Note that you probably need to install the GitHub version of meshio:

pip install https://github.com/nschloe/meshio/archive/master.zip

[prince2110]

Finally figured out this little secret which I was missing out !! :)

And then I get this:

Thanks a lot guys for helping me !!

[kenza1994]

Hello,

I'm on a project where I need to do what you did in your example.

So, I took a cylinder that I cut into two mesh zones (zone1 = material1 and zone2 = material2). I output an .inp file

After that, I converted it to .vtk file with meshio by adding the -s command

However, when I open my .vtk file with paraview, nothing is displayed.

So, I tried to take yours, it opens and I can visualize the meshes (I can't see the colors of the different materials).

That why, I don't understand why on my side I can't see my files on paraview or see the colors.

How do you do it and which versions do you use?

Is there, something to modify in this .inp file to see the two distinct materials as well as the geometry and the mesh on paraview?

My file .inp is attached.

job.zip