Can't managed to modify mesh from k file

[Thmbzh]

Hi everyone,
I have some issues to modify the mesh of my .k file.
I succeed into modifying the tuto one and changing parameters, but when using my file, nothing happens.
I tried to change every parameters with different values but nothing happens with the mesh.
I might going wrong with the definition of my parameters but hard to know.
I run a 2D axysym mesh.

Thanks in advanced
Stay safe at home !

[mtezzele]

Dear @Thmbzh, are you able to read the coordinates of your mesh? I mean, maybe the problem is not in the deformation part, but in the parsing part.

[Thmbzh]

Thank you for your kick answer.
I'm sorry, I'm not familiar with this kind of tools, how can I see if I can read the coordinates ?
I don't think the problem is from the script itself, I used the tuto one. So it is about the mesh indeed.

[mtezzele]

The parse method returns the mesh points. Is this variable empty or not?

[Thmbzh]

This is the script I use

from pygem import FFDParameters, FFD, KHandler, VtkHandler
import vtkplotter as vtkp

"params = FFDParameters()"
"params.read_parameters(filename='/home/thomas/PyGeM/tests/test_datasets/parameters_test_ffd_pipe_unv_C0.prm')"

"%cat ../tests/test_datasets/parameters_test_ffd_pipe_unv_C0.prm"
k_handler = KHandler()
mesh_points = k_handler.parse('/home/thomas/PyGeM/TEST_PIECE/exemple_1.k')

"free_form = FFD(params, mesh_points)"
"free_form.perform()"
"new_mesh_points = free_form.modified_mesh_points"
"k_handler.write(new_mesh_points, 'test_pipe_mod_C0.k')"

"%cat ../tests/test_datasets/parameters_test_ffd_pipe_unv_C1.prm"

params = FFDParameters()
params.read_parameters(filename='/home/thomas/PyGeM/TEST_PIECE/parameters_test.prm')

free_form = FFD(params, mesh_points)
free_form.perform()
new_mesh_points = free_form.modified_mesh_points
k_handler.write(new_mesh_points, 'test_C1.k')

[mtezzele]

Just check the content of the mesh_points variable, that is the result of the parse method. If it is empty then you have to modify the parse method in order to properly read your file. We tested only on one file (the one present in the tests folder). It was an external contribution so we do not have a direct experience with k files.

[Thmbzh]

I somehow succeed to have a deformation. It seems it didn't like to have z = 0 for 2D elements, worked with z = -1.0
I'm truly sorry if I didn't tried enough variables, didn't want to bother you with such a little thing.
Thanks for your support, I'll try to give you any feedback with this type of mesh !