Le béton se présente comme un matériau constitué de granulats, jouant le rôle d'inclusions, enchâssés dans une matrice correspondant à la pâte de ciment hydraté. A l'échelle de la pâte, l'hydratation du ciment génère un milieu multiphasique, constitué d'un squelette solide et de pores remplis ou partiellement remplis d'eau selon leur taille. Le composant principal du squelette solide est le gel de C-S-H, les autres composants étant de nature cristalline. La qualification de gel du composant C-S-H est liée à sa nanostructure dont la schématisation la plus admise consiste en une phase aqueuse adsorbée, en sandwich avec des feuillets solides de nature cristalline. Il est bien admis que la structure du C-S-H est à l'origine de son comportement viscoélastique et donc de celui du béton. Ce comportement viscoélastique peut s'expliquer par un réarrangement de sa nanostructure sous l'effet des contraintes mécaniques appliquées à l'échelle macroscopique. La modélisation macroscopique du fluage du béton ne permet pas d'expliquer la variabilité du fluage d'une formulation de béton à une autre. En effet, les paramètres des modèles macroscopiques ne peuvent être identifiés que par l'analyse de résultats expérimentaux obtenus par des essais réalisés sur des éprouvettes de béton. Ces paramètres ne sont valables que pour une formulation donnée. L'identification de ces paramètres conduit donc à des programmes expérimentaux très coûteux et ne fournit pas suffisamment d'informations sur la sensibilité des paramètres macroscopiques à la variabilité des caractéristiques mécaniques et morphologies des constituants. Dans ce travail, on suppose qu'il existe une échelle microscopique à laquelle les mécanismes moteurs du fluage ne sont pas impactés par la formulation du béton. A cette échelle, celle du C-S-H, les propriétés viscoélastiques peuvent être considérées avoir un caractère intrinsèque. L'influence de la formulation ne concerne alors que les concentrations des différents hydrates. Trois approches, analytiques, semi-analytiques et numériques sont alors proposées pour estimer, par une homogénéisation multi-échelle, les propriétés viscoélastiques macroscopiques du béton à partir des propriétés de ses constituants ainsi qu'à partir de sa microstructure. Ces approches sont basées sur l'extension des schémas d'homogénéisation élastique au cas viscoélastique au moyen du principe de correspondance qui utilise la transformée de Laplace-Carson. Les propriétés effectives sont alors déterminées directement dans l'espace de Carson. Par la suite, celles dans l'espace temporel sont obtenues par la transformée inverse. Les approches proposées apportent des solutions aussi bien dans un cadre général que sous certaines hypothèses restrictives : coefficient de Poisson viscoélastique microscopique ou macroscopique constant, module de compressibilité constant. Sur le plan théorique, deux schémas d'homogénéisation ont été étudiés : le schéma de Mori-Tanaka, basé sur le problème de l'inclusion d'Eshelby, et le schéma auto-cohérente généralisé basé sur la neutralité énergique de l'inclusion. Les résultats obtenus montrent que sous ces hypothèses restrictives, le spectre macroscopique se présente comme une famille de sous ensembles de temps caractéristiques bornés par les temps caractéristiques microscopiques. Par ailleurs, les propriétés thermodynamiques, de croissance monotone et de concavité, des fonctions de retard macroscopiques ne sont préservées par l'homogénéisation que sous certaines conditions de compatibilité des spectres microscopiques. Sur le plan pratique, les méthodes développées ont été appliquées pour construire la complaisance de fluage propre macroscopique du béton en connaissant les données communes de toutes sortes de bétons et celles correspondant à une formulation donnée. Les résultats expérimentaux disponibles sont alors exploités pour analyser le caractère intrinsèque des propriétés viscoélastiques à l'origine du fluage du béton / Concrete can be considered as a multiscale heterogeneous material in which elastic inclusions, i.e., aggregates, are embedded in a viscoelastic matrix which corresponds to the hardened cement paste. At the cement paste scale, the hydration process leads to a partially saturated porous multiphase medium whose solid skeleton is mainly constituted of the C-S-H gel, the other hydrates being crystalline solids. The gel-like properties of the C-S-H are essentially related to the existence of strongly adsorbed structural water in the interlayer nanospace. It is well admitted that this feature is the microscopic origin of the creep mechanism of C-S-H and therefore of concrete at the macroscale. This viscoelastic behavior can be explained by a rearrangement of the nanostructure of the C-S-H component due to a macroscopically applied loading. Therefore, the macroscopic modeling of concrete creep does not allow explaining the variability of creep with respect to the concrete mix. Indeed, the identi?cation of the viscoelastic parameters of macroscopic models requires performing creep tests at the concrete scale. Thus, the validity of the identi?ed parameters is limited to the concrete mix under consideration and the identi?cation process does not permit to carry out the influence of the microstructure and the mechanical constituent properties on the macroscopic creep. In this work, we assume that there exists a microscopic scale at which the driving creep mechanisms are not affected by the concrete mix. At this scale, the C-S-H one, the viscoelastic properties can be considered to be intrinsic. The concrete mix will essentially influence the volume fraction of the hydrates. Analytical, semi-analytical and numerical approaches are developed in order to derive, by multi-scale homogenization, the macroscopic viscoelastic properties of concrete on the basis of the knowledge of the microscopic properties, of the constituents together with the knowledge of the mix-dependent microstructure. These approaches consist in extending results of elastic homogenization schemes to viscoelasticity with the help of the correspondence principle using the Laplace-Carson transform. Therefore, the effective properties are obtained in a straightforward manner in the Carson transform space, which requires the transform inversion in order to determine the time evolution of these properties. The proposed models permit to answer to this requirement, both under some restrictive assumptions (constant Poisson's ratio at either microscopic or macroscopic scale, constant bulk microscopic modulus) and in a general way. At the theoretical level, two homogenization schemes are investigated : the Mori-Tanaka scheme based on the Eshelby solution and the Generalized Self Consistent scheme, based on the energetic neutrality condition of the inclusion. The obtained results show that under the aforementioned restrictive conditions, the macroscopic creep spectrum presents as a family of sets of retardation times which are bounded by two successive microscopic characteristic times. Furthermore, the monotonically increasing and concavity features of the macroscopic creep compliance, which derive from thermodynamic considerations, are not always preserved, except if some compatibility conditions are ful?lled by the microscopic spectrums. From the practical point of view, these methods developed are used to establish the complaisance creep function of concrete from the knowledge of the common data of all types of concrete, as well as the other parameters that characterize a given formulation. Hence, by a back-analysis of creep tests obtained on cement paste or on concrete, the intrinsic basic creep properties of the gel C-S-H are estimated

L’augmentation de la puissance électrique des composants électroniques pose le problème de la dissipation de la chaleur générée. Les boîtiers électriques doivent permettre la dissipation de cette chaleur en conservant une isolation électrique. La solution retenue pour évacuer la chaleur par transfert thermique consiste en matériaux composites dont les renforts par leur structuration vont améliorer la conductivité thermique. Des composites à matrice polymère ont été choisis pour leur aptitude de mise en forme. La conductivité thermique et l’isolation électrique sont assurées par des charges céramiques. Les méthodes d’homogénéisation donnent des pistes d’amélioration du comportement de composites en fonction des propriétés de leurs constituants, de leur géométrie et de leur distribution. Elles fournissent ainsi une formulation optimisée de matériaux répondant à certaines caractéristiques issues de cahiers des charges émanant du partenaire industriel (Institut de Soudure). La conductivité thermique attendue des composites impose une forte fraction volumique de charges pour compenser le caractère isolant de la matrice polymère. Des méthodes d’homogénéisation ont été développées pour prédire la conductivité thermique effective pour de forts taux de charges (supérieur à 20%) et des contrastes élevés de conductivité thermique. La présence d’une interphase provenant d’incompatibilités fortes entre les composants doit également être modélisée / The increase of electronic components in the integrated circuits and the required electrical power set the question of the dissipation of the heat generated. The electrical box must favor the heat dissipation while maintaining electrical insulation. The solution chosen to transfer the heat is to develop composite materials whose reinforcements by their structure will improve the thermal conductivity. Polymer-based composite materials were chosen for their building ability. Thermal conductivity and electrical insulation are insured by ceramic reinforcements. The homogenization methods allow to improve the composites’ design according to the properties of their constituents, their geometry and their distribution. They thus provide an optimized formulation of materials satisfying the characteristics emanating from the industrial partner (‘Institut de Soudure’). The expected thermal conductivity of the composites imposes a high volume fraction of reinforcements to counterbalance the insulating polymer matrix. Homogenization methods have been developed to provide predictions of effective thermal conductivity for high (greater than 20%) reinforcement rates and high thermal conductivity contrasts. The presence of an interphase resulting from strong physico-chemical incompatibilities between the components must also be modeled

Cette thèse est consacrée au développement d'une méthode numérique visant à modéliser des écoulements dans des milieux poreux fissurés, ainsi qu'à déterminer leur perméabilité effective à partir des avancements théoriques récents. En parallèle, elle a été aussi l'occasion de continuer sur la voie théorique et d'obtenir de nouveaux résultats sur ce plan.Les équations générales gouvernant l'écoulement dans de tels matériaux sont rappelées, et plus particulièrement, la conservation de la masse à l'intersection entre fissures est établie explicitement [132, 139]. Des solutions générales du potentiel sont proposées sous la forme d'une équation intégrale singulière décrivant l'écoulement dans et autour des fissures entourées par une matrice infinie soumise à un champ lointain [136, 139]. Ces solutions représentent le champ de pression dans le milieu infini en fonction de l'infiltration dans les fissures qui prennent en compte complètement l'interaction et l'intersection entre fissures. En considérant le problème d'une fissure super-conductrice, des solutions analytiques de l'écoulement ont été développées et ces solutions sont utilisées comme cas de référence pour valider la solution numérique. De plus, la solution théorique de ce problème dans le cas 3D permet de comparer le modèle d'écoulement de Poiseuille dans une fissure elliptique d'épaisseur nulle et le modèle d'inclusion ellipsoïdale aplatie soumise à l'écoulement de Darcy [140]. Des outils numériques ont été développés en se basant sur la méthode des équations intégrales singulières afin de résoudre les équations générales du potentiel [132, 180]. Cela permet, d'une part, de modéliser l'écoulement stationnaire dans un domaine poreux contenant un grand nombre de fissures et, d'autre part, de proposer une solution semi-analytique de l'infiltration dans une fissure isolée dépendant de la perméabilité de la matrice, de la conductivité de la fissure et de la variable géométrique de la fissure. Cette dernière est l'élément essentiel pour déterminer de la perméabilité effective d'un milieu poreux fissuré en utilisant des schémas d'homogénéisation. Ce modèle auto-cohérent révèle un seuil de percolation qui est alors applicable pour l'estimation de la perméabilité effective d'un matériau contenant un grand nombre de fissures. L'approche par sous-structuration permet l'extension de la solution générale du potentiel, écrite pour un domaine infini, à celle d'un domaine fini [181]. Une solution analytique de l'écoulement dans et autour d'une fissure partiellement saturée est établie, fondée sur la solution analytique pour la fissure super-conductrice. Celle-ci est alors utilisée pour estimer la perméabilité effective d'un milieu poreux fissuré non-saturé [141]. Le modèle de la perméabilité effective est appliqué dans le contexte du stockage géologique du CO2 en vue d'étudier le comportement d'une zone de faille constituée par un noyau argileux et des zones fissurées. La pression d'injection provoque l'augmentation de la pression interstitielle dans le réservoir. Cette surpression engendrée dans le réservoir peut affecter la perméabilité de zones fissurées ce qui conduit à des phénomènes hydromécaniques couplés. Les résultats de simulations numériques permettent d'évaluer, d'une part, le risque de la remontée de la saumure à l'aquifère supérieur, et d'autre part, le risque de l'initiation d'une rupture sur le plan de la faille / This thesis aims to develop a method for numerical modelling of fluid flow through fractured porous media and for determination of their effective permeability by taking advantage of recent results based on formulation of the problem by Singular Integral Equations. In parallel, it was also an occasion to continue on the theoretical development and to obtain new results in this area. The governing equations for flow in such materials are reviewed first and mass conservation at the fracture intersections is expressed explicitly. Using the theory of potential, the general potential solutions are proposed in the form of a singular integral equation that describes the steady-state flow in and around several fractures embedded in an infinite porous matrix under a far-field pressure condition [136, 139]. These solutions represent the pressure field in the whole body as functions of the infiltration in the fractures, which fully take into account the fracture interaction and intersections. Closed-form solutions for the fundamental problem of fluid flow around a single fracture are derived, which are considered as the benchmark problems to validate the numerical solutions. In particular, the solution obtained for the case of an elliptical disc-shaped crack obeying to the Poiseuille's law has been compared to that obtained for ellipsoidal inclusions with Darcy's law [140].The numerical programs have been developed based on the singular integral equations method to resolve the general potential equations [132, 180]. These allow modeling the fluid flow through a porous medium containing a great number of fractures. Besides, this formulation of the problem also allows obtaining a semi-analytical infiltration solution over a single fracture depending on the matrice permeability, the fracture conductivity and the fracture geometry. This result is the important key to upscalling the effective permeability of a fractured porous medium by using different homogeneisation schemes. The results obtained by the self-consistent scheme have been in particular established. The multi-region approach can be used to extend the general potential solution written for the infinite domain to that for a finite domain [181]. A closed-form solution for flow in and around a single partially saturated fracture, surrounded by an infinite matrix subjected to a far-field condition, is also derived combining the solutions for a superconductive fracture and for an imprevious fracture. This solution is then employed to estimate the effective permeability of unsaturated fractured porous media [141].The effective permeability model is applied to study the hydromechanical behaviour of a fault zone constituted by a clay core surrounded by fractured zones in the context of CO2 geological storage. The pressure injection induces an overpressure in the reservoir that may affect the permeability of the fractured zones leading to complexe coupled hydromechanical phenomena. The simulation results allow evaluating the risk of leakage of the reservoir brine to higher aquifers as well as the risk of fault reactivation