Kirchhoff-Love shell implementation in FEniCS

Hi. I am looking for possible examples or demos or tutorials to implement the Kirchhoff-Love theory for thin shells via FEniCS. All I could find was the Naghdi problem of fenics-shell, which is a bit heavier for my current level of understanding.

I am a newbie in FEniCS, variational formulation, UFL etc…So any kind of guidance and code snippets is deeply appreciated.

Thank you in advance :slight_smile:

Hi,
Kirchhoff-Love shells is not the easiest to implement as it requires C^1-continuity of the deflection field and FEniCS does not provide elements ensuring this continuity. A first remedy is to resort to a DG formulation as detailed here Clamped Kirchhoff-Love plate — fenics-shells

A second approach is to relax this continuity and use a Reissner-Mindlin theory instead which would be very close to Kirchhoff-Love in the thin plate limit, provided that shear locking issues have been taken care of. You can find many implementations of that in fenics-shells or also here for plates:
https://comet-fenics.readthedocs.io/en/latest/demo/reissner_mindlin_quads_reduced_integration/reissner_mindlin_quads.py.html
https://comet-fenics.readthedocs.io/en/latest/demo/reissner_mindlin_crouzeix_raviart/reissner_mindlin_CR.html

Finally, if you look for curved shells in general, fenics-shells has indeed a few implementations using an explicit reference mapping. You can also look at this demo which is able to handle non-manifold meshes:
https://comet-fenics.readthedocs.io/en/latest/demo/linear_shell/linear_shell.html

5 Likes

Dear bleyerj

Thanks for the reply. I had actually seen all these links but your answer has put them in the right perspective for me ! will start by trying out these tutorials.

can you please give more insights and examples into this statement?