Numerical simulation of a fine-tunable Föppl-von Kármán model for foldable and bilayer plates

A numerical scheme is proposed to identify low energy configurations of a F\"oppl-von K\'arm\'an model for bilayer plates. The dependency of the corresponding elastic energy on the in-plane displacement $u$ and the out-of-plane deflection $w$ leads to a practical minimization of the functional via a decoupled gradient flow. In particular, the energies of the resulting iterates are shown to be monotonically decreasing. The discretization of the model relies on $P1$ finite elements for the horizontal part $u$ and utilizes the discrete Kirchhoff triangle for the vertical component $w$. The model allows for analysing various different problem settings via numerical simulation: (i) stable low-energy configurations are detected dependent on a specified prestrain described by elastic material properties, (ii) curvature inversions of spherical and cylindrical configurations are investigated, (iii) elastic responses of foldable cardboards for different spontaneous curvatures and crease geometries are compared.