Gelöste Aufgaben/LM01

Aufgabenstellung

Wie Minimum-Prinzipe "funktionieren" und wie die Mathematik dazu aussieht, untersuchen wir hier an einem klassischen Beispiel.

Lageplan

Gesucht ist die Lage der Knoten von vier starren Kettengliedern (Länge ℓ, Masse m) im Erdschwerefeld.

Lösung mit Maxima

Lorem Ipsum ....

tmp

Die Aufgabe ist ein einfaches Beispiel für Lösungsansätze mit "Lagrange-Multiplikatoren".

Wir arbeiten mit dem Prinzip vom Minimum der Potentiellen Energie - also dem Potential der Gewichtskraft der Kettenglieder. Besonders charmant - im Vergleich zu Ansätzen mit dem Kräfte-Gleichgewicht ist hier das "Wegfallen" der Schnittkräfte - die brauchen wir hier nicht explizit angeben.

Header

Text

📑 toggle code

1+1

tmp

Koordinaten

Jeder Verbindungspunkt eines Kettengliedes hat zwei Koordinaten U_j, V_j.

Die Arbeitsfunktion der Gewichtskraft ist

$A=\sum _{n=0}^{N}mg\frac{U_{i-1}-U_i}{2}$.

Im Gleichgewicht  hat diese Arbeitsfunktion ein Minimum, allerdings müssen dabei für die starren Kettenglieder die Bedingungen

${(U_i-U_{i-1})}^2+{(V_i-V_{i-1})}^2={\ell }^2$

erfüllt sein!

Declarations

Text

📑 toggle code

1+1

tmp

Die Lagrange-Funktion für die Kette lautet also

$\begin{array}{lll}\mathrm{\Lambda }=& & \sum _{n=0}^{N}mg\frac{U_{i-1}-U_i}{2}\\ & +& \sum _{n=0}^N{\lambda }_i\cdot ({(U_i-U_{i-1})}^2+{(V_i-V_{i-1})}^2-{\ell }^2)\end{array}$

mit den Lagrangeschen Multiplikatoren λ_i.

Hier ist N =4, die Randbedingungen für die Verschiebungen lauten also

$U_0=0,U_4=0,V_0=0\text{ und }{V}_N=L$

und wir wählen hier

$L=3\ell$.

Mit der Abkürzung

$\gamma =mg\ell$ und $U_i=\ell {\tilde{U}}_i,V_i=\ell {\tilde{V}}_i,$

ist dann

$\begin{array}{lll}\mathrm{\Lambda }=& & \frac{({\tilde{U}}_3+{\tilde{U}}_2)\gamma }{2}+\frac{{\tilde{U}}_3\gamma }{2}+\frac{({\tilde{U}}_2+{\tilde{U}}_1)\gamma }{2}+\frac{{\tilde{U}}_1\gamma }{2}\\ & +& ({(3-{\tilde{V}}_3)}^2+{\tilde{U}}_3^2-1){\lambda }_4\\ & +& ({({\tilde{V}}_3-{\tilde{V}}_2)}^2+{({\tilde{U}}_3-{\tilde{U}}_2)}^2-1){\lambda }_3\\ & +& ({({\tilde{V}}_2-{\tilde{V}}_1)}^2+{({\tilde{U}}_2-{\tilde{U}}_1)}^2-1){\lambda }_2\\ & +& ({\tilde{V}}_1^2+{\tilde{U}}_1^2-1){\lambda }_1\end{array}$.

Die Gleichgewichtsbedingungen unter Nebenbedingungen lauten dann

$\begin{array}{ll}\frac{\partial \mathrm{\Lambda }}{\partial U_i}& =0\\ \frac{\partial \mathrm{\Lambda }}{\partial V_i}& =0\\ \frac{\partial \mathrm{\Lambda }}{\partial {\lambda }_i}& =0\end{array}$.

Equilibrium Conditions

Text

📑 toggle code

1+1

tmp

Die Gleichungen, die wir nun lösen müssen sind:

$\begin{array}{ll}0=& -2{\tilde{U}}_2{\lambda }_2+2{\tilde{U}}_1{\lambda }_2+2{\tilde{U}}_1{\lambda }_1+1\\ 0=& -2{\tilde{V}}_2{\lambda }_2+2{\tilde{V}}_1{\lambda }_2+2{\tilde{V}}_1{\lambda }_1\\ 0=& {\tilde{V}}_1^2+{\tilde{U}}_1^2-1\\ 0=& -2{\tilde{U}}_3{\lambda }_3+2{\tilde{U}}_2{\lambda }_3+2{\tilde{U}}_2{\lambda }_2-2{\tilde{U}}_1{\lambda }_2+1\\ 0=& -2{\tilde{V}}_3{\lambda }_3+2{\tilde{V}}_2{\lambda }_3+2{\tilde{V}}_2{\lambda }_2-2{\tilde{V}}_1{\lambda }_2\\ 0=& {\tilde{V}}_2^2-2{\tilde{V}}_1{\tilde{V}}_2+{\tilde{U}}_2^2-2{\tilde{U}}_1{\tilde{U}}_2+{\tilde{V}}_1^2+{\tilde{U}}_1^2-1\\ 0=& 2{\tilde{U}}_3{\lambda }_4+2{\tilde{U}}_3{\lambda }_3-2{\tilde{U}}_2{\lambda }_3+1\\ 0=& 2{\tilde{V}}_3{\lambda }_4-6{\lambda }_4+2{\tilde{V}}_3{\lambda }_3-2{\tilde{V}}_2{\lambda }_3\\ 0=& {\tilde{V}}_3^2-2{\tilde{V}}_2{\tilde{V}}_3+{\tilde{U}}_3^2-2{\tilde{U}}_2{\tilde{U}}_3+{\tilde{V}}_2^2+{\tilde{U}}_2^2-1\\ 0=& {\tilde{V}}_3^2-6{\tilde{V}}_3+{\tilde{U}}_3^2+8\end{array}$.

Und die sind nichtlinear! Dafür brauchen wir eine Lösungsroutine - wir verwenden das Newton–Verfahren.

Wir finden:

$\begin{array}{lll}{\tilde{U}}_1=& & 0.807\\ {\tilde{V}}_1=& & 0.59\\ {\lambda }_1=& -& 0.929\\ {\tilde{U}}_2=& & 1.22\\ {\tilde{V}}_2=& & 1.5\\ {\lambda }_2=& -& 0.603\\ {\tilde{U}}_3=& & 0.807\\ {\tilde{V}}_3=& & 2.41\\ {\lambda }_3=& -& 0.603\\ {\lambda }_4=& -& 0.929\end{array}$.

✔ Lösungen nichtlinearer Gleichungen:

... sind eine Wissenschaft für sich. Selbst für scheinbar triviale Probleme wie dieses laufen wir für N>4 schon in große numerische Probleme hinein. Der Trick ist häufig, als Startpunkt der Iteration Werte vorzugeben, die einigermaßen in der Nähe der erwarteten Lösung liegen. Sonst verabschiedet sich der Newton-Algorithmus oft mit einem "convergance error". Oder wir finden Lösungen, die zwar Gleichgewichtsbedingungen sind aber offensichtlich instabil sind.

Solving

Text

📑 toggle code

1+1

tmp

Lösungskonstellation

Hier erhalten wir jedoch eine Lösung - und so sieht sie aus;

Post-Processing

Text

📑 toggle code

1+1

---

Links

• ...

Literature

• ...