La démarche actuelle des industries aérospatiales est de diminuer le coût de lancement des engins spatiaux par une réduction de la masse des composants. Dans l’optique de cette démarche, de nouveaux matériaux sont élaborés et permettent de satisfaire aux exigences de densification, de dissipation thermique et de réduction de masse des équipements électroniques embarqués dans les satellites. Cette thèse est une contribution à l’étude de deux de ces matériaux, l’AlSi CE9F et une nuance d’alumine cofrittée à température, destinés à réaliser des boitiers hybrides de protections de composants électroniques, initialement conçus en Kovar. Les objectifs sont d’affiner les connaissances sur propriétés mécaniques des deux matériaux et de mettre en place des règles de conceptions propres à leurs comportements mécaniques. En effet, l’AlSi CE9F et l’alumine ont un comportement à rupture fragile. La détermination de leurs résistances à la rupture est alors réalisée dans le cadre de la théorie de Weibull. Des séries d’essais de flexion quatre points et trois points sont effectués. Elles permettent d’identifier les paramètres de Weibull des deux matériaux à température ambiante et de mettre en évidence les effets de volume. L’étude expérimentale est poursuivie sur l’AlSi CE9F afin de déterminer l’influence de la température sur ses propriétés mécaniques à travers deux approches. La première s’intéresse à une variation monotone de la température et la seconde à des cycles thermiques entre -50 et 125°C. Si la première étude ne montre qu’une faible évolution du module d’élasticité, la seconde démontre que les cycles thermiques contribuent à l’amélioration de la résistance à la rupture de l’AlSi CE9F. Cette augmentation de la contrainte à la rupture se traduit également par une évolution de sa microstructure. Dans un second temps, un modèle de Weibull est numériquement mis en place à partir des paramètres identifiés et du critère de la contrainte équivalente de Freudenthal. Ce critère est analysé et validé à travers l’étude de trois éprouvettes en AlSi CE9F à chargements complexes. Le modèle validé est enfin utilisé pour décrire le comportement mécanique de deux composants dans différentes configurations de sollicitation, réalisés respectivement en alumine HTCC et en AlSi CE9F. Une méthodologie de dimensionnement est alors mise en place et permettra de disposer de nouvelles règles de conception équivalentes à celles existant sur les matériaux classiques. / The aerospace companies currently want to decrease the price of spacecraft launching with a reduction of the mass components. New materials were recently developed to satisfy the rising requirements of thermal dissipation, densification and weight decrease of on-board electronic equipment intended to satellite. This thesis is a contribution to the characterization of two of these innovative materials: AlSi CE9F and a grade of alumina HTCC. These materials are designed to manufacture hybrid boxes for computing chips, originally made in Kovar. The objectives are to improve the mechanical properties knowledge of these materials and to develop a know-how design specific to their mechanical behaviours. Indeed, AlSi CE9F and alumina have brittle fracture behaviour. The strength analysis is also realized in connection with the Weibull theory. The Weibull’s parameters are identified from the four points and three points bending strength and the volume effects are highlighted. The experimental study is completed by the analysis of the temperature influence on the mechanical properties of AlSi CE9F through two approaches. The first one considers a monotonic variation of temperature and shows a minor evolution of the elastic modulus. The second one proves that thermal cycles between -50 and 125°C improve the strength value of AlSi CE9F. This increase is also reflected by an evolution of its microstructure. Secondly, a Weibull’s model is numerically established based on identified parameters and the Freudenthal’s equivalent stress criterion. The Freudenthal’s criterion is analysed and confirmed through the study of complex loading samples made in AlSi CE9F. The confirmed model is finally used to describe the mechanical behaviours of two components respectively made in AlSi CE9F and alumina HTCC, thoroughly in several loading configurations. A design methodology is developed and will bring new rules in modelling and design, closed to those existing in conventional materials.

Le projet ANR ICARE, dans lequel s'inscrit cette thèse, vise au développement de méthodes pour l'analyse de structures complexes et de grande taille. Le défi scientifique consiste à investiguer des zones très localisées, mais potentiellement critiques vis-à-vis de la tenue mécanique d'ensemble. Classiquement, sont mis en œuvre aux échelles globale et locale des représentations, discrétisations, modèles de comportement et outils numériques adaptés à des besoins de simulation gradués en complexité. Le problème global est traité avec un code généraliste dans le cadre d'idéalisations topologiques (formulation plaque, simplification géométrique) et comportementale (homogénéisation) ; l'analyse locale quant à elle demande la mise en œuvre d'outils spécialisés (routines, codes dédiés) pour une représentation fidèle de la géométrie et du comportement.L'objectif de cette thèse consiste à développer un outil efficace de couplage non-intrusif pour la simulation multi-échelles / multi-modèles en calcul de structures. Les contraintes de non-intrusivité se traduisent par la non modification de l'opérateur de rigidité, de la connectivité et du solveur du modèle global, ce qui permet de travailler dans un environnement logiciel fermé. Dans un premier temps, on propose une étude détaillée de l'algorithme de couplage global/local non-intrusif. Sur la base d'exemples et de cas-test représentatifs en calcul de structures (fissuration, plasticité, contact...), on démontre l'efficacité et la flexibilité d'un tel couplage. Aussi, une analyse comparative de plusieurs outils d'optimisation de l'algorithme est menée, et le cas de patchs multiples en interaction est traité. Ensuite le concept de couplage non-intrusif est étendu au cas de non-linéarités globales, et une méthode de calcul parallèle par décomposition de domaine avec relocalisation non-linéaire est développée. Cette méthode nous a permis de paralléliser un code industriel séquentiel sur un mésocentre de calcul intensif. Enfin, on applique la méthode de couplage au raffinement de maillage par patchs d'éléments finis. On propose un estimateur d'erreur en résidu explicite adapté au calcul de solutions multi-échelles via l'algorithme de couplage. Puis, sur la base de cet estimateur, on met en œuvre une procédure non-intrusive de raffinement local de maillage. Au travers de ces travaux, un outil logiciel de couplage non-intrusif a été mis au point, basé sur l'échange de données entre différents codes de calcul (protocole Message Passing Interface). Les développements effectués sont intégrés dans une surcouche Python, dont le rôle est de coupler plusieurs instances de Code_Aster, le code d'analyse de structures développé par EDF R&D, lequel sera utilisé dans l'ensemble des travaux présentés. / This PhD thesis, part of the ANR ICARE project, aims at developing methods for complex analysis of large scale structures. The scientific challenge is to investigate very localised areas, but potentially critical as of mechanical systems resilience. Classically, representation models, discretizations, mechanical behaviour models and numerical tools are used at both global and local scales for simulation needs of graduated complexity. Global problem is handled by a generic code with topology (plate formulation, geometric approximation...) and behaviour (homogenization) simplifications while local analysis needs implementation of specialized tools (routines, dedicated codes) for an accurate representation of the geometry and behaviour. The main goal of this thesis is to develop an efficient non-intrusive coupling tool for multi-scale and multi-model structural analysis. Constraints of non-intrusiveness result in the non-modification of the stiffness operator, connectivity and the global model solver, allowing to work in a closed source software environment. First, we provide a detailed study of global/local non-intrusive coupling algorithm. Making use of several relevant examples (cracking, elastic-plastic behaviour, contact...), we show the efficiency and the flexibility of such coupling method. A comparative analysis of several optimisation tools is also carried on, and the interacting multiple patchs situation is handled. Then, non-intrusive coupling is extended to globally non-linear cases, and a domain decomposition method with non-linear relocalization is proposed. Such methods allowed us to run a parallel computation using only sequential software, on a high performance computing cluster. Finally, we apply the coupling algorithm to mesh refinement with patches of finite elements. We develop an explicit residual based error estimator suitable for multi-scale solutions arising from the non-intrusive coupling, and apply it inside an error driven local mesh refinement procedure. Through this work, a software tool for non-intrusive coupling was developed, based on data exchange between codes (Message Passing Interface protocol). Developments are integrated into a Python wrapper, whose role is to connect several instances of Code_Aster, the structural analysis code developed by EDF R&D, which will be used in the following work.

Le procédé de rétreint fait partie de la famille des techniques de mise en forme sans enlèvement de matière. La déformation du lopin est obtenue par chocs successifs d'un ensemble de matrices disposées autour de la pièce. Ce procédé est généralement utilisé pour la réduction de section de tubes ou de barres. Dans un contexte industriel, la maîtrise des paramètres procédé et la compréhension des phénomènes sous-jacents est indispensable pour non seulement limiter le temps de développement de nouveaux produits mais aussi diminuer le nombre de rebus des références qui posent problème. L'objectif de la thèse est de répondre à ces attentes en exploitant les possibilités offertes par la simulation numérique. Dans un premier temps, les alliages étudiés ont été caractérisés mécaniquement. En effet grâce à des essais de traction et d'impact de Taylor et à l'emploi d'une méthode d'identification par analyse inverse, les paramètres de la loi d'écoulement de Johnson-Cook ont été déterminés. Dès lors différents modèles, utilisant le code Abaqus/Explicit, ont été mis en place (du 2D axisymétrique au 3D en passant des modèles en reprenant que le quart de la pièce) et la validation s'est faite grâce à des campagnes d'essais réalisés sur site. Enfin, un progiciel développé en C++ sera livré à l'industrie. Il intègre différentes fonctionnalités comme la prédiction des efforts de forge par calcul analytique. Tous ces outils numériques et analytiques ont permis de mieux comprendre le procédé en termes de chemin d'écoulement de la matière, de distribution des contraintes et déformations et de profil d'évolution d'évolution de l'effort tout au long de la mise en forme. ABSTRACT : The shrinking process is part of the family of formatting techniques without removing material. / The shrinking process is part of the family of formatting techniques without removing material.

Ce travail de thèse est consacré à l'étude des milieux granulaires humides. On étudie principalement le cas de faibles teneurs en eau. L'influence du phénomène de capillarité sur la rétention d'eau et sur le comportement mécanique du milieu est analysée. A chaque état d'équilibre, la distribution d'eau dans le milieu granulaire est obtenue en appliquant une succion homogène dans un volume élémentaire représentatif. La méthode multi-échelles utilisée se nourrit d'expérimentations/modélisation à l'échelle des interactions entre les grains, et produit des simulations en éléments discrets du comportement macroscopique qui sont comparées aux résultats d'expérimentation. A l'échelle des interactions capillaires, une étude expérimentale du pont liquide pilotée en succion permet la validation d'un modèle basé sur l'approximation toroïdale du profil du pont liquide. A l'échelle macroscopique, les courbes de rétention d'eau simulées sont proches des courbes de rétention d'eau expérimentales réalisées sur des milieux modèles composés de billes de verre. Enfin, le comportement mécanique macroscopique est simulé. On note une fragilisation du matériau en fonction de la succion dans le domaines des faibles teneurs en eau qui s'explique par une diminution de la densité des ponts liquides. La prise en compte d'une rugosité des grains permet une meilleure description de la transition entre un état humide et un état sec. L'analyse des contraintes montre la pertinence du tenseur des contraintes associé aux interactions attractives dans l'étude de l'évolution de la contrainte à la rupture (compression simple) et de la compressibilité (compression oedométrique) en fonction de la succion. / This work deals with study of humid granular media. Weak water percent media are mainly studied. The influence of capillarity on water retention and mechanical behaviour is analyzed. At each state of equilibrium, water distribution is obtained by applying a homogeneous suction in the entire elementary representative volume. The multi-scale approach used, is based on experimental/modelling at the local scale (interaction between grains) and produce discrete elements simulation of the macroscopic behaviour which are compared to experimental results. At the scale of the interactions, an experimental study of the liquid bridge piloted by suction validates a model based on the toroidal approximation of the liquid bridge shape. At the macroscopic scale, simulated water retention curves are near of experimental water retention curves made on a model sample composed of glass beads. Then, the mechanical behaviour is simulated. We notice a loss of cohesi on in function of suction because of a reduction of the liquid bridge density. The introduction of a roughness allows a best representation of the transition between the non saturation state and the dry state. The analysis of the stress reveal the pertinence of the stress tensor associated with attractives interactions in the evolution of rupture stress (simple compression) and compressibility (oedometric compression) in function of suction.

L'évaluation des risques naturels et du risque sismique en particulier a pris une importance croissante pour les autorités publiques ces dernières années, entraînant une adaptation importante des outils utilisés jusqu'à présent dans l'étude de la vulnérabilité du bâti. Il est nécessaire de mieux évaluer ces risques afin de mieux y faire face. De nouvelles approches doivent être proposées pour répondre à ces objectifs dont les approches numériques font partie. Parmi les méthodes numériques développées pour tenir compte de la spécificité du bâti maçonné figurent les approches par éléments discrets comme la méthode NSCD (« Non Smooth Contact Dynamics »). Les travaux présentés dans cette thèse visent à mieux appréhender le comportement mécanique et dynamique des ouvrages maçonnés en utilisant ce type de méthode implémenté dans le code calcul LMGC90, afin d'affiner la compréhension et la caractérisation de la vulnérabilité sismique des édifices en maçonnerie. La méthodologie que nous cherchons à mettre en place vise à traiter la vulnérabilité sismique de bâtiments spécifiques ou au contraire d'ensembles de bâtiments de même typologie. / The perception and the evaluation of natural risks and more particularly of the seismic hazard have never been so important and highlighted by the authorities. This phenomenon leads to a deep adjustment of tools used up to now in the assessment of the vulnerability of buildings. New approaches as numerical approaches have to be proposed so as to best estimate the risk. Among all the numeriacl methods developed to take into account specificities of masonry buildings, approaches by discrete elements methods as the NSCD method (“Non Smooth Contact Dynamics”) seem to be particularly well adapted. This PhD thesis aims to best understand the mechanical behaviour of masonry structures by using this discrete method implemented in the software LMGC90 so as to give an accurate estmation of the seismic vulnerability of buildings. The very purpose of this work is to develop a calculation tool able to deal with the vulnerabilityof both specific and common masonry buildings.

Dans le contexte industriel actuel, l'utilisation du modèle de forces cohésives en milieu solide revêt d'une importance capitale dans des nombreux domaines d'applications dont en particulier le domaine du transport, du nucléaire ou du génie-civil. En conséquence, la présente étude propose une modélisation numérique de la rupture fragile d'une plaque en acier en présence des forces cohésives. Le logiciel Cast3M est utilisé pour la mise en œuvre numérique. Le modèle de type Dugdale est utilisé et permet de calculer dans le cas d'une plaque élasto-plastique la distribution de contraintes à la pointe de fissure. Le modèle suppose, qu'à la pointe de fissure sur une longueur cohésive Lzc, les contraintes cohésives σc (x1), suivant l'axe Ox1 restent constante avec une amplitude égale à σc. Le domaine de rupture élastique et le domaine de ruine plastique a été prisent en compte. L'approche peut être considérée comme une extension de la mécanique linéaire de ruptures dans laquelle est introduite la notion de contrainte critique, ce qui permet en particulier de mieux rendre compte des effets d'échelle. A partir du principe de superposition des contraintes en mode I de rupture, les calcules dans la configuration déformée de la plaque fissurée a été étendue à d'autres géométries plus simples : une plaque avec fissure sans force cohésive, mais avec charge extérieure appliquée, une plaque avec une fissure purement cohésive et une plaque avec une fissure partiellement cohésive et sans charge extérieure. Ces différentes formes géométriques ont été prises en compte dans la modélisation numérique ainsi que la création de la fissure et le démarrage de la propagation. Le paramètre énergétique calculé selon la formule de l'intégrale J de Rice (1968) est également prise en compte et est comparé avec les résultats analytiques de Ferdjani et al. (2007). / In the context of industrial applications, the linear fracture mechanics theory is not sufficient to account for various aspects of the crack propagation and it becomes necessary to use more sophisticated models like cohesive force models. The goal of the present thesis is to develop such a model in order to account for all the process of the crack propagation, from the nucleation in a sound structure to the final failure of that structure. Specifically, we use Dugdale model which contains both the concept of critical stress and of internal length. The presence of a critical stress allows us to account for the nuckleation of a crack even if the body is initially sound, in contrast with Griffith theory which require the existence of preexisting cracks. The presence of an internal length allows us to account for pertinent scale effects, once again in contrast with Griffith theory in which the predicted scale effects are only correct for large specimen, not for small specimen. This approach is used to describe the full process of cracking in several cases, either by purely analytical methods when the geometry is simple enough, or by the finite element method with the code Cast3M in the case of complex geometry.

Ce travail de thèse concerne la prédiction des structures en béton soumises à des chargements extrêmes. Il s'intéresse plus particulièrement au comportement du béton sous fort confinement où la contrainte peut atteindre des niveaux de l'ordre du giga Pascal. La modélisation de ce comportement doit être capable de reproduire la compaction irréversible. Pour ce faire, deux lois de comportement élasto-plastique - endommageable ont été développées et implantées dans un code aux éléments discrets. Les paramètres utilisés dans ces lois sont calibrés par les simulations des essais de traction/compression uniaxial, des essais hydrostatiques et triaxiaux. Une fois les paramètres calibrés, la loi montrant le meilleur agrément avec l'expérience a été choisie pour la prédiction de la réponse du béton sous différents niveaux de confinement. Les résultats du modèle sont analysés non seulement à l'échelle macroscopique mais également à l'échelle de l'élément discret. La nécessité de prendre en compte une loi d'interaction de type élasto-plastique-endommageable est aussi montrée. La deuxième partie du travail de thèse développe une méthode de couplage entre le modèle éléments discrets et un modèle d'écoulement compressible en tenant compte des mécanismes physiques fondamentaux d'interaction entre l'écoulement interne et les particules solides d'un matériau poreux. Le problème d'écoulement est résolu par une méthode en volumes finis, où le volume est discrétisé en tétraèdres issus d'une triangulation régulière de Delaunay. Notre modèle est une adaptation aux fluides compressibles d'un modèle développé initialement pour les écoulements incompressibles. Ce couplage a été utilisé pour simuler le comportement triaxial des bétons humides et saturés sous différents niveaux de confinement. Les résultats nous montrent une bonne reproduction du comportement non-drainé du béton saturé sous faible confinement. Pour fort confinement, les simulations ne se rapprochent des résultats expérimentaux qu'au prix d'une compressibilité du fluide plus faible que celle de l'eau. Par ailleurs, la contrainte effective était une variable pertinente pour décrire le comportement du béton humide par un état limite intrinsèque indépendant du degré de saturation. / This thesis work deals with the predicting of concrete structures submitted to some extreme loadings, and, more particularly, focuses on behavior of concrete under a high-confining pressure. At this range of pressures, irreversible compaction of the material occurs and needs to be considered. Doing so, two elasto-plastic-damaged constitutive laws have been developed and implanted into a discrete element numerical code. Local parameters to be used in these constitutive laws are identified by simulating reference uniaxial traction/compression tests and triaxial compression tests. Once these parameters have been obtained, the law showing the best agreement with the experimental data has been chosen to predict the reponse of concrete sample for triaxial compressive tests at different levels of confinement. The numerical results have been analyzed not only at macroscopic scale but also at discrete element scale. The need of a constitutive law taking into account the elasto-plastic-damaged behavior has been also proved. The second objective of the thesis work was to develop a fluid flow – coupled discrete element model by considering fundamental physical mechanisms of the interaction between the internal fluide flow and the solid particles of a porous material. The flow problem is solved by the finite volume method, where the volume is discretized into tetrahedra issue of a regular Delaunay triangulation. Our model is an adaptation for elastic fluids of a model originally developed for incompressible flows. The developed fluid-flow coupled discrete element has been used to simulate the undrained triaxial behavior of concrete under different levels of confinement. The results show a good reproduction of undrained behavior of saturated concrete under low confinement. For high confinement, the simulations only resemble the experimental results when the fluid compressibility is lower than that of water. Moreover, the effective stress was a relevant variable to describe the behavior of the wet concrete by an intrinsic limit state independent of the degree of saturation.

Dans une problématique d'analyse de la durabilité des structures en béton armé, la quantification de la localisation des déformations et des propriétés des fissures sont deux points clés. Ce travail présente une méthode permettant, dans le cadre de la mécanique des milieux continus, d'améliorer la description de l'évolution de la localisation des déformations. En se basant sur une approche continue du problème, l'évolution des nonlinéarités dans le béton est décrite au travers d'un modèle d'endommagement régularisé. Pour améliorer la description de la localisation des déformations, une modification de la méthode de régularisation nonlocale intégrale sur les variables internes est proposée. L'influence de l'état de contrainte sur les interactions nonlocales est introduite dans la régularisation, afin de prendre en compte la dégradation de la structure ainsi que l'influence des conditions aux limites sur les interactions nonlocales. Cette méthode, implantée dans le code aux éléments finis Cast3M, est validée sur différents cas tests analysant l'évolution des nonlinéarités de l'enclenchement de l'endommagement jusqu'à la rupture et permet notamment de résoudre des pathologies identifiées pour la méthode nonlocale originale. La comparaison avec des résultats expérimentaux montre également la capacité du modèle à décrire l'évolution de la fissuration dans une structure. Le modèle développé est ensuite utilisé pour analyser le comportement de structures en béton armé et sert de base pour introduire une description de la fissuration dans une modélisation simplifiée de type poutre multifibre. A partir de calcul 3D sur des éléments en béton armé utilisant le modèle développé, une loi uniaxiale est identifiée pour déterminer la fissuration dans une fibre en fonction de l'énergie dissipée par le modèle d'endommagement. Une comparaison avec des résultats expérimentaux est effectuée et montre la capacité de cette approche simplifiée à estimer la fissuration. / For the durability analysis of reinforced concrete structures, the modelling of strain localization and the estimation of cracking properties are hot topics. This work introduces a method allowing, in the framework of continuous mechanics, to improve the description of the evolution of strain localization. Based on a continuous description of the problem, the evolution of nonlinearities in concrete is described with a regularized damage model. In order to improve the description of strain localization, a modification of the nonlocal integral regularization method is proposed. The influence of the stress state on the nonlocal interactions is introduced in the regularization method, in order to take into account the degradation of the structure (decrease of the bearing capacities) as well as the influence of free boundary conditions. This method, implemented in the finite element code Cast3M, is validated against several cases of study, by analyzing the evolution of nonlinearities from damage initiation up to failure. It allows solving several pathologies pointed out for the original nonlocal method. The comparison with experimental results shows also the capacity of the proposed model to describe the evolution of cracking in a structure. Then, the model developed is used to analyse the behaviour of reinforced concrete structures and to develop a method to quantify cracking in a multifiber beam element modelling. From 3D calculation on reinforced concrete element with the new nonlocal model developed, a uniaxial law is identified in order to estimate cracking as a function of the energy dissipated by the damage model. A comparison with experimental data is performed and shows the potentiality of this simplified approach to estimate cracking.

Le galetage est une opération de finition et de traitement mécanique de surface, particulièrement utilisée dans le secteur du décolletage. Afin de valoriser pleinement l'opération, les industriels ont besoin d'estimer les caractéristiques des pièces traitées; il s'agit de la rugosité, de la dureté superficielle et de la distribution des contraintes résiduelles. Pour répondre à cette problématique, différentes modélisations ont été construites par le passé. Les deux plus abouties s'appuient principalement sur l'hypothèse d'équivalence entre la sollicitation de galetage et la sollicitation d'indentation. Dans ce rapport, une nouvelle approche de la modélisation de l'opération de galetage est proposée. Elle repose sur une étude fine des phénomènes physiques mis en jeu par la sollicitation de galetage. Dans un premier temps, le problème mécanique posé par la modélisation de l'opération de galetage est résolu sous l'hypothèse de comportement purement élastique. Dans ce cas, l'équivalence entre la sollicitation de galetage et la sollicitation d'indentation est bien re-démontrée. Dans un second temps, l'effet de la plasticité sur cette équivalence est étudié. Un phénomène de formation et écoulement de bourrelet, particulièrement influent sur la sollicitation de galetage, est alors mis en évidence. Ce constat invalide totalement l'hypothèse d'équivalence galetage - indentation formulée précédemment. Ainsi, il est démontré que pour modéliser fidèlement les phénomènes induits par l'opération, il est indispensable de considérer le contact roulant libre du galet sur un solide en trois dimensions. En tirant partie de cette étude, une nouvelle méthode de modélisation de l'opération de galetage par la technique des Éléments Finis est alors proposée. Une attention particulière est portée sur la définition d'une méthode de caractérisation mécanique préalable du matériau, notamment adaptée à la problématique par l'utilisation d'un essai de compression spécifique. En ce qui concerne l'état géométrique de la surface, les résultats numériques sont en bon accord avec les mesures expérimentales. De plus, pour les grandeurs mécaniques (dureté superficielle et contraintes résiduelles), une correspondance qualitative est obtenue. Finalement, il apparaît que ce modèle apporte une contribution notable à la compréhension de l'opération, ouvrant la voie à une modélisation plus fidèle. / Roller burnishing is a cost effective surface enhancement process. To value this process, industry needs a tool to estimate the roughness, the hardness and the residual stress distribution of a burnished workpiece. To achieve this goal, different models have been built in the past. The two most comprehensive are based on the hypothesis of equivalence between normal contact and rolling contact. In this study, a new approach of the modeling of roller burnishing is proposed. For this, physical phenomena that occur around the contact area have been precisely studied. Firstly, the mechanical problem that is posed is solved analytically assuming a perfectly elastic behavior. In this case, equivalence between normal contact and rolling contact is proved. Secondly, the effect of plasticity on that result is studied. A phenomenon of accumulation and flow of material, which effect appears to be particularly important on the mechanical and geometrical characteristics of the burnished surface is highlighted. This observation imposes to reject the hypothesis of equivalence between normal contact and rolling contact. Consequently, it is conclude that any model can't be reliable until it considers the rolling contact of the roller on the workpiece. On this basis, a new finite element model of roller burnishing is established. Finally, numerical results are compared to experimental measurements. Concerning geometrical state, a quantitatively good correlation is observed while a qualitatively agreement is obtained for the two mechanical indicators. It can be concluded that this new EF model give a better understanding of the mechanics of roller burnishing process and will give the opportunity to improve parameters process or to adapt parameters according to the wished mechanical and geometrical characteristics of the workpiece.

En Guinée, la disponibilité des ressources en eau est soumise à de graves problèmes de gestion avec des conséquences importantes pour les populations, notamment à Kankan où la plupart des puits tarissent pendant la saison sèche, et où la fourniture régulière de l’eau à partir de la rivière Milo n'est pas assurée dans tous les quartiers de la ville. Le Milo et son bassin versant représente donc un enjeu important pour les habitants de cette ville. Le Milo (480 km de long) qui draine un bassin versant 13 810 km2 est le plus important des affluents rive droite du Niger. Pour évaluer l’impact du bassin versant du Milo et de la ville de Kankan sur la qualité des eaux et des sédiments de ce cours d’eau, 30 échantillons de sédiments de fond ont été prélevés au cours d’une première campagne (juin-juillet 2013) sur l’ensemble des cours d’eau du bassin versant, tandis que durant une deuxième campagne (avril 2014-mai 2015), un suivi régulier a été mis en place sur le Milo en amont et en aval de Kankan. Au cours de cette dernière, 232 échantillons d’eau du Milo (dont 116 utilisés pour l’isotopie), et 26 échantillons de sédiments de fond ont été prélevés sur les stations de Bordo (amont de Kankan) et de Karifamoriah (aval de Kankan). Les concentrations en éléments majeurs et traces, terres rares, COP, COD, pH, alcalinité, l’azote organique, les isotopes de l’oxygène et de l’hydrogène, ainsi que la micro-granulométrie des sédiments ont été mesurés. Les résultats obtenus durant le cycle hydrologique 2014-2015 ont permis d’estimer le flux de matières exportées en solution par le Milo à Kankan à 47 863 t.an-1 et le flux de matières particulaires à 76 759 t.an-1. 92% de ce tonnage est exporté durant la période des hautes eaux. Ces flux permettent d’estimer des vitesses moyennes d’érosion physique (8 t.km-2.an-1) et d’altération chimique (5 t.km-2.an-1) du bassin relativement faibles. Le flux spécifique de CO2 consommé par l’altération chimique est lui aussi relativement faible (76.103 mole.km-2.an-1), mais il reste dans la moyenne des flux de CO2 mesurés sur les bassins couverts de sols latéritiques. Le degré de contamination des sédiments et sa variation spatio-temporelle ont été évalués grâce au calcul du facteur d’enrichissement (FE), en normalisant les concentrations en éléments traces par rapport à Al, Sc et Ti et en prenant le PAAS et l’UCC comme matériaux de référence. La quasi-totalité des éléments traces proviennent essentiellement des processus d’altération et sont d’origine naturelle. Toutefois, un enrichissement modéré à significatif (5 < FE < 20) a été mis en évidence pour Zr et Hf (enrichissements naturels) et As et Sb (impacts anthropiques). Finalement, l’apport anthropique en éléments traces, quand il existe, reste modéré. Les profils de concentrations en terres rares ont permis de montrer que la signature géochimique des sédiments de fond du Milo est homogène sur l’ensemble du bassin et caractéristique des sols latéritiques, non perturbée par la ville de Kankan et elle est proche des signatures du PAAS et de l’UCC. Pour l’ensemble des éléments traces et des terres rares la contribution anthropique reste faible et les fractions disponibles (phases labiles) sur les sédiments (extraction à l’EDTA) sont également faibles, confirmant ainsi que l’essentiel de ces éléments rentrent principalement dans la composition des phases résiduelles.