Dear Kostas,
first, thank a lot for the feedback, the convergence is much cleaner now!
Thanks
also for sharing this script here, it contains several features i had
never use so far, it's quite interesting/intriguing.
However, i am afraid i did not properly explained what my issue was.
Somehow,
i was able to make the solver converge even without enabling some
options for the solve and the result was qualitatively correct, as far
as i could tell.
But it turned out there is a large
difference between the pressure predicted by the Hertz theory and the
result i get and i was wondering whether i am doing something wrong. I
tried to list the potential suspects in the first email
Regarding your remarks:
-Yes,using
a punctual load bother me a lot, which is why i wonder if i am supposed
to do that to reproduce Hertz analytical case.
-Regarding
the unit, i wanted to stay in the small deformation assumption as i was
not sure the analytical solution is supposed to hold for large
deformation. Therefore, here, the Saint Venant Kirchhoff "becomes" a
mere Hooke's law.
Do you have any idea how to reproduce the analytical results?
Thanks
to your help regarding the parameters fiddling, i was at least able to
reproduce a different case where there is no ponctual force (surfacic
load instead) with a half-disk and a different analytical formula
and...It worked (see the results enclosed for the pressure comparison on
the contact boundary)!
Before, it failed to converge.
However, this formula is not really well documented at all, i would
rather have a comparison with the classical formula.
Thanks in advance for your help,
David.
PS: sorry for the duplicate, i forgot to put the mailing list in copy previously