L'étude réalisée au cours de cette thèse s'intéresse au fluage secondaire de gainages de Zircaloy-4 détendu dans des conditions représentatives d'un Accident de Perte de Réfrigérant Primaire (APRP) dans le circuit primaire d'un réacteur à eau pressurisée. Elle s'intègre dans le cadre du projet ANR PERFROI. Ce travail s'articule autour de trois axes que sont la conception d'essais de fluage à haute température et en pression interne sur des gaines, leur réalisation et exploitation, la modélisation de ces essais et la caractérisation du comportement au fluage secondaire de cet alliage, et l'étude du fluage secondaire d'éprouvettes de Zr-4 pré-oxydées. Le banc d'essai est conçu pour induire un chargement thermique hétérogène le long des éprouvettes tubulaires testées en utilisant un chauffage par induction. Le fluage à haute température des métaux est très dépendant de la température. Ce gradient thermique couplé à une pression interne mène alors à une déformation hétérogène de l'éprouvette dans l'axe du tube. En réalisant plusieurs chargements en pression interne, plusieurs conditions thermo-mécaniques sont appliquées en un seul essai. Les essais sont instrumentés de manière à mesurer les distributions cinématique et thermique dans une région de 20 mm située entre les spires de l'inducteur. Deux méthodes d'analyse d'images numériques sont adaptées dans ce but et ont été validées par des essais dédiés. Les essais sont modélisés sous le logiciel commercial ABAQUS6.11-2. Un recalage d'une loi de fluage secondaire est réalisé, par des calculs éléments finis intégrés dans un solveur de Gauss-Newton. Cette démarche est validée avec un essai virtuel puis appliquée à 9 essais réalisés sur des échantillons de Zr-4 vierge. Les résultats permettent de déterminer à la fois l'influence de la température et du chargement mécanique de manière fine pour la gamme de température et pression souhaitée. Ils mettent en évidence un changement brutal de comportement à une température de 810 °C, correspondant au début de la transition de la phase alpha vers un matériau biphasé. De plus, les exposants de Norton identifiés dans le domaine alpha pour des contraintes inférieures à 25 MPa sont de l'ordre de l'unité, et entre 4 et 5 pour les contraintes supérieures. Ces ordres de grandeurs sont respectivement associés dans la littérature à des mécanismes de fluage par diffusion et dislocation. Au-dessus de 810 °C, les caractéristiques associées dans la littérature à la superplasticité sont mis en évidence : l'exposant de Norton est compris entre 2,5 et 3,5 et l'énergie d'activation est supérieure à 300 kJ/mol. Des essais sont enfin menés sur des éprouvettes seulement pré-hydrurées, puis pré-oxydées et pré-hydrurées. Les mesures réalisées pendant ces essais sont associées à des métallographies pour mettre en évidence l'effet des couches de zircone et de phase alpha enrichie en oxygène sur la déformation de la gaine à 800 °C. / The present work investigates the steady-state creep behavior of Stress Relieved Annealed Zircaloy-4 claddings under thermo-mechanical conditions simulating a Loss-Of-Coolant Accident (LOCA) in a primary loop of Pressurized Water Reactors. It is part of the PERFROI project focusing on the flow blockage within fuel rods bundle and its potential impact on the cool ability of a nuclear core reactor. The present work includes three main elements that are (a) the design, the performance and the processing of an innovative experiment, (b) the modeling of the tests and the determination of steady-state creep constitutive laws at temperatures ranging from 750 to 850 °C and (c) the preliminary study of the corrosion effect on the Zr-4 claddings creep behavior. An experiment addressing to heterogeneous thermal conditions was designed using an induction heating device. The creep behavior of metals at high temperatures is known to be very sensitive to the temperature. Under an internal pressurization of the cladding this thermal gradient induces its heterogeneous deformation. Using this configuration and three successive internal pressure loadings with a single experiment a lot of thermal mechanical conditions are available. Both thermal and kinematics full field measurements are performed into a 20 mm gauge length located between the induction coils using respectively near infra-red thermography and two dimensional digital image correlation. These two key digital image methods are assessed using dedicated experiments. The tests are modeled using the ABAQUS commercial software. The expected creep behavior law is updated using finite element model calculations included into a Gauss-Newton solver. The method was validated using virtual experiments and then applied to the actual tests. A high consistency between the experiments is obtained and is in agreement with available literature data. Norton exponent and activation energy maps are plotted into the thermal mechanical condition range of interest. These parameters highlight changes of dominant deformation mechanisms in the alpha-phase domain at a 26 MPa von Mises stress and of micro-structure at an 810 °C temperature (at which phase transformation is expected). In the mixed phase domain, the material parameters are similar to those associated in the literature data with super plasticity deformation mechanism. Finally, experiments are performed on pre-oxidized samples at 800 °C. The rupture of the outer zirconia and oxygen enriched alpha layers is observed during the test at a 1.1\% hoop strain.

Ce travail a pour objectif de développer des techniques théoriques et expérimentales pour la détermination des propriétés thermophysiques et terme source. Deux formes de comportement temporel pour le terme source ont été étudiées : un constant et un qui varie dans le temps. La variation dans le temps a été considérée comme une pulse carrée ou une variation sinusoïdale. Deux formes d’échauffement ont été utilisées : une résistance électrique et un laser diode. Pour l’acquisition des données une caméra de thermographie par infrarouge a été utilisée. La stratégie nodale a été utilisée pour contourner le problème des grosses quantités de données générées par la caméra. Le problème direct a été résolu par différences finies, et deux approches pour la solution du problème inverse ont été utilisées, en fonction du comportement temporel du terme source. Les deux approches sont basées sur des méthodes d’inférence statistiques dans une approche Bayésienne, avec la méthode de Monte Carlo via les Chaînes de Markov pour le terme source constant, et le filtre de Kalman pour le problème dont le terme source varie dans le temps. Des manipulations contrôlées ont été faites dans un échantillon avec des propriétés thermophysiques déterminées par des méthodes classiques dans la littérature. / This work deals with the development of new theoretical and experimental techniques for the efficient estimation of thermophysical properties and source-term in micro and macro-scale. Two kinds of source term were studied: a constant and a time varying source term. The time wise variation of the source term had a sinusoidal and a pulse form. Two devices were used for the sample heating: An electrical resistance and a laser diode. For the data acquisition, an infrared camera was used, providing a full cartography of properties of the medium and also non-contact temperature measurements. The direct problem was solved by the finite differences method, and two approaches were used for the solution of the inverse problem, depending on the time varying behavior of the source term. Both approaches deal with the parameters estimation within the Bayesian framework, using the Markov Chain Monte Carlo (MCMC) method via the Metropolis Hastings (MH) algorithm for the constant source term, and the Kalman filter for the time-varying source term. The nodal strategy is presented as a method to deal with the large number of experimental data problems. Experiments were carried out in a sample with well-known thermophysical properties, determined by classical methods.

Ce travail porte sur l'identification de modèles de zones cohésives. Ces modèles, proposés initialement dans les années 60 sont maintenant de intensivement utilisés dans les simulations numériques pour rendre compte de l'initiation et de la propagation de fissures pour différents matériaux et structures.L'identification de ces modèles reste encore aujourd'hui une problématique délicate. Les développements récents de techniques d'imagerie permettent d'accéder à des champs de mesures locales (e.g. déformation et température, …). On se propose dans ce travail d'utiliser la richesse des informations issues de ces techniques d'imagerie pour mettre en place une procédure d'identification qui prenne en compte à la fois le développement de la localisation (effet de structure) mais aussi la nature des différentes irréversibilités mises en jeu (comportement thermo-mécanique). On s'intéresse à des comportements élasto-plastiques endommageables de matériaux ductiles. L'endommagement est associé à un comportement cohésif de l'interface entre les éléments volumiques supposés purement élasto-plastiques.La procédure d'identification comporte deux étapes. La première consiste à caractériser la forme et les paramètres de la loi cohésive sur des essais de traction standard à partir d'une analyse des champs mécaniques localement développés. La seconde étape consiste à vérifier la cohérence thermo-mécanique du modèle identifié en confrontant les mesures calorimétriques déduites des champs de température avec les prévisions du modèle identifié.Cette méthode est appliquée avec succès sur différents matériaux (acier Dual Phase et cuivre). Une attention particulière est portée sur la caractérisation de la longueur caractéristique qui est nécessairement introduite dans l'identification. On montre que cette longueur peut être estimée au regard des différents paramètres introduits dans les traitements d'images.Cette méthode est appliquée sur différents matériaux (acier et cuivre). Une attention particulière est portée sur la caractérisation de la longueur caractéristique qui est nécessairement introduite dans l'identification. On montre que cette longueur peut-être corrélée à l'échelle d'identification des processus d'endommagement sous-jacents. Ainsi les modèles cohésifs identifiés sont fournis au modélisateur avec l'échelle physique à laquelle ils résument l'endommagement volumique du matériau. / This work deals with the identification of cohesive zone models. These models were intially proposed in the 1960s. They are now more and more frequently used in numerical simulations to account for crack initiation and propagation in different materials and structures.The identification of these models still remains a delicate issue. The recent developments in imaging techniques now allow reaching local measurement fields (e.g. strain, temperature,…). We propose here to use the large amount of information given by these techniques to set up an identification procedure accounting for either the localization development (structural effect) and also the character of the different irreversibility sources encountered (thermo-mechanical behavior). We study damageable elasto-plastic ductile materials. Damage is associated to a cohesive behavior of the interface between volumic elements supposed to remain purely elasto-plastic.The identification procedure involves two steps. The first one consists in characterizing the shape and the parameters of the cohesive zone on tensile tests by analyzing the mechanical fields locally developed. The second one consists in checking the thermo-mechanical consistency of the identified model by confronting the calorimetric measurements deduced from temperature fields with the previsions of the identified model.This method is applied on different materials (Dual Phase steel and copper). A specific caution is conferred to the characterization of the characteristic length necessarily introduced by the identification. It is shown that this length can be estimated regarding the different parameters introduced in the image processing.

Cette thèse porte sur le diagnostic de peintures murales masquées par une couche de chaux ou de peintures avec des repeints.Trois méthodes originales sont présentées et le principe validé expérimentalement.D’abord une nouvelle méthode d’excitation de couches profondes par le rayonnement infrarouge lointain λ>20 μm, qui présente une meilleure efficacité par rapport à la méthode classique d’excitation par conduction thermique. Le fait de filtrer les courtes longueurs d’onde élimine l’échauffement en surface du matériau. La semi-transparence de la couche de chaux dans l’infrarouge lointain permet d’éclairer directement les couches de peintures et d’obtenir des thermogrammes révélant les motifs masqués.Nous avons aussi étudié la transmission des couches de peintures dans l’infrarouge lointain par spectrométrie à transformée de Fourier. Les peintures présentent vers 30 μm dans l’infrarouge lointain des zones d’absorption caractéristiques des groupements fonctionnels. Aussi la mesure directe de l’emission totale dans l’infrarouge lointain à température ambiante de peintures permet leur identification. Nous avons conçu un système qui pourrait être déplacé sur le terrain.Nous avons enfin introduit une nouvelle approche pour le diagnostic des peintures murales masquées par une couche de chaux basée sur la mesure du temps de vol des photons balistiques rétrodiffusés collectés par une caméra à balayage de fente d’une résolution de 2 ps. / This study deals with the diagnosis of hidden paintings by a layer of lime or another paintings.Three original methods are presented and experimentally validated the principle. First, a new excitation method of the deeper layers based on far infrared radiation λ>20 μm, which has a higher efficiency compared to the conventional method of thermal conduction excitement. The fact that the filter removes short wavelength in the heating surface of the material. The semi-transparent layer of lime in the far infrared can illuminate layers of paint directly and get thermograms revealing hidden patterns.We also studied the transmission of the layers of paint in the far infrared by Fourier transform spectroscopy. The painting present in the far infrared 30> μm areas obsorption witch caracteristic of functional groups. Also the measurement of total emission in the far infrared at room temperature allows their identification.Finally, we introduced a new approach for diagnosing murals hidden by a layer of lime based on measuring the time of flight of photons backscattered ballistic collected by a streak camera with a resolution of 2 ps.

Ce mémoire aborde la question de la production d’observables intrinsèques au comportement thermomécanique des matériaux pour mieux en formuler les lois d’états. Ces observables sont les sources de chaleur thermomécaniques, activées par sollicitation mécanique. Ces sources peuvent être reconstruites dans l’espace et le temps par inversion de mesures de champs de température obtenus par thermographie IR. Nous présentons essentiellement deux méthodes développées lors de ce travail de thèse qui reposent sur des approches spectrales réduites (dont la décomposition sur Modes de Branche) et des inversions séquentielles (méthode de Beck) ou itératives (Gradient Conjugué). Concernant cette dernière, nous proposons d’y adjoindre une régularisation efficace en s’inspirant de techniques de filtrage par TSVD. S’agissant de matériaux sujets aux instabilités plastiques (PolyEthylène Haute Densité) pour lesquels les vitesses locales peuvent être non négligeables, l’inversion des mesures en température nécessite que l’on considère un opérateur d’advection-diffusion, qui impose alors l’apport d’une connaissance supplémentaire : le champ de vitesses locales. Celui-ci est mesuré par corrélation d’images 3D et nous détaillons le travail expérimental mené ainsi que les résultats obtenus sur des essais de traction pilotés par vidéo-extensométrie. Nous montrons que pour des essais quasi-statiques à vitesses relativement élevées, les effets d’advection sont généralement négligeables. Nous montrons également en quoi la richesse des informations thermomécaniques (Sources) et cinématiques (Taux de déformation, vitesses) permet de mieux comprendre la dynamique de l’instabilité plastique. Enfin nous critiquons les résultats obtenus sur la reconstruction de source par confrontation des deux algorithmes développés et par une analyse physique des phénomènes observés / This work concerns the way intrinsic observables can be produced, which are related to the thermomechanical behavior of materials and necessary for better formulation of state laws. These observables are Thermomechanical Heat Sources (THS) which are activated through mechanical excitation. These sources can be reconstructed both in space and time by the inversion of measured temperature fields obtained through IR thermography. We develop two main methods in this work which rely on spectral reduced approaches (one of them being the decomposition on Branch Modes) and both on a sequential inversion (Beck’s method) and an iterative one (Conjugated Gradient). Regarding the latter, we suggest to combine the standard approach with an efficient regularization method which comes from the filtering techniques based on TSVD. As we are concerned with materials which can be subjected to plastic instabilities (High Density PolyEthylene) for which local velocities of matter displacement can be non negligible, the inversion of the measurements must be performed with the advection-diffusion operator of heat transfer. It is then necessary to obtained additional knowledge: the velocity field. This one is measured by 3D Digital Image Correlation and we detail the experimental work we have carried out, which are based on tensile tests monitored with video-extensometry. We show that for quasi-static tests at relatively high strain rates, the advective effects are generally negligible. We also show the richness of the information brought by this dual thermomechanical (heat sources) and kinematical (strain-rates, velocities) information. It allows for a better understanding of the plastic instability (necking) dynamics. Lastly, we criticize the obtained results on THS reconstruction by the confrontation between the two algorithms and by a physical analysis of the observed phenomena

Ces travaux traitent de l'utilisation de la thermographie infrarouge dans le cadre du contrôle non destructif des ouvrages de génie civil.Une première partie, traite de l’étude des paramètres influençant la mesure in situ, de la capacité de la technique à détecter la variation de propriétés intrinsèques au matériau, et de la détection de délaminations. Les résultats présentés sont issus de mesures expérimentales appuyées par une approche numérique aux éléments finis. Dans un premier temps, une étude sur la capacité de la thermographie infrarouge à détecter des variations de porosité ou de teneur en eau a été réalisée. Dans un second temps, des travaux sont menés sur la détermination des seuils de détection des délaminations en fonction des conditions d’exposition. Ils ont montré un seuil de détection correspondant à un rapport de 2, entre l’extension latérale du défaut et sa profondeur, pour un ensoleillement direct, et à un rapport de 3,3 pour un ouvrage soumis uniquement aux variations de température de l’air. La réalisation d’un suivi temporel combiné à l’étude de l’évolution des gradients temporels de température permet d’améliorer ces seuils de détection. Enfin, une étude originale sur le pontage présent au niveau des délaminations, a montré la prédominance de l’influence de celui‐ci sur la profondeur des délaminations.Une deuxième partie porte sur les travaux réalisés dans le cadre du projet ANR‐SENSO. Ils traitent de la combinaison des résultats issus de différentes techniques de CND dans le cadre plus large de l'amélioration des diagnostics pour la gestion du patrimoine. / This thesis focuses on the use of infrared thermography as a tool for non destructive testing ofbuildings. Mainly, the application is on civil engineering projects.The first part includes identifying the parameters that can affect this in situ technique. Thisparticularly deals with the infrared thermography capacity to detect intrinsic property variations, anddelamination detection. Combination of experiments on concrete slabs and numerical simulationsare used. In a first step, a study on the capacity of thermography to detect porosity and watercontent variation was conducted. In a second step a study on the thresholds for detectingdelaminations based on exposure conditions is carried out. As an outcome, the threshold that hasbeen detected corresponds to a ratio of 2 between the lateral extension of the defect and its depthto direct sunlight; while a ratio of 3,3 if it is exposed to air temperature variations. This studysuggests that a time monitoring combined with the study of the evolution of temporal temperaturegradients can improve the detection limits. Finally, an original study showed the predominance ofthe influence of bridging on the depth of delamination.The second part tackles the works carried out during the ANR project SENSO. Results fromdifferent non destructive tools were coupled for the purpose of improving diagnosis in the assetmanagement.

Les travaux que nous avons menés concernent la modélisation de l’interaction laser-matière couplée à de la thermographie infrarouge, pour la détection de défauts débouchants et sous-jacents dans des matériaux métalliques. A terme, ce procédé de thermographie se veut être une alternative possible à la magnétoscopie et au ressuage dans le domaine des contrôles non destructifs. Ces travaux se sont déroulés en trois phases. Dans un premier temps, nous nous sommes intéressés à la caractérisation des propriétés optiques et thermophysiques des matériaux, pour alimenter le modèle de simulation numérique mais également pour les besoins des mesures de validation du modèle numérique. La deuxième partie de l’étude a permis la mise en place d’un modèle de simulation numérique représentatif de la réalité avec une approche en deux étapes : une première consistant à calculer des champs de températures éléments finis et une seconde à modéliser la caméra de thermographie infrarouge représentatif de celle utilisée dans le projet. Des solutions innovantes ont été utilisées pour obtenir un modèle qui intègre tous les paramètres, avec des temps de calculs faibles. Ce modèle a été ensuite validé par des mesures comparatives des champs de températures numériques et expérimentaux. Enfin, nous nous sommes appuyés sur ce modèle pour rechercher les paramètres optimaux de ce procédé, grâce à la méthode des plans d’expériences numériques. / The aim of this work is the detection of open and subjacent defects in metallic materials using laser-material interaction coupled with infrared thermography. This process is a possible alternative for magnetic particles testing and dye penetrant testing in the field of non-destructive testing. This work is divided into three main parties. At first, we have been interested in the characterization of optical and thermophysical properties of materials we used, in order to have good boundary conditions and also for the needs of temperatures fields measurements for validation. The second part concern the development of a numerical simulation model with two step approach: the first involves modelling temperature field with finite element method and the second details the infrared camera modelling. Innovative solutions were used to obtain model which integrates all parameters, with low time calculation. After validated temperature field, this model allows us to look for the process optimal parameters through numerical experimental designs.

Les critères de localisation et d’endommagement sont généralement basés sur un cadre dissipatif et ce travail s’intéresse aux couplages thermomécaniques accompagnant les micromécanismes de déformation. Il repose en partie sur des données expérimentales obtenues précédemment dans le laboratoire par Bodelot pour un polycristal d’acier A316L. Ce travail tire profit d'une combinaison de techniques différentes, en particulier de mesures in situ de champs cinématiques et thermiques ainsi que de l’Orientation Imaging Microscopy, de la profilométrie et d’une micrographie de surface. Différents outils ont été développés afin (1) d'identifier automatiquement les systèmes de glissement activés, (2) d’estimer l’émissivité de la surface permettant ainsi une détermination des champs thermiques avec une précision de 30 mK, (3) de projeter les champs bruts expérimentaux sur la microstructure et (4) de permettre la modélisation du polycristal et de ses conditions aux limites thermomécaniques réelles dans un cadre de plasticité cristalline dans le code EF Abaqus. Il a notamment été montré que les variations de température fournissent une estimation précise et aisée de la limite d'élasticité macroscopique ainsi que la détermination de la contrainte de cisaillement critique à l'échelle granulaire. En outre, les mesures cinématiques ont permis l'identification des systèmes de glissement activés. Des bilans énergétiques expérimentaux et numériques ont été réalisés et une grande influence de l'hétérogénéité polycristalline sur les mécanismes de stockage d’énergie a été soulignée. Les méthodes proposées contribueront à améliorer les critères d’endommagement basés sur un cadre dissipatif / Strain localization and damage criteria of materials and structures are commonly based on a dissipative framework and this work investigates the thermomechanical couplings accompanying the deformation micromechanisms. It is partly based on experimental data obtained previously in the laboratory by Bodelot for a A316L austenitic stainless steel polycrystal. This work takes profit of a multi-technique approach combining, in particular, in-situ kinematic and thermal fields measurements as well as Orientation Imaging Microscopy, profilometry and surface micrography. Different tools have been developed (1) to automatically identify the activated slip systems directly from the surface micrography, (2) to approach the surface emissivity field allowing an accurate determination of the thermal fields with a 30 mK precision, (3) to project raw experimental fields on the microstructure and (4) to allow the modeling of the polycrystal aggregate and its real thermomechanical boundary conditions by using a crystal plasticity framework within the Abaqus FE code. It has notably been shown that the temperature variations provides an easy and accurate estimation of the macroscopic yield stress at the specimen scale as well as the determination of the Critical Resolved Shear Stress at the intragranular scale. In addition, the local kinematic measurements allow the in-situ identification of the activated slip systems. Experimental and numerical energy balances have been conducted and a great influence of the polycrystalline heterogeneity on the energy storage mechanism has been underlined. The proposed methods would help improving physical based dissipative criteria for damage analysis

L'objectif de cette thèse est de mener une caractérisation expérimentale complète du comportement mécanique et thermomécanique d'un polyuréthane thermoplastique compact et sous forme moussée, utilisé dans la fabrication de pièces antivibratoires dans l'industrie automobile. Cette caractérisation est effectuée en s'appuyant sur des techniques d'imagerie quantitative, la thermographie infrarouge en particulier. D'un point de vue mécanique, la dépendance des différents phénomènes - effet Mullins, hystérésis mécanique, déformation résiduelle et variation de volume - à la densité initiale et à la vitesse de chargement a été établie. Sur le plan thermique et calorimétrique, la signature de ces différents phénomènes a été caractérisée et les bilans énergétiques menés ont permis d'estimer la dissipation mécanique en fonction des conditions d'essais et de la densité initiale du matériau. Les résultats de ces travaux constituent une base de données expérimentale pour l'identification de paramètres constitutifs et l'enrichissement de lois de comportement. / This PhD thesis deals with the full mechanical and thermomechanical characterization of thermoplastic polyurethane, in its compact and foam states, dedicated to antivibratory systems for automotive. This characterization has been carried out by using quantitative imaging techniques, especially infrared thermography. From a mechanical point of view, the dependency of the different phenomena involved in the material deformation - Mullins effect, mechanical hysteresis, residual deformation and volume change - on the initial density and loading rate has been established. The thermal and calorimetric signature of these phenomena has been characterized and the mechanical dissipation was estimated from energy balances. The results obtained provide an experimental database for the identification of constitutive parameters and the enrichment of behavior laws.

L'objectif de la thèse est d'analyser le comportement mécanique de matériaux granulaires composites bidimensionels en terme de textures granulaires en utilisant deux approches : étude expérimentale par "thermoelastic stress analysis" et étude numérique par dynamique moléculaire. Les systèmes granulaires composites sont préparés à l'aide de cylindres en polyoxyméthylène (POM) et polyéthylène haute densité (PEHD), présentant un rapport de rigidité de 4 entre eux. Différents rapports de diamètres et de nombres de particules sont considérés. Les résultats expérimentaux et numériques sont en bon accord à l'échelle macroscopique. En particulier, le réseau fort (qui est ici caractérisé par des contraintes hydrostatiques supérieures à la valeur moyenne) contient moins de 50% des particules, et présente une distribution décroissance exponentielle quel que soit le type de particules considéré pour l'analyse (particules souples, particules rigides, toutes les particules). De plus, la distribution des contacts entre particules rigides (contacts POM-POM) est anisotrope et tend à s'organiser dans le sens de la direction du chargement extérieur appliqué, tandis que les autres types de contact agissent principalement pour maintenir le système en équilibre. / The main objective of our dissertation is to analyze the mechanical behavior of two-dimensional composite granular materials through the granular textures. Thermoelatic stress analysis experiments and molecular dynamics simulations are used for this purpose. The composite granular systems are prepared from polyoxymethylene (POM) and high-density polyethylene (HDPE) cylinders with a stiffness ratio of about 4 between them. Different configurations in terms of ratios of diameter size and ratio of particle numbers are systematically investigated. Experimental and numerical results are good correlated at the macroscopic scale. In particular the strong network, which is here characterized by hydrostatic stresses higher than the mean value, consists of less than 50% of all particles, and exhibits an exponential decay whatever the type of particles considered for the analysis (soft, stiff, or both types). In addition, the contact distributions between stiff particles (POM-POM contacts) is anisotropic with an effort to arrange parallel to the direction of the external applied load, whereas the other types of contacts just act to sustain the granular system in equilibrium.