Multigrid methods for stress-based contact problems in linear elasticity

Staff - Faculty of Informatics

Date: 14 September 2021 / 08:30 - 10:00


You are cordially invited to attend the PhD Dissertation Defence of Gabriele Rovi on Tuesday 14 September 2021 at 08:30 on Zoom (Meeting-ID: 895 8808 9373 code: 078302)

Problems of contact mechanics arise in many engineering applications. The contact conditions make the problem constrained and are the main challenge to be tackled. Different weak forms can be used for the modeling of contact problems. The primal weak form solves for the only displacement in H_1. After the finite element (FE) discretization, the problem can be solved by means of the monotone multigrid (MMG) method, a solver for constrained problems with optimal complexity. However, the primal formulation is affected by locking and can compute the stress, a physical quantity of primary interest, only by means of differentiation of the displacement. In contrast to the primal formulation, the stress-based formulations are not affected by locking, as they are based on the stresses as main unknowns. Thus they are quite attractive for nearly incompressible and incompressible materials. In this thesis, we study the first-order system least squares (FOSLS) and the dual formulations for the Signorini problem, i.e., a unilateral contact problem. In the first approach, an energy functional subject to only box-constraints has to be minimized. In the second case, also global equality constraints must be enforced and an LBB (Ladyzhenskaya-Babuska-Brezzi) condition must be satisfied. This thesis extends the MMG method for the primal formulation to the stress-based formulations applied to the Signorini problem for nearly incompressible and incompressible materials. However, in the stress-based formulations, the stress belongs to the space H_div and therefore special care is needed for the FE discretization and the corresponding solution methods. Linear multigrid methods which work for H_div spaces have been already investigated. To the author's knowledge, this thesis is the first attempt to generalize the MMG from the primal formulation to the stress-based ones. To this purpose, we generalize the Arnold-Falk-Winther patch smoother for H_div-regular problems, so that the contact constraints are solved locally. We show several numerical experiments that illustrate the performance of our new multigrid method for both the FOSLS and the dual cases.

Dissertation Committee:
- Prof. Rolf Krause, Università della Svizzera italiana, Switzerland (Research Advisor)
- Prof. Cesare Alippi, Università della Svizzera italiana, Switzerland (Internal Member)
- Prof. Michael Multerer, Università della Svizzera italiana, Switzerland (Internal Member)
- Prof. Christian Hesch, Universität Siegen. Germany (External Member)
- Prof. Gerhard Starke, University of Duisburg, Germany (External Member)



Festa in corte

Academy of Architecture, Faculty of Communication, Culture and Society, Faculty of Biomedical Sciences, Faculty of Economics, Faculty of Informatics