Dear Yves
Thanks a lot for your advice.
> It seems that on some distorted quadrilateral elements, the inversion of the geometric transformation converges well, but outside the element (you can trace the inversion in the function geotrans_inv_convex::invert_nonlin in bgeot_geotrans_inv.cc to see it) which means that there are possibly several solutions to the inversion.
Thanks. I will check it's part. I will ask again if I could not trace it.
> For the moment, I do not see how to solve that. Appart to shift a bit the transformation such as
> md.add_interpolate_transformation_from_expression("Proj1", mesh1, mesh2, "[X(1);-X(2)-1e-2]")
I will try this.
Best regards Tetsuo
2021年2月16日(火) 17:51 Yves Renard <Yves.Renard@insa-lyon.fr>:
Dear Tetsuo,
It seems that on some distorded quadrilateral elements, the inversion of the geometric transformation converges well, but outside the element (you can trace the inversion in the function geotrans_inv_convex::invert_nonlin in bgeot_geotrans_inv.cc to see it) which means that there are possibly several solutions to the inversion. For the moment, I do not see how to solve that. Appart to shift a bit the tranformation such as
md.add_interpolate_transformation_from_expression("Proj1", mesh1, mesh2, "[X(1);-X(2)-1e-2]")
Best regards,
Yves
On 13/02/2021 11:04, Tetsuo Koyama wrote:
Dear GetFEM community
I want to solve the contact problem between two cylinders. The mesh is like the following.This is a 1/4 model with symmetry conditions on the left and having Dirichlet condition [0.-0.1] on the top and [0.0, 0.0] on the bottom.
I used this tutorial `Example of wheel in contact (Assembly between two meshes, transformations, use of fixed size variables)` as a guide to create a program like the attached script "hertz.py", but an error occurred:
RuntimeError: (Getfem::InterfaceError) -- Error in getfem_generic_assembly_compile_and_exec.cc, line 1433 virtual int getfem::ga_instruction_interpolate_base::exec():
No valid element for the transformation. Probably transformation failed
I think my projection is correct.md.add_interpolate_transformation_from_expression("Proj1", mesh1, mesh2, "[X(1);-X(2)]")
Do you have any idea how to solve this problem?Thanks a lot for reading.
Best regards Tetsuo
-- Yves Renard (Yves.Renard@insa-lyon.fr) tel : (33) 04.72.43.87.08 INSA-Lyon 20, rue Albert Einstein 69621 Villeurbanne Cedex, FRANCE http://math.univ-lyon1.fr/~renard ---------
