Summary: | Les simulations en mécanique du solide présentent des difficultés comme le traitement de l'incompressibilité ou les non-linéarités dues aux grandes déformations, aux lois de comportement et de contact. L'objectif principal de ce travail est de proposer des méthodes éléments finis capables de gérer l'incompressibilité en grandes déformations en utilisant des éléments de faible ordre. Parmi les approches de la littérature, les formulations mixtes offrent un cadre théorique intéressant. Dans ce travail, une formulation mixte à trois champs (déplacements, pression, gonflement) est introduite. Dans certains cas, cette formulation peut être condensée en formulation à deux champs. Cependant, il est connu que le problème discret obtenu par une approche éléments finis de type Galerkin n'hérite pas automatiquement de la condition de stabilité “inf-sup” du problème continu : les éléments finis utilisés, et notamment les ordres d'interpolation doivent être choisis de sorte à vérifier cette condition de stabilité. Cependant, il est possible de s'affranchir de cette contrainte en ajoutant des termes de stabilisation à la formulation EF Galerkin. Cette approche permet entre autres d'utiliser des ordres d'interpolation égaux. Dans ce travail, des éléments finis stables de type P2/P1 sont utilisés comme référence, et comparés à une formulation P1/P1, stabilisée soit avec une fonction bulle, soit avec une méthode VMS (Variational Multi-Scale) basée sur un espace sous-grille orthogonal à l'espace EF. Combinées à un modèle grandes déformations basé sur des déformations logarithmiques, ces approches sont d'abord validées sur des cas académiques puis sur des cas industriels. === Simulations in solid mechanics exhibit difficulties as dealing with incompressibility or nonlinearities due to finite strains, constitutive laws and contact. The basic motivation of our work is to propose efficient finite element methods capable of dealing with incompressibility in finite strain context, and using low order elements. Among the approaches in the literature, mixed formulations offer an interesting theoretical framework. In this work, a three-field mixed formulation (displacement, pressure, volumetric strain) is investigated. In some cases, this formulation can be condensed in a two-field formulation. However, it is well-known that the discrete problem given by the Galerkin finite element technique, does not inherit the “inf-sup” stability condition from the continuous problem: the finite elements used, and in particular the interpolation orders must be chosen so as to satisfy this stability condition. However, it is possible to circumvent it, by adding terms stabilizing the FE Galerkin formulation. The latter approach allows the use of equal order interpolation. In this work, stable finite elements of type P2/P1 are used as reference, and compared to a P1/P1 formulation, stabilized with a bubble function, or with a VMS method (Variational Multi-Scale) based on a sub-grid-space orthogonal to the FE space. Combined to a finite strain model based on logarithmic strain, these approaches are first validated on academic cases and then on industrial cases.
|