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.

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.

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.

Les forces de contact dans un milieu granulaire se répartissent de façon très inhomogène en un réseau de "chaînes de forces" qui supporte la plus grosse partie des contraintes. Il est primordial de bien comprendre l'influence du désordre des forces à l'échelle du grain sur les propriétés d'équilibre mécanique d'un milieu granulaire à l'échelle macroscopique.<br /><br />La mesure de forces dans un milieu granulaire est délicate à cause d'un couplage fondamental entre les déformations du capteur et la mobilisation des forces de friction entre grains ou entre grains et paroi. Cependant, en définissant proprement le protocole de mesure, nous avons pu améliorer de façon très significative la reproductibilité des résultats en comparaison des mesures que l'on trouve dans la littérature.<br /><br />Nous nous sommes intéressé aux liens qui existent entre la structure de l'empilement granulaire et la répartition des contraintes. Sous l'effet d'un cisaillement ou de vibrations, l'équilibre d'une colonne granulaire dans un silo évolue considérablement ainsi que la structure de l'empilement comme le révèlent des mesures de densité moyenne et locale. Sous le sommet d'un tas de sable formé par écoulement des grains en avalanches, j'ai observé très clairement un minimum ou "trou" de pression, alors que la pression est maximum si les grains sont déposés en couches horizontales. Nous avons aussi mesuré les fluctuations résiduelles de la pression en fonction de la taille des grains ou de la hauteur de remplissage du silo et ai observé que leur dépendance avec la taille des grains montrent une régression statistique anormale en comparaison de celle déduite de la distribution des forces à l'échelle du grain.<br /><br />La plupart des observations sont qualitativement et quantitativement bien reproduites par le modèle OSL dans lequel les contraintes se propagent selon deux directions dont l'une s'identifie à la direction moyenne des chaînes de forces. La notion de propagation anisotrope de forces permet de comprendre la distribution de pression sous un tas ou la forme de la courbe de saturation de la pression dans un silo, y compris les oscillations de la pression en présence d'une surcharge.<br /><br />Nous avons proposé en outre un modèle de durcissement d'arches qui permet d'expliquer l'existence d'un mode d'écoulement fragmentaire après renversement d'un tube rempli de grains et met en évidence le rôle non négligeable de l'élasticité des parois.

L'analyse de la pression exercée par un écoulement granulaire sur un obstacle nécessite la connaissance du coefficient de traînée. Néanmoins, les études bibliographiques expérimentales ne fournissent pas les données internes de l'écoulement intervenant dans son calcul. C'est pourquoi, nous proposons la détermination du coefficient de traînée par l'application d'outils numériques tridimensionnels, basés sur la méthode des éléments discrets (MED), et validés par des études expérimentales. L'expérience est un canal d'écoulement de laboratoire, intégralement modélisé par la MED. Les lois locales de contact intègrent un comportement normal élastique, hystérétique associé à un critère de glissement tangentiel. La validation du modèle intégral repose sur la similarité entre les résultats expérimentaux et numériques des paramètres physiques de l'écoulement et de l'effort d'impact sur l'obstacle. Cette validation est obtenue sans calibration importante des paramètres numériques. Le modèle intégral est ensuite optimisé, en termes de temps de calcul et de possibilités d'études, sous la forme d'un canon granulaire, dans lequel n'est simulée que la zone d'interaction entre l'écoulement et l'obstacle. Les propriétés de l'écoulement : vitesse, épaisseur, densité, sont contrôlées car intégrées en tant que paramètres initiaux numériques. L'influence de caractéristiques associées à l'écoulement, comme le nombre de Froude, à la taille des grains, ou à la forme et la taille de l'obstacle, sur le coefficient de traînée, est alors analysée. Enfin, une extension micromécanique est proposée dans l'analyse de l'influence d'une adhésion et d'une cohésion intégrées aux lois locales de contact.

L'utilisation d'un matériau granulaire au sein de structures industrielles afin de diminuer les niveaux vibratoires est une solution alternative aux revêtements viscoélastiques, qui ont une dépendance forte aux conditions de température. Pour l'industrie navale et aéronautique l'enjeu est ainsi d'améliorer la fiabilité et le confort. Les niveaux d'interaction entre les grains dépendent des paramètres de contrôle (niveaud'accélération), mais aussi des caractéristiques des matériaux constitutifs des grains, du taux d'humidité dufluide environnant, de la géométrie et des dimensions des grains, ainsi que des conditions de confinement.Pour une accélération donnée, indépendante du point de la structure, la pertinence du choix d'unmatériau granulaire par rapport à une même masse indéformable est mise en évidence par l'étude de lavibration d'une cavité rigide montée sur un oscillateur partiellement remplie. La variation paramétrique destypes de confinements dans la cavité ainsi que le matériau constitutif des grains permet d'extraire lesparamètres influents dans la dissipation d'énergie. Le développement d'un modèle analytique à constanteslocalisées permet de simuler l'énergie dissipée par le système via un nombre réduit de coefficients représentantla dynamique vibratoire de la matière en grains. La confrontation de ce modèle avec différentesexpérimentations permet de valider son efficacité et son caractère prédictif dans la dissipation de l'énergievibratoire d'un système dynamique.Pour une structure modale en vibration, où le niveau d'accélération est dépendant du point de lastructure, l'optimisation du positionnement des amas de grain est montrée. / The use of granular material in industrial structures to reduce vibration levels is an alternative toviscoelastic surfacing, which is highly dependent on temperature conditions. For the naval and aeronauticalindustry, the challenge is to improve reliability and comfort. The levels of interaction between the grainsdepend on the control parameters (acceleration level), but also on the characteristics of the constituentmaterials of the grains, the moisture content of the surrounding fluid, the geometry and dimensions of thegrains, as well as their containment conditions.For a given acceleration, independent of the point of structure, the relevance of the choice of a granularmaterial with respect to the same non-deformable mass is demonstrated by the study of the vibration of a rigidcavity mounted on an oscillator partially filled. The parametric variations of the types of confinement in thecavity as well as the constituent material of the grains make it possible to extract the influent parameters inthe energy dissipation. The development of an analytical model with localized constants allows to simulate theenergy dissipated by the system via a reduced number of coefficients representing the vibratory dynamics ofthe granular matter. The comparison of this model with experiments makes it possible to validate itseffectiveness and its predictive character in the dissipation of the vibratory energy of a dynamic system.For a modal structure, where the acceleration level is dependent on the point of the structure, theoptimization of the grain cluster positioning is shown.

Les simulations numériques à l'échelle du pore sont fréquemment utilisées pour étudier le comportement des écoulements multiphasiques largement rencontrées dans phénomènes naturels et applications industrielles. Dans ce travail, la morphologie de structures liquides et l'action capillaire sont examinées à l'échelle des pores par la méthode de Boltzmann sur réseau (LBM) à plusieurs composants selon le modèle de Shan-Chen. Les résultats numériques obtenus sont en bon accord avec les solutions théoriques. Les simulations numériques sont étendues à microstructures complexes au-delà du régime pendulaire.La LBM a été utilisée pour modéliser l'écoulement multiphasique à travers un milieu poreux idéalisé dans des conditions de drainage primaire quasi-statique. Les simulations LBM ont fourni une excellente description du déplacement de l'interface fluide-fluide à travers les grains. Pendant le drainage, les simulations LBM sont capables de reproduire la déconnexion d'une phase dans le milieu granulaire sous la forme de ponts pendulaires ou structures liquides complexes. Malheureusement, le temps de calcul nécessaire pour ce type de simulations est assez élevé. Afin d’optimiser les ressources de calcul, nous présentons un modèle 2D (modèle Throat-Network) basé sur des solutions analytiques pour décrire l'écoulement biphasique à travers un ensemble de disques dans un temps de calcul très réduit, donc le modèle 2D est susceptible de remplacer les simulations LBM lorsque les ressources de calcul sont limitées. L'approche souligne l'importance de simuler le problème a l'échelle de la gorge du pore pour obtenir les relations volume - pression capillaire locales. Le modèle Throat-Network est un point de départ pour le modèle hybride proposé pour résoudre les problèmes en 3D. Le modèle hybride combine l’efficacité de l’approche réseau de pores et la précision du LBM à l’échelle des pores. Le modèle hybride est basé sur la décomposition de l’échantillon en petits sous-domaines, dans lesquels des simulations LBM sont effectuées pour déterminer les propriétés hydrostatiques principales (pression capillaire d'entrée, courbe de drainage primaire et morphologie du liquide pour chaque gorge du pore). Malgré la réduction significative des temps de calcul obtenus avec le modèle hybride, le temps n’est pas négligeable et les modélisations numériques d'échantillons de grandes tailles ne sont pas réalistes. Les approximations données par les méthodes Incircle et MS-P, qui prédisent les propriétés hydrostatiques, sont comparées à celles de LBM et du modèle hybride. / Numerical simulations at the pore scale are a way to study the behavior of multiphase flows encountered in many natural processes and industrial applications. In this work, liquid morphology and capillary action are examined at the pore-scale by means of the multicomponent Shan-Chen lattice Boltzmann method (LBM). The accuracy of the numerical model is first contrasted with theoretical solutions. The numerical results are extended to complex microstructures beyond the pendular regime.The LBM has been employed to simulate multiphase flow through idealized granular porous media under quasi-static primary drainage conditions. LBM simulations provide an excellent description of the fluid-fluid interface displacement through the grains. Additionally, the receding phase trapped in the granular media in form of pendular bridges or liquid clusters is well captured. Unfortunately, such simulations require a significant computation time. A 2D model (Throat-Network model) based on analytical solutions is proposed to mimic the multiphase flow with very reduced computation cost, therefore, suitable to replace LBM simulations when the computation resources are limited. The approach emphasizes the importance of simulating at the throat scale rather than the pore body scale in order to obtain the local capillary pressure - liquid content relationships. The Throat-Network model is a starting point for the a hybrid model proposed to solve 3D problems. The hybrid model combines the efficiency of the pore-network approach and the accuracy of the LBM at the pore scale to optimize the computational resources. The hybrid model is based on the decomposition of the granular assembly into small subsets, in which LBM simulations are performed to determine the main hydrostatic properties (entry capillary pressure, capillary pressure - liquid content relationship and liquid morphology for each pore throat). Despite the reduction of computation time, it is still not negligible and not affordable for large granular packings. Approximations by the Incircle and the MS-P method, which predict hydrostatic properties, are contrasted with the results provided by LBM and the hybrid model. Relatively accurate predictions are given by the approximations.

Carbone Savoie est une société française spécialisée dans les solutions cathodiques en carbone et graphite pour l’électrolyse de l’aluminium primaire. La fabrication de ces produits requiert un savoir-faire important en raison de la complexité de procédé. Le malaxage, première étape du procédé, permet le mélange intime des particules carbonées et d’un liant. Les grains carbonés de tailles variées comme le coke calciné, anthracite calciné, et le graphite constituent un empilement multimodal d’une dizaine de microns à quelques centimètres. Le liant actuel est liquide à haute température, se caractérise par un haut rendement en carbone, est cokéfiable mais toxique. La législation européenne sur les produits chimiques (REACH) l’a identifié comme « substance extrêmement préoccupante », le risque associé étant une interdiction de son utilisation. Dans ce contexte, Carbone Savoie a lancé un projet de recherche visant à substituer cette matière première dangereuse. Le premier objectif a été de déterminer des compositions possibles pour ce nouveau liant ainsi que leur comportement rhéologique. Une première composition a montré un comportement newtonien à pH basique. Un additif à point de gélification, dont le comportement a fait l’objet d’une étude, a été utilisé pour une deuxième composition. La meilleure méthode de mélange fut ensuite sélectionnée puis son efficacité avec le nouveau liant validée. Une troisième composition, qui est une suspension, a montré un comportement rhéofluidifiant. Cette dernière s’est avérée être la meilleure candidate pour substituer le brai. Le deuxième objectif du travail a été la détermination des conditions de malaxage avec le nouveau liant. L’étude du mélange granulaire carboné est basée sur une méthode de suivi d’intensité de malaxage. Grâce à la courbe d’intensité enregistrée en fonction du ratio massique liquide/solide, les différents états de l’eau dans le mélange en fonction de la saturation en liant du système ont été définis : régime pendulaire, funiculaire, capillaire et goutte. Le régime capillaire est particulièrement intéressant car il s’agit du régime d’extrusion des cathodes. L’influence sur le procédé (cinétique de mouillage, étendue du régime capillaire) de la viscosité du liant a été étudiée : une viscosité élevée permet une plus grande tolérance du procédé vis-à-vis des variations de l’empilement granulaire des matières sèches d’une formulation à l’autre ou encore la difficulté de produire des poudres de granulométrie fixe. Enfin, une corrélation entre les résultats obtenus et l’échelle industrielle a été trouvée. Les propriétés finales d’une cathode peuvent être reliées à la qualité de la pâte en sortie de malaxeur grâce à la détermination de la cohésion via une cellule de cisaillement annulaire. / Carbone Savoie is a French company specialized in the production of carbon and graphite cathodic solutions for the primary aluminum industry. The production of cathodes requires an important know-how because of its complex manufacturing process. The first step of the process consists in mixing the carbonaceous raw materials with a binder. The raw materials (calcined coke and anthracite or graphite) have different size and represent a multimodal carbon aggregate (from tens of microns to a few centimeters). The current binder is liquid at high temperature, has a high carbon yield and is coking but poisonous. The European regulations on chemical products (REACH) identified the pitch binder as “substance of very high concern”. Therefore, that product could become prohibited. Within this context, Carbone Savoie initiated a research project in order to substitute that hazardous material. First of all, the project aims at determining possible compositions for the new binder and to describe their rheological behavior. A first formulation had a Newtonian behavior for high pH values. An additive with a gel point which behavior has been studied was used for a second formulation. The best mixing process was then chosen and its efficiency was tested with this new binder. A third formulation which was a suspension showed a shear-thinning behavior. It appeared to be the best candidate to substitute pitch as a binder for the production of cathode blocks. The second goal of the project was the determination of the mixing conditions with this new binder. The study of the mixing of the carbon aggregate with the binder is based on a tracking of the mixing intensity. Thanks to the plot of the intensity-liquid/solid ratio curve, different states were defined for the system according to the binder content (pendular state, funicular state, capillary state and solid/liquid dispersion state). The capillary state is of a major concern because the paste is extruded in that state. The impact on the process (wetting kinetics, width of the capillary state) of the rheological behavior of the binder has been studied: a high viscosity allows a better tolerance to the variability in the size of the carbon aggregate from a formulation to another and also to the trouble of producing powders with very small particles sizes. Finally, a correlation between these results and the industrial scale has been found. Final properties of cathodes blocks can be linked with the paste quality just after the mixing step thanks to the determination of the paste cohesion with a ring shear tester.

Un matériau granulaire confiné dans une colonne s'appuie sur les parois de celle-ci pour s'opposer aux forces extérieures qu'il subit. Il en résulte un couplage entre propriétés mécaniques du matériau et propriétés de surface de la colonne. L'étude de ce couplage constitue l'objet de cette thèse. Dans une première expérience, nous avons étudié le poids exercé par un milieu granulaire constitué de billes à la base d'une colonne. Nous avons défini une procédure expérimentale telle que le matériau est au seuil de glissement en tout point aux parois. Le poids mesuré sature avec la hauteur de billes ; la masse de saturation se met à l'échelle avec le rayon au cube, et diminue lorsque la densité et la friction aux parois augmentent. Lorsqu'on ajoute un surpoids à la surface supérieure du milieu granulaire, on observe un maximum, suivi d'une lente décroissance de la masse pesée en fonction de la masse versée. Ce maximum est plus élevé lorsque la densité et la friction aux parois augmentent. Les résultats obtenus sont comparés à différentes théories mécaniques des matériaux granulaires. Le modèle de Janssen décrit parfaitement le matériau granulaire non contraint par un surpoids. Mais l'ajout d'un surpoids met en défaut ce modèle, ainsi que le modèle OSL et l'élasticité homogène isotrope, même si cette dernière reproduit qualitativement tous les phénomènes observés. Dans une deuxième expérience, nous avons poussé la colonne granulaire vers le haut. La force de résistance à la poussée croît alors exponentiellement avec la hauteur de billes, ce que nous expliquons par le modèle de Janssen adapté à la poussée. A basse vitesse, nous observons un mouvement de « stick-slip », caractérisé par une force maximale de résistance à la poussée avant glissement qui croît fortement lorsque la vitesse décroît et lorsque l'humidité augmente. Nous expliquons ce phénomène par un vieillissement de la friction statique aux parois. A haute vitesse, nous observons une transition hystérétique vers un régime de glissement continu caractérisé par une force de résistance qui croît avec la vitesse et le taux d'humidité. Les phénomènes observés sont différents lorsqu'on change le matériau de la colonne, ce qui prouve l'influence primordiale des parois sur les effets observés. La mesure directe du coefficient de friction bille-paroi permet de retrouver les phénomènes observés. La poussée d'un système de billes de verre polydisperses dans une colonne en PMMA révèle un stick-slip régulier à très basse vitesse, qui devient irrégulier à plus haute vitesse, la distribution des glissements s'élargissant vers les petites valeurs.

Cette étude a pour objectifs la caractérisation expérimentale des frottements saccadés dans les matériaux granulaires modèles constitués des billes de verre monodisperses en compression triaxiale drainée et l'identification des paramètres de contrôle. Cinq paramètres macroscopiques caractérisent ces frottements saccadés : la chute de déviateur et la contraction volumique, l'intermittence de déformation axiale, le module d'Young et le coefficient de Poisson. Les frottements saccadés affectent simultanément le déviateur et la déformation volumique. Le comportement macroscopique est globalement contractant tandis que le matériau tend vers un état limite critique en grandes déformations, à la manière des sables lâches. Cependant il présente localement, dans les phases de blocage qui suivent immédiatement les ruptures temporaires, le comportement dilatant des sables denses, qui obéit à une relation contrainte-dilatance linéaire et unique. Les frottements saccadés disparaissent au-delà d'une vitesse critique d'écrasement axial, qui dépend de la contrainte de confinement et de la taille des grains. Le module d'Young dynamique par propagation d'ondes varie avec la contrainte de confinement selon une loi de puissance. Le module d'Young quasi-élastique au départ des phases de blocage est constant à l'intérieur du domaine élastique, de même que le coefficient de Poisson. Le suivi par granulométrie laser et par analyse d'images des matériaux après un ou plusieurs essais triaxiaux permet de suivre l'évolution de la taille moyenne et de la forme des grains. Tandis que les instabilités par saccade disparaissent suite à un certain nombre d'essais, on observe, simultanément à une légère diminution du volume moyen, l'apparition progressive de populations d'objets non sphériques par une fusion des grains analogue au frittage.