L'étude du comportement du béton en température est un problème majeur qui vise in fine à évaluer le niveau de sécurité des structures sous des sollicitations thermiques sévères, lors d'incendies par exemple. Pour cela, de nombreux modèles sont développés dans un cadre de couplage thermo-hydro-mécanique (THM), pour prendre en compte les différents processus physico-chimiques et mécaniques mobilisés par ces sollicitations et conditionnant la stabilité du matériau béton.Cependant, ces modélisations ne prennent souvent pas en compte explicitement la nature hétérogène du matériau. En effet, le béton est un matériau constitué d'inclusions noyées dans une matrice cimentaire pouvant être schématisée comme un milieu poreux ouvert partiellement saturé en eau. Les inclusions sont caractérisées par leurs natures minéralogiques, leurs morphologies et leurs tailles classées en fuseaux granulaires. Cette hétérogénéité introduit une complexité supplémentaire : la nécessité de prendre en compte la microstructure pour quantifier l'effet de l'incompatibilité (thermique, hydrique et mécanique) inclusion-matrice sur le comportement THM du béton. Ce travail constitue une première avancée dans ce sens. A ce titre, une modélisation élément fini multi-échelle tridimensionnelle (3D) est développée permettant d'affecter des comportements spécifiques à la matrice et aux inclusions. Pour la matrice, siège des transports de masse dans son réseau poreux connecté, une approche thermo-hydrique à trois fluides (eau, vapeur et air sec) est adoptée et est couplée à une loi de comportement poro-mécanique endommageable. Les inclusions (granulats) sont considérées hydriquement inertes, une approche thermo-mécanique avec endommagement est alors retenue.Une analyse, par simulations numériques, des effets de la nature minéralogique des granulats (calcaires ou silico-calcaires) de leurs distributions et de leurs morphologies a été menée sur des configurations 2D et 3D. Les effets étudiés ont notamment concerné l'influence de ces paramètres sur les fluctuations locales des champs de température, de pression de gaz et d'endommagement au regard de la dispersion des mesures expérimentales. L'analyse est limitée à l'échelle mésoscopique, celle où les granulats de taille caractéristique supérieure à 5 mm sont considérés, le reste (stable et pâte de ciment) étant une matrice homogène. Enfin, cette analyse a mis en évidence le besoin de mettre en œuvre une approche expérimentale cohérente avec une analyse mutli-échelle, à la fois pour la caractérisation des propriétés (thermiques, hydriques et mécaniques) de chaque constituant et pour l'étude des évolutions des champs lors des changements d'échelles. Un protocole expérimental a été définit et des premiers résultats d'essais sont présentés et analysés au travers de résultats obtenus dans la partie modélisation / The investigation of the behavior of heated concrete is a major research topic which concerns the assessment of safety level of structures when exposed to high temperatures, for instance during a fire. For this purpose, several modeling approaches were developed within thermo-hydro-mechanical (THM) frameworks in order to take into account the involved physic-chemical and mechanical processes that affect stability of heated concrete. However, existing models often do note account explicitly for the heterogeneity of the material : concrete is composite material that may be schematized as an assembly of inclusions (aggregates) embedded in a cementitious matrix (cement paste). This latter may be described as a partially saturated open porous medium. The aggregates are characterized by their mineralogical nature together with their morphology and size distribution. The material heterogeneity bring an additional complexity : the need to take into account the microstructure in order to quantify the effect of matrix-inclusion thermal, hygral and mechanical incompatibilities on the THM behavior of concrete. This work is a first step in this direction. For this purpose, a three-dimensional (3D) multi-scale finite element model is developed. It allows affecting specific behaviors to matrix and inclusions. For the former, where mass transports occur within the connected porous network, a three-fluids approach (liquid water, vapor and dry air) is adopted and is coupled to a poro-mechanical damage based approach. For inclusions (aggregates) no hygral component arises a pure thermo-mechanical model is considered. The developed model is then used to investigate, either by 2D or 3D numerical simulations, effects of mineralogical nature, morphology and distribution of aggregates. Studied effects have mainly concerned the influence of these parameters on local fluctuations of simulated temperature, gas pressure and damage fields with regard to experimentally observed dispersion. The analysis is here limited to the mesoscale, at which only aggregates with a characteristic size above 5 mm are meshed while the remaining inclusions together with the cement paste are considered to be a homogeneous matrix. Finally, the numerical analysis carries out the need to perform an experimental campaign that is consistent with a multi-scale approach of the THM behavior of concrete : an experimental campaign that allows to identify thermal, hygral and mechanical properties of each concrete constituent and that permit to assess evolution of fields during upscalling. An experimental protocol is then elaborated for this purpose and some obtained results are presented and analyzed with regard to results obtained in the modeling part

La réduction du volume des glaciers et la fusion printanière plus précoce de la neige causée par le réchauffement climatique provoquent des variations du cycle hydrologique à la fois pour les têtes de bassin versant, mais aussi pour les zones situées plus à l’aval. Afin de prédire correctement l’amplitude des changements possibles futurs et d’envisager une gestion adaptée, une bonne connaissance de l’interaction entre les glaciers, le climat et les écoulements hydriques est nécessaire. L’objectif de cette étude est d’évaluer l’effet de la variabilité climatique sur le fonctionnement hydro-glaciologique et ses conséquences sur la disponibilité de l’eau du bassin versant de l’Arve (Alpes françaises) depuis 1960. Ce bassin s’étend sur une surface de 1958 km2 et est composé de cinq bassins versants emboités (Arveyron d’Argentière, Arveyron de la Mer de Glace, Arve au Pont des Favrands, à Sallanches and au Bout du Monde), tous influencés par la fusion glaciaire et nivale mais dans différentes proportions étant donnée la large gamme d’extension de couverture glaciaire s’étalant de 5 á 53%. Ce travail est basé sur des longs jeux de données glaciologiques, météorologiques, hydrologiques et de couverture de neige qui sont issues soit de mesures ponctuelles dans l’espace soit de données obtenues par télédétection.L’analyse des tendances a été réalisée sur des données hydrologiques et météorologiques des cinq bassins versants emboités. Pour cela, le cycle saisonnier du débit est ajusté en utilisant une fonction mathématique de type “modèle à pic asymétrique”. Les changements observés des débits ont été reliés aux variables météorologiques ainsi que à l’évolution de la couverture glaciaire. Les résultats indiquent un comportement contrasté entre les bassins versants selon les taux d’englacements, avec une tendance croissante des valeurs de débit dans les bassins versants fortement englacés (couverture de glacier >30%) et une décroissante pour les moins englacés. La sensibilité du cycle hydrologique au changement climatique futur a été évaluée. Pour le milieu du 21e siècle, on prévoit que le volume annuel écoulé serait réduit de 16% pour l’Arveyron d’Argentière et de 31% pour l’Arveyron de la Mer de Glace. Pour la période estivale, la quantification détaillée de chaque terme de l’équation du bilan hydrologique, ainsi que leurs incertitudes, sur les bassins versants de l’Arveyron d’Argentière et de l’Arveyron de la Mer de Glace-Leschaux a permis de souligner l’importance des transferts d’eau souterraine pour représenter et prédire le comportement hydro-glaciologique d’un bassin versant donné. Deux model d’écoulement distribués de type degré-jour couples à un modèle de routage hydrologique à réservoir linéaire ont était utilisé sur le bassin versant de l’Arveyron d’Argentière sur la période 1960–2009. La calibration est effectuée autant sur la base des données de débit qu’avec une approche multicritère avec les données de débit, de couverture neigeuse et du bilan de masse annuel, à pas de temps journalier. Les résultats montrent l’aptitude d’utiliser un modèle classique degré-jour pour simuler le comportement hydro-glaciologique et la production d’eau sous-glaciaire d’un bassin versant fortement glaciaire. Pour la période 1960–2004, une valeur de Kling Gupta Efficiency de 0.85 entre le débit simulé et observe à était obtenu. La calibration multicritère semble réduire les incertitudes des simulations. / Glacier recession and the anticipation of spring snow melt driven by a warming climate could lead to changes in the hydrological cycle affecting not only the headwater catchments but also the areas downstream. In order to correctly predict the magnitude of future possible changes and to consider appropriate strategies of water management, a good understanding of the interaction between glaciers, climate and hydrology is needed. The aim of this study is to assess the effect of climate variability on the hydro-glaciological behaviour and its consequence on water availability in the Arve River catchment (French Alps) since 1960. It covers 1958 km2 and is composed by five nested catchments (Arveyron d’Argentière, Arveyron de la Mer de Glace, Arve at Pont des Favrands, Arve at Sallanches and Arve at Bout du Monde), all influenced by glacier and snow melt but characterized by various percentages of glacier cover ranging from 5 to 53%. This research is based on a long dataset of in situ or remote sensing glaciological, meteorological, hydrological and snow cover area data.Trend analyses are performed on the hydrological and meteorological data at all the considered sites. The seasonal cycle of each catchment is fitted using a mathematical function, namely the asymmetric peak model, and changes in the discharge are related to observed changes in the meteorological variables and the glaciers’ evolution. Results point out a contrasting behaviour among the catchments characterized by different glacier covers, showing an increasing trend on the discharge values in highly glacierized catchments (with a glacier cover >30%) and a decrease in the low glacierized ones. The sensitivity of the seasonal cycle to the future climate is evaluated. In the mid-21st century the annual runoff would be reduced by 16% for Arveyron d’Argentière and 31% for Arveyron de la Mer de Glace. Over the summer season, a detailed quantification of each term of the hydrological balance equation, as well as their uncertainties, on the Argentière and Mer de Glace-Leschaux drainage basins allows to underline the importance of considering the groundwater transfers to represent and predict the hydro-glaciological behaviour of a considered catchment. Two different distributed temperature index melt models coupled with a linear reservoir discharge model are used on the Arveron d’Argentière catchment over the 1960–2009 period. The calibration is carried out against discharge only and with a multi- criteria approach considering the discharge, the snow cover area and the glacier-wide annual mass balance values at daily time step. Results demonstrate the suitability of the use of a classical degree day model in simulating the hydro-glaciological behaviour and the subglacial water production of a highly glacierized catchment. A KGE of 0.85 is obtained between the observed and simulate discharge values over the 1960–2004 period. The use of a multi-criteria approach seems to reduce the simulation uncertainties.

In a hydro-dominated system, such as New Zealand, the continual improvement and development of effective optimization and simulation software to inform decision making is necessary for effective resource management. Stochastic Constructive Dual Dynamic Programming (SCDDP) is a technique which has been effectively applied to the New Zealand system for optimization and simulation. This variant of Dynamic Programming (DP) allows optimization to occur in the dual space reducing the computational complexity and allows solutions from a single run to be formed as price signal surfaces and trajectories. However, any application of this method suffers from issues with computational tractability for higher reservoir numbers. Furthermore, New Zealand specific applications currently provide limited information on the system as they all use the same two-reservoir approximation of the New Zealand system. This limitation is of increasing importance with the decentralization of the New Zealand electricity sector. In this thesis we develop this theory with respect to two key goals: • To advance the theory surrounding SCDDP to be generalizable to higher reservoir numbers through the application of the point-wise algorithm explored in R. A. Read, Dye, S. & Read, E.G. (2012) to the stochastic case. • To develop at least two new and distinct two-reservoir SCDDP representations of the New Zealand system to provide a theoretical basis for greater flexibility in simulation and optimization of hydro-thermal scheduling in the New Zealand context.

Der weltweit steigende Energiebedarf stellt die Menschheit vor immer größere Herausforderungen. Im Angesicht des Klimawandels und der begrenzten Verfügbarkeit von fossilen Energieträgern liegt eine besondere Verantwortung bei der Entwicklung von erneuerbaren Energieressourcen. Dabei spielt die Geothermie eine besondere Rolle, da sie zur Deckung der Grundlast geeignet ist. Zu der Bewertung von geothermischen Potentialen leisten die Geowissenschaften einen großen Beitrag. Das geothermische Potential eines Standorts hängt vor allem von der Art und Ausbreitung von Störungszonen ab, da sie wesentlich den Grundwasserfluss kontrollieren. In dieser Arbeit werden verschiedene geowissenschaftliche Methoden miteinander kombiniert, um die Wirksamkeit von Störungszonen zu bewerten. Strukturgeologische Kartierungen und hydrogeologische Felduntersuchungen erklären die Ausbreitung und hydraulische Funktion von Störungszonen im Untersuchungsgebiet. Geohydrochemische Untersuchungen geben Auskunft über Fluid- und Gesteinszusammensetzung und deren Wechselwirkungen. Numerische Simulationen des Gebietes zeigen, dass sowohl Störungszonen als auch Fluideigenschaften wichtig für die Verteilung von Druck und Temperatur im Reservoir sind. Das Untersuchungsgebiet ist das Hochenthalpy-Geothermiefeld Lahendong in Sulawesi (Indonesien). Ein Kraftwerk produziert hier 80 MW Elektrizität. Die ersten Untersuchungen im Gebiet starteten in den 1970er Jahren. Jedoch sind Genese und Verteilung der thermalen Wässer noch nicht komplett verstanden. Das aktuelle konzeptionelle Modell zeigt eine Aufteilung in mehrere Reservoirbereiche. Die verschiedenen Bereiche sind durch horizontal impermeable Störungszonen voneinander getrennt. Den nördlichen Bereich kennzeichnen saure Wässer und den südlichen Teil pH-neutrale Wässer. Auch die Produktivität variiert stark zwischen den beiden Reservoirbereichen. In vertikaler Richtung sind Störungszonen jedoch durchlässig, was dazu führt, dass heiße Quellen entlang von Störungen oder deren Kreuzungspunkten auftreten. Die Reservoirgesteine in Lahendong sind basaltische Andesite, Tuffe und vulkanische Brekkzien. Die Permeabilitätsverteilung der Störungszonen wird durch die Ausbreitung von Rissen kontrolliert. Risse sind vor allem in der Bruchzone der Störung verbreitet, was zu einer hydraulischen Durchlässigkeit parallel zur Störung führt. Die Unterteilung des Reservoirs in Lahendong wird durch diese Rissverteilung bestimmt. Der nördliche saure Teil des Reservoirs ist durch höhere Produktivitätsraten, Gasaustritte an der Oberfläche und stark alterierte und geklüftete Gesteine im Untergrund charakterisiert. Der südliche Abschnitt ist heißer und hat weniger stark alterierte Gesteine. Die beobachteten Reservoireigenschaften wurden von hydrochemischen und hydraulischen Modellierungen bestätigt. Der Grundwasserfluss mit Neubildung und Austritten aus dem Reservoir ist im Gelände und im Modell vor allem durch Störungszonen kontrolliert. Jedoch ist der Grundwasserfluss auch durch den Aggregatzustand des Wassers beeinflusst. Für die Modellierung von 2-Phasen-Ausbreitung müssen die Permeabilitäten angepasst werden, um gleiche Temperatur- und Druckbedingungen zu modellieren. Der Hauptbeitrag dieser Arbeit liegt in der Demonstration, dass eine systematische strukturgeologische Analyse für das Verständnis von Grundwasserfluss in geothermischen Reservoiren unentbehrlich ist. Es wurde bestätigt, dass die Kombination von tektonischen, hydrogeologischen und geohydrochemischen Informationen den wichtigsten Beitrag für das Verständnis von Grundwasserströmungen leistet. Die Grundwasserströmung ist der wichtigste Faktor für die Wahl des richtigen Standortes für Produktion und Injektion in geothermischen Feldern. Eine detaillierte Analyse gewährleistet eine nachhaltige Nutzung des Feldes und verringert Risiken, wie schwach produzierende Bohrungen oder die Produktion von stark korrosiven Wässern. Auf dieser Grundlage kann eine Felderweiterung geplant werden, wie es in Lahendong angedacht ist.

L'étude de sûreté d'un stockage souterrain de déchets radioactifs passe, du point de vue de la géomécanique, par la caractérisation de la zone endommagée induite par le creusement du puits d'accès aux cavités de stockage. L'objet de cette thèse est donc de simuler la réponse hydro-mécanique des argilites de l'Est au creusement du puits d'accès au laboratoire souterrain de Meuse/Haute-Marne, en vue d'une comparaison avec les mesures qui seront acquise in situ en 2005 (expérimentation REP).<br/> Nous avons d'abord analysé l'ensemble des essais mécaniques qui ont été réalisés entre 1995 et 2001. Nous observons que chaque série d'essais a sa propre cohérence, que ce soit en termes de paramètres élastiques, de résistance à la rupture ou de capacité au fluage. Mais il existe de fortes différences entre les séries. Ces différences, qui sont dues au protocoles expérimentaux et non au matériau, nous ont permis de mettre en évidence trois résultats importants : les phases de resaturation imposées par plusieurs laboratoires ont dégradé les propriétés mécaniques des argilites de l'Est, l'existence d'une contrainte effective au sens de Biot n'est pas évidente pour ces argilites, et leur résistance à la rupture augmente nettement avec le séchage.<br/>D'après ces essais, nous avons distingué trois mécanismes de déformation irréversible : compaction, pré-rupture, rupture. Nous les avons décrits chacun par une loi de comportement élasto-plastique simple, fondée sur le critère de Mohr-Coulomb ou de Drücker-Prager, avec un écrouissage linéaire. On obtient alors une loi de comportement complète pour les argilites de l'Est en considérant les trois mécanismes simultanément (plasticité multi-critère).<br/>Enfin, les différentes modélisations que nous avons réalisées permettent de tirer les conclusions suivantes :la prise en compte des déformations irréversibles avant-rupture (compaction et pré-rupture) a pour effet de réduire considérablement l'extension de la zone en rupture qui apparaît autour du puits (par rapport au cas où le comportement est élastique jusqu'à la rupture) ; la diffusion des pressions interstitielles modifie peu la réponse mécanique du massif (en particulier, les déformations plastiques créées lors du déchargement mécanique n'évoluent pas au cours du drainage du massif) ; l'anisotropie des contraintes naturelles horizontales, dont par ailleurs nous confirmons l'existence, influence fortement la zone endommagée autour du puits (polarisation en azimut) ; enfin nous montrons qu'il n'est pas indispensable de simuler le détail des opérations de creusement du puits.

Ce mémoire est consacré à l'étude de la rétention d'eau et de la consolidation des sols non saturés dans un cadre thermo-hydro-mécanique. Il est organisé en deux parties qui concernent respectivement "l'effet de la déformation et de la température sur le comportement hydrique" et "l'effet de la température et de la succion sur le comportement mécanique" de sols non saturés. Dans la première partie, on montre la pertinence du concept de surface caractéristique d'un sol par opposition à la courbe de rétention dont la portée s'avère limitée dans le cas de milieux très déformables. Ce concept est illustré expérimentalement dans le cas d'un sable limoneux argileux. La modélisation mathématique de la surface caractéristique s'appuie sur une importante campagne expérimentale comprenant la mesure d'environ 240 triplets (indice des vides, teneur en eau, succion). Par ailleurs, un modèle de comportement thermo-hydrique est proposé afin de déterminer la surface caractéristique et la courbe de rétention à une température donnée. Ce modèle est validé dans le cas de deux matériaux (une céramique à base de terre cuite et un sable limoneux argileux) testés par l'auteur ainsi que dans le cas d'autres matériaux de la littérature. Enfin, on présente une application à la détermination de la perméabilité des sols non saturés prenant en compte la déformation et la température. Dans la deuxième partie, l'effet de la température et de la succion sur le comportement mécanique est examiné à partir d'essais de consolidation du limon de Sion. Ces essais sont réalisés à différents niveaux de température et de succion. Pour chacun de ces essais, les indices de gonflement et de compression ainsi que la contrainte de préconsolidation sont mesurés. On détermine ainsi l'influence de la température et de la succion sur ces paramètres essentiels du comportement mécanique. Enfin, on propose une fonction d'approximation modélisant l'évolution de la contrainte de préconsolidation en fonction de la température et de la succion.

The present work deals with the behavior of ancient railway sub-structure in France. A statistical study was firstly undertaken on problems occurred in the whole ancient French railway network. The analysis evidenced the particular importance of sub-grade quality for the performance of the sub-structure and the track geometry. Afterwards, an ancient railway line in the West of France was investigated. The analysis showed that the degradation speed of this line was correlated with different parameters such as the nature of sub-grades and the thickness of different layers. An increase trend of degradation speed with the increase in interlayer thickness was identified. The interlayer has a positive impact since it reduces the train-induced stress applied to the sub-grade. The hydro-mechanical behavior of interlayer soil under different conditions (water content, fines content, stress, number of cycles) was investigated. A set of triaxial tests and infiltration tests were performed for this purpose. By analyzing the shear strength properties, the permanent axial strain and the resilient modulus of interlayer soil, we found that the water content and the fines content must be considered together. Adding more fines into the interlayer presents a positive impact under unsaturated conditions thanks to the suction effect, but a negative impact under saturated conditions. The infiltration column tests with drying/wetting cycles showed that the hydraulic conductivity of interlayer soil is governed by fines fraction but did not change significantly with fines content. In order to study the mechanism of interlayer creation and mud pumping, a physical model of 550 mm inner diameter was developed. Soil samples representing the ancient French railway substructure with a ballast layer overlying an artificial silt layer (mixture of crushed sand and kaolin were tested. The effects of monotonic and cyclic loadings, water content and dry unit mass of sub-soil were investigated. It was found that the pore water pressure developed in the sub-soil and the sub-soil stiffness are the key factors for the migration of fine particles or the creation of interlayer/mud pumping. Water is the necessary condition, but it is the soil compressibility that governs the phenomenon to occur

M. Tech. (Chemical Engineering, Faculty of Engineering and Technology), Vaal University of Technology / Acid mine drainage (AMD) is a serious problem associated with mining activities, and it has the potential to contaminate surface and ground water. The aim of this study was to evaluate the performance of constructed wetland and photocatalysis in treating AMD. Three identical unvegetated upflow constructed wetlands packed with natural zeolite (clinoptilolite) and coarse silica sand were made of a cylindrical plastic pipe, and the slurry photocatalyst was prepared using quartz material. A hydro-alcohol thermal method was used to prepare an anatase core-void-shell TiO2 photocatalyst. The results showed that the three unvegetated upflow constructed wetlands (CW) had relatively similar percentage removal of heavy metals despite their varying concentrations within the AMD. The removals were: Fe (86.54 - 90.4%); Cr (56.2 - 64.5%); Mg (56.2 - 67.88%); Ca (77.1 - 100%); and 100% removal was achieved for Be, Zn, Co, Ni, and Mn. The removal of sulphate was also 30%. Heavy metals concentration in CW packing material was significantly higher in the outlet of the constructed wetlands than in the inlet. The adsorption isotherms revealed that the experimental data fitted the Langmuir Isotherms better, which suggested a monolayer coverage of heavy metals on the surface of the adsorbents; thermodynamic studies showed that the nature of adsorption taking place was physical; the kinetics models showed that the adsorption was first order reaction. A higher photocatalytic reduction (62%) of Cr(VI) was obtained at pH 2, 30 mg/l Cr(VI) initial concentration, and three hours of irradiation time. It was also found that the presence of Fe(III) enhanced the reduction of Cr(VI). The core-void-shell TiO2 photocatalyst showed a better activity than the commercial P25 Degussa for the reduction of Cr(VI) to Cr(III). The kinetic studies showed that the reduction of Cr(VI) was first order reaction. Photocatalytic reduction of Cr(VI) in real AMD sample was achieved only for the Douglas North Discharge (DND) sample (68%), and the Fe(III) reduction was found to be 83%. Therefore it was concluded that the combination of constructed wetland and UV/ TiO2 photocatalysis employing anatase core-void-shell TiO2 as a photocatalyst has a potential to reduce the toxicity of Cr(VI)-laden acid mine drainage.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal.

Les sols argileux peuvent présenter des variations importantes de volume lors d’un changement hydrique, tel une période de sécheresse. Ce phénomène de retrait-gonflement des sols argileux est à l’origine de nombreux préjudices aux bâtis ce qui le situe, après les inondations, au second rang des aléas naturels français en matière de coût des dommages : 5 milliard d’euros entre 1988 et 2007. Les maisons individuelles sont les plus touchées par ce phénomène qui provoque un tassement différentiel du bâti à l’origine de sa dégradation (fissures dans les éléments de maçonnerie non armée en particulier). L’objectif de cette thèse est d’étudier la vulnérabilité des bâtis vis-à-vis de l’aléa retrait-gonflement à travers une analyse de l’interaction sol-structure. Ce travail s’est porté sur les trois grands volets suivants : 1 - le comportement hydromécanique des sols argileux, 2 - l’interaction sol-structure, 3 - l’évaluation des dommages structuraux. Le comportement hydromécanique des sols argileux a été modélisé par le concept de surface d’état et intégré dans un modèle d’interaction sol-structure inspiré du modèle de Winkler, avec des approches de plus en plus complexes : cas unidimensionnel tout d’abord, 2D ensuite en assimilant la structure à une poutre, puis en 3D en représentant la structure par une plaque. Ces modèles d’interaction sol-structure permettent de calculer la déflexion relative d’un bâti, en fonction d’une valeur de succion imposée et des propriétés mécaniques ou hydromécaniques de la structure et du sol d’assise. La déflexion relative du bâti est alors comparée à des valeurs seuils afin d’évaluer le dommage.La dernière partie de cette thèse consiste en une étude de faisabilité pour le développement de courbes de vulnérabilités adaptées à la problématique du retrait gonflement. Ces courbes représentent la moyenne des dommages d’un type de bâti conçu sur un type de sol gonflant pour une valeur de succion imposée. Ces courbes ont été développées à partir d’une typologie des sols gonflants et d’une typologie des bâtis issue des bases de données des maisons sinistrées. La moyenne des dommages est calculée par la méthode de Monte-Carlo, en prenant en compte la variabilité des paramètres du bâti. / Clayey soils can present large volumetric deformations in response of water content change. This phenomenon of shrinkage and swelling of clayey soils is recognized as a costly natural hazard throughout the world. In France, this cost is reported between 1988 and 2007 to 5 million Euros, ranking in second class of the French natural hazards in terms of cost of damage after the floods. The individual masonry buildings with shallow foundations are the most affected by this phenomenon that causes a differential settlement of the building leading to cracks in facades and structural elements, especially in unreinforced masonry elements.The objective of this thesis was to study the vulnerability of the building, against the shrink-swell hazard through an analysis of soil-structure interaction. This work was focused on three major points: - Hydro-mechanical behavior of clay soils - Soil-structure interaction - Evaluation of structural damage.The hydro-mechanical behavior of clay soils was considered by the concept of state surface and integrated into a soil-structure interaction model, based on the Winkler model, with increasingly complex approaches: 1D, 2D and 3D by modeling the building behavior respectively by a spring system, the beam element and the plate element. These models of soil-structure interaction are able to calculate the relative deflection of the building, according to the values of imposed suction, mechanical properties of the structure and hydro-mechanical properties of the soil. The relative deflection of the building is then compared to threshold values of the classical damage categories to assess the building damage.Finally a feasibility study was conducted on the development of vulnerability curves adapted to the problem of shrinkage – swelling of clayey soils. These curves represent the average of structural damage versus suction, for each type of building. These curves were developed based on a classification of expansive soils and a building typology outcome from databases of affected buildings. The average damage is calculated by the Monte-Carlo method, taking into account the variability of the building parameters.