Le concassage entre trois appuis (multipoints) d'un assemblage cohésif de 5000 polyèdres convexes formant une tesselation de Voronoï cylindrique est simulé en 3D à l'aide du code de simulation numérique discrète (DEM) décrit dans l'article DESCANTES Y. (2022, https://hal.science/hal-03649287v1). La vidéo de gauche permet de suivre l'initiation et la propagation de fractures à la surface de l'éprouvette, tandis que la vidéo de droite montre la tesselation et en simultané les liens cohésifs qui se rompent (en rouge) au sein de l'assemblage au fil des sollicitations. Les grains étant rigides, les fractures sont initiées et se propagent nécessairement à leur périphérie.

Cette simulation s'inspire des simulations numériques réalisées à la suite de la thèse d'A. Neveu et décrites dans l'article Artoni et al (2019, https://hal.science/hal-01881801v1). Les calculs sont adimensionnés, de sorte que le diamètre de l'éprouvette vaut 1 et sa hauteur 0,5. Les appuis sont fixes, sauf l'appui du haut qui se déplace verticalement vers le bas à vitesse constante. Les polyèdres sont soumis à des forces cohésives visco-élastiques et à des forces de contact exercées par leurs voisins immédiats, par contre la gravité est nulle. Le frottement de Coulomb vaut ici 0. 

Le code ayant été parallélisé à l'aide d'OPENMP, le temps de calcul est de l'ordre de 10h en CPU-multithreading (~ 40 threads).