Le travail présenté dans cette thèse porte sur l'estimation du flot optique par méthodes variationnelles en petits et en grands déplacements. Nous proposons un modèle basé sur la combinaison locale-globale à laquelle nous ajoutons la prise en compte des variations de la luminosité. La particularité de ce manuscrit réside dans l'utilisation de la méthode des éléments finis pour la résolution des équations. En effet, cette méthode se fait pour le moment très rare dans le domaine du flot optique. Grâce à ce choix de résolution, nous proposons d'implémenter un contrôle local de la régularisation ainsi qu'une adaptation de maillage permettant d'affiner la solution au niveau des arêtes de l'image. Afin de réduire les temps de calcul, nous parallélisons les programmes. La première méthode implémentée est la méthode parallèle en temps appelée pararéel. En couplant un solveur grossier et un solveur fin, cet algorithme permet d'accélérer les calculs. Pour pouvoir obtenir un gain de temps encore plus important et également traiter les séquences en haute définition, nous utilisons ensuite une méthode de décomposition de domaine. Combinée au solveur massivement parallèle MUMPS, cette méthode permet un gain de temps de calcul significatif. Enfin, nous proposons de coupler la méthode de décomposition de domaine et le pararéel afin de profiter des avantages de chacune. Dans une seconde partie, nous appliquons tous ces modèles dans le cas de l'estimation du flot optique en grands déplacements. Nous proposons de nous servir du pararéel afin de traiter la non-linéarité de ce problème. Nous terminons par un exemple concret d'application du flot optique en restauration de films. / The work presented in this thesis focuses on the estimation of the optical flow through variational methods in small and large displacements. We propose a model based on the combined local-global strategy to which we add the consideration of brightness intensity variations. The particularity of this manuscript is the use of the finite element method to solve the equations. Indeed, for now, this method is really rare in the field of the optical flow. Thanks to this choice of resolution, we implement an adaptive control of the regularization and a mesh adaptation to refine the solution on the edges of the image. To reduce computation times, we parallelize the programs. The first method implemented is a parallel in time method called parareal. By combining a coarse and a fine solver, this algorithm speeds up the computations. To save even more time and to also be able to handle high resolution sequences, we then use a domain decomposition method. Combined with the massively parallel solver MUMPS, this method allows a significant reduction of computation times. Finally, we propose to couple the domain decomposition method and the parareal to have the benefits of both methods. In the second part, we apply all these models to the case of the optical flow estimation in large displacements. We use the parareal method to cope with the non-linearity of the problem. We end by a concrete example of application of the optical flow in film restoration.

Ce travail de thèse porte sur le développement d’outils de modélisation pour le contrôle non destructif (CND) par courants de Foucault (CF). Il existe actuellement une tendance à la mise en œuvre de capteurs souples qui représentent une solution pertinente pour inspecter des pièces ayant une surface complexe. L’objectif principal de cette thèse est l’élaboration de techniques permettant la modélisation de ce type de capteurs au sein de la méthode des éléments finis (MEF).Lors de la modélisation d’un capteur souple avec la MEF, trois problématiques se manifestent. La première concerne le maillage des milieux fins qui apparaissent dans ce type de configuration (distance capteur-pièce contrôlée, bobine plate, revêtement fin...). Le maillage de ces régions par des éléments simpliciaux peut poser des problèmes numériques (éléments déformés quand un maillage grossier est considéré et grand nombre d’inconnues quand un maillage fin est utilisé). La deuxième problématique concerne le déplacement du capteur. Si les différents sous-domaines géométriques (air, pièce, capteur...) sont convenablement remaillés pour chaque position du capteur, le temps total pourra être pénalisant. La troisième problématique relative à la modélisation d’un capteur souple porte sur l’imposition du courant dans des bobines déformées.Une comparaison de différentes approches nous a conduit à retenir la méthode overlapping, qui permet de considérer simultanément des milieux fins avec maillages non-conformes. Cette méthode permet d’effectuer le recollement de deux surfaces non planes et/ou de géométries différentes. La méthode overlapping a été implantée dans deux formulations duales (magnétique et électrique) en 2D et 3D intégrées dans le code de calcul DOLMEN (C++) du LGEP. La méthode overlapping a été validée pour la prise en compte de plusieurs types de régions minces (air, milieux conducteurs, milieux magnétiques, bobines plates...). La modélisation des capteurs souples nécessite aussi l’implantation d’une technique permettant d’imposer correctement les courants dans un inducteur de forme arbitraire. Une technique a été sélectionnée et programmée, d’abord pour des bobines classiques (volumiques) déformées puis pour des bobines plates flexibles. Afin de valider les développements réalisés, différentes configurations de test ont été traitées et les résultats comparés avec des solutions de référence analytiques ou expérimentales. / The objective of this thesis is the development of modeling tools for eddy current testing (ECT). Currently there is a tendency to use flexible sensors which represent a viable solution for inspecting parts with a complex surface. The main objective of this thesis is the development of techniques for taking into account this kind of sensors within the finite element method (FEM).When modeling a flexible sensor with the FEM, three issues have to be considered. The first one is related to the mesh of thin regions that appear in this kind of configuration (sensor-inspected part distance, thin coating...). The meshing of these regions with simplicial elements can cause numerical problems (distorted elements when a coarse mesh is considered and high number of unknowns when a fine mesh is used). The second issue is related to the displacement of the sensor on the part surface. If the different subdomains (air, part, sensor ...) are properly remeshed for each position of the sensor, the required time can be penalizing. The third issue, related to the modeling of a flexible sensor, is the computation of the current in distorted sensor coils.A comparison of different approaches has led to select the overlapping element method, which allows to simultaneously consider the thin regions with nonconforming meshes. This method allows to perform the connection of two surfaces which can be non-planar and/or have different geometries. The overlapping method has been implemented in two dual formulations (magnetic and electric) in 2D and 3D and integrated in the computation code (C++) DOLMEN of LGEP. The overlapping method has been validated for several kinds of thin regions (air, conductive regions, magnetic regions, flat coils...). The modeling of flexible sensors also requires the establishment of a technique for properly imposing the current in an inductor of arbitrary shape. A technique has been selected and implemented for conventional (volumic) distorted coils but also for flexible flat coils. Different test configurations have been considered in order to validate the developments and the results have been compared with analytical references or experimental solutions.

Les réfractaires électrofondus, qui constituent l’objet de ce travail, appartiennent au système alumine-zircone-silice. Ils sont obtenus par coulée dans des moules à des températures supérieures à 2000°C, rendant très difficile toute instrumentation. De nombreux phénomènes intrinsèques au matériau interviennent lors du refroidissement qui suit la coulée. Parmi ces derniers, cette recherche a essentiellement porté sur la transformation de phase (de tétragonale à monoclinique) de la zircone et aux phénomènes associés (gonflement, plasticité,…).A partir d’essais mécaniques à haute température réalisés en laboratoire, les lois de comportements thermiques et mécaniques ont été caractérisées et modélisées en cours de transformation de la zircone. La plasticité à très bas seuil de contrainte observée a, en particulier, été décrite par une vitesse de déformation dérivée du modèle de Leblond, une fonction de rendement de type Cam-clay sans consolidation et une fonction de rendement définissant l’avancement de la transformation en fonction de la température. Après implémentation dans un code de calcul par éléments finis et validation par confrontation avec des résultats d’essais sous contraintes multiaxiales, ce modèle a été assemblé aux autres composantes du comportement mécanique (fluage, élasticité,…), pour décrire l’ensemble des phénomènes thermomécaniques observés lors du refroidissement.Parallèlement, des coulées de blocs en laboratoire, instrumentées par des thermocouples et des capteurs d’émission acoustique, ont permis de reconstruire par simulation numérique l’évolution du champ de température à l’intérieur de la dalle au cours du refroidissement. L’enthalpie de solidification et celle associée à la transformation de phase, préalablement quantifiée par ATD, ont été prises en compte. L’application du modèle mécanique complet, associant toutes les composantes du comportement, a permis de calculer l’évolution du champ de contraintes généré par les gradients thermiques en fonction du temps et, en particulier, de mettre en évidence le rôle essentiel joué par la transformation de phase sur la relaxation des contraintes. / Fused-cast refractories, which are concerned by this work, belong to the alumina-zirconia-silica system. They are obtained by casting in molds at temperatures higher than 2000°C, that make very difficult any instrumentation. Many phenomena intrinsic to the material occur during cooling-down after casting. Among these latter, this research essentially focused on the phase transformation (from tetragonal to monoclinic) of zirconia and the associated phenomena (swelling, plasticity,...).From high temperature mechanical tests performed in laboratory, the thermal and mechanical behavior laws were characterized and modeled during the zirconia transformation. Plasticity at very low stress threshold was observed. A Leblond type model has been extended by introducing a Cam-clay yield function without consolidation. In this model, the progress of the transformation is controlled by the evolution of the temperature. This model was complemented by other components of the mechanical behavior (creep, elasticity, ...). It has been validated by experimental tests under multiaxial loadings that replicate the main thermomechanical phenomena observed during cooling.In parallel, blocks casted in laboratory conditions, instrumented with thermocouples and acoustic emission sensors, allowed a numerical simulation of the change in temperature field within the block during cooling-down. This simulation took into account the solidification enthalpy and the enthalpy associated to the phase transformation, previously quantified by DTA. The implementation of the complete mechanical model integrating all the behavior components led to a calculation of the stress field changes generated by thermal gradients as a function of time and, in particular, to highlight the essential role played by the phase transformation on stress relaxation.

Ce travail porte sur le problème des partitions minimales, à l'interface entre théorie spectrale et optimisation de forme. Une introduction générale précise le problème et présente des résultats, principalement dûs à B. Helffer, T. Hoffmann-Ostenhof et S. Terracini, qui sont utilisés dans le reste de la thèse.Le premier chapitre est une étude spectrale asymptotique du laplacien de Dirichlet sur une famille de domaines en dimension deux qui tend vers un segment. L'objectif est d'obtenir une localisation des lignes nodales dans la limite des domaines minces. En appliquant les résultats de Helffer, Hoffmann-Ostenhof et Terracini, on montre ainsi que les domaines nodaux des premières fonctions propres forment des partitions minimales.Le deuxième chapitre étudie les valeurs propres de certains opérateurs de Schrödinger sur un domaine plan avec condition au bord de Dirichlet. On considère des opérateurs qui ont un potentiel électrique nul et un potentiel magnétique d'un type particulier, dit d'Aharonov-Bohm, avec des singularités en un nombre fini de points appelés pôles. On démontre que les valeurs propres dépendent continuement des pôles. Dans le cas de pôles distincts et éloignés du bord, on prouve que cette dépendance est analytique lorsque la valeur propre est simple. On exprime de plus une condition suffisante pour que la fonction qui aux pôles associe une valeur propre présente un point critique. On utilise alors la caractérisation magnétique des partitions minimales pour montrer que l'énergie minimale est une valeur critique d'une de ces fonctions.Le troisième chapitre est un article écrit en collaboration avec Virginie Bonnaillie-Noël. Il porte sur une famille d'exemples, les secteurs angulaires de rayon unité et d'ouverture variable, dont on tente de déterminer les partitions minimales. On applique pour cela les théorèmes généraux rappelés dans l'introduction afin de déterminer les partitions nodales qui sont minimales. On s'intéresse plus particulièrement aux partitions minimales en trois domaines. En appliquant les idées du deuxième chapitre, on montre que pour certaines valeur de l'angle, il n'existe aucune partition minimale qui soit symétrique par rapport à la bissectrice du domaine. D'un point de vue quantitatif, on obtient des encadrements précis de l'énergie minimale.Le quatrième chapitre consiste en l'étude des partitions minimales de tores plats dont on fait varier le rapport entre longueur et largeur. On utilise une méthode numérique très différente de celle du troisième chapitre, basée sur un article de B. Bourdin, D. Bucur et É. Oudet. Elle consiste en une relaxation suivie d'une optimisation par un algorithme de gradient projeté. On peut ainsi tester des résultats théoriques antérieurs. Les résultats présentés suggèrent de plus la construction explicite de familles de partitions (en liaison avec des pavages du tore) qui donnent une nouvelle majoration de l'énergie minimale.Un dernier chapitre de perspectives présente plusieurs applications possibles des méthodes décrites dans la thèse. / This work is concerned with the problem of minimal partitions, at the interface between spectral theory and shape optimization. A general introduction gives a precise statement of the problem and recall results, mainly due to B. Helffer, T. Hoffmann-Ostenhof and S.Terracini, that are used in the rest of the thesis.The first chapter is an asymptotic spectral study of the Dirichlet Laplacian on a familly of two-dimensional domains converging to a line segment. The aim is to localize the nodal lines when the domains become very thin. With the help of the results of Helffer, Hoffmann-Ostenhof, and Terracini, we then show that the nodal domains of the first eigenfunctions give minimal partitions.The second chapter studies the eigenvalues of some Schrödinger operators on a domain with Dirichlet boundary conditions. We consider operators that have no electric potential and a so-called Aharonov-Bohm magnetic potential, which has singularities at a finite number of points called poles. We prove that the eigenvalues are continuous functions of the poles. When the poles are distinct and far from the boundary, we prove that this function is analytic, assuming the eigenvalue is simple. We also give a sufficient condition for the function to have a critical point. Using the magnetic characterization of minimal partitions, we show that the minimal enery is a critical value for one of these functions.The third chapter in an article written in collaboration with Virginie Bonnaillie-Noël. It studies minimal partitions for sectors of unit radius with a variable angular opening. We apply the general results presented in the introduction, together with numerical computations, to determine nodal partitions that are minimal. We focus on partitions into three domains. Using ideas from the second chapter, we show that, for some values of the angle, there is no minimal partition that is symmetric with respect to the bisector. Form a quantitative point of view, we obtain precise bounds on the minimal energy.The fourth chapter studies the minimal partitions of flat tori in function of the ratio between width and length. We use a numerical method that is quite different from chapter three, and is based on an article by B. Bourdin, D. Bucur, and É. Oudet. It consists in a relaxation of the problem, followed by optimization with the help of a projected gradient algorithm. The results shown here additionally suggest explicit families of partitions, which consist in tilings of tori by polygons, that give upper bounds on the minimal energy. In the last chapter we consider several possible applications of the methods described in the thesis.

Une dégradation des propriétés mécaniques des arcs surperélastiques en alliages à mémoire de forme à base NiTi, utilisés dans les traitements orthodontiques, a été observée après absorption d'hydrogène. L’effet de l’hydrogène a été étudié, dans un premier temps, sur le comportement global des arcs à l’aide des essais de traction et dans un deuxième temps, sur le comportement local à l’aide des essais de nanoindentation. Pour modéliser le comportement des AMFs après absorption d'hydrogène, une première approche a été proposée, en introduisant la dépendance des paramètres de la transformation martensitique à la concentration moyenne d'hydrogène, observée dans les courbes contrainte-déformation obtenues, dans un modèle existant dédié aux AMFs. Une deuxième approche consiste à proposer un modèle thermomécanochimique couplé. Pour ce faire, un nouveau potentiel thermodynamique est défini en introduisant la déformation chimique due à la présence de l'hydrogène et l'interaction de ce dernier avec les variantes de martensite. Des forces thermodynamiques sont déduites de ce potentiel et équilibrées en faisant intervenir les phénomènes dissipatifs : mécanique, thermique et chimique. L'exploitation du modèle proposé a nécessité le développement d'un élément fini spécifique adoptant la concentration d'hydrogène comme un degré de liberté supplémentaire qui prend en compte le couplage complet et la formulation proposée de l'équilibre thermique et chimique / A degradation of the mechanical properties of NiTi-based shape memory alloys superelastic archs, used in orthodontic treatments, was observed after hydrogen absorption. The effect of hydrogen was first investigated on the global behaviour of archs using tensile tests and secondly on the local behaviour using nanoindentation tests. A first approach to model the behavior of AMFs after hydrogen absorption has been proposed, by introducing the dependence of martensitic transformation parameters on the average hydrogen concentration, observed in the stress-strain curves obtained, in an existing dedicated model to SMA. A second approach is to propose a coupled thermomechanical model. A new thermodynamic potential is defined by introducing the chemical strain due to the presence of hydrogen and the interaction of the latter with martensite variants. Thermodynamic forces are deduced from this potential and balanced by involving mechanical, thermal and chemical dissipative phenomena. The exploitation of the proposed model required the development of a special finite element adopting the hydrogen concentration as an additional degree of freedom that takes into account the thermomechanical coupling and the proposed formulation of the thermal and chemical equilibrium

Nous présentons dans cette thèse nos travaux sur la réduction d'ordre appliquée à des simulations d'aérothermie. Nous considérons le couplage entre les équations de Navier-Stokes et une équations d'énergie de type advection-diffusion. Les paramètres physiques considérés nous obligent à considéré l'introduction d'opérateurs de stabilisation de type SUPG ou GLS. Le but étant d'ajouter une diffusion numérique dans la direction du champs de convection, afin de supprimer les oscillations non-phyisques. Nous présentons également notre stratégie de résolution basée sur la méthode des bases réduite (RBM). Afin de retrouver une décomposition affine, essentielle pour l'application de la RBM, nous avons implémenté une version discrète de la méthode d'interpolation empirique (EIM). Cette variante permet de la construction d'approximation affine pour des opérateurs complexes. Nous utilisons notamment cette méthode pour la réduction des opérateurs de stabilisations. Cependant, la construction des bases EIM pour des problèmes non-linéaires implique un grand nombre de résolution éléments finis. Pour pallier à ce problème, nous mettons en oeuvre les récents développement de l'algorithme de coconstruction entre EIM et RBM (SER). / We present in this thesis our work on model order reduction for aerothermal simulations. We consider the coupling between the incompressible Navier-Stokes equations and an advection-diffusion equation for the temperature. Since the physical parameters induce high Reynolds and Peclet numbers, we have to introduce stabilization operators in the formulation to deal with the well known numerical stability issue. The chosen stabilization, applied to both fluid and heat equations, is the usual Streamline-Upwind/Petrov-Galerkin (SUPG) which add artificial diffusivity in the direction of the convection field. We also introduce our order reduction strategy for this model, based on the Reduced Basis Method (RBM). To recover an affine decomposition for this complex model, we implemented a discrete variation of the Empirical Interpolation Method (EIM) which is a discrete version of the original EIM. This variant allows building an approximated affine decomposition for complex operators such as in the case of SUPG. We also use this method for the non-linear operators induced by the shock capturing method. The construction of an EIM basis for non-linear operators involves a potentially huge number of non-linear FEM resolutions - depending on the size of the sampling. Even if this basis is built during an offline phase, we usually can not afford such expensive computational cost. We took advantage of the recent development of the Simultaneous EIM Reduced basis algorithm (SER) to tackle this issue.

Electromagnetic coupling phenomena between overhead power transmission lines and other nearby structures are inevitable, especially in densely populated areas. The undesired effects resulting from this proximity are manifold and range from the establishment of hazardous potentials to the outbreak of alternate current corrosion phenomena. The study of this class of problems is necessary for ensuring security in the vicinities of the interaction zone and also to preserve the integrity of the equipment and of the devices there present. However, the complete modeling of this type of application requires the three- -dimensional representation of the region of interest and needs specific numerical methods for field computation. In this work, the modeling of problems arising from the flow of electrical currents in the ground (the so-called conductive coupling) will be addressed with the finite element method. Those resulting from the time variation of the electromagnetic fields (the so-called inductive coupling) will be considered as well, and they will be treated with the generalized PEEC (Partial Element Equivalent Circuit) method. More specifically, a special boundary condition on the electric potential is proposed for truncating the computational domain in the finite element analysis of conductive coupling problems, and a complete PEEC formulation for modeling inductive coupling problems is presented. Test configurations of increasing complexities are considered for validating the foregoing approaches. These works aim to provide a contribution to the modeling of this class of problems, which tend to become common with the expansion of power grids. / Les phénomènes de couplage électromagnétique entre les lignes aé- riennes de transmission d\'énergie et des structures voisines sont inévitables, surtout dans les zones densément peuplées. Les effets indésirables découlants de cette proximité sont variés, allant de l\'établissement des tensions dangereuses à l\'apparition de phénomènes de corrosion dus au courant alternatif. L\'étude de cette classe de problèmes est nécessaire pour assurer la sécurité dans les alentours de la zone d\'interaction et aussi pour préserver l\'intégrité des équipements et des dispositifs présents. Cependant, la modélisation compl ète de ce type d\'application implique la représentation tridimensionnelle de la région d\'intérêt et nécessite des méthodes numériques de calcul de champs spécifiques. Dans ces travaux, des problèmes liés à la circulation de courants électriques dans le sol (ou de couplage dit conductif) seront abordés avec la méthode des éléments finis. D\'autres problèmes résultants de la variation temporelle des champs électromagnétiques (ou de couplage dit inductif) seront aussi considérés et traités avec la méthode PEEC (Partial Element Equivalent Circuit) généralisée. Plus précisément, une condition limite particulière sur le potentiel électrique est proposée pour tronquer le domaine de calcul dans l\'analyse par éléments finis des problèmes de couplage conductif et une formulation PEEC complète pour la modélisation de problèmes de couplage inductif est présentée. Des con gurations tests de complexités croissantes sont considérées pour valider les approches précédentes. Ces travaux visent ainsi à apporter une contribution à la modélisation de cette classe de problèmes, qui tendent à devenir communs avec l\'expansion des réseaux électriques.

Cette thèse s'intéresse à la génération, à l'utilisation et à l'évaluation des maillages hex-dominants, composés d'hexaèdres et de tétraèdres, dans la cadre de la simulation numérique par la méthode des éléments finis. Les éléments finis hexaédriques sont souvent préférés aux éléments tétraédriques car ils offrent un meilleur ratio entre précision et temps de calcul dans un certain nombre de situations. Cependant, si la génération automatique de maillages tétraédriques est aujourd'hui un domaine bien maîtrisé, ce n'est pas le cas de la génération de maillages hexaédriques alignés avec le bord, qui reste un problème largement ouvert. En l'absence de progrès significatifs, les approches actuelles se contentent de maillages hex-dominants afin de tirer parti des performances supérieures des hexaèdres et de la flexibilité géométrique des tétraèdres, qui rend possible le maillage automatique. Dans une première partie, nous développons des algorithmes robustes pour la génération de maillages hex-dominants à partir de champs de directions, notamment pour l'isolement et le remplissage des régions difficiles à mailler (singularités et autres dégénérescences). Dans la seconde partie, nous essayons de déterminer dans quelles situations et dans quelle mesure les maillages hexaédriques, et hex-dominants générés précédemment, sont plus intéressants que les maillages tétraédriques. Ceci implique spécifiquement d'étudier plusieurs manières d'effectuer des simulations par éléments finis avec les maillages hybrides, dont une approche où nous utilisons des contraintes de continuité pour maillages non-conformes. Pour mesurer l'influence du maillage sur l'approximation des solutions, nous proposons une nouvelle méthode d'échantillonnage pour calculer très efficacement des distances globales entre solutions éléments finis définies sur des domaines compliqués / This thesis focuses on generation, usage and evaluation of hex-dominant meshes, which are made of hexaehedra and tetrahedra, in the context of the finite element method. Hexahedron finite elements are often preferred to tetrahedron elements because they offer a better compromise between accuracy and computation time in certain situations. However, if tetrahedral meshing is a well mastered subject, it is not the case of hexahedral meshing. Generating hexahedral meshes with elements aligned to the borders is still an open and difficult problem. Meanwhile, current automated approaches can use hex-dominant meshes in order to take advantage of both hexahedron accuracy and geometrical flexibility of tetrahedra. In the first part, we develop robust algorithms for the generation of hex-dominant meshes with elements aligned with the borders. Specifically, we propose a method to extract and fill the areas where hexahedral meshing is difficult (singularities and degeneracies). In the second part, we try to identify and to quantify the advantages of hexahedral and hex-dominant meshes over tetrehedral ones. This requires to study various ways to apply the finite element method on hybrid meshes, including one in which we propose to use continuity constraints on hexahedral-tetrahedral non-conforming meshes. To measure the impact of meshes on the finite element accuracy, we develop a new sampling method which allows to compute efficiently global distances between finite element solutions defined on complicated 3D domains

Le forgeage à chaud est un procédé de formage des métaux utilisé pour former des matériaux qui sont difficiles à former à froid ainsi que pour réaliser des géométries complexes. La réduction de la limite d’élasticité à haute température et une augmentation subséquente de l’aptitude à la mise en forme constituent le principal mécanisme à l’origine du procédé. Les méthodes numériques constituent un moyen efficace de prédire les états de contrainte / déformation du produit à différentes étapes de la mise en forme. Bien que les méthodes classiques soient suffisamment précises pour fournir une représentation appropriée du procédé, elles ont tendance à être coûteuses en ressources informatiques. Cela limite leur utilisation dans des cas concret, en particulier pour des études d’optimisation du procédé. L’approche pseudo inverse (API), développée dans le contexte du forgeage à froid 2D axisymétrique, fournit une estimation rapide des champs de contrainte et de contrainte dans le produit final pour une forme initiale donnée. Dans ce travail, l’API est étendue pour inclure les effets thermiques et visco-plastiques dans le procédé de forgeage ainsi que dans le cas général 3D. Les résultats sont comparés aux codes commerciaux disponibles basés sur les approches classiques pour montrer l’efficacité et les limites de l’API. Les résultats obtenus indiquent que l’API est un outil assez efficace pouvant être utilisé à la fois pour des simulations 2D et 3D du forgeage à chaud. / Hot forging is a metal forming process used to form difficult-to-form materials as well as to achieve complex geometries. The reduction of yield stress at high temperatures and a subsequent increase in formability is the primary mechanism that drives the process. Numerical methods provide an efficient means to predict the material yield and the stress/strain states of the product at different stages of forming. Although classical methods are accurate enough to provide a suitable representation of the process, they tend to be computationally expensive. This limits its use in practical cases especially for process optimization. Pseudo Inverse Approach (PIA) developed in the context of 2D axisymmetric cold forming, provides a quick estimate of the stress and strain fields in the final product for a given initial shape. In this work, the PIA is extended to include the thermal and viscoplastic effects on the forging process as well as to the general 3D case. The results are compared with commercially available software, based on the classical approaches, to show the efficiency and the limitations of PIA. The results obtained indicate that PIA is a quite effective tool that can be used for both 2D and 3D simulations of hot forging.

Le but du travail présenté est la compréhension puis la simulation numérique des mécanismes d'érosion des augets de turbine Pelton par impacts répétés de gouttes d'eau dans le but de prédire la durée de vie des composants. Tout d'abord, les phénomènes de propagation d'ondes dans les milieux fluide et solide sont étudiés. Cela permet de mettre en lumière l'évolution temporelle et la distribution spatiale des pressions de contact, et l'apparition de microjets par éjection supersonique du fluide au contact. Les études expérimentales de l'érosion par gouttes d'eau traduisent un dommage basé sur la fissuration par fatigue. Des simulations numériques en dynamique rapide couplées fluide-structure sont alors effectuées. Le domaine solide est discrétisé par la Méthode des éléments Finis (MEF), et le domaine fluide par la méthode Smoothed Particle Hydrodynamics (SPH), qui est une méthode particulaire (sans maillage) particulièrement adaptée aux grandes distorsions et au suivi des surfaces libres. L'analyse des états de contraintes vient corroborer la nature cyclique de l'endommagement. La simulation d'érosion est alors réalisée à l'aide de critères de fatigue multiaxiaux. Le choix se porte vers un premier critère général de l'American Society of Mechanical Engineers (ASME), utilisant les valeurs principales des différences de contraintes au cours du temps. Le second choix concerne un critère à plan critique : le critère de Dang Van 2. Il traite séparément la contrainte hydrostatique et le cisaillement alterné maximal local. Ces critères permettent de définir les régions érodées du solide au bout d'un nombre d'impact donné, ce qui fait de cette démarche une méthode prédictive. Une étude paramétrique pour différentes tailles de gouttes et vitesses d'impact est ensuite réalisée, puis on évalue l'influence de la présence d'une couche de coating. / The goal of this work is the comprehension and the numerical simulation of erosion caused by repeated droplet impact on Pelton turbine buckets, to predict the lifetime of these components. First, waves propagation phenomenon inside fluid and solid domains are presented, which allows determining the time evolution and spatial distribution of contact pressure, and the birth of lateral microjets by supersonic ejection of the fluid on the contact. Experimental studies of erosion by droplet impact highlight a fatigue cracking-based erosion mechanism. Then, coupled FSI computation are performed. The solid subdomain is discretized by the Finite Element Method (FEM), and the fluid subdomain by the Smoothed Particle Hydrodynamics (SPH), which is a particle method (meshless) effectively recommended for large distortions and free surface tracking. Stress analysis confirms the cyclic nature of the damage mechanism, and erosion simulation is performed using multiaxial fatigue criteria. The first selected criterion is a general one from the American Society of Mechanical Engineers (ASME) using principal values of stress differences over time. The second one is the Dang van 2 criterion, belonging to the family of critical plane criteria. This criterion considers separately the effects due to hydrostatic stress on one hand, and the ones induced by maximum local shear on the other. These two criteria are used to building the equivalent eroded zones of the solid subdomain for a given number of impacts, which allows to qualify this procedure as a predictive predictive. Finally, a parametric study for different droplet sizes and velocites is computed, and the effects of a coating layer are investigated.