La dégradation géométrique de la voie ferrée sur les Lignes à Grande Vitesse (LGV) est un phénomène qui entraîne des coûts de maintenance importants pour assurer une qualité de circulation et de sécurité. Cette dégradation géométrique est imputable, en partie au tassement de la couche de ballast qui constitue l'une des parties de la voie ferrée.Le tassement du ballast est un phénomène difficile à estimer et à prédire car sa nature granulaire et ses caractéristiques particulières entraînent des variabilités de ses propriétés mécaniques. Ce travail de thèse est consacré à l'étude du comportement mécanique du ballast, pour le développement d'un modèle prédictif du tassement sur voie ferrée à partir de la caractérisation mécanique initiale du matériau (grâce à l'utilisation du pénétromètre léger Panda) et des sollicitations auxquelles la voie est soumise. Grâce à l'étude paramétrique développée sur un banc d'essais à échelle réelle, nous avons constitué et validé un modèle de prédiction qui est basé sur une loi de relaxation logarithmique.La modélisation numérique par éléments discrets en utilisant la méthode de la Dynamique des Contacts permet d'étudier les propriétés mécaniques du matériau ballast à l'échelle des grains. L'analyse des déformations transitoires montre une dépendance claire de la déformation moyenne par rapport à la contrainte appliquée et au rapport d'aspect en raison de l'action du frottement aux frontières. Les fluctuations de ces déformations, sont très importantes et semblent évoluer avec la déformation moyenne. Finalement, la pertinence de la caractérisation du ballast à partir de l'essai Panda a été vérifiée à partir d'une étude paramétrique sur les mécanismes d'enfoncement grâce à l'utilisation de cette approche numérique discrète. / Railway track degradation on high-speed lines is a phenomenon which causes high maintenance costs to ensure quality traffic and safety.This geometric degradation is due in part to the vertical strain (settlement) of the ballast layer, which compose the track.The ballast settlement is a difficult phenomenon to estimate and predict because its granular nature and characteristics of this layer induces variability of its mechanical properties.This work is devoted to the study of the mechanical behavior of ballast, to develop a predictive model of track settlement from the initial mechanical characterization of the material (by means a light penetrometer Panda) and the loading on track. By means several test on a full-scale model of railway track, we established and validated a predictive model based on a logarithmic relaxation law. Numerical modeling by discrete element method using the Dynamic Contacts allows to study the mechanical properties of ballast at grain-scale.Transient deformation analysis shows a clear dependence of the average deformation with respect to the stress and aspect ratio as a result of the frictional feedback at the boundaries. Fluctuations of these deformations are significant and seem to evolve with the average deformation. Finally, the relevance of the ballast characterization using the Panda test has been verified by means a parametric study on mechanisms during the penetration process by discrete element method.

Les matériaux à particules hautement déformables sont des formes complexes de matière avec de nombreuses applications en chimie, pharmacie, cosmétique et agro-alimentaire. L’effet conjugué du désordre et des grandes déformations des particules conduit à des propriétés mécaniques nouvelles par rapport aux matériaux à particules indéformables. En particulier, la compressibilité et la résistance au cisaillement sont contrôlées par une combinaison de réarrangements et de changement de forme des particules. Dans ce travail de thèse, nous avons développé une approche numérique originale pour la simulation de ces systèmes. Pour permettre aux particules de se déformer indéfiniment, nous avons modélisé chaque particule par un agrégat de particules primaires sans frottement qui interagissent par une force d’attraction de type Lennard-Jones et une contrainte de non-interpénétration. La dissipation d’énergie par collisions inélastiques entre les particules primaires confère un caractère plastique aux déformations des particules. Nous avons utilisé ce modèle pour étudier les propriétés de compaction et de cisaillement de ces systèmes. Nos résultats ont permis de mettre en évidence le caractère non-linéaire de la compressibilité lorsque la compacité progressivement augmente au-delà de celles des assemblages de particules indéformables. Sous cisaillement, un état critique est atteint avec une dilatance contrôlée par la pression de confinement. Dans cet état, nous avons exploré les distributions des formes des particules, les textures et les distributions des forces pour différentes valeurs de la pression. Nous avons également comparé la compressibilité simulée par l’approche développée avec celle obtenue par la Méthode de Points Matériels (MPM) en utilisant des particules élastiques. / Soft-particle materials are complex forms of matter that occur in numerous applications in chemical, pharmaceutical, cosmetic and food products. Joint effects of disorder and large particle deformations lead to novel mechanical properties that differ from those of rigid-particle materials. In particular, the compressibility and shear resistance depend on both particle rearrangements and their shape change. In this doctoral work, we developed an original approach for numerical simulation of these systems. To allow the particles to deform without breakage, each particle is modeled as an aggregate of frictionless primary particles interacting via a Lennard-Jones attraction force and impenetrability constraints. Energy dissipation by inelastic collisions between primary particles leads to the plastic nature of particle deformations. This model was used to investigate the compaction and shear behavior of soft-particle systems. We find that the compressibility is strongly nonlinear as the packing fraction increases beyond that of a random close packing of rigid particles. In continuous shearing, a critical state is reached with a dilatancy that depends on the confining pressure. In this state, we investigate the shear resistance, distributions of particle shapes, fabric properties and inter-particle forces as a function of the confining pressure. We also compare our results with those obtained by using the Material Point Method (MPM) with elastic particles.

Cette étude porte sur l'influence de la forme des particules sur le comportementmécanique des matériaux granulaires, et les mécanismes de transfert de charge quis'y développent, notamment dans les cas des conduits enterrés. La géométrie desparticules (polygones de forme convexe ou assemblage de particules de forme concaveconstitués de plusieurs disques superposés et indissociables) a été caractérisée par uncoefficient de forme α. Cette étude est basée sur une approche numérique par élémentsdiscrets. Des simulations numériques de l'essai de compression biaxiale montrent queles caractéristiques macroscopiques ou géométriques de l'échantillon granulaire, telque l'angle de frottement macroscopique, la compacité, ou la nature des bandes decisaillement, dépendent fortement du coefficient de forme α et de la convexité ou nonconvexité des grains.Les mécanismes de transfert de charge au dessus d'un conduit souple ont été étudiésexpérimentalement (rouleaux bidimensionnels en condition de déformation plane) etnumériquement (MED). Les expérimentations réalisées montrent que la présence duconduit à peu d'in_uence sur le comportement macroscopique de l'assemblée granulairelors d'une sollicitation biaxiale. Les résultats du modèle numèrique convergentavec les rèsultats expérimentaux et mettent en évidence la présence des mécanismes detransfert de charge au dessus du conduit dont les intensités dépendent du coefficientde forme α. / This study was devoted to the in_uence of grain shape on the mechanical behaviourof granular materials and its e_ect on load transfer over underground pipes. Shapeof convex polygons and concave clumps of discs was generalised with a geometricalparameter α. In the study a Discrete Element Modeling (DEM) approach was used.Biaxial compression of granular assemblies revealed that mechanical and geometricalproperties like porosity, macroscopic friction or shear localisation depends both on αand on grain (non-)convexity.The intergranular load transfer over a _exible pipe was studied both experimentally(2D rods in a plain strain apparatus) and numerically (DEM). The experimentsshowed that the pipe has no signi_cant impact on the macroscopic behaviour of theassembly. The numerical model complied with the experiments and revealed that thearching e_ect in a sheared granular medium exists above the pipe and is magni_edwith the increase of α of the grains.

Nous nous intéressons à l'écoulement de milieux granulaires cohésifs soumis à des vibrations horizontales en vue d'une application au remplissage des moules de presse. L'effet des vibrations est analysé pour un grain ou un milieu granulaire déposé sur un plan incliné et un milieu granulaire contenu dans une conduite verticale de section constante. Pour un grain polyédrique, trois régimes sont identifiés : stick, stick-slip et slip. Un modèle basé sur le bilan des forces appliquées à un grain a permis de préciser les limites de ces trois domaines. La vitesse des grains tend vers une asymptote lorsque l'inclinaison est inférieure à 0,7 fois l'angle de frottement entre les grains et le plan. Par ailleurs, l'inertie du grain conduit à une réduction de l'amplitude de son oscillation lorsque l'accélération augmente. La comparaison modèle-expérience permet d'estimer le coefficient de frottement. Pour un milieu granulaire de faible épaisseur, nous retrouvons les mêmes tendances. Les milieux granulaires épais sont le lieu de forts cisaillements provoquant leur étalement très rapide. En conduite verticale, une compétition s'établit entre les déplacements horizontal et vertical des particules. Selon l'inertie, les vibrations appliquées pourront conduire soit à une dilatance soit à une compaction. Cette dernière provient du déplacement de la conduite et du milieu granulaire en sens opposé pendant une partie d'une période de vibration. Si le phénomène de compaction atteint le centre de la conduite, il se forme des arches qui bloquent l'écoulement. Les observations faites à l'échelle mésoscopique permettent d'interpréter les vitesses d'écoulement mesurées à l'échelle macroscopique. / With the target to improve the feeding step during the moulding process, we study the flow of cohesive granular matter submitted to vibration. Vibration effects are first analyzed on the sliding motion of a single particle on inclined plane and also on granular matter sample deposited on an inclined plane or in a vertical funnel. For a single particle, three regimes are identified: stick, stick-slip and slip regimes. A simple model based on the movement equation of one grain allows to determinate the limits between the three regimes. The grain velocity reaches asymptote when the plane inclination is below 0.7 times the friction angle between the grain and the substrate. Otherwise, when the acceleration increases, the grain inertia leads to reduce the transverse amplitude of the grain oscillation motion. The comparison between experimental and numerical results allows to estimate a value for the friction coefficient parameter. For a granular sample, we find the same tendencies. Thick granular layers are submitted to high shear, causing a fast spreading. In vertical funnel, there is a competition between vertical and horizontal motions. Depending on granular inertia, submitted vibrations could lead to dilatancy due to the shear or lead to compaction. This compaction is due to an opposite motion direction between the sample and the funnel during one part of the period vibration. If the compaction reaches the center of the funnel, arches are formed and jammed the flow. Observations realized at mesoscopic scale allow to interpret flow velocities measured at macroscopic scale.

Le travail présenté dans cette thèse vise à apporter une compréhension fine des couplages hydromécaniques et plus particulièrement du phénomène d'effondrement capillaire. L'utilisation de matériaux granulaires modèles observés à petites échelles grâce au microtomographe à rayons X constitue le point de départ de cette étude. Les résultats portent dans un premier temps sur la capacité à reproduire et observer le phénomène d'effondrement capillaire à l'échelle de l'éprouvette, appelée aussi échelle macroscopique. Une des premières problématiques importantes de ce projet a donc porté sur la reconstitution d'éprouvettes suffisamment lâches permettant d'observer le phénomène effondrement. Le potentiel d'effondrement, alors mesuré au simple œdomètre sur ces éprouvettes reconstituées, décroit avec la charge mécanique appliquée. En contrôlant le volume d'eau au cours de l'imbibition, les éprouvettes sont sujettes à une succession d'incréments de déformation rapides et ceci quelle que soit la dynamique de mouillage. En parallèle, l'utilisation de la technique d'analyse non destructive par microtomographie aux rayons X nous a permis de quantifier à une échelle mésoscopique (échelle de quelques grains), l'évolution de matériaux granulaires sous chargement hydrique et mécanique. Cela a été possible par l'utilisation de la technique de corrélation d'images numériques (DIC) associée à des outils d'analyse adaptés aux matériaux non saturés. Ces outils révèlent que l'hétérogénéité locale des déformations et des teneurs en eau augmente au cours de l'imbibition principalement à cause de l'effet de la gravité sur la répartition de l'eau dans l'éprouvette. L'hétérogénéité initiale des éprouvettes induite par la méthode de préparation est aussi identifiée. Cependant, malgré la présence de ces hétérogénéités, nous montrons que les mesures réalisées à l'échelle de l'éprouvette restent représentatives du comportement mésoscopique du matériau au cours de l'effondrement / This thesis deals with hydromechanical couplings in soils and more particularly focusing on the capillary collapse phenomenon. Although thoroughly studied at macroscopic scale in fine soils, the mechanisms governing capillary collapse in granular soils have been less observed.At first, the results concern the capacity to reproduce and to observe the phenomenon of capillary collapse at the scale of the specimen, called also the macroscopic scale.Thus, one of the first important objectives of this project concerned the preparation method of the specimens, which have to be loose enough to observe this phenomenon. The collapse potential, determined here by the one-dimensional collapse test, decreases with the mechanical stress. Moreover, when the quantity of water is controlled, incremental strains of specimens arise whatever the imbibition dynamic.To complement existing data, an experimental program was carried out using X-ray computed tomography (X-ray CT), a relevant tool to obtain 3D informations at grain scale. In combination with Volumetric Digital Image Correlation (V-DIC), X-ray CT also makes it possible to get information at a mesoscopic scale, at which the use of continuum mechanics still remains relevant. In this work, both microscopic and mesoscopic analyses have been carried out on a partially saturated sand subjected to gradual controlled imbibition in an oedometric cell.Local heterogeneities of strains and water contents appear during the imbibition, which are principally explained by the effect of gravity and by an inhomogeneity of the initial density.However, in spite of the existence of those localisations, by comparing the analyses at the mesoscopic and at the macroscopic scale, we show that measurements performed at the scale of the specimen are representative of the local behaviour of the material

Cette thèse expérimentale porte sur les ondes ultrasonores se propageant dans un milieu granulaire sec, sous contrainte. Dans une première partie, on s'intéresse à la propagation des ondes cohérentes de grande longueur d'onde. La mesure des vitesses de ces ondes en fonction de la contrainte et de la fréquence, montre que les propriétés élastiques du milieu résultent du contact de type Hertz-Mindlin, et que le milieu est non dispersif. La mesure d'atténuation indique un mécanisme dominant du type diffraction de Rayleigh. La propagation d'une onde de Biot à travers l'air dans un empilement sous gravité a été aussi mise en évidence. Dans une seconde partie, on s'est intéressé à la diffusion multiple d'onde acoustique. Nous avons déduit le coefficient de diffusion et le temps inélastique d'absorption en fonction de la contrainte et de la fréquence d'excitation. Le libre parcours moyen de transport a été déterminé en fonction du rapport de la longueur d'onde sur la taille de grains, montrant deux régimes distincts à basse et à haute fréquence. La dernière partie, porte sur le comportement non linéaire hystérétique d'un milieu granulaire, par méthode de résonance. Les modes fondamentaux d'ondes cohérentes ont été étudiés en fonction de l'amplitude d'excitation et de la contrainte de confinement. Un ramollissement élastique de 10%, pour l'onde longitudinale, et de 20%, pour l'onde transversale, sont observés, en bon accord avec un modèle de milieu effectif basé sur la théorie de Mindlin. Les effets d'interaction irréversible son-matière ont été examinés, tels que la compaction intermittente et le processus de cicatrisation. Enfin, on observe que l'atténuation de l'onde transversale dépend conjointement de l'amplitude d'excitation et de la fréquence / This experimental work deals with the ultrasound propagation in a dry granular medium under stress. The first chapter of the thesis concerns the coherent wave propagation in the long-wavelength limit. The measurement of sound velocity as a function of the applied stress and wave frequency confirm that the nonlinear elasticity of the medium originates from the Hertz-Mindlin contact and show that the medium is free of velocity-dispersion. On the other hand, the amplitude measurement reveals that the dominant mechanism of attenuation is due to the Rayleigh-like scattering. Additionally, the airborne Biot wave is observed and investigated. The second chapter is devoted to the diffusively scattered ultrasound. We particularly determine the diffusion coefficient and the inelastic absorption time as a function of stress and frequency. The mean transport free path was investigated versus the ratio of the wavelength to the grain size, which shows two distinct behaviours at low and high frequencies, respectively. The third and last chapter is focused on the hysteretic nonlinear behaviour studied by the resonance method. The fundamental modes of coherent waves are investigated as a function of the applied stress and excitation amplitude. The softening of elastic modulus is observed and can rises to 10% for the longitudinal wave and 20% for the transversal wave, which are in good agreement with the model prediction based on the Mindlin friction theory. Furthermore, the irreversible sound-matter interaction is examined; the intermittent compaction and the slow healing process are observed. Finally, the attenuation of transversal waves is shown to depend simultaneously on wave amplitude and frequency

Depuis les dernières années, le gouvernement provincial a investi des montants record dans le réseau routier québécois, notamment pour la construction de nouvelles routes et l’entretien des chaussées. En plus d’être doté d’une grande superficie territoriale, le Québec connaît un climat rigoureux qui complexifie la conception des chaussées. Les changements de saison entraînent des répercussions néfastes quant à la performance des chaussées, dont l’accumulation de déformations permanentes (ɛp) dans les structures granulaires. Bien que la déformation permanente puisse se développer dans l’ensemble de la structure routière, seule la couche de fondation a été considérée dans cette étude. Afin d’améliorer les connaissances sur la déformation permanente à long terme des matériaux granulaires sollicités aux changements saisonniers nordiques, quatre sources granulaires de fondation (MG-20) typique du Québec (grauwacke, gneiss granitique, calcaire et basalte) ont été soumises à des essais de déformation permanente en laboratoire. Ces essais ont été réalisés à l’aide de deux appareils : la presse hydraulique à chargement répété et le simulateur à charge roulante. Les résultats obtenus ont été modélisés selon le modèle de Desden. Ce modèle mathématique est composé de quatre paramètres qui décrivent la phase de post compaction (a), le taux de déformation permanente (b) et la phase de rupture (c et d). Les principaux résultats de déformation permanente montrent que la composition minéralogique, la granulométrie, le degré de saturation (%SR) et la variation des conditions environnementales influencent le comportement plastique d’un matériau granulaire (MG). En effet, les résultats montrent que la composition minéralogique peut provoquer une augmentation de 62% du paramètre b pour une même teneur en eau et que de passer d’une granulométrie grossière à une granulométrie plus sableuse peut entrainer une augmentation moyenne de 30,75% du paramètre b. De plus, les résultats montrent qu’une variation du degré de saturation peut entraîner une augmentation de la déformation permanente de 147%. Les essais réalisés en condition de gel et de dégel ont également montré que le cycle de gel et dégel a un effet direct sur le paramètre de post compaction. Un essai au simulateur à charge roulante a été effectué dans le but de valider les résultats obtenus en condition triaxiale. Les résultats montrent que le simulateur provoque à l’échantillon une déformation permanente plus importante que celle enregistrée en condition triaxiale. Cependant, la phase de post compaction reste très similaire, peu importe l’instrumentation utilisée. Cette étude constitue un premier essai au Canada quant à l’application d’une charge roulante mobile en contexte de gel et dégel. Deux relations de corrélation entre les paramètres a et b du modèle de Dresden et les propriétés de base des matériaux granulaires ont été développées. La première relation développée est une relation de type linéaire qui relie le paramètre de post compaction (a) et la variable complexe nfopt et ce. La seconde relation développée est également de type linéaire et lie le paramètre b et le degré de saturation, le pourcentage de fracturation, le coefficient d’uniformité et la masse volumique sèche du matériau. L’intégration des deux relations de corrélation dans la relation mathématique de Sweere a permis de développer une relation de prédiction de déformation permanente à long terme.

Le travail présenté dans ce mémoire a pour objectif d'aboutir à une meilleure caractérisation expérimentale du comportement non drainé du sable d'Hostun à l'échelle locale, en relation avec les différentes méthodes de reconstitution utilisées. On vise en particulier à mettre en évidence la façon dont la microstructure influence les propriétés des sols granulaires, notamment vis-à-vis de l'instabilité de liquéfaction. Le point clé de ce travail consiste en l'utilisation de capteurs piézoélectriques pour mesurer les vitesses d'onde de cisaillement au sein d'éprouvettes triaxiales reconstituées par trois techniques différentes à savoir le damage humide, la pluviation à sec et la pluviation sous eau. Les résultats obtenus mettent en évidence l'influence considérable du mode de fabrication sur le comportement mécanique du sable. Ils montrent que pour une même densité initiale, et à même pression de confinement, le comportement mécanique observé est très différent selon les trois modes de préparation. Les valeurs de vitesses d'onde de cisaillement enregistrées sont différentes d'un mode de déposition à un autre. Concernant l'anisotropie, le matériau préparé par damage humide présente un certain degré d'anisotropie inhérente, qui est toutefois minimal. En revanche, la pluviation à sec et la pluviation sous eau induisent des structures initiales plus isotropes. / This work aims to achieve a better experimental characterization of the undrained behavior of the Hostun HN31 sand at the local scale, depending on the different procedures used for sample reconstitution. We particularly highlight the influence of the microstructure on the mechanical behavior of sand especially in relation to liquefaction instability. The key point of this work consists of the use of piezoelectric transducers to measure shear wave velocities within triaxial sand specimens. These specimens have been prepared by different methods, namely: moist tamping, dry pluviation and water pluviation. The obtained results highlight the considerable influence of the deposition mode on the mechanical behavior of sand. They show that for the same initial density, and at the same confining pressure, the mechanical behavior is very different from the three preparation methods. Furthermore, the values of shear wave velocity are different from one deposition mode to another. Concerning the anisotropy, the material prepared by moist tamping presents an inherent degree of anisotropy, which it is however minimal. On the other hand, dry pluviation and water pluviation induce more isotropic initial structures. Moreover, the initial anisotropy is less important than the anisotropy induced by loading. These experimental results indicate that the initial state of a granular medium can not be described by only couple of scalar variables of initial density and confining pressure: a variable that can reflect the initial microstructure and an evolution law of this variable are necessary to model the experimental observations.

Cette thèse traite de la modélisation et des méthodes numériques pour la simulation d'écoulements de fluides non-Newtoniens, et particulièrement, de matériaux granulaires. Une application de ce travail concerne les poudres de couverture utilisées pour protéger thermiquement le métal de l'air dans le procédé de coulée en source d'alliages métalliques. Ces poudres sont conditionnées dans des sacs disposés dans la lingotière, qui brûlent suite aux fortes chaleurs engendrées, et permettant son écoulement sur la surface du métal. Ainsi, la simulation numérique apparaît comme un puissant outil pour l'optimisation du procédé, et notamment, de l'étalement de ces poudres.Dans ce travail, une formulation éléments finis a été proposée pour modéliser l'écoulement multiphasique des matériaux granulaires dans un formalisme de la mécanique des milieux continus. Les équations associées sont résolues via des schémas numériques stabilisés, couplés avec la méthode Level-Set pour capturer et suivre le profil du matériau granulaire au cours de la simulation. Dans un premier temps, les outils numériques ont été testés sur des cas d'écoulements de fluides de Bingham, où les fortes non-linéarités sont traitées par une méthode de régularisation. Puis la formulation est étendue aux écoulements de granulaires secs, dont le comportement piezzo-dépendent est traduit par la loi mu(I). Le modèle a été validé sur des cas d'effondrement de colonnes de grains, et une étude de sensibilité aux conditions aux limites et constantes physiques du modèle est proposée.Enfin, des cas industriels de chutes de poudres sur substrats solide et métal fondu ont été menés, amenant à des premières pistes pour l'optimisation du procédé de coulée en lingotière. / This thesis is devoted to the modeling and numerical methods for the simulation of non-Newtonian flows, and focuses particularly on granular materials flows. This work is applied to molten powders aiming to ensure metal thermal protection from the air in ingot casting process of metallic alloys. These powders are conditionned into bags disposed into the mold, which burn due to high temperatures, and allowing the powder spreading onto the metal surface. Thus, numerical simulation appears as a powerful tool for the process optimization, and especially, for the powder spreading.In this work, a finite element formulation has been proposed for the modeling of granular multiphase flows, by a continuum approach. The associated equations are solved using stabilized numerical schemes, coupled with the Level-Set method to capture and follow the granular profile during the simulation. First, the numerical tools have been implemented for Bingham flows, by using regularization a method. Then, the formulation was extended to dry granular flows, by the use of the mu(I) rheology constitutive model for describing its pressure-dependent behavior. The model has been validated on granular collapses, and a sensitivity analysis to boundary conditions and physical constants has been proposed.Finally, industrial cases of powder chutes ontoboth solid and liquid metla substrates have been conducted, leading to preliminary solutions for the optimization of ingot casting process.

Les méthodes par éléments discrets constituent une alternative par rapport aux méthodes<br />par éléments finis pour la modélisation du comportement d'un milieu divisé : sol, milieux<br />granulaires, etc . Dans ce mémoire l'objectif est d'étudier un matériau granulaire soumis à des sollicitations<br />cycliques : le ballast des voies ferrées soumis au passage des trains. Nous présentons à la<br />suite d'une étude bibliographique sur le comportement du ballast, la méthode de résolution NonSmooth<br />Contact Dynamics, et nous exposons l'ensemble des développements nécessaires dans le<br />cas bidimensionnel et tridimensionnel : un algorithme de détection entre polygones et entre polyèdres<br />convexes, le paramétrage, la résolution du problème de contact frottant. Nous déterminons<br />ensuite les paramètres optimaux pour des simulations de chargement cyclique. Dans une dernière<br />partie, nous présentons ensuite un ensemble de résultats montrant l'aptitude des méthodes par éléments<br />discrets à décrire l'évolution d'un massif granulaire sous chargement cyclique. L'analyse<br />d'une comparaison entre une expérience et la simulation met évidence le comportement particulier<br />de ce type de système, qui s'assimile à une couche mince granulaire où des micro-structures particulières<br />guident la réponse mécanique du milieu. Dans le cas tridimensionnel, des comportements<br />mécaniques connus ont été retrouvés, à titre d'exemple, en analysant l'évolution de structures<br />maçonnées. Nous présentons une étude de la résistance latérale de la voie ballastée, en considérant<br />des formes de grains issues de la digitalisation de grains réels, pour laquelle on retrouve un<br />comportement identifié expérimentalement.