Differences

This shows you the differences between two versions of the page.

--- a_mechanik:drehschwingungen:gruppenseiten:gruppe316:start [28 January 2021 18:40] – 1. annapaul
+++ a_mechanik:drehschwingungen:gruppenseiten:gruppe316:start [29 January 2021 20:37] (current) – annapaul
@@ Line 6: / Line 6: @@
 Diese Wiki-Seite behandelt den Home-Lab Versuch " Drehschwingung".
 ====== Versuchsdurchführung ======
-===== Messung 1 =====
+===== Allgemeiner Aufbau =====
-Bei der 1. Messung wurde eine
+Zunächst wurde eine Gitarrensaite fest im das gewünschte Objekt gewickelt und anschließend noch mit einem Faden fest darüber gewickelt damit dies hält. Diese Konstruktion wurde mit Hilfe von Klebestreifen an einem Tisch befestigt. Dies wurde auf dem Tisch zusätzlich noch mit Büchern beschwert, da ein Klebestreifen als nicht genügend Sicherung angesehen wurde.
-===== Messung 2 =====
+Der Aufbau ist in folgendem Bilde ersichtlich
+{{ :a_mechanik:drehschwingungen:gruppenseiten:gruppe316:stricknadelbild.jpg?400 |}}
-=====Messung 3 Luftwiderstand =====
+Das gewünschte Objekt wird um einen Winkel ausgedreht und die Zeit die dieses für eine Schwingung braucht wird mit Hilfe der Augen und einer Stoppuhr bestimmt.
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:luftwid.jpg?400 |}}
+===== Messung 1 =====
+Bei der 1. Messung wurde eine 16g schwere Sticknadel, mit einem Radius von 4,5mm, an einer Gitarrensaite befestigt und bei unterschiedlichen Längen mit jeweils einem ungefährenen Winkel von 90° aus der Ruhelage verschoben. Dort wird diese Losgelassen und die Zeit von drei schwingungen wurde mit Hilfe einer Stopuhr gemessen. Dies wurde bei jeder Länge 5 Mal wiederholt.
+{{ :a_mechanik:drehschwingungen:gruppenseiten:gruppe316:stricknadelgrafik.jpg?400 |}}
-Die Messung 3 wurde Identisch zur Messung 2 durchgeführt, jedoch wurde um den Einfluss des Luftwiderstandes zu Untersuchen, an den Besen ein Starkes Blatt Papier geklebt (130g/m^2) (siehe Foto).  Der Besen wurde dann auch 5 mal bei gleichbleibenden Winkel fallen gelassen und mit dem gleichen Winkel ohne das Papier verglichen.
+===== Messung 2 =====
+Für die 2. Messung wird Messung 1 wiederholt, jedoch nun mit einem Band anstelle einer Gitarrensaite.
+{{ :a_mechanik:drehschwingungen:gruppenseiten:gruppe316:band.jpg?400 |}}
-===== Zeitmessung =====
+=====Messung 3 und 4 =====
-Um den Fehler bei der Zeitmessung zu minimieren, wurde das Experiment per Video aufgezeichnet. Somit muss keine Schrecksekunde berücksichtigt werden und es kann eine genaue Angabe über die Fallzeit gemacht werden. Die Videos wurden in einem Schnittprogramm hochgeladen
+Hier wurden verschiedene Objekte bei einer Länge mehrfach Schwingen lassen.
+Das 1. Objekt war ein Topfdeckel.
+{{ :a_mechanik:drehschwingungen:gruppenseiten:gruppe316:topfdeckelbild.jpg?400 |}}
-Zur Zeitmessung wurde ein Video aufgenommen, welches im Nachhinein analysiert wurde. Hierbei wurde Video An in 25 fps und Video Ay in 30 fps geschossen. Bei den Videos wurden die genauen Frames bestimmt, bei denen der Besen losgelassen wurde und bei welchem er auf dem Boden aufschlug. Durch die Differenz wurde die Zeit bestimmt, die der Besen zum Fallen brauchte.
+Das 2. Objekt ein Teddybär.
-Für die Unsicherheit wurden jeweils 3fps angenommen. Ein Frame zu Anfang beim Loslassen, eins beim Aufschlagen und eins fürs menschliches Versagen. Dies wirkt sich bei den Unterschiedlichen Videos natürlich anders auf die Sekunden aus. Bei Video An entspricht dies +- 0,12 sec und bei dem Video Ay +- 0,1 sec.
+{{ :a_mechanik:drehschwingungen:gruppenseiten:gruppe316:teddy.jpg?400 |}}
-Im Anschluss folgt ein screenshot um besser sehen zu können wie genau die Zeit gemessen wurde.
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:zeitmessung.png?600 |}}
-Ein Ausschnitt der Excel Tabelle von Messung 1 Ay. Hierbei ist zu bemerken, dass bei der Zeit (aa:bb:cc:dd) aa für die Stunden, bb für die Minuten, cc für die Sekunden und dd für die fps steht.
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:bspexceltabay.png?400 |}}
 ====== Messdaten ======
-In der ersten Tabelle befinden sich die Messwerte zur ersten Messung mit einem Besen der Länge l=1,2m. Für die Fallzeiten wurde bereits der Mittelwert aller 5 Messungen bestimmt, der auch später in der Auswertung verwendet wird. Die zweite Tabelle enthält die Daten zur Messung mit einem Besen der Länge l=1,34m.
-^ Winkel in Grad  ^ Mittelwert der Fallzeit in Sekunden ^ Simulierte Fallzeit^
-| 5              | 1,21 | 1,03|
-| 7              | 1,06 | 0,95|
-| 10             | 0,98 | 0,85|
-| 12             | 0,88 | 0,79|
-| 17             | 0,75 | 0,69|
-| 23             | 0,66 | 0,6|
-| 27             | 0,63 | 0,56|
-| 32             | 0,55 | 0,51|
-| 37             | 0,50 | 0,46|
-| 42             | 0,47 | 0,43|
-| 50             | 0,40 | 0,37|
-| 54             | 0,38 | 0,35|
-| 59             | 0,35 | 0,31|
-^ Winkel in Grad  ^ Mittelwert der Fallzeit in Sekunden ^ Simulierte Fallzeit^
-| 9,9              | 1,08 | 0,9|
-| 17,2             | 0,84 | 0,73|
-| 22,3             | 0,76 | 0,65|
-| 25,4             | 0,73 | 0,61|
-| 33,3             | 0,59 | 0,53|
-| 35,6             | 0,56 | 0,5|
-| 42,4             | 0,57 | 0,45|
-| 46,6             | 0,49 | 0,42|
-| 48,4             | 0,48 | 0,41|
-====== Computerprogramm ======
-===== Zeitschrittverfahren =====
-Das Computerprogramm zum numerischen Lösen der Differentialgleichung mit dem Zeitschrittverfahren wurde in einem Mathematica Notebook geschrieben. Hierbei wurden die Besenlänge l=1.45m und der Anfangswinkel phi0= 0.25 rad gewählt. Zunächst wird in "bewegung" das Zeitschrittverfahren angewandt, um den Winkel in Abhängigkeit von der Zeit zu erhalten und somit den Bewegungsverlauf der Besenspitze. Der erste Wert dieser Liste gibt uns somit den Startwinkel und der letzte Wert die Fallzeit. Dies wird eingebunden in "weitere", um das erste Tupel aus Anfangswinkel und Fallzeit zu erhalten. Weiterhin wird eine Liste erstellt, die die Fallzeit für größer werdende Winkel angibt und schließlich kann das ganze mit dem in Mathematica eingebauten Befehl ListPlot gezeichnet werden. Zum genaueren Untersuchen der Werte ist es jedoch von Vorteil, sich die Listen ausgeben zu lassen.
-<code>
-l = 1.45;
-\[CurlyPhi]0 = 0.25;
-g = 9.81;
-zeitschritt = 0.005;
-\[Tau] = Sqrt[(2 l)/(3 g)];
-bgl[\[CurlyPhi]_] := Sin[\[CurlyPhi]]/\[Tau]^2;
-bewegung := Module[{liste, t, \[CurlyPhi], \[CurlyPhi]1},
-   t = t0; t0 = 0; \[CurlyPhi] = \[CurlyPhi]0;
-   t = t0; \[CurlyPhi]1 = 0;
-   liste = {{\[CurlyPhi]0, t0}};
-   While[\[CurlyPhi] < \[Pi]/
-, \[CurlyPhi]1 = \[CurlyPhi]1 + zeitschritt*bgl[\[CurlyPhi]];
-    \[CurlyPhi] = \[CurlyPhi] + zeitschritt*\[CurlyPhi]1;
-    t = t + zeitschritt;
-    liste = Append[liste, {\[CurlyPhi], t}]]; liste];
-weitere =
-  Module[{\[CurlyPhi]0, lüst},
-   lüst = {{First[First[bewegung]], Last[Last[bewegung]]}};
-   \[CurlyPhi]0 = 0.25;
-   While[\[CurlyPhi]0 < \[Pi]/2,
-    lüst = Append[
-      lüst,
-      {\[CurlyPhi]0, Last[Module[
-         {liste, t, \[CurlyPhi], \[CurlyPhi]1},
-         liste = {t0};
-         \[CurlyPhi] = \[CurlyPhi]0;
-         t = t0;
-         \[CurlyPhi]1 = 0;
-         While[
-          \[CurlyPhi] < \[Pi]/2,
-          \[CurlyPhi]1 = \[CurlyPhi]1 + zeitschritt*bgl[\[CurlyPhi]];
-          \[CurlyPhi] = \[CurlyPhi] + zeitschritt*\[CurlyPhi]1;
-          t = t + zeitschritt;
-          liste = Append[liste, t] ];
-         liste]]}]; \[CurlyPhi]0 = \[CurlyPhi]0 + 0.01]; lüst];
-ListPlot[weitere,
- AxesLabel -> {"Anfangswinkel \!\(\*SubscriptBox[\(\[CurlyPhi]\), \(0\
-\)]\) in rad", "Fallzeit t in s"}, ImageSize -> Large,
- GridLines -> {{.2, .4, .6, .8, 1.0, 1.2, 1.4,
-.6}, {.1, .2, .3, .4, .5, .6, .7, .8}}, PlotLabel -> "Fallzeit"]
-</code>
-Die entstehende Grafik sieht dann wie folgt aus:
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:fallzeit_numerisch.png?600 |}}
-Um schließlich Aussagen treffen zu können über den Einfluss der Zeitschritte auf die simulierte Fallzeit, kann man sich entweder die gesamte Liste Ausgeben lassen, z.B bei einem Zeitschritt von 0.03s:
-<code>
-In: bewegung
-Out: {{0.25, 0}, {0.25226, 0.03}, {0.256799, 0.06}, {0.263658,
-.09}, {0.272897, 0.12}, {0.284598, 0.15}, {0.298864,
-.18}, {0.315819, 0.21}, {0.33561, 0.24}, {0.358409,
-.27}, {0.384413, 0.3}, {0.413841, 0.33}, {0.446943,
-.36}, {0.483992, 0.39}, {0.525291, 0.42}, {0.57117,
-.45}, {0.621986, 0.48}, {0.678125, 0.51}, {0.739993,
-.54}, {0.808019, 0.57}, {0.882649, 0.6}, {0.964333,
-.63}, {1.05352, 0.66}, {1.15065, 0.69}, {1.25612, 0.72}, {1.37027,
-.75}, {1.49337, 0.78}, {1.62558, 0.81}}
-</code>
-Da man an dieser Stelle aber die Fallzeit erhalten möchte, reicht es sich den letzten Wert anzusehen:
-<code>
-In: Last[bewegung]
-Out: {1.62558, 0.81}
-</code>
-Hierbei fällt jedoch auf, dass für größere Zeitschritte der letzte Wert meist einen Anfangswinkel enthält, der größer als 90° ist. Für Winkel größer als 90° sollte das Zeitschrittverfahren jedoch bereits enden. Je kleiner also die Zeitschritte sind, desto präziser wird die Fallzeit bestimmt.
-=====Einfluss Besenlänge und Anfangswinkel=====
-Eine weitere Visualisierung der Bewegungs außerhalb des Zeitschrittverfahrens ist mit folgendem separaten Code ersichtlich:
-<code>
-g := 9.81;
-\[Tau][l_] := Sqrt[(2 l)/(3 g)];
-sol[\[CurlyPhi]0_, l_] :=
- NDSolve[{\[CurlyPhi]''[t] ==
-    Sin[\[CurlyPhi][t]]/\[Tau][l]^2, \[CurlyPhi][
-] == \[CurlyPhi]0, \[CurlyPhi]'[0] == 0}, {\[CurlyPhi][
-    t], \[CurlyPhi]'[t], \[CurlyPhi]''[t], \[CurlyPhi]'''[t] }, {t, 0,
-}]
-Plot[Evaluate[{\[CurlyPhi][t], \[CurlyPhi]'[t], \[CurlyPhi]''[
-     t], \[CurlyPhi]'''[t]} /. sol[0.25, 1.45]], {t, 0, 1},
- AxesLabel -> {"x", "y"}, GridLines -> Automatic,
- PlotLabel ->
-  "Numerische Lösung der Bewegungsgleichung und ihre Ableitungen",
- PlotLegends -> Automatic, ImageSize -> Large]
-</code>
-Der sich ergebende Plot ist gut geeignet, um qualitative Aussagen über die Fallzeit zu treffen, denn die Nullstelle der dritten Ableitung gibt diese an. Dies hängt natürlich damit zusammen, dass die Beschleunigung beim Aufprall auf den Boden am größten ist und die Funktion der Beschleunigung, also die zweite Ableitung, dort ihr Maximum hat. Daran lässt sich der Einfluss der Stablänge zeigen.
-Der folgende Plot gehört zum obigen Code mit der Besenlänge l=1.45m:
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:screenshot_155_.png?600 |}}
-Die Legende gilt auch für die weiteren Plots. Die nachfolgenden beschreiben die Bewegung eines Besens der Länge l=1.2m, 1.0m, 0.8m:
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:l_1.2m.png?600 |}}
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:l_1m.png?600 |}}
+Messung 1
+^ Länge in cm  ^ Mittelwert der Periodendauer in Sekunden ^ Standardabweichung der Periodendauer^
+| 55          | 1,19 |0,027|
+| 43,5             | 1,07 | 0,032|
+| 37          |0,97 | 0,041|
+| 29             | 0,86 | 0,026|
+| 25             | 0,79 | 0,023|
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:l_0.8m.png?600 |}}
-Man erkennt, dass die Nullstelle des orangenen Graphen immer weiter nach links wandert, je kürzer der Besen ist. Je länger also ein Besenstiel ist, desto länger ist auch die Fallzeit.
-Mit dem gleichen Verfahren lässt sich auch zeigen, dass die Fallzeit für kleinere Winkel größer ist.
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:phi_0.5.png?600 |}}
+Messung 2
+^ Länge in cm  ^ Mittelwert der Periodendauer in Sekunden ^ Standardabweichung der Periodendauer^
+| 28            | 11,64 | 0,760|
+| 23            | 10,38 | 0,205|
+| 18            | 9,25 | 0,109|
+| 13,5            | 7,99 | 0,152|
-In einer reibungsfreien Umgebung hat die Masse keinen Einfluss auf die Fallzeit, sie wird aus der Bewegungsgleichung rausgekürzt.
-=====Beschleunigung der Stabspitze=====
+Messung 3 und 4
-Es wird weiterhin im selben Notebook gearbeitet wie im vorherigen Abschnitt. Hier soll gezeigt werden, dass die Beschleunigung der Stabspitze eines beliebig langen Besens zu einem beliebigen Anfangswinkel mit einer größeren Beschleunigung zu Boden fällt als eine frei fallende Punktmasse. Die Herleitung der Formel erfolgt in unserem Versuchsbericht im Unterkapitel "Beschleunigung".
+^ Objekt ^ Länge in cm ^ Mittelwert der Periodendauer in Sekunden ^ Standardabweichung der Periodendauer^
-Mit dem folgenden Code wurde eine Manipulate Umgebung erstellt, in der die Besenlänge und der Anfangswinkel variiert werden können. Der Plot stellt die Tangentialbeschleunigung dar. Für unsere Betrachtung ist nur das erste Maximum von Bedeutung, alle weiteren, die eventuell auftreten, können vernachlässigt werden. Außerdem wurde eine Linie eingezeichnet, die den Wert der Fallbeschleunigung markiert.
+| Topfdeckel | 43             | 19,28 | 0,081|
-<code>
+| Teddy | 42            | 4,32 | 0,126|
-Manipulate[
- Plot[{Evaluate[
-/2*g*Sin[\[CurlyPhi][t]] /. sol[\[Alpha], \[Lambda]]], 9.81}, {t,
-, 1.5},
-  AxesLabel -> {"t", "\!\(\*SubscriptBox[\(\[Alpha]\), \(T\)]\)"},
-  GridLines -> Automatic,
-  PlotLabel -> {"Tangentialbeschleunigung Besenspitze für l=", \
-\[Lambda],
-    " und \!\(\*SubscriptBox[\(\[CurlyPhi]\), \(0\)]\)=", \[Alpha],
-    ""}, PlotLegends ->
-   "\!\(\*SubscriptBox[\(\[Alpha]\), \
-\(T\)]\)=\!\(\*FractionBox[\(3\), \(2\)]\)g*Sin[\[CurlyPhi][t]]"], {\
-\[Alpha], 0.01, \[Pi]/2}, {\[Lambda], 0.1, 2.0}]
-</code>
-Nun wird für einige beliebige Werte exemplarisch gezeigt, dass die Stabspitze mit einer höheren Beschleunigung fällt als eine Punktmasse es tun würde.
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:0.25_1.45.png?600 |}}
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:0.4_0.25.png?600 |}}
+======Fehler ======
+{{ :a_mechanik:drehschwingungen:gruppenseiten:gruppe316:gfehler.jpg?400 |}}
+Hierbei:
+r=0,0003m; u(r)= 0,00025m
+l=0,20m; u(l)=0,005m
+R=0,045m; u(R) als 0 angenommen, Angabe des Herstellers
+m=0,016kg; u(m)=0,0005kg
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:0.7_0.25.png?600 |}}
+Daraus folgt, dass der Fehler für G 23,9*10^9 Pa ist
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:1.45_0.5.png?600 |}}
-{{ :a_mechanik:kippender_besenstiel:gruppenseiten:gruppe316:1.45_1.06.png?600 |}}
+{{ :a_mechanik:drehschwingungen:gruppenseiten:gruppe316:ifehler.jpg?400 |}}
+Für den Topf gilt folgendes:
+r=0,0003m; u(r)= 0,00025m;
+l=0,43m; u(l)=0,005m;
+T=9,28s; u(T)=0,036s;
+Daraus folgt, dass der Fehler für den Topfdeckel 0,112kg*(m^2) ist
+Für den Teddy gilt folgendes:
+r=0,0003m; u(r)= 0,00025m;
+l=0,42m; u(l)=0,005m;
+T=4,32s; u(T)=0,056s;
+Daraus folgt, dass der Fehler für den Teddy 0,003kg*(m^2) ist

Tools

menus and quick search

quick search

site status

Page Tools

meta data for this page

Differences