Ce travail a porté sur l'effet de la présence d'un fluide interstitielle sur le comportement des roches granulaires saturées dans la phase pré et post localisation. A partir d'essais de compression triaxiale et isotrope sur un grés de Fontainebleau d'une porosité de 20%, nous avons caractérisé à la fois le comportement en conditions drainés et non drainés de cette roche. Les résultats issus des essais isotropes ont mis en évidence l'influence de la pression de confinement et de la surpression interstitielle sur les paramètres poroélastque de la roche. La notion de la contrainte effective de Biot ainsi que la détermination analytique du coefficient de Skempton sont discutées. En ce qui concerne les essais de compression triaxiale, une attention particulière a été portée sur l'évolution de la surpression interstitielle issue des essais non drainés et le phénomène de localisation dans les deux conditions de drainage (drainé/ non drainés). Nous sommes parvenus à montrer par exemple que le seuil de localisation ne semble pas être affecté par la condition de drainage. Le résultas des essais triaxiaux ainsi que la prédiction de la localisation des déformations sont bien reproduits à l'aide d'un modèle de comportement avec élasticité non-linéaire et plasticité couplée à l'endommagement. L'étude de la localisation a été complétée par des observations microscopiques des bandes de cisaillement. L'analyse de la microstructure des bandes de cisaillement, par analyse d'images obtenues au microscope électronique à balayage, nous a permis d'étudier l'évolution de la porosité, de la surface spécifique et de la taille des grains à l'intérieur de ces bandes. En particulier, les changements de porosité et la structure des bandes en zones contractantes et/ou dilatantes sont discutés en fonction de la contrainte de confinement et le rôle de la surpression interstitielle est discuté. Nous présentons également une évaluation quantitative de l'évolution de la perméabilité en fonction des modifications de la microstructure en utilisant le modèle de Walsh et Brace (1984). Ce document se termine par une étude consacrée à l'étude de bande de cisaillement à grande échelle que constituent les failles, ou une analogie est observée entre la structure des bandes de cisaillement dans les granulaires obtenues sur le grès de Fontainebleau et la structure des zones de failles. Les travaux présentés ici s'inscrivent dans le cadre de notre participation au projet européen DGLAB-CORINTH sur les risques sismique dans le Golfe de Corinth.

L'objectif technique majeur de cette thèse se rapporte au développement d'un outilde simulation numérique fiable et robuste adapté à la maîtrise et la connaissance dufonctionnement mécanique des ouvrages géotechniques. Cet outil de simulation doitnotamment permettre une maîtrise des scénarii les plus pénalisants pour leur tenue,notamment sous l'aléa sismique.La qualité des résultats d'un tel outil dans le cadre de la méthode des éléments finisest fonction du modèle de comportement utilisé, de la qualité de l'intégration dumodèle et de sa résolution numérique. Dans le cadre de ce travail, le modèleélastoplastique de comportement cyclique des sols de l'ECP (dit de Hujeux) a étéintroduit dans Code_Aster selon un schéma d'intégration implicite, permettant ainsiune représentation fine et précise des phénomènes mis en jeu durant l'aléasismique. Ce modèle possède également l'avantage d'être adapté au comportementde différents types de sols dans le cadre des milieux poreux sous l'hypothèse despetites déformations. Après validation sur des chemins de chargement variés, lemodèle est à présent utilisé pour la simulation de la construction par couches, de lamise en eau et de la tenue sismique de barrages en terre.Le caractère non standard et adoucissant de ce modèle conduit à mettre en oeuvredes techniques de régularisation pour résoudre le problème de dépendancepathologique des résultats aux maillages lors de l'apparition de modes de ruinelocalisée. Le modèle de second gradient de dilatation est utilisé en complément aumodèle de comportement de Hujeux pour contrôler la largeur des bandes delocalisation apparaissant sur les structures étudiées. La prise en compte d'unecinématique enrichie permet de rendre objectives aux maillages les réponses desstructures durant leur phase d'adoucissement mais n'instaure pas l'unicité dessolutions aux problèmes posés suite aux instabilités. Dans le cadre d'essaisbiaxiaux drainés sur des matériaux dilatants exprimés au sein des milieux dusecond gradient de dilatation, il est apparu une dépendance de la largeur desbandes de cisaillement à l'état de contraintes initial en plus des propriétés desmatériaux.A partir d'un exemple analytique d'une bande de matériau cisaillée, cettedépendance a pu être exprimée, comparée et maîtrisée en fonction des paramètresde régularisation par rapport aux résultats des simulations numériques. L'extensionde cette approche à un cas de stabilité de pente sous chargement d'une fondationsuperficielle a ensuite été entrepris. La dépendance des largeurs de bandes à l'étatde contraintes initial est apparue comme un élément clé de la maîtrise du couplageentre le modèle de second gradient de dilatation et les modèles de type Cam-Clay

Cette étude a pour objectifs la caractérisation expérimentale des frottements saccadés dans les matériaux granulaires modèles constitués des billes de verre monodisperses en compression triaxiale drainée et l'identification des paramètres de contrôle. Cinq paramètres macroscopiques caractérisent ces frottements saccadés : la chute de déviateur et la contraction volumique, l'intermittence de déformation axiale, le module d'Young et le coefficient de Poisson. Les frottements saccadés affectent simultanément le déviateur et la déformation volumique. Le comportement macroscopique est globalement contractant tandis que le matériau tend vers un état limite critique en grandes déformations, à la manière des sables lâches. Cependant il présente localement, dans les phases de blocage qui suivent immédiatement les ruptures temporaires, le comportement dilatant des sables denses, qui obéit à une relation contrainte-dilatance linéaire et unique. Les frottements saccadés disparaissent au-delà d'une vitesse critique d'écrasement axial, qui dépend de la contrainte de confinement et de la taille des grains. Le module d'Young dynamique par propagation d'ondes varie avec la contrainte de confinement selon une loi de puissance. Le module d'Young quasi-élastique au départ des phases de blocage est constant à l'intérieur du domaine élastique, de même que le coefficient de Poisson. Le suivi par granulométrie laser et par analyse d'images des matériaux après un ou plusieurs essais triaxiaux permet de suivre l'évolution de la taille moyenne et de la forme des grains. Tandis que les instabilités par saccade disparaissent suite à un certain nombre d'essais, on observe, simultanément à une légère diminution du volume moyen, l'apparition progressive de populations d'objets non sphériques par une fusion des grains analogue au frittage.

L'objectif technique majeur de cette thèse se rapporte au développement d'un outilde simulation numérique fiable et robuste adapté à la maîtrise et la connaissance dufonctionnement mécanique des ouvrages géotechniques. Cet outil de simulation doitnotamment permettre une maîtrise des scénarii les plus pénalisants pour leur tenue,notamment sous l'aléa sismique.La qualité des résultats d'un tel outil dans le cadre de la méthode des éléments finisest fonction du modèle de comportement utilisé, de la qualité de l'intégration dumodèle et de sa résolution numérique. Dans le cadre de ce travail, le modèleélastoplastique de comportement cyclique des sols de l'ECP (dit de Hujeux) a étéintroduit dans Code_Aster selon un schéma d'intégration implicite, permettant ainsiune représentation fine et précise des phénomènes mis en jeu durant l'aléasismique. Ce modèle possède également l'avantage d'être adapté au comportementde différents types de sols dans le cadre des milieux poreux sous l'hypothèse despetites déformations. Après validation sur des chemins de chargement variés, lemodèle est à présent utilisé pour la simulation de la construction par couches, de lamise en eau et de la tenue sismique de barrages en terre.Le caractère non standard et adoucissant de ce modèle conduit à mettre en oeuvredes techniques de régularisation pour résoudre le problème de dépendancepathologique des résultats aux maillages lors de l'apparition de modes de ruinelocalisée. Le modèle de second gradient de dilatation est utilisé en complément aumodèle de comportement de Hujeux pour contrôler la largeur des bandes delocalisation apparaissant sur les structures étudiées. La prise en compte d'unecinématique enrichie permet de rendre objectives aux maillages les réponses desstructures durant leur phase d'adoucissement mais n'instaure pas l'unicité dessolutions aux problèmes posés suite aux instabilités. Dans le cadre d'essaisbiaxiaux drainés sur des matériaux dilatants exprimés au sein des milieux dusecond gradient de dilatation, il est apparu une dépendance de la largeur desbandes de cisaillement à l'état de contraintes initial en plus des propriétés desmatériaux.A partir d'un exemple analytique d'une bande de matériau cisaillée, cettedépendance a pu être exprimée, comparée et maîtrisée en fonction des paramètresde régularisation par rapport aux résultats des simulations numériques. L'extensionde cette approche à un cas de stabilité de pente sous chargement d'une fondationsuperficielle a ensuite été entrepris. La dépendance des largeurs de bandes à l'étatde contraintes initial est apparue comme un élément clé de la maîtrise du couplageentre le modèle de second gradient de dilatation et les modèles de type Cam-Clay / The main technical objective of this PhD thesis deals with the development of a soilbehavior numerical tool. It should be robust, efficient and adapted to model themechanical behavior of geotechnical structures (e.g. embankment dam) under theworst loading scenarii such earthquakes.In the finite element method, the quality results is directly linked to the soilconstitutive model, the integration scheme and the numerical resolution. In this PhDThesis, the ECP elastoplastic soil model is introduced in Code_Aster through animplicit scheme. An implicit scheme ensures to respect the theoretical formulation ofthe model. The ECP constitutive model is one of the models available in theliterature to represent the behaviour of different kinds of soils under cyclic loadingsand it is used since the 80's by hydraulic engineers at EDF. It is expressed in termsof effective stresses and infinitesimal strains. The developments are validated forlaboratory tests in a large scale of loading paths. On the other hand, the study of asand embankment was performed and compared to the results obtained with thefinite element software GEFDyn developed at ECP.The ECP model is based on a non-associated flow rule and it is able to reproduce asoftening behavior. When shear bands occur in the structure, these properties leadto a pathological sensitivity of the results depending on the mesh size. Therefore, aregularization technique has to be used to circumvent this problem and to obtainobjective results with respect to the mesh. The second gradient of dilation model isthus chosen to be coupled to the ECP model and in this way, to ensure a spatiotemporalindependence of results. However, this mesh independence still evolves ina potential domain of solutions, when instability occurs. The simulations of drainedbiaxial tests on laboratory samples show a dependence of shear bands thickness inregard to the initial stress state and material properties.An analytical problem of a dilatant shear band is used to extract the key factors.These theoretical solutions are compared and validated to numerical responses,which are in good accordance. A bearing capacity problem was also solved todemonstrate the potential of the method. The conclusion of this work establishes themain role of initial stress state over the shear band thickness in the context of thesecond gradient of dilation model and the models based on Cam-Clay approach.

Cette recherche vise à analyser et à identifier, par des essais, le comportement au fluage d'une argile de type kaolinite sous un état de contraintes triaxial. Une investigation fine de l’état microstructural, avant et après fluage, a par ailleurs été réalisée. Les résultats macroscopiques ont permis de constater que des phénomènes de dilatance et de contractante s’amplifient pendant le fluage. Les résultats montrent également que la position, sur le plan des contraintes (p'-q), du niveau de contrainte vis-à-vis des différents domaines volumétriques, contrôle le développement de la déformation au fluage et du taux de déformation. Les résultats de MEB ont montré que l'évolution microstructurale de l'argile, après le chargement mécanique, dépend fortement de l'histoire du chargement. / The objective of this research is to analyze the creep behavior of a typical clay along triaxial tests, moreover, to analyze the microstructural mechanisms of creep. Analysis on the macroscopic results ascertained that both dilatancy and contractancy phenomena could occur during creep. The magnitude of the dilatancy/contractancy during creep was guided by the test conditions: stress level and the over consolidation ratio, which specifically governed the direction of the volumetric strain variations. The results of SEM indicated that the microstructural evolution of the clay after the mechanical loading depend on the stress history. Afterwards, the structural evolution in creep phase depends on the structural pattern developed in monotonic loading

Les argiles raides sont généralement rencontrées à grande profondeur (plus de 300 mètres). Ces matériaux se caractérisent par leur teneur en eau relativement faible. De nombreux projets industriels justifient l'intérêt récent porté à ces matériaux. La thèse traite en particulier d'argiles raides susceptibles de répondre à des sollicitations par des déformations élastiques, plastiques et viscoplastiques. C'est la prise en compte simultanée d'une rupture engendrant des déformations plastiques et des déformations différées de types viscoplastiques qui font l'originalité de cette thèse. Cette nouvelle approche constitue une amélioration indispensable dans la description du comportement mécanique d'argiles raides. Ce travail est présenté en trois parties : la première partie est consacrée à l'étude expérimentale et à la modélisation du comportement mécanique des argiles raides. Dans cette partie, nous décrivons les matériaux considérés, et les essais réalisés. Nous discutons les résultats et élaborons un modèle viscoplastique avec rupture. La deuxième partie, « Application du modèle VPR au calcul d'un tunnel de section circulaire », est consacrée à l'élaboration d'une solution originale semi analytique et d'un algorithme implanté dans GEOMEC91. La troisième partie, « Application du modèle VPR au calcul 2d axisymétrique des tunnels soutenus et étude paramétrique », montre très clairement l'influence du modèle VPR sur la convergence du tunnel au moment de la pose du soutènement, et par conséquent sur les sollicitations de ce dernier. Les résultats des calculs montrent une forte influence de la cohésion à court terme sur la convergence du tunnel.

Cette étude a pour objectifs la caractérisation expérimentale des frottements saccadés dans les matériaux granulaires modèles constitués des billes de verre monodisperses en compression triaxiale drainée et l'identification des paramètres de contrôle. Cinq paramètres macroscopiques caractérisent ces frottements saccadés : la chute de déviateur et la contraction volumique, l'intermittence de déformation axiale, le module d'Young et le coefficient de Poisson. Les frottements saccadés affectent simultanément le déviateur et la déformation volumique. Le comportement macroscopique est globalement contractant tandis que le matériau tend vers un état limite critique en grandes déformations, à la manière des sables lâches. Cependant il présente localement, dans les phases de blocage qui suivent immédiatement les ruptures temporaires, le comportement dilatant des sables denses, qui obéit à une relation contrainte-dilatance linéaire et unique. Les frottements saccadés disparaissent au-delà d'une vitesse critique d'écrasement axial, qui dépend de la contrainte de confinement et de la taille des grains. Le module d'Young dynamique par propagation d'ondes varie avec la contrainte de confinement selon une loi de puissance. Le module d'Young quasi-élastique au départ des phases de blocage est constant à l'intérieur du domaine élastique, de même que le coefficient de Poisson. Le suivi par granulométrie laser et par analyse d'images des matériaux après un ou plusieurs essais triaxiaux permet de suivre l'évolution de la taille moyenne et de la forme des grains. Tandis que les instabilités par saccade disparaissent suite à un certain nombre d'essais, on observe, simultanément à une légère diminution du volume moyen, l'apparition progressive de populations d'objets non sphériques par une fusion des grains analogue au frittage. / The objectives of this study are the experimental characterisatino of the stick-slip instabilities in a model granular material and the identification of relevant control parameters. As monodisperse glass beads are subjected to drained triaxial compression tests, five macroscopic parameters characterize the stick-slip phenomenon: the deviator drop, the jump in volumetric contraction, the intermittence of the axial strain, Young’s modulus and Poisson’s ratio. The stick-slip events simultaneously affect the deviator stress and the volumetric strain. While the global material behavior is that of a loose sand, gradually contracting and hardening as it approaches its large strain critical state, its response in the “stick” phases immediately following the “slip” instabilities is similar to that of dense, dilatant sands, with a unique, linear stress-dilatancy relationship. Stick-slip events disappear beyond a critical axial strain rate, depending on the confining stress and on the grain diameter. The Young modulus associated to wave propagation varies with the confining stress according to a power law. The quasi-elastic modulus measured at the beginning of the stick phase is constant inside the elastic domain, as well as the Poisson ratio. The evolution of grain size and shape after one or several triaxial tests is monitored by laser granulometry and image analysis. The gradual vanishing of stick-slip events, on repeating the tests, is likely related to the global decreasing trend of average particle volume and to the formation of non-spherical objects, apparently by some phenomenon analogous to sintering.