Résolution du contact frottant déformable-déformable par lagrangien augmenté et GCR préconditionné

Ce document traite le problème de l’élasticité en grandes déformations avec contact frottant. Dans le premier chapitre, nous présentons un rappel des lois auxquelles nous ferons référence : lois de comportement matériau, conditions aux limites, contact et frottement selon la loi de Coulomb. Le coeur de la thèse commence par la discrétisation des conditions de bord. L’imposition en moyenne de la non-pénétration et de la loi de Coulomb est explicitée dans le but de se ramener en dimension finie. C’est dans ce domaine que les algorithmes sont développés. Les conditions d’optimalité sont donc calculées pour le problème déjà discrétisé et sont linéarisées par la suite selon le paradigme des contraintes actives. Suivant cette stratégie, l’état de contact est fixé à chaque itération de la méthode de Newton, transformant les conditions de contact en contraintes d’égalité. Le système de conditions d’optimalité ainsi obtenu peut alors être linéarisé. La résolution de ce nouveau système linéaire constitue l’autre partie importante de cette thèse. Les méthodes de pénalisation et d’Uzawa sont d’abord exposées ; elles serviront de point de comparaison. Nous construisons par la suite un algorithme basé sur un GCR préconditionné par une factorisation faisant apparaître un complément de Schur. Le tout est présenté d’abord dans le cas plus simple du contact sans frottement, puis en ajoutant le frottement. Des résultats numériques, parfois académiques (des cubes) et parfois industriels (des lamelles de pneus), viennent finalement appuyer les convictions que non seulement ces méthodes sont rapides, mais qu’elles sont aussi robustes.This thesis addresses the problem of large deformation elasticity with frictional contact. In the first chapter, we recall the laws to which we refer: material behavior laws, boundary conditions, contact and friction according to Coulomb’s law. The heart of the thesis begins with the discretization of boundary conditions. Mortar formulations for non-penetration and Coulomb’s law are explained in order to set the problem in a finite dimensional space. Optimality conditions are calculated for the discretized problem and are subsequently linearized according to the paradigm of active set strategy. Under this approach, the contact state is set at each iteration of Newton’s method, transforming the contact conditions in equality constraints. The system of optimality conditions thus obtained can then be linearized. The resolution of this new linear system is another important part of this thesis. Methods of penalization and Uzawa are first exposed; they serve as a comparison point. Thereafter, we introduce a preconditioned GCR algorithm based on a factorization using a Schur complement. The whole problem is presented first in the simplest case of contact without friction, and, in the last chapter, with friction. Numerical results, sometimes academic (cubes) and sometimes industrial (slides tires) finally support the belief that not only these methods are efficient, but they are also robust