Numerical approximation of the solution of an obstacle problem modelling the displacement of elliptic membrane shells via the penalty method
Abstract
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.
- Publication:
-
arXiv e-prints
- Pub Date:
- May 2022
- DOI:
- 10.48550/arXiv.2205.11293
- arXiv:
- arXiv:2205.11293
- Bibcode:
- 2022arXiv220511293M
- Keywords:
-
- Mathematics - Analysis of PDEs
- E-Print:
- The paper went through a thorough restyling, as an easier way of proving the regularity result was found. We also added new proofs on the convergence of a finite element scheme and numerical simulations. arXiv admin note: text overlap with arXiv:2107.07673