Cette thèse porte sur l'étude de l'écoulement et du transfert thermique de polymère dans les deux procédés de rotomoulage: rotomoulage réactif et rotomoulage de poudres thermoplastiques. L'écoulement du polymère réactif au cours du rotomoulage a été modélisé par l'écoulement d'un fluide réactif à surface libre sous gravité à l'intérieur d'un cylindre horizontal en rotation autour de son axe principal. Dans ce cas, la méthode lagrangienne ''Smoothed Particle Hydrodynamics (SPH)'' a été appliquée en utilisant une expression adaptée de la viscosité variable du polymère réactif. À l'état de poudre, l'écoulement et le transfert thermique de polymère au cours du rotomoulage des poudres thermoplastiques a été modélisés en utilisant un code de calcul basé sur la ''Méthode des Éléments Discrets (DEM)''. Les régimes d'écoulement, la ségrégation, le mélange et le transfert thermique des poudres dans un moule en rotation ont été simulés. Chaque fois que cela était possible des résultats des modèles ont été validés expérimentalement avec succès. / This Ph.D research work is about the study of polymer flow and heat transfer in both processes of rotational molding: reactive rotational molding and rotational molding of thermoplastic powders. The flow of the reactive polymer in the reactive rotational molding was modeled by a free surface flow under gravity inside a horizontal rotating cylinder. The meshless and lagrangian method ''Smoothed Particle Hydrodynamics (SPH)'' was applied to simulate the flow by using an adapted variable viscosity expression. In powder state, the flow and heat transfer of polymer during rotational molding of thermoplastic powders are modeled using a solver based on the ''Discrete Element Method (DEM)''. In which, the flow regime, the segregation, the mixture and the heat transfer of polymer powders inside a horizontal rotating cylinder was simulated. Results of the model have been validated experimentally with success.

L'agglomération des fines entraîne une augmentation importante de la demande en eau des mélanges cimentaires et constitue un frein à la valorisation des sous-produits contenant des proportions élevées de fines. Une modélisation discrète de la microstructure des agglomérats de fines particules est développée. Le modèle intègre les forces de contact avec prise en compte du frottement et du roulement, les forces physico-chimiques et les forces hydrodynamiques calculées avec l'approximation de drainage libre. Le processus d'agglomération, de fragmentation et de restructuration sous cisaillement d'agglomérats fractals de fines particules minérales est étudié en vue de quantifier l'eau immobilisée dans les flocs. Un agglomérat initial isolé et lâche est soumis à un taux de cisaillement constant. Lorsque les forces hydrodynamiques sont prépondérantes par rapport aux forces de cohésion, l'agglomérat mère se fragmente en agglomérats secondaires jusqu'à atteindre un état stationnaire, où le nombre, la taille et la structure des flocs peuvent être considérés comme constants. L'influence du gradient de vitesse, du contexte physico-chimique, de la taille et de la polydispersité des particules sur les caractéristiques microstructurales des flocs à l'équilibre est étudiée. Il est montré que le comportement d'une suspension diluée monodisperse peut être décrit par un unique nombre adimensionnel faisant intervenir la force de cohésion maximale, la taille des particules et le gradient de vitesse. / The agglomeration of fines leads to a significant increase in water demand of cement mixtures and hinders the valorization of by-products containing high fines proportions. A discrete modeling of the microstructure of fine particles agglomerates is developed. The model takes into account contact forces, physico-chemical, and hydrodynamic forces computed with the free draining approximation. The agglomeration, fragmentation and restructuring process of fractal agglomerates of fine mineral particles under shear is studied to quantify the entrapped water in the flocs. A loose initial agglomerate isolated is submitted to a constant shear rate. When the hydrodynamic forces predominate over the adhesive forces, the parent agglomerate is broken up into secondary agglomerates until reaching a steady state, where the number, size and structure of flocs can be considered constant. The influence of the velocity gradient, the physico-chemical context, particles size and their polydispersity on the microstructural characteristics of flocs at steady state is studied. It is shown that the behavior of a monodisperse dilute suspensions can be described by a dimensionless single number involving maximum cohesion force, particle size and velocity gradient.

La technique de changement d’échelles a été largement développée dans la littérature pour décrire le comportement global des milieux granulaires en prenant en compte leurs propriétés locales. Cette technique considère classiquement deux échelles : l’échelle macroscopique du volume élémentaire représentatif et l’échelle microscopique du contact entre particules. Le défi majeur de ce changement d’échelles "micro-macro" réside dans la définition de la déformation macroscopique : en effet, si la contrainte macroscopique peut être clairement définie à partir des forces de contact, il a été montré qu’il n’était pas approprié de déduire la déformation macroscopique à partir de la cinématique aux contacts. Dans ce cadre, ce travail propose d’introduire une troisième échelle dite mésoscopique. Cette échelle, à laquelle peuvent être définies à la fois contrainte et déformation, est intermédiaire entre les échelles microscopique et macroscopique et permet de palier au défi majeur mentionné ci-dessus. Elle est définie au niveau d’arrangements locaux de particules, appelés sous-domaines, et sa pertinence est étudiée sur la base d’échantillons numériques composés de particules circulaires puis sphériques, simulés par la méthode des éléments discrets. Les milieux bidimensionnels sont géométriquement représentés par un graphe de particules composé de sous-domaines fermés, encore appelés cellules de vide, dont la frontière est constituée de branches connectant les centres de particules en contact : l’échelle mésoscopique est donc définie au niveau de ces cellules de vide fermées. A cette échelle locale, on décrit tout d’abord la structure du milieu en termes de densité et de texture puis l’on définit les variables statique et cinématique locales du milieu en termes de contrainte et déformation. De fortes hétérogénéités des milieux granulaires en termes de structure, déformation et contrainte sont mises en évidence à l’échelle mésoscopique, avec de plus une structuration des hétérogénéités de contraintes et de déformations et une forte corrélation entre ces deux quantités. Concernant les milieux tridimensionnels, une partition en cellules de vide fermées est impossible du fait de la complexité de la structure 3D de ces milieux. On propose donc une méthode de partition du milieu basée sur la distribution des vides en son sein. La méthode consiste en premier lieu en une subdivision du milieu en tétraèdres, par une partition de Delaunay, puis en une association de tétraèdres voisins selon un critère prédéfini en vue de la création de sous-domaines, non fermés, mais au rôle analogue aux sous-domaines fermés de l’étude 2D. Le critère d’association proposé est basé sur le rapport entre la taille des constrictions (vide sur chaque face des tétraèdres) et la taille des pores au voisinage de chaque constriction. Cette méthode d’association constitue donc l’étape préliminaire à l’extension au cas tridimensionnel des résultats obtenus dans le cas bidimensionnel. / The technique of change of scales has been extensively developed in the literature to describe the global behaviour of granular materials taking into account their local properties. This method usually considers two scales : the macroscopic scale at the level of representative elementary volume and the microscopic scale at the level of contact between particles. The major difficulty of this "micro-macro" change of scales lies in the definition of the macroscopic strain : indeed, the macroscopic stress is clearly defined from contact forces while it is not appropriate to derive the macroscopic strain from the kinematics at contacts. In this framework, this work proposes to introduce a third scale called mesoscopic scale. This scale, at which both stress and strain can be defined, is intermediate between macroscopic and microscopic scales and allows to overcome the major difficulty mentioned above. The mesoscopic scale is defined at the level of local arrangements of particles, called sub-domains, and its relevance is studied on numerical 2D and 3D materials composed of circular then spherical particles, simulated with the discrete element method. Bidimensional media are geometrically represented by a particle graph composed of closed sub-domains, also called void cells, whose border is constituted by the branches joining the centers of particles in contact : the mesoscopic scale is thus defined at the level of these closed void cells. At this local scale, we fist describe the structure of the medium in terms of density and fabric ; we define then the static and kinematic variables in terms of stress and of strain. Strong heterogeneities of granular media in terms of structure, stress and strain are highlighted at this scale, with a structuration of heterogeneities of stress and strain and a significant correlation between these two quantities. Concerning tridimensional media, a partition into closed void cells is impossible, because of the complexity of the 3D structure of these media. We propose a partition method based on the distribution of voids inside the medium. This method consists first in subdividing the medium into tetrahedrons by a Delaunay partition and then in associating neighbouring tetrahedrons, according to a criterion to be defined. This allows us to form sub-domains which are not closed but which play a role analogous to the role of closed sub-domains in the 2D study. The proposed association criterion is based on the ratio between the size of constriction (void on each face of tetrahedrons) and the size of pores around each constriction. This partition method constitutes a preliminary step for an extension of the results obtained in the bidimensional case to the tridimensional case.

Ce travail de recherche a pour objectif d’étudier le comportement au cisaillement des matériaux granulaires homogènes et hétérogènes par l’essai triaxial. Les travaux sont réalisés sur des billes de verre et comportent à la fois des essais triaxiaux de laboratoire et des simulations numériques en DEM (Discret Element Method). D’un point de vue expérimental, les essais sont tout d’abord réalisés sur des échantillons homogènes (billes monodisperses). Ensuite, on étudie des systèmes de billes monodisperses comportant une inclusion de grande taille et des mélanges de billes bidisperses pour représenter les cas hétérogènes. D’un point de vue numérique, une condition limite cylindrique paroi rigide, basée sur la formule de Lamé, est intégrée dans un code aux éléments discrets existant. Une série de procédures est proposée afin de modéliser l’essai triaxial dans des conditions comparables aux essais. Les résultats numériques et expérimentaux sont ensuite comparés, à la fois pour les systèmes monodisperses et pour les milieux hétérogènes. Le modèle numérique permet de reproduire de manière très satisfaisante la courbe déviatorique dans toutes les conditions expérimentales utilisées. Il permet également de décrire correctement les déformations volumiques de l’échantillon, malgré la simplicité de la condition aux limites utilisée. L’outil numérique développé constitue ainsi un complément intéressant aux essais expérimentaux. Il permet en effet d’étudier l’influence de paramètres micromécaniques sur le comportement global des matériaux. / This research work aims at studying the shear behavior of homogeneous and heterogeneous granular materials by triaxial test. The work is performed on glass beads both in laboratory tests and by numerical simulations in DEM (Discrete Element Method). From an experimental point of view, tests are first performed on homogeneous samples (monodisperse beads), then expands to heterogeneous cases like monodisperse systems comprising a large inclusion and bidisperse mixtures. From a numerical point of view, a cylindrical rigid wall boundary condition, based on the Lamé formula is integrated into an existing discrete elements code. A series of procedures is proposed to model the triaxial test in conditions similar to experimental tests. The numerical and experimental results are compared both for monodisperse and heterogeneous systems. The numerical model can reproduce deviatoric curves very satisfactorily in all experimental conditions. It also allows describing correctly the volumetric strains of the sample despite the simplicity of the boundary condition used. The developed numerical tool is an interesting complement to experiment triaxial tests. It makes it possible to study the influence of micromechanical parameters on the overall behavior of materials.

La dysplasie acétabulaire est une pathologie causant de l'instabilité articulaire et une réduction de la superficie de chargement au niveau de la hanche. Celle-ci augmente le risque de développement précoce d'ostéoarthrose, causant de la douleur et menant à la perte de la fonction articulaire. La dysplasie est conventionnellement traitée par ostéotomie péri-acétabulaire, mais les effets précis de cette correction sur la biomécanique intra-articulaire restent à être élucidées. À l'aide de la méthode de simulation par éléments discrets, la distribution de pression articulaire a été évaluée de manière pré et post opératoire chez 9 patients dysplasiques. Différentes mesures orthopédiques évaluant le degré de correction acétabulaire ont été acquises et plusieurs scores numériques évaluant la surcharge articulaire ont été calculés. Des corrélations entre ces mesures et ces scores ont ensuite été effectués. Il a été noté que les mesures d'angle centre-externe latéral, d'index acétabulaire et d'index d'extrusion de la tête fémorale sont corrélés avec l'augmentation de l'aire de contact moyen durant la marche (p < 0.05, R² de 0.22 à 0.50). L'influence significative de la méthode de reconstruction cartilagineuse (distributions d'épaisseurs moyennes saine vs. dysplasique) et du profil de marche utilisé (phase d'appui moyen sain vs. dysplasique) a été remarqué auprès de la majorité des métriques évaluées. Ces résultats, ainsi que les facteurs qui les influencent, peuvent servir à assister la planification chirurgicale, effectuer des validations in-vitro d'optimisation d'ostéotomie péri-acétabulaire et, par la suite, guider la conception de nouveaux traitements. / Acetabular dysplasia causes articular instability and the reduction of the loading area in the hip. This causes the early onset of osteoarthritis which leads to pain and the eventual loss of articular function. Acetabular dysplasia is commonly treated by periacetabular osteotomy, but the precise effects of this correction on intra-articular biomechanics is not yet well known. Using discrete element analysis, pressure distribution has been calculated in the dysplastic hips of 9 patients before and after periacetabular osteotomy. A variety of radiographic measurements and numerical metrics have been evaluated in order to quantify the degree of acetabular correction and excessive loading. Correlations between these measures and metrics have been established in order to reveal any causal effects. In particular, the corrections of the lateral center-edge angle, acetabular index and femoral head extrusion index are correlated to the rise of average contact area during gait (p-value < 0.05, R² = 0.22 to 0.50). Additionally, the influence of cartilage layer reconstruction method (population-based healthy vs. dysplastic thicknesses) and loading parameters (healthy vs. dysplastic gait profiles) was significant for a majority of evaluated metrics. These results and influencing factors can help to assist with surgical planning, provide a basis for in-vitro validation of surgical optimization and act as a standard for designing new orthopedic treatments for acetabular dysplasia.

Les anodes en carbone font partie de la réaction chimique de réduction de l'alumine, qui est consommée lors du procédé d'électrolyse Hall-Héroult. Le comportement des agrégats secs de coke en tant que composant principal des anodes de carbone (environ 85 %) a un rôle clé et exceptionnel dans leurs propriétés finales. L'analyse de défaillance des agrégats de coke sec permet non seulement de mieux comprendre les mécanismes de déformation des matériaux granulaires sous charge compressive, mais peut également identifier les causes potentielles de défauts structurels des anodes en carbone, telles que les fissures horizontales. Dans ce travail, il sera montré qu'un mode de défaillance particulier peut être responsable de la génération de fissures dans les anodes en carbone. Le comportement de rupture des agrégats de coke n 'est pas seulement affecté par les paramètres du processus de compactage, tels que la pression de confinement et la vitesse de déformation axiale, mais il dépend également fortement de la distribution granulométrique et de la forme des particules de coke. La méthode des éléments discrets (DEM) est utilisée pour modéliser le comportement micromécanique des agrégats de coke sec pendant le processus de compactage. De plus, le critère de travail de second ordre est utilisé pour analyser la rupture des éprouvettes de granulats de coke. Les résultats révèlent que l'augmentation de la pression de confinement augmente la probabilité du mode de diffusion de la rupture dans l'éprouvette. D'autre part, l'augmentation de la vitesse de déformation augmente les chances du mode de localisation de la déformation de la rupture dans l'éprouvette. De plus, les résultats indiquent que l'utilisation de fines particules ainsi que la diminution de la sphéricité des particules de coke augmenteront la plage de stabilité des agrégats de coke. De plus, en utilisant l'analyse des évaluations de contour de micro-déformation pendant le processus de compactage, il est montré que, à la fois en ajoutant des particules fines aux agrégats de coke et en diminuant la sphéricité des particules de coke, la possibilité de créer une bande de compression dans le coke agrégats est réduit. Étant donné que la présence des bandes de compactage dans la pâte d'anode crée une zone sujette à la génération de fissures horizontales, les résultats de cette étude pourraient conduire à la production d'anodes en carbone avec moins de défauts structurels. / Carbon anodes are part of the chemical reaction of alumina reduction, that is consumed during the Hall-Héroult electrolysis process. The behavior of dry coke aggregates as the main component of carbon anodes (about 85 %) has an exceptional key role in their final properties. The failure analysis of dry coke aggregates not only leads to a better understanding of the deformation mechanisms of granular materials under compressive loading but also can also identify the potential causes of structural defects in carbon anodes, such as horizontal cracks. In this work, it will be shown that a particular failure mode can be responsible for the crack generation in the carbon anodes. The failure behavior of the coke aggregates is not only affected by the compaction process parameters, such as the confining pressure and axial strain rate, but it is also strongly dependent on the size distribution and shape of coke particles. The discrete element method (DEM) is employed to model the micro-mechanical behavior of the dry coke aggregates during the compaction process. In addition, the second-order work criterion is used to analyze the failure of the coke aggregate specimens. The results reveal that increasing the confining pressure enhances the probability of the diffusing mode of the failure in the specimen. On the other hand, the increase of the strain rate augments the chance of the strain localization mode of the failure in the specimen. In addition, the results indicate the fact that the use of fine particles as well as decreasing the sphericity of coke particles will increase the stability range of the coke aggregates. Moreover, by using the analysis of micro-strain contour evaluations during the compaction process, it is shown that, both by adding fine particles to the coke aggregates and by decreasing the sphericity of coke particles, the possibility of creating a compression band in the coke aggregates is reduced. Since the presence of the compaction bands in the anode paste creates an area that is prone to horizontal crack generation, the results of this study could lead to the production of carbon anode with fewer structural defects.

Les structures de sphères creuses appartiennent à la famille des matériaux cellulaires qui ont récemment été étudiés pour leurs multiples propriétés. Dans le cas de cette thèse, le but des sphères creuses est de dissiper l’énergie d’impact d’un oiseau sur un cockpit d’avion. Elles sont développées dans le cadre du projet SAMBA (Shock Absorber Material for Birdshield Application) afin d’optimiser leur énergie spécifique absorbée (J/kg).Dans un premier temps, des essais quasi-statiques (v = 5 mm/min) et dynamiques (v = 2 m/s) de compression uni-axiale sont conduits à température ambiante sur une seule sphère creuse de diamètre 30 mm. Une propagation rapide de fissures macroscopiques est observée. Le formalisme de la Mécanique Élastique Linéaire de la Rupture (MELR) est utilisé pour estimer le taux de restitution d’énergie critique dynamique GIdc du matériau constitutif. La position du sommet de fissure est mesurée pendant la propagation de fissure à l’aide d’une caméra rapide. La Méthode des Éléments Discrets (DEM) permet de simuler la rupture dynamique en implémentant une technique de relâchement des nœuds. Le taux de restitution d’énergie GIdc peut être estimé à partir de l’histoire (position et temps) du sommet de fissure. Le modèle numérique montre que les structures sphériques dissipent une proportion importante de l’énergie par des effets dynamiques. A une même vitesse de propagation, plus l’épaisseur de coque est fine, plus les effets inertiels générés par la rupture sont importants et ce pour une même vitesse de propagation.Le modèle numérique DEM est ensuite employé pour reproduire la rupture dynamique sur une sphère creuse à l’aide d’un critère en contrainte seule ou un critère mixte en contrainte – énergie. Les bons résultats obtenus démontrent la capacité de la DEM à représenter la propagation de fissures en régime dynamique.Finalement, des essais numériques et expérimentaux multi-sphères sont réalisés afin évaluer le comportement des sphères creuses au sein d’un assemblage. / Hollow sphere structure (HSS) belongs to cellular solids that have been studied recently for its multiples properties. In our case, HSS aims to absorb soft impacts energy on an airliner cockpit. HSS is investigated through the SAMBA (Shock Absorber Material for Bird-shield Application) project because of its promises in term of specific energy dissipated (J/kg) during impact.First of all, quasi-static and dynamic (v = 5 mm/min to v = 2 m/s) uniaxial compression tests are conducted at room temperature on a single sphere (D = 30 mm). Rapid crack propagation (RCP) is observed to be predominant at macroscopic scale. The formalism of Linear Elastic Fracture Mechanics (L.E.F.M.) is therefore used to estimate the dynamic energy release rate GIdc . The crack tip location is measured during the crack propagation using a high speed camera. The Discrete Element Method (DEM) is used to simulate the dynamic fracture by implementing the node release technique. The dynamic energy release rate can be determined using an experimentally measured crack history. In spherical structures the numerical results reveal a high proportion of energy dissipated through inertial effects as well as a dependence of the thickness of the hollow sphere over the range of 0.04 mm to 1.2 mm.The DEM model Is then employed to reproduce the RCP according to two failure criterions: a stress criterion and a coupled stress-energy criterion. It reveals to be an interesting way to model the mechanical behavior of brittle materials.Eventually, experimental and numerical multi-spheres tests are performed to evaluate the behavior of brittle hollow spheres within an assembly.

Les systèmes de stockage d'énergie thermique sont des éléments centraux de divers types de centrales de production d'énergie fonctionnant à l'aide de sources d'énergie renouvelables. Ces systèmes de stockage thermique de type solide/fluide peuvent être considérés comme une solution rentable dans les centrales solaires à concentration. Un tel dispositif est constitué d'un réservoir rempli d'un lit granulaire à travers lequel circule un fluide caloporteur. Cependant, dans de tels dispositifs, le réservoir pourrait être soumis à une accumulation de contraintes thermiques au cours des cycles de chargement et de déchargement en raison de la dilatation thermique différentielle entre le milieu granulaire et la paroi du réservoir. Cette thèse a été consacrée à l'étude du comportement thermomécanique du lit granulaire à l'intérieur d'un réservoir de stockage thermique de type régénérateur solide/fluide. Pour atteindre cet objectif, deux approches ont été suivies : les approches numériques et expérimentales. Un modèle numérique basé sur la méthode des éléments discrets a été développé pour décrire le comportement du réservoir sous cyclage thermique. L'évolution des contraintes appliquées aux parois de la cuve au cours des cycles thermiques, est étudiée en tenant compte à la fois des charges thermiques et mécaniques, ainsi que de la cinématique du milieu granulaire à l'échelle des particules (c'est-à-dire des éléments discrets) et de la déformation thermo-mécanqiue du réservoir. Des simulations ont été effectuées pour différentes configurations thermiques (chauffage/refroidissement homogène du réservoir sur sa hauteur ou gradient vertical de température) et différents cas de condition limite (paroi rigide, paroi du réservoir avec un coefficient de dilatation thermique supérieur au lit de roche ou inversement). Le comportement du réservoir dépend des conditions thermiques et conditions aux limites imposées. De plus, un prototype de thermocline a été conçu et construit au sein du laboratoire au CEA, visant à étudier l'accumulation de contraintes sur les cycles thermiques. Le dispositif expérimental, appelé ESPERA, est équipé de dispositifs de mesure de force installés à différentes hauteurs sur la paroi du réservoir. Les dispositifs de mesure de force ont été développés et calibrés au CEA. Leur sensibilité a également été testée en développant un autre banc expérimental, P'tit- Pousse. Des essais expérimentaux ont été réalisés permettant de prouver l'accumulation des contraintes au cours des cycles. Finalement, une comparaison entre les résultats numériques et les premiers résultats expérimentaux obtenus en fin de thèse est proposée. / Thermal Energy Storage (TES) systems are central elements of various types of power plants operated using renewable energy sources. Packed bed TES can be considered as a cost effective solution in concentrated solar power plants (CSP). Such a device is made of a tank filled with a granular bed through which a heat-transfer fluid circulates. However, in such devices, the tank might be subjected to an accumulation of thermal stresses during cycles of loading and unloading due to the differential thermal expansion between the filler and the tank wall. This research was devoted to investigate the thermo-mechanical behavior of the granular bed inside a packed bed TES tank. To achieve this objective, two approaches were undertaken in this work, i.e. numerical and experimental. A numerical model was defined to describe the tank's behavior under thermal cycling based on the discrete element method (DEM). The evolution of tank wall stresses over thermal cycles, taking into account both thermal and mechanical loads, as well as the kinematics of the granular material at the particles scale (i.e. discrete elements), are studied here. The deformability of the tank itself under thermo-mechanical loads is also included in the numerical model. Simulations were performed for different thermal configurations (i.e. the tank is heated homogeneously along its height or with a moving vertical gradient of temperature) and different boundary condition cases (i.e. rigid wall, tank wall with a higher thermal expansion coeficient than the bed or inversely). The behavior of the tank is dependent on the imposed thermal and boundary conditions. In addition to this, a thermocline prototype was designed and constructed at the CEA laboratory, aiming at studying the stress accumulation over the thermal cycles. The experimental setup, called ESPERA, is equipped with force measurement devices, installed at different height positions on the wall. The force measurement devices were developed and calibrated at the CEA. Their sensitivity was also tested using a different developed setup, P'tit-Pousse. Experimental tests were carried out proving the stress accumulation over the cycles. Eventually, a comparison between the numerical results and the lately-obtained preliminary experimental measurements is proposed.

Le renforcement de sols par inclusions rigides est une technique de fondation dans laquelle un matelas granulaire de transfert de charge est mis en \oe uvre entre l'ouvrage et le sol support. Les reports de charges se développant dans le matelas permettent d'augmenter la charge de l'ouvrage transmise au réseau d'inclusions mis en place dans le sol support. Une importante réduction des tassements peut ainsi être obtenue. L'absence de recommandations relatives au dimensionnement et à la réalisation de de ce type de renforcement de sol a conduit à la mise en place du projet national ASIRi. Une des étapes essentielles à la compréhension du fonctionnement de ce type d'ouvrage réside dans les mécanismes de transfert de charge agissant dans le matelas granulaire qui demeurent mal connus et complexes. Afin d'apporter des éléments de compréhension relatifs à ces mécanismes, une étude expérimentale mettant en jeu des matelas granulaire soumis à une sollicitation simplifiée a été réalisée sur modèle réduit. Dans un second temps, l'étude expérimentale a été complétée par des modélisations numériques utilisant la Méthode des Éléments Discrets. Enfin, une étude paramétrique et numérique portant sur le comportement de matelas de transfert de charge sur inclusions rigides est proposée. Elle se base sur des configurations d'ouvrages usuelles et permet d'évaluer l'influence de différentes solutions techniques sur les reports de charge.

L'objectif de ce travail de thèse a été de modéliser avec une méthode aux éléments discrets la fracturation en tension, et plus particulièrement la fragmentation dynamique, sur des matériaux fragiles pour de hautes vitesses de déformation. La fragmentation est un phénomène irréversible, non linéaire et aléatoire. Elle intervient dans de nombreux domaines de la vie courante, quelque soit l'échelle considérée. La modélisation numérique de ce phénomène permettrait une prédiction de certains paramètres statistiques de la fragmentation, comme le nombre de fragments, la taille des fragments, la distribution de la taille des fragments, etc. Pour cette thèse, la Méthode des Éléments Discrets (DEM) s'est révélée être un excellent moyen pour simuler la fracturation en raison de sa nature discrète. Toutefois, une bonne méthode de simulation numérique ne suffit pas à elle seule pour modéliser la fragmentation dynamique. Un critère de rupture doit également être inséré, afin d'introduire un endommagement. Ce critère de rupture s'écrit au niveau d'un lien entre deux particules et il engendre un dommage, en faisant décroître la contrainte locale jusqu'à l'obtention d'une fissuration discrète. Dans un premier temps, un critère de rupture de Camacho-Ortiz |24| a été introduit dans une méthode aux éléments discrets. Ce critère se traduit par un endommagement en fonction d'une ouverture de fissure. Ce premier critère a donné de bons résultats comparé à ceux de |69, 88, 97, 143-147| sur la convergence des paramètres de la fragmentation sur des cas simples, mais nécessite un grand nombre de particules. Dans un second temps, afin d'envisager la modélisation de la fragmentation sur des cas plus complexes en trois dimensions à de hautes vitesses de déformation, un second critère de rupture a été introduit. Ce critère de rupture s'appuie sur une approche physique différente, qui prend en compte l'hétérogénéité des matériaux fragiles avec leurs défauts susceptibles d'évoluer et de provoquer une rupture locale. Pour cela, il fait intervenir une loi probabiliste de Weibull afin d'introduire des défauts par élément de volume. Ce critère a été développé par C. Denoual, P. Forquin et F. Hild |29, 33, 42-44|. Tout d'abord, ce second critère de rupture a été testé sur des cas simples en obtenant une convergence des paramètres statistiques de la fragmentation avec un nombre environ $10$ fois moins important de particules que pour la première méthode. Un cas plus complexe en trois dimensions de modélisation de barre d'Hopkinson en trois dimensions a permis de tester de manière " qualitative " la méthode.