CONVERGENT FEM FOR A MEMBRANE MODEL OF LIQUID CRYSTAL POLYMER NETWORKS\ast

We design a finite element method for a membrane model of liquid crystal polymer networks. This model consists of a minimization problem of a nonconvex stretching energy. We discuss properties of this energy functional such as lack of weak lower semicontinuity. We devise a discretization with regularization, propose a novel iterative scheme to solve the nonconvex discrete minimization problem, and prove stability of the scheme and convergence of discrete minimizers. We present numerical simulations to illustrate convergence properties of our algorithm and features of the model.