L'objet de cette thèse est d'étudier un problème de la magnétostatique tridimensionnel. On propose trois formulations mixtes couplant une méthode d'éléments finis pour tenir compte du milieu hétérogène et une méthode éléments de frontière pour le milieu extérieur homogène. Pour la méthode intégrale on a utilisé les équations de Calderon, l'opérateur de Neumann-Dirichlet ou d'autres opérateurs intégraux. L'utilisation des éléments d'arête de Nédélec pour le champ magnétique, et les éléments de face de Raviart pour l'induction magnétique permet d'utiliser des méthodes éléments finis conformes. Des résultats numériques ont permis de valider ces méthodes. La deuxième partie a porté sur la comparaison de diverses discrétisations pour l'opérateur de Poincaré-Steklov. Ces méthodes ont été comparées sur une formulation de la magnétostatique. Enfin, on propose des formulations discontinues du problème de la magnétostatique avec des conditions aux limites. On montre que ces formulations sont consistantes et des estimations d'erreur sont obtenues.

Le caoutchouc industriel est présent dans de nombreuses applications, de la plus basique à la plus technique. Son mode de défaillance prépondérant est la rupture due à la propagation de fissures par fatigue : les sollicitations mécaniques ainsi que l'atmosphère extérieure provoquent dans un premier temps l'apparition d'une amorce de fissure. Sous l'effet des sollicitations mécaniques, cette fissure se propage jusqu'à rupture de la pièce. L'objectif de ce travail est de faciliter la simulation numérique de la propagation de fissures dans les élastomères. Pour cela, on utilise la méthode des éléments finis étendus X-FEM. Cette méthode a été développée afin de limiter le recours au remaillage dans le cadre de la fissuration des métaux. En outre, elle permet d'enrichir l'approximation éléments finis par des fonctions provenant de la physique du problème. La première partie de ce travail consiste à adapter cette méthode à la mécanique non-linéaire de la rupture. On s'intéresse en particulier au choix de la formulation de résolution ainsi qu'à la recherche de fonctions d'enrichissement adaptées. Dans un deuxième temps, on s'intéresse à l'enrichissement des formulations mixtes pour la gestion de la contrainte d'incompressibilité. Des stratégies ont été développées afin de préserver la stabilité de ces formulations. Ces enrichissements permettent la vérification de la condition inf-sup dans le cas de trous, des inclusions et des fissures sous l'hypothèse des petites perturbations. Enfin, dans une troisième partie, on détaille l'application du concept de forces configurationnelles comme critère directionnel pour la propagation de fissures 2D et 3D.

Les travaux de recherche de cette habilitation se situent dans le domaine de l'Analyse Numérique des Equations aux Dérivées Partielles et portent sur la modélisation, la discrétisation, l'analyse a priori et a posteriori de schémas et la simulation numérique de différents problèmes issus de la mécanique. Un fil conducteur de ces travaux est l'utilisation et l'étude des méthodes d'éléments finis (conformes, non-conformes, mixtes, de Galerkin discontinus, stabilisés) et des formulations mixtes. Les domaines d'application abordés sont la mécanique des solides élastiques, l'ingénierie pétrolière et la mécanique des fluides, newtoniens et non-newtoniens. Ainsi, des problèmes d'élasticité linéaire, comme la discrétisation de deux modèles de plaque mince en flexion munie de conditions aux limites physiques, ont été considérés. Des écoulements anisothermes dans les milieux poreux, décrits par les équations de Darcy-Forchheimer avec un bilan d'énergie exhaustif dans les cas mono et multi-phasique, ainsi qu'un couplage thermo-mécanique puits - réservoir pétrolier ont aussi été étudiés, dans le cadre d'une collaboration industrielle avec Total. Enfin, plusieurs questions en mécanique des fluides ont été abordées, comme la discrétisation robuste des équations de Stokes par une méthode de Galerkin discontinue en lien avec les éléments finis non-conformes, le traitement des conditions aux limites non-standard pour les équations de Navier-Stokes, la modélisation hiérarchique multi-dimensionnelle des écoulements fluviaux à surface libre, la simulation réaliste des écoulements de liquides polymères et la stabilité des schémas numériques par rapport aux paramètres physiques, en particulier pour le modèle de Giesekus.

Les travaux de recherche présentés dans cette HDR sont répartis en trois thématiques complémentaires. La première concerne la modélisation des instabilités locales dans les matériaux métalliques (localisation, striction ...). Il s'agit de développer des outils théoriques et numériques de prédiction de ces phénomènes et leur validation au travers de courbes limites de formage pour des aciers ferritiques ou dual-phase. Deux approches complémentaires ont été adoptées : la première est micromécanique se basant sur la plasticité cristalline couplée à des techniques de transition d'échelles, tandis que la seconde est phénoménologique utilisant des modèles de comportement avancés que nous avons couplé à de l'endommagement. Le critère de bifurcation de Rice, qui est relié à la perte d'ellipticité des équations gouvernant le problème aux limites, a été particulièrement analysé au travers de ces deux approches. L'influence sur l'apparition d'une localisation de plusieurs paramètres matériau ou liés à la microstructure ainsi que l'effet de certains mécanismes déstabilisants ont été analysés. Il est mis en évidence, notamment, le rôle tout à fait essentiel de la formation de points de vertex aux points courants de la surface de plasticité dans la détection d'une localisation de la déformation plastique en bande de cisaillement. Le deuxième axe de recherche est relatif aux instabilités de type structure (flambage, plissement ...). Ces phénomènes se manifestent le plus souvent en présence de structures minces ou relativement élancées. Dans ce cadre, nous proposons une approche basée sur l'étude de stabilité des évolutions quasi statiques. Pour des matériaux sensibles à la vitesse de déformation, l'absence de bifurcation ainsi que d'états d'équilibre nous amènent naturellement à poser le problème comme celui de la stabilité de trajectoires quasi statiques. Pour les matériaux élasto-plastiques, cette approche est justifiée par le fait que bien souvent nous sommes en présence d'une évolution quasi statique, pour un trajet de chargement donné, même si chaque point de cette évolution représente un état d'équilibre. Le critère de stabilité unifié que nous proposons est donné par la positivité de la seconde variation de l'énergie totale et est valable pour des solides visco-élastiques, visco-plastiques ou élasto-plastiques. Plus récemment, nous avons étendu ce critère à des modèles à gradients. Enfin, la comparaison de ce critère avec les résultats existants relatifs au flambage plastique nous a permis de rediscuter le critère de non-bifurcation de Hill en relation avec le choix du modèle de plasticité approprié. Il est mis en évidence, de nouveau, le rôle important joué par le développement de points de vertex sur la surface de plasticité lors de la prédiction du flambage avec la théorie du module tangent élasto-plastique. Le troisième axe de recherche concerne le développement d'éléments finis de type coques volumiques. L'idée de ce nouveau concept est de combiner les avantages des éléments coques et 3D pour formuler des éléments particulièrement adaptés à la simulation des structures minces. Ces nouveaux éléments de coques volumiques ont l'avantage de prendre en compte la flexion suivant une direction choisie tout en conservant la formulation classique des éléments volumiques. Ils permettent ainsi de modéliser des structures comportant des parties minces et des zones épaisses, sans les difficultés habituelles de raccord de maillages contenant des éléments coques et 3D. Leur utilisation en emboutissage est très prometteuse, en particulier pour des problèmes où les effets dans l'épaisseur sont d'importance majeure. Étant sous-intégrés, pour améliorer leurs performances, les modes de hourglass générés sont alors efficacement stabilisés par des techniques récentes. Les modes de verrouillages (membrane et cisaillement transverse) sont éliminés par des techniques de projection pouvant se mettre sous le formalisme " Assumed Strain Method ".