Au cours du temps, les cavités souterraines sont soumises à un vieillissement et plusieurs types de dégradation peuvent apparaitre. Les anciennes exploitations souterraines, parfois constituées d’un ou plusieurs niveaux, n’ont, sans doute, pas été conçues pour être stables à long terme. Elles ont été réalisées à une époque où n'existaient pas d'enjeux en surface, de zones de travaux, ce qui permettait d’éviter de se préoccuper des mouvements de sol induits. Elles ont pu quelquefois être totalement ou partiellement remblayées, mais pas de manière systématique. L’effondrement d'une cavité souterraine engendre la déconsolidation des niveaux supérieurs des terrains de recouvrement. Ces mécanismes peuvent provoquer en surface deux types de désordres : un affaissement ou un fontis. L'affaissement et le fontis peuvent provoquer des graves dommages aux structures et aux infrastructures en surface mais aussi mettre en péril la sécurité des populations. Les travaux de cette thèse s’articulent autour de l’étude du foisonnement, du fontis et de sa propagation dans les terrains de recouvrement.Les objectifs de cette thèse sont doubles : tout d’abord il s’agit d’étudier le foisonnement de la roche lors d’un effondrement de toit de carrières souterraines ; ensuite il s’agit de modéliser la propagation du fontis dans les terrains de recouvrement et ainsi hiérarchiser les paramètres associés à ce phénomène.La première partie de cette thèse repose sur une étude bibliographique qui récapitule les méthodes d’exploitation, méthodes d’analyse de stabilité de carrières souterraines, méthodes de prévision de la hauteur d’effondrement et estimation de foisonnement. A l’issue de cette synthèse bibliographique l’étude s’est focalisée sur les carrières souterraines à faible profondeur exploitées par chambres et piliers. Ainsi, la modélisation numérique par la méthode des éléments discrets (MED) a été choisie pour analyser l’instabilité des toits de carrières souterraines.La seconde partie s’intéresse au développement d’un modèle numérique qui a pour objectifs : d’une part, le développement d’un Programme de Discrétisation des Massifs Rocheux (PDMR) qui constitue le préprocesseur du logiciel STTAR3D et le développement d’un code permettant le calcul du coefficient de foisonnement des débris de l’effondrement. D’autre part, l’implémentation des lois de comportement sur STTAR3D.La troisième partie consiste à déterminer, d’une part les caractéristiques physico-mécaniques d’échantillons prélevés dans la carrière de la Brasserie (Paris-France), qui a été choisie pour une tester le modèle développé et d’autre part, les deux paramètres de la loi de comportement utilisée pour modéliser les contacts à savoir  et µ.Enfin, la dernière partie de ce travail est constituée des simulations numériques pour lesquelles les paramètres de la loi de comportement mesuré expérimentalement ont été introduits dans STTAR3D. Dans la première étude numérique menée, on s’intéresse à l’effet de la hauteur de chute, du rayon de l’ouverture initiale du fontis et du degré de fracturation sur le foisonnement des décombres, ainsi qu’à l’effet de la variation du foisonnement sur la hauteur de l’effondrement et sur l’affaissement. Dans un second temps, on réalise un modèle de la carrière de la Brasserie dont on calcule le comportement par simulation numérique afin d’obtenir l’affaissement en surface et la hauteur de l’effondrement qui sont comparés aux observations in-situ / Over time, the underground cavities are subjected to aging and several types of degradation can occur. The old underground cavities have probably not been designed to be stable over the long term. They have sometimes been totally or partially backfilled, but not in a systematic way. The collapse of a mine causes deconsolidation of the upper levels of the overburden. These mechanisms can cause two types of disorders on the surface: subsidence or sinkhole. Subsidence and sinkhole can cause severe damage to structures and infrastructures in surface, but also jeopardize the safety of the population.The work of this thesis revolves around the study of rock's bulking, sinkhole and its spread in the overburden. The aims of this thesis are twofold: firstly, to study the bulking of rock during the roofs mine collapse; Secondly, modeling the spread of the sinkhole in the overburden and thus to prioritize the parameters associated with this phenomenon.The first part of this thesis is a bibliographical study which summarizes the methods of exploitation, methods of analysis of stability of underground quarries, methods of prediction of the height of collapse and estimation of the bulking factor. At the end of this bibliographic synthesis, the study focused on shallow underground quarries operated by rooms and pillars. Thus, numerical modeling using the discrete element method (MED) was chosen to analyze the instability of roofs of underground quarries.The second part focuses on the development of a numerical model with the following objectives: on the one hand, the development of a Rock Mass Discretization Program (RMDP) which constitutes the preprocessor of the STTAR3D software and the development of a Code allowing calculation of the bulking factor of the rubble of collapse. On the other hand, implementation of the behavior laws on STTAR3D.The third part consists of determining, on the one hand, the physicals and mechanicals characteristics of samples taken from the quarry of the Brasserie (Paris-France), which was chosen to test the model developed. On the other hand, determining of parameters of the behavior law used for modeling the contacts, namely “” and “μ”.Finally, the last part of this work is made of numerical simulations for which the parameters of the behavior law measured experimentally have been introduced in STTAR3D. In the first numerical study, we investigate the effect of fall height, the radius of the initial opening of the sinkhole and the fracturing degree on the bulking of the rubble, as well as the effect of variation of the bulking on the collapse height and on the subsidence. In a second step, a model of the Brasserie’s mine is realized, the behavior of which is studied by numerical simulation in order to obtain the subsidence on the surface and the collapse height, which are compared with the in-situ observations

Les vagues courtes sont omniprésentes à la surface des océans, avec des longueurs de quelques dizaines de mètres au mètre typiquement. Connaitre leurs directions de propagation en mer est important à plusieurs titres, notamment pour la compréhension de la dynamique de l'état de mer, des échanges air-mer ou de la dérive de particules en surface. Ces distributions directionnelles sont étudiées ici au regard des progrès récents réalisés en techniques d'instrumentation. L'analyse du bruit sismo-acoustique enregistré en grandes profondeurs permet d'extraire un comportement quasi-universel qui dépend indirectement de cette distribution à travers ladite intégrale de recouvrement. Il est cohérent avec des observations directes du champ de vagues obtenues à partir de reconstructions tridimensionnelles de la surface de l'océan. Alors que la direction de propagation des vagues longues s'aligne avec celle du vent, les vagues courtes s'en détachent d'autant plus à mesure que leurs échelles diminuent (bimodalité).La comparaison de ces observations avec les prédictions d'un modèle numérique de vagues, basé sur l'environnement WAVEWATCH®III, permet de constater que ces modèles sont qualitativement valables mais encore quantitativement incorrects. Une des possibilités explorées pour corriger cet effet est la prise en compte de sources de vagues courtes à ±90° de la direction du vent, qui pourraient être associées au déferlement des vagues longues. Une telle source à elle seule n'explique pas les formes des distributions directionnelles observées. D'autres mécanismes pourraient intervenir que de futures investigations pourront tenter de clarifier. / Short surface gravity waves are ubiquitous at the ocean surface, with lengths from a few tens of meters to a meter typically.Knowing their propagation directions at sea is important in several respects, especially for the understanding of sea-state dynamics, airsea interactions and particles surface drift.Their directional distributions are here investigated in the light of the recent progress made in instrumentation techniques. The analysis of ocean bottom seismo-acoustic noise records allows for the extraction of a quasi-universal behavior which indirectly depends on this distribution through the socalled overlap integral. It is coherent with direct observations of the wave field obtained from tri-dimensional reconstructions of the ocean surface elevation field. While the propagation direction of long waves aligns with the wind direction, short waves progressively detach from it towards small scales (bimodality).Comparing those measurements with the predictions of a spectral numerical wave model, based on WAVEWATCH®III environment, allows to realize that they provide qualitatively correct but quantitatively incorrect predictions. One of the possibilities here explored to correct for it, is by accounting for the sources of energy at ±90° to the wind direction, which could be associated with the breaking of long waves. This source term on its own does not explain the shapes of the observed directional distributions. Other mechanisms could come into play that future investigations will help clarify.

Les pieds d’aubes de soufflantes de turboréacteurs font face à des sollicitations de type fretting. Il en résulte deux types d'endommagements: (i) l'amorçage et la propagation de fissures, (ii) l'usure des surfaces en contact. OBJECTIF: Afin de fournir les outils permettant de répondre à la problématique industrielle, une méthode semi-analytique permettant de traiter des problèmes de contacts élasto-plastiques revêtus et/ou hétérogènes est développé à partir d'éléments existants et de solutions analytiques et numériques novatrices. METHODE: La structure est simplifiée en supposant un contact entre deux massifs élastiques semi-infinis. Des solutions analytiques donnant pour: + les contributions élémentaires de chargements normaux et tangentiels constants sur une surface rectangulaire + les contributions élémentaires de déformations plastiques supposées constantes sur un volume parallélépipédique + les contributions élémentaires de déformations d’incompatibilité liées à un problème hétérogène (inclusions, revêtements, endommagement,...) supposées constantes sur un volume parallélépipédique. >> les déplacements en surface ou les contraintes dans le volume. Les déplacements en surface ou les contraintes dans le volume sont alors exprimés en utilisant des produits de convolution discrets entre des coefficients d'influence et la source surfacique (chargements en surface) ou volumique (déformations plastiques ou d'incompatibilité). Le problème normal et le problème tangentiel en glissement total ou en glissement partiel peuvent alors être résolus, en prenant en compte les effets plastiques et hétérogènes. L'algorithme d'optimisation sous contrainte utilisé (contact elastique) est celui développé par L.Gallego tandis que la base du solver plastique utilisé a été développé par C. Jacq (contact elasto-plastique sans frottement) RESULTATS: De nouvelles solutions analytiques sont obtenues pour le calcul des déplacements résiduels tangentiels. / Fan blade roots of jet engines are operating under fretting sollicitations. This results in two types of damage: (i) initiation and crack propagation, and (ii) wear of contact surfaces. OBJECTIVE: A semi-analytical method has been developed in order to solve the industrial problem. This method is fast-computing and able to deal with elastic-plastic coated contacts and / or heterogeneous, and has bee developed from existing items and innovative analytical and numerical solutions. METHOD: The structure is simplified by assuming two semi-infinite elastic solid in contact, and analytical solutions for: + Contributions of elementary normal and tangential loads constant over a rectangular area + Contributions of elementary plastic strains constant on a parallelepiped + Contributions basic incompatibility of deformations related to a heterogeneous (inclusions, coatings, damage ,...) constant on a parallelepiped. >> on surface displacements and stresses in the volume. Surface displacements and stresses in the volume are then expressed using discrete convolution products (2D or 3D-FFT) between the coefficients of influence and the surfacic (loading surface) or volumic (plastic deformation or inconsistency) source. The normal and tangential problems in both the gross-slip and stick-slip regime can be solved, taking into account the effect of plasticity and heterogeneities. The contact algorithm used (elastic contact) has been developed by L. Gallego while the base of the plastic solver has been developed by C. Jacq (frictionless elastic-plastic contact) RESULTS: New analytical solutions are given for calculating the residual tangential displacements.

La durée de vie en fonctionnement des pièces mécaniques peut être augmentée par la présence de contraintes résiduelles de compression. Inhérentes à la plupart des procédés de fabrication, contraintes résiduelles jouent un rôle sur la tenue en service des éléments mécaniques. La connaissance et la maîtrise de ces contraintes résiduelles et des procédés associés sont donc fondamentales. La mise en place de méthodes numériques performantes pour prédire ces contraintes résiduelles permettra à terme d'éviter de nombreux et coûteux essais et d'étudier l'influence des principaux paramètres. Cette thèse présente le développement et l'application de méthodes de calcul semi analytiques (SA) à la modélisation de procédés mécaniques de mise en compression des surfaces. Les méthodes SA, initialement développées pour la simulation des contacts élasto-plastiques, ont l'avantage de temps de calcul très significativement réduits par rapport aux méthodes numériques conventionnelles. Cette méthode est d'abord utilisée pour la simulation d'un procédé connu sous le nom de galetage, avec un chargement de type contact roulant. Enfin elle est utilisée pour la simulation d'impacts uniques puis répétés, et les développements associés présentés. Le contact roulant entre deux massifs élasto-plastiques, sans considération du frottement, est d'abord étudié. L'influence de la plasticité et du type d'écrouissage (isotrope ou cinématique), de la géométrie des massifs en contact (sphéres ou ellipsoides) et du type de chargement (indentation ou roulement) sur la pression de contact et les déformations plastique sont analysés. La simulation d'impacts est ensuite abordée. La méthode développée est tout d'abord validée numériquement, puis confrontée à l'expérience. Trois matériaux ont été plus particulièrement étudiés : 316L, AA 7010, et l'Inconel 600. Les dimensions des impacts ainsi que les déformations générées, mesurées par stéréo corrélation sont utilisées pour valider expérimentalement la méthode. Le procédé de billage ultrasonore a été tout particulièrement étudié. Dans un premier temps nous nous sommes attachés à la décrire la cinématique des billes dans une chambre fermée, mises et maintenues en mouvement par une sonotrode. L'utilisation de formules analytiques pour l'estimation des coefficients de restitution, lors des nombreux chocs entre les billes ou avec les parois de la chambre, a permis d'affiner le calcul des vitesses moyennes d'impact en fonction des paramètres du procédé. La méthode SA est ensuite utilisée pour déterminer le champ de déformations plastiques induits par les impacts. Enfin une méthode de projection est proposée pour permettre in fine la prédiction des contraintes résiduelles pour des structures fines ou épaisses. / The life time of mechanical components can be increased by the presence of compressive residual stresses. Inherent to most production processes, residual stresses play a critical role on the mechanical parts behaviour. The knowledge and mastering of residual stresses and linked processes are thus fundamental. The development of efficient numerical methods to predict these residual stresses will allow to save costly experiments and to study the influence of the main parameters. This PhD presents the development and application of semi analytical methods (SAM) to the modelling of mechanical processes of compressive residual stress generation. The SAMs, initially developed for the simulation of elasto-plastic contacts, have the advantage of significantly low computation time compared to classical numerical methods. This method is first used to simulate the low plasticity burnishing process, with a rolling loading. Then, it is used for the simulation of impacts, first unique and then repeated. The frictionless rolling contact between two elasto-plastic bodies is first studied. The influence of plasticity, of the hardening model (isotropic or kinematic), of the geometry of the bodies in contact (spheres or ellipsoids) and of the loading type (indentation or rolling) on the contact pressure and plastic strains are analysed. Impacts simulation is then addressed. The developed method is first validated numerically then confronted to experimentations. Three materials have been particularly studied: 316L, AA 7010 and Inconel 600. The impacts dimensions and the generated strains, measured by digital image correlation, are used to validate experimentally the method. The ultrasonic shotpeening process has been specifically studied. The description of the kinematics of the shots put in movement by a sonotrode in a closed peening chamber has first been studied. The use of analytical formulae for the estimation of the coefficients of restitution, during the numerous impacts between shots and with the chamber’s walls, allowed refining the calculation of the average impact velocity as a function of the process parameters. The SAM is the used to determine the plastic strain field induced by the impacts. At last a projection method is proposed to finally determine the residual stress field in thick or thin structures.

Ce travail porte sur la morphodynamique des bancs de galets alternés dans une rivière de montagne aménagée. Le site expérimental est un tronçon long de 8 km de la rivière Arc située dans les Alpes françaises. Une analyse des cartes historiques, des photos aériennes et des relevés topographiques montre que le lit de la rivière a évolué d’un style en tresse vers un lit rectiligne avec des bancs alternés contraints entre deux digues. Au milieu des années 1990, le lit de la rivière a été aplani pour former des profils trapézoïdaux. Cependant, les bancs alternés ont réapparu en moins de dixans. Une analyse des facteurs contrôlant leur formation est fournie. Les observations de terrain et l’analyse numérique 1D suggèrent que le motif de bancs alternés représente un état de quasi équilibre du tronçon. D’autre part, la topographie d’un sous-tronçon incluant un banc de galets aété mesurée 18 fois entre 2006 et 2012. Une méthodologie de construction des MNT basée sur leslignes directrices et incluant une estimation de l’erreur a été proposée pour une application aux litsde rivières. La technique de différenciation des MNT a été utilisée pour détecter les changements morphologiques du banc. Ainsi, la morphodynamique du sous-tronçon a été étudiée en la reliant aux principaux processus qui la contrôlent tels que l’hydrologie, les formes des hydrogrammes, la répartition des contraintes de cisaillement, les apports amont en sédiments et l’énergie du coursd’eau. La forme des hydrogrammes de crue ainsi que les apports en sédiments semblent être lesprincipaux facteurs influençant le bilan sédimentaire sur le site étudié. En outre, des modélisationsnumériques 2D hydro-sédimentaires ont été utilisées pour appréhender la dynamique du banc. Lessimulations ont révélé l’importance de l’interaction entre les sédiments fins et les galets sur ladynamique globale du banc. / The current research investigates the morphodynamics of gravel alternate bars in an engineeredmountainous river. The experimental site corresponds to an 8 km long reach of the Arc River, in the French Alps. An analysis of historical maps and recent aerial photographs and topographic campaigns shows that the original braided river bed has evolved to a straight bed with alternate bar system constrained within two embankments following two centuries of man-made engineeringprojects. In the mid 1990s the river bed was flattened to fit a trapezoidal cross-sectional design. However, alternate bars re-appeared in less than ten years. An analysis of the controlling factors of bar formation and evolution is provided. Field observation and 1D numerical analysis suggest that the alternate-bar pattern represents a quasi-equilibrium state of the river bed. A singlegravel bar within the reach was monitored 18 times from 2006 to 2012. A methodology of DTM construction based on breaklines is proposed for river beds and an error estimation is presented. DTM differencing technique was used to detect morphologic changes of the gravel bar. Therefore, the gravel bar and main channel dynamics were analyzed relative to main controlling processes suchas extreme hydrologic events, sediment supply, hydrograph shape, bed shear stress distribution andriver energy. Sediment supply and hydrograph shape appeared to be the main factors influencing the study site sediment budget. Moreover, 2D numerical modelling of water flow and sediment transport was performed. Sensitivity analysis was conducted to yield results in agreement with observations. Simulations revealed the relevance of the interaction between fine sediments and gravels on the bar dynamics.

Les procédés additifs, auxquels appartient la fusion laser de poudres métalliques, ont la capacité de créer des structures à géométries complexes, avec la possibilité d'intégrer des formes creuses, par exemple des canaux de refroidissement assurant un contrôle thermique optimum. Ce procédé permet de fabriquer des pièces réelles à partir de poudres métalliques, par fusion du matériau, couche par couche, en accord avec le modèle CAO. Au cours du procédé, de nombreux cycles thermiques et d'importants gradients thermiques se produisent dans la pièce au cours de sa fabrication. Ces gradients de température induisent des déformations plastiques hétérogènes et de ce fait des contraintes résiduelles. Ces contraintes peuvent nuire à la qualité de la pièce obtenue, par exemple sa résistance mécanique. Ces travaux ont pour objectifs de proposer un modèle numérique, s’appuyant sur la méthode des éléments finis afin d'étudier l'apparition des contraintes résiduelles lors du procédé de fusion laser de poudres métalliques. Le logiciel multiphysique ABAQUS® a été utilisé pour effectuer les analyses thermiques et mécaniques. La technique « d'ajout et de suppression des éléments » a été utilisée afin de simuler la fusion et la solidification de la matière au cours du procédé. Les propriétés mécaniques dépendantes de la température de l'acier maraging, utilisé dans notre cas, ont été obtenues à l’aide d’essais expérimentaux de caractérisations et intégrées dans le modèle. Les calculs sont réalisés de manière découplée, dans un premier temps le calcul thermique est effectué, puis les résultats sont utilisés pour réaliser le calcul mécanique et finalement prédire les champs de contraintes. Dans le cadre de ce travail, une méthode originale s'appuyant sur la technique de mesure des contraintes résiduelles par enlèvement de couches successives a été mise au point pour mesurer ces contraintes en direct au cours du procédé. Les résultats renseignent sur le niveau et la distribution des contraintes dans la pièce créée et le support. Deux paramètres ont été testés afin d'étudier leur influence sur le niveau des contraintes résiduelles : le temps d’étalement de la poudre entre deux couches successives et la hauteur des couches. Le modèle numérique paramétrable permet d'analyser les effets de paramètres liés au procédé sur la répartition des contraintes résiduelles dans les pièces fabriquées. Les résultats montrent que la variation de l'épaisseur du support n'affecte pas la répartition des contraintes dans la pièce créée. Le préchauffage du support à une température de 800°C réduit les contraintes résiduelles L'étude de quelques trajectoires laser montre leurs influences sur la répartition des déformations plastiques cumulées ainsi que la hauteur des couches de poudres ou de la forme du support (embase, colonnes). / The Selective Laser Melting process, belonging to Additive processes , have the ability to create structures with complex geometries , with the possibility of including cavities, such as cooling channels providing optimum temperature control. This process enables the manufacture of three-dimensional parts from metal powders by melting the material , layer by layer, in agreement with the CAD model. In the process , high temperatures and thermal gradients cycles occur in the part during the process. These temperature gradients induce heterogeneous plastic strain and residual stresses. These residual stresses may affect the quality of the part obtained, for example the fatigue life. This work aims to propose a numerical model , based on the finite element method to study the appearance of residual stresses during laser melting process of metallic powders . The ABAQUS® Multiphysics software was used to perform the thermal and mechanical analyzes. The movement of the laser beam and the resolution of the thermal problem can predict the evolution of the temperature in the part and support. The "birth and death elements" technique was used to simulate the melting and solidification of the material during the process. Dependent mechanical properties of the temperature of the maraging steel used in this case were obtained using experimental testing and characterization and were established in the model. The calculations are decoupled : initially thermal calculation is performed and the results are used to perform mechanical calculations and finally predict the residual stress fields. In this work, a novel method based on the technique of measuring residual stresses by removing layers was developed to measure these stresses directly in the process. The results provide information on the level and distribution of stresses in the created part and support. Two parameters were tested to study their influence on the level of residual stress : time to spread the powder between two successive layers and layer height. The model is used to analyze the effects of process parameters related to the distribution of residual stresses in the manufactured parts. The results show that the variation of the thickness of the support does not affect the distribution of stresses in the part created. Preheating the substrate to a temperature of 800 °C reduces the residual stresses. The study of some laser strategies shows their influence on the distribution of plastic strain thus the height of the layers of powder or in the form of support (base, columns).

Dans le contexte environnemental actuel, l’amélioration des performances énergétiques des transmissions mécaniques par engrenage est un réel challenge. De part, leur compacité et leur arrangement axisymétrique, les transmissions mécaniques de type réducteur à trains épicycloïdaux sont de plus en plus répandues dans divers applications (éolien, aéronautique,…). Il est généralement admis que les principales sources de dissipation de puissance dans de telles transmissions sont dues au frottement aux dentures, au mode de lubrification (barbotage ou injection de lubrifiant), au piégeage d’huile entre les dents lors de l’engrènement et à la ventilation des mobiles. Il est également à noter que les pertes de puissance générées par une transmission ne peuvent être découplées de son comportement thermique. En effet, les échauffements locaux dans la transmission ont un impact sur les propriétés du lubrifiant, qui ont elles-mêmes une influence sur les pertes de puissance. A partir de ce postulat, le travail présenté dans ce manuscrit propose un modèle numérique permettant de quantifier les différents postes de pertes de puissance générées dans un réducteur épicycloïdal, lubrifié par injection d’huile, pour une application aéronautique. / In the current environmental context, gearbox efficiency has become a major issue. Because of their compactness and axi-symmetric arrangement, planetary gearboxes are widely used in several applications (such as wind, aerospace…). It is generally accepted that total power losses in such gear transmissions can be decomposed into the contributions of the friction between teeth, the lubrication process (oil splash or jet), the oil trapping during meshing and tooth windage. It is noticed that power losses produced by a mechanical transmission cannot be dissociated from its thermal behavior. Indeed, the local warmings impact lubricant physical properties and these last have influence on power losses in return. Based on this postulate, the purpose of this work is to develop a numerical model in order to estimate the power losses generated into a jet lubricated planetary gear train, for an aeronautical application.

De nos jours, la modélisation numérique constitue un outil indispensable pour le prototypage de convertisseurs électromagnétiques. Les matériaux magnétiques jouent un rôle essentiel dans la conversion de l’énergie, il est donc nécessaire de maîtriser leur comportement et leur représentation. L’objectif de ce travail s’inscrit dans ce cadre et s’attache à élaborer des lois réalistes de comportement de matériaux afin de les inclure dans des simulateurs de circuits. Concernant le comportement statique, le modèle de Jiles-Atherton a été implémenté puis adapté, simplifié et modifié afin d’en améliorer la précision et l’implémentation. La modélisation dynamique du matériau a été effectuée grâce au modèle DWM élaboré au laboratoire Ampère. Ce modèle intègre les effets dynamiques excédentaires grâce à une loi « dynamique de matériau » implémentée au sein de l’équation de diffusion magnétique. Ce modèle a été ensuite homogénéisé afin d’en améliorer son implémentation future dans un simulateur de circuit. Chacun des différents modèles a été testé et validé sur plusieurs échantillons. / Nowadays, numerical modeling is an indispensable tool for the prototyping of electromagnetic converters. Magnetic materials play an essential role into the energy conversion so it is necessary to control their behavior as well as their modeling. The objective of this work is to develop realistic laws of material behavior for circuit simulators use. Regarding the static behavior, the Jiles-Atherton model has been implemented and adapted, simplified and modified to improve accuracy and implementation. Dynamic modeling of the material was performed using the model DWM developed into the Ampere laboratory. This model incorporates the excedentary dynamic effects thanks to a "dynamical material law" implemented into the magnetic diffusion equation. Then this model was homogenized to improve its future implementation in a circuit simulator. Each of the different models has been tested and validated on several samples.

Le fretting et l'usure sont des problèmes récurrents dans le domaine de l'aéronautique. Les contacts aube/disque au niveau du compresseur ou de la turbine haute pression des moteurs d'avion, par exemple, sont soumis à d'importantes sollicitations de fretting à de fortes températures. L'enjeu des industriels est d'optimiser la durée de vie de ces composants et d'être capable de prévoir l'amorçage de fissures. Afin d'améliorer la tenue des pièces, des revêtements sont utilisés pour les protéger. Leurs propriétés mécaniques et matériaux ont un impact direct sur le contact et la durée de vie. Les choix de matériaux, du nombre de couches, de l'épaisseur, de l'ordre sont donc primordiaux. De par leur composition (fibres, mono-cristaux), leur élaboration (extrusion) ou leur mode de déposition, l'hypothèse de considérer des matériaux homogènes isotropes s'avère trop réductrice. L'anisotropie est un paramètre important à prendre en compte au niveau du dimensionnement. Les matériaux composites sont de plus en plus utilisés dans l'aé-ronautique. Dans cette optique, cette thèse a pour objectif l'étude du comportement des matériaux homogènes anisotropes, en s'intéressant à l'influence des principaux paramètres mécaniques caractéristiques d'un matériau afin de mieux appréhender leurs effets. On s'attardera sur le module de Young (ou module d'élasticité), le module de Coulomb (ou module de cisaillement) et le coefficient de Poisson, et leurs valeurs selon les différentes directions. Comme attendu, le module de Young dans la direction normale au contact joue un rôle prépondérant dans la détermination du profil de pression. Néanmoins, l'influence du module de Young dans le plan tangent au contact n'est pas à négliger, il modifie aussi la forme de l'aire de contact. L'orientation du matériau par rapport au contact est par conséquent un paramètre à prendre en considération, il peut directement atténuer ou accentuer l'effet du module de Young dans une direction priviligiée. Les module de Coulomb et coefficient de Poisson ont aussi été analysés. Il en résulte qu'ils influent significativement sur le contact. Ces résultats se confirment dans le cas d'un massif revêtu, à la différence que les effets du revêtement et du substrat peuvent se compenser. L'impact des propriétés du revêtement sera d'autant plus important que celui-ci sera épais. L'échelle du contact par rapport aux matériaux utilisés importe aussi sur les profils de pression. Une comparaison entre le modèle anisotrope homogène et un modèle isotrope hétérogène a été réalisée. A l'échelle mesoscopique, le composite est composé d'une matrice avec des fibres qui induisent des pics de pression alors qu'à l'échelle macroscopique, le matériau composite est perçu comme un matériau homogène, les profils de pression sont lissés. / Fretting and wear are recurrent problems in the field of aeronautics. Contacts the blade / disk at the compressor or high-pressure turbine aircraft engines, for example, are subjected to high stresses at high temperatures. The challenge for manufacturers is to maximize the lifetime of these components and be able to predict crack initiation. To improve handling parts, coatings are used to protect them. Materials and their mechanical properties have a direct impact on the contact and the lifetime. The choice of materials, number of layers, the thickness of the order are therefore essential. By their composition (fibers, single crystals), elaboration (extrusion) or their mode of deposition, the hypothesis to consider homogeneous isotropic materials is too simplistic. Anisotropy is an important parameter to take into account in the design. Composite materials are increasingly used in the aeronautic. In this context, this thesis aims to study the behavior of anisotropic homogeneous materials, focusing on the influence of the main parameters mechanical characteristics of a material to better understand their effects. We focus on the Young's modulus (or modulus of elasticity), the module Coulomb (or shear modulus) and the Poisson's ratio, and values ​​in different directions. As expected, the Young's modulus in the direction normal to the contact plays an important role in determining the pressure profile. However, the influence of Young modulus in the plane tangent to the contact is not to neglect it also alters the shape of the contact area. The orientation of the material with respect to the contact is therefore a parameter to take into consideration, it can directly reduce or enhance the effect of Young's modulus in a direction privileged. The module Coulomb and Poisson's ratio were also analyzed. As a result they significantly affect the contact. These results are confirmed in the case of a coated solid, unlike the effects of coating and substrate can compensate. The impact properties of the coating will be even more important than it is thick. The scale of the contact relative to the materials used is also important on the pressure profiles. A comparison between the model anisotropic homogeneous and isotropic heterogeneous model have been performed. At mesoscopic scale, the composite is composed of a matrix with fibers that induce pressure peaks while at the macroscopic level, the composite material is seen as a homogeneous material, the pressure profiles are smoothed.

L'aquifère des basaltes du Dalha de la région de Dikhil, au Sud de la République de Djibouti, a fait l'objet d'une étude basée sur une approche pluridisciplinaire (hydrochimique, isotopique, hydrodynamique, modélisation 2D des écoulements) afin de comprendre le fonctionnement de ce système aquifère dans le cadre d'une gestion de ressources en eaux.L'étude hydrochimique a permis de mettre en évidence, d'une part, les processus d'acquisition de la minéralisation des eaux suivant les principales directions d'écoulement, et d'autre part, des faciès hydrochimiques complexes sur la base de la dominance respective des cations et anions. Trois faciès hydrochimiques sont prédominants, à savoir : Na-SO4-HCO3 représentant les eaux à minéralisation modérée, Na-HCO3-SO4 caractéristique des eaux à faible minéralisation et Na-Mg-Cl-SO4 qui correspond aux eaux à forte minéralisation et circulant principalement dans les aquifères des basaltes du Dalha et stratoïdes. Les résultats de l'analyse des isotopes stables de l'environnement (d18O et d2H) sont en parfait accord avec ceux issus de l'hydrochimie, et indiquent que les eaux souterraines de la zone d'étude sont d'origine météorique.D'un point de vue hydrogéologique, les résultats des pompages d'essais exécutés sur les forages captant l'aquifère du Dalha et l'aquifère stratoïde, ont révélé le fort degré d'hétérogénéité des ces aquifères.Le modèle numérique d'écoulement bidimensionnel de l'aquifère des basaltes de Dalha, en régime permanent, a confirmé l'hétérogénéité de l'aquifère mis en évidence par les pompages d'essais et a permis d'établir un bilan de flux qui indique que l'exploitation de l'aquifère a atteint une limite critique / The Dalha basalts aquifer of the region of Dikhil, in the South West of the Republic of Djibouti, was the object of a study based on a multidisciplinary approach (hydrochemical, isotopic, hydrodynamics, groundwater modeling) to understand the functioning of this aquifer system within the framework of a sustainable management of its water resources.The hydrochemical study pointed out, on one hand, the processes of mineralization along the main water flow directions, and on the other hand, complex hydrochemical facies on the basis of the respective dominance of cations and anions. Three hydrochemical facies are dominant, namely: Na-SO4-HCO3 representing waters with moderate mineralization, Na-HCO3-SO4 facies characteristic of waters with low mineralization and Na-Mg-Cl-SO4 which corresponds to waters with strong mineralization and circulating mainly in Dalha and stratiform basalts aquifers. Results of environmental stable isotopes (18O and 2H) are in perfect agreement with those of hydrochemistry, and indicate that groundwater in the study area are of meteoric origin. In hydrogeological point of view, the results of pumping tests performed on wells capturing these aquifers, revealed the high degree of heterogeneity of these aquifers.The numerical groundwater modeling of Dalha basalts aquifer performed in steady state has confirmed the heterogeneity of the aquifer highlighted by pumping tests and helped to establish a water balance which indicates that the exploitation of the aquifer has reached a critical limit beyond which overexploitation may induce an irremediable degradation of the resource.