D'un point de vu fondamental, la structure micromagnétique d'un point de Bloch est prédite depuis plus de 50 ans, mais représente cependant une singularité topologique dans le cadre de la théorie du micromagnétisme. Par conséquent, une description purement micromagnétique du point deBloch s'avère difficile. Ce manuscrit décrit les outils mathématiques et un ensemble d'algorithmes permettant de combiner un modèle d'Heisenberg classique avec des simulations micromagnétiques sur la base des éléments finis. A l'aide de ces algorithmes, nous pouvonsrigoureusement étudier les propriétés caractéristiques d'un point de Bloch d'une paroi de domaine de type vortex localisée dans un nanofil cylindrique ferromagnétique. Cette thèse décrit le pinning/depinning des points de Bloch au réseau atomique ainsi que les différents modes depropagation détectables pour des parois de domaine sous l'influence d'un champ magnétique externe, qui peut conduire le système à des vitesses supermagnoniques. / The fundamental micromagnetic structure of Bloch points was predicted by micromagnetic theory more than 50 years ago, but represents a topological singularity within the theory of micromagnetism. This complicates a pure micromagnetic description. This thesis describes thenecessary mathematical background and a set of algorithms to combine a classical atomistic Heisenberg model with micromagnetism on the basis of the finite element method. By means of those algorithms the characteristic properties of Bloch points in vortex domain walls inferromagnetic solid cylindrical nanowires can be studied rigorously. ln addition to the pinning/depinning of Bloch points at the atomistic lattice the thesis reports on the different modes of propagation detectable for vortex domain walls in that system under the influence of an externalmagnetic field, which can drive the system of the domain wall and the Bloch point with supermagnonic velocities.

Les travaux présentés dans ce mémoire s'inscrivent dans l'étude et l'amélioration des modèles d'endommagement continus régularisés (non locaux), l'objectif étant d'étudier la transition entre un champ d'endommagement continu défini sur l'ensemble d'une structure et un modèle discontinu de fissuration macroscopique.La première étape consiste en l'étude semi-analytique d'un problème unidimensionnel (barre en traction) visant à identifier une famille de lois d'interface permettant de basculer d'une solution non homogène obtenue avec un modèle continu à gradient d'endommagement vers un modèle discontinu de fissuration cohésive. Ce passage continu / discontinu est construit de telle sorte que l'équivalence énergétique entre les deux modèles soit assurée, et reste exacte quelque soit le niveau de dégradation atteint par le matériau au moment où cette transition est déclenchée.Cette stratégie est ensuite étendue au cadre 2D (et 3D) éléments finis dans le cas de la propagation de fissures rectilignes (et planes) en mode I. Une approche explicite basée sur un critère de dépassement d'une valeur « critique » de l'endommagement est proposée afin de coupler les modèles continus et discontinus au sein d'un même calcul quasi-statique par éléments finis. Enfin, plusieurs résultats de simulations menées avec cette approche couplée sont présentés. / The present work deals with the study and the improvement of regularized (non local) damage models. It aims to study the transition from a continuous damage field distributed on a structure to a discontinuous macroscopic failure model.First, an analytical one-dimensional study is carried out (on a bar submitted to tensile loading) in order to identify a set of interface laws that enable to switch from an inhomogeneous solution obtained with a continuous gradient damage model to a cohesive zone model. This continuous / discontinuous transition is constructed so that the energetic equivalence between both models remains ensured whatever the damage level reached when switching.This strategy is then extended to the bi-dimensional (and tri-dimensional) case of rectilinear (and plane) crack propagation under mode I loading conditions, in a finite element framework. An explicit approach based on a critical damage criterion that allows coupling both continuous and discontinuous approaches is then proposed. Finally, results of several simulations led with this coupled approach are presented.

Avec l'avancé des méthodes numériques, de nouvelles méthodes dites « sans maillage » sont apparues pour remédier à certaines limitations de la méthode des éléments finis. Ces méthodes ont la particularité de n’employer aucun maillage prédéfini : elles utilisent un ensemble de nœuds dispersés dans le domaine considéré et sur ses frontières. L’objectif de cette étude est de montrer l’intérêt de l’application des méthodes sans maillage basées sur les fonctions de base radiale pour la simulation des procédés de mise en forme en général et de poinçonnage rapide en particulier. Une attention particulière sera portée sur la méthode RPIM qui offre l’avantage de proposer une interpolation nodale. La démarche proposée dans ce document consiste: à rappeler succinctement les principes essentiels des méthodes sans maillage en précisant leurs avantages par rapport aux méthodes classiques, à présenter et à mettre en œuvre la technique numérique de la méthode sans maillage RPIM avec un calibrage des paramètres caractéristiques, et enfin, à traiter des exemples numériques de procédés de mise en forme avec amorçage et propagation de fissure qui confirmeront ces avantages. / In recent years, new methods named Meshfree methods have been developed to surmount limitations of the finite element method. The main characteristic of these methods is to not employ any pre-defined mesh: they use a set of nodes scattered within the problem domain as well as sets of nodes scattered on the boundaries of the domain. A particular attention is paid to the RPIM method, which proposes a nodal interpolation. The followed steps are: a presentation of « Meshfree methods » and their advantages compared to the traditional methods, an introduction to the meshfree RPIM method with a calibration of its associated parameters, and finally, a discussion on results obtained with the RPIM in forming processes exhibiting initiation and propagation of a crack showing the interest of the approach.

Depuis plus de 150 ans, le presse-garnitures (ou presse-étoupes) est le moyen le plus simple de réaliser l’étanchéité d’un robinet tout en permettant le mouvement de la tige pour ouvrir/fermer le robinet. C’est aussi la principale cause de fuites sur un robinet. Malgré l’utilisation de garnitures en graphite expansé qui réduit de manière significative le nombre de fuites, des problèmes comme l’augmentation de l’effort de manœuvre de la tige y ont été aussi associés. L’amélioration de l’étanchéité tout en réduisant l’effort de manœuvre de la tige est constamment recherchée sur un presse-garnitures de robinet afin d’optimiser industriellement. Pour atteindre cet objectif, le presse-garnitures est au cœur des études scientifiques et industrielles depuis 1957. Toutes les études identifiées dans la littérature ne permettent pas de maîtriser l’étanchéité et l’effort de manœuvre puisqu’elles ne prennent pas en compte le troisième corps présent dans le contact tige-garnitures lors des manœuvres de la tige. Les travaux de cette thèse proposent alors d’étudier le comportement tribologique du contact tige-garnitures. La démarche adoptée nécessite d’identifier les conditions de contact, ce qui est difficile par une instrumentation directe. Il faut donc reconstituer ces conditions par le couplage entre une approche expérimentale et une approche numérique. L’approche expérimentale est utilisée pour reproduire les conditions du contact tige-garnitures sur un banc d’essai simulant le fonctionnement d’un robinet à l’échelle 1:1. Afin d’améliorer la compréhension du comportement tribologique, des caractérisations tribologiques ont été réalisées sur la tige et les garnitures en graphite expansé après essai. Ainsi, elles ont permis de définir le circuit tribologique du contact tige-garnitures, qui contrôle à la fois l’étanchéité autour de la tige et l’effort de manœuvre. En effet, les débits du troisième corps contrôlent l’étanchéité tandis que le mécanisme d’accommodation de vitesse du troisième corps par cisaillement contrôle l’effort de manœuvre. Afin de compléter l’analyse expérimentale précédente, une approche numérique a été utilisée par l’intermédiaire d’un modèle par éléments finis. Le modèle permet d’évaluer la qualité d’étanchéité d’un presse-garnitures, vis-à-vis de la pression du fluide à étancher. De plus, l’influence du frottement aux interfaces tige-garnitures et boîte-garnitures sur l’étanchéité et l’effort de manœuvre est mise en évidence par l’intermédiaire d’une étude de sensibilité. Les travaux de cette thèse permettent l’appréciation de l’importance de prendre en compte l’influence du troisième corps dans un contact. Dans le cas d’un presse-garnitures, le troisième corps dans le contact tige-garnitures modifie des conditions de contact et celles-ci influent à leur tour l’étanchéité autour de la tige et l’effort de manœuvre. / For over 150 years, compression packing has been the simplest method of sealing a valve while allowing movement of the stem to move to open or close the valve. It is also the major source of valve leakage. Despite the use of flexible graphite packing, which reduces significantly the number of valve leakage incidents, it leads to other problems such as the increasing force required to cycle the stem. Improvement in quality of sealing while decreasing force required to cycle the stem is required in a valve compression packing in order to optimize the industrial system. To achieve this objective, compression packing has been a topic of scientific and industrial research since 1957. All published research does not allow to control the sealing quality and the force required the cycle the stem as it takes no account of the third body in the stem-packing contact while the stem moves. Thus, the work of this thesis designed to investigate the tribological behaviour of the stem-packing contact. The chosen approach requires identification of the contact conditions, which is practically impossible using local measurement. Hence, the contact conditions should be reconstructed using a coupled experimental and numerical investigation. The experimental investigation is carried out on a test rig, which simulates a full-scale operated valve to reproduce the conditions of the stem-packing contact. In order to improve the understanding of the tribological behaviour, the tribological characterization is made of the stem and the flexible graphite packing after the test. Hence, they helped to define the tribological circuit of the stem-packing contact, which controls both the sealing quality around the stem and the force required to cycle the stem. Indeed, the third body flow controls the sealing quality while the third body velocity shearing mode controls the force required to cycle the stem. In order to complete the experimental analysis, numerical investigation is carried out using a finite element model. The model allows to evaluate the sealing quality of a compression packing, with the regards to the fluid pressure. In addition, the influence of friction at the stem-packing and housing-packing interfaces is highlighted through a parametric analysis. The work in this thesis helps to raise awareness about the importance of taking into account the influence of the third body in a contact. In the case of compression packing, the third body changes the contact conditions which in turn, influence the sealing quality around the stem and the force required to cycle the stem

Le développement actuel de l’industrie vise à prévoir l’intégrité des structures dans le temps ou dans le cas de sollicitation extrême. Les risques liés à la propagation des fissures dans le cas de chocs ou d’impacts sont encore difficiles à prévoir. Les codes de calcul dans ce domaine regroupent plusieurs méthodes de simulation au sein d’un même code de calcul. Afin de présenter les différentes méthodes numériques mises en oeuvre, ce mémoire a été découpé en trois parties distinctes. Dans la première partie, nous présentons la bibliographie, puis notre apport aux méthodes de simulation numérique en l’appliquant au cas de la propagation de fissure dynamique et enfin les résultats obtenus à partir des méthodes proposées. Nous comparons ces simulations à des résultats expérimentaux ou à des simulations 2D trouvés dans la littérature. À travers la bibliographie, nous présenterons la théorie de la mécanique de la rupture pour arriver à un critère de propagation de fissure adapté à la dynamique transitoire. Ce critère a déjà été utilisé pour la fissuration dynamique en 2 dimensions. Nous décrirons la méthode des éléments finis étendus utilisée jusqu’ici principalement en quasi-statique. Nous donnerons les avantages mais aussi les limites de mise en oeuvre de cette méthode, notamment à travers le choix des enrichissements et de l’intégration des éléments coupés par la fissure. La méthode des level-sets est ensuite présentée : elle permet de décrire et faire évoluer la fissure indépendamment de la structure. On met en évidence le besoin de robustesse pour faire évoluer la fissure en dynamique explicite. La seconde partie est consacrée au développement et à l’extension de la méthode en 3D. Après avoir rappelé le critère de propagation en 3D fragile et avec plasticité, on cherche à proposer des schémas d’intégration spatiale plus économiques. Une nouvelle stratégie de propagation des level-sets basé sur la géométrie est proposée pour la dynamique explicite 3D. Enfin dans la troisième partie, nous appliquerons les méthodes à des cas de propagation de fissure bidimensionnelle puis tridimensionnelle. Nous simulerons dans un premier temps des cas 2D en mode I puis en mode mixte, afin de vérifier que l’on arrive à résultats proches des cas déjà simulés en 2D. Pour terminer par des simulations de propagation tridimensionnelle de fissure avec arrêt et redémarrage de la fissure. Tous ces développements on été implémentés dans le code de calcul de dynamique explicite EUROPLEXUS, co-propriété du CEA et de la Commission Européenne. / The current development of the industry focus on structural integrity over time or in the case of extremes stresses. Risks related to the cracks propagation in the event of shocks or impacts are still difficult to predict. Computing codes in this area groups several methods of simulation within the same computer software. To present the various numerical methods used, this thesis was divided into three distinct parts. In the first part we present the literature. Then, in second part, our contribution to the numerical simulation methods are presented by applying it to the case of dynamic crack propagation. Finally the results obtained from the proposed methods are described. We compare these simulations with experimental results or 2D simulations found in the literature. Through the first part, we present the theory of fracture mechanics to reach a criterion of crack propagation adapted to the transient dynamics. This criterion has been used for dynamic cracks in two dimensions. We describe the extended finite element method mainly used for quasi-static problems. We give the advantages but also the limits of this method: the choice of enrichment and the integration method are particularly important. The level-sets method is then presented: it allows to describe and develop the crack regardless of the structure. It highlights the need of robustness due to explicit dynamics scheme. The second part is devoted to the development and extension of the method in 3D. After reminding the propagation criterion in 3D, we try to offer more economic patterns of spatial integration. A new strategy of level-sets propagation based on geometrical approach is proposed for the explicit dynamic and applied in 3D. In the third part, we apply the methods to the case of two-dimensional crack propagation and three-dimensional. We initially simulate 2D mode I then mixed mode, to ensure that we arrive at results close to earlier 2D simulations. To finish, we present three-dimensional simulations of crack propagation with stopping and restarting crack. All these developments have been implemented in the computing software EUROPLEXUS , co-owned by the CEA and the European Commission.

Le travail présenté concerne le développement d’un cadre d’applications pour le traitement d’écoulements dispersés, tenant compte de l’effet de la coagulation sur la distribution des tailles des particules. Nous explorons également quelques techniques de validation expérimentale. Les modèles sont valables pour un écoulement incompressible et isotherme, avec des particules qui ont un temps de relaxation faible en comparaison à celui du fluide. Pour la phase dispersée, une méthode eulérienne est utilisée, ce qui permet d’extrapoler la vitesse des particules de celle du fluide. La distribution des tailles est modélisée à l’aide du « Direct Quadrature Method of Moments ». Cette approche permet de résoudre des équations de transport pour les poids et les abscisses d’une approximation de la distribution à l’aide des fonctions Dirac delta. L’effet de la coagulation est pris en compte à l’aide d’un noyau de collisions qui utilise la vitesse instantanée du fluide. Toutes les équations de transport sont résolues à l’aide de la méthode des éléments finis. Pour le fluide, les stabilisations « Streamline Upwind » et « Pressure Stabilized Petrov-Galerkin » sont utilisées ensemble avec une stabilisation grad-div. Afin de limiter le temps de calcul pour une simulation directe, une formulation utilisant un traitement explicite des termes d’advection est proposée. Avec l’apparition de gradients élevés, les équations de transport pour les particules nécessitent une stabilisation supplémentaire. Tout le travail est disponible dans le projet de logiciel libre Coolfluid 3, en utilisant un langage spécifique permettant une implémentation directe pour des modèles en éléments finis. Le code qui en résulte ressemble à la forme variationnelle des équations utilisées. Le programme est générique en termes de dimensions spatiales et de type d’éléments. Une première validation utilise des résultats trouvés dans la littérature comme référence. La précision des méthodes est vérifiée à l’aide des vortices Taylor-Green. Pour l’écoulement et les concentrations des particules, une simulation directe d’un canal turbulent est effectuée. Le noyau de coagulation est vérifié à l’aide de particules de différentes tailles qui tombent à travers un vortex Burgers. Finalement, quelques techniques de validation expérimentale sont utilisées dans une cellule d’essai. La technique « Particle Image Velocimetry » est utilisée pour les vitesses du fluide, tandis que la distribution des tailles est mesurée à l’aide du « Phase Doppler Anemometry » et « Multiple Wavelength Light Extinction ». La technique d’extinction de lumière est capable de produire des distributions des tailles qui peuvent être comparées facilement avec les résultats numériques. / The present work deals with the development of a framework for the modeling of dispersed flows, including the effect of coagulation on the particle size distribution. We also explore some techniques for experimental validation. Models are developed for incompressible, isothermal flow containing particles that have a small relaxation time compared to the fluid time scale. For the dispersed phase, an equilibrium Eulerian approach is used, extrapolating the particle velocity from the fluid velocity. The size distribution is modeled using the Direct Quadrature Method of Moments. In practice, this results in solving transport equations for the weights and abscissa of a Dirac delta approximation of the size distribution. To model the effect of coagulation, a collision kernel that makes use of the resolved instantaneous velocity is developed. All transport equations are solved using the Finite Element Method. For the fluid, the Streamline Upwind and Pressure Stabilized Petrov-Galerkin method are used, with additional grad-div stabilization. To decrease the solution time for DNS, a segregated formulation with an explicit advection term is proposed. The particle transport equations require cross-wind diffusion in addition to the streamline upwind stabilization when large gradients occur. All work is available in the open source Coolfluid 3 framework, using an Embedded Domain Specific Language we developed for the implementation of finite element models. The resulting code closely resembles the variational form of the equations and is generic in terms of element type and the number of spatial dimensions. A first validation uses literature results as reference. Correctness and accuracy of the methods are verified using the Taylor-Green vortex flow. For the fluid and particle concentration, direct numerical simulation of a turbulent channel flow is performed. The particle coagulation kernel is tested using particles of different sizes falling through a Burgers vortex. Finally, some experimental validation techniques are used on a small test chamber. Particle image velocimetry is used for the fluid motion, while the size distributions are measured using Phase Doppler Anemometry and Multiple Wavelength Light Extinction. The light extinction technique was found to produce size distributions that could provide valuable reference data for our particle model.

Les mécanismes produisant les signaux observés par la méthode de Polarisation Provoquée (PP) sur une large gamme de fréquences (entre 1 mHz et 20 kHz) ne sont pas encore complètement identifiés. Deux sujets particuliers ont été abordés dans cette thèse. L'origine du signal observé à haute fréquence (>1 kHz) a été analysé en effectuant des mesures de la résistivité complexe sur de l'eau à différentes salinités. Les résultats montrent des écarts sur la phase par rapport à la réponse attendue à haute fréquence. Ils dépendent du type d'électrode de mesures et de la conductivité du milieu. Un modèle basé sur un circuit électrique équivalent a été proposé pour modéliser ces effets. Nous avons aussi exploré le mécanisme responsable de la polarisation en présence de grains semi-conducteurs en analysant la dépendance de temps de relaxation. La réponse spectrale d'un milieu sableux saturé a été étudiée en variant la concentration et le type de l'électrolyte, la taille, le type et la quantité de grains semi-conducteurs insérés. En utilisant la méthode des éléments finis, nous avons entrepris une simulation numérique 2D basée sur la solution des équations de Poisson-Nernst-Planck dans le domaine temporel et spectral. Les résultats expérimentaux sont conformes à ceux issus de la simulation numérique et montrent une décroissance comparable du temps de relaxation avec l'augmentation de la concentration de l'électrolyte. Finalement, une campagne géophysique de terrain sur un site paléo-miniers contenant des grains semi-conducteurs complète l'approche de laboratoire. Des mesures de PP temporelles permettent de délimiter les zones de scories sur le site et d'en estimer le volume. / The physical mechanisms responsible for the induced polarization response over the frequency range (from 1 mHz to 20 kHz) are not completely understood. In particular, within the framework of this thesis, two subjects have been addressed. The origin of the signal observed at high frequency (HF) (>1 kHz) was analyzed by carrying out Spectral IP measurements on tap water samples. A phase deviation from the expected response has been observed at HF. The resulted deviation in phase appears to be dependent on the measuring electrode type (potential electrodes) and conductivity of the medium. A model based on an equivalent electrical circuit and designed to represent HF response, has been proposed to correct these effects. The mechanism responsible for the polarization in a medium containing semi-conductor grains has been investigated by analyzing the dependence of the relaxation time. We carried out experimental measurements on a sand medium containing different types of semi-conductors. The spectral response was studied by varying the concentration and type of the electrolyte, the size and content of semi-conductor grains. By using the finite element method, a 2D numerical simulation based on Poisson-Nernst-Planck equations was performed in time and frequency domains. The experimental results are qualitatively in accordance with numerical simulation. It showed a comparable decrease in the relaxation time when increasing the electrolyte concentration. Finally, field measurements on a paleo-mining site containing semi-conductor grains have been acquired. Time-domain IP measurements allowed us to define the zones of slag in the site and led to estimate the slag volume.

La corrosion localisée des alliages d'aluminium a pu être reliée qualitativement au rôle de la microstructure, des phases intermétalliques (IM) en particulier, qui provoquent la sensibilisation à la corrosion localisée des alliages d’aluminium contenant par exemple du cuivre. Mais un effort de modélisation est nécessaire si on veut aller vers de nouvelles méthodes d’essais.L’objectif de ce travail est d'étudier et de modéliser, sur des systèmes métallurgiques très simplifiés, les effets de couplage entre réactions électrochimiques et chimiques qui ne sont pas pris en compte dans l'analyse des essais conventionnels.L’endommagement lors de l’amorçage de la corrosion microstructurale a pu être simulé numériquement par résolution, par la méthode des éléments finis, de l’équation de Nernst-Planck décrivant les phénomènes de transport et les réactions des espèces chimiques en solution.Le comportement en milieu non chloruré de trois combinaisons où le couple IM (cathode)/matrice Al (anode) peut être assimilé à un couple Cuivre/Aluminium a été étudié, ceci de façon à rendre compte de la dissolution sélective qui fait que les phases S ou θ, présentes dans les alliages au cuivre, subissent un enrichissement en cuivre en surface. Pour valider le modèle « AlOH3 » on s’est donc appuyé sur la mesure de la vitesse de corrosion de l’anode pour trois combinaisons testées appelées systèmes modèles.En terme de mécanismes, la confrontation de l’ensemble de ces expériences aux résultats de la simulation numérique démontre que l’amorçage de la dissolution de l’anode nécessite un changement de pH local, et donc que la dépassivation locale de la surface d’aluminium par un effet chimique est l’élément déclenchant du couplage électrique. Pour aller vers une application industrielle et simuler successivement, l’amorçage et la propagation de la corrosion intergranulaire, il reste pour l’amorçage à intégrer les effets d’interaction entre phases et pour la propagation à construire à partir de ce modèle numérique, un modèle valable en milieu confiné et désaéré qui devra satisfaire aux règles physiques imposées par la dimension du système que constitue le joint de grain / In literature, localized corrosion of aluminum alloys has been qualitatively related to the role of microstructure, intermetallic phases (IM) in particular, which induces sensitization to localized corrosion of aluminum alloys containing copper, for example. But modeling effort is needed to go further in the development of new testing methods.The objective of this work is to study and model on very simplified metallurgical systems, the effects of coupling between chemical and electrochemical reactions that are not taken into account in the analysis of conventional tests.The damage induced by the microstructural corrosion initiation has been numerically simulated by resolution, by the finite element method, of the Nernst-Planck equation describing the mass transport phenomena and reactions of chemical species in solution.Behavior in non-chlorinated medium of three combinations where the couple IM (cathode) / Al matrix (anode) has been assumed to be close from a Copper / Aluminium couple has been studied, mimicking the selective dissolution of S or θ phases present in copper alloys leading to a surface copper enrichment. Validation of the model "AlOH3" defined in this study, was based on the measurement of the corrosion rate of the anode for the three model systems which were tested.In terms of mechanisms, the comparison of the experimental results of the numerical simulation clearly highlights that the initiation of the dissolution of the anode requires a change in the local pH, and therefore the local depassivation of the surface aluminum by a chemical effect is the triggering element of the electrical coupling.To transfer this numerical approach to an industrial application and simulate successively, the initiation and the propagation of intergranular corrosion, it remains, for initiation to introduce the effects of interaction between phases and for propagation to build, from this numerical model, a robust model in confined and deaerated environment which must comply with the rules imposed by the physical size of the system as defined by the grain boundary

Les enveloppes métalliques en acier non allié choisies par l’Andra comme surconteneurs pour le stockage des déchets radioactifs seront principalement soumises à la corrosion généralisée. L’objectif de ce travail est de démontrer qu’il est possible de simuler sur le long terme cette vitesse de corrosion de façon à bien dimensionner les surconteneurs pour qu’ils restent étanches et confinent la radioactivité des déchets qu’ils contiennent pendant plusieurs siècles.La modélisation par la méthode des éléments finis, basée sur la résolution numérique de l’équation de Nernst-Planck en potentiel libre a donc été utilisée pour prévoir le comportement électrochimique de ces aciers sur des durées inaccessibles à des expériences en laboratoire. A partir d’études paramétriques (influence des constantes de cinétique électrochimique, des constantes cinétiques associées aux réactions chimiques, de la température) et en se basant sur plusieurs comparaisons entre expériences de laboratoire et simulations numériques, il a été possible de déterminer l’ensemble des paramètres et des règles de fonctionnement nécessaires pour construire et utiliser un modèle numérique prédictif le plus complet possible.Dans la simulation à long terme, qui constitue le résultat majeur du travail, on constate qu’après épuisement de l’oxygène dissous initialement piégé, la vitesse de corrosion est contrôlée par le recouvrement de la surface métallique dû à l’apparition de produits de corrosion (magnétite, sidérite). Ce phénomène induit une diminution de la porosité totale, identifiée comme déterminante dans le ralentissement de la vitesse de corrosion. / Unalloyed steels selected by Andra for nuclear waste storage would be mainly affected by general corrosion. The aim of this study is to demonstrate that it is possible to simulate the long term corrosion rate in order to estimate the thickness of the containers expected to maintain the confinement for several centuries. Modeling by the finite element method, based on the resolution of Nernst-Planck equation in free potential conditions has been used to predict the electrochemical behavior of these steels impossible to reach at long time through laboratory experiments. From parametric studies (influence of electrochemical kinetic constants, kinetic constants dealing with chemical reactions, temperature) and in the light of several comparisons between laboratory and numerical experiments, it was possible to determine all the parameters and all the rules to build and use the most completed predictive numerical model.The main result of this study is a long term prediction model. After a first step it was established that dissolved oxygen initially trapped is consumed, whereas the corrosion rate is controlled by the fractional coverage of the surface due to the formation of corrosion products (magnetite, siderite). As a consequence, the decrease of the total porosity can be identified as a key parameter in the reduction of the corrosion rate.

La méthode du contour est une technique récente pour la mesure de contraintes résiduelles dans des structures mécaniques. Elle permet d'obtenir un champ de contraintes résiduelles en 2D. Dans la méthode du contour, l'échantillon est découpé en deux parties. Les contraintes résiduelles normales aux faces découpées se détendent, ce qui cause leur déformation. Le profil de déformation des faces coupées est mesuré, puis utilisé pour le calcul de contraintes résiduelles. La méthode du contour est basée sur le principe de superposition élastique. Sa théorie assume que la détente de contraintes résiduelles est entièrement élastique. Cependant, en pratique, la détente élastique de contraintes résiduelles peut être accompagnée d'une détente plastique, ce qui induit des erreurs dans les contraintes mesurées. Ce projet se concentre sur l'étude et le contrôle des erreurs de plasticité dans la méthode du contour. Les résultats de cette recherche fournissent d'importantes informations sur la façon dont les déformations plastiques affectent les résultats de la méthode du contour. Par ailleurs une nouvelle technique de coupe pour l'atténuation des erreurs de plasticité a été développée. Des outils de calcul et des instructions sont proposés pour l'estimation des erreurs de plasticité. Enfin des instructions pour la minimisation des erreurs de plasticité sont proposées et appliquées à la mesure de contraintes résiduelles dans un échantillon. / The contour method has emerged as a promising technique for residual stress measurement in relatively large, thick and complex engineering components. The method involves making a cut in the sample of interest, measuring the subsequent relaxed deformation of the cut face and using this profile to back-calculate the original residual stress field by FE modelling. The method is based on the theory of elasticity in that the stress relaxation during test specimen cutting is assumed to be entirely elastic. However, when measuring residual stresses close to the material yield stress, plasticity can occur and affect the measurements.The main aim of this thesis was to develop methods of mitigating and estimating plasticity-induced errors in the contour method. The outcomes of this research provide a valuable insight into how accumulation of plasticity affects the performance of the contour method. A novel cutting strategy that aims at mitigating the plasticity-induced errors has been developed. Furthermore, procedures are developed to estimate the plasticity-induced errors. Finally guidelines are proposed and applied to a case study for mitigating the plasticity-induced errors in the contour method.