Cette thèse porte sur l'étude de l'écoulement et du transfert thermique de polymère dans les deux procédés de rotomoulage: rotomoulage réactif et rotomoulage de poudres thermoplastiques. L'écoulement du polymère réactif au cours du rotomoulage a été modélisé par l'écoulement d'un fluide réactif à surface libre sous gravité à l'intérieur d'un cylindre horizontal en rotation autour de son axe principal. Dans ce cas, la méthode lagrangienne ''Smoothed Particle Hydrodynamics (SPH)'' a été appliquée en utilisant une expression adaptée de la viscosité variable du polymère réactif. À l'état de poudre, l'écoulement et le transfert thermique de polymère au cours du rotomoulage des poudres thermoplastiques a été modélisés en utilisant un code de calcul basé sur la ''Méthode des Éléments Discrets (DEM)''. Les régimes d'écoulement, la ségrégation, le mélange et le transfert thermique des poudres dans un moule en rotation ont été simulés. Chaque fois que cela était possible des résultats des modèles ont été validés expérimentalement avec succès. / This Ph.D research work is about the study of polymer flow and heat transfer in both processes of rotational molding: reactive rotational molding and rotational molding of thermoplastic powders. The flow of the reactive polymer in the reactive rotational molding was modeled by a free surface flow under gravity inside a horizontal rotating cylinder. The meshless and lagrangian method ''Smoothed Particle Hydrodynamics (SPH)'' was applied to simulate the flow by using an adapted variable viscosity expression. In powder state, the flow and heat transfer of polymer during rotational molding of thermoplastic powders are modeled using a solver based on the ''Discrete Element Method (DEM)''. In which, the flow regime, the segregation, the mixture and the heat transfer of polymer powders inside a horizontal rotating cylinder was simulated. Results of the model have been validated experimentally with success.

L'agglomération des fines entraîne une augmentation importante de la demande en eau des mélanges cimentaires et constitue un frein à la valorisation des sous-produits contenant des proportions élevées de fines. Une modélisation discrète de la microstructure des agglomérats de fines particules est développée. Le modèle intègre les forces de contact avec prise en compte du frottement et du roulement, les forces physico-chimiques et les forces hydrodynamiques calculées avec l'approximation de drainage libre. Le processus d'agglomération, de fragmentation et de restructuration sous cisaillement d'agglomérats fractals de fines particules minérales est étudié en vue de quantifier l'eau immobilisée dans les flocs. Un agglomérat initial isolé et lâche est soumis à un taux de cisaillement constant. Lorsque les forces hydrodynamiques sont prépondérantes par rapport aux forces de cohésion, l'agglomérat mère se fragmente en agglomérats secondaires jusqu'à atteindre un état stationnaire, où le nombre, la taille et la structure des flocs peuvent être considérés comme constants. L'influence du gradient de vitesse, du contexte physico-chimique, de la taille et de la polydispersité des particules sur les caractéristiques microstructurales des flocs à l'équilibre est étudiée. Il est montré que le comportement d'une suspension diluée monodisperse peut être décrit par un unique nombre adimensionnel faisant intervenir la force de cohésion maximale, la taille des particules et le gradient de vitesse. / The agglomeration of fines leads to a significant increase in water demand of cement mixtures and hinders the valorization of by-products containing high fines proportions. A discrete modeling of the microstructure of fine particles agglomerates is developed. The model takes into account contact forces, physico-chemical, and hydrodynamic forces computed with the free draining approximation. The agglomeration, fragmentation and restructuring process of fractal agglomerates of fine mineral particles under shear is studied to quantify the entrapped water in the flocs. A loose initial agglomerate isolated is submitted to a constant shear rate. When the hydrodynamic forces predominate over the adhesive forces, the parent agglomerate is broken up into secondary agglomerates until reaching a steady state, where the number, size and structure of flocs can be considered constant. The influence of the velocity gradient, the physico-chemical context, particles size and their polydispersity on the microstructural characteristics of flocs at steady state is studied. It is shown that the behavior of a monodisperse dilute suspensions can be described by a dimensionless single number involving maximum cohesion force, particle size and velocity gradient.

L'objectif de ce travail de recherche est le développement d'un outil numérique capable de simuler le comportement d'infrastructures sensibles soumises à des charges dynamiques extrêmes sous l'effet d'aléas naturels ou humains tels que les impacts aériens. L'étude propose pour ce faire une modélisation 3D par éléments discrets, capable de décrire des états de destructions avancés en obtenant des macro-fissures et des fragments réalistes grâce à la nature discontinue du modèle.Dans un premier temps, on a étudié de manière exhaustive l'influence des paramètres de création du maillage sur les caractéristiques dudit maillage et sur le comportement macroscopique du béton. Ensuite, on a transformé le modèle de charge dynamique en une modélisation plus réaliste de l'énergie de rupture dynamique en contrôlant l'augmentation de distance limite maximum à l'interaction. Par ailleurs, on a défini une condition de ratio entre la taille des éléments discrets et celle des éléments finis pour s'assurer d'un contact correct entre eux.La procédure d'identification des paramètres du modèle est réalisée en simulation numérique avec des essais en laboratoire: Compression et traction quasi-statiques, essai tri-axial à haut confinement, écaillage dynamique. Enfin, la fiabilité de l'approche est vérifiée sur trois essais différents d'impact violent:1) Perforation et pénétration de projectiles à tête ogivale dans des cibles cylindriques confinées (CEA-Gramat)2) Essais d'impact tranchant de projectile avec une géométrie homothétique particulière sur des dalles en béton (Erzar) / The objective of this work is development of a numerical tool capable to simulate sensitive infrastructures subjected to severe dynamic loadings due to natural or manmade hazards, such as aircrafts impacts. This study proposes a 3D discrete element method able to predict advance damage states obtaining realistic macro-cracks and materials fragments due to its discontinue nature.We thoroughly studied the influence of mesh creation parameters on the mesh characteristics and on the macroscopic concrete behaviour. Then, we modified the constitutive model for dynamic loading with a more realistic modelling of the dynamic fracture energy by controlling the increase of the maximum distance limit at the interaction scale. In addition, we defined a condition for ratio between the size of finite and discrete elements for proper contact between them.The identification procedure for the parameters of the constitutive model is analysed with simulation on laboratory test: quasi-static compression and tension, high confinement triaxial and dynamic spalling. Finally, the reliability of our approach is verified on three different types of hard impact test: 1) perforation and penetration of ogive-nosed projectiles on confined cylindrical targets (CEA-Gramat); 2) edge-on impact tests of projectiles with a particular homothetic geometry on concrete tiles (Erzar); 3) drop-weight impact on a reinforced concrete beam (University of Toronto).

Les ouvrages type merlons, construits pour se prémunir contre les chutes de blocs rocheux, sont dimensionnés en pratique au moyen de méthodes empiriques. Afin d'optimiser la géométrie de ces constructions leur modélisation a été entreprise. La méthode aux éléments discrets est utilisée en raison des fortes restructurations du remblai constitutif pouvant intervenir durant un impact. A partir d'un état de l'art, les caractéristiques mécaniques du remblai sont évaluées. La calibration des paramètres quasi-statiques du matériau modèle est alors réalisée par simulations d'essais triaxiaux. Les paramètres gérant le comportement dynamique sont calibrés à leur tour par simulations d'impacts à énergies modérées. Le remblai modèle obtenu est validé par la modélisation d'essais d'impacts expérimentaux dont la gamme d'énergies correspond à celle rencontrée pour les merlons. Il est alors utilisé pour simuler des impacts sur merlons. Plusieurs tailles d'ouvrages sont étudiées, répondant chacune à des gammes d'énergies différentes. Une analyse paramétrique portant sur divers aspects du bloc et de l'ouvrage est réalisée. Elle permet de définir les variables principales gérant le phénomène d'impact. La variation couplée de ces dernières est réalisée afin d'établir leur influence sur la capacité du merlon à contenir le bloc.

Dans la fabrication de lentilles asphériques pour la photolithographie, l'étape delissage est critique. C'est aujourd'hui le seul procédé qui peut corriger les défauts de hautesfréquences spatiales responsables de diffusions parasites, de diminutions de transmittance etde contraste. Cette opération doit préserver la forme asphérique basse fréquence tout enlissant les défauts de hautes fréquences. Un tel comportement peut être obtenu pour des outilscombinant une couche continue flexible pour le suivi basse fréquence et une couche de poixfractionnée pour le lissage de défauts hautes fréquences. Les buts principaux de cette étudesont de prédire l'efficacité de lissage et le suivi de forme de différents outils et ensuite dedéterminer l'outil permettant le meilleur compromis. A cette fin, un modèle multi-échelles estdéveloppé. A l'échelle de l'outil entier, une étude paramétrique par éléments finis permet dedéterminer les caractéristiques de la couche flexible ainsi que la force appliquée optimale afind'obtenir l'homogénéité de la pression à l'échelle de la forme asphérique globale. A l'échelle dela couche de poix, la Méthode par Eléments Discrets est utilisée pour investiguer l'interfaceoutil-pièce. Un modèle basé sur le concept de la poutre cohésive viscoélastique est développé,prenant en compte la réponse fréquentielle de la couche de poix. La comparaison avecl'Analyse Mécanique Dynamique montre la capacité de la DEM à modéliser le comportementviscoélastique. L'opération de lissage est ensuite modélisée par DEM et analytiquement. Lesdonnées expérimentales obtenues par la méthode de la Densité Spectrale de Puissancemettent en évidence l'impact des propriétés de la poix sur l'efficacité de lissage. Lacomparaison entre les résultats numériques et expérimentaux montre que le modèle DEMdonne des résultats prometteurs pour la modélisation du lissage de défauts.

La technique de changement d’échelles a été largement développée dans la littérature pour décrire le comportement global des milieux granulaires en prenant en compte leurs propriétés locales. Cette technique considère classiquement deux échelles : l’échelle macroscopique du volume élémentaire représentatif et l’échelle microscopique du contact entre particules. Le défi majeur de ce changement d’échelles "micro-macro" réside dans la définition de la déformation macroscopique : en effet, si la contrainte macroscopique peut être clairement définie à partir des forces de contact, il a été montré qu’il n’était pas approprié de déduire la déformation macroscopique à partir de la cinématique aux contacts. Dans ce cadre, ce travail propose d’introduire une troisième échelle dite mésoscopique. Cette échelle, à laquelle peuvent être définies à la fois contrainte et déformation, est intermédiaire entre les échelles microscopique et macroscopique et permet de palier au défi majeur mentionné ci-dessus. Elle est définie au niveau d’arrangements locaux de particules, appelés sous-domaines, et sa pertinence est étudiée sur la base d’échantillons numériques composés de particules circulaires puis sphériques, simulés par la méthode des éléments discrets. Les milieux bidimensionnels sont géométriquement représentés par un graphe de particules composé de sous-domaines fermés, encore appelés cellules de vide, dont la frontière est constituée de branches connectant les centres de particules en contact : l’échelle mésoscopique est donc définie au niveau de ces cellules de vide fermées. A cette échelle locale, on décrit tout d’abord la structure du milieu en termes de densité et de texture puis l’on définit les variables statique et cinématique locales du milieu en termes de contrainte et déformation. De fortes hétérogénéités des milieux granulaires en termes de structure, déformation et contrainte sont mises en évidence à l’échelle mésoscopique, avec de plus une structuration des hétérogénéités de contraintes et de déformations et une forte corrélation entre ces deux quantités. Concernant les milieux tridimensionnels, une partition en cellules de vide fermées est impossible du fait de la complexité de la structure 3D de ces milieux. On propose donc une méthode de partition du milieu basée sur la distribution des vides en son sein. La méthode consiste en premier lieu en une subdivision du milieu en tétraèdres, par une partition de Delaunay, puis en une association de tétraèdres voisins selon un critère prédéfini en vue de la création de sous-domaines, non fermés, mais au rôle analogue aux sous-domaines fermés de l’étude 2D. Le critère d’association proposé est basé sur le rapport entre la taille des constrictions (vide sur chaque face des tétraèdres) et la taille des pores au voisinage de chaque constriction. Cette méthode d’association constitue donc l’étape préliminaire à l’extension au cas tridimensionnel des résultats obtenus dans le cas bidimensionnel. / The technique of change of scales has been extensively developed in the literature to describe the global behaviour of granular materials taking into account their local properties. This method usually considers two scales : the macroscopic scale at the level of representative elementary volume and the microscopic scale at the level of contact between particles. The major difficulty of this "micro-macro" change of scales lies in the definition of the macroscopic strain : indeed, the macroscopic stress is clearly defined from contact forces while it is not appropriate to derive the macroscopic strain from the kinematics at contacts. In this framework, this work proposes to introduce a third scale called mesoscopic scale. This scale, at which both stress and strain can be defined, is intermediate between macroscopic and microscopic scales and allows to overcome the major difficulty mentioned above. The mesoscopic scale is defined at the level of local arrangements of particles, called sub-domains, and its relevance is studied on numerical 2D and 3D materials composed of circular then spherical particles, simulated with the discrete element method. Bidimensional media are geometrically represented by a particle graph composed of closed sub-domains, also called void cells, whose border is constituted by the branches joining the centers of particles in contact : the mesoscopic scale is thus defined at the level of these closed void cells. At this local scale, we fist describe the structure of the medium in terms of density and fabric ; we define then the static and kinematic variables in terms of stress and of strain. Strong heterogeneities of granular media in terms of structure, stress and strain are highlighted at this scale, with a structuration of heterogeneities of stress and strain and a significant correlation between these two quantities. Concerning tridimensional media, a partition into closed void cells is impossible, because of the complexity of the 3D structure of these media. We propose a partition method based on the distribution of voids inside the medium. This method consists first in subdividing the medium into tetrahedrons by a Delaunay partition and then in associating neighbouring tetrahedrons, according to a criterion to be defined. This allows us to form sub-domains which are not closed but which play a role analogous to the role of closed sub-domains in the 2D study. The proposed association criterion is based on the ratio between the size of constriction (void on each face of tetrahedrons) and the size of pores around each constriction. This partition method constitutes a preliminary step for an extension of the results obtained in the bidimensional case to the tridimensional case.

Ce travail de recherche a pour objectif d’étudier le comportement au cisaillement des matériaux granulaires homogènes et hétérogènes par l’essai triaxial. Les travaux sont réalisés sur des billes de verre et comportent à la fois des essais triaxiaux de laboratoire et des simulations numériques en DEM (Discret Element Method). D’un point de vue expérimental, les essais sont tout d’abord réalisés sur des échantillons homogènes (billes monodisperses). Ensuite, on étudie des systèmes de billes monodisperses comportant une inclusion de grande taille et des mélanges de billes bidisperses pour représenter les cas hétérogènes. D’un point de vue numérique, une condition limite cylindrique paroi rigide, basée sur la formule de Lamé, est intégrée dans un code aux éléments discrets existant. Une série de procédures est proposée afin de modéliser l’essai triaxial dans des conditions comparables aux essais. Les résultats numériques et expérimentaux sont ensuite comparés, à la fois pour les systèmes monodisperses et pour les milieux hétérogènes. Le modèle numérique permet de reproduire de manière très satisfaisante la courbe déviatorique dans toutes les conditions expérimentales utilisées. Il permet également de décrire correctement les déformations volumiques de l’échantillon, malgré la simplicité de la condition aux limites utilisée. L’outil numérique développé constitue ainsi un complément intéressant aux essais expérimentaux. Il permet en effet d’étudier l’influence de paramètres micromécaniques sur le comportement global des matériaux. / This research work aims at studying the shear behavior of homogeneous and heterogeneous granular materials by triaxial test. The work is performed on glass beads both in laboratory tests and by numerical simulations in DEM (Discrete Element Method). From an experimental point of view, tests are first performed on homogeneous samples (monodisperse beads), then expands to heterogeneous cases like monodisperse systems comprising a large inclusion and bidisperse mixtures. From a numerical point of view, a cylindrical rigid wall boundary condition, based on the Lamé formula is integrated into an existing discrete elements code. A series of procedures is proposed to model the triaxial test in conditions similar to experimental tests. The numerical and experimental results are compared both for monodisperse and heterogeneous systems. The numerical model can reproduce deviatoric curves very satisfactorily in all experimental conditions. It also allows describing correctly the volumetric strains of the sample despite the simplicity of the boundary condition used. The developed numerical tool is an interesting complement to experiment triaxial tests. It makes it possible to study the influence of micromechanical parameters on the overall behavior of materials.

Les anodes en carbone font partie de la réaction chimique de réduction de l'alumine, qui est consommée lors du procédé d'électrolyse Hall-Héroult. Le comportement des agrégats secs de coke en tant que composant principal des anodes de carbone (environ 85 %) a un rôle clé et exceptionnel dans leurs propriétés finales. L'analyse de défaillance des agrégats de coke sec permet non seulement de mieux comprendre les mécanismes de déformation des matériaux granulaires sous charge compressive, mais peut également identifier les causes potentielles de défauts structurels des anodes en carbone, telles que les fissures horizontales. Dans ce travail, il sera montré qu'un mode de défaillance particulier peut être responsable de la génération de fissures dans les anodes en carbone. Le comportement de rupture des agrégats de coke n 'est pas seulement affecté par les paramètres du processus de compactage, tels que la pression de confinement et la vitesse de déformation axiale, mais il dépend également fortement de la distribution granulométrique et de la forme des particules de coke. La méthode des éléments discrets (DEM) est utilisée pour modéliser le comportement micromécanique des agrégats de coke sec pendant le processus de compactage. De plus, le critère de travail de second ordre est utilisé pour analyser la rupture des éprouvettes de granulats de coke. Les résultats révèlent que l'augmentation de la pression de confinement augmente la probabilité du mode de diffusion de la rupture dans l'éprouvette. D'autre part, l'augmentation de la vitesse de déformation augmente les chances du mode de localisation de la déformation de la rupture dans l'éprouvette. De plus, les résultats indiquent que l'utilisation de fines particules ainsi que la diminution de la sphéricité des particules de coke augmenteront la plage de stabilité des agrégats de coke. De plus, en utilisant l'analyse des évaluations de contour de micro-déformation pendant le processus de compactage, il est montré que, à la fois en ajoutant des particules fines aux agrégats de coke et en diminuant la sphéricité des particules de coke, la possibilité de créer une bande de compression dans le coke agrégats est réduit. Étant donné que la présence des bandes de compactage dans la pâte d'anode crée une zone sujette à la génération de fissures horizontales, les résultats de cette étude pourraient conduire à la production d'anodes en carbone avec moins de défauts structurels. / Carbon anodes are part of the chemical reaction of alumina reduction, that is consumed during the Hall-Héroult electrolysis process. The behavior of dry coke aggregates as the main component of carbon anodes (about 85 %) has an exceptional key role in their final properties. The failure analysis of dry coke aggregates not only leads to a better understanding of the deformation mechanisms of granular materials under compressive loading but also can also identify the potential causes of structural defects in carbon anodes, such as horizontal cracks. In this work, it will be shown that a particular failure mode can be responsible for the crack generation in the carbon anodes. The failure behavior of the coke aggregates is not only affected by the compaction process parameters, such as the confining pressure and axial strain rate, but it is also strongly dependent on the size distribution and shape of coke particles. The discrete element method (DEM) is employed to model the micro-mechanical behavior of the dry coke aggregates during the compaction process. In addition, the second-order work criterion is used to analyze the failure of the coke aggregate specimens. The results reveal that increasing the confining pressure enhances the probability of the diffusing mode of the failure in the specimen. On the other hand, the increase of the strain rate augments the chance of the strain localization mode of the failure in the specimen. In addition, the results indicate the fact that the use of fine particles as well as decreasing the sphericity of coke particles will increase the stability range of the coke aggregates. Moreover, by using the analysis of micro-strain contour evaluations during the compaction process, it is shown that, both by adding fine particles to the coke aggregates and by decreasing the sphericity of coke particles, the possibility of creating a compression band in the coke aggregates is reduced. Since the presence of the compaction bands in the anode paste creates an area that is prone to horizontal crack generation, the results of this study could lead to the production of carbon anode with fewer structural defects.

La dysplasie acétabulaire est une pathologie causant de l'instabilité articulaire et une réduction de la superficie de chargement au niveau de la hanche. Celle-ci augmente le risque de développement précoce d'ostéoarthrose, causant de la douleur et menant à la perte de la fonction articulaire. La dysplasie est conventionnellement traitée par ostéotomie péri-acétabulaire, mais les effets précis de cette correction sur la biomécanique intra-articulaire restent à être élucidées. À l'aide de la méthode de simulation par éléments discrets, la distribution de pression articulaire a été évaluée de manière pré et post opératoire chez 9 patients dysplasiques. Différentes mesures orthopédiques évaluant le degré de correction acétabulaire ont été acquises et plusieurs scores numériques évaluant la surcharge articulaire ont été calculés. Des corrélations entre ces mesures et ces scores ont ensuite été effectués. Il a été noté que les mesures d'angle centre-externe latéral, d'index acétabulaire et d'index d'extrusion de la tête fémorale sont corrélés avec l'augmentation de l'aire de contact moyen durant la marche (p < 0.05, R² de 0.22 à 0.50). L'influence significative de la méthode de reconstruction cartilagineuse (distributions d'épaisseurs moyennes saine vs. dysplasique) et du profil de marche utilisé (phase d'appui moyen sain vs. dysplasique) a été remarqué auprès de la majorité des métriques évaluées. Ces résultats, ainsi que les facteurs qui les influencent, peuvent servir à assister la planification chirurgicale, effectuer des validations in-vitro d'optimisation d'ostéotomie péri-acétabulaire et, par la suite, guider la conception de nouveaux traitements. / Acetabular dysplasia causes articular instability and the reduction of the loading area in the hip. This causes the early onset of osteoarthritis which leads to pain and the eventual loss of articular function. Acetabular dysplasia is commonly treated by periacetabular osteotomy, but the precise effects of this correction on intra-articular biomechanics is not yet well known. Using discrete element analysis, pressure distribution has been calculated in the dysplastic hips of 9 patients before and after periacetabular osteotomy. A variety of radiographic measurements and numerical metrics have been evaluated in order to quantify the degree of acetabular correction and excessive loading. Correlations between these measures and metrics have been established in order to reveal any causal effects. In particular, the corrections of the lateral center-edge angle, acetabular index and femoral head extrusion index are correlated to the rise of average contact area during gait (p-value < 0.05, R² = 0.22 to 0.50). Additionally, the influence of cartilage layer reconstruction method (population-based healthy vs. dysplastic thicknesses) and loading parameters (healthy vs. dysplastic gait profiles) was significant for a majority of evaluated metrics. These results and influencing factors can help to assist with surgical planning, provide a basis for in-vitro validation of surgical optimization and act as a standard for designing new orthopedic treatments for acetabular dysplasia.

Les structures de sphères creuses appartiennent à la famille des matériaux cellulaires qui ont récemment été étudiés pour leurs multiples propriétés. Dans le cas de cette thèse, le but des sphères creuses est de dissiper l’énergie d’impact d’un oiseau sur un cockpit d’avion. Elles sont développées dans le cadre du projet SAMBA (Shock Absorber Material for Birdshield Application) afin d’optimiser leur énergie spécifique absorbée (J/kg).Dans un premier temps, des essais quasi-statiques (v = 5 mm/min) et dynamiques (v = 2 m/s) de compression uni-axiale sont conduits à température ambiante sur une seule sphère creuse de diamètre 30 mm. Une propagation rapide de fissures macroscopiques est observée. Le formalisme de la Mécanique Élastique Linéaire de la Rupture (MELR) est utilisé pour estimer le taux de restitution d’énergie critique dynamique GIdc du matériau constitutif. La position du sommet de fissure est mesurée pendant la propagation de fissure à l’aide d’une caméra rapide. La Méthode des Éléments Discrets (DEM) permet de simuler la rupture dynamique en implémentant une technique de relâchement des nœuds. Le taux de restitution d’énergie GIdc peut être estimé à partir de l’histoire (position et temps) du sommet de fissure. Le modèle numérique montre que les structures sphériques dissipent une proportion importante de l’énergie par des effets dynamiques. A une même vitesse de propagation, plus l’épaisseur de coque est fine, plus les effets inertiels générés par la rupture sont importants et ce pour une même vitesse de propagation.Le modèle numérique DEM est ensuite employé pour reproduire la rupture dynamique sur une sphère creuse à l’aide d’un critère en contrainte seule ou un critère mixte en contrainte – énergie. Les bons résultats obtenus démontrent la capacité de la DEM à représenter la propagation de fissures en régime dynamique.Finalement, des essais numériques et expérimentaux multi-sphères sont réalisés afin évaluer le comportement des sphères creuses au sein d’un assemblage. / Hollow sphere structure (HSS) belongs to cellular solids that have been studied recently for its multiples properties. In our case, HSS aims to absorb soft impacts energy on an airliner cockpit. HSS is investigated through the SAMBA (Shock Absorber Material for Bird-shield Application) project because of its promises in term of specific energy dissipated (J/kg) during impact.First of all, quasi-static and dynamic (v = 5 mm/min to v = 2 m/s) uniaxial compression tests are conducted at room temperature on a single sphere (D = 30 mm). Rapid crack propagation (RCP) is observed to be predominant at macroscopic scale. The formalism of Linear Elastic Fracture Mechanics (L.E.F.M.) is therefore used to estimate the dynamic energy release rate GIdc . The crack tip location is measured during the crack propagation using a high speed camera. The Discrete Element Method (DEM) is used to simulate the dynamic fracture by implementing the node release technique. The dynamic energy release rate can be determined using an experimentally measured crack history. In spherical structures the numerical results reveal a high proportion of energy dissipated through inertial effects as well as a dependence of the thickness of the hollow sphere over the range of 0.04 mm to 1.2 mm.The DEM model Is then employed to reproduce the RCP according to two failure criterions: a stress criterion and a coupled stress-energy criterion. It reveals to be an interesting way to model the mechanical behavior of brittle materials.Eventually, experimental and numerical multi-spheres tests are performed to evaluate the behavior of brittle hollow spheres within an assembly.