Numerical approximation of the solution of Koiter's model for an elliptic membrane shell subjected to an obstacle via the penalty method

Xin Peng,Paolo Piersanti,Xiaoqin Shen

2024-03-12

Abstract:This paper is devoted to the analysis of a numerical scheme based on the Finite Element Method for approximating the solution of Koiter's model for a linearly elastic elliptic membrane shell subjected to remaining confined in a prescribed half-space. First, we show that the solution of the obstacle problem under consideration is uniquely determined and satisfies a set of variational inequalities which are governed by a fourth order elliptic operator, and which are posed over a non-empty, closed, and convex subset of a suitable space. Second, we show that the solution of the obstacle problem under consideration can be approximated by means of the penalty method. Third, we show that the solution of the corresponding penalised problem is more regular up to the boundary. Fourth, we write down the mixed variational formulation corresponding to the penalised problem under consideration, and we show that the solution of the mixed variational formulation is more regular up to the boundary as well. In view of this result concerning the augmentation of the regularity of the solution of the mixed penalised problem, we are able to approximate the solution of the one such problem by means of a Finite Element scheme. Finally, we present numerical experiments corroborating the validity of the mathematical results we obtained.

Numerical Analysis,Analysis of PDEs

What problem does this paper attempt to address?

The paper aims to address the numerical approximation problem of a linearly elastic elliptical membrane shell under specific obstacle conditions. Specifically, the study focuses on the behavior of a linearly elastic elliptical membrane shell constrained within a predetermined half-space in the absence of friction. This problem is formulated as a set of variational inequalities involving a fourth-order elliptic operator.

The main contributions of the paper are as follows:
1. **Uniqueness and Variational Inequalities**: It is proven that the considered problem has a unique solution and satisfies a set of variational inequalities governed by a fourth-order elliptic operator.
2. **Penalty Method Approximation**: It is demonstrated that the solution to the problem can be approximated using the penalty method.
3. **Enhanced Regularity**: It is proven that the solution to the penalty problem exhibits higher regularity at the boundary.
4. **Mixed Variational Form**: A mixed variational form of the penalty problem is provided, and it is proven that its solution also exhibits higher regularity at the boundary.
5. **Finite Element Method**: Based on the above results, a convergent finite element scheme is proposed to approximate the solution.
6. **Numerical Experiments**: The effectiveness of the mathematical results is validated through numerical experiments.

Overall, the goal of the paper is to conduct a precise numerical analysis of the behavior of a linearly elastic elliptical membrane shell under specific obstacle conditions in the absence of friction and to propose effective numerical methods to solve such problems.