Geometrical mixed finite element methods for fourth order obstacle problems in linearised elasticity

TL;DR

This paper introduces a new mixed finite element method for fourth order obstacle problems in linearised elasticity, demonstrating convergence and addressing geometric constraints with numerical validation.

Contribution

It develops a novel mixed finite element approach for fourth order variational problems, including obstacle problems and shallow shells, with convergence analysis and implementation insights.

Findings

Convergence of the proposed finite element method for biharmonic obstacle problems.

Implementation challenges for curved shallow shells requiring additional constraints.

Simplified implementation for flat shallow shells without symmetry constraints.

Abstract

This paper is devoted to the study of a novel mixed Finite Element Method for approximating the solutions of fourth order variational problems subjected to a constraint. The first problem we consider consists in establishing the convergence of the error of the numerical approximation of the solution of a biharmonic obstacle problem. The contents of this section are meant to generalise the approach originally proposed by Ciarlet \& Raviart, and then complemented by Ciarlet \& Glowinski. The second problem we consider amounts to studying a two-dimensional variational problem for linearly elastic shallow shells subjected to remaining confined in a prescribed half-space. We first study the case where the parametrisation of the middle surface for the linearly elastic shallow shell under consideration has non-zero curvature, and we observe that the numerical approximation of this model via a mixed Finite Element Method based on conforming elements requires the implementation of the additional constraint according to which the gradient matrix of the dual variable has to be symmetric. However, differently from the biharmonic obstacle problem previously studied, we show that the numerical implementation of this result cannot be implemented by solely resorting to Courant triangles. Finally, we show that if the middle surface of the linearly elastic shallow shell under consideration is flat, the symmetry constraint required for formulating the constrained mixed variational problem announced in the second part of the paper is not required, and the solution can thus be approximated by solely resorting to Courant triangles. The theoretical results we derived are complemented by a series of numerical experiments.