Ce travail concerne l'analyse de la vulnérabilité des structures de protection et des ouvrages sensibles en béton soumis à des actions dynamiques sévères (impacts, explosions) dues à des risques anthropiques d'origine accidentelle ou non. L'objet est la mise au point d'outils prévisionnels de simulation capables de décrire de manière objective le comportement dynamique du béton. Pour cela, une approche numérique novatrice reposant sur la méthode des Eléments Discrets (MED) est développée. Une première partie de cette thèse concerne la simulation des essais quasi-statiques de compression et traction uniaxiales. Une loi de transfert de moment (LTM) a été introduite pour pallier au problème de fragilité en compression simple. Ensuite, la procédure d'identification des paramètres du modèle modifié a été optimisée pour bien reproduire le comportement macroscopique du béton. Enfin, le modèle a été validé en représentant correctement le comportement quasi-statique de plusieurs types de béton. La deuxième partie du travail traite la simulation des essais de traction dynamique du béton aux barres de Hopkinson. Les résultats ont montré la nécessité de prendre l'effet de vitesse de déformation dû au matériau pour bien reproduire le comportement expérimental. Ensuite, Les paramètres du modèle permettant de reproduire cet effet de vitesse ont été identifiés. Enfin, des essais avec des taux de déformation très élevés ont été simulés et les résultats numériques ont été en accord avec le comportement observé expérimentalement. / This work concerns the analysis of the vulnerability of sensitive concrete structures subjected to severe dynamic actions such as impacts due to natural hazards or human factors. The object is to develop a numerical tool that can describe objectively the dynamic behaviour of concrete. Then, a 3D discrete element method (DEM) was developed and used to perform the analysis. The first part of this thesis focuses on the simulation of quasi-static uniaxial compression and traction tests. A moment transfer law (MTL) was introduced to overcome the problem of brittle compressive behavior. Then, the identification procedure of the modified DEM model has been optimized in order to reproduce very well the macroscopic behaviour of concrete. Finally, the model has been validated by representing properly the real quasi-static behavior of different types of concrete. The second part of the study deals with the simulation of the dynamic Hopkinson traction bar tests of concrete. The results showed that a local rate effect has to be introduced to reproduce the strain rate dependency, which would then be a material-intrinsic effect. Then, the parameters of the model have been identified. Finally, simulations were run at high strain rates and showed consistent results with respect to experimental behaviour.

Les condensateurs multicouches en céramique (Multilayer Ceramic Capacitors, MLCCs) sont des composants passifs clés de l'électronique moderne. Les MLCCs sont constitués d'une alternance d'électrodes métalliques et de couches diélectriques de céramique. Les puces ultraminces sont composées de mélanges de couches micrométriques métalliques et céramiques et d'additifs de céramique de taille nano. Un certain nombre de défauts tels que des fissures, des délaminations des couches et des discontinuités au sein de l'électrode, peuvent survenir dans la fabrication de ces MLCCs ultraminces. Un dispositif expérimental à rayons X (TXM, Transmission X-ray Microscope) avec une résolution spatiale de 30 nm au synchrotron APS (Advanced Photon Source, USA) a été utilisé pour caractériser un volume cylindrique représentatif de Ø 20 µm × 20 µm extrait d'une puce 0603 (1,6 mm × 0,8 mm) au nickel (Ni) + titanate de baryum (BaTiO3, ou BT) avant et après frittage sous argon hydrogéné (2%). La tomographie 3D de la microstructure montre que les discontinuités de l'électrode finale sont liées à des hétérogénéités initiales dans les couches d'électrodes. La radiographie in-situ aux rayons X pendant le frittage (vitesse de chauffage de 10 °C/mn, température de maintien à 1200 °C pendant 1 heure, puis refroidissement à 15 oC/min) d'un volume représentatif d'électrode au palladium (+ baryum-néodyme-titanate) confirme bien que les discontinuités dans l'électrode proviennent de l'hétérogénéité initiale de la poudre, qui est lié à la nature du compactage d'un matériau particulaire. La discontinuité se produit à l'étape précoce du cycle de frittage. A ce stade, l'électrode métallique commence à fritter tandis que le matériau diélectrique peut être considéré comme un substrat inerte qui contraint le frittage de l'électrode.Des études corrélatives utilisant un FIB-SEM (Focused Ion Beam Scanning Electron-microscopie) en tomographie à haute résolution (5 × 5 × 5 nm3) ont été effectuées sur des échantillons MLCC à vert et frittés. Elles confirment que la résolution de la nanotomographie X est suffisante pour étudier l'évolution des hétérogénéités. Cependant la tomographie par FIB permet à la nanotomographie X d'être réinterprétée avec plus de précision. D'autre part, le FIB fournit les paramètres des particules pour les simulations DEM.La méthode des éléments discrets (DEM) a été utilisée pour simuler la microstructure du système multicouche lors du frittage. Tout d'abord, le frittage de la matrice de nickel avec inclusions BT a été simulé en utilisant le code dp3D. Nous avons pu montrer que la vitesse de densification de la matrice diminue avec l'augmentation la fraction volumique d'inclusions et avec la diminution de la taille des inclusions. Pour une fraction volumique donnée, et une taille d'inclusions donnée, une meilleure dispersion des inclusions conduit à un retard plus marqué de la densification du frittage de la matrice de nickel.Le co-frittage de multicouches de BT/Ni/BT a été simulé en tenant compte des informations collectées à partir de la tomographie FIB-SEM à résolution élevée (taille des particules, distribution de taille, hétérogénéités, et pores). On constate que les discontinuités d'électrodes proviennent des hétérogénéités initiales dans le comprimé à vert et se forment au début de frittage sous contrainte. Ces résultats de simulation sont en bonne correspondance avec les observations expérimentales. Une étude paramétrique indique que les discontinuités d'électrodes peuvent être minimisées par l'homogénéisation de la compacité, par l'augmentation de l'épaisseur des électrodes et par l'utilisation d'un chauffage rapide.A partir des résultats expérimentaux et des simulations DEM, une conclusion générale peut être avancée: la discontinuité finale provient de l'hétérogénéité initiale dans les couches d'électrodes et survient à un stade précoce de frittage lorsque les couches diélectriques contraignent les couches d'électrodes

L'objectif de cette thèse est d'étudier le comportement à la rupture des matériaux fragiles avec une approche capable de relier l'amorçage de fissures à leurs propagations. On adopte la méthode des éléments discrets (DEM), avec le matériau composé d'un ensemble bidimensionnel régulier de particules en contact. Ceci qui nous permet de dériver une expression pour le facteur d'intensité de contraintes en fonction des forces et des déplacements relatifs des deux contacts adjacents à la pointe d'une fissure. Un critère de rupture classique, basé sur la ténacité du matériau, est ensuite adopté pour l'analyse de la propagation des fissures, représenté par la perte des forces de contact entre les particules. La vérification de la formulation est faite par la comparaison des simulations numériques à des solutions classiques de la mécanique de la rupture en mode I, mode II et mode mixte. Ensuite, on étend l'application du critère discret au comportement d'un matériau sain (non pré-fissuré), soumis à des contraintes homogènes, pour lequel on retrouve un critère de rupture du type Rankine. Le résultat final est un modèle discret simple totalement compatible avec les approches continues usuelles, et qui ne demande pas d'essais de calibration, typiques des approches discrètes conventionnelles.

Ce travail concerne l'analyse de la vulnérabilité des structures de protection et des ouvrages sensibles en béton soumis à des actions dynamiques sévères (impacts, explosions) dues à des risques anthropiques d'origine accidentelle ou non. L'objet est la mise au point d'outils prévisionnels de simulation capables de décrire de manière objective le comportement dynamique du béton. Pour cela, une approche numérique novatrice reposant sur la méthode des Eléments Discrets (MED) est développée. Une première partie de cette thèse concerne la simulation des essais quasi-statiques de compression et traction uniaxiales. Une loi de transfert de moment (LTM) a été introduite pour pallier au problème de fragilité en compression simple. Ensuite, la procédure d'identification des paramètres du modèle modifié a été optimisée pour bien reproduire le comportement macroscopique du béton. Enfin, le modèle a été validé en représentant correctement le comportement quasi-statique de plusieurs types de béton. La deuxième partie du travail traite la simulation des essais de traction dynamique du béton aux barres de Hopkinson. Les résultats ont montré la nécessité de prendre l'effet de vitesse de déformation dû au matériau pour bien reproduire le comportement expérimental. Ensuite, Les paramètres du modèle permettant de reproduire cet effet de vitesse ont été identifiés. Enfin, des essais avec des taux de déformation très élevés ont été simulés et les résultats numériques ont été en accord avec le comportement observé expérimentalement. / This work concerns the analysis of the vulnerability of sensitive concrete structures subjected to severe dynamic actions such as impacts due to natural hazards or human factors. The object is to develop a numerical tool that can describe objectively the dynamic behaviour of concrete. Then, a 3D discrete element method (DEM) was developed and used to perform the analysis. The first part of this thesis focuses on the simulation of quasi-static uniaxial compression and traction tests. A moment transfer law (MTL) was introduced to overcome the problem of brittle compressive behavior. Then, the identification procedure of the modified DEM model has been optimized in order to reproduce very well the macroscopic behaviour of concrete. Finally, the model has been validated by representing properly the real quasi-static behavior of different types of concrete. The second part of the study deals with the simulation of the dynamic Hopkinson traction bar tests of concrete. The results showed that a local rate effect has to be introduced to reproduce the strain rate dependency, which would then be a material-intrinsic effect. Then, the parameters of the model have been identified. Finally, simulations were run at high strain rates and showed consistent results with respect to experimental behaviour.

L'objectif principal de ce travail de thèse est d'améliorer la prise en compte des peuplements de taillis dans les logiciels d'analyse trajectographique.Dans un premier temps, un modèle numérique permettant de créer des peuplements virtuels de taillis à l'échelle du versant a été développé et validé sur la base d'inventaires forestiers réalisés dans des taillis.Deux modèles numériques permettant de simuler des impacts de blocs sur des franc pieds et des cépées ont été également développés en utilisant la Méthode des Éléments Discrets (MED).Ces modèles ont été calibrés par des essais d'impact sur des tiges de hêtre.Ils permettent d'intégrer l'influence du houppier et du système racinaire, de modéliser explicitement le contact entre le bloc et les tiges impactées, et d'intégrer les non-linéarités matérielles (rupture des tiges, délaminage) se développant dans le tronc au cours de l'impact.Ces travaux ont conduit à la construction d'un modèle trajectographique MED permettant de simuler la propagation d'un bloc dans une forêt de taillis à l'échelle du versant. Finalement, le rôle protecteur de différents peuplements de taillis contre l'aléa de chute de bloc a été caractérisé à l'aide de ce modèle. / This research work is dedicated to improve the integration of coppice stands in rockfall analyses.First, a model was built to create virtual coppice stands. This model was validated using field inventories in coppice stands.Two numerical models were developed to simulate impacts of blocks on single trees and coppice stools using the Discrete Elements Method (MED).These models were calibrated using laboratory impact tests on beech stems.They account for the influence of the root system and of the crown on the tree dynamic response, the explicit modeling of the contact between the block and the impacted stem and the non-linearity evolution into the trunk during impact.Finally, a DEM rockfall software was developed to model rockfall propagation in coppice stands.The protective role against rockfall hazard of different coppice stands was characterized with this model.

L’objectif des travaux présentés dans ce mémoire de thèse est de développer une modélisation numérique de la compaction des poudres composées de particules sécables dans le cadre de la méthode de Dynamique des Contacts en vue d’application au procédé de fabrication du combustible nucléaire. Les particules sont modélisées comme des agrégats cohésifs de fragments potentiels (cellules) de formes polyédriques irréguliers. A l’aide de ce modèle de cellules liées (Bonded Cell Method), nous avons réalisé une étude paramétrique de la résistance des particules par rapport aux paramètres géométriques et mécaniques du modèle. Nos résultats révèlent deux régimes et une mise à l’échelle en loi de puissance de la résistance à la compression en fonction de l’adhésion normale et du rapport entre l’adhérence tangentielle et l’adhésion normale entre cellules. Nous avons optimisé les paramètres du modèle pour la compaction uni-axiale des assemblages d’un grand nombre de particules sécables. Les simulations ont permis d’identifier les mécanismes de compaction et de rupture des particules, et de caractériser l’évolution de la texture et des tailles et formes des fragments. Les résultats obtenus montrent clairement que le processus de compaction est fortement non-linéaire en raison notamment de l’évolution de l’étalement granulométrique qui contrôle la texture et la transmission des contraintes. Enfin, nous avons mené une étude systématique de l’effet de la polydispersité de taille dans le cas de particules sphériques. / The goal of this PhD work is to develop a numerical modeling approach of the compaction of powders composed of crushable particles in the framework of the Contact Dynamics method in view of application to the manufacture process of nuclear fuel. The particles are modeled as cohesive aggregates of potential fragments (cells) of irregular polyhedral shape. Using this Bonded Cell Method, we performed a parametric investigation of the strength of particles with respect to the geometrical and mechanical model parameters. Our results reveal two regimes and a power-law scaling of the compressive strength as a function of the ratio between tangential adherence and normal adhesion between cells. We optimized the model parameters for the uniaxial compression of packings of a large number of crushable particles. The simulations allow us to identify the mechanisms of compaction and fracture of particles, and to characterize the evolution of texture and the sizes and shapes of fragments. Our simulation results clearly show that the compaction process is strongly nonlinear as a consequence of the evolution of particle size distribution, which controls the texture and stress transmission. We also perform a systematic analysis of the effect of size polydispersity in the case of spherical particles.

Le formage de matériaux multiphasés comprend des mécanismes complexes en lien avec la rhéologie,la morphologie et la topologie des phases.Du point de vue numérique,la modélisation de ces phénomènes en résolvant les équations aux dérivées partielles (EDP) décrivant le comportement continu des phases n'est pas trivial.En effet,de nombreuses discontinuités associées aux phases se déplacent et peuvent interagir.Ces phénomènes peuvent être conceptuellement déclicats à intégrer au modèlecontinu et coûteux en termes de calcul.Dans cette thèse,la méthode des éléments discrets (DEM) est utilisée pour modéliser phénoménologiquement les déformations finies inélastiques dans les multi-matériaux.Les lois d'interactions attractive-répulsive sont imposées à des particules fictives,dont les ré-arrangements collectifs modélisent les déformations irréversibles de milieux continus.Le comportement numérique des empilements de particules est choisi pourreproduire des traits caractéristiques de la viscoplasticité parfaite etisochore:contrainte d'écoulement,sensibilité à la vitesse de déformation,conservation du volume.Les résultats d'essais de compression de bi-matériaux simples,simulés avec la DEM,sont comparés à la méthode des éléments finis (FEM) et sont en bon accord.Le modèle est entendu pour pouvoir supporter des sollicitations de traction.Une méthode de détection de contacts et d'auto-contacts d'objets physiques estproposée,basée sur l'approximation locale des surfaces libres.Les capacités de la méthodologie globale sont testées sur des mésostructurescomplexes,obtenues par tomographie aux rayons X.La compression à chaud d'un composite métallique dense est modélisée.La co-déformation peut être observées à l'échelle spatiale des phases.Deux cas de matériaux ``poreux'' sont considérés.Premièrement la simulation de la compression puis traction d'alliagesd'aluminium présentant des pores.Ces pores proviennent du coulage du matériau,leur fermeture et ré-ouverture mécanique est modélisée,y compris la coalescence à grande déformation.Deuxièmement la simulation de la compression de mousse métallique de faibledensité.Typiquement utilisée dans le but d'absorber de l'énergie mécanique,la compression jusqu'à densification provoque de nombreuses interactions entreles bras de matière. / Forming of multiphase materials involves complex mechanisms linked with therheology,morphology and topology of the phases.From a numerical point of view,modeling such phenomena by solving the partial differential equation (PDE) system accounting for thecontinuous behavior of the phases can be challenging.The description of the motion and the interaction of numerous discontinuities,associated with the phases,can be conceptually delicate and computationally costly.In this PhD,the discrete element method (DEM) is used to phenomenologically model finite inelastic strain inmulti-materials.This framework,natively suited for discrete phenomena,allows a flexible handling of morphological and topological changes.Ad hoc attractive-repulsive interaction laws are designed betweenfictitious particles,collectively rearranging to model irreversible strain in continuous media.The numerical behavior of a packing of particles can be tuned to mimic keyfeatures of isochoric perfect viscoplasticity:flow stress, strain rate sensitivity, volume conservation.The results for compression tests of simple bi-material configurations,simulated with the DEM,are compared to the finite element method (FEM) and show good agreement.The model is extended to cope with tensile loads.A method for the detection of contact and self-contact events of physicalobjects is proposed,based on a local approximation of the free surfaces.The potential of the general methodology is tested on complex mesostructuresobtained by X-ray tomography.The high temperature compression of a dense metallic composite is modeled.The co-deformation can be observed at the length scale of the phases.Two cases of ``porous'' material are considered.Firstly,the simulation of the compression and the tension of aluminum alloys with poresis investigated.These pores stem from the casting of the material,their closure and re-opening is modeled,including the potential coalescence occurring at large strain.Secondly,the compression of a metallic foam,with low relative density,is modeled.Typically used in energy absorption applications,the compression up to densification involves numerous interactions between thearms.

Un modèle numérique décrivant le comportement de fibres souples en suspension dans un écoulement de fluide en régime inertiel a été développé au moyen d'un couplage entre la méthode des éléments discrets et la méthode des volumes finis. Chaque fibre est discrétisée en plusieurs éléments de type poutre permettant de prendre en compte une déformation (flexion, torsion, allongement) et un mouvement de corps rigide. Les équations du mouvement des fibres sont résolues au moyen d'un schéma explicite du second ordre (temps et espace). Le mouvement de la phase fluide est décrit par les équations de Navier-Stokes, qui sont discrétisées et résolues au moyen d'un schéma aux volumes finis non structurés, d'ordre 4 (temps et espace). Le couplage entre la phase solide (discrète) et la phase fluide (continue) est obtenue par une pseudo méthode IBM (Immersed Boundary Method) dans laquelle l'effort hydrodynamique est calculé analytiquement. Plusieurs modèles de force hydrodynamique issus de la littérature sont analysés et leur validité ainsi que leurs limites sont identifiées. Pour des nombres de Reynolds (Re) correspondant au régime inertiel (0.01 < Re < 100, Re défini à l'échelle de la fibre), des formulations non-linéaires de la force hydrodynamique exercée par un écoulement uniforme sur un cylindre infini sont utilisées. Le couplage a aussi été utilisé pour des fibres rigides en écoulement de Stokes, en utilisant l'expression de la force de traînée issue de la théorie des corps élancés (`slender body theory'). Une expression du moment hydrodynamique par unité de longueur est obtenu à partir de simulations numériques par volumes finis de l'écoulement autour d'un cylindre élancé.Le modèle développé a été validé par comparaison avec plusieurs résultats expérimentaux et analytiques, du régime de Stokes (pour des fibres rigides) jusqu'aux régimes inertiels. Dans le cas du régime de Stokes, des simulations numériques du cisaillement de suspensions de fibres semi-diluées ont été réalisées. Le modèle développé permet de capturer les interactions hydrodynamiques et non-hydrodynamiques entre les fibres. Les interactions élasto-hydrodynamiques pour $Re$ fini ont été validées dans deux cas. Dans le premier cas, la flèche d'une fibre encastrée-libre dans un écoulement uniforme a été obtenu par calcul numérique et le résultat validé par comparaison aux résultats expérimentaux de la littérature. Dans le second cas, la conformation de fibres élancées et très déformables dans un écoulement turbulent homogène et isotrope a été obtenu par calcul numérique et le résultat validé par comparaison aux résultats expérimentaux de la littérature. Deux études numériques ont été réalisées pour étudier l'effet de la présence de fibres en suspension sur la turbulence au sein du fluide suspensif. Le modèle numérique a permis de reproduire le phénomène de réduction/amplification de la turbulence dans un écoulement en canal ou en conduite, dû à l'évolution microstructurale de la phase fibreuse. / A numerical model describing the behavior of flexible fibers under inertial flows was developed by coupling a discrete element solver with a finite volume solver.Each fiber is discretized into several beam segments, such that the fiber can bend, twist and rotate. The equations of the fiber motion were solved usinga second order accurate explicit scheme (space and time). The three dimensional Navier-Stokes equations describing the motion of the fluid phase was discretizedusing a fourth th order accurate (space and time) unstructured finite volume scheme. The coupling between the discrete fiber phase and the continuous fluid phasewas obtained by a pseudo immersed boundary method as the hydrodynamic force on the fiber segments were calculated based on analytical expressions.Several hydrodynamic force models were analyzed and their validity and short-comings were identified. For Reynolds numbers (Re) at the inertial regime(0.01 < Re < 100, Re defined at the fiber scale), non linear drag force formulations based on the flow past an infinite cylinder was used. For rigid fibers in creeping flow, the drag force formulation from the slender body theory was used. A per unit length hydrodynamic torque model for the fibers was derived from explicit numerical simulations of shear flow past a high aspect ratio cylinder. The developed model was validated against several experimental studies and analytical theories ranging from the creeping flow regime (for rigid fibers) to inertial regimes. In the creeping flow regime, numerical simulations of semi dilute rigid fiber suspensions in shear were performed.The developed model wasable to capture the fiber-fiber hydrodynamic and non-hydrodynamic interactions. The elasto-hydrodynamic interactions at finite Reynolds was validated with against two test cases. In the first test case, the deflection of the free end of a fiber in an uniform flow field was obtained numerically and the results were validated. In the second test case the conformation of long flexible fibers in homogeneous isotropic turbulence was obtained numerically and the results were compared with previous experiments. Two numerical studies were performed to verify the effects of the suspended fibers on carrier phase turbulence and the numerical model was able to reproduce the damping/enhancement phenomena of turbulence in channel and pipe flows as a consequence of the micro-structural evolution of the fibers.

L’objectif de cette thèse est d’étudier le comportement à la rupture des matériaux fragiles avec une approche capable de relier l’amorçage de fissures à leurs propagations. On adopte la méthode des éléments discrets (DEM), avec le matériau composé d'un ensemble bidimensionnel régulier de particules en contact. Ceci qui nous permet de dériver une expression pour le facteur d'intensité de contraintes en fonction des forces et des déplacements relatifs des deux contacts adjacents à la pointe d’une fissure. Un critère de rupture classique, basé sur la ténacité du matériau, est ensuite adopté pour l'analyse de la propagation des fissures, représenté par la perte des forces de contact entre les particules. La vérification de la formulation est faite par la comparaison des simulations numériques à des solutions classiques de la mécanique de la rupture en mode I, mode II et mode mixte. Ensuite, on étend l’application du critère discret au comportement d’un matériau sain (non pré-fissuré), soumis à des contraintes homogènes, pour lequel on retrouve un critère de rupture du type Rankine. Le résultat final est un modèle discret simple totalement compatible avec les approches continues usuelles, et qui ne demande pas d’essais de calibration, typiques des approches discrètes conventionnelles. / The objective of this thesis is to study the fracture behavior of brittle materials by an approach which relates crack initiation to crack growth. We adopt the discrete element method (DEM) and we represent the material by a 2D regular set of particles in contact. This allows us to derive an expression for the stress intensity factor as a function of the forces and relative displacements of two adjacent contacts at the crack tip. A classical failure criterion, based on the material’s toughness, is then adopted for the analysis of crack propagation, represented by the loss of contacts forces between particles. The formulation is verified by the comparison of numerical simulations to classical solutions of fracture mechanics in mode I, mode II and mixed mode. Afterwards, we apply our discrete criterion to uncracked materials under homogenous stress conditions, obtaining a Rankine like behavior.. The work results in a simple discrete model which is totally compatible to continuum mechanics, where no calibration tests are required, in contrast to most of discrete approaches.

Ce travail de thèse a pour objet d’analyser le processus de dissipation d’énergie dans les amortisseurs par milieux granulaires. Des études de nature expérimentale, analytique et numérique ont été menées afin de pouvoir détecter les paramètres clefs agissant sur la dissipation; un modèle expérimental minimal a été présenté en premier lieu afin de mettre en évidence l’efficacité des milieux granulaire en tant qu’amortisseurs de vibrations. Un second modèle expérimental plus élaboré a été développé, avec de multiples protocoles expérimentaux, pour mener une étude paramétrique et détecter leurs impacts sur l’évolution du facteur de perte du système. On montre que le coefficient de perte ne dépend pas du matériau des particules ou leur nombre, mais dépend fortement de la masse totale des grains embarqués et sur l’amplitude du signal vibrant. Nos mesures montrent aussi la contribution de l'écoulement visqueux de l'air entourant les grains au facteur de perte globale des amortisseurs.La partie analytique à son tour a permis de retrouver le comportement obtenu expérimentalement par le billet du développement des équations du mouvement du système, celle des énergie cinétique et énergie dissipée afin de proposer enfin une équation maitresse qui décrit l’évolution du facteur de perte réduit au sein de notre système. Pour atteindre une plus grande précision, une modélisation du système granulaire par la méthode des éléments discrets (DEM) a permis de retrouver les mêmes conclusions et ainsi valider les constatations expérimentales et le modèle analytique proposé. / This thesis aims to analyze the process of energy dissipation in particle dampersExperimental, analytical and numerical studies have been conducted in order to identify key parameters influencing the dissipation; minimal experimental model was introduced first to highlight the efficiency of granular media as shock absorbers of vibrations. A second more sophisticated experimental model was developed, with multiple experimental protocols, to conduct a parametric study and detect their impact on the evolution of the system loss factor. It is shown that the loss coefficient is independent of the particle material or their number, but depends strongly on the total mass of embedded grains and on the amplitude of the vibrating signal. Our measurements also show the contribution of viscous flow of the air surrounding the grains to the overall loss factor.The analytical part in turn led to the discovery The behaviour obtained experimentally by the development of the equations of motion of the system, that of kinetic energy dissipated and energy to finally offer a mistress equation which describes the evolution of the loss factor reduced within our system. To achieve greater accuracy, a model of the granular system by the discrete element method (DEM) allowed to find the same conclusions and thus validate the experimental findings and the proposed analytical model.