Gelöste Aufgaben/Kw26

Aufgabenstellung

Ein Aggregat der Masse 2m ist durch zwei Gummifedern wie skizziert gelagert und wird senkrecht durch seine Gewichtskraft und eine äußere Last F=2 H sin(Ω t) belastet.

Für dieses System mit einem nichtlineares Federelement soll die Lösung als Anfangswertproblem gefunden und mit der linearen Näherungslösung verglichen werden.

Gesucht sind

die Lösung der zugeordneten linearisierten Bewegungsgleichung,
die numerische Lösung des nichtlinearen Problems.

Für die Gummifeder gibt der Lieferant folgende eine Kennlinie an:

Gegeben:

u₀,
m, g, H, Ω'

Lösung mit Maxima

Header

Gummi-Elemente zeigen eine deutliche Nichtlinearität bei großen Auslenkungen. Wie man diese Nichtlinearität bei dynamischer Belastung berücksichtigt, zeige ich Ihnen hier.

Declarations

Wir berücksichtigen nur die vertikale Bewegung des Aggregats, dafür nehmen wir die Koordinate u:

Als System-Parameter setzten wir an:

$H = m g γ$
(Kraft-Amplitude als Vielfaches der Gewichtskraft) und
$\begin{array}{l} Ω = \frac{2 π}{T_{Ω}}, \\ T_{Ω} = η T_{0} \end{array}$
(Erreger-Periodendauer ist ein Vielfaches der Schwingungsdauer des zugeordneten linearen Systems).

Spring Characteristic

Die Gummi-Element Federkennlinie approximieren wir über die punktsymmetrische Funktion

\frac{K (u_{t})}{m g} = c_{1} \cdot (\frac{u_{t}}{u_{0}}) + c_{3} \cdot {(\frac{u_{t}}{u_{0}})}^{3}

Die beiden Konstanten c1 und c2 erhalten wir aus den Bedingungen

\frac{K (u_{0})}{m g} = 1 und \frac{K (2 u_{0})}{m g} = 3

,

die wir aus der Kennlinie ablesen zu

c_{1} = \frac{5}{6}, c_{3} = \frac{1}{6}

So sieht die Kennlinie dann aus:

Equilibrium Conditions

Die Bewegungsgleichungen des Systems bekommen wir direkt aus dem Kräftegleichgewicht am Körper (mit dem Prinzip von d'Alembert) zu

2 m \ddot{u} + 2 K (u_{t}) = 2 m g + F (t)

Damit die "homogene Lösung" der Bewegung (der Einschwingvorgang in die periodische Lösung) nach endlicher Zeit weg ist, mogeln wir noch etwas Dämpfung hinzu:

2 m \ddot{u} + 2 b \dot{u} + 2 K (u_{t}) = 2 m g + F (t)

.

Wir brauchen auch die linearisierte Bewegungsgleichung. Dazu entwickeln wir die Federkennlinie um die stationäre Ruhelage u = u₀ des Aggregats in eine Taylor-Reihe bis zum linearen Glied, setzen also

u = u_{0} + u_{l} (t) mit u_{l} ≪ u_{0}

und finden die linearisierte Form der Bewegungsgleichung zu

2 m \cdot \ddot{u_{l}} + 2 \cdot b \cdot {\dot{u}}_{l} + \frac{8 \cdot m g}{3 \cdot u_{0}} \cdot u_{l} = F (t)

.

Dimensionless Representation

Dimensionslos machen wir die Bewegungsgleichung mit der Referenz-Auslenkung u₀ und der Periodendauer T₀ der zugeordneten linearen Bewegungsgleichung, also

$u = μ \cdot u_{0}$ und
$t = τ \cdot T_{0}$

Umschreiben der linearen und nichtlinearen Bewegungsgleichungen in dimensionslose Form liefert mit

b = 2 D m ω_{0},

die nichtlineare Form

$\frac{d^{2}}{d τ^{2}} \cdot μ + 4 \cdot π \cdot D \cdot (\frac{d}{d τ} \cdot μ) + \frac{π^{2}}{2} \cdot (5 μ + μ^{3}) = 3 \cdot π^{2} \cdot (1 + γ s i n (\frac{2 \cdot π \cdot τ}{η}))$ und

die linearisierte Form

$\frac{d^{2}}{d τ^{2}} \cdot μ_{l} + 4 D \cdot π \cdot (\frac{d}{d τ} \cdot μ_{l}) + 4 \cdot π^{2} \cdot μ_{l} = 3 \cdot π^{2} \cdot s i n (\frac{2 \cdot π \cdot τ}{η}) \cdot γ$ .

Solving

Teil A: Partikulare Lösung der linearen Bewegungsgleichung

Von der linearen Bewegungsgleichung wollen wir nur die partikulare Lösung - also die Lösung nach Abklingen des Einschwingvorgangs (dafür brauchten wir die Dämpfung!). Wir wählen als Ansatz

μ_{l} (τ) = U_{s} \cdot s i n (\frac{2 \cdot π \cdot τ}{η}) + U_{c} \cdot c o s (\frac{2 \cdot π \cdot τ}{η})

(Ansatz vom Typ der rechten Seite)

und erhalten für die Koeffizienten

\begin{array}{l} U_{c} = - \frac{3 \cdot η^{3} \cdot D \cdot γ}{8 \cdot η^{2} \cdot D^{2} + 2 \cdot η^{4} - 4 \cdot η^{2} + 2}, \\ U_{s} = \frac{3 \cdot η^{4} \cdot γ - 3 \cdot η^{2} \cdot γ}{16 \cdot η^{2} \cdot D^{2} + 4 \cdot η^{4} - 8 \cdot η^{2} + 4} \end{array}

.

Die Lösung plotten wir unten.

Teil B: Numerische Lösung des Anfangswertproblems

Für die numerische Lösung brauchen die Bewegungsgleichung als Differentialgleichung erster Ordnung als

\begin{array}{ll} \frac{d μ}{d τ} & = ν, \\ \frac{d ν}{d τ} & = - 4 \cdot π \cdot D \cdot ν - \frac{π^{2}}{2} (5 μ + μ^{3}) + 3 π^{2} (1 + γ \sin (\frac{2 \cdot π \cdot t}{η})) \end{array}

.

Die lösen wir nun als Anfangswertproblem

\underline{\dot{q}} = \underline{f} (\underline{q})

mit

\underline{q} = (\begin{array}{c} μ \\ ν \end{array})

.

Als Anfangsbedingungen wählen wir die stationäre Ruhelage des Aggregats, also

{\underline{q}}_{0, 1} = (\begin{array}{l} 1 \\ 0 \end{array})

.

Und wir führen die numerische Integration mit drei Parameter-Sätzen durch:

Überkritische Anregung:

$η = 0.5, γ = 1, D = 0.1$ Hier ist η = 1/2, also erfolgt die Anregung mit einer Periodendauer, die halb so groß ist, wie die Periodendauer des linearisierten Systems - oder: die Erregerfrequenz ist doppelt so hoch wie die Eigenfreuqenz.

Anregung in der Nähe der Resonanzfrequenz:

$η = 1.1, γ = 1, D = 0.1$ Mit η = 1.1 regen wir in der Nähe der Erregerfrequenz des linearen Systems an - etwas darunter, die Phase des Systems ist also fast noch konphas mit der Erregung.

Unterkritische Anregung:

$η = 2.0, γ = 1, D = 0.1$ Hier ist Erregerfrequenz nur hald so groß, wie die Eigenfrequenz des linearen Systems, die Phase des Systems ist also konphas zur Erregung.

Post-Processing

Für die dreiAnfangsbedingungen erhalten wir die Auslenkungen U(τ):

#	Auslenkung	Geschwindigkeit
1)
2)
3)

Für das nichtlineare System können wir uns den Einschwingvorgang im Phasendiagramm ansehen:

Und stellen fest:

nur für die überkritische Lösung ist die Auslenkungen des Systems so klein, dass eine Linearisierung gerechtfertigt ist.

Für die Parameter-Sätze 3 und besonders 2 entstehen Schwingungen mit hohen Amplituden, die wir in jedem Fall als nichtlineares Problem lösen müssen.

Und im eingeschwungenen Zustand können wir die Lösung für lineares und nichtlineares System einfach vergleichen. Die Lösungen des linearen Systems sind immer Ellipsen, die Lösungen des nichtlinearen Systems sind deutlich "kunstvoller":

🧨 {{{title}}}:

{{{text}}}

Links

...

Literature

...

Gelöste Aufgaben/Kw26

Inhaltsverzeichnis

Aufgabenstellung