Mathematisch gesehen, warum funktioniert das Zusammenfassen von Massenmatrix und Lastvektor?


13

Ich weiß, dass Menschen häufig konsistente Massenmatrizen durch gewürfelte Diagonalmatrizen ersetzen. In der Vergangenheit habe ich auch einen Code implementiert, bei dem der Lastvektor nicht auf FEM-konsequente Weise, sondern auf konzentrierte Weise zusammengesetzt wird. Aber ich habe nie darüber nachgedacht, warum wir das überhaupt dürfen.

Was ist die Intuition hinter Klumpen, die es einem erlaubt, sie auf Masse- und Lastvektoren anzuwenden? Was ist die mathematische Rechtfertigung dafür? In welchen Situationen ist das Zusammenfassen nicht erlaubt / keine gute Annäherung für Masse- und Lastvektoren?

Antworten:


16

Bei der Finite-Elemente-Methode werden die Matrixeinträge und die Einträge auf der rechten Seite als Integrale definiert. Wir können diese im Allgemeinen nicht genau berechnen und Quadratur anwenden. Es gibt jedoch viele Quadraturformeln, die man wählen kann, und man wählt sie oft so, dass (i) der durch die Quadratur verursachte Fehler in der gleichen Größenordnung wie der aufgrund der Diskretisierung oder zumindest nicht wesentlich schlechter ist, und (ii) Die Matrix hat bestimmte Eigenschaften, die sich als zweckmäßig herausstellen.

Ein Beispiel für diese Arbeitsweise ist das Massenklumpen: Wenn man eine bestimmte Quadraturformel wählt (nämlich die mit Quadraturpunkten an den Interpolationspunkten des finiten Elements), ist die resultierende Massenmatrix zufällig diagonal. Dies ist sehr praktisch für die rechnerische Implementierung und der Grund, warum diese Quadraturformeln verwendet werden. Es ist auch der Grund, warum es "funktioniert": Diese spezielle Wahl der Quadraturformel hat immer noch eine einigermaßen hohe Ordnung.


Geniale Antwort, wie immer. Mich würde auch Ihre Meinung zum zweiten Teil der Frage sehr interessieren, wenn eine Verdichtung nicht zulässig ist / eine schlechte Annäherung vorliegt , wenn mir etwas in den Sinn kommt.
Anton Menshov

2
@AntonMenshov: Es scheint schwierig (vielleicht unmöglich?) Zu sein, eine gute Annäherung über das Zusammenfassen von Elementen höherer Ordnung zu erhalten, da in diesem Fall das Zusammenfassen von Elementen höherer Ordnung (z. B. diagonal) einer Quadratur niedrigerer Ordnung entspricht Polynome.
Paul

@WolfgangBangerth: Ich glaube ich verstehe jetzt. Es ist also so, als würde man Newton-Cotes-Regeln für die Integration anstelle der Gaußschen Quadratur verwenden. Da jede Lagrange-Interpolationsfunktion Einheitswerte an einem bestimmten Knoten hat, führt das Migrieren der Quadraturpunkte zu den Knoten dazu, dass nur die diagonalen Terme ungleich Null werden (zumindest für lineare Elemente).
Paul

1
2n3

1
Der wichtige Punkt ist, dass für Elemente höherer Ordnung die "konzentrierte" Massenmatrix über bestimmte Quadraturformeln definiert werden muss. Die ursprüngliche Form - von der der Begriff "zusammengefasst" stammt - fügte die Einträge außerhalb der Diagonale zur Diagonale hinzu, aber das funktioniert nur, wenn sie alle positiv sind. Wenn Sie die Gaußsche Quadratur anwenden, gilt dies für Elemente niedrigster Ordnung, jedoch nicht für Elemente höherer Ordnung.
Wolfgang Bangerth

5

Diagonale Matrizen haben offensichtliche Vorteile bei der Beschleunigung numerischer Berechnungen, und die Antwort von Wolfgang Bangerth ist eine gute Erklärung für die Berechnung einer diagonalen Massenmatrix, sie beantwortet jedoch nicht die Frage des OP "Warum funktioniert das ? " Im Sinne von "Warum ist das?" es ist eine gute Annäherung an die Physik, die Sie modellieren ".

Konzeptionell können Sie die Reaktion eines Elements in drei Teile unterteilen: Translationsbewegung eines starren Körpers, starre Drehung um den Massenmittelpunkt des Elements und Verformung des Elements.

12vTMvv

einein3ein5

Daher brauchen Sie wirklich nur eine "gute" Annäherung an die starren Körperteile der Bewegung, dh 6 DOFs, und tatsächlich konvergiert eine gute Annäherung an nur die KE aus der Starrkörper- Translation , dh 3 DOFs, wie die Elementgröße ist reduziert.

Die diagonalen Terme der Elementmatrix enthalten mehr als genug unabhängige Parameter, um diese 3 oder 6 KE-Terme mit ausreichender Genauigkeit darzustellen. Tatsächlich können Sie für Elemente höherer Ordnung Massendiagonale-Massenmatrizen verwenden, bei denen die diagonalen Terme für die Mittelknoten Null sind.

Es ist zu beachten, dass dies eine völlig andere Situation ist als die potentielle Energie des Elements, bei der die Beiträge der Translation und Rotation des starren Körpers Null sind und das Einzige, was zählt, die Dehnungsenergie darstellt, die der Deformation des Elements entspricht . Eine diagonale Steifigkeitsmatrix wäre daher keine realisierbare Idee!


5

Neben den anderen Antworten gibt es Szenarien, in denen Fehler in der Massenmatrix keinen Einfluss auf das gewünschte Ergebnis haben.

K(u) u=f(u)u^K(u) u+C(u) u˙+M u¨=f(u)MCu˙=u¨=0M

MM-1

1 Obwohl das Denken über dynamisches physikalisches Verhalten mit einer "richtigen" Massenmatrix natürlich einfacher ist - z. B. kann der Drehimpuls durch konzentrierte Massenmatrizen nicht ordnungsgemäß erhalten werden.


1
Cu˙Cu˙Mu¨Cu˙C

u¨
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.