Mathématiquement, pourquoi l'agrégation de matrice de masse / vecteur de charge fonctionne-t-elle?

Je sais que les gens remplacent souvent des matrices de masse cohérentes par des matrices diagonales groupées. Dans le passé, j'ai également implémenté un code dans lequel le vecteur de charge est assemblé de manière groupée plutôt que de manière cohérente FEM. Mais je n'ai jamais cherché à savoir pourquoi nous sommes autorisés à le faire en premier lieu.

Quelle est l'intuition derrière l'agrégation qui permet de l'appliquer aux vecteurs de masse et de charge? Quelle en est la justification mathématique? Dans quelles situations l'agrégation n'est-elle pas autorisée / n'est-elle pas une bonne approximation pour les vecteurs de masse et de charge?

Dans la méthode des éléments finis, les entrées de matrice et les entrées de droite sont définies comme intégrales. Nous ne pouvons, en général, pas les calculer exactement et appliquer la quadrature. Mais il existe de nombreuses formules de quadrature que l'on pourrait choisir, et on les choisit souvent de manière à ce que (i) l'erreur introduite par la quadrature soit du même ordre que celle due à la discrétisation, ou du moins pas sensiblement pire, et (ii) la matrice possède certaines propriétés qui s'avèrent pratiques.

Le regroupement des masses en est un exemple: si l'on choisit une formule de quadrature particulière (à savoir celle avec des points de quadrature situés aux points d'interpolation de l'élément fini), alors la matrice de masse résultante se trouve être diagonale. C'est assez pratique pour l'implémentation informatique, et la raison pour laquelle les gens utilisent ces formules en quadrature. C'est aussi la raison pour laquelle cela "fonctionne": ce choix particulier de formule en quadrature a toujours un ordre relativement élevé.

Les matrices diagonales ont des avantages évidents pour accélérer les calculs numériques, et la réponse de Wolfgang Bangerth est une bonne explication sur la façon de calculer une matrice de masse diagonale, mais elle ne répond pas à la question du PO "pourquoi ça marche " au sens de "pourquoi c'est une bonne approximation de la physique que vous modélisez ".

Conceptuellement, vous pouvez séparer la réponse d'un élément en trois parties: le mouvement de translation d'un corps rigide, la rotation rigide autour du centre de masse de l'élément et la déformation de l'élément.

La fonction de base de la matrice de masse de l'élément est de représenter l'élément KE sous une forme quadratique (c'est-à-dire où sont les vitesses nodales).$\frac 1 2 v^T M v$$v$

À mesure que la taille de l'élément diminue, la contribution à la KE de la rotation du corps rigide diminue plus rapidement que la contribution de la translation (pour un élément solide avec une taille linéaire typique de , la masse est proportionnelle à mais les moments d'inertie sont proportionnelle à ) et la contribution de la déformation de l'élément est négligeable (au moins pour les problèmes de petites déformations élastiques).$a$$a^3$$a^5$

Par conséquent, vous n'avez vraiment besoin que d'une "bonne" approximation des parties du corps rigides du mouvement, c'est-à-dire 6 DOF, et en fait une bonne approximation uniquement du KE de la translation du corps rigide , c'est-à-dire 3 DOF, convergera lorsque la taille de l'élément est réduit.

Les termes diagonaux de la matrice d'éléments contiennent plus qu'assez de paramètres indépendants pour représenter ces 3 ou 6 termes KE avec une précision suffisante. En fait, pour les éléments d'ordre supérieur, vous pouvez utiliser des matrices de masse de masse diagonale où les termes diagonaux pour les nœuds médians sont nuls.

Notez qu'il s'agit d'une situation complètement différente de l'énergie potentielle de l'élément, où les contributions de la translation et de la rotation du corps rigide sont nulles, et la seule chose qui compte est de représenter l'énergie de déformation correspondant à la déformation de l'élément . Une matrice de rigidité diagonale ne serait donc pas une idée réalisable!

En plus des autres réponses, il existe des scénarios dans lesquels les erreurs dans la matrice de masse n'ont aucune influence sur le résultat souhaité.

Étant donné un problème de déformation non linéaire de la forme qui a une solution unique , on pourrait envisager de résoudre ce problème comme le problème dynamique avec la matrice de masse et la matrice d'amortissement grâce à une discrétisation temporelle appropriée. Il est facile de voir que le problème d'origine est résolu lorsqu'un état de repos a été atteint ( ). Surtout, n'a aucune influence sur ce résultat (tant que la convergence vers le résultat unique est atteinte).$K(u) \space u = f(u)$$\hat{u}$$K(u) \space u + C(u) \space \dot{u} + M \space \ddot{u} = f(u)$$M$$C$$\dot{u} = \ddot{u} = 0$$M$

La justification de l'agrégation de masse a alors plus à voir avec la vitesse de convergence plutôt qu'avec, par exemple, l'erreur de quadrature ¹ : la facilité d'inverser une diagonale traduit par des itérations de pas de temps rapides avec un schéma d'intégration de temps approprié (où est le seul matrice inverse).$M$$M^{-1}$

^{1 Bien que le raisonnement sur le comportement physique dynamique soit bien sûr plus facile avec une matrice de masse «correcte» - par exemple, le moment angulaire peut être incorrectement conservé par des matrices de masse groupées.}