Finite element approximation of large bending isometries

A finite element scheme for the approximation of large isometric deformations with minimal bending energy is devised and analyzed. The convergence to a stationary point and energy decreasing property of an iterative algorithm for the numerical solution of the scheme is proved. Numerical experiments illustrate the performance of the iteration and show that the discretization leads to accurate approximations for large vertical loads and compressive boundary conditions.