Ce travail apporte une contribution à la modélisation numérique du forgeage à chaud à grande vitesse. Dans ce cadre, une formulation dynamique partiellement explicite en accélération/pression est proposée. Tout d'abord, plusieurs formulations implicites vitesse/pression (quasi-statique et dynamique) sont présentées et discutées. Leur solution est considérée comme référence par la suite. La formulation dynamique partiellement explicite est ensuite présentée dans le contexte viscoplastique. La vitesse courante étant supposée connue, le traitement du comportement non-linéaire du matériau est évité. Ainsi, le système associé au problème est linéaire. La discrétisation éléments finis est basée sur le mini-élément P1 + /P1. Son application à une approche dynamique mixte accélération/pression est examinée et quelques approximations sont effectuées afin de simplifier la formulation. Cette nouvelle approche a été implémentée avec succès dans une version du code FORGE2. Des comparaisons avec les approches implicites de référence ont été réalisées pour des tests de traction et de compression et prouvent la pertinence de notre approche. Une amélioration concernant la perte de volume est constatée. De plus, une réduction du CPU d'au moins 20% est obtenue lorsque le même pas de temps est utilisé dans les cas implicite et partiellement explicite. La procédure est généralisée au cas élasto-viscoplastique. Afin de simplifier la résolution du problème, une nouvelle stratégie est explorée. Elle consiste à résoudre un système linéaire uniquement en pression. L'accélération est calculée ensuite explicitement en fonction de la pression. Les estimations réalisées nous font envisager des gains significatifs de temps de calcul en 2D et en 3D. La formulation partiellement explicite proposée représente donc une alternative pour la simulation du forgeage à chaud.

Nous présentons dans cette contribution les techniques que nous avons mises en oeuvre pour paralléliser un code éléments finis 2D dédié à la simulation du forgeage de pièces axisymétriques. Les modèles de comportement conduisent à résoudre des équations de type Stokes généralisé, exprimées sous forme mixte en vitesse et pression. La discrétisation spatiale est effectuée par une méthode éléments finis originale basée sur une stabilisation du MINI-élément P1+P1<br />Cette approche mène à des systèmes linéaires symétriques non définis positifs que l'on peut inverser avec un solveur itératif. L'introduction de préconditionneurs par factorisation incomplète LDL(0) ainsi que l'optimisation de la résolution non-linéaire nous permet de concurrencer une méthode directe sur un maillage de plus de 3000 noeuds.<br />Une stratégie de parallélisation SPMD couplée avec un solveur itératif avec préconditionnement diagonal aboutit à, un solveur parallèle simple et efficace, ne dépendant ni de la partition ni du nombre de domaines. Différentes stratégies sont envisagées pour développer des factorisations incomplètes parallèles. Un préconditionneur additif de Schwarz est notamment proposé. Celui-ci est construit à partir des matrices locales, complétées sur leur diagonale aux interfaces et avec un coefficient de sur-relaxation. Des résultats sur des simulations industrielles sont donnés par une machine parallèle à mémoire partagée. Ceux-ci, obtenus sur des problèmes 2D et 3D, prouvent la pertinence de notre approche.<br />Les stratégies développées permettent ainsi de réduire de manière significative les temps de simulation de la majorité des cas industriels. Elles permettent aussi d'élargir les champs d'application des codes de calcul à des simulations industrielles très complexes ou avec des maillages de plus de 15000 noeuds en 2D

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.

Dans ce travail nous nous intéressons à la résolution du problème d'évolution de Navier-Stokes incompressible totalement discrétisé en temps et en espace, en dimension deux par une méthode à deux grilles. Dans un premier temps, nous étendons la méthode à deux grilles, appliquée par V. Girault et J.-L. Lions au problème de Navier-Stokes instationnaire semi-discrétisé au problème de Navier-Stokes totalement discrétisé en temps (par un schéma d'ordre un) et en espace (par une méthode d'éléments finis d'ordre un). Dans la première étape, le problème non-linéaire est discrétisé en espace et en temps sur une grille grossière de pas d'espace H avec un pas de temps Delta t. Puis dans la deuxième étape, le problème, linéarisé autour de la vitesse u_H calculée à l'étape précédente, est discrétisé en espace sur une grille fine de pas d'espace h et le même pas de temps. L'idée de la méthode à deux grilles est que, sous des hypothèses adéquates, la contribution de u_H à l'erreur dans le terme non-linéaire en espace, est mesurée en norme L^2 en espace et en temps et a un ordre plus élevé que si elle était mesurée en norme H^1. Dans un deuxième temps, vu que le but est de gagner en ordre de convergence de l'erreur totale du schéma ainsi qu'en complexité nous étudions un schéma à deux grilles d'ordre deux en temps du problème totalement discrétisé en temps et en espace de Navier-Stokes. Nous présentons les résultats suivants: dans le cas de la résolution du schéma d'ordre un en temps, si h = H^2 = Delta t, alors l'erreur globale de l'algorithme à deux grilles est de l'ordre de h. Dans le cas du schéma d'ordre deux en temps, si h^2 = (Delta t)^2 = H^3, alors l'erreur globale de l'algorithme à deux grilles est de l'ordre de h^2: résultats identiques à ceux de la résolution directe du problème non-linéaire sur une grille fine.