Ce travail présente une caractérisation expérimentale de la localisation de la déformation et la rupture sur un matériau synthétique analogue de la roche dans différentes conditions de chargement axisymétrique. L’évolution de la localisation de la déformation a été caractérisée par la technique de corrélation d’images bidimensionnelle. Deux gammes d’essais en compression et en extension ont été réalisées à différentes pressions de confinement Pc. Profitant de la technique DIC et de la caractérisation détaillée des propriétés constitutives de GRAM1 par Mas et Chemenda (2014; 2015), nous montrons que l'initiation de la localisation des bandes de déformation est précédée par la déformation élastique uniforme puis la déformation élasto-plastique. La localisation de la déformation est initiée dans le régime de dilatance (la dilatance β est positive) et d’adoucissement (le module d’écrouissage h est négatif). Au cours de l'évolution du chargement, la déformation à l'intérieur de la bande de localisation de la déformation devient compactive (β <0) et est accompagnée par un durcissement du matériau (h> 0). Ceci provoque dans un premier temps l’élargissement de la bande puis l’apparition de nouvelles bandes qui se forment successivement. La formation de ces nouvelles bandes provoque un ralentissement de la déformation dans les bandes préexistantes, ce qui densifie progressivement le réseau. Dans les essais en extension axisymétrique, les fractures se forment perpendiculaires à la contrainte principale mineure σ3 sous un état de contrainte extensif à faible Pc puis compressif quand Pc augmente. Nos résultats présentent une transition continue de la fracture extensive à la fracture cisaillante avec une augmentation de la contrainte compressive. Ceci est en accord avec les résultats obtenus pour les essais en extension réalisées sur des roches naturelles (Ramsey et Chester, 2004 ; Bobich 2005) / This work presents an experimental characterization of the strain localization and the rupture on a granular rock analogue material under different conditions of axisymmetric loading. The evolution of the strain localization was characterized by the two-dimensional image correlation technique. Two series of triaxial compression and extension tests were carried out at different confining pressure Pc. Taking advantage of the DIC technique and detailed characterization of GRAM1’s constitutive properties by Mas and Chemenda (2014, 2015), we show that the initiation of deformation localization bands is preceded by the uniform elastic and then elasto-plastic deformation. The strain localization is initiated in the dilatant regime (positive dilatancy factor β) and strain-softening (negative hardening modulus h). During the band evolution, the deformation within it becomes compactive (β<0) and is accompanied by the material hardening (h>0), which causes the band to widen and new bands to form successively. The formation of new bands causes a slowing down the deformation within the prior bands, which makes the progressively densified band network to continuously evolve. In axisymmetric extension tests, the fractures are formed perpendicular to the minor principal stress σ3 in an extensive stress state at low Pc then compressive when Pc increases. Our results show a continuous transition from the extensive fracture to shear fracture with an increase of compressive stress. This is suitable of the results obtained for extension tests performed on natural rocks (Ramsey and Chester, 2004; Bobich, 2005)

Les risques liés à la propagation de fissures sous impact sont encore très difficiles à estimer. La détermination de critères de rupture dynamique uniquement à partir de résultats expérimentaux reste délicate. Ainsi la première étape pour valider des lois de propagation de fissures sous impact passe par le développement d'outils de simulation numérique. Depuis les années 1970, de nombreux codes de calcul mécanique ont été dédiés à l'étude de la propagation de fissures, notamment dans le cas du phénomène de fatigue. La principale difficulté consiste dans la nécessité de suivre la géométrie de la fissure au cours du temps. Ces dernières années, des méthodes alternatives basées sur la partition de l'unité ont permis une description implicite des discontinuités mobiles. C'est le cas de la méthode des éléments finis étendue (X-FEM) qui paraît particulièrement adaptée à la simulation de la propagation dynamique de fissures sous chargement mixte où les trajets de fissures ne sont pas connus a priori. Si ces outils numériques permettent maintenant de représenter l'avancée dynamique d'une fissure, les résultats numériques doivent être comparés à des résultats expérimentaux pour s'assurer que les lois introduites sont physiquement fondées. Notre objectif est donc de développer conjointement des techniques expérimentales fiables et un outil de simulation numérique robuste pour l'étude des phénomènes hautement transitoires que sont l'initiation, la propagation, l'arrêt et le redémarrage de fissures sous impact.<br />Des expériences de rupture dynamique ont donc été réalisées sur du Polyméthacrylate de méthyle (PMMA) durant lesquelles la mixité du chargement varie et des arrêts et redémarrages de fissures se produisent. Deux bancs d'essais différents ont été utilisé, le premier basé sur la technique des barres de Hopkinson (ou barres de Kolsky), le second mettant en jeu un vérin rapide. Le PMMA étant transparent, la position de la fissure au cours de l'essai a été acquise grâce à des caméras rapides mais aussi en utilisant un extensomètre optique (Zimmer), habituellement dédié à la mesure de déplacements macroscopiques d'un contraste noir/blanc. L'utilisation de cet extensomètre pour suivre la fissure au cours de l'essai a permis d'obtenir une localisation très précise de la pointe de la fissure en continu, permettant ainsi l'étude des phases transitoires de propagation. Afin d'étudier le même phénomène dans des matériaux opaques comme les aluminiums aéronautiques (Al 7075), des techniques de corrélation d'images numériques ont été employées en mouchetant les éprouvettes impactées. De nouveaux algorithmes ont été développés afin de traiter les images issues d'une caméra ultra-rapide (jusqu'à 400 000 images par seconde).<br />Plusieurs géométries ont été envisagées afin d'étudier différents cas de propagation dynamique : initiation en mode I pur, initiation en mode mixte, propagation, arrêt, redémarrage, interaction entre deux fissures, influence d'un trou sur le trajet d'une fissure, branchement dynamique de fissures. Ces expériences ont ensuite été reproduites numériquement afin de valider les algorithmes et les critères de rupture choisis.

Le développement de méthodes d’identification de paramètres de lois de comportement de matériaux est devenu primordial pour avoir accès à la connaissance complète du comportement. En effet, les méthodes de mesure optiques, comme la Corrélation d’Images Numériques, permettent d’obtenir les quantités cinématiques de la relation de comportement sous forme de champs de vecteurs. En revanche, les contraintes ne sont généralement pas mesurables et il est nécessaire d’identifier les paramètres de la loi de comportement du matériau considéré pour y avoir accès. Plusieurs méthodes ont vu le jour et permettent de répondre à cette problématique mais la plupart d’entre elles supposent une homogénéité du matériau. Ce mémoire traite de l’application de certaines de ces méthodes, notamment la méthode de l’écart à l’équilibre (MEQ) et la méthode de recalage de modèle éléments finis (MREF), dans des matériaux hétérogènes à microstructure complexe où les propriétés mécaniques évoluent spatialement dans le volume. L’objectif est d’identifier ces propriétés mécaniques locales qui régissent la cinématique mesurée de tels matériaux dans le cadre de l’élasticité linéaire isotrope. Dans un premier temps, les deux méthodes citées sont décrites, implémentées et comparées sur des cas simulés en 2D. La MREF est préférée à la MEQ car plus robuste vis-à-vis des incertitudes de mesure. Basée sur un formalisme itératif, une parallélisation de l’algorithme a été opérée pour diminuer le coût en temps de la méthode. Des expérimentations dans le plan sur des éprouvettes en polyuréthane où les hétérogénéités sont maîtrisées ont permis de valider la méthode. Enfin, deux applications en 3D sur un matériau en mousse polyuréthane et un composite à base de fibres de bois démontrent l’intérêt d’une telle méthode pour l’identification de propriétés mécaniques locales. La mise en évidence d’une relation entre les propriétés locales identifiées et les propriétés locales de la microstructure de ces matériaux est effectuée. / The development of identification methods of material constitutive equation parameters has become fundamental to completely know the mechanical behavior. Indeed, optical methods, such as Digital Image Correlation, allows to get kinematics quantities of the constitutive equation as vectors fields. But, stresses are usually not available experimentally and one has to identify constitutive equation parameters to compute them. Several methods have been developed and answer to that problematic but most of them suppose the materials as homogeneous. This memoir is about the application of some of these methods, such as the equilibrium gap method (EGM) and the finite element model updating method (FEMU), in the case of heterogeneous materials with complex structures where mechanical properties vary spatially in the volume. The objective is to identify these local mechanical properties which rule the measured kinematics of such materials considering the isotropic linear elasticity. Firstly, both methods are detailed, implemented and compared on 2D simulated cases. The FEMU method is preferred because it is more robust in the presence of noisy data. Based on an iterative process, a parallelisation of the algorithm is achieved in order to reduce the cost of the method. In-plane experiments on polyurethane samples where heterogeneities are controlled have validated the method. Finally, two 3D applications on a polyurethane foam material and a wood-based fibrous composite have demonstrated the interest of this approach to identify local mechanical properties. The highlighting of a relationship between identified local properties and microstructural properties of these materials is made.

Les outillages de mise en forme en acier martensitique de type AISI H11 sont des pièces critiques dont le comportement en service est étroitement lié à leurs structures internes et à leur évolution. Les conditions des sollicitations lors de la mise en oeuvre du procédé est souvent à l'origine de modifications microstructurales en surface, à savoir la morphologie des lattes de martensite, les orientations cristallographiques, l'état d'écrouissage interne ou encore le profil de surface. Ces aspects peuvent éventuellement altérer les performances mécaniques de l'acier AISI H11. Afin d'appréhender et d'optimiser le comportement mécanique de celui-ci, une approche multi-échelle est mise en oeuvre dans ce travail. Celle-ci s'articule autour d'une investigation expérimentale et d'un traitement numérique. L'étude expérimentale s'attache à reproduire, à l'échelle du laboratoire, des surfaces équivalentes à celles issues lors des procédés de mise en oeuvre des outillages. Des techniques de caractérisation spécifiques, à savoir le MEB, l'EBSD, la nanoindentation ou encore l'altimétrie permettent de mettre en évidence un gradient de la stéréologie du matériau en surface et sous-surface. Les hétérogénéités locales induites concernent la morphologie des lattes de martensite, les orientations cristallographiques, l'état d'écrouissage interne mais également le profil de surface. Des essais mécaniques in-situ associés à la technique de corrélation d'images numériques sont réalisés pour des chargements monotones quasi-statiques et cycliques de type traction-traction. Une investigation des champs mécaniques locaux en surface est ainsi effectuée, elle permet d'analyser les schémas de localisations des déformations non linéaires liés aux artéfacts stéréologiques. Le traitement numérique s'intéresse à une modélisation multi-échelle, et plus particulièrement à des calculs par la méthode des éléments finis sur des microstructures virtuelles générées par tesselations de Voronoï. Celles-ci sont effectuées de manière à reproduire les structures martensitiques et considèrent des relations d'orientations spécifiques (de type Kurdjumov-Sachs) à l'issue du traitement thermique entre les lattes de martensite et le grain austénitique parent. Les équations constitutives du modèle de plasticité cristalline (élasto-viscoplastique) de Méric-Cailletaud sont implantées dans le code de calcul par éléments finis Abaqus dans le cadre de l'hypothèse des petites perturbations (HPP) et de la théorie des transformations finies. La formulation du modèle dans le contexte de la théorie des transformations finies est effectuée dans le cadre d'une description spatiale où la notion de dérivée objective est considérée. Celle-ci consiste en celle d'Oldroyd ou de Truesdell de manière à ce qu'une telle formulation soit équivalente à une description lagrangienne. Le traitement numérique a permis de reproduire de manière qualitative les schémas de localisation en surface mise en évidence lors de l'investigation expérimentale. L'influence des divers paramètres stéréologiques, évoqués ci-dessus, sur les champs mécaniques locaux a été analysée. De par cette approche, il a été possible de mettre en évidence certains mécanismes élémentaires, notamment les effets d'interaction et de surface. Enfin, il a été constaté que la prise en compte des rotations des réseaux cristallins par la théorie des transformations finies permet de relâcher certaines zones de localisation des champs mécaniques autour d'artéfacts stéréologiques. / AISI H11 martensitic tool steels are critical mechanical components that behaviour during service is drastically linked to their internal structures and their possible evolution. Their manufacture processes are often at the origin of microstructural changes at the surface, namely the morphology of martensitic laths, the crystallographic orientations, the internal hardening state and the surface profile These aspects can potentially alter the mechanical performance of AISI H11 martensitic steel. In order to get better insight into and optimize its mechanical behaviour, a multi-scale approach involving an experimental investigation and a numerical treatment is taken in this work.The experimental investigation focuses to reproduce, at the laboratory scale, equivalent surfaces to those resulting from tool steels manufacture processes. Specific characterization techniques, namely SEM, EBSD, nanoindentation and altimetry enable to highlight a stereology gradient of the material in surface and sub-surface. The induced local heterogeneities consist in morphology of martensitic laths and crystallographic orientations, internal hardening state and surface profile. In-situ mechanical tests with digital image correlation technique (DIC) are carried out for monotonous quasi-static and tension-tension cyclic loads. An investigation of the local mechanical fields at the surface is thus performed and allows to analyze the localizations schemes of nonlinear strains which are related to stereological artifacts.The numerical treatment is focused on a multi-scale modelling, and more particularly on finite element calculations on virtual microstructures which are generated by Voronoi tesselations. The latters are carried out such that to reproduce martensitic structures and consider a specific orientation relationship between martensitic laths and parent austenitic grains (i.e. Kurdjumov-Sachs) after the heat treatment. The constitutive equations of the (elasto-viscoplastic) crystal plasticity of Méric-Cailletaud are implemented in the finite element code Abaqus in the context of the small strain assumption and the finite strain theory. The formulation of the model in the context of finite strain theory is is given a spatial description where the notion of objective derivative, namely the so called one of Oldroyd or Truesdell, is used in such a way that such formulation is equivalent to a Lagrangian description.The numerical treatment has allowed to qualitatively reproduce the localization patterns at the surface which have been highlighted in the experimental investigation. The influence of the different stereological parameters mentioned above on the local mechanical fields was analyzed. By this approach, it was possible to highlight some elementary mechanisms including interaction and surface effects. Finally, it was found that the inclusion of lattice rotations via the theory of finite strain allows to release certain areas of mechanical fields localization that are related to stereological artifacts.

Ce travail s’inscrit dans le contexte d’une meilleure prédiction de la durée de vie de structures soumises au risque d’amorçage et de propagation de fissures multiples de corrosion sous contrainte (CSC). Ainsi, en développant une méthodologie expérimentale originale basée sur des mesures conjointes de corrélation d’images numériques (DIC), d’émission acoustique (EA) et de bruit électrochimique (EN), et leur analyse, ce travail vise à contribuer à l’amélioration de la compréhension des mécanismes impliqués dans le développement de colonies de fissures courtes de CSC qui interagissent entre elles, et à la simulation du comportement de cette colonie. Le choix de conditions optimales de traitement thermique d’un alliage base-Ni et de pH d’une solution de polythionate, a permis la maitrise des paramètres géométriques et morphologiques de colonies de fissures intergranulaires, qui ont été investiguées par DIC grâce à l’optimisation d’une préparation de surface adaptée. Les différentes phases de propagation d’une fissure unique et d’une colonie de fissures ont été identifiées, de même que les mécanismes associés, par des expérimentations et analyses réalisées en 2D et en 3D. Cette approche expérimentale innovante a permis de poser les bases d’une approche numérique puis de la valider. Un focus particulier a été porté sur l’EN au travers d’une analyse critique des perturbations engendrées par le bruit de l’instrumentation et par l’asymétrie du système d’étude. Les limitations de la technique pour son application à l’étude quantitative de la corrosion sous contrainte ont été évaluées sur la base des résultats de l’étude. Une transposition de la démarche expérimentale à hautes pressions et températures est proposée comme perspective à court terme de ce travail, qui permet également d’envisager la prise en compte de différents modes de propagation des fissures en lien avec la microstructure du matériau dans l’approche numérique. / This work concerns with the current needs of enhancing the tools used for predicting the remnant lifetime of structures subjected to the risk of initiation and propagation of multiple stress corrosion cracks (SCC). The approach consists in developing an original experimental methodology based on joint measurements of digital image correlation (DIC), acoustic emission (EA) and electrochemical noise (EN). The final objective is to contribute to both the understanding of the mechanisms involved in the development of interacting short stress corrosion cracks and to the modeling of the colony behavior. The choice of optimal conditions for the heat treatment of a Nickel-base alloy and for the pH of a polythionate solution allowed controlling the morphological parameters of intergranular cracks colonies, which were investigated by DIC owning to an optimized suitable surface treatment. The different propagation stages of a single crack and some colonies were identified, together with the involved mechanisms, through experiments and analyses performed in 2D and 3D. This innovative experimental approach allowed settling the basements of the numerical modeling and validating it. A particular attention was focused on EN measurements through a critical analysis of the perturbations generated by the instrumental noise and the asymmetry of the studied system. The limitations of the technique for its application to the quantitative study of SCC were evaluated on the basis of the present results. A transposal of the experimental approach towards high temperature and pressure conditions of test was finally proposed as a short-term prospect of this work, also allowing considering other modes of crack propagation linked to the material microstructure in the numerical approach.

Cette thèse porte sur l'étude de milieux granulaires constitués de particules molles. Elle s'appuie d'une part sur la méthode de corrélation d'images numériques (DIC) et d'autre part sur des simulations associant la méthode des éléments finis (FEM) et la méthode de Dynamique des Contacts (CD) permettant de tenir compte de la déformabilité élevée des particules et des interactions entre particules. Pour une large gamme de matériaux complexes (élastique, plastique et mousse solide), l'étude préliminaire de particules uniques comprimées radicalement valide la technique de DIC pour les mesures des champs de déformation, et pour la détermination de l'énergie de déformation pour des matériaux élastiques. Des grandeurs macroscopiques et micro-structurales tels que la compacité, la coordinence, l'énergie élastique, la densité de la probabilité de force de contact ainsi que de la densité d'énergie sont mesurées, et comparées entre les résultats numériques et expérimentaux en particulier dans le régime où les déformations des particules sont d'amplitudes finies. L'accord quantitatif entre ces approches permet de valider à la fois la méthode expérimentale pour l'étude de milieux granulaires à particules déformables, et d'autre part l'approche numérique. Forts de ces résultats, nous avons mené une étude numérique de la compression uniaxiale d'un assemblage de cylindres Néo-Hookéens. L'effet du frottement sur les paramètres macroscopiques et mico-structuraux est déterminé, y compris lorsque la compacité du système est proche de 100% / This thesis deals with the study of granular media composed of soft particles. It relies on the one hand on the method of Digital Image Correlation(DIC) and on the other hand, on simulations coupling the Finite Element Method and the Contact Dynamics method, taking into account the high deformability of particles and interactions between particles. For a wide range of complex materials (elastic, plastic and solid foam), the preliminary study of radially compressed single particles validates the DIC technique for strain field measurements, and for the determination of the strain energy for elastic materials. Macroscopic and micro-structural quantities such as packing fraction, coordination, elastic energy, probability density of the contact force as well as strain energy density, are measured and compared between the numerical and experimental results in particular in the regime where the deformations of the particles are large. The quantitative agreement between these approaches makes it possible to validate both the experimental method for the study of granular media with deformable particles, and the numerical approach. Based on these results, we conducted a numerical study of the uniaxial compression of a Neo-Hookean cylinder assembly. The effect of friction on the macroscopic and mico-structural parameters is determined, even when the packing fraction of the system is close to 100%.

Le comportement de matériaux subissant des chargements dynamiques à hautes énergies (impacts, explosions \dots) est généralement sensible à la vitesse de déformation (viscoplastique) et/ou à l'endommagement. Les procédures classiques de caractérisation des modèles de comportement associés utilisent des essais statiquement déterminés nécessitant le respect d'hypothèses limitatrices. Ainsi, le traitement de champs mécaniques hétérogènes est impossible, ce qui limite l'exploitation des essais à de faibles gammes de déformations et oblige à réaliser plusieurs tests à vitesse de déformation constante pour la caractérisation de la viscoplasticité. Enfin, ces restrictions empêchent de profiter pleinement de la quantité importante d'information accessible par mesures de champs. Une solution est alors de recourir à des essais statiquement indéterminés, permettant d'exploiter des champs hétérogènes. Parmi les outils disponibles, la Méthode des Champs Virtuels (MCV) présente des avantages certains par rapport aux méthodes classiques de recalage par éléments finis. L'accent est mis dans ces travaux sur le développement de la MCV pour la caractérisation du modèle viscoplastique de Johnson-Cook. Un des grands atouts de la MCV est qu'elle permet de caractériser la partie viscoplastique du modèle à partir d'un unique essai en conditions de chargement dynamiques, grâce à son exploitation statiquement indéterminée de champs de déformations et de vitesses de déformations hétérogènes. Une perspective à court terme, dont la faisabilité est vérifiée sur des tests numériques, est l'identification simultanée par la MCV de paramètres de modèles élastoplastiques endommageables (Lemaitre).

La compréhension du comportement musculaire est un champ d'investigation primordial dans divers domaines tels que le sport, le crash automobile et la médecine. Peu de techniques in vivo permettent aujourd'hui de caractériser les propriétés mécaniques du muscle. Aussi, l'objectif de ce travail est de développer une méthode d'identification in vivo à partir demesures de champs de déplacement. La méthode d'identification est composée de 3 étapes interdépendantes. La première consiste à réaliser une IRM de la cuisse pour segmenter manuellement différents tissus musculaires (quadriceps, ischios, gracilis et sartorius) et le tissu adipeux. Un comportement Néo-Hookéen est choisi pour modéliser le comportement hyperélastique (C10, D). Dans un second temps, un dispositif expérimental de compression est développé pour mesurer le champ de déplacement in vivo à partir des techniques d'imagerie échographique et de Corrélation d'Images Numériques. Finalement, une méthode inverse est mise en œuvre pour identifier les paramètres C10 et D de chaque tissu. Par le biais d'un exemple numérique, l'erreur sur les paramètres identifiés est évaluée. Les cartographies des champs de déplacement expérimentaux confirment les observations qualitatives obtenues sur les images échographiques et sont validées par les champs de déformation obtenus par approximation diffuse. Une faible erreur d'identification (C10<3%, D<7%) est obtenue à partir de l'exemple numérique, et les paramètres mécaniques identifiés sont en accord avec la littérature. Ces résultats valident la démarche inverse mise en œuvre qui permettra, à terme, de suivre l'évolution des pathologies et de mener des simulations prédictives.

L'objectif dans cette recherche expérimentale est de contribuer à la compréhension des phénomènes de l'endommagement et de la fissuration dans les sols argileux saturés sous chargement hydrique (dessiccation) et mécanique (chemins de contrainte). L'étude mécanique est axée sur l'identification des propriétés élastiques du matériau endommagé sous chargement triaxial à différents niveaux et chemins de contrainte. Les principaux points traités étaient, l'influence de la contrainte de consolidation (p'ic), du degré de surconsolidation (OCR) et du chemin de contrainte sur la vitesse d'onde de compression. Les résultats montrent que les propriétés élastiques sont affectées et évoluent avec la déformation. La contractance semble augmenter la vitesse d'onde de compression à mesure que p'ic croît et densifie le matériau. En revanche, la dilatance semble plutôt diminuer la vitesse d'onde de compression. Par ailleurs, différents chemins triaxiaux conduits à un niveau de contrainte donné ont montré que les vitesses d'onde sont très proches si les chemins sont normalement consolidés. Les vitesses sont proches également si les chemins sont surconsolidés. En revanche, à un même état de contrainte, la vitesse d'onde est différente si on compare un chemin normalement consolidé avec un chemin surconsolidé. La deuxième partie de l'étude porte sur une analyse approfondie des mécanismes d'amorçage et de propagation des fissures sous dessiccation libre et contrôlée d'une argile initialement saturée. L'analyse via la corrélation d'images numériques (CIN) a permis, en premier lieu, une caractérisation multiéchelle du phénomène de retrait opéré au cours du séchage. Différents phénomènes ont été identifiés : le phénomène de relaxation, caractérisé par un changement de signe de la déformation locale extension / compression, le mécanisme de fissuration par extension (mode I), identifié par la cartographie des directions des déformations principales. L'étude a été complétée par des essais de traction indirecte par flexion 3 points réalisés sur des poutrelles de kaolin K13 initialement soumises à différentes succions. L'effet de la succion initiale sur la résistance à la traction a donc été observé. Les résultats montrent que plus la succion initiale est élevée, plus la pente de la courbe de variations de la résistance à la traction est élevée. D'autre part, la déformation de la zone tendue de la poutrelle atteint sa valeur critique lorsque la contrainte atteint la résistance maximale à la traction de l'argile / The objective of this experimental research is to contribute to understanding phenomena of damage and cracking in saturated clay soils under hydric loading (desiccation) and mechanical loading (stress paths). Mechanical study focuses on identifying the elastic properties of the damaged material under triaxial loading at different levels and stress paths. The main issues were the influence of; the consolidation stress (p'ic), the overconsolidation ratio (OCR) and the stress path on the compression wave velocity. The results show that the elastic properties are affected and evolve with strain. Contractancy appears to increase the compression wave velocity as p'ic increase and densifies the material. However, dilatancy seems rather decrease the compression wave velocity. Moreover, various triaxial paths performed at a given stress level showed that the wave velocities are very close if the paths are normally consolidated. Velocities are also close if the paths are overconsolidated. However, at the same stress state, the wave velocity is different when comparing a normally consolidated path with an overconsolidated path. The second part of the study focuses on a thorough analysis of the cracks' initiation and propagation mechanisms under free and controlled desiccation on initial saturated clay. Analysis via the digital image correlation (DIC) allowed, first, a multi-scales characterization of the shrinkage phenomenon operated during drying. Different phenomena were identified; the relaxation phenomenon characterized by a change of the local strain's sign extension / compression, the cracking mechanism by extension (mode I) identified by the map of the directions of the principal strains. The study was completed with indirect tensile tests by bending performed on kaolin k13 beams initially submitted to different suctions. The effect of the initial suction on the tensile strength has been observed. The results show that more the initial suction is high, the greater the slope of the curve of changes in the tensile strength is high. On the other hand, strain of the tension zone of the beam reaches its critical value when the stress reaches the maximum tensile strength of the clay

Le travail réalisé durant cette thèse a été motivé par l'étude des mécanismes microscopiques à l'origine du fluage dans les matériaux granulaires.%En particulier, on cherche à explorer des techniques expérimentales et numériques pour l'étude d'un tel phénomène.Dans une première partie, on cherche à mesurer les déplacements des grains dans un matériau granulaire par observations en micro-tomographie X. Une telle identification ne peut être efficacement réalisée pour des phénomènes rapides avec les méthodes classiques de corrélation d'images numériques. Une nouvelle méthode nommée emph{corrélation discrète des projections numériques} qui contourne cette difficulté est développée dans cette thèse. Cette méthode, basée sur la corrélation des projections de tomographie, permet de mesurer les déplacements avec un nombre réduit de projections (100 fois moins que les méthodes classiques), ce qui diminue énormément le temps d'acquisition nécessaire pour la mesure. La méthode, appliquée à des données expérimentales, donne une précision comparable à celles des méthodes classiques tandis que le temps d'acquisition nécessaire est réduit à quelques minutes. Une étude portant sur l'analyse des sources d'erreurs affectant la précision des résultats est également présentée.Le but de la deuxième partie est de réaliser des simulations numériques pour fournir une caractérisation de l'essai oedométrique. Différents assemblages de billes de verre légèrement poly-disperses interagissant à travers des contacts élastiques de Hertz-Mindlin et frottement de Coulomb ont été utilisés. Ces simulations ont permis d'étudier l'évolution de certains paramètres structuraux du matériau modèle, préparant ainsi le terrain pour de futures études sur le fluage. Il a été particulièrement souligné que les contacts élastiques utilisés dans ces simulations ne reproduisent pas l'irréversibilité des déformations observée dans les expériences sur des sables. Cependant, l'irréversibilité est bien visible sur le nombre de coordination et l'anisotropie. Alors que les paramètres élastiques peuvent exprimer la réponse pour des petits incréments de déformations, la compression oedometrique est belle et bien anélastique, principalement à cause de la mobilisation du frottement. Le rapport entre les contraintes horizontales et verticales (coefficient du sol au repos) n'est particulièrement constant que lorsque l'anisotropie de structure est instaurée dans l'état initial de l'assemblage. Il est par ailleurs relié à l'anisotropie interne de la structure par une formule simple. Finalement, les coefficients du tenseur élastique dépendent principalement du nombre de coordination et son anisotropie est plus liée à l'anisotropie des contacts qu'à celle des forces / The present work is motivated by the study of creep in granular materials at the microscopic scale.The first part of this thesis deals with displacement measurements by microtomography. Classical digital image correlation fails to catch time-dependent (possibly fast) phenomena such as short-term creep. A new method named emph{Discrete Digital Projection Correlation} is developed to overcome this limitation. This method requires very few projections (about 100 times less than classical methods) of the deformed state to perform the correlation and retrieve grain displacements. Therefore, the acquisition time is remarkably reduced, which allows to study time-dependent phenomena.The method is tested on experimental data. While its accuracy compares favorably to that of conventional methods, it only requires acquisition times of a few minutes. The origins of measurement errors are tracked by numerical means, on simulated grain displacements and rotations.The second part is a numerical simulation study, by the Discrete Element Method (DEM), of oedometric compression in model granular materials, carried out with a simple model material: assemblies of slightly polydisperse spherical beads interacting by Hertz-Mindlin contact elasticity and Coulomb friction. A wide variety of initial states are subject to compression, differing in density, coordination number and fabric anisotropy. Despite apparently almost reversible strains, oedometric compression proves an essentially anelastic and irreversible process,due to friction, with important internal state changes affecting coordination number and anisotropy. Elastic moduli only describe the response to very small stress increments about well equilibrated configurations. The ratio of horizontal stress to vertical stress (or coefficient of earth pressure at rest, commonly investigated in soil mechanics) only remains constant for initially anisotropic assemblies. A simple formula relates it to force and fabric anisotropy parameters, while elastic moduli are mainly sensitive to the latter. Further studies of contact network instabilities and rearrangements should pave the way to numerical investigations of creep behavior