Contact Problem

Hi there,

As a initial step I am trying to solve the linear elastic problem and my geometry is coming in contact with its own surfaces upon loading. This contact is not correct as the surfaces are penetrating into each other.

I am trying to over come this by defining the contact by penalty formulation which is provided in this explanation.

https://comet-fenics.readthedocs.io/en/latest/demo/contact/penalty.html

But I am finding it hard to understand as this problem of mine is a generic geometry and I do not know where the contacts happen.

If anyone can help me with this I would be thankful.

Attached is my output image.

Very Similar question to How to implement the penalty method to avoid self-contact - #2 by dokken