On Fri, May 6, 2022 at 8:45 AM Konstantinos Poulios <logari81@googlemail.com> wrote:
Hi,Your script looks very nice and clean. It just needs a bit of fiddling with the augmentation parameter and the solver. The attached version produces this resultJust some general comments:- Point constraints/loads make almost never sense (unless you know beforehand the reaction force at a constrained point will be zero)- Maybe there is some issue with your units, I needed to increase the load to a very large value to get some visible deformations- The Saint-Venant Kirchhoff material law, that we often use for demonstrations, is a "toy" material law. If you do anything with moderate/large strains, use some proper hyperelastic law- The augmentation parameter should be in the same order of magnitude like the modulus of elasticity, from 1-2 orders of magnitude below to 1-2 orders of magnitude above.Best regardsKostasOn Thu, May 5, 2022 at 8:04 PM David Danan <daviddanan9021@gmail.com> wrote:Dear Getfem-users,My issue might be utterly simple, but i cannot find what i am doing wrong...I am trying to compare the solution provided by Getfem to the analytical solution of an Hertzian contact in 2D (namely, a deformable ball against a perfectly rigid plane) but, so far, i do not get correct results.Please, find enclosed my script. I use gmsh python api to generate the mesh but, if you do not want to use it, you can comment the line to generate the mesh in the script and juste use the mesh also enclosed.Here are some of my concerns:-Are the formula used in the script for the analytical contact pressure/contact surface area etc... correct in 2D?-Is it correct to use a punctual nodal force in such a case? Does it make sense in this context?-Is it correct to compare the analytical pressure to the lagrange multipliers?I have also tried to use a testcase provided in Code_Aster documentation and i got much better results, although not perfect but within the tolerance given in the documentation below.
https://www.code-aster.org/doc/default/en/man_v/v1/v1.01.246.pdfIt might be because these formula are for 3D only but i am not sure yet.Thanks in advance for your help,kind regards,David.
