Dans cette étude, l’intérêt est porté sur l’effet des conditions opératoires sur l’enrobage en voie sèche de grosses particules « hôtes » par de fines particules « invitées » et aussi sur la modélisation de cet enrobage selon la méthode par éléments discrets (DEM) afin de mieux comprendre les phénomènes mis en jeu. Dans ce travail, les matériaux choisis comme particules hôtes sont les Suglets® et les particules invitées sont en stéarate de magnésium (StMg). Ces deux éléments sont introduits dans un mélangeur à fort taux de cisaillement, le Cyclomix. Les propriétés du produit final, comme la coulabilité, la mouillabilité et le degré d’avancement de l’enrobage, ont été caractérisés. La variation des propriétés est étudiée en fonction de la durée de traitement dans le Cyclomix pour diverses vitesses de rotation, taux de remplissage et rapport de taille de particules hôte et invitée. La coulabilité a été améliorée avec la durée du mélange ou la vitesse de rotationLe degré d’avancement présente une tendance semblable pour différentes conditions opératoires. Sa variation est représentée par une loi exponentielle empirique en fonction du temps de mélange, paramétrée par une constante ajustable. Cette constante permet d’estimer l’efficacité de l’enrobage. La simulation des mouvements de particules dans le mélangeur par DEM a permis d’obtenir des informations sur la position, la vitesse des particules, et d’autres paramètres énergétiques. Les champs de vitesse réelle ou numérique liés aux mouvements de particules, analysés par PIV (Particule Image Velocimetry), sont analogues. La constante d’enrobage dépend de la vitesse de rotation simulée et peut être ainsi prédite par DEM. / Investigations of the effect of the operating conditions on the dry coating and the motion of particles in a high shear mixer by the DEM simulation have been done. Big Suglets® and small Magnesium Stearate materials have been chosen for the dry coating process. The treatment has been carried out in a Cyclomix, a high shear mixer, at different mixing time, rotational speeds, filling ratio of the samples and particle size ratio to observe the effect of the operating conditions on the dry coating. A conversion ratio has been introduced to quantify the degree of coating and measured. The flowability has been improved when increasing the operation time in the mixer. Higher speeds of rotation can improve the flowability more rapidly. However, the flowability doesn’t seem to be sensitive to the filling ratio. The curves of conversion ratio versus time exhibit the same kind of variations for different operating conditions. It is possible to approximate this tendency by an exponential function in which a characteristic parameter “coating rate constant” is introduced, linked to the efficiency of the dry coating process, since the conversion ratio shows a linear relation with a flowability index and the wettability angle. Simulation of the particle motion in the mixer has been carried out by a Discrete Element Method: different parameters characterizing the location, the velocities, and the attachment of the particles have been derived. The velocity fields of the real and simulated particle motions, compared by Particle Image Velocity (PIV), are quite similar, validating the DEM method and allowing the prediction of the coating process.

L’augmentation de la vitesse ferroviaire impose l’amélioration des performances de freinage liés aux matériaux de friction utilisés dans les freins à disques et à semelles. Cependant, l'usure reste un point limitant en terme de performance mais aussi de sécurité. Elle a pour conséquence de détériorer les matériaux (aspects mécano-thermiques) mais aussi de contaminer le contact entre roue et rail par des particules de semelles de frein pouvant isoler électriquement celui-ci (aspects mécano-électriques). Le procédé permettant de localiser les trains, ayant pour principe la conduction électrique entre roue et rail, est alors compromis.Ces problématiques sont gouvernées par les différents contacts (roue/rail, roue/semelle, ...) qui s'inscrivent dans le concept de triplet tribologique. Un triplet est composé des corps en contact (premiers corps), de l'interface (troisième corps) et du mécanisme qui les maintient en contact. Afin de comprendre les phénomènes se produisant au sein des contacts, une analyse multi-physique locale et un découplage des paramètres (mécaniques, thermiques, électriques) s'avèrent nécessaires. Toutefois, cela est difficile à établir expé- rimentalement, une modélisation numérique par éléments discrets est alors pertinente pour palier à ces besoins.L'écoulement dynamique d'un troisième corps constitué de particules de freinage et la dégradation d'un matériau de freinage réel (suite à des méthodes d'imageries) sont modélisés. Les résultats obtenus ont permis d'appréhender les compétitions existant entre paramètres physiques d'une part (mécanique, électrique, thermique) et paramètres d'éch-elles d'autre part (local, global). Ces compétitions tendent à minimiser la résistance élec- trique lorsque les constituants sont en proportion adéquat (meilleure conduction électrique entre roue et rail) mais aussi à équilibrer les débits du circuit tribologique, ce qui conduit à protéger les matériaux (usure contrôlée). En outre, les analyses numériques ont permis de mettre en évidence de nouveaux paramètres liés à ``l'étalement'' du troisième corps et d'appréhender les points clés permettant de se rapprocher d'un contact réel. / The increase of the rail speed imposes the improvement of the braking performances related to the friction materials used in brake pad and brakes shoes. However, wear remains a limiting point in terms of performance but also safety. It has the effect of damaging the materials (mechano-thermal aspects) but also to contaminate the contact between wheel and rail by particles of brake shoes that can electrically isolate it (mechanical-electrical aspects). The method for locating the trains, whose principle is the electrical conduction between wheel and rail, is then compromised.These issues are governed by the different contacts (wheel/rail, wheel/shoes, ...) that fit into the concept of tribological triplets. These are composed of the bodies in contact (first bodies), the interface (third body) and the mechanism that keeps them in contact. In order to understand the phenomena occurring within the contacts, a local multi-physics analysis and a decoupling of the parameters (mechanical, thermal, electrical) are necessary. However, this is difficult to establish experimentally, numerical modeling by discrete elements method is then relevant for these needs.The dynamic flow of a third body consisting of braking particles and the degradation of a real braking material (following imaging methods) are modeled. The results obtained made it possible to understand the competitions existing between physical parameters on the one hand (mechanical, electrical, thermal) and scale parameters on the other hand (local, global). These competitions tend to minimize the electrical resistance when the constituents are in adequate proportion (better electrical conduction between wheel and rail) but also to balance the flows of the tribological track, which leads to protect the materials (controlled wear). In addition, the numerical analyzes made it possible to highlight new parameters related to "spreading" of the third body and to apprehend the key points making it possible to approach a real contact.

Neste trabalho, apresenta-se um estudo a respeito de problemas de interação Fluido-Partícula e dos métodos e formulações usados para resolvê-los. Além disso, propõe-se um novo método para resolver problemas acoplados de mecânica dos fluidos e mecânica das partículas. A ideia é baseada em trabalhos anteriores de (Gomes e Pimenta 2015) e (Campello 2016), e o objetivo é desenvolver um modelo computacional adequado e eficiente para simular problemas envolvendo fluidos carregados com partículas sólidas. O problema de fluido é resolvido por uma abordagem Euleriana de elementos finitos, usando elementos mistos locais de velocidade-pressão, os quais satisfazem a condição de compatibilidade LBB. O sistema de equações não linear obtido é resolvido iterativamente pelo método de Newton-Raphson. Uma característica importante é que a malha de fluidos permanece fixa durante a passagem do fluxo, assim como nas abordagens Eulerianas clássicas. O problema das partículas, por sua vez, é resolvido com uma abordagem Lagrangiana de elementos discretos, em que os contatos de partícula com partícula e partícula com paredes (limites fixos) são livremente permitidos e resolvidos. A influência do fluido no movimento das partículas é representada por meio de forças e momentos, que são calculados a partir do escoamento e impostos às partículas de maneira acoplada, iterativa e explícita. As interfaces entre o fluido e as partículas são tratadas por meio da técnica de fronteiras imersas, onde as condições de contorno do fluido nos contatos com partículas são impostas através da interpolação de funções descontínuas e constantes de multiplicadores de Lagrange ao longo das interfaces. É adotado um método explícito, interativo e escalonado para conseguir a convergência dentro de cada passo de tempo do problema. Para ilustrar o potencial do método proposto, são apresentadas simulações de escoamentos carregados de partículas. / In this work, a study is presented on Fluid-Particle interaction problems and the methods and formulations used to solve them. In addition, is proposed a new method for solving coupled problems of fluid mechanics and particle mechanics. The idea is based on previous works by (Gomes e Pimenta 2015) and (Campello 2016), and the goal is to develop an efficient computational model suited to simulate problems involving flowing fluid media laden with solid particles. The fluid problem is resolved by an Eulerian finite element approach using local element velocity-pressure pairs satisfying the LBB compatibility condition, with the resulting nonlinear system of equations being iteratively solved by the Newton-Raphson method. As an important feature, the fluid mesh remains fixed during the flow, just as in classical Eulerian approaches. The particles´ problem, in turn, is resolved in a Lagrangian discrete element approach, wherein both particle-to-particle and particle-to-wall (fixed boundaries) contacts are freely permitted and resolved. The influence of the fluid on the motion of the particles is represented by means of forces and moments, which are computed from the fluid flow and imposed on the particles in a coupled, iterative and explicit way. The fluid-particles´ interfaces are treated by means of immersed boundary technique, in which the fluid interface conditions with the (nonmatching) particles´ boundaries are imposed through discontinuous piecewise constant Lagrange multipliers interpolating functions along the interfaces. An explicit, staggered and interactive scheme is adopted to achieve convergence within each time step of the problem. In order to illustrate the potentialities of the proposed scheme, particle-laden fluid flow simulations are presented.

A presente pesquisa tem como objetivo realizar um estudo teórico e desenvolver simulações computacionais envolvendo a dinâmica de sistemas gás-sólido. O foco principal do trabalho é a modelagem do particulado através da análise das forças de contato entre partículas de materiais granulares utilizando modelos contínuos baseados na mecânica dos solos e na teoria cinética dos escoamentos granulares (sistemas grandes com muitas partículas, formulação Euleriana - Volumes Finitos) e modelos discretos baseados nas características físicas dos materiais (sistemas intermediários e número limitado de partículas, formulação Lagrangeana - Método dos Elementos Discretos). Investigam-se os modelos existentes na literatura com intuito de melhorar os modelos contínuos e discretos baseados na interação entre as partículas que caracterizam a dinâmica do particulado em sistemas gás-sólido. Propõe-se uma nova abordagem para a determinação do coeficiente de rigidez da mola baseada em uma equivalência entre os modelos lineares e não-lineares. Utiliza-se o código fonte MFIX para realizar simulações computacionais da dinâmica de sistemas gás-sólido, analisando o processo de fluidização, mistura e segregação de partículas, influência das correlações de arrasto, e análise das forças de contato entre as partículas através do novo método para a determinação do coeficiente de rigidez da mola . Os resultados obtidos são comparados com dados numéricos e experimentais da literatura. / The purpose of the present study is to perform a theoretical study and develop numerical simulations involving dynamic in gas-solid systems. The focus of the work is the modeling of particulate matter using continuous models based on soil mechanics and the kinetic theory of granular flows (large systems with many particles, Eulerian formulation - Finite Volume) and discrete models based on physical characteristics of the particles (intermediate systems and limited number of particles, Lagrangian formulation - Discrete Element Method). It is proposed a new approach to determine the normal spring stiffness coefficient of the linear model through the numerical solution for the overlap between particles in non-linear models. The linear spring stiffness is determined using an equivalence between the linear and the non-linear models. It is used the MFIX computational code to perform numerical simulations of the dynamics of gas-solid systems. It is analyzed the processes of fluidization, mixing and particle segregation and the influence of drag correlations. The proposed approach for normal spring stiffness coefficient is applied in the numerical simulations of two problems: single freely falling particle and bubbling fluidized bed. The results were compared with numerical and experimental data from literature.

Cette thèse porte sur l'étude numérique des écoulements fluide-particules rencontrés dans l'industrie. Ces travaux se situent dans le cadre de la compréhension des phénomènes qui se déroulent dans des tambours tournants et réacteurs à lit fixe en présence de particules de forme non convexe. En effet, la forme des particules influence de manière importante la dynamique de ces milieux. A cet effet, nous nous sommes servis de la plateforme numérique parallèle Grans3D pour la dynamique des milieux granulaires et PeliGRIFF pour les écoulements multiphasiques. Dans la première partie de cette thèse, nous avons développé une nouvelle stratégie numérique qui permet de prendre en compte des particules de forme arbitrairement non convexe dans le solveur Grains3D. Elle consiste à décomposer une forme non convexe en plusieurs formes convexes quelconques. Nous avons nommé cette méthode « glued-convex ». Le modèle a été validé avec succès sur des résultats théoriques et expérimentaux de tambours tournants en présence de particules en forme de croix. Nous avons aussi utilisé le modèle pour simuler le chargement de réacteurs à lits fixes puis des lois de corrélation sur les taux de vide ont été déduites de nos résultats numériques. Dans ces travaux, nous avons aussi testé les performances parallèles de nos outils sur des simulations numériques à grande échelle de divers systèmes de particules convexes. La deuxième partie de cette thèse a été consacrée à l'extension du solveur PeliGRIFF à pouvoir prendre en compte la présence de particules multilobées (non convexes) dans des écoulements monophasiques. Une approche du type Simulation Numérique Directe, basée sur les Multiplicateurs de Lagrange Distribués / Domaine Fictif (DLM/FD), a alors été adoptée pour résoudre l'écoulement autour des particules. Une série d'études de convergence spatiale a été faite basée sur diverses configurations et divers régimes. Enfin, ces outils ont été utilisés pour simuler des écoulements au travers de lits fixes de particules de forme multi-lobée dans le but d'étudier l'influence de la forme des particules sur l'hydrodynamique dans ces lits. Les résultats ont montré une consistance avec les résultats expérimentaux disponibles dans la littérature. / Non convex granular media are involved in many industrial processes as, e.g., particle calcination/drying in rotating drums or solid catalyst particles in chemical reactors. In the case of optimizing the shape of catalysts, the experimental discrimination of new shapes based on packing density and pressure drop proved to be difficult due to the limited control of size distribution and loading procedure. There is therefore a strong interest in developing numerical tools to predict the dynamics of granular media made of particles of arbitrary shape and to simulate the flow of a fluid (either liquid or gas) around these particles. Non-convex particles are even more challenging than convex particles due to the potential multiplicity of contact points between two solid bodies. In this work, we implement new numerical strategies in our home made high-fidelity parallel numerical tools: Grains3D for granular dynamics of solid particles and PeliGRIFF for reactive fluid/solid flows. The first part of this work consists in extending the modelling capabilities of Grains3D from convex to non-convex particles based on the decomposition of a non-convex shape into a set of convex particles. We validate our numerical model with existing analytical solutions and experimental data on a rotating drum filled with 2D cross particle shapes. We also use Grains3D to study the loading of semi-periodic small size reactors with trilobic and quadralobic particles. The second part of this work consists in extending the modelling capabilities of PeliGRIFF to handle poly-lobed (and hence non-convex) particles. Our Particle Resolved Simulation (PRS) method is based on a Distributed Lagrange Multiplier / Fictitious Domain (DLM/FD) formulation combined with a Finite Volume / Staggered Grid (FV/SG) discretization scheme. Due to the lack of analytical solutions and experimental data, we assess the accuracy of our PRS method by examining the space convergence of the computed solution in assorted flow configurations such as the flow through a periodic array of poly-lobed particles and the flow in a small size packed bed reactor. Our simulation results are overall consistent with previous experimental work.

L’objectif de ce travail est d’élaborer et de caractériser pour des applications biomédicalesun matériau auto-enchevêtré à base d’une seule fibre d’alliage à mémoire deforme de type Nickel-Titane. Nous avons optimisé un procédé de fabrication consistantà enchevêtrer et figer un ressort par des traitements thermiques. Les échantillonsont été caractérisés en compression et traction, avec suivi par caméra optique ettomographie in-situ. Les structures obtenues sont homogènes, isotropes, superélastiquesà température ambiante jusqu’à des déformations d’au moins 30%, et peuventdevenir ferroélastiques avec un effet mémoire d’au moins 16% par un traitement thermiqueadditionnel. Leur comportement en compression est consolidant puis dilatantet en traction, légèrement auxétique. Une comparaison avec des milieux similairesconstitués de fils ductiles et viscoélastiques, ainsi qu’avec des simulations par élémentsdiscrets sur des milieux élastiques sans frottement, montre que les propriétésmécaniques des structures auto-enchevêtrées sont contrôlées par leur architecturesingulière, à mi-chemin entre milieux continus et discrets. / The aim of this work is to process and characterize for biomedical applications,self-entangled structures made of a single NiTi shape memory fiber. We have optimizeda processing route consisting in entangling and shape-setting a spring bythermomechanical treatments. The samples were characterized in compression andtension, using optical and x-ray tomographic observations. The structures thus obtainedare homogeneous, isotropic, superelastic at room temperature up to strains ofat least 30%, and can become ferroelastic with a shape memory effect up to at least16% strain by an additional heat treatment. The mechanical behavior in compressionis first consolidating and then dilating, while in tension, the samples are slightlyauxetic. A comparison with similar media made of ductile and viscoelastic fibers,as well as with discrete element simulations on friction-free elastic fibers, show thatthe mechanical properties of these self-entangled structures are controlled by theirunique architecture, in-between continuous and discrete media.

L’objectif de ce travail est d’élaborer et de caractériser pour des applications biomédicalesun matériau auto-enchevêtré à base d’une seule fibre d’alliage à mémoire deforme de type Nickel-Titane. Nous avons optimisé un procédé de fabrication consistantà enchevêtrer et figer un ressort par des traitements thermiques. Les échantillonsont été caractérisés en compression et traction, avec suivi par caméra optique ettomographie in-situ. Les structures obtenues sont homogènes, isotropes, superélastiquesà température ambiante jusqu’à des déformations d’au moins 30%, et peuventdevenir ferroélastiques avec un effet mémoire d’au moins 16% par un traitement thermiqueadditionnel. Leur comportement en compression est consolidant puis dilatantet en traction, légèrement auxétique. Une comparaison avec des milieux similairesconstitués de fils ductiles et viscoélastiques, ainsi qu’avec des simulations par élémentsdiscrets sur des milieux élastiques sans frottement, montre que les propriétésmécaniques des structures auto-enchevêtrées sont contrôlées par leur architecturesingulière, à mi-chemin entre milieux continus et discrets. / The aim of this work is to process and characterize for biomedical applications,self-entangled structures made of a single NiTi shape memory fiber. We have optimizeda processing route consisting in entangling and shape-setting a spring bythermomechanical treatments. The samples were characterized in compression andtension, using optical and x-ray tomographic observations. The structures thus obtainedare homogeneous, isotropic, superelastic at room temperature up to strains ofat least 30%, and can become ferroelastic with a shape memory effect up to at least16% strain by an additional heat treatment. The mechanical behavior in compressionis first consolidating and then dilating, while in tension, the samples are slightlyauxetic. A comparison with similar media made of ductile and viscoelastic fibers,as well as with discrete element simulations on friction-free elastic fibers, show thatthe mechanical properties of these self-entangled structures are controlled by theirunique architecture, in-between continuous and discrete media.

Les situations où deux fluides non miscibles sont présents dans un matériau granulaire déformable sont largement rencontrées dans la nature et dans de nombreux domaines de l'ingénierie et de la science. Comprendre l'évolution de tels systèmes multiphases nécessite la connaissance de toutes les phases, leur distribution et interactions. Un modèle micro-hydromécanique couplé est présenté dans cette thèse sur la base de travaux précédents, visant à simuler le drainage quasi-statique de matériaux granulaires déformables. Il combine une approche de type réseau de pores et la méthode des éléments discrets (DEM) pour les fluides et les grains respectivement. Un critère local de mouvement d'interfaces fluides est établi, afin d'approximer au mieux le rôle de la géométrie porale sur les phénomènes capillaires et notamment les forces exercées sur les grains solides à l'intérieur de chaque pore. Une attention particulière est dédiée aux événements de piégeage du fluide drainé et à l'invasion préférentielle le long des bords du domaines. Le modèle est valide par la comparaison avec des résultats expérimentaux (courbes de rétention d'eau). Nous appliquons le modèle pour étudier deux questions: (1) les effets de taille finie et à la question du volume élémentaire représentatif (REV); (2) le paramètre de contrainte effective de Bishop et la relation entre contrainte effective macroscopique contrainte de contact micromécanique. Finalement, une extension du modèle au régime pendulaire est présentée et des premiers résultats sont présentés et discutés. / The situation of two immiscible fluids through a deformable granular material is widely encountered in nature and in many areas of engineering and science. To understand the physical evolution of the multiphase system is of great importance for the applications. It requires the knowledge of all component phases, their distribution and interactions. A pore-scale coupled hydromechanical model is presented in this thesis based on previous work, aiming at simulating the quasi-static drainage of a deformable granular materials. The model combines a pore network approach and the discrete element method (DEM) for the fluids and grains, respectively. A local criterion for determining the local movements of the fluids interfaces established to approximate the role of the local pore geometry on capillarity and namely on the forces exerted on the solid grains inside each pore. Special attentions have been paid to the entrapment events of the receding fluid and to the preferential invasion along the boundaries. The model is validated through comparisons with experimental results (water retention curves). We apply the model for examining two issues: (1) finite size effects and the concept of representative elementary volume (REV); (2) Bishop's effective stress parameter and to the relationship between macro-scale effective stress and micro-scale contact stress. Finally, an extension to the pendular regimes is proposed and first results are presented and analyzed.

Les matériaux granulaires multiphasiques occupent une place très importante dans notre environnement qui suscitent un grand intérêt de nombreuses communautés scientifiques, notamment celles de la mécanique des sols ou de la géotechnique. Le caractère divisé permet les milieux granulaires multiphasiques d'avoir un comportement mécanique global qui trouve leur origine, leur distribution et interactions entre les phases de composition. Un modèle de couplage hydromécanique est présenté dans ce travail de thèse pour l'application à la modélisation microscopique des couplages hydromécaniques dans les matériaux granulaires saturés. Le modèle numérique est basé sur un couplage de la méthode des éléments discrets (DEM) avec une formulation en volumes finis, à l'échelle des pores (PFV), du problmème de l'écoulement d'un fluide visqueux incompressible. Le solide est modélisé comme un arrangement de particules sphériques avec des interactions de type élasto-plastique aux contacts solide-solide. On considère un écoulement de Stokes incompressible, en supposant que les forces inertielles sont négligeables par rapport aux forces visqueuses. La géométrie des pores et leur connectivité sont définies sur la base d'une triangulation régulière des sphères, qui aboutit à un maillage tétrahédrique. La définition des conductivités hydrauliques à l'échelle des pores est un point clef du modèle, qui se rapproche sur ce point à des modèles de type pore-network. Une importance particulière réside dans les lois d'interactions fluide-solide permettant de déterminer des forces de fluide appliquées sur chacune des particules, tout en assurant un coût de calcul acceptable pour la modélisation en trois dimensions avec plusieurs millieurs des particules. Des mesures de perméabilités sur des assemblages bidisperses de billes de verre sont présentées et comparées aux prédictions du modèle et aux formules empiriques/semi-empiriques dans la littérature, ce qui valide la définition de la conductivité locale et met en évidence le rôle de la distribution granulométrique et la porosité. Une approche numérique pour analyser l'interaction mécanique fluide-solide et les mécanismes d'érosion interne est finalement présentée. / Multiphase granular materials occupy a very important place in our environment that are of great interest to many scientific communities, including those of soil mechanics or geotechnical engineering. The divided nature allows multiphase granular media to have a global mechanical behaviour which originates from all component phases, their distribution and interactions. Acoupled hydromechanical model is presented in this work for the application to microscopic modeling of coupled hydromechanical in saturated granular materials. The numerical model uses a combination of the discrete element method (DEM) with a pore-scale finite volume (PFV) formulation of flow problem of an incompressible viscous fluid. The solid is modeled as an assembly of spherical particles, where contact interactions are governed by elasto-plasticrelations. Stokes flow is considered, assuming that inertial forces are small in comparison with viscous forces. Pore geometry and pore connections are defined locally through regular triangulation of spheres, from which a tetrahedral mesh arises. The definition of pore-scale hydraulic conductivities is a key aspect of this model. In this sense, the model is similar to a pore-network model. The emphasis of this model is, on one hand the microscopic description of the interaction between phases, with the determination of the forces applied on solid particles by the fluid, on the other hand, the model involves affordable computational costs, that allow the simulation of thousands of particles in three dimensional models. Permeability measurements on bidispersed glass beads are reported and compared with model predictions and empirical formulas/semi-empirical in the literature, validating the definition of local conductivities and bringing out the role of particle size distribution and porosity. A numerical approach to analyze the fluid-solid mechanical interaction and mechanisms of internal erosion is finally presented.

Le cadre général de ce travail concerne la modélisation du comportement mécanique sous impact de merlons de protection contre les chutes de blocs rocheux. Les ouvrages ciblés sont à technologie cellulaire, et sont composés de gabions à enveloppe grillagée remplis de matériaux granulaires (pierres concassées ou mélange de sable et de déchiquetas de pneus). L'absence de recommandations relatives au dimensionnement de ce type de structure de protection a conduit à la mise en place du projet ANR REMPARe couplant expérimentations in situ et modélisations numériques. Dans le cadre de ce projet, une méthode numérique innovante couplant les approches discrète et continue a été développée. Localement, la méthode des éléments discrets est employée pour conserver la précision de modélisation dans les zones fortement sollicitées proches de l'impact. Le remblai technique aval, peu sollicité, a été discrétisé par une méthode continue aux différences finies. La jonction entre ces deux zones fait l'objet d'un couplage basé sur des critères de conservation d'énergie et de condition cinématique en déplacement. Un modèle de comportement pour les gabions a été proposé et calibré à partir d'essais expérimentaux réalisés à l'échelle de la cellule. Des simulations d'impacts sur un merlon de protection prototype ont été menées puis validées par comparaison avec les mesures issues des expérimentations en vraie grandeur. Afin d'améliorer la compréhension des mécanismes engendrés, l'étude paramétrique portant sur les propriétés des matériaux et les caractéristiques d'impact a permis d'évaluer leur influence sur le mécanisme d'impact. / This study focuses on numerical modelling of the mechanical behaviour of rockfall protection structure submitted to impact. The concerned structures are constituted by cell assembly composed by wire mesh cages filled with granular material (crushed stone, mixture of sand and shredded tires). The lack of design norms led to the set up of the ANR REMPARE project which is based on real scale experiment and numerical modelling. In this project, an innovative numerical method, combining continuous and discrete approaches, has been developed. Locally, the discrete element method is used to maintain the accuracy of modelling in the highly stressed areas near the impacted zone. Downstream, the embankment, lowly stressed, is simulated by a finite difference method. The link between the two domains is ensured by a combined method which is based on energy conservation and respect of border kinematic conditions. A model for the geocells is proposed and calibrated from experimental tests carried out at the cell scale. Impacts on rockfall protection prototype are simulated and validated by comparison with measurements obtained on real scale experiments. Parametric studies on material's properties and characteristics of impact are carried out to evaluate their influence on the impact mechanisms.