Méthode des éléments discrets

Laméthode des éléments discretsest une méthode crée en 1979 par Cundallpour résoudreles équations de la dynamiquedans un milieu granulaire. Cette méthode consiste à définir le domaine à l’aide de particules liées entre elles par des lois d’interaction.
Cette méthode a été développée après pour résoudre les problèmes de discontinuité dans un milieu continu, et, plus précisément, pour visualiser la propagation de failles. Les scientifiques ont souhaité développer cette méthode parce que la méthode la plus connue,la méthode des éléments finis, atteint ses limites pour étudier des discontinuités.
Plusieurs méthodes ont été générées par la méthode principale de Cundall, elles sont classées en trois familles :
Consiste à modéliser le domaine par des poutres, ici des poutres d’Euler-Bernouilli, associées à des lois élastiques. Cette méthode donne des bons résultats pour visualiser l'amorçage et la propagation de fissures.
Cette méthode est une méthode de Cundall modifiée : la matière est modélisée par des sphères qui représentent des points de masse, liées entre elles par des ressorts. Les particules de forme sphériques ou circulaires présentent l’avantage de permettre un calcul très simple des contacts. Cette méthode est très efficace pour modéliser lecomportement élastique des matériaux.
C’est une combinaison entre les deux modèles particulaires et lattice, permettant de conjuguer les avantages des deux. Il est construit par des sphères reliées entre elles par des poutres[FY1]
La méthode des éléments discrets permet de modéliser un ou plusieurs solides rigides ou déformables en interaction par contact. Il est alors possible de déterminer les déplacements et les rotations de chacun des solides au cours de la simulation.
Il est également possible d'utiliser cette méthode pour évaluer l'endommagementde solideshétérogènes, c'est-à-dire l'amorçage et la propagation de fissures dans le milieu. Pour cela, un critère de rupture est imposé au niveau des liaisons cohésives qui relient les éléments discrets entre eux. Si le critère est dépassé, la liaison est rompue et n'a plus d'incidence dans le reste de la simulation.
Enfin, cette méthode peut aussi être utilisée pour simuler plusieurs types d'écoulement granulaire. Elle permet une étude détaillée de la dynamique des écoulement à l'échellemicroscopique, ce qui est souvent difficile voire impossible à faire expérimentalement.
Cependant, cette méthode a le défaut de nécessiter de grands temps de calcul et une puissance de calcul assez élevée. En effet, pour obtenir des résultats suffisamment précis, il est nécessaire de travailler un grand nombre d'éléments discrets pour modéliser le ou les solide(s). Mais plus le nombre d'éléments est élevé, plus les performances nécessaires et le temps de calcul augmentent. Ce dernier peut alors monter jusqu'à plusieurs heures, voire plusieurs jours. C'est la raison pour laquelle cette méthode n'est pas souvent utilisée dans l'industrie. Il est cependant possible de réduire ce temps en utilisant desprocesseurs graphiques, ou GPU, ou en réalisant la simulation sur des clusters de calculs.
L'utilisation de la méthode des éléments discrets nécessite de vérifier que cette approche permet bien de retrouver le comportement mécanique relatant de lamécanique des milieux continus. Il est alors nécessaire de déterminer et de quantifier les lois d'interactions à l'échelle microscopique, c'est-à-dire de paramétrer correctement les liaisons cohésives reliant les éléments discrets entre eux, en intégrant des paramètres microscopiques qui peuvent être évaluées grâce à desessais de traction/compressionsimple,flexionsimple etc. Le but est de trouver les paramètres qui permettent d'obtenir le même comportement que celui obtenu expérimentalement, ou avec d'autres méthodes comme la méthode deséléments finis.
Pour donner un exemple, dans le modèle de liaison hybride "lattice - particulaire" (ou poutre), les éléments discrets sont reliés entre eux par des poutres le plus souvent cylindriques, de typeBernouilli. Cela permet de décrire la géométrie de la poutre avec uniquement deux paramètres : la longueur de la poutre, et son rayon. On peut ensuite ajouter des paramètres mécaniques, comme lemodule de Young, lecoefficient de Poisson, la limite de rupture etc, et des paramètres dynamiques comme ladensitédes éléments discrets.