Le présent travail est consacré à l'étude du comportement thermo-hydro-mécanique linéaire et non linéaire des roches poreuses de type argilites par approche de changement d'échelle. A partir des observations microstructurales de ces matériaux, un modèle conceptuel a été proposé. Dans ce modèle, le volume élémentaire représentatif du milieu hétérogène est composé d'une phase matricielle argileuse contenant des inclusions sphériques de minéraux de quartz et de calcite et des inclusions ellipsoïdales aplaties représentant l'espace poreux. Dans un premier temps, le procédé de la modélisation a été exploité par la détermination des propriétés effectives isotropes et isotrope transverses des argilites : la conductivité thermique et les propriétés thermo-hydro-mécanique croisées. En outre, de nombreuses études numériques ont mis en évidence l'influence de la morphologie de l'espace poreux, de la minéralogie et des schémas d'estimation sur les résultats prédictifs. Dans un deuxième temps, nous avons modélisé le comportement mécanique non linéaire (élasto-plastique, élasto-viscoplastique) des roches argileuses. La comparaison entre les simulations numériques et les résultats expérimentaux disponibles (essai de compression triaxiale, essai de fluage) a confirmé la validation du modèle développé / The present work deals with the linear and non-linear thermo-hydro-mechanical behaviour of porous rocks such as the argillite by the multiscale modelling approach. Based on microstructure observations, a conceptual model was proposed. In this model, the representative elementary volume of a heterogeneous medium is composed of an argillaceous matrix containing spherical inclusions of minerals quartz and calcite and ellipsoidal inclusions representing the pore space. In a first step, the process of modelling has been exploited by determining the isotropic and transversely isotropic effective properties of the argillite: thermal conductivity and thermo-hydro-mechanical properties. Furthermore, many numerical studies have highlighted the influence of the morphology of the pore space, of the mineralogy and of the estimate schemes to the predictive results. In a second step, we modelled the non linear mechanical behaviour (elasto-plastic, elasto-viscoplastic) of argillaceous rocks. The comparison between numerical simulations and available experimental results (triaxial compression test, creep test) confirmed the validation of the model developed

Les structures en béton des centrales nucléaires peuvent être soumises à des contraintes thermo- hydriques modérées, caractérisées par des températures de l’ordre de la centaine de degrés aussi bien en conditions de service qu’accidentelles. Ces contraintes peuvent être à l’origine de désordres importants notamment la fissuration qui a pour effet d’accélérer les transferts hydriques dans la structure. Dans le cadre de l’étude de la durabilité de ces structures, le modèle THMs a été développé au Laboratoire d’Etude du Comportement des Bétons et des Argiles (LECBA) du CEA Saclay pour simuler le comportement du béton face à des sollicitations couplées thermo-hydro-mécaniques. Dans cette thèse on s’est intéressé à l’amélioration dans le modèle THMs d’une part de l’estimation des paramètres mécaniques et hydromécaniques du matériau en conditions partiellement saturées et en présence de fissuration et d’autre part de la description de la fissuration. La première partie a été consacrée à la mise au point d’un modèle basé sur une description multi-échelle de la microstructure des matériaux cimentaires, en partant de l’échelle des principaux hydrates (portlandite, ettringite, C-S-H, etc.) jusqu’à l’échelle macroscopique du matériau fissuré. Les paramètres investigués sont obtenus à chaque échelle de la description par des techniques d’homogénéisation analytiques. Dans la seconde partie on s’est attaché à décrire numériquement de façon précise la fissuration notamment en termes d’ouverture, de localisation et de propagation. Pour cela une méthode de réanalyse éléments finis/éléments discrets a été proposée et validée sur différents cas-test de chargement mécanique. Enfin la procédure a été mise en œuvre dans le cas d’un mur chauffé et une méthode de recalcul de la perméabilité a été proposée permettant de montrer l’intérêt de la prise en compte de l’anisotropie du tenseur de perméabilité lorsqu’on s’intéresse à l’étude des transferts de masse dans une structure en béton fissurée. Mots clés : matériaux cimentaires, homogénéisation, modélisation multi-échelle, microfissures, éléments discrets, éléments finis, chargements thermo-hydro-mécaniques. / The concrete structures of nuclear power plants can be subjected to moderate thermo-hydric loadings characterized by temperatures of the order of hundred of degrees in service conditions as well as in accidental ones. These loadings can be at the origin of important disorders, in particular cracking which accelerate hydric transfers in the structure. In the framework of the study of durability of these structures, a coupled thermo-hydro-mechanical model denoted THMs has been developed at Laboratoire d’Etude du Comportement des Bétons et des Argiles (LECBA) of CEA Saclay in order to perform simulations of the concrete behavior submitted to such loadings. In this work, we focus on the improvement in the model THMs in one hand of the assessment of the mechanical and hydromechanical parameters of the unsaturated microcracked material and in the other hand of the description of cracking in terms of opening and propagation. The first part is devoted to the development of a model based on a multi-scale description of cement-based materials starting from the scale of the main hydrated products (portlandite, ettringite, C-S-H etc.) to the macroscopic scale of the cracked material. The investigated parameters are obtained at each scale of the description by applying analytical homogenization techniques. The second part concerns a fine numerical description of cracking. To this end, we choose to use combined finite element and discrete element methods. This procedure is presented and illustrated through a series of mechanical tests in order to show the feasibility of the method and to proceed to its validation. Finally, we apply the procedure to a heated wall and the proposed method for estimating the permeability shows the interest to take into account an anisotropic permeability tensor when dealing with mass transfers in cracked concrete structures. Keywords : cement-based materials, homogenization, multi-scale modelling, microcracks, discrete elements, finite elements, thermo-hydro-mechanical loadings.

Le comportement des roches-réservoirs souterraines est un sujet important pour de nombreuses applications liées à la production d’énergie (extraction d’hydrocarbures, séquestration de CO2, ...). L'une des principales questions posées est celle de l'effet des déformations sur les propriétés de transfert hydraulique de la roche, en particulier en conditions saturées. En effet, la déformation des géomatériaux est rarement homogène en raison de conditions aux limites complexes et de sa tendance intrinsèque à se localiser. Cette non-uniformité spatiale de la déformation produit un champ de perméabilité hétérogène. Cela remet en question la validité (a) des méthodes traditionnelles d'analyse macroscopique et (b) des mesures établies principalement loin des zones de déformation localisée. Ainsi, pour améliorer la caractérisation des géo-matériaux, il est crucial d’avoir des mesures locales de la perméabilité, et de connaître la relation entre la déformation et la perméabilité, qui gouverne leur comportement hydraulique.Cette thèse porte sur l’étude du couplage hydromécanique des roches par tomographie aux neutrons et aux rayons X, ainsi que sur le développement de nouvelles méthodes d'analyse. Même si le recours à l'imagerie par rayons X en géosciences devient de plus en plus accessible, la détection directe des fluides a été très limitée en raison du faible contraste air/eau dans les géomatériaux. Contrairement aux rayons X, les neutrons sont très sensibles à l’hydrogène présent dans l'eau. La radiographie par neutrons permet donc d'obtenir des images où la détection du fluide est bien plus facile. De plus, les neutrons sont sensibles aux isotopes, ce qui veut dire que l’eau lourde et celle normale, qui ont des propriétés physico-chimiques proches, peuvent être distinguées avec une grande précision. Il faut noter que l’imagerie aux neutrons pour les roches est un domaine expérimental qui est essentiellement inexploré, ou limité à des études 2D d'échantillons secs, avec peu ou pas de contrôle sur les conditions aux limites.Dans le cadre de ce travail, nous avons conçu une nouvelle cellule triaxiale, avec un contrôle asservi, pour effectuer des expériences d'écoulement de fluides multiples dans un échantillon de roche saturé et chargé mécaniquement avec acquisition des données neutroniques en 4D. Une autre originalité du projet est l'utilisation d'installations d'imagerie neutroniques à haute performance (CONRAD-2 au Helmholtz Zentrum à Berlin et NeXT à l'Institut Laue-Langevin à Grenoble), profitant de la technologie de pointe et des lignes de faisceaux les plus puissantes du monde. Cela a permis d'acquérir des données à une fréquence optimale pour notre étude.Ce travail présente les résultats de plusieurs campagnes expérimentales couvrant une série de conditions initiales et de conditions aux limites relativement complexes. Pour quantifier le couplage hydromécanique local, nous avons appliqué un certain nombre de procédures de post-traitement standard et nous avons également développé un ensemble de méthodes de mesure, par exemple pour suivre le front d’eau et déterminer les cartes de vitesse 3D. Les résultats montrent que la vitesse du front d’eau entraîné par imbibition dans un échantillon sec est augmentée à l’intérieur d'une bande de cisaillement compactante, tandis que la vitesse d’écoulement entraîné par la pression est réduite dans les échantillons saturés, quelque soit la réponse volumétrique de la bande de cisaillement (compactante / dilatante). La nature des données 3D et des analyses s'est révélée essentielle dans la caractérisation du comportement mécanique complexe des échantillons et de la vitesse d'écoulement qui en résulte.Les résultats expérimentaux obtenus contribuent à la compréhension de l'écoulement dans les matériaux poreux sous chargement, garantissent la pertinence de l'analyse et permettent d’etablir une méthode expérimentale pour d'autres campagnes hydromécaniques in-situ. / The behaviour of subsurface-reservoir porous rocks is a central topic in resource engineering industry and has relevant applications for hydrocarbon and water production or CO2 sequestration. One of the key open issues is the effect of deformations on the hydraulic properties of the host rock, specifically in saturated environments. Deformation in geomaterials is rarely homogeneous because of the complex boundary conditions they undergo as well as for their intrinsic tendency to localise. This non uniformity of the deformation yields a non uniform permeability field, meaning that the traditional macroscopic analysis methods are outside their domain of validity. These methods are in fact based on measurements taken at the boundaries of a tested sample, under the assumption of internal homogeneity. At this stage, our understanding is in need of direct measurements of the local fluid permeability and its relationship with localiseddeformation.This doctoral dissertation focuses on the acquisition of such local data about the hydro mechanical properties of porous geomaterials in full-field, adopting neutron and x-ray tomography, as well as on the development of novel analysis methods. While x-ray imaging has been increasingly used in geo-sciences in the last few decades, the direct detection of fluid has been very limited because of the low air/water contrast within geomaterials. Unlike x-rays, neutrons are very sensitive to the hydrogen in the water because of their interaction with matter (neutrons interact with the atoms’ nuclei rather than with the external electron shell as x-rays do). This greater sensitivity to hydrogen provides a high contrast compared to the rock matrix, in neutron tomography images that facilitates the detection of hydrogen-rich fluids. Furthermore, neutrons are isotope-sensitive, meaning that water (H 2 0) and heavy water (D20), while chemically and hydraulically almost identical, can be easily distinguished in neutron imaging.The use of neutron imaging to investigate the hydromechanical properties of rocks is a substantially under-explored experimental area, mostly limited to 2D studies of dry, intact or pre-deformed samples, with little control of the boundary conditions. In thiswork we developed a new servocontrolled triaxial cell to perform multi-fluid flow experiments in saturated porous media, while performing in-situ loading and acquiring 4-dimensional neutron data.Another peculiarity of the project is the use of high-performance neutron imaging facilities (CONRAD-2, in Helmholtz Zentrum Berlin, and NeXT-Grenoble, in Institut Laue-Langevin), taking advantage of the world’s highest flux and cutting edge technology to acquire data at an optimal frequency for the study of this processes. The results of multiple experimental campaigns covering a series of initial and boundary conditions of increasing complexity are presented in this work.To quantify the local hydro-mechanical coupling, we applied a number of standard postprocessing procedures (reconstruction, denoising, Digital Volume Correlation) but also developed an array of bespoke methods, for example to track the water front andcalculate the 3D speed maps.The experimental campaigns performed show that the speed of the water front driven by imbibition in a dry sample is increased within a compactant shear band, while the pressure driven flow speed is decreased in saturated samples, regardless of the volumetric response of the shear band (compactant/dilatant). The 3D nature of the data and analyses has revealed essential in the characterization of the complex mechanical behaviour of the samples and the resultant flow speed.The experimental results obtained contribute to the understanding of flow in porous materials, ensure the suitability of the analysis and set an experimental method for further in-situ hydromechanical campaigns.

Cette thèse présente une plateforme d'assimilation de données issues de l'imagerie vidéo et intégrée au modèle numérique d'évolution de profil de plage 1DBEACH. Le manque de jeux de données bathymétriques haute-fréquence est un des problèmes récurrents pour la modélisation morphodynamique littorale. Pourtant, des relevés topographiques réguliers sont nécessaires non seulement pour la validation de nos modèles hydro-sédimentaires mais aussi dans une perspective de prévision d'évolution morphologique de nos plages sableuses et d'évolution de la dynamique des courants de baïnes en temps réel. Les récents progrès dans le domaine de l'imagerie vidéo littorale ont permis d'envisager un moyen de suivi morphologique quasi-quotidien et bien moins coûteux que les traditionnelles campagnes de mesure. En effet, les images dérivées de la vidéo de type timex ou timestack rendent possible l'extraction de proxys bathymétriques qui permettent de caractériser et de reconstruire la morphologie de plage sous-jacente. Cependant, ces méthodes d'inversion bathymétrique directes sont limitées au cas linéaire et nécessitent, selon les conditions hydrodynamiques ambiantes, l'acquisition de données vidéo sur plusieurs heures voire plusieurs jours pour caractériser un état de plage. En réponse à ces différents points bloquants, ces travaux de thèse proposaient l'implémentation puis la validation de méthodes d'inversion bathymétrique basées sur l'assimilation dans notre modèle de différentes sources d'observations vidéo disponibles et complémentaires. A partir d'informations hétérogènes et non redondantes, ces méthodes permettent la reconstruction rapide et précise d'une morphologie de plage dans son intégralité pour ainsi bénéficier de relevés bathymétriques haute fréquence réguliers. / This thesis presents data-model assimilation techniques using video-derived beach information to improve the modelling of beach profile evolution.The acquisition of accurate and recurrent nearshore bathymetric data is a difficult and challenging task which limits our understanding of nearshore morphological changes. This is particularly true in the surf zone which exhibits the largest degree of morphological variability. In addition, surfzone bathymetric data are crucial from many perspectives such as numerical model validation, operational rip current prediction or real-time nearshore evolution modelling. In parallel, video imagery recently arose as a low-cost alternative to direct measurement in order to daily monitor beach morphological changes. Indeed, bathymetry proxies can be extracted from video-derived images such as timex or timestacks. These data can be then used to estimate underlying beach morphologies. However, simple linear depth inversion techniques still suffer from some restrictions and require up to a 3-day dataset to completely characterize a given beach morphology. As an alternative, this thesis presents and validates data-assimilation methods that combine multiple sources of available video-derived bathymetry proxies to provide a rapid, complete and accurate estimation of the underlying bathymetry and prevent from excessive information.

Lorsqu'un hélicoptère vol au dessus de l'eau le constructeur se doit d'anticiper l'amerrissage d'urgence (ditching) pour la sécurité de l'équipage. La solution adoptée consiste en des systèmes de flottabilités gonflables. Ces systèmes pour EUROCOPTER ont deux buts : amortir l'impact et assurer la flottabilité le temps que l'appareil soit évacué. Des normes définissent les limites du ditching en terme de vitesses d'impact et en terme d'état de mer. Lors du dimensionnement de l'appareil et des flottabilités pour l'impact des centaines de calculs sont a réitérer pour converger vers la solution optimale.Cette thèse s'inscrit dans un projet d'optimisation d'un logiciel maison voué au Ditching. Elle fait suite à celle de N.Malleron(2009) focalisée sur le cas vertical. Elle apporte des contributions à la compréhension de cette IFS complexe via différentes approches expérimentales et l'exploration de différents modèles que l'on tente de corréler avec les résultats expérimentaux. / When operating on offshore conditions, helicopters may encounteremergency situations during which an immediate sealanding isrequired. In this case, the aircraft may experience a water impact event.Floatability systems are inflated to cushion the impact and toensure that the floatation is maintained while the crew evacuates.The dimensioning procedure of helicopter's floatability systemsremains very challenging. In this process, accurate prediction of thebehavior and the global loads on the helicopters floaters during ditchingis crucial. To obtain a better understanding of the floatability system,and to be able to develop more accurate models, the hydrodynamicalcharacteristics of the floaters, as flexible bodies, need to be furtherinvestigated. Experiments assessing the hydrodynamics characteristicsof these particular types of deformable bodies are rare. To this aim, anexperimental campaign has been carried out to identify the nature ofloads (global and local) and potential deformations. Numericalinvestigations are proposed to complete our understanding ofphenomena observed during experimental tests.An academic experimental campaign have been led in the flume of Ecole Centrale de Marseille. A baloon was drop from 0.5 to 2.4m with various inner pressure. Deformations and kinematics are measured and analysed before comparisons with a semi-analytical modelisation. Some results with the in-house software DITCHER are also proposed in this thesys.

Energy is central to economic development. It has been established that there is a clear correlation between energy consumption and living standards. Nigeria is a country of very industrious and enterprising people. However, due to non availability of adequate energy in the country, especially in the remote, off grid locations, the entrepreneurial inclination of the average Nigerian living in these locations has been largely stunted. Over the years, successive governments in the country, in realisation of the pivotal role of energy in national development, have explored various options to improve energy supply and availability, but the situation has not experienced any remarkable improvement. This has forced many businesses and households to resort to self provision through generators, often at exorbitant costs. This research work addresses the challenge of energy in remote, off grid locations by appraising the techno economic potential of renewable energy, using Obudu Ranch as a case study. This ranch is the foremost tourism resort in Nigeria, and has played host to a number of international events over the years. Presently, electricity is being generated through the use of diesel powered generating sets. The adjoining communities are currently without electricity, although a few of the residents have acquired generators for self provision, mostly for their domestic use. Aside the high cost associated with this, the discharge of noxious contaminants into the atmosphere is undesirable. The research entailed a working collaboration with some notable Non Governmental Organisations (NGOs) that have done extensive ground work in the area for access to some secondary data, as well as a number of corporate and governmental agencies that are relevant to the study. Further, the ranch was visited to establish hands-on, the existing renewable energy sources. A trade-off of these sources was carried out with reference to a number of relevant evaluation parameters to identify the most suited option for addressing the energy challenge. A comparative analysis of this selected source was then made to establish its techno economic potential against the existing source of power generation- diesel powered generating sets, which currently costs R1.5 million annually in running expenses. The findings from this research have established that a Renewable Energy source (mini hydro) is a more cost effective option than the diesel powered gen set, providing a 43% reduction in cost of energy generation, and a 42% reduction in the life cycle cost over the five year of analysis, compared to the status quo. In addition, it is also more environmentally friendly. Conclusively, the findings and recommendations of this research effort, if well implemented, will be beneficial to the ranch, the adjoining communities and other relevant stakeholders. / Thesis (M.Ing. (Development and Management Engineering))--North-West University, Potchefstroom Campus, 2009.

Energy is central to economic development. It has been established that there is a clear correlation between energy consumption and living standards. Nigeria is a country of very industrious and enterprising people. However, due to non availability of adequate energy in the country, especially in the remote, off grid locations, the entrepreneurial inclination of the average Nigerian living in these locations has been largely stunted. Over the years, successive governments in the country, in realisation of the pivotal role of energy in national development, have explored various options to improve energy supply and availability, but the situation has not experienced any remarkable improvement. This has forced many businesses and households to resort to self provision through generators, often at exorbitant costs. This research work addresses the challenge of energy in remote, off grid locations by appraising the techno economic potential of renewable energy, using Obudu Ranch as a case study. This ranch is the foremost tourism resort in Nigeria, and has played host to a number of international events over the years. Presently, electricity is being generated through the use of diesel powered generating sets. The adjoining communities are currently without electricity, although a few of the residents have acquired generators for self provision, mostly for their domestic use. Aside the high cost associated with this, the discharge of noxious contaminants into the atmosphere is undesirable. The research entailed a working collaboration with some notable Non Governmental Organisations (NGOs) that have done extensive ground work in the area for access to some secondary data, as well as a number of corporate and governmental agencies that are relevant to the study. Further, the ranch was visited to establish hands-on, the existing renewable energy sources. A trade-off of these sources was carried out with reference to a number of relevant evaluation parameters to identify the most suited option for addressing the energy challenge. A comparative analysis of this selected source was then made to establish its techno economic potential against the existing source of power generation- diesel powered generating sets, which currently costs R1.5 million annually in running expenses. The findings from this research have established that a Renewable Energy source (mini hydro) is a more cost effective option than the diesel powered gen set, providing a 43% reduction in cost of energy generation, and a 42% reduction in the life cycle cost over the five year of analysis, compared to the status quo. In addition, it is also more environmentally friendly. Conclusively, the findings and recommendations of this research effort, if well implemented, will be beneficial to the ranch, the adjoining communities and other relevant stakeholders. / Thesis (M.Ing. (Development and Management Engineering))--North-West University, Potchefstroom Campus, 2009.

Les tunnels profonds sont souvent construits dans les roches sédimentaires et métamorphiques stratifiées qui présentent habituellement des propriétés anisotropes en raison de leur structure et des propriétés des constituants. Le présent travail vise à étudier les tunnels profonds dans un massif rocheux anisotrope élastique en portant une attention particulière sur les effets des couplages hydromécaniques par des approches analytiques et numériques. Une solution analytique pour un tunnel creusé dans un massif rocheux anisotrope saturé est développée en tenant compte du couplage hydro-mécanique dans le régime permanent. Cette solution analytique est utilisée pour réaliser une série d’études paramétriques afin d'évaluer les effets des différents paramètres du matériau anisotrope sur le comportement du tunnel. Dans un deuxième temps la solution analytique est élargie pour décrire le comportement du tunnel pendant la phase transitoire hydraulique. Afin de compléter ces études analytiques qui prennent en compte seulement un couplage unilatéral (dans le sens que seul le comportement hydraulique influence le comportement mécanique et pas l’inverse) de l’analyse numérique avec un couplage complet, ont été réalisés. Une application de la solution analytique sur la méthode de convergence-confinement est aussi bien abordée qui peut prendre en compte l'influence du front de taille du tunnel sur le travail du soutènement ainsi que sur le massif. La solution obtenue peut servir comme un outil de dimensionnement rapide des tunnels en milieux poreux en le combinant avec des approches de dimensionnement comme celle de convergence confinement. / Deep tunnels are often built in the sedimentary and metamorphic foliated rocks which exhibits usually the anisotropic properties due to the presence of the discontinuity. The analysis of rock and liner stresses due to tunnel construction with the assumption of homogeneous and isotropic rock would be incorrect. The present thesis aims to deal with the deep tunnel in anisotropic rock with a particular emphasis on the effects of hydraulic phenomenon on the mechanical responses or reciprocal effects of hydraulic and mechanical phenomena by combining analytical and numerical approach. On that point of view, a closed-formed solution for tunnel excavated in saturated anisotropic ground is developed taking into account the hydromechanical coupling in steady-state. Based on the analytical solution obtained, parametric studies are conducted to evaluate the effects of different parameters of the anisotropic material on the tunnel behavior. The thesis considers also to extend the analytical solution with a time-dependent behavior which takes into account the impact of the pore pressure distribution on mechanical response over time, i.e., one way coupling modeling. In addition, some numerical analysis based on fully-coupled modeling, i.e., two ways coupling, are conducted which are considered as the complete solution for the analytical solution. An application of the closed-form solution on convergence-confinement method is as well referred which can take into account the influence of the tunnel face on the work of the support as well as the massif. The obtained solution could be used as a quick tool to calibrate tunnels in porous media by combining with design approaches such as convergence-confinement method.

Nous présentons deux approches différentes pour décrire le couplage hydromécanique des géomatériaux. Dans une approche de type phénoménologique nous traitons le milieu poreux comme un milieu continu équivalent dont les interactions entre la phase fluide et le squelette solide constituent le couplage du mélange à l'échelle macroscopique. En caractérisant le comportement de chaque phase nous arrivons à décrire le comportement couplé du milieu couplé saturé.Nous utilisons cette approche pour modéliser des essais expérimentaux faits sur un cylindre creux pour une roche argileuse (argile de Boom). Les résultats expérimentaux montrent de façon claire que le comportement de cette roche est fortement anisotrope. Nous avons choisi de modéliser ces essais en utilisant une lois de comportement élasto-plastique pour laquelle la partie élastique est transversalement isotrope.Le problème aux conditions aux limites étudié met en évidence des déformations localisées autour du forage intérieur. Afin de décrire de façon objective le développement de ces bandes de cisaillement nous avons considéré un milieu continu local de type second gradient qui permet d'introduire une longueur interne. De ce fait nous avons pu étudier le problème d'unicité en montrant qu'un changement de la discrétisation temporelle du problème aux limites peut conduire à des solutions différentes.Dans la deuxième approche étudiée nous caractérisons la microstructure du matériau avec des grains et un réseau de canaux pour la phase fluide. À l'aide d'un processus numérique d'homogénéisation nous arrivons à calculer numériquement la contrainte du mélange et le flux massique. Cette méthode d'homogénéisation numérique a été implémentée dans un code aux éléments finis afin d'obtenir des résultats macro. Une validation de l'implentation est proposée pour des calculs en mecanique pure et en hydromécanique. / We present two different approaches to describe the hydromechanical behaviour of geomaterials. In the first approach the porous media is studied through an equivalent continuum media where the interaction between the fluide and solid phases caracterize the coupling behaviour at the macroscale.We take into account this approach to model experimental tests performed over a hollow cylinder sample of clay rock (Boom Clay), considered for nuclear waste storage. The experimental results clearly show that the mechanical behaviour of the material is strongly anisotropic. For this reason we chose an elasto-plastic model based on Drucker-Prager criterion where the elastic part is characterized by cross anisotropy.The numerical results of boundary value problem clearly show localised strains around the inner hollow section. In order to regularize the numerical problem we consider a second gradient local continuum media with an enriched kinematic where an internal lenght can be introduced making the results mesh independent. The uniqueness study is carried out showing that changing the temporal discretization of the problem leads to different solutions.In the second approach we study the hydromechanical behaviour of a porous media that it is characterised by the microstructure of the material. The microstructure taken into account is composed by elastic grains, cohesives interfaces and a network of fluid channels. Using a periodic media a numerical homogenization (square finite element method) is considered to compute mass flux, stress and density of the mixture. In this way a pure numerical constitutive law is built from the microstructure of the media. This method has been implemented into a finite element code (Lagamine, Université de Liège) to obtain results at the macroscale. A validation of this implementation is performed for a pure mechanical boundary value problem and for a hydromechanical one.

Le mouvement de Séchilienne est situé à une trentaine de kilomètre au sud-est de l'agglomération grenobloise sur la bordure Sud du Massif de Belledonne. Il affecte le versant rocheux en rive droite de la vallée alpine de la Romanche façonné par les nombreux cycles de glaciation et déglaciation et qui lors de la dernière phase glaciaire a été complètement recouvert par le glacier. Le raidissement des pentes et l'approfondissement du pied du versant associés au phénomène de décompression affectant le versant et correspondant au relâchement des contraintes liées à la disparition progressive du glacier seraient donc les facteurs prédominants dans l'initiation et l'évolution du mouvement de Séchilienne. Les récentes contraintes temporelles obtenues sur l'escarpement sommital ont cependant montré une initiation différée du mouvement par rapport à la déglaciation et suggèrent une initiation lors de l'Optimum climatique de l'Holocène suite à une modification brusque du climat. L'initiation et l'évolution de grands mouvements de versants rocheux en domaine montagneux peuvent être le résultat de la contribution de différents facteurs géologiques (la fracturation, la lithologie ou le régime tectonique), géomorphologiques (relief, pente) mais également de facteurs de forçage externe tels que les séismes, le pergélisol, les variations hydrologiques, le retrait du glacier ou encore les effets climatiques. Ces facteurs étant tous potentiellement actifs sur le mouvement de Séchilienne, l'identification des facteurs contrôlant l'initiation et l'évolution du mouvement et pouvant expliquer le délai observé entre la déglaciation et l'initiation du mouvement est donc délicate.Nous avons donc pour cela mené une étude multidisciplinaire permettant à partir (1) des données géophysiques et d'essais en laboratoire de déterminer en profondeur les limites du mouvement et d'évaluer le volume déstabilisé ainsi que son état de déstructuration, (2) d'une analyse de la structure détaillée du versant en surface et en profondeur de préciser l'origine des fractures observées mais également leur rôle de contrôle sur le mouvement (3) d'une analyse géomorphologique des objets présents à la surface du versant, d'identifier de discrétiser et de sélectionner les structures glaciaires et gravitaires sur lesquels nous avons par la suite appliqué la méthode de datation permettant ainsi d'obtenir le modèle cinématique du retrait du glacier de la Romanche le long du versant de Séchilienne et de l'initiation et de l'évolution du mouvement (4) d'une modélisation thermique du versant basée sur la courbe d'évolution des températures reconstruite depuis la dernière phase glaciaire quaternaire de caractériser la géométrie et la persistance du pergélisol depuis la dernière phase glaciaire jusqu'à la période actuelle.La confrontation de l'ensemble de ces résultats et leur replacement dans le contexte climatique régnant depuis la dernière phase glaciaire il ya 21000ans jusqu'à l'actuelle permet ainsi de proposer un modèle conceptuel de l'initiation et de l'évolution du mouvement de Séchilienne et du rôle des différents facteurs étudiés. Le mouvement se serait donc initié suite à la décompression du versant en réponse au retrait du glacier. La mise en place et le développement du pergélisol au sein du versant aurait cependant limité la propagation du mouvement en renforçant mécaniquement le versant. Son action ne c'est cependant pas limitée à renforcer le versant mais aurait également sur une échelle de temps plus importante favoriser l'endommagement et l'affaiblissement mécanique de la zone affecté imposant ainsi au mouvement sa géométrie actuelle. / The Sechilienne landslide is located thirty kilometers at the Southeast of Grenoble on the Southern edge of the massif of Belledonne. It affects the rock slope on the right bank of the Romanche alpine valley shaped by many cycles of glaciation and deglaciation and which was completely covered by the glacier during the last glacial phase. The oversteepening of the slopes and the deepening of the foot of the slope associated with the phenomenon of debutressing affecting the slope and corresponding to the lateral stress release resulting from ice melting would be so predominant factors in the initiation and development of the Sechilienne landslide. However the recent temporal constraints obtained on the scarp showed a delayed motion relative to the initiation and deglaciation suggest initiation during the Holocene climatic Optimum after a sudden climate change. The initiation and evolution of large movements of rock slopes in mountainous area may be the result of the contribution of different geological factors (fracturing , lithology or tectonic regime) , geomorphology (relief, slope) but also external forcing factors such as earthquakes , permafrost , hydrological changes , the glacier disappearance or climatic effects. These factors are all potentially active in the Sechilienne landslide, so the identification of factors controlling the initiation and development of the movement and which can explain the delay observed between deglaciation and the initiation of movement is tricky.Therefore, we conducted a multidisciplinary study from (1) geophysical data and laboratory tests to determine the depth limits of the mass movement and to assess the volume and the breakdown state, (2) an analysis of the detailed structures of the slope at the surface and in depth to identify the origin of the observed fractures and also their role in the control of the movement (3) a geomorphological analysis of objects on the surface of the slope to discretize , to identify and to select the glacial and gravitational structures on which then we applied the datation method to obtain the kinematic model of the Romanche glacier retreat along the slope of Séchilienne and of the initiation and the evolution of the movement (4) a thermal modeling of the slope based on the temperature curve proposed since the last quaternary glacial stage to characterize the geometry and the persistence of permafrost since the last glacial phase to the current period.The comparison of all these results and their replacement in the climatic conditions prevailing since the last glacial stage 21000 years ago to the present, allows to propose a conceptual model for the initiation and the evolution of Séchilienne landslide and the role of the various studied factors. So, the landslide would be initiated by slope debutressing in response to the glacier retreat. However the establishment and the development of the permafrost within the slope would limit the propagation of the movement by reinforcing mechanically the slope.