IncompressibleLinearElasticity from tet-vs-hex not working

[ateixeira0163]

Hi! Great work!
As the title said, I'm wanting to simulate and incompressible material and using NeoHookean with a high poisson ratio (0.48) seems to yield bad results. I've then seen that there's an option as IncompressibleLinearElasticity that was even used in the paper comparing tet and hex meshes, but the example showed doesn't work anymore and just changing from NeoHookean to IncompressibleLinearElasticity with nu = 0.5 yields a error with IPC : initial distance is 0, returning toi=0! until the whole thing crashes.

I imagine this is due to the updates made to the code, but, is there an example using IncompressibleLinearElasticity that I could see/use?
Thank you very much!

[danielepanozzo]

Afaik we do not have any IPC example that uses the IncompressibleLinearElasticity. We have not used that feature in a while. A possible workaround, which is easy to try, is to use higher-order basis. Quadratic or cubic basis behave much better than linear ones for a high poisson ratio.

[ateixeira0163]

This is done by the ['space']['discr_order'] parameter?

[teseoch]

Yes

[teseoch]

If you are still struggling maybe you can share the whole log in debug or trace