Dans ce travail, nous proposons une approche essentiellement expérimentale. Dans une première partie, nous exposons les résultats obtenus par imagerie rapide et mesures électriques corrélées dans le temps, concernant l'étude de la transition entre des électrodes de cuivre dans l'air et l'argon, et entre des électrodes de tungstène dans l'argon pour des pressions comprises entre 100 mbar et 1 bar. Dans aucun cas une transition à l'arc n'est observée sans constriction quasi-totale de la colonne positive. Dans le cas d'électrodes en tungstène dans l'argon, la vitesse d'extension de la zone filamentaire depuis l'anode vers la cathode a été vue dépendante de la pression et de l'intensité. La durée d'une transition spontanée a été linéairement corrélée à la distance inter-électrode, suggérant un mécanisme propagatif à vitesse constante dans la colonne. Plus la pression est élevée, plus la durée totale des transitions est longue ; par ailleurs, à pression constante, cette durée de transition spontanée décroît à intensité croissante, pour tendre vers une durée limite, correspondant à la durée constante de changement de structure de la zone cathodique. Dans une seconde partie, nous nous intéressons à l'échauffement de surfaces d'anodes et de cathodes en cuivre dans l'air et dans l'argon soumises à des arcs électrique pour plusieurs gammes de courant (de 30 à 120 A). Deux méthodes de mesures sont utilisées. Enfin, nous utilisons les mesures de températures de surface pour proposer deux modèles simples de bilan de puissance à la surface des électrodes (flux de chaleur). L'un considère uniquement le chauffage dû au spot d'arc et l'autre le chauffage supplémentaire de la colonne d'arc.

Les systèmes de positionnement et de navigation par satellite (GNSS) ont connu un formidable essor au cours des deux dernières décennies au point d'être aujourd'hui un élément technologique universel. Le recours au positionnement précis par satellite est intégré à notre vie quotidienne: puces GPS intégrées dans les téléphones portables, systèmes de radionavigation utilisés pour le guidage d'engins agricoles, l'approche aérienne et maritime. Les stations GNSS déployées dans les réseaux géodésiques permanents participent au maintien des références géodésiques nationales et internationales et à la prévision météorologique. Ces applications scientifiques requièrent une mise en œuvre rigoureuse des équipements et un traitement sophistiqué des mesures. Une limitation de précision commune vient la présence d'erreurs dans les mesures dues à la traversée de l'atmosphère. L'amélioration des méthodes de correction de la propagation des signaux radiofréquence dans l'atmosphère a permis au cours des dix dernières années d'atteindre une précision sub-centimétrique en positionnement vertical et de faire des GNSS une technique de référence pour le sondage de la vapeur d'eau dans la troposphère. Ce manuscrit rassemble les résultats des recherches que j'ai développées au LATMOS et à l'IGN depuis une douzaine d'années dans le domaine du positionnement précis par GPS, de la télédétection de la vapeur d'eau atmosphérique par GPS et lidar Raman, et de l'étude du cycle de l'eau atmosphérique au moyen de systèmes d'observations et de modèles de prévision météorologique.

Au cours du 20ème} siècle, les enregistrements marégraphiques suggèrent une hausse du niveau de la mer de 1.8 mm/an. Plus récemment, les observations spatiales indiquent une hausse de 3.3 mm/an sur la période 1993-2009. Cette augmentation au cours du temps est attribuée au réchauffement global de la planète enregistré depuis plusieurs années maintenant. Durant cette thèse, nous analysons les observations et les causes de la hausse moyenne globale du niveau de la mer. Nous estimons les variations stériques du niveau de la mer grâce aux données du projet international Argo et, les variations de masse des océans liés aux apports d'eau des continents à l'aide des mesures de la mission GRACE. Une autre étude se concentre sur la variabilité du stock d'eaux continentales des plus grands bassins hydrologiques de la planète, à l'aide des données GRACE, et l'impact de cette composante à la hausse du niveau de la mer. Puis, nous analysons l'impact de la variabilité interannuelle du stock d'eaux continentales aux variations du niveau de la mer sur diverses périodes. Enfin, nous étudions le bilan des contributions climatiques à la hausse observée du niveau de la mer sur la période altimétrique totale et pour les années récentes. Dans une deuxième partie, nous étudions les variations régionales du niveau de la mer. Nous établissons dans un premier temps, les causes de la variabilité régionale du niveau de la mer et nous interprétons le signal résiduel issu de la différence entre le niveau de la mer observé et l'expansion thermique des océans. Nous étudions ensuite les variations passées des structures régionales du niveau de la mer sur les dernières décennies (1950-2003). Le but est de reconstruire les variations passées du niveau de la mer en 2-D et ainsi d'avoir une connaissance plus approfondie de l'évolution régionale du niveau moyen des mers. Ce travail a pour objectif ultime de contraindre les modèles climatiques couplés utilisés par l'IPCC (Intergovernmental Panel on Climate Change) pour prédire l'évolution future du niveau de la mer au cours du 21ème siècle. Cette analyse nous permet de détecter un signal basse fréquence dans la variabilité régionale du niveau de la mer que nous observons non seulement dans les données in situ mais aussi dans les modèles climatiques couplés.

Des quantités importantes de gaz réactifs et de particules sont émises en Afrique équatoriale par l'activité humaine et naturellement par les écosystèmes forestiers, les éclairs et les sols. Les différents processus de transport qui caractérisent cette région pendant la saison de la mousson (e.g. Jets d'Est Africain et Tropical, convection profonde) peuvent conduire à la redistribution de ces émissions en dehors du continent Africain et influencer ainsi aussi bien le bilan régional que global de l'ozone (O3). Le but de mes travaux de thèse est d'améliorer notre compréhension sur l'impact de ces émissions et de ces processus de transport sur la composition chimique de l'atmosphère en Afrique équatoriale. Les données mesurées pendant la campagne AMMA qui a eu lieu en été 2006 en Afrique de l'ouest ainsi que les données du programme MOZAIC sont analysées et confrontées aux résultats du modèle global LMDz_INCA. Cela permet d'évaluer les performances de ce modèle au dessus de l'Afrique équatoriale d'une part et, d'autre part, d'analyser l'influence des différentes émissions et des différents processus de transport qui caractérisent cette région pendant la saison de la mousson sur le bilan régional de l'ozone. Une description et une validation des résultats du modèle LMDz_INCA avec des observations satellitaires sont présentées dans la première partie de cette thèse. La deuxième partie est consacrée à l'impact du transport convectif et des oxydes d'azote (NOx) émis par les éclairs sur l'ozone et ses précurseurs. L'analyse de l'impact des émissions par les feux de biomasse depuis l'Afrique centrale est ensuite présentée. La dernière partie traite de la contribution des émissions africaines et asiatiques au bilan de l'ozone en Afrique équatoriale.

La réduction des impacts environnementaux des entreprises s'appuie notamment sur de nombreux instruments d'évaluation environnementale (ACV, Bilan Carbone, etc.) censés guider leurs actions. Or, malgré leur importante diffusion, la question de leur performativité, c'est-à-dire de leur capacité à transformer effectivement les pratiques managériales, reste posée. Autrement dit, nombre de ces outils d'évaluation ne deviennent pas ipso facto des instruments de gestion. Dans les travaux en ingénierie de l'environnement, les difficultés d'appropriation des instruments d'évaluation sont associées aux propriétés intrinsèques de ces outils (ergonomie, robustesse des modèles scientifiques) et font l'objet d'un travail de révision permanente. Dans cette thèse, nous explorons les dimensions organisationnelles et cognitives de cette appropriation, en particulier l'émergence de communautés épistémiques et pratiques structurées autour des instruments et qui constituent le milieu réceptif à leur déploiement.Issue d'un travail en CIFRE dans un cabinet de conseil en environnement, cette thèse étudie, à partir de six cas empiriques, les processus et les conditions dans lesquels les instruments enclenchent des dynamiques d'action collective. Au-delà du niveau élémentaire des instruments, nous étudions les dispositifs qui sont mis en œuvre dans le cours d'actions collectives ainsi que leurs effets inattendus. La thèse met en évidence l'existence de dispositifs plus ou moins élaborés qui influent sur les capacités d'apprentissage des entreprises. La caractérisation de ces dispositifs permet d'identifier de nouvelles figures d'acteurs et des communautés intermédiaires, plus ou moins intégrées à l'entreprise, jouant un rôle clé dans les dynamiques d'action collective.

L'étude expérimentale effectuée sur la précipitation du SrMoO₄ a montré qu'il présente plusieurs morphologies (bipyramide et fuseaux) possibles sous forme d'agglomérats multiéchelles. Les paramètres opératoires: concentration initiale en SrMoO₄, température et vitesse d'agitation ont une influence sur la morphologie et la taille des cristaux et agglomérats. La précipitation du SrMoO₄ a été suivie en phase liquide par conductimétrie dans un réacteur discontinu. En fin de précipitation, la forme des cristaux et des agglomérats a été observée au MEB/FEG et la distribution de taille de particules (DTP) mesurée.La poudre a été analysée par DRX, BET et porosimétrie Hg. L'objectif de ce travail est le développement d'un modèle d'agglomération couplé avec la nucléation et la croissance cristalline de façon à déterminer la morphologie des agglomérats obtenus. Le logiciel Parsival a été utilisé pour résoudre les bilans de population et ainsi obtenir les DTP et les profils de désursaturation.

Le fuligule milouin (Aythya ferina) est une espèce de canard plongeur répandue dans l’ensemble du Paléarctique. Montrant des signes de déclin de ses effectifs nicheurs en Europe depuis les années 2000, son statut de conservation IUCN a été réévalué de « Préoccupation Mineure » à « Vulnérable ». Etant une espèce chassée, la question du maintien de ces prélèvements cynégétiques est légitime. L’objectif de ce travail de thèse est donc de comprendre son fonctionnement démographique et les mécanismes sous-jacents à sa dynamique de population, pour proposer des mesures de gestion et de conservation. Pour cela, nous nous sommes attachés à estimer la tendance de ce déclin sur sa voie de migration nord-ouest européenne à partir des données de recensements de la mi-janvier. Puis, nous avons étudié la connectivité migratoire de cette voie de migration avec les autres voies préssuposées, afin de mieux comprendre l’origine de ce déclin. Enfin, nous avons estimé deux paramètres démographiques clés, la survie et le succès de nidification, pour les intégrer dans un cadre matriciel permettant d’estimer un taux de croissance asymptotique et comprendre quels sont les paramètres démographiques clés sur lesquels doivent reposer les actions de gestion. Il ressort de ce travail de thèse que le déclin du fuligule milouin est probalement dû surtout à une diminution de sa productivité en Europe de l’Est et en Russie. Le contrôle et la gestion de la productivité par l’Homme étant limités, uniquement des travaux de gestion autour des habitats liés à la nidification sont envisageables. Ces travaux pourraient être facilement menés en Europe mais beaucoup plus difficilement en Sibérie, principale aire de reproduction des milouins caractérisée par une vaste superficie et un fort isolement. Les taux de survie des oiseaux bagués en France sont plus faibles que pour les oiseaux bagués dans d’autres pays, potentiellement du fait d’une pression de chasse plus élevée. Amener le taux de survie en France au même niveau que dans les pays voisins (Suisse et Grande-Bretagne), par une modification de la réglementation sur la chasse, serait de nature à ramener le taux de croissance à l’équilibre. La mise en place de la gestion adaptative pourrait faciliter les objectifs de conservation fixés en adaptant annuellement ces prélèvements aux effectifs présents et aux connaissances dont on dispose. / The Common Pochard (Aythya ferina) is a regular diving duck species in the western Palearctic. However, a worrisome decline of its wintering population led to an up-listed IUCN status from “Least Concern” to “Vulnerable”. This species is still hunted in Europe despite this decline. Hence, one may wonder about the sustainability of its harvest. The aim of this work was to understand the population dynamics of this species, and the drivers of these mechanisms. For this purpose, we assessed the declining trend in northwestern Europe using the mid-january censuses. Then, we studied the migratory connectivity with the two others flyways, in order to better understand the origin of the decline. Finally, we assessed two main demographic parameters (survival rate and nesting success) and combined these into a matrix population model. This model allowed us to assess an asymptotic growth rate and to determine the key demographic parameters on which management actions should focus. The main results of this thesis indicate that the decrease in productivty in Europe and in Russia could have been the main reason of the decline. However, given the limited human action to improve Pochard productivity, only breeding habitat improvement could be considered. Such improvements could be easily considered in Europe, but not in Siberia, the main breeding region characterized by a large area and a strong geographic isolation. Survival rates were lower in France than in neighbouring countries, possibly owing to a greater hunting pressure. A more moderate hunting pressure could lead to higher survival rates, and a balanced growth rate. Setting up an adaptative harvest management scheme could help reaching defined management goals, by annually adapting hunting quotas to current knowledge and assessment of Pochard population size.

Les arcs électriques de défaut se produisant entre deux câbles endommagés peuvent donner lieu à un phénomène appelé " arc tracking ", qui se caractérise par le maintien et la propagation de cet arc le long des câbles, entraînant leur ablation progressive. Dans un réseau aéronautique, un tel défaut peut avoir de graves conséquences lorsqu'il apparaît au sein d'un toron de câbles, d'autant plus qu'il est susceptible de se transférer à une structure avoisinante. Un dispositif expérimental a été développé pour l'étude de ce phénomène. Il est constitué d'une alimentation conçue selon des spécifications propres au domaine aéronautique, permettant de générer un signal continu ou alternatif pendant une durée réglable. Afin de se rapprocher des conditions de vol (pression réduite à haute altitude), les essais sont effectués au sein d'une enceinte fermée, soit à 1 bar, soit à 95 mbar. Trois types de câbles sont testés, un en aluminium et deux en cuivre. En régime continu, ils sont reliés au + et au - de l'alimentation, et à deux phases en régime alternatif. Une plaque en aluminium, représentant la structure de l'avion et connectée au neutre du générateur en régime alternatif, est placée en regard de la partie dénudée des câbles, là où se produit l'arc, permettant un éventuel transfert de l'arc depuis les câbles vers cette plaque. Chaque test s'accompagne de l'acquisition systématique des signaux de courant et de tension, et deux capteurs de flux radiatif sont positionnés autour de l'arc. Une caméra rapide permet la visualisation du comportement de l'arc au cours de l'essai. Enfin, les échantillons de câbles sont pesés avant et après chaque test. Un bilan de puissance est réalisé pour les différentes configurations testées (courant, pression, câbles), afin d'estimer la répartition de la puissance totale entre les électrodes (fusion et vaporisation des câbles, conduction, rayonnement) et la colonne de plasma (rayonnement, conduction, convection). La part transférée aux électrodes est estimée au moyen de la valeur de la chute de tension aux électrodes et du courant. Une partie de cette puissance aux électrodes est à l'origine de la fusion et de la vaporisation des câbles (le matériau isolant comme le métal). Celle-ci est estimée grâce à un calcul thermodynamique à partir de la masse de câble ablatée, en considérant que la totalité a été fondue et que 1% a été vaporisé. La puissance transmise à la colonne d'arc est en partie rayonnée, mais le spectre de longueurs d'onde inférieures à 200 nm est absorbé dès les premiers microns d'air, et une estimation est effectuée à partir du calcul coefficient d'émission nette. Une partie importante de l'étude a porté sur le transfert de l'arc à la plaque en aluminium, en fonction de la distance entre les câbles et cette plaque et selon un critère de température en face arrière, mesurée par thermographie infrarouge. Il a ainsi été montré qu'à pression réduite l'arc est beaucoup plus diffus, et a tendance à s'accrocher à la plaque en aluminium à des distances environ deux fois plus grandes qu'à pression atmosphérique. Cependant, la puissance totale disponible étant inférieure d'environ 20%, les dégâts à la plaque en aluminium sont moins importants. / Fault electric arcs occurring between two damaged wires may lead to a phenomenon called "arc tracking", which consists in a maintaining and propagating arc along the cables, causing their progressive ablation. In an aeronautical network, such a fault may have serious consequences if it happens within a wiring harness, especially as it is likely to transfer to a nearby structure. An experimental setup has been developed to study this phenomenon. It consists of a power supply designed in accordance with aeronautical specifications, which enables to generate a direct or alternating signal during an adjustable time. In order to reproduce flight conditions in a better way (low pressure at high altitude), tests are carried out within a closed chamber, whether at 1 bar pressure, or at 95 mbar pressure. Three types of cables are tested, one mainly made of aluminium and two mainly made of copper. In direct regime, they are connected to the + and the - poles of the generator, and to two phases in alternating regime. An aluminium plate, representing the aircraft structure and connected to the neutral in alternating regime, is placed in front of the stripped part of the cables, where the arc occurs, enabling a potential transfer of the arc from the cables toward the plate. Every test involves a systematic acquisition of the current and voltage signals, and two radiative flux sensors are positioned around the arc. A fast camera enables a viewing of the arc behaviour during the test. Finally, the sample is weighed before and after each test. A power balance is realised for the various tested configurations (current, pressure, cables), in order to estimate the distribution of the total power between the electrodes (fusion and vaporisation of the cables, but also conduction and radiation) and the plasma column (radiation, conduction, convection). The amount of power that is transferred to the electrodes is estimated with the value of the electrode voltage drop and the current. Part of this power makes the cables melt and vaporise (the insulating material as well as the metal), it is estimated with a thermodynamic calculation based on the ablated mass of cable, considering that the whole has been molten, and 1% vaporised. The amount of power transmitted to the arc column is partly radiated, but the wavelength spectrum emitted below 200 nm is absorbed through the very first microns in air, and estimation is made involving the calculation of the net emission coefficient. A significant part of the study was related to the transfer of the arc toward the aluminium plate, as a function of the distance between the cables and that plate, and according to a temperature criterion at the back side, by means of infrared thermography. Thus, it was shown that at low pressure, the arc is much less localised, and tends to transfer to the aluminium plate at distances around twice bigger than at atmospheric pressure. However, the total available power is around 20% smaller, and the damage caused to the aluminium plate is less severe.

Cette thèse porte sur l’impact du type de sol sur le fonctionnement biogéochimique d’un écosystème forestier typique d’Europe centrale et occidentale. Pour cette étude, un dispositif expérimental a été mis en place par l’INRA et l’ANDRA au sein d’une hêtraie (Fagus sylvatica) dans la forêt de Montiers-sur-Saulx (Meuse, France) afin de suivre les cycles biogéochimiques entre les différents compartiments de l’écosystème (l’atmosphère, la canopée, la litière et le sol). L’intérêt particulier du site expérimental de Montiers vient du contraste important entre les types de sols pour un même peuplement : rendisol, calci-brunisol et alocrisol/brunisol. Pour la première fois, les stocks et les flux de l’eau et des éléments ont été mesurés et comparés in situ et sur le long terme (entre janvier 2012 et décembre 2015) sur des sols forestiers différents, toutes les autres conditions du site égales par ailleurs (climat, apports atmosphériques, âge et structure du peuplement). Les stocks totaux et échangeables dans le sol et les stocks dans les différents compartiments de la biomasse végétale (branches, tronc, racines grosses et fines) ont été déterminés par le biais d’analyses chimiques et de modélisations. Les flux élémentaires ont été calculés à partir des prélèvements mensuels des solutions de l’écosystème (dépôt atmosphérique, pluviolessivat, écoulement de tronc, solutions libres et liées du sol) et des chutes de litière. Le modèle hydrique BILJOU© a été utilisé pour estimer les composantes du bilan hydrique de l’écosystème. Les résultats montrent que la réserve utile du sol, les stocks des éléments dans le sol, les compositions chimiques des solutions du sol, la stratégie de colonisation racinaire et la production de biomasse pérenne diffèrent fortement en fonction du type de sol. Contrairement à ce que l’on pouvait supposer compte tenu des différences marquées entre les signatures chimiques des solutions des trois sols, les compositions foliaires du peuplement et les flux des éléments hors sol en solution (échange de la canopée, apport au sol) et sous forme solide (chute de litière) ne diffèrent pas ou peu entre les trois stations expérimentales pour la majorité des éléments. La part biologique des cycles minéraux est globalement prédominante et le renouvellement des racines fines représente un flux de recyclage généralement supérieur à la chute de litière aérienne. Nous montrons également que la quantité d’eau transpirée par la canopée, directement liée à la taille de la réserve utile du sol, est le facteur de contrôle principal de la productivité du peuplement à l’échelle annuelle. Des processus d’adaptation du peuplement aux conditions physico-chimiques du sol semblent réduire, voire compenser entièrement, les facteurs secondaires de contrôle de la productivité du peuplement (notamment la disponibilité des nutriments dans le sol). Ainsi un enjeu significatif pour les gestionnaires forestiers pourrait être l’adaptation des pratiques sylvicoles à des parcelles de gestion davantage basées sur les propriétés physiques des sols, et en particulier les profondeurs d’apparition de la roche-mère et de colonisation racinaire / This thesis deals with the impact of soil type on the biogeochemical functioning of a typical forest ecosystem of Central and Western Europe. For this study, a strongly instrumented experimental site was implemented in a beech stand (Fagus sylvatica) within the state forest of Montiers-sur-Saulx (Meuse, France), in order to monitor the biogeochemical cycling between the different compartments of the ecosystem (the atmosphere, the canopy, the forest floor and the soil). The particular value of the Montiers experimental site resides in the strong contrast between soil types under the same beech stand: Rendzic Leptosol, Eutric Cambisol and Dystric Cambisol. For the first time, stocks and fluxes of water and elements were measured and compared in situ and over the long term (from January 2012 to December 2015) on different forest soil types, all other site conditions being equal (climate, atmospheric inputs, stand age and structure). Total and exchangeable pools in the soil and stocks in the different compartments of the vegetation (branches, trunk, fine and coarse roots) were determined via chemical analysis and modelisation. Fluxes of elements were calculated from monthly sampling of the ecosystem’s solutions (atmospheric deposition, throughfall, stemflow, gravitational and bound soil solutions) and of litter fall. The water balance model BILJOU© was used to estimate the different components of the water budget. The results show that the soil water holding capacity, the stocks of elements in the soil, the chemical composition of soil solutions, the rooting strategy and the perennial biomass production differ strongly between soil types. Contrary to what might have been expected in regard to the marked differences between the chemical signatures of the soils solutions in the three soils, the foliar elemental composition of the beech stand and the aboveground fluxes of elements in solution (canopy exchange, stand deposition) and in solid state (litter fall) do not differ significantly between the three experimental stations for the majority of the studied elements. The biological part of the mineral cycles is overall predominant and the recycling through fine roots turnover is generally higher than litter fall. We also show that the quantity of water transpired by the forest canopy, directly linked to the soil water holding capacity, is the primary control factor of the annual stand productivity. Stand adaptation mechanisms to physico-chemical soil properties seemed to reduce, or even entirely compensate for, secondary factors controlling the stand productivity (in particular nutrient availability in the soil). Thus a significant issue for forest managers might be to further adapt forestry practices to management units based on soil physical properties, especially depth to the bedrock and rooting depth

L’étude du comportement des roches nécessite de comprendre leur réponse sous diverses sollicitations. L’étude énergétique de l’endommagement des roches est indispensable pour prédire les phénomènes dynamiques. Ces phénomènes sont dus au développement de fissures dans les roches soumises à des fortes contraintes initiales et induites. La fissuration est une forme de dissipation d’énergie qui permet de rétablir l'équilibre du milieu. L’objectif de la thèse est de modéliser la fissuration dans les milieux rocheux dans la perspective d’étudier le comportement des ouvrages souterrains à grande profondeur. Le développement des modèles capables de représenter la fissuration, la coalescence des fissures et leur interaction avec des fractures préexistantes est indispensable. De la littérature, il ressort deux principales approches théoriques et numériques de modélisation de la fissuration : continue et discrète. Une synthèse critique de ces approches nous a conduit à retenir l’approche discrète et plus particulièrement le code Yade dans le cadre de cette thèse. Ce code permet de simuler explicitement la fissuration avec ou sans fractures pré-existantes. Des développements ont été effectués afin de tenir compte des différentes formes d’énergie intervenant dans le comportement des roches soumises à des sollicitations. En particulier une corrélation entre l’énergie de fissuration numérique et l’activité microsismique observée sur des échantillons en laboratoire en compression. Les différentes composantes énergétiques explicitées puis implémentées dans Yade sont : le travail externe, l’énergie potentielle, l'énergie élastique, l’énergie de frottement, l’énergie de fissuration, l'énergie cinétique et l’énergie dissipée en amortissement. La validation de l’approche énergétique a été réalisée grâce à la simulation des essais de laboratoire. L’évolution des différentes composantes énergétiques a permis de vérifier que le bilan des énergies est correctement évalué. Le bilan énergétique a également été vérifié à échelle de structures en simulant l’excavation souterraine d’un Mine-by Experiment (URL Manitoba). L’extension de la zone endommagée induite par l’excavation prédite numériquement a été comparée à celle observée in-situ autour du Mine-by Experiment. Il a été constaté que l’endommagement prédit est similaire à celui observé dans les directions des contraintes initiales mineure et majeure. Par ailleurs la formulation énergétique permet d’étudier numériquement les processus de fissuration des roches. Wassermann (2006) a réalisé des essais de compression uniaxiale et triaxiale sur des échantillons de minerai de fer lorrain. Nous avons modélisé ces essais. La comparaison qualitative des événements acoustiques et des énergies de fissuration issues respectivement des essais et des simulations numériques a montré des tendances similaires. Par contre, la comparaison d’un point de vue quantitatif a montré que le nombre des événements acoustiques numériques est plus important que celui déterminé expérimentalement. L’énergie dissipée par fissuration numérique est également plus importante que celle obtenue sur les essais. Cette différence est expliquée par le fait que les capteurs du dispositif expérimental ne détectent pas tous les événements acoustiques. Les résultats obtenus pourront permettre de mieux comprendre les cinétiques des phénomènes dynamiques dans les ouvrages souterrains profonds. Une autre application a consisté à modéliser un pilier de mine de fer de Joeuf (Lorraine). Le modèle numérique montre deux modes de fissuration dans le pilier : (a) écaillage en peau de pilier, (b) deux bandes de rupture s’initiant du mur et toit du pilier pour se propager vers le cœur du pilier. Ce travail offre de bonnes perspectives pour mieux comprendre la propagation de la fissuration à plus grande échelle mais aussi de progresser dans la recherche de corrélation entre la géomécanique et la géophysique / The study of rock mass behavior requires the understanding of their response under various loadings. The study of rock damage from an energetic point of view is essential in order to predict dynamic phenomena. These phenomena are due to the development of cracks in rocks subjected to strong initial and induced stresses. Fracturing is a form of energy dissipation that restores the balance of the involved medium. The aim of the thesis is to model rock cracks and study the behavior of underground structures at great depths. The development of models able to simulate the fracturing, the coalescence of cracks and their interaction with pre-existing fractures is essential. In the literature, there are two main theoretical and numerical approaches for crack modeling: continuous and discrete. A detailed analysis of these approaches has led us to choose the discrete approach and more particularly the code Yade. This code enables to simulate explicitly cracks propagation with or without pre-existing fractures. Developments have been made to evaluate the different forms of energy involved in rock behavior. In particular, a correlation between the cracks energy determined numerically and the microseismic activity observed in laboratory samples has been performed. The various energy components developed and then implemented in Yade are: external work, potential energy, elastic energy, friction energy, cracks energy, kinetic energy and damping energy. Validation of the energy approach was carried out by simulating laboratory tests. The evolution of the various energy components permits to verify that the energy balance is correctly evaluated. The energy balance was also verified at a structure scale by simulating the underground excavation of a Mine-by Experiment (URL Manitoba). The extension of the damaged zone induced by excavation and predicted by numerical simulations was compared with that observed in-situ around the Mine-by Experiment. It has been found that the predicted and the observed damage are similar in the directions of initial minor and major initial stresses. In addition, the energy formulation enables to study numerically the fracturing process of rocks. Wassermann (2006) performed uniaxial and triaxial compression tests on samples of iron ore from Lorraine. We have modeled these tests. The qualitative comparison of acoustic events and cracks energies determined from tests and numerical simulations showed similar trends. On the other hand, the quantitative comparison showed that the number of numerical acoustic events is greater than the number of experimental acoustic events. Also, the energy dissipated by cracks determined numerically is greater than the energy measured in the tests. This difference is explained by sensors accuracy of the experimental device, which are not able to detect all the generated acoustic events. The results obtained will allow us to better understand the dynamic phenomena in the deep underground structures. Another application consisted in modeling an iron ore pillar of Joeuf (Lorraine). The numerical model shows two modes of cracking in the pillar: (a) flaking of pillar wall, (b) two breaking bands initiating from the wall and the roof of the pillar to propagate towards his core. This provides good perspectives for better understanding cracks propagation at a larger scale, also to progress in the understanding of the correlation between geomechanics and geophysics