L’expansion actuelle du secteur de la haute vitesse ferroviaire vers des vitesses d’exploitation et chargements plus importants conduit à la nécessité d’améliorer les systèmes de freinage garantissant leurs performances malgré l’énergie supplémentaire à dissiper.La méthodologie actuelle de développement de ces matériaux étant principalement basée sur un retour d’expérience essai/erreur, l’objectif vise donc à développer une méthodologie alternative incorporant des modèles théoriques et numériques pour la conception de matériaux de friction. Cependant, les difficultés s’avèrent multiples avec des interactions complexes entre les aspects tribologiques, thermiques, mécaniques, et physico-chimiques, mais également entre les différentes échelles allant de la surface au système. Ce travail se concentre donc sur les aspects thermomécaniques dans le but de développer un modèle numérique prédictif intégrant une évolution du matériau pouvant affecter les performances de freinage. Dans un premier temps, une méthodologie originale est développée pour caractériser le matériau de friction soumis à des sollicitations réalistes et permettant d’identifier le comportement du matériau. La caractérisation fut réalisée pour différents historiques de freinage. Des modèles de comportement thermo-élasto-plastique ont été ainsi établis. La deuxième étape de ce travail porte sur le développement d’un modèle d’éléments finis de freinage intégrant les comportements matériaux précédemment identifiés ainsi que leur évolution avec l’historique de chargement. Les résultats illustrent l’impact de l’évolution du matériau de friction sur les performances de freinage. / The railway high speed sector is currently expanding with an increase of the maximum operating speed and load leading to the necessity of improving the capacity of the disc braking to ensure performances despite the additional energy to dissipate. The current methodology used for the development of these material is based on trial-error experimental feedback. The goal is to develop alternative methodologies, with theoretical and numerical models for designing these components. Difficulties are the complex interactions between tribological, thermal, mechanical, chemical effects and interactions between scales from the surface to the system. This work is focused on the thermomechanical aspects with the challenge of developing a realistic numerical model and considering the material evolution assuming that this evolution affects the braking performances. In the first step, an original methodology has been developed to characterize the friction material submitted to realistic solicitations and allowing identifying bulk material behavior. Characterization has been done for different braking loading histories. Corresponding thermo-elastic-plastic behavior models have been proposed. The second step of this work is the development of a finite element model of the brake system, including the sintered material models previously identified and the evolution with loading history. The numerical results illustrate the impact of the friction material evolution regarding the braking performances. They are also compared to experimental tests carried out on a real scale braking bench.

Cette étude a été effectuée à la suite du développement de lois et d'outils applicables à la modélisation numérique du comportement élastoviscoplastique de matériaux cristallins : des lois de comportement de monocristaux, des lois de transition d'échelle pour les modèles d'homogénéisation, un code de calcul Eléments Finis adapté au calcul parallèle et un programme de génération de microstructures polycristallines 3D. Disposant de ces éléments, nous avons étudié le comportement de polycristaux 3D en petites déformations, aux échelles macroscopique, intergranulaire et intragranulaire.<br />Le milieu polycristallin est décrit par des polyèdres de Voronoï, donnés sous la forme d'un fichier de voxels (L. Decker, D. Jeulin, ENSMP). L'implémentation de la méthode FETI dans le code EF ZéBuLoN (F. Feyel, S. Quilici, ENSMP-ONERA) permet la résolution en parallèle de problèmes à très grand nombre de degrés de liberté. Ainsi nous avons accès à un nombre illimité de réalisations de microstructures et nous pouvons faire figurer suffisamment d'éléments dans un maillage pour que soit possible la description des champs intragranulaires dans un polycristal 3D.<br />Pour commencer nous montrons les spécificités de notre approche par rapport aux travaux de modélisation de la plasticité cristalline. La première partie de l'exploitation des outils a consisté à analyser la sensibilité des résultats aux données de la modélisation (nombre d'éléments, nombre de grains, réalisation de microstructure . . . ) afin d'établir une configuration de calcul valable pour des simulations sur un Volume Elémentaire Représentatif de polycristal isotrope. En seconde partie nous mettons en évidence l'hétérogénéité de comportement inter- et intragranulaire et l'apport de la méthode par rapport à une démarche autocohérente. Ceci est complété par une analyse de l'influence des joints de grain et des conditions aux limites sur la réponse d'un essai en traction simple, aux différentes échelles de la modélisation. Nous caractérisons ainsi un effet local et un effet moyen pour tous les grains, en fonction de la distance à un joint ou à un bord. En annexe sont donnés les résultats de simulations obtenus avec un modèle non-local des milieux de Cosserat (S. Forest, ENSMP) qui ont permis de quantifier un effet de taille de grain sur le comportement effectif de polycristaux.

L'objectif de la présente thèse et de développer un cadre d'étude et de simulation apte à prédire les propriétés effectives de diffusion dans les matériaux cimentaires saturés. Le principe suivi, fondé sur les approches multi-échelles séquencées, permettra de modéliser chaque échelle d'intérêt, depuis le gel des C-S-H, jusqu'au matériau béton, et ainsi via des simulations de prédire les coefficients de diffusion effectifs.Le cadre d'étude est donc fondé sur les points clés suivant :- Des morphologies de type "matrice - inclusion" sont générées via des excursions de champs aléatoires corrélés tri-dimensionnels. Cette approche, combinée avec quelques opérations morphologiques, permet la génération de morphologies inclusionnaires 3D suivant des distributions de tailles imposées. Il est ensuite possible de faire évoluer ces morphologies, depuis un état initial, en fonction d'un modèle prescrit (par exemple en termes de fraction volumiques), via le seuil de l'excursion.- Les cibles morphologiques de chaque échelle considérées tentent de refléter les résultats expérimentaux et modèles les plus communément admis. En supplément, à l'échelle de la pâte de ciment, un modèle "simple" d'hydratation est développé afin de prédire l'évolution de celle-ci dans le temps en termes de fractions volumiques. Le choix a été fait d'utiliser pour fondement un modèle simple (Jennings & Tennis), incorporant des paramètres essentiels, tels que la composition minéralogique du ciment ou le rapport e/c. Des modifications mineures ont été apporté au modèle, afin de prendre en compte de récents résultats expérimentaux, concernant principalement la fin de l'hydratation ainsi que sa cinétique.- Le cadre d'homogénéisation a été construit pour être utilisable à toutes les échelles considérées, et il permet ainsi de prédire les propriétés effectives de diffusion avec des considération énergétiques d'équivalence entre échelles. De façon analogue à la théorie de Hill pour l'élasticité, le cadre d'étude assure une séparation d'échelle cohérente vis-à-vis des propriétés effectives de diffusion.- Enfin, d'un point de vue numérique, les simulations utilisent les méthodes E-FEM, couplées à l'utilisation de maillages non-adaptés sur lesquels les morphologies sont "projetées". De cette façon, les discontinuités géométriques (discontinuité de propriétés matériaux) sont sont assurés par les élément finis enrichis, et ne requièrent pas d'éléments joint, ou de "moyennage" de propriétés.Les résultats de chaque échelle sont utilisés comme paramètres d'entrés pour les échelles suivantes. Depuis l'échelle des hydrates (C-S-H), jusqu'à l'échelle du béton, les coefficients de diffusion effectifs sont ainsi estimés. / The aim of the present Ph.D. is to develop a methodology and a simulation framework able to predict as accurately as possible effective properties regarding diffusion in saturated cement based materials. A sequenced multiscale framework is developed to perform numerical homogenization on diffusive properties, and thus predict effective properties. The general outline is to model every scale in concrete-like materials, from the C-S-H gel, to the concrete scale, and with the help of a proper framework, be able to predict effective properties. The methodology relies on the following key points for every considered scale :- Inclusion-matrix like morphologies are generated through level set methods applied to multi-dimensional correlated random fields. This approach, together with few morphological operations allows to generate 3D morphologies with given particle size distributions (PSD), or pore size distributions. Then, it is also possible, from the initial state, to make those morphologies evolve according to any descriptive model (e.g. in terms of volume fraction), via the level set.- The morphological description of each scale is chosen according to the most acknowledged experimentations and models. In addition, at the cement paste scale, a hydration model was developed to predict the evolution of the cement paste through time, in terms of volume fractions. The choice was made to keep a simple model (Jennings & Tennis), with few but essential input parameters (mineral cement composition, w/c ratio, PSD). It has only been modified to include recent experimental results, such as the maximal hydration degree, or the hydration rates.- The upscaling framework was built to suit every considered scale, and allows to assess effective properties based on energy conservation principles. Similar to Hill's theory for elasticity, the framework is build to ensure to proper scale separation regarding effective diffusion coefficients.- Lastely, the numerical context is based on E-FEM techniques, which allows to use unstructured meshes, on which the morphologies are projected. This way, any material discontinuity in dealt within the FE, and thus does not require averaging properties, or joint elements. Results from every scale will be used as inputs at the higher scales, and effective diffusion coefficients are thus estimated.

Le présent travail propose et développe une méthodologie complète d’identification des propriétés mécaniques locales d’un matériau hétérogène à l’échelle des phases constitutives. Il s’agit d’une combinaison de compétences très diverses, à la fois en théorique, en simulation numérique et en expérimentation. Plus précisément, la partie théorique concerne la détermination des relations nano-micro pour le module d’indentation homogénéisé par des techniques de changement d’échelle; la partie numérique basée sur la théorie du calcul à la rupture est réalisée en vue de trouver de ces dernières relations applicables pour la dureté; et la dernière partie est effectuée pour récupérer les propriétés homogénéisées par la voie expérimentale à l’aide de la technique de nano-indentation. L’étude expérimentale de la thèse est pour l’objectif de déterminer des propriétés d’indentation de différents échantillons de pâte de ciment. Un programme expérimental complet, est développé, qui permet de caractériser des phases principales à l’échelle micrométrique de ce matériau, parmi lesquelles nous nous intéressons surtout à celles plus importantes correspondantes à des phases de la matrice C-S-H. La modélisation du problème lié à l’enfoncement d’une pointe d’indentation dans un matériau est étudiée. Pour cela, la première voie, basée sur l’approche cinématique du calcul à la rupture, consiste à tenter de construire des mécanismes de ruine analytiquement, puis à les faire évoluer en fonction du changement de la géométrie initiale, afin d’obtenir la charge de ruine correspondante. La seconde voie consiste ensuite à suivre la même approche, mais en construisant numériquement ces mécanismes de ruine. La charge obtenue dépend naturellement des paramètres de critères retenus, que l’on détermine grâce à la combinaison avec les résultats expérimentaux. Les critères de résistance de Von-Mises et de Tresca valables pour des matériaux purement cohérents ainsi que celui de forme elliptique sont examinés dans ce travail / The present work proposes and develops a complete methodology for identifying the local mechanical properties of a heterogeneous material at the scale of the constitutive phases. It is a combination of very diverse skills in theory, in numerical simulation and in experimentation. More precisely, the theoretical part concerns the determination of the nano-micro relations for the indentation module; the numerical part based on the yield design theory is carried out to find the last relations applicable for the hardness; and the last part is performed to obtain homogenized properties by the experimental way using the nano-indentation technique. The experimental study of the thesis is for the purpose of determining indentation properties of different cement paste samples. A complete experimental program, is developed, which allows characterizing the main phases at the micrometric scale of this material, among which we are mainly interested in the C-S-H matrix phases. The modeling of the problem related to the penetration of an indentation point into a material is studied. For this, the first way, based on the kinematic approach of the yield design theory, consists in trying to construct ruin mechanisms analytically, then to make them evolve according to the change of the initial geometry, in order to obtain the corresponding ultimate load. The second way is then to follow the same approach, but by building numerically these ruin mechanisms. The obtained load depends naturally on the retained criteria parameters, which are determined by the combination with the experimental results. The Von-Mises and Tresca strength criteria for purely coherent materials as well as the elliptical one are examined in this work

En sciences et dans l'industrie pour limiter le nombre de maquettes expérimentales dans les projets R&D afin de mieux comprendre et de bien interpréter les phénomènes ultrasonores complexes observés sur site, des simulations de contrôles ultrasonores sont effectuées. Ces simulations sont d'autant plus réalistes que la description des structures à contrôler est précise, en particulier au niveau des constantes d'élasticité et d'atténuation intrinsèque. Les objectifs de cette étude sont d'améliorer d'une part les connaissances sur l'influence des caractéristiques métallurgiques des matériaux anisotropes et hétérogènes sur la propagation ultrasonore et d'autre part les performances des codes de calcul (logiciel ATHENA d'EDF) qui nécessitent de disposer des données d'entrée pertinentes, notamment en ce qui concerne les constantes d'élasticité et l'atténuation ultrasonore. Cette étude est dédiée à la caractérisation des matériaux à gros grains, comme les aciers austéno-ferritiques moulés par centrifugation et les soudures en acier inoxydable austénitique ou en alliages à base nickel. Un système expérimental unique permettant de mesurer les constantes d'élasticité et l'atténuation en incidence oblique à été mis au point. Le point fort de ce dispositif est qu'il permet de travailler au-delà de l'angle critique longitudinal et donc de mesurer les propriétés d'atténuation des ondes transversales. Les constantes d'élasticité sont déduites des vitesses ultrasonores à partir d'un processus d'optimisation basé sur la résolution de problème inverse. Nous avons montré les potentialités d'algorithmes d'optimisation globaux tels que les algorithmes génétiques moins susceptibles de converger vers des minima locaux de la fonction à minimiser. Les résultats obtenus à partir des mesures expérimentales sont en accord avec la littérature. Des résultats de l'atténuation des ondes longitudinales et transversales par décomposition du faisceau en spectre d'ondes planes sont présentés.

Prévoir où, quand et comment les matériaux cassent est une problématique qui occupe scientifiques et ingénieurs depuis des siècles. Ce problème est rendu complexe par le fait que la concentration des contraintes en pointe de fissure lie intimement le comportement observé à l'échelle macroscopique aux inhomogénéités de microstructure à des échelles très fines. Ceci induit une dynamique de fissuration erratique, composée d'événements d'endommagement rapides et imprévisibles séparés de périodes calmes (e.g. dynamique des tremblements de terre le long des failles). Par essence, ces aspects statistiques ne peuvent pas être traités avec l'approche de la mécanique des milieux continus traditionnels. Dans un premier temps, nous tentons d'appréhender ce problème au travers d'une expérience modèle qui consiste à faire propager une fissure dans une roche artificielle dont nous contrôlons la microstructure. La vitesse de chargement du système de fracture est réglable sur une large gamme de valeurs. La vitesse de fissuration et l'énergie mécanique sont enregistrées en temps réel. En parallèle, l'émission acoustique associée aux événements de fracture ainsi que leur localisation sont mesurées via des capteurs piézoélectriques, puis analysées comme cela est communément fait en sismologie. Ces expériences nous permettent de caractériser quantitativement la dynamique intermittente de la fissuration. Elles montrent qu'un certain nombre des lois empiriques observées en géophysique sur la sismicité (loi de Richter-Gutenberg, d'Omori, de Voight, d'Utsu...) se retrouvent dans notre système modèle. Dans un deuxième temps, nous adressons ce problème théoriquement et numériquement, en identifiant le phénomène de fracture dans les matériaux hétérogènes avec celui de la propagation d'une ligne élastique sur un potentiel aléatoire 2D. Ceci permet de déterminer quantitativement, en termes de vitesse de chargement, de tailles des hétérogénéités, de propriétés du matériau, et de géométrie de structure, quand la dynamique de fissuration est régulière et compatible avec l'approche ingénieur des milieux continus, et quand elle devient erratique et nécessite une approche statistique. Dans ce dernier cas, nous caractérisons la statistique de de cette dynamique et relions celle-ci aux paramètres de l'expérience.

En sciences et dans l’industrie pour limiter le nombre de maquettes expérimentales dans les projets R&D afin de mieux comprendre et de bien interpréter les phénomènes ultrasonores complexes observés sur site, des simulations de contrôles ultrasonores sont effectuées. Ces simulations sont d’autant plus réalistes que la description des structures à contrôler est précise, en particulier au niveau des constantes d’élasticité et d’atténuation intrinsèque. Les objectifs de cette étude sont d’améliorer d’une part les connaissances sur l’influence des caractéristiques métallurgiques des matériaux anisotropes et hétérogènes sur la propagation ultrasonore et d’autre part les performances des codes de calcul (logiciel ATHENA d’EDF) qui nécessitent de disposer des données d’entrée pertinentes, notamment en ce qui concerne les constantes d’élasticité et l’atténuation ultrasonore. Cette étude est dédiée à la caractérisation des matériaux à gros grains, comme les aciers austéno-ferritiques moulés par centrifugation et les soudures en acier inoxydable austénitique ou en alliages à base nickel. Un système expérimental unique permettant de mesurer les constantes d’élasticité et l’atténuation en incidence oblique à été mis au point. Le point fort de ce dispositif est qu’il permet de travailler au-delà de l’angle critique longitudinal et donc de mesurer les propriétés d’atténuation des ondes transversales. Les constantes d’élasticité sont déduites des vitesses ultrasonores à partir d’un processus d’optimisation basé sur la résolution de problème inverse. Nous avons montré les potentialités d’algorithmes d’optimisation globaux tels que les algorithmes génétiques moins susceptibles de converger vers des minima locaux de la fonction à minimiser. Les résultats obtenus à partir des mesures expérimentales sont en accord avec la littérature. Des résultats de l’atténuation des ondes longitudinales et transversales par décomposition du faisceau en spectre d’ondes planes sont présentés. / In industry, to limit the number of experimental models in R&D projects, to better understand and to well interpret the complex ultrasonic phenomena observed du ring controls on site, simulations of ultrasonic controls are carried out. These simulations are all the more realistic as the description of structures to control is accurate, especially in terms of elastic constants, and intrinsic attenuation. The objectives of this study are firstly to improve knowledge about the influence of the metallurgical properties of anisotropie and heterogeneous materials on the ultrasonic propagation and secondly the performance of the computation codes (software ATHENA EDF) which need to have the relevant inputs, particularly as regards the elastic constants and ultrasonic attenuation. This study is dedicated to the characterization of coarse materials such as austenitic-ferritic steel centrifugally cast and the welding in steel austenitic stainless or in alloy nickel-based. A unique experimental system for measuring the elastic constants and attenuation at oblique incidence has been developed. The strong point of this device is that it can work beyond the longitudinal critical angle and thus to measure the attenuation properties of transversal waves. The elastic constants are deduced from ultrasonic speed from an optimization process based on the resolution of Inverse problems. We have shown the potential of global optimization algorithms such as genetic algorithms Jess likely to converge to local minima of the function to minimize. The results obtained from experimental measurements are in agreement with literature. Results of the attenuation of the longitudinal and transverse waves by beam decomposition into spectrum of plane waves are represented.

La méthode d'homogénéisation en calcul à la rupture permet de calculer les propriétés de résistance à l'échelle macroscopique d'un matériau hétérogène périodique (matériaux composites, sols renforcés) à partir de la solution d'un problème de calcul à la rupture défini sur la période de base. Une méthode numérique spécifique destinée à résoudre un tel problème est mise au point. Elle est fondée sur la mise en oeuvre de l'approche cinématique du calcul à la rupture et conduit : 1) à la recherche du minimum d'une fonction convexe d'un nombre fini de paramètres scalaires pour laquelle un algorithme de résolution est proposé. 2) à une évaluation par l'extérieur du convexe de résistance macroscopique qui se révèle d'autant plus précise que le nombre de paramètres de minimisation est élevé. On procède à de nombreuses applications de cette méthode à des problèmes de "contrainte plane" (plaques renforcées ou perforées chargées dans leur plan) et son extension prometteuse à des problèmes en "déformation plane" est esquissée.

La première partie du mémoire est consacrée à une étude bibliographique. Après un rapide aperçu historique des travaux portant sur les milieux à microstructure, une étude de la théorie des matériaux simples de W. Noll est présentée, ainsi qu'un exposé des travaux de G. Capriz sur les matériaux à microstructure. On introduit ensuite les notions d'univers d'observation "profond", de microplacement et de macroplacement dans cet univers d'un ensemble de particules. Les propriétés de ces placements conditionnent les définitions des milieux monoparticulaires, multiparticulaires, multiconstituants ou multiphasiques que l'on donne. On décrit les propriétés cinématiques et dynamiques de ces milieux. On propose l'introduction d'un champ de fonctions de phase pour traiter les problèmes de changements de phases. Deux approches complémentaires pour obtenir les équations d'équilibre d'un milieu multiparticulaire à microstructure sont ensuite mises en oeuvre, en s'intéressant plus particulièrement aux milieux de Cosserat multiparticulaires (distribution continue d'orientation de trièdres au même point de l'espace physique). Le chapitre IV consiste en un exposé complet d'une méthode de changement d'échelle apparentée aux outils classiques de l'homogénéisation que nous nommons globalisation. Cette méthode permet, par défocalisation d'un milieu à microstructure "classique", de déduire les équations d'équilibre macro des équations d'équilibre du milieu initial. L'intérêt réside dans l'identification physique des grandeurs constitutives du milieu macro à partir des grandeurs constitutives du milieu micro. La démarche est menée à terme dans le cadre général des microstructures sur lesquelles l'action du groupe des rotations de IR3 est transitive. La fin du mémoire apporte quelques indications sur l'ébauche d'une méthode de détermination du comportement d'un milieu multiparticulaire à partir de la globalisation des principales équations de bilan du milieu micro.

La durabilité des structures en béton est désormais intégrée dans la démarche de conception des ouvrages de Génie Civil. En effet, quel que soit le type de sollicitation (mécanique, thermique, hydrique) une fissuration est susceptible de se produire risquant d'impacter la durée de vie de l'ouvrage par la pénétration d'agents agressifs. L'aptitude au service peut elle-même être affectée pour les structures où une étanchéité est requise (enceinte de confinement de centrales nucléaires, réservoirs de gaz naturel liquéfié, barrages, stockages des déchets radioactifs ou de CO2, etc.). Dans ce contexte industriel, la prédiction du débit de fuite traversant des éléments composés de matériaux à base cimentaire est donc un enjeu scientifique et industriel majeur. Pour parvenir à cet objectif de simulation numérique, il est nécessaire de mettre en place un couplage hydro-mécanique. L'anisotropie de la fissuration induite par les sollicitations mécaniques complexes conduit à un tenseur de perméabilité macroscopique anisotrope. La détermination de ce tenseur est un enjeu important dans l'objectif de mener des calculs à l'échelle macroscopique avec des modèles phénoménologiques. De plus, les calculs de perméabilité sont un moyen de comparer les volumes fissurés obtenus par les différents modèles mécaniques. La modélisation de la fissuration pour les matériaux quasi-fragiles hétérogènes à l'échelle mésoscopique tels que le béton est complexe et suivant les approches utilisées, les résultats peuvent fortement varier. C'est pourquoi l'étude numérique proposée dans la thèse comporte une comparaison entre deux approches mécaniques : - une première basée sur une modélisation mécanique de type E-FEM (Embedded Finite Element Method) [Benkemoun et al., 2010] - - une seconde basée sur une modélisation mécanique d'endommagement [Mazars, 1984] régularisée en énergie de fissuration [Hillerborg et al., 1976]. Le travail numérique associé à cette thèse consiste donc à développer un modèle couplant de manière faible un modèle mécanique à un modèle de transfert en 3D à l'échelle mésoscopique. En se basant sur le concept de « double porosité », la perméabilité du milieu fissuré est vue comme la combinaison d'une perméabilité diffuse et isotrope (liée au réseau poreux initial du béton et à son degré de saturation) et d'une perméabilité « discrète » et orientée au sein des fissures (le calcul de cette dernière étant basé sur les ouvertures de fissures données par le modèle mécanique et sur les équations de la mécanique des Navier-Stokes en régime permanent). La comparaison des résultats obtenus sur différents résultats expérimentaux issus de la littérature (un tirant traversé par de l'eau [Desmettre et Charron, 2011] et un élément structurel traversé par de l'air sec [Nahas et al., 2014]) permet de comparer la pertinence des deux modèles mécaniques utilisés ainsi que l'approche utilisée pour estimer le débit traversant des éléments en béton fissurés. / The durability of concrete structures is nowadays fully integrated in the civil engineering constructions design process. Whatever the loading is (mechanical, thermic, hydric), cracks may appear and impact the structure lifespan by the infiltration of aggressive agents. The serviceability can be directly impacted for the structures playing an air/water tightness role (containment building nuclear power plants, liquefied natural gas storage tanks, dams, radioactive waste disposal, etc.). The prediction of the flow going through elements composed of a cementitious material is therefore a major scientific and industrial issue. To achieve this goal, a hydro-mechanical coupling must be implemented. The anisotropic cracking induced by complex mechanical loadings leads to an anisotropic macroscopic permeability tensor. This tensor computation is an important issue dealing with phenomenological models for macroscopic problems. The cracking modelling of quasi-brittle materials, heterogeneous at the mesoscopic scale like concrete, is complex and different mechanical approaches can lead to various results. Therefore, permeability calculations are an elegant way to examine cracking patterns obtained with several mechanical models. Consequently, this study compares two mechanical approaches: - the first one is based on an Embedded Finite Element Method (E-FEM) mechanical model [Benkemoun et al., 2010] - - the second one is based on a damage mechanical model [Mazars, 1984] regularised by the fracture energy of the material [Hillerborg et al., 1976]. This thesis presents a hydro-mechanical approach weakly coupling a mechanical model with a permeation model in 3D at the mesoscopic scale. This work is based on the “double porosity” concept splitting the permeability into two parts: the first one is isotropic and corresponds to flows within the porosity of the material- the second one, based upon a set of cracks with different orientations and openings, is anisotropic. For the latter, each crack is a path for mass flow according to the fluid laws considering two infinite planes. In order to check this approach relevance, numerical results are compared to experimental results extracted from the literature (an experiment where water goes through a specimen made of a steel reinforcing bar covered with concrete under load [Desmettre et Charron, 2011] and a device where dry air goes through a structural element made of reinforced concrete [Nahas et al., 2014]). The computation of the flow going to those cracked concrete elements helps to understand the presented approach efficiency and the differences between the two used mechanical models.