Le ballast est la couche granulaire, formée de grains rocheux irréguliers de taille centimétrique, qui supporte la voie ferrée par l'intermédiaire des traverses. Son rôle est crucial dans la transmission des efforts statiques et dynamiques imposés par la circulation des trains vers le sol supportant la voie. Pour comprendre le phénomène de tassement du ballast, qui croît fortement avec la vitesse des trains, et qui entraîne de coûteuses opérations de maintenance, il est important de bien connaître l’influence de la morphologie des grains sur le comportement mécanique du ballast. C’est l’objectif de cette thèse, dans laquelle est développée une technique innovante de caractérisation morphologique fondée sur la décomposition orthogonale aux valeurs propres (« POD »). Cette technique, appliquée à des données expérimentales détaillées (scans 3D) permet d’identifier un choix optimal de fonctions de forme. Il est alors possible de générer numériquement des assemblages granulaires représentatifs du ballast avec différents degrés de raffinement selon le nombre de fonctions de base retenu. Pour passer aux applications mécaniques, les grains, ainsi modélisés par un nombre limité de coefficients, doivent ensuite par une opération de maillage être assimilés à des polyèdres convexes dont le comportement mécanique est simulé par la méthode des éléments discrets (dynamique des contacts, avec le logiciel LMGC90). On détermine alors des caractéristiques de base comme l'angle de frottement interne et la compacité dans l'état critique de cisaillement quasi-statique maintenu. Ces propriétés sont stabilisées à partir d’un nombre assez faible de fonctions de forme (de l’ordre de 10). Dans le domaine de la description morphologique, l’approche fondée sur la POD ouvre d’intéressantes perspectives pour améliorer les caractérisations usuelles. Elle permet de différencier les ballasts selon leur degré d’usure ou la carrière dont ils proviennent. Elle conduit également à une méthode d’identification des faces et des courbures des grains / The ballast is the granular layer formed of irregular centimeter-sized rock grains, which supports the railway through the sleepers. Its role is crucial in the transmission of static and dynamic forces imposed by the movement of trains towards the ground supporting the track. To understand the phenomenon of ballast settlement, which increases significantly with the speed of the trains, and which requires costly maintenance operations, it is important to understand the influence of the grain morphology on the mechanical behavior of the ballast. Such is the objective of this thesis. First, an innovative method for particle shape characterization is developed, based on the Proper Orthogonal Decomposition (POD). When applied to highly accurate experimental data (3D scans of ballast grains), this technique enables the identification of an optimal set of shape functions, such that the shape information is condensed in a minimum number of coefficients. It becomes then possible to numerically generate granular assemblies that are representative of the morphology of ballast grains with different degrees of refinement according to the number of retained shape functions. In order to investigate their mechanical properties, the grains modeled by a limited number of coefficients have to be subjected to a meshing procedure and turned into convex polyhedra, suitable for numerical simulations of mechanical behaviours via the discrete element method (contact dynamics, using the LMGC90 software). We determine some fundamental characteristics such as the internal friction angle and the solid fraction in the critical state of quasistatic steady shear flow. These properties are stabilized after including a rather small number of shape functions in to the model (on the order of 10).The method also opens promising avenues for morphological characterization applications. It proves able to distinguish ballast grains with varying degrees of wear, or produced in different quarries. It also leads to a procedure for identifying faces and curvatures

Nous étudions la rhéologie d'un ensemble de grains sphériques et frottants, par la simulation numérique, à l'échelle des grains, d'écoulements de cisaillement sous une contrainte normale P contrôlée, en présence d'un liquide interstitiel. En faible teneur, ce liquide se présente sous forme de ménisques intergranulaires qui transmettent des forces capillaires attractives ; s'il sature l'espace intergranulaire, on s'intéresse alors à l'écoulement de Stokes de la suspension dense ainsi constituée, où dominent les forces visqueuses. Les assemblages de grains secs constituent un système de référence aux propriétés mécaniques bien connues, en particulier l'approche de l'état critique de la mécanique des sols dans la limite quasi-statique. L'effet des ménisques capillaires qui joignent les grains en présence d'un liquide en faible saturation (régime pendulaire) est étudié pour les taux de cisaillement allant du régime quasi statique au régime inertiel. La rhéologie est caractérisée par le frottement interne apparent, la compacité de l'assemblage, les différences de contraintes normales et diverses variables internes, fonctions de deux paramètres de contrôle adimensionnés : le nombre inertiel I et la pression réduite P*, qui compare les forces de confinement à l'adhésion dans les contacts. Notre étude concerne les états homogènes, ce qui exclut les états de cisaillement localisés observés à faible P*, de l'ordre de 0,1. Le coefficient de frottement interne augmente de 0.35 (cas sec) à 0.9 environ pour P*=0.4, tandis que la compacité décroît de 0.59 à 0.52. L'important effet des forces capillaires sur la rhéologie, sensible pour des P* de plusieurs unités, est relié à la texture anisotrope des contacts et des ponts liquides. Lorsque P* décroît, nombre de contacts cohésifs sont maintenus pour des intervalles de déformation de plusieurs unités, survivant aux effets de rotation et de cisaillement de l'écoulement, et forment des amas percolants dans le système entier. Les résultats sont modérément sensibles à la saturation dans le régime pendulaire, mais fortement affectés par l'hystérèse de la conformation des ménisques. En présence de forces visqueuses et non plus capillaires, une version simplifiée de la dynamique stokésienne est adoptée dans laquelle les forces de lubrification entre proches voisins, supposées dominantes, sont les seules interactions hydrodynamiques. La rhéologie est fortement influencée par les contacts intergranulaires directes, qu'autorise la coupure à courte distance de la singulérité de lubrification du fait de la rugosité de surface des particules. Le même état critique que celui des grains secs est approché dans la limite quasi-statique. Nous discutons de lois rhéologiques exprimées en fonction du nombre visqueux qui remplace alors le nombre inertiel, et de la divergence de la viscosité effective à l'approche de la compacité critique en écoulement permanent, ou de la compacité maximale des assemblages aléatoires pour les configurations isotropes désordonnées / We numerically simulate the shear flow of dense assemblies of 3D frictional spherical grains under a fixed normal stress P in steady-state, either in the presence of a small amount of an interstitial liquid, which gives rise to capillary menisci and attractive forces, or in the fully saturated state, when the mechanical properties of suspensions in Stokes flow are controlled by hydrodynamic and contact forces. Dry grain assemblies are used as a reference system for which the rheological properties - in particular the approach to the critical state – are rather well known and can be measured with good accuracy. A non-saturating wetting fluid creates capillary attractive intergranular forces, the effects of which on the rheology are investigated in the pendular state, with shear rates ranging from quasistatic to inertial regimes. The system behavior is characterized by the dependence of internal friction coefficient, solid fraction, normal stress differences and internal state parameters on two dimensionless control parameters: the inertial number, I and the reduced pressure, P*, comparing confining forces to contact tensile strength. We focus on steady homogeneous flows, excluding localized flow patterns which we observe to occur for low P* (of order 0.1). The apparent internal friction coefficient increases to 0.9 in the quasistatic limit for P*=0.4, from its dry value 0.35, while solid fraction decreases from 0.59 to 0.52. We relate the significant effect of capillary forces on the macroscopic behavior of the system, up to P* values of several unities, to fabric anisotropy parameters of contact and distant interactions. As P* decreases, many cohesive contacts are observed to survive the tumbling motion associated to the shear flow, and their average age exceeds the reciprocal shear rate. Corresponding clusters of grains with enduring capillary bonds gather a large proportion of grains and percolate through the sample. The results are shown to be moderately sensitive to saturation within the pendular range, yet rather strongly affected by the hysteretic nature of liquid bridges. In the presence of viscous forces, assuming lubrication effects to dominate the hydrodynamic interactions, we adopt a simplified version of the (overdamped) Stokesian dynamics approach, in which hydrodynamic interactions only couple close neighbours. Rheological properties are strongly influenced by direct intergranular contacts and friction, which are permitted due to a very small distance lubrication cutoff modeling surface asperities. The same critical state as in the dry case is approached in the quasistatic limit. We discuss expressions of rheological laws involving the viscous number instead of the inertial number, and the divergence of effective viscosities in steady flow and in isotropic random suspensions as either the critical state or the random close packing solid fraction are approached