Can this be done in NGsolve?
Hi,
yes, you can. The space is called HDivDivSurface in NGSolve.
An example for Koiter shells is here:
https://docu.ngsolve.org/latest/i-tutorials/unit-6.1.4-shells/shell.html#
Best,
Joachim
Hi,
if you want to solve HHJ on surfaces without the other shell stuff as membrane and shearing terms, you can find a simple example attached
https://ngsolve.org/media/kunena/attachments/889/hhj_surface.py
Best
Michael
Attachment: hhj_surface.py
Great, thanks. And this can be adapted for curved surface triangulations, using Lagrange polynomials for the triangles?
Yes, correct, geometry is approximated by arbitrary order polynomially curved trigs or quads, and the HHJ spaces is mapped to the curved surface using double-sided Piola transform. All nasty derivatives are computed.