Numerical Approximation of the Solution of an Obstacle Problem Modelling the Displacement of Elliptic Membrane Shells via the Penalty Method

In this paper we establish the convergence of a numerical scheme based, on the Finite Element Method, for a time-independent problem modelling the deformation of a linearly elastic elliptic membrane shell subjected to remaining confined in a half space. Instead of approximating the original variational inequalities governing this obstacle problem, we approximate the penalized version of the problem under consideration. A suitable coupling between the penalty parameter and the mesh size will then lead us to establish the convergence of the solution of the discrete penalized problem to the solution of the original variational inequalities. We also establish the convergence of the Brezis–Sibony scheme for the problem under consideration. Thanks to this iterative method, we can approximate the solution of the discrete penalized problem without having to resort to nonlinear optimization tools. Finally, we present numerical simulations validating our new theoretical results.