Ce travail présente une nouvelle méthode de prédiction de l'effondrement d'un terrain minier suite à une exploitation intensive. L'approche développée a été appliquée pour le cas de la concession de Drouville, appartenant à la Compagnie des Salins du Midi et des Salines de l'Est.<br/>Le gisement y était exploité par cavités de dissolution interconnectées par facturation hydraulique. Des développements incontrôlés de la dissolution ont abouti à la jonction de deux cavités voisines et à la formation d'une grande cavité SG4/5 dont l'effondrement a été provoqué pour des raisons de sécurité.<br/>La rétro-analyse de cet effondrement a été le point de départ de la recherche qui est présentée, dont l'objectif final était la mise au point d'une méthode d'exploitation sûre permettant de récupérer la plus grande partie du gisement sans risquer d'évolution ultérieure.<br/>Afin d'atteindre det objectif, nous avons utilisé la méthodes des éléments finis pour résoudre le problème mécanique.<br/>Le modèle développé intègre des lois rhéologiques décrivant le comportement des matériaux rencontrés. Les paramètres de ces lois ont été obtenus par rétro-analyse directe (méthode d'essai et d'erreur / méthode de minimisation d'erreur).<br/>L'application du modèle pour des simulations bidimensionnelles et tridimensionnelles nous a permis d'identifier les zones de forte plastification et de rupture. De plus, il nous a été possible d'anticiper la réaction du massif rocheux durant le processus d'effondrement.<br/>En utilisant cette technique conjointement avec une analyse détaillée du comportement mécanique du massif, une nouvelle méthode d'exploitation intensive par élargissement progressif de la zone d'effondrement a été proposée.

A travers cette étude, on s'attache à montrer que les saumures contenues dans les roches halitiques des formations salifères du Bassin bressan, apparaissent en quantités non négligeables et peuvent jouer un rôle important dans les mécanismes de déformation de ces roches, ceci tant à l'échelle humaine d'un projet souterrain que lors de la formation des dômes ou diapirs. Dans la première partie, au cours de l'étude pétrographique et pétrophysique sur des carottes provenant de sondages profonds, nous analysons l'espace poreux des principaux faciès des roches halitiques en adaptant les méthodes de mesure classiques. Ceci permet d'estimer la teneur in situ en saumures intercristallines. On caractérise ensuite les inclusions fluides dans les cristaux de halite et une méthode de quantification de la teneur en saumures intracristallines est développée. La deuxième partie de l'étude concerne l'analyse du rôle potentiel de ces saumures dans les mécanismes de déformation. Une synthèse bibliographique est effectuée, puis certains aspects sont précisés à l'aide d'essais expérimentaux réalisés avec une microcellule de fluage permettant l'observation au microscope, en continu, d'une lame épaisse de matériau sous contrainte et en lui imposant un gradient de température. A l'échelle humaine, on montre notamment que la déformation cataclastique peut dans certains cas être prépondérante et peut fortement augmenter la perméabilité. L'ajout de saumures ou l'éclatement thermique des inclusions fluides peut faciliter la rupture de la roche. La thermomigration des inclusions fluides intracristallines a également été analysée. A l'échelle des temps géologiques, on conforte, par des observations et des essais, l'hypothèse du rôle essentiel joué par les mécanismes de dissolution-cristallisation lorsqu'ils s'associent aux mécanismes de déformation plastique pour provoquer la recristallisation dynamique des roches halitiques.

Ce travail est une contribution à la caractérisation de l'espace poreux des roches sédimentaires, par l'étude des équilibres capillaires dans une vingtaine d'échantillons provenant de divers horizons géologiques et géographiques; ce en relation avec une connaissance fine de la texture de ces roches. Les caractéristiques pétrographiques des roches étudiées ont été définies grâce à la combinaison de divers moyens classiques : observation visuelle, lames minces, photographies au microscope électronique à balayage (MEB) et diffraction des rayons X. Quant aux caractéristiques pétrophysiques, 3 ont été mesurées : il s'agit de la porosité, de la perméabilité, et du facteur de formation. L'analyse de la texture a été envisagée à partir de deux techniques basées sur l'étude de la répartition d'un couple de fluides non miscibles (air-eau et mercure-vapeur) dans les réseaux poreux. Il s'agit de la méthode des états restaurés et de porosimétrie au mercure. Elles permettent l'obtention de courbes de pression capillaire en fonction de la saturation (courbe Pf et courbe porosimétrique). De ces courbes on déduit plusieurs paramètres qui sont en relation avec la géométrie de l'espace poreux, tel que la décomposition de la porosité en deux classes : la macroporosité (rayon d'accès > 1µm) et la microporosité (rayon d'accès < 1µm), ainsi que le rayon moyen d'accès de chacune d'elles. Les échantillons ont été ensuite répartis dans deux catégories : les monoporeux et les biporeux. La combinaison des mesures pétrophysiques et des paramètres issus des courbes de pression capillaire a permis l'obtention de nombreuses corrélations entre les différents facteurs. La confrontation entre l'observation des échantillons au MEB et les courbes de pression capillaire assimilées à des hyperboles d'équation Y = A/X a montré que le facteur A est en relation directe avec le mode de la répartition des accès des pores. D'autre part la comparaison entre les courbes Pf et porosimétrique montre que l'on peut déduire l'une de l'autre moyennant quelques manipulations complémentaires simples.

4 chapitres : 1. Quelques notions sur les modélisations de la houle. 2. Le problème élastique. 3. Stabilité d'un massif de sol cohérent. - 4. Stabilité d'un massif de sol granulaire.

Ce travail, contribue à mieux comprendre les phénomènes complexes lors des tirs de mine afin d'approcher une conception optimisée. La première partie, concernant la détonation, est destinée à la compréhension des modèles numériques fondés sur la théorie de la propagation de détonation. La seconde partie, concernant la modélisation de la fracturation des roches, s'adresse à la simulation des tirs. Une nouvelle loi de comportement dynamique des roches (type élastique avec endommagement anisotrope) est proposée et est intégrée dans RADIOSS, un code de calcul des problèmes dynamiques transitoires. Ce modèle permet également de calculer les dimensions représentatives des fragments ; les calculs avec diverses modélisations (en 2-D et en 3-D, détonation/structure) sont effectués. La dernière partie est consacrée à l'optimisation des tirs à l'aide des traitements numériques avancés. La connaissance des mécanismes de la fragmentation des roches est approfondie notamment sur l'influence de certains paramètres de conception. Les exemples sur la recherche des règles par l'apprentissage automatique montrent l'efficacité de cette approche. L'exemple sur la recherche d'une conception optimisée montre la puissance et la perspective pratique de l'approche par algorithmes génétiques couplés avec l'apprentissage automatique.

A l'opposé du dimensionnement de fondations d'ouvrages sur les sols, régi par des règles de calcul validées, le dimensionnement de fondations au rocher est mal maîtrisé. Jusqu'à présent, ceci n'a pas posé de problème majeur, car le dimensionnement de fondations se trouvait plutôt limité par la résistance du béton que par celle du massif rocheux. Mais la construction d'ouvrages d'art de plus en plus majestueux nécessite aujourd'hui une meilleure connaissance de ce domaine de la mécanique des roches. La première partie recense les méthodes les mieux adaptées pour déterminer la déformabilité et la résistance d'un massif rocheux. Elle présente les différentes méthodes de dimensionnement de fondations superficielle et semi-profonde au rocher. Enfin, elle analyse les textes réglementaires français et étrangers. Il existe peu de méthodes de dimensionnement de fondations qui tiennent compte du caractère discontinu des massifs rocheux ; le cas de fondations soumises à des efforts latéraux et à des moments renversants n'est quasiment pas traité. La deuxième partie consiste en l'analyse numérique du comportement de fondations superficielle et semi-profonde sur massif rocheux à l'aide du code de calcul par éléments distincts UDEC. Les résultats montrent que les discontinuités du massif rocheux ont un rôle primordial sur le comportement sous effort incliné de fondations semi-profondes. Selon leur pendage, les discontinuités peuvent diminuer considérablement la charge limite admissible. Quelques modèles analytiques simples à peu de blocs sont ensuite proposés, afin d'estimer - pour un dimensionnement préliminaire - la charge limite de rupture d'un puits marocain sous effort latéral. Les résultats numériques obtenus sont bidimensionnels; les résultats tridimensionnels peuvent être estimés à l'aide de coefficients de transfert 2D/3D, qu'il serait intéressant de vérifier à l'aide d'un code de calcul par éléments distincts tridimensionnel. Enfin, il serait nécessaire de valider les résultats de modélisation numérique sur un site réel.

Parmi les techniques de renforcement des sols par inclusions, les micropieux offrent un très vaste champ d'application pour les fondations, les reprises en sous-oeuvre et la stabilisation des pentes et talus. La présente recherche qui fait partie du projet national FOREVER (Renforcement des Sols par Micropieux) a pour objet de modéliser le comportement des micropieux isolés ou en groupe, plus particulièrement développer des outils numériques permettant le calcul. Pour l'étude du comportement des micropieux isolés sous charge axiale, en prenant en compte deux modèles pour la mobilisation du frottement latéral et la résistance de pointe (modèle basé sur l'utilisation des résultats pressiométriques MENARD et modèle basé sur le mécanisme de cisaillement pur d'anneaux concentriques avec une nouvelle définition pour le rayon d'influence), trois micropieux d'essai en vraie grandeur réalisé sur le chantier du TGV Paris-Bruxelles (1994) sont analysés. Un modèle, dit modèle "Multizones", a été présenté et développé pour l'analyse non linéaire des micropieux isolés et chargés latéralement. Ce modèle basé sur le concept de rotule plastique, calcule l'expansion de la zone plastifiée le long du micropieu. Il est introduit dans le programme PILATE-LCPC, donnant le nouveau programme "PILATEP". Les résultats de PILATEP sont comparés avec ceux obtenus par la méthode des éléments finis (programme CESAR-LCPC) ainsi qu'avec des résultats d'essais en modèle réduit par MAHMOUD et BURLY (1994). Un nouveau modèle "hybride" de type "facteur z" et "facteur y" pour le calcul des groupes et des réseaux de micropieux est présenté. Ce modèle combine, d'une part, la méthode des fonctions de transfert de charge pour analyser chacun des micropieux du groupe (modèles 't-z' et 'p-y'), et d'autre part, les concepts du continuum élastique pour déterminer l'interaction entre les micropieux du groupe, appelée "effet de groupe". Ce modèle est introduit dans le programme GOUPIL-LCPC, donnant naissance au programme "GOUPEG". Afin de valider ce logiciel, on compare les résultats de GOUPEG aux résultats expérimentaux du groupe de trois micropieux en vraie grandeur (à Rueil-Malmaison, 1994). Dans un réseau, les micropieux sont chargés axialement et latéralement. Dans ce cas, le couplage entre les efforts axiaux et latéraux est étudié. Enfin, les résultats de GOUPEG sont comparés, avec succès, à ceux obtenus par les méthodes du continuum ainsi que par la méthode des éléments finis (programme CESAR-LCPC).

Ce mémoire porte sur l'étude du comportement et de la rupture d'une roche (grès rouge de Wildmoor) de réservoir lors de l'exécution et de la mise en production d'un forage pétrolier. Divers chemins de sollicitations, en conditions axisymétriques (essais de compression uniaxiale et triaxiale et essais d'extension), ont permis de caractériser le comportement mécanique de la roche. Les essais avec cycles de charge-décharge ont permis de quantifier l'endommagement du matériau. La loi de comportement développée tient compte de la dépendance en contraintes des paramètres élastiques du matériau et de l'endommagement. Le comportement plastique de la roche est décrit par un critère de Mohr-Coulomb non-linéaire tenant compte de la dépendance en pression des paramètres plastiques (coefficients de frottement et de dilatance) de la roche. Il comporte une phase d'écrouissage isotrope en frottement, suivie d'une phase de radoucissement cinématique de la cohésion. Ce modèle de comportement est appliqué à l'étude de la localisation des déformations dans le cadre d'une théorie de plasticité non-coaxiale. L'application à l'étude de stabilité d'un forage pétrolier a comporté le développement d'un programme de calcul du champ de contraintes autour d'un puits en cours de creusement et de production. L'étude des conditions d'apparition de différents modes d'instabilités (formation de bandes de cisaillement et instabilité de surface) a mis en évidence l'influence de la dilatance de la roche, de la profondeur du forage et du gradient de pression lors de la mise en production du puits.

L'objet de cette recherche est de modéliser le comportement Thermo-Hydro-Mécanique des massifs rocheux fracturés en vue du stockage de déchets radioactifs. A cette fin, nous avons été amené à proposer une méthodologie de modélisation des massifs rocheux fracturés qui a été appliquée à un site souterrain réel (Site EDF de Nouvelle Romanche). Cette méthodologie consiste, dans un premier temps, à rechercher les VER hydraulique et mécanique. Au dessus du plus grand de ces VER, les développements effectués dans un programme d'éléments finis permettent de modéliser toutes les fractures de manière explicite. On détermine les propriétés mécaniques homogénéisées en condition drainée et non drainée en simulant des essais triaxiaux représentatifs du massif soumis à un chargement. Ces simulations permettent d'analyser l'évolution des propriétés hydrauliques homogénéisées avec l'état de contraintes. La réalisation des expériences de type drainé et non drainé permet de discuter de la pertinence de l'assimilation d'un milieu fracturé à un milieu poreux. Ces simulations ont mis en évidence, en particulier, le rôle du comportement en cisaillement des fractures sur les propriétés homogénéisées mécaniques et hydrauliques. Du point de vue de la thermique, tant que le mode principal de propagation de la chaleur est la conduction, les propriétés thermiques homogénéisées sont peu différentes des propriétés thermiques de la matrice rocheuse.

Dans les massifs rocheux fissurés, il est très difficile de réaliser des ouvrages sans avoir fait, au préalable, une étude de sa fracturation et de son comportement mécanique et hydraulique. L'étude d'un massif rocheux fissuré doit comporter les étapes suivantes : - Relevé systématique des discontinuités du massif, - Étude statistique de ses discontinuités pour identifier les différentes familles qui découpent le massif, - Étude du comportement mécanique des discontinuités, - Étude générale de stabilité, - Étude hydraulique. Notre recherche est relative à l'étude générale de stabilité. Les méthodes classiques de calcul de stabilité, basées sur l'analyse des forces en situation d'équilibre limite et ne prenant pas en compte la déformabilité du massif et la répartition des contraintes existantes dans celui-ci, ne donnent pas entière satisfaction. La méthode que nous avons développée est une méthode dite par "éléments distincts". Elle permet de modéliser les massifs rocheux fissurés par un assemblage de blocs définis par le système de discontinuités et de calculer leurs interactions réciproques. L'interaction entre les blocs se transmet à travers leurs contacts (ou joints) suivant les lois de comportement mécanique de ces derniers. On peut exprimer les contraintes au niveau des joints en fonction de leurs rigidités et des déplacements des blocs. La minimisation de l'énergie potentielle totale de la structure étudiée permet de produire un système d'équations où sont exprimés les forces exercées sur le système et les déplacements des blocs par l'intermédiaire d'une matrice de rigidité. Il est souvent nécessaire de considérer des lois de comportement mécanique des joints non-linéaires, ce qui conduit à un système d'équations non-linéaires dont la résolution se fait d'une manière incrémentale itérative. Par cette méthode, on peut alors étudier les déformations causées par la réalisation d'une excavation dans le massif, rendre compte des ruptures locales qui peuvent évoluer suivant des mécanismes de ruptures progressives. Dans cette recherche, nous avons également modélisé les systèmes d'ancrages par boulons ou tirants actifs et par armatures passives scellées dans le massif. Le comportement mécanique des ancrages est considéré comme élasto-plastique avec écrouissage cinématique.