Dans le cadre des études de faisabilité du stockage des déchets radioactifs de haute activité et à vie longue en formation géologique profonde, nous sommes amenés à effectuer des modélisations des phénomènes de couplage THM à court et à long terme engendrés par l’excavation et le stockage. Grâce à sa faible perméabilité qui impose un transport par diffusion et sa capacité de rétention, les argilites du Callovo-Oxfordien ont été choisies comme potentielles barrières géologiques pour le stockage. Au vu de la complexité des phénomènes mises en jeu, la thèse se limitera à l’étude de son comportement mécanique à long terme, volet essentiel à l’étude des couplages précités. Pour la caractérisation du comportement différé des argilites, deux mécanismes de fluage sont envisagés: la déformation viscoplastique de la matrice argileuse et la propagation sub-critique de microfissures. Dans ce travail, les deux mécanismes de fluage sont étudiés respectivement. Comme les approches micromécaniques permettent de conduire à une meilleure description du comportement en relation avec les aspects microstructuraux de la roche et d'éviter ainsi le recours à une identification des paramètres pour chaque zone géomecanique, dans ce travail, nous avons mis en oeuvre une modélisation multi-échelle construite à partir des comportements des constituants et des données microstructurales de l'argilite. / Claystones are extensively investigated in France in the framework of feasibility studies for geological repository of nuclear waste. Under the coordination of the French Agency for radioactive waste management (ANDRA), laboratory tests and numerical modeling have been performed for the characterization of thermo-hydromechanical properties of claystones. Due to the very long life of radioactive waste, the characterization of long term hydromechanical behaviors of claystones is a crucial issue. For this purpose, laboratory creep tests associated with microstructural analysis have been performed on claystones with different degrees of saturation. It is found that two main mechanisms are at the origin of macroscopic time dependent deformation. The first one is related to viscoplastic flow of clay matrix and the second one is due to subcritical propagation of microcracks in clay matrix. In this thesis, two creep mechanisms are studied respectively. As micromechanical approaches can lead to a better description of behavior in relation to the microstructure of the rock, in this work, we implemented a multi-scale model based on the behaviors of each component and microstructure of claystones.

Au cours de ce travail, nous nous sommes attachés à prédire les caractéristiques rhéologiques de suspensions de particules dans des situations non newtoniennes. Nous avons tout d'abord étudié expérimentalement le comportement rhéologique de suspensions bidisperses de particules non colloïdales et non browniennes immergées dans un fluide à seuil dans les situations où les interactions entre particules sont hydrodynamiques (mécaniques). Nous avons montré que l'influence de la polydispersité des particules sur le module élastique et le seuil de contrainte de la suspension pouvait être décrit en utilisant un modèle d'empilement développé pour les matériaux granulaires secs. Deux relations analytiques simples permettant d'estimer le module élastique et la contrainte seuil d'une suspension non newtonienne polydisperse en fonction des caractéristiques des constituants ont alors été proposées et validées. Nous nous sommes ensuite intéressés à la prédiction des caractéristiques globales de suspensions contenant une forte concentration de particules de taille différente dans le cadre d'une approche par changement d'échelle. Pour cela, nous avons mis en uvre une approche par motifs morphologiques permettant la prise en compte de paramètres géométriques tels que la distribution de tailles des inclusions ou les distances entre particules voisines. Les résultats de ce modèle ont été validés par comparaison avec des données expérimentales de la littérature. Nous avons ensuite montré que les modèles morphologiques permettaient également de décrire l'évolution des modules complexes d'un matériau thixotrope subissant une transition de phase liquide-solide due à la croissance et la percolation d'un réseau de particules attractives. Dans la dernière partie de cette étude, nous avons étudié expérimentalement l'influence de l'histoire de cisaillement sur la structuration et la déstructuration d'une suspension colloïdale thixotrope

Au cours de ce travail, nous nous sommes attachés à prédire les caractéristiques rhéologiques de suspensions de particules dans des situations non newtoniennes. Nous avons tout d'abord étudié expérimentalement le comportement rhéologique de suspensions bidisperses de particules non colloïdales et non browniennes immergées dans un fluide à seuil dans les situations où les interactions entre particules sont hydrodynamiques (mécaniques). Nous avons montré que l'influence de la polydispersité des particules sur le module élastique et le seuil de contrainte de la suspension pouvait être décrit en utilisant un modèle d'empilement développé pour les matériaux granulaires secs. Deux relations analytiques simples permettant d'estimer le module élastique et la contrainte seuil d'une suspension non newtonienne polydisperse en fonction des caractéristiques des constituants ont alors été proposées et validées. Nous nous sommes ensuite intéressés à la prédiction des caractéristiques globales de suspensions contenant une forte concentration de particules de taille différente dans le cadre d'une approche par changement d'échelle. Pour cela, nous avons mis en uvre une approche par motifs morphologiques permettant la prise en compte de paramètres géométriques tels que la distribution de tailles des inclusions ou les distances entre particules voisines. Les résultats de ce modèle ont été validés par comparaison avec des données expérimentales de la littérature. Nous avons ensuite montré que les modèles morphologiques permettaient également de décrire l'évolution des modules complexes d'un matériau thixotrope subissant une transition de phase liquide-solide due à la croissance et la percolation d'un réseau de particules attractives. Dans la dernière partie de cette étude, nous avons étudié expérimentalement l'influence de l'histoire de cisaillement sur la structuration et la déstructuration d'une suspension colloïdale thixotrope / In this work, we aimed to predict the rheological properties of suspensions of particles in non-Newtonian fluids. We first studied experimentally the rheological behaviour of bidisperse suspensions of non colloidal particles and non-Brownian immersed in a yield stress fluid in situations where interactions between particles are hydrodynamic (mechanical). We showed that the influence of the polydispersity on the elastic modulus and yield stress of the suspension could be described using a packing model developed for dry granular materials.Two simple closed form estimates for both the elastic modulus and the yield stress of non-Newtonian polydisperse suspension as a function of the characteristics of the constituents were proposed and validated. This study was then extended in a micromechanical analytical framework to predict the macroscopic characteristics of suspensions containing highly concentration of particles with different size. Our ne w model, called a "Morphologically Representative Pattern - Based Approach" (MRP-based approach), successfully predicts the behaviour of such materials by taking into account some geometrical parameters such as the particle size distribution or the distance between the nearest-neighbour particles. This model has been validated by a comparison with the experimental data in the literature. We then showed that the MRP-based approach also allowed describing the evolution of complex modules of a thixotropic material undergoing a liquid-solid transition due to the growth and the percolation of an attractive particle network. In the last part of this work, we studied experimentally the influence of the flow history on the structure of a thixotropic colloidal suspension

Le stockage en formation géologique profonde est l'une des solutions possibles pour la gestion à très long terme des déchets radioactifs de Haute Activité et à Vie Longue (HAVL). Son étude nécessite d'acquérir une bonne connaissance du comportement et des propriétés des formations géologiques potentielles ainsi que de leur évolution au cours du temps sous l'effet des sollicitations induites par une éventuelle installation de stockage.<br />La formation géologique hôte sera soumise à des chargements mécanique et thermique dus respectivement au creusement des ouvrages de stockage et au dégagement de chaleur des colis de déchets stockés. Ces chargements thermomécaniques vont engendrer une modification de la distribution des contraintes dans la couche hôte et les déformations des ouvrages ainsi que l'extension des zones endommagées (EDZ) pourraient provoquer des instabilités locales et globales. <br />Les objectifs de ce travail consistent en l'élaboration de techniques de calculs pour optimiser la modélisation numérique du comportement thermoélastique des ouvrages à grande échelle et en l'évaluation de la perturbation thermomécanique induite par le stockage sur la formation géologique hôte. Dans cette optique, après une présentation de l'état des connaissances sur les aspects thermomécaniques des roches liés au stockage souterrain, des modélisations numériques 2D et 3D du comportement thermoélastique d'ouvrage individuel et d'un réseau de galeries ont été réalisées par une approche discrète. Toutefois, cette approche classique est pénalisante pour l'étude du comportement global d'installations de stockage. Pour pallier à cela, une approche de modélisation numérique, basée sur la technique d'homogénéisation des milieux périodiques, a été proposée. Des formulations ainsi que des procédures numériques ont été élaborées pour déterminer le comportement thermoélastique effectif d'une structure hétérogène équivalente. Le modèle, obtenu par cette méthode, a été validé avec des méthodes d'homogénéisation existantes telles que le modèle auto-cohérent, ainsi que les bornes de Hashin-Shtrikman. La comparaison du comportement thermoélastique effectif avec le comportement thermoélastique réel de référence a montré une bonne cohérence des résultats. Pour une application au stockage profond, les coefficients thermoélastiques effectifs d'un réseau de galeries circulaires ont pu être déterminés en 2D pour différentes dimensions de l'entraxe des galeries.

On s'intéresse plus particulièrement à la modélisation du comportement d'un matériau hétérogène méso-fissuré (béton, roche,...), soumis à une sollicitation thermo-hydro-mécanique avec prise en compte du couplage géométrique. Pour conduire cette étude, on s'appuie notamment sur les approches micro-mécaniques du comportement des milieux méso-fissurés non saturés développées depuis quelques années au Laboratoire des Matériaux et des Structures du Génie Civil - Ur Navier - Université Paris Est. Le milieu fissuré non saturé traité ici est constitué d'une matrice solide homogène élastique linéaire et de fissures connectées saturées par deux fluides immiscibles : un liquide et un gaz séparés par une surface capillaire. La fissure est traditionnellement considérée comme une cavité ellipsoïdale (cas 3D) ou elliptique (cas 2D) dont le rapport d'aspect tend vers zéro. Deux morphologies typiques de matériau sont considérés dans ce travail : la situation où les fissures sont toutes orientées dans la même direction et la situation où les fissures possèdent des orientations aléatoires. Dans une première étape, on rappelle brièvement les résultats disponibles concernant la modélisation des fissures non saturées par des cavités ellipsoïdales aplaties. A la fin de cette première partie, on complète les résultats déjà disponibles en étudiant l'influence de l'histoire de chargement sur la réponse de matériau. Dans une deuxième étape, on s'attache à valider une partie des résultats obtenus en utilisant une description des efforts capillaires dans les fissures par une précontrainte homogène en seréférant aux solutions analytiques exactes disponibles dans la littérature permettant de décrire le comportement d'une fissure isolée au sein d'une matrice élastique. Dans une troisième étape, on s'intéresse aux phénomènes de propagation des fissures en condition non saturée. Les lois de propagation sous critique et le phénomène de branchement des fissures sont également prises en compte dans cette approche. La dernière partie de la thèse concerne l'influence de la température sur le comportement des milieux poreux non saturés

La généralisation des charges non linéaires dans l'industrie cause de nombreuses perturbations dans les réseaux électriques. Ainsi, les courants harmoniques augmentent les pertes dans les transformateurs, par le biais des effets pelliculaire et de proximité. Dès lors, ceux-ci sont assujettis à un sur échauffement. Afin d'éviter une défaillance diélectrique ainsi qu'un vieillissement prématuré, la température du point chaud ne doit pas excéder celle spécifiée par la classe. Néanmoins, pour une charge donnée il est difficile de prévoir le comportement électromagnétique et thermique d'un transformateur. L'objectif de cette thèse est donc de caractériser ce fonctionnement pour se prémunir lors de la phase de conception des nuisances ultérieures. <br /><br />Une première modélisation s'appuie sur une méthode analytique et axisymétrique permettant le calcul de la distribution des densités de courant dans les enroulements. Puis, après une homogénéisation des conducteurs, une méthode de résolution par éléments finis aboutit à l'obtention de la température en régime permanent. Enfin, la méthode est vérifiée sur un transformateur de type sec et de puissance modérée, associé à un banc expérimental innovant. <br /><br />Un second modèle propose une nouvelle approche analytique et numérique. Tout d'abord, une méthode des éléments finis 3D (MEF) permet le calcul du champ magnétique. Ce calcul non linéaire est alors associé aux équations de circuit. Puis, une MEF thermique 3D fournie la distribution de température. Finalement, le modèle couplé est appliqué à un transformateur de distribution de petite puissance et instrumenté pour vérifier le suréchauffement engendré par des courants non sinusoïdaux.

On s'intéresse plus particulièrement à la modélisation du comportement d'un matériau hétérogène méso-fissuré (béton, roche,...), soumis à une sollicitation thermo-hydro-mécanique avec prise en compte du couplage géométrique. Pour conduire cette étude, on s'appuie notamment sur les approches micro-mécaniques du comportement des milieux méso-fissurés non saturés développées depuis quelques années au Laboratoire des Matériaux et des Structures du Génie Civil - Ur Navier - Université Paris Est. Le milieu fissuré non saturé traité ici est constitué d'une matrice solide homogène élastique linéaire et de fissures connectées saturées par deux fluides immiscibles : un liquide et un gaz séparés par une surface capillaire. La fissure est traditionnellement considérée comme une cavité ellipsoïdale (cas 3D) ou elliptique (cas 2D) dont le rapport d'aspect tend vers zéro. Deux morphologies typiques de matériau sont considérés dans ce travail : la situation où les fissures sont toutes orientées dans la même direction et la situation où les fissures possèdent des orientations aléatoires. Dans une première étape, on rappelle brièvement les résultats disponibles concernant la modélisation des fissures non saturées par des cavités ellipsoïdales aplaties. A la fin de cette première partie, on complète les résultats déjà disponibles en étudiant l'influence de l'histoire de chargement sur la réponse de matériau. Dans une deuxième étape, on s'attache à valider une partie des résultats obtenus en utilisant une description des efforts capillaires dans les fissures par une précontrainte homogène en seréférant aux solutions analytiques exactes disponibles dans la littérature permettant de décrire le comportement d'une fissure isolée au sein d'une matrice élastique. Dans une troisième étape, on s'intéresse aux phénomènes de propagation des fissures en condition non saturée. Les lois de propagation sous critique et le phénomène de branchement des fissures sont également prises en compte dans cette approche. La dernière partie de la thèse concerne l'influence de la température sur le comportement des milieux poreux non saturés / The main topic of my work is the development of a micromechanical model for the behaviour of unsaturated mesocracks in media (concrete, rock...) in which the thermo-hydro-mechanical loadingsand thermo-hydro-mechanical couplings are taken into account. For this, we used the micromechanical approach model of behaviour of cracked porous media recently developed at LMSGC. My thesis is focused on the equilibrium configurations of a porous material whose pore space is saturated by a vapour and a liquid phase. The behaviour of an elastic medium containing unsaturated mesocracks is studied in the framework of a micromechanical approach. The cracks are filled by two immiscible fluids, namely a liquid and a gas, separated by a capillary interface. Furthermore, it is assumed that the set of cracks constitutes a connected network ; the capillary pressure is uniform over a representative elementary volume. The cracks are modelled as flat oblate spheroid cavities. Several geometrical configurations of cracks in porous media are considered in the framework of Eshelby-based homogenization methods (parallel cracks, randomly oriented cracks). First, a previously developed model showed that when coupling between the deformation of the cracks and the capillary forces is taken into account, there is no more a one-to-one relationship between the loading parameters and the state-variables. Thus, we describe the loading history prescribed to the material in order to compute its response. Second, we validate these results referring to the exact solutions available in the literature to describe the behaviour of a unsaturated crack within an elastic matrix. Third, the description of crack propagation in unsaturated media is considered in the framework of linear elastic fracture mechanics. The phenomenon of subcritical crack growth due to stress corrosion cracking is taken into account in this approach. Mixed mode fracture in the plane is also examined. Finally, we are interested in the influence of the temperature on the behavior of unsaturated porous media in the framework of the micromechanical approach

Pour l’évaluation de la résistance ultime des ouvrages l’ingénieur de génie civil fait appel à différentes méthodes plus ou moins empiriques, dont de nombreuses manuelles, du fait de la lourdeur excessive des méthodes par éléments finis non-linéaires mises en œuvre dans les logiciels de calcul à sa disposition. Le calcul à la rupture, théorisé par J. Salençon, indique la voie de méthodes rigoureuses, tout à fait adaptées à cette problématique, mais dont la mise en œuvre systématique dans un logiciel a longtemps buté sur l’absence de méthodes numériques efficaces. Ce verrou de mathématique numérique a été levé récemment (Algorithme de point intérieur).Dans ce contexte l’objectif de la présente thèse est de mettre au point les méthodes permettant d’analyser, au moyen du calcul à la rupture, la capacité ultime d’éléments en béton armé tridimensionnels. Les deux approches du calcul à la rupture, que sont les approches statique et cinématiques, seront mises en œuvre numériquement sous la forme d’un problème d’optimisation résolu à l’aide d’un solveur mathématique dans le cadre de la programmation semi définie positive (SDP).Une large partie du travail sera consacré à la modélisation des différents matériaux constituant le béton armé. Le choix du critère pour modéliser la résistance du béton sera discuté, tout comme la méthode pour prendre en compte le renforcement. La méthode d’homogénéisation sera utilisée dans le cas de renforcement périodique et une adaptation de cette méthode sera utilisée dans le cas de renforts isolés. Enfin, les capacités et le potentiel de l’outil développé et mis en œuvre au cours de cette thèse seront exposés au travers d’exemples d’application sur des structures massives / To evaluate the load bearing capacity of structures, civil engineers often make use of empirical methods, which are often manuals, instead of nonlinear finite element methods available in existing civil engineering softwares, which are long to process and difficult to handle. Yield design (or limit analysis) approach, formalized by J. Salençon, is a rigorous method to evaluate the capacity of structures and can be used to answer the question of structural failure. It was, yet, not possible to take advantage of these theoretical methods due to the lack of efficient numerical methods. Recent progress in this field and notably in interior point algorithms allows one to rethink this opportunity. Therefore, the main objective of this thesis is to develop a numerical model, based on the yield design approach, to evaluate the ultimate capacity of massive (3D) reinforced concrete structural elements. Both static and kinematic approaches are implemented and expressed as an optimization problem that can be solved by a mathematical optimization solver in the framework of Semi-Definite Programming (SDP).A large part of this work is on modelling the resistance of the different components of the reinforced concrete composite material. The modelling assumptions taken to model the resistance of concrete are discussed. And the method used to model reinforcement is also questioned. The homogenization method is used to model periodic reinforcement and an adaptation of this technique is developed for isolated rebars. To conclude this work, a last part is dedicated to illustrate the power and potentialities of the numerical tool developed during this PhD thesis through various examples of massive structures