Quelltext der Seite Hüpfender Tennisball

Ein Tennisball wird aus einer bestimmten Höhe über einem Betonboden losgelassen. Die nachfolgende Bewegung soll mit Hilfe eines eindimensionalen Modells simuliert werden.

==Spezifikation==
[[Bild:HuepfenderTennisball.jpg|thumb|]]Die Tennisregeln der International Tennis Federation (ITF) schreiben für Bälle folgendes vor:
*Die äussere Hülle des Balles muss gleichförmig und nahtlos, seine Farbe weiss oder gelb sein.
*Der Ball muss diesen Anforderungen entsprechen und eine Masse von mehr als 56,7&nbsp;g und weniger als 58,5&nbsp;g aufweisen.
*Der Durchmesser muss mehr als 6,35&nbsp;cm und weniger als 6,67&nbsp;cm betragen.
*Es gibt mehrere spezifizierte Ballarten. Jeder Ball muss eine Sprunghöhe von mehr als 134,62&nbsp;cm und weniger als 147,32&nbsp;cm aufweisen, wenn er aus einer Höhe von 254&nbsp;cm auf eine ebene, harte Fläche, z.&nbsp;B. Beton, fallen gelassen wird.
Aus diesen Regeln werden folgende Spezifikationen abgeleitet
*Masse 0.058 kg
*Radius 0.0325 m
*Widerstandsbeiwert 0.4
*Stosszahl 0.55

==Modell==
[[Bild:HuepfenderTennisball_SD.jpg|thumb|Systemdiagramm]]Die [[Impulsbilanz]] bildet das Rückgrat dieses Modells zur eindimensionalen [[Translationsmechanik]]. Die positive Richtung zeige nach oben. Deshalb entzieht das Gravitationsfeld der Kugel mit der konstanten Rate ''F<sub>G</sub> = m*g'' Impuls. Die Luft entzieht der Kugel mit zunehmender Geschwindigkeit immer mehr Impuls. Ist die Geschwindigkeit negativ, führt sie Impuls zu. Um dieses Verhalten sauber zu modellieren, formulieren wir den [[Luftwiderstand]] mit dem Produkt aus der Geschwindigkeit und ihrem Betrag

::<math>F_W=\frac{\varrho_L}{2}c_WA\left|v\right|v</math>

Der statische [[Auftrieb]], der auch in Luft wirkt, ist hier nicht berücksichtigt bzw. mit der Masse des Balls verrechnet worden.

Die Geschwindigkeit eines Körper ist '''dynamisch''' durch den Quotienten aus [[Impuls]] und [[Masse]] festgelegt. Weil die Geschwindigkeit '''kinematisch''' als Änderungsrate des Ortes definiert ist, kann die Strecke mittels einer "Rohr-Topf-Anordnung" berechnet werden. In diesem Fall muss das Rohr (flow) in den Topf (Reservoir) hinein orientiert sein. Weil die positive Richtung nach oben zeigt, nimmt auch die Höhe bis zum Aufschlag ab.

Um den Aufschlag zu modellieren, wird zuerst eine positive Strecke (<tt>Delta_z</tt>)

:<tt>Delta_z = if z>0 then 0 else -z</tt>

Die eigentliche Aufschlagdynamik basiert auf einer starken Feder, deren Federkonstante beim Entspannen um die Stosszahl kleiner ist als beim Zusammendrücken

:<tt>FF = if v<0 then D*Delta_z else Stosszahl*D*Delta_z </tt>

==Simulation==
[[Bild:HuepfenderTennisball_h_v.png|thumb|Höhe und Geschwindigkeit]]Würde man den Tennisball im Vakuum fallen lassen, wäre die Beschleunigung betragsmässig gleich der Gravitationsfeldstärke. Die Geschwindigkeit des Balls würde jede Sekunde um 9.81 m/s abnehmen. Seine Endgeschwindigkeit wäre dann bei einer Fallhöhe von 2.54 m gleich

:<math>v_e=\sqrt{2gh}</math> = 7.06 m/s

was eine Fallzeit von

:<math>t=\frac{v_e}{g}</math> = 0.7196 s

ergibt. Die Simulation ergibt eine Endgeschwindigkeit von knapp 6.96 m/s bei einer Fallzeit von 0.724 s. Der Unterschied rührt von der Wirkung des Luftwiderstandes.

Weil der Boden als harte, teilelastische Feder (Federkonstante 10 kN/m) modelliert worden ist, muss der Zeitschritt der Simulation sehr klein gewählt werden. Nimmt man einen Zeitschritt von 10 Mikrosekunden (0.00001 s), müssen 500'000 Schritte ausgeführt werden, wozu ein Laptop aus dem Jahre 2007 0.28 Sekunden benötigt.

==Boden modellieren==
[[Bild:HuepfenderTennisball_FF_Dz.png|thumb|Federkennlinie]] Die Graphik zeigt das Kraft-Verformungs-Diagramm der teilelastischen Feder, mit welcher der Boden modelliert worden ist. Nimmt man die Hälfte der im ersten Aufschlag erzeugten Verformung von 0.016 m und dividiert diese 8 mm durch die Aufprallgeschwindigkeit von 7 m/s, erhält man für die Bremszeit etwa 1.1 Millisekunden. Diese Abschätzung zeigt, dass der Zeitschritt bis etwa auf eine Millisekunde angehoben werden darf, ohne dass die Genauigkeit der Rechnung leidet. In diesem Fall muss der Rechner "nur" 5000 Schritte ausführen.

Der Aufschlag am Boden könnte mit einem realistischeren Modell simuliert werden. Nur sind solche Modelle häufig komplex, rechenintensiv und nur unter grossem Aufwand zu [[Validierung|validieren]].


[[Kategorie:MMS]] [[Kategorie:Trans]] [[Kategorie:TransMod]] [[Kategorie:Modelle]]