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 <[email protected]>:
>
> 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
>
> [image: mesh.png]
>
>
> --
>
> Yves Renard ([email protected]) tel : (33) 04.72.43.87.08
> INSA-Lyon
> 20, rue Albert Einstein
> 69621 Villeurbanne Cedex, FRANCE
> http://math.univ-lyon1.fr/~renard
>
> ---------
>
>
