Masse-Feder-System (Computergrafik)

Das Masse-Feder-System (kurz MFS, englisch Mass-Spring-System), in der Literatur auch Masse-Feder-Modell (kurz MFM), ist ein Verfahren in der Computergrafik, um das Deformationsverhalten von auf Polygonnetzen basierenden Objekten zu bestimmen. Es kann sowohl auf ebene als auch auf räumliche Objekte angewandt werden.

Die Gründe für die weite Verbreitung des Masse-Feder-Systems sind das leichtverständliche physikalische Konzept, die einfache Implementierung und die relativ niedrigen Rechenanforderungen.

Geschichte

Das Masse-Feder-System basiert auf den physikalischen Ansätzen des Hookeschen Gesetzes.

Wegen der geringen Rechenleistung damaliger Computer wurde es anfänglich nur für Flächen genutzt (sowohl zweidimensionale, als auch im dreidimensionalen Raum eingebettete Flächen, zum Beispiel ein Tischtuch, das über einen Tisch fällt). Deshalb wurde es bei seinen ersten Verwendungen ab dem Jahr 1981 zunächst nur zur Haut- und Stoffverhaltenssimulation verwendet. Erst 1989 wurden die ersten MFS-Volumenmodelle veröffentlicht.

Heute ist es ein gängiges Verfahren in der Computerspielbranche und in der Medizin zur Simulation operativer Eingriffe.

Aufbau

Zwei Masseknoten, die mit einer Feder und einem Dämpfer verbunden sind.

In einem Masse-Feder-System werden Objekte durch ein Netz aus Masseknoten, mechanischen Federn und Dämpfungsgliedern dargestellt.

Ablauf der Verschiebung eines Punktes in einem Masse-Feder-System zu verschiedenen Zeitpunkten.

Die Knoten eines Polygonnetzes bekommen jeweils eine Masse. Die Kanten werden durch eine mechanische Feder, die mit einem Dämpfungsglied parallel geschaltet ist, dargestellt. Somit hat sie eine Federkonstante, eine Ruhelänge und eine Abklingkonstante. Das Dämpfungsglied ist notwendig, wenn das Objekt nach der Deformation nicht schwingen, sondern eine Ruhelage finden soll. In Abhängigkeit von den Masseknoten wird der Reibungskoeffizient errechnet.

Die Masseknoten sind regelmäßig im ganzen Objekt verteilt und stellen es so dar. Zur Verbindung der Masseknoten werden Federn zwischen sie gesetzt. Sind zwei Masseknoten miteinander verbunden, nennt man sie benachbart. Die Ausgangsform des Objekts liegt vor, wenn alle Federn in Ruhelage sind ($E_{\mathrm {pot} }=0$), also deren Ruhelänge haben. Wird nun ein Masseknoten verschoben, liegt an den benachbarten Federn eine potentielle Energie an ($E_{\mathrm {pot} }\neq 0$), die nun durch Verschieben der benachbarten Massepunkten minimiert wird. Eine solche Verschiebung wirkt sich auf das ganze Netz aus, wodurch sich dieses verschiebt.

Um einer Verschiebung im Raum entgegenzuwirken, können Punkte „befestigt“ werden (an einer absoluten Position oder an einer relativen Position zu einem anderen Massepunkt). Dadurch kann unter anderem eine plastische Deformation erzielt werden.

Das Masse-Feder-System kann, um physikalisches Verhalten realistischer zu machen, um eine Kollisionserkennung ergänzt werden, damit die Massepunkte auf Kollisionen untereinander und auf Kollisionen mit anderen Objekten reagieren können.

Berechnung

In der folgenden Herleitung wird von einem $d$-dimensionalen Raum ausgegangen.

Masse-Feder-Systeme haben ihren Ursprung in der klassischen Mechanik. Sie basieren auf der Idee, dass ein Körper $M$ durch einen einzelnen Masseknoten $m$ repräsentiert werden kann. Weitergehend kann dieser Masseknoten in $n$ kleinere Masseknoten $m_{i}$, verteilt über den ganzen Körper, unterteilt werden, sodass gilt:

$m=\sum \limits _{i=1}^{n}m_{i}$

Nach dem Hookeschen Gesetz ist das Verhalten einer Feder mit parallelem Dämpfungsglied folgendermaßen definiert:

${\vec {F}}_{s}=-k\cdot {\vec {s}}-\mu \cdot {\dot {\vec {s}}}$

wobei $k$ die Federkonstante, ${\vec {s}}$ die Auslenkung der Feder zur Ruhelage und $\mu$ der Reibungskoeffizient ist.

Newtons Definition zum Verhalten einer bewegten Masse ist:

${\vec {F}}=m\cdot {\vec {a}}$

${\vec {F}}$ ist die Kraft, die dem Körper der Masse $m$ die Beschleunigung ${\vec {a}}$ erteilt.

Auf Basis dieser letzten beiden Gleichungen kann die Bewegungsgleichung eines einzelnen Masseknotens $i$ hergeleitet werden.

Auf den Masseknoten $i$ mit der Masse $m_{i}$ wird bei der Verschiebung um ${\vec {s}}_{i}\in \mathbb {R} ^{d}$ die Kraft ${\vec {F}}_{i}$ ausgeübt:

${\vec {F}}_{i}=m_{i}\cdot {\ddot {{\vec {s}}_{i}}}$

Die Auslenkung ${\vec {s}}_{ij}$ wird bei zwei Masseknoten $i$ und $j$ in einem Masse-Feder-System mit der Distanz voneinander, abzüglich der Länge der Feder im Ruhezustand $l_{0,ij}$, berechnet:

${\vec {s}}_{ij}={\frac {{\vec {x}}_{j}-{\vec {x}}_{i}}{\begin{vmatrix}{\vec {x}}_{j}-{\vec {x}}_{i}\end{vmatrix}}}\cdot ({\begin{vmatrix}{\vec {x}}_{j}-{\vec {x}}_{i}\end{vmatrix}}-l_{0,ij})$

${\vec {x}}_{i}$ ist die Position des $i$-ten Masseknotens. Mit ${\frac {{\vec {x}}_{j}-{\vec {x}}_{i}}{\begin{vmatrix}{\vec {x}}_{j}-{\vec {x}}_{i}\end{vmatrix}}}$ wird der Richtungsvektor als Einheitsvektor von Masseknoten $i$ zu Masseknoten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle j} berechnet. Der Betrag von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \begin{vmatrix}\vec x_j - \vec x_i\end{vmatrix} - l_{0,ij}} ist der Abstand zwischen Masseknoten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle i} und Masseknoten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle j} minus der Ruhelänge der Feder, was den Längenskalar darstellt. Das Skalarprodukt Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec s_{ij}} ist die Auslenkung des Masseknotens Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle i} in Abhängigkeit vom Masseknoten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle j} .

Das führt zu folgender Kraft, die auf den einzelnen Masseknoten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle i} mit dem benachbarten Knoten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle j} wirkt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_{s,ij} = -k_{ij} \cdot \vec s_{ij} - \mu_{ij} \dot{\vec s_i} }

Die externen Kräfte, die auf den Masseknoten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle i} wirken, können als Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_{\mathrm{ext},i}} zusammengefasst werden:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_{\mathrm{ext}, i} = \vec F_{G, i} + \vec F_{e, i} }

mit dem Vektor der äußeren Kraft Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_{\mathrm{ext}, i}} , dem Gravitationsvektor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_{G, i}} und dem Vektor der anderen externen Einwirkungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_{e, i}} .

Fasst man alle auf den Masseknoten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle i} wirkenden Kräfte zusammen, erhält man folgende Differenzialgleichung als Bewegungsgleichung, wobei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle N} die Menge der benachbarten Masseknoten ist:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_{i} = \sum \limits_{j \in N} \vec F_{s,ij} + \vec F_{\mathrm{ext}, i}}

Es gilt Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_i = m_i \ddot{\vec s_i}} :

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle m_i \ddot{\vec s_i} = \sum \limits_{j \in N}\vec F_{s, ij} + \vec F_{\mathrm{ext}, i} }

Die Bewegungsgleichung für das gesamte System erhält man, indem man alle Verschiebungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle s_i} mit dem Vektor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec s = (\vec s_1, \vec s_2, \dots, \vec s_n)^T} , die Massen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle m_i} mit der Matrix Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle M} , die Reibungskoeffizienten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \mu_{ij}} mit der Matrix Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle D} und die Federkonstanten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle k_{ij}} der einzelnen Masseknoten durch die Matrix Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K} beschreibt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec F_\mathrm{ext} = M \cdot \ddot{\vec s} + D \cdot \dot{\vec s} + K \cdot \vec s }

Die Matrizen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle D} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K} sind beide symmetrisch, da Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \mu_{ij} = \mu_{ji}} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle k_{ij} = k_{ji}} . Daraus folgt, dass jeweils nicht alle Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle n^2} Komponenten, sondern nur Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{n^2+3n}{2}} Komponenten berechnet werden müssen. Außerdem können die Komponenten auf der Hauptdiagonalen der Matrizen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle D} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K} (also die Komponenten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \mu_{ii}} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle k_{ii}} ) vernachlässigt werden, da es keine Feder gibt, die an beiden Enden mit demselben Masseknoten verbunden ist. Die Anzahl der zu berechnenden Komponenten verringert sich somit auf Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{n^2+n}{2}} .

Siehe auch

• Federpendel
• Harmonischer Oszillator

Weblinks

Commons: Masse-Feder-Systeme – Sammlung von Bildern, Videos und Audiodateien

• Johan Jansson: Beispiele zum Masse-Feder-System