meta data for this page
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
a_mechanik:kippender_besenstiel:gruppenseiten:gruppe342:start [ 4 January 2021 18:24] – [Vorüberlegungen] lukaskoepp | a_mechanik:kippender_besenstiel:gruppenseiten:gruppe342:start [21 January 2021 12:22] (current) – [Ergebnisse] lukaskoepp | ||
---|---|---|---|
Line 10: | Line 10: | ||
Wir erhalten folgende Bewegungsgleichung für den Winkel $\phi$ des Körpers in Bezug auf das Drehmoment $J$ | Wir erhalten folgende Bewegungsgleichung für den Winkel $\phi$ des Körpers in Bezug auf das Drehmoment $J$ | ||
- | $ J \ddot{\phi} = F_G $ | + | $ J \ddot{\phi} = F_G \cdot \frac{l}{2} = m \cdot g \cdot sin(\phi) \cdot \frac{l}{2}$ |
+ | |||
+ | Für J ohne Luftreibung gilt: | ||
+ | |||
+ | $ J = \int_{0}^{l} x^2 \cdot \frac{m}{l} \,dx = \frac{m \cdot l^2}{3} $ | ||
+ | |||
+ | Es ergibt sich nun die Differentialgleichung $ \ddot{\phi} = \frac{3 \cdot g}{2 \cdot l} \cdot sin(\phi) $ und wir erkennen, dass die Masse bei einer reibungsfreien Bewegung irrelevant ist. | ||
+ | |||
+ | Die Länge $l$ tritt antiproportional in der Gleichung auf, ist also die Länge größer so wird auch die Fallzeit größer | ||
+ | |||
+ | Gleiches gilt hier für das Ausbalancieren eines Besens. Ist die Masse nicht homogen verteilt, so drehen wir den Besen, dass der massereiche Anteil weit entfernt ist. Dies erhöht das Trägheitsmoment und der Besen kommt nicht so schnell aus der Balance. | ||
+ | |||
+ | Der Anfangswinkel $\phi_0$ spielt eine entscheidene Rolle. Für $\phi_0 = 0$ befinden wir uns auf einem instabilen Extremum, was theoretisch eine unendlich lange Fallzeit haben kann. Real wird der Besenstiel irgendwann aus der Ruhe gebracht. Durch Erhöhung vom Anfangswinkel erreichen wir natürlcih kürzere Fallzeiten, bis wir bei $\phi_0 = \pi/2$ beim Boden ankommen. | ||
====== Messung ====== | ====== Messung ====== | ||
Line 16: | Line 28: | ||
- | Zuerst betrachten wir einen Besen. Dieser hat die Länge $l_{Besen} = 1,30 \pm 0,01m$. | + | Zuerst betrachten wir einen Besen. Dieser hat die Länge $l_{Besen} = 1,30 \pm 0,03m$. |
- | Wir messen die Starthöhe und bestimmen $\phi_0 = \arctan{s/h}$, wobei $s$ die Entferneung entlang der Küchenzeile und $h$ bei die Höhe der Kückenzeile ist. | + | Wir messen die Starthöhe und bestimmen $\phi_0 = \arctan{s/h_k}$, wobei $s$ die Entferneung entlang der Küchenzeile und $h_k$ bei die Höhe der Kückenzeile ist. |
- | Ab einer gewissen Entfernung reicht der Besen nicht mehr bis zur Küchenplatte. Nun messen wir die Höhe $h$ bei einem festen Abstand $s_{fest}$ und bestimmen den Winkel erneut als $\phi_0 = \arctan{s/h}$. | + | Ab einer gewissen Entfernung reicht der Besen nicht mehr bis zur Küchenplatte. Nun messen wir die Höhe $h$ bei einem festen Abstand $s_{Mark}$ und bestimmen den Winkel erneut als $\phi_0 = \arctan{s_{Mark}/h}$. |
- | Unsicherheiten $u(h) = 0,03m; u(T) = 0.1s???$ | + | {{: |
- | Wir lassen den Stab mit zwei Fingern los, sodass er keinen Schwung bekommt. Die Zeit wird mit der akkustischen Stopuhr gemessen. Dabei klopfen wir beim Loslassen auf die Küchenzeile. | + | Unsicherheiten: |
- | Höhe Küchenzeile | + | $ u(s$ bzw. $h) = 0,03m, u(T) = 0.1s $ |
- | Entfernung Markierung $s_{Mark} = 0,90 \pm 0,01m$. | + | $ u(\phi_0) = \sqrt{(\frac{h^2 \cdot u(s)}{h \cdot (h^2 + s^2)})^2 + (\frac{- s\cdot h^2 \cdot u(h)}{h^2 \cdot (h^2 + s^2)})^2} = \frac{0.03m}{\sqrt{h^2 + s^2}} \leq 2^° $ |
+ | |||
+ | Den Fehler von einem Satz von Messpunkten $T_i$ erhalten wir aus dem Standardfehler der als Standabweichung$/ | ||
+ | |||
+ | Wir lassen den Stab mit zwei Fingern los, sodass er keinen Schwung bekommt. Die Zeit wird mit der akkustischen Stoppuhr gemessen. Dabei klopfen wir beim Loslassen auf die Küchenzeile. | ||
+ | |||
+ | Höhe Küchenzeile $h_k = 0,92 \pm 0,03m$. | ||
+ | |||
+ | Entfernung Markierung $s_{Mark} = 0,90 \pm 0,03m$. | ||
^ ^ Fallzeit/s Messung Nr. ||||| | ^ ^ Fallzeit/s Messung Nr. ||||| | ||
Line 60: | Line 80: | ||
| 0,320 | 0,227 | 0,220 | 0,235 | 0,203 | 0,237 | | | 0,320 | 0,227 | 0,220 | 0,235 | 0,203 | 0,237 | | ||
+ | ===== Messung: Einfluss der Luftreibung ===== | ||
- | ====== Computerprogramm ====== | + | Wir hängen ein Blatt Papier mit Fläche $A=10,5\cdot 18,6\ (cm)^2 |
- | ===== Messung ===== | + | Wir messen die Fallzeiten |
- | <code c [enable_line_numbers=" | + | | Fallzeit/ |
- | import numpy as np | + | und bestimmen die Unsicherheit $u(T)$ erneut über die Standardabweichung. |
- | import matplotlib.pyplot as plt | + | |
- | h_kueche | + | ===== Ergebnisse ===== |
- | s_mark | + | |
- | s1 = [0.160, 0.365, 0.665, 0.930] | ||
- | h1 = [0.655, 0.320] | ||
- | winkel = [np.arctan(x/ | + | {{: |
- | print(' | + | |
- | t_laengen1 | + | Der Übersichtlichkeit halber zeichnen wir die numerischen Graphen kontinuierlich, |
- | | + | |
- | | + | |
- | | + | |
- | t_laengen1 = [sum(x)/ | + | |
- | t_hoehen1 = [[0.315, 0.299, 0.315, 0.304, 0.292], | + | Wir sehen, dass unsere Vorhersage nicht ganz stimmen. Die Diskussion möglicher Fehlerquellen erfolgt weiter unten. |
- | [0.255, 0.197, 0.268, 0.241, 0.262]] | + | |
- | t_hoehen1 = [sum(x)/ | + | |
- | t1 = t_laengen1 + t_hoehen1 | + | |
- | print(' | + | Wir sehen außerdem, dass die Luftreibung die Fallzeit des Besenstiels durch das Anhängen eines Blatt Papiers verlängert. Aufgrund einer Messung zu einem späteren Zeitpunkt können wir hier aber nicht ganz sicher sein, ob es nicht kleine Veränderungen z.B. in der Zeitstoppung gab. |
- | | + | ====== Weitere Messunsicherheiten ====== |
- | # Winkel wie gehabt | + | |
- | t_laengen2 = [[0.723, 0.667, 0.708, 0.681, 0.637], | + | Wir betrachten weitere Störquellen. |
- | | + | |
- | | + | |
- | | + | |
- | t_laengen2 = [sum(x)/ | + | |
- | t_hoehen2 = [[0.286, 0.307, 0.333, 0.283, 0.303], | ||
- | [0.227, 0.220, 0.235, 0.203, 0.237]] | ||
- | t_hoehen2 = [sum(x)/ | ||
- | t2 = t_laengen2 + t_hoehen2 | ||
- | x_err, y_err = 0, 0 | + | In Bezug auf die Schallreflexion mit $c \approx 300 m/s$ haben wir bei unserer Küche mit ca. 3m Länge bereits bei Reflexion von zum Beispiel $R \approx 80 \% $ schon nach $t=0.1s$ eine Dämpfung auf $0,9^{10} \approx 10 \% $ der Ausgangsleistung. Außerdem vermuten wir, dass die Erkennung des Anstiegs des Schallpegels über einen Treshhold nicht von der Relflexion beeinflusst wird, bzw. im Aufschlag des Besen der ersten Schallpeak |
- | plt.errorbar(winkel, t1, xerr=x_err, yerr=y_err, marker=' | + | |
- | plt.errorbar(winkel, t2, xerr=x_err, yerr=y_err, ls='', | + | |
- | plt.xlabel(' | + | |
- | plt.ylabel(' | + | |
- | plt.legend() | + | |
- | plt.show() | + | |
- | </ | + | Die Reflexivität der Wand ist vermutlich deutlich kleiner. |
+ | Beim Erzeugen des ersten Audio-Peaks durch Klopfen auf den Boden und gleichzeitigen Loslassens kann jedoch leicht ein Zeitfehler entstehen, da man dazu neigt, das Klopfen mit dem Aufschlag zu synchronisieren. Dies könnte verantwortlich sein für die systematisch zu geringe gemessene Fallzeit. | ||
+ | ====== Luftreibung ====== | ||
- | ====== Balder lutscht Eier ====== | + | Wir betrachten als Modell für die Luftreibung eine Form der Newton Reibung für turbulente Fluide. Wir erhalten als Reibungsmodell also |
- | ====== Weitere Messunsicherheiten ====== | + | $ F_R = ½ A c_w \rho_{Luft} v^2$ |
- | Wir betrachten weitere Störquellen. | + | wobei $A$ die Querschnittsfläche, |
+ | $ v = S/T = \pi/2 \cdot l/2 \cdot 1/T $ | ||
- | In Bezug auf die Schallreflexion | + | Wir erhalten für den Besen mit $r_B \approx |
- | Die Reflexivität der Wand ist vermutlich deutlich kleiner. | + | $ F_{R, Besen} \approx ½ r l (\pi/2 \cdot l/2 \cdot 1/T)^2 \approx 0,030 N \\ F_{R, Eisen} \approx 0,005 N $ |
- | Beim Erzeugen des ersten Audio-Peaks | + | Verglichen mit der Gravitation sind wir also um mindestens zwei Größenordnungen kleiner, selbst wenn die effektive Gravitation $g_{eff} \leq 9.81$ etwas kleiner als der Ortsfaktor durch durch den $\sin{\phi}$ ist. Wir betrachten den Einfluss der Luftreibung zumindest für die im Versuch verwendeten Geometrien als irrelevant. |
- | Außerdem ist im Moment des Loslassens ein zusätzliches Anschubsen nicht auszuschließen. | + | Wir betrachten noch einmal das Anhängen eines Blatt Papiers an den Besen. Für eine lange Rechteckplatte erhalten wir hier $c_w = 2$. Wir haben $A_{Papier} = 0,02\ m^2$ gegeben. Damit errechnet sich mit |
- | ====== Diese Seiten ====== | + | $F_{R, Papier} \approx |
- | Diese Seite und ihre Unterseiten sind Ihr Bereich im APwiki für die Bearbeitung | + | |
- | des Heim-Versuchs " | + | |
- | im Präsenzpraktikum das Heft hat. Das heißt, es ist Ihre Logbuch für das, was | + | |
- | Sie konkret experimentell und bei der Programmierung durchführen. | + | |
- | Legen Sie Fotos ab, notieren Sie Messwerte, laden sie ihr Programm hoch. Form | + | für den Besenstiel ein Gesamtwert von $F \approx F_{R, Besen} + F_{R, Papier} \approx 0,11 N$, was theoretisch einen kleinen Einfluss auf die Fallzeit haben könnte. |
- | und Formatierung sind dabei zweitrangig. | + | ====== Computerprogramm ====== |
- | Damit dieser Bereich diese Aufgabe erfüllen kann, haben wir ihn mit speziellen | + | ===== Numerisches Verfahren ===== |
- | Zugriffsrechten ausgestattet: | + | ==== 1 Bestimmung der Fallzeit ==== |
- | - Ihre Gruppe hat das exklusive Schreibrecht für diese Seite. | + | |
- | - Die Seite ist nur für Ihre Gruppe, die Tutoren und die Praktikumsleitung einsehbar. | + | |
- | Unten auf dieser Seite finden Sie einen Abschnitt " | ||
- | findet die Kommunikation mit Ihrem Tutor statt. Sie oder er wird Ihnen dort | ||
- | Rückmeldung zu Ihrem Versuchsbericht geben. | ||
- | Hier im Wiki gibt es [[: | + | <code python |
- | Formatierung ihres Versuchsberichts mit Latex]]. Den Versuchsbericht geben Sie | + | import numpy |
- | dann im Ilias ab. | + | import matplotlib.pyplot as plt |
- | < | ||
- | den Start erleichtern. Sie können es nach Belieben löschen und durch Ihre | ||
- | eigenen inhalte ersetzen. </ | ||
- | ===== Computerprogramm ===== | + | def evaluate(length, |
- | Dokumentieren Sie hier im Wiki das Programm, das Sie für die Lösung der Bewegungsgleichung des Besenstiels geschrieben haben. Dafür eignet sich dafür besonders gut die Umgebung < | + | tau_quadrat |
+ | phi = phi0 | ||
+ | v_phi = 0 | ||
+ | time = 0 | ||
+ | while phi < numpy.pi | ||
+ | a_phi = numpy.sin(phi) / tau_quadrat | ||
+ | v_phi += timestep * a_phi | ||
+ | phi += timestep * v_phi | ||
+ | time += timestep | ||
+ | return time | ||
- | Außerdem ist es möglich einen Link zum Download des präsentierten Programm-Codes anzuzeigen. Dazu geben Sie in dem einleitenden code-Tag einen Dateinamen an. Der Download bezieht sich unmittelbar auf das Im Editor eingetragene Programmstück. Ein getrennter Upload ist nicht nötig. | ||
- | Beispiel: | + | # testing the evaluation |
- | < | + | prec_test |
- | #include < | + | phi0s_test = [0 + prec_test * i for i in range(1, round(2 * numpy.pi / prec_test) + 1)] |
- | int main() | + | sols_test = [evaluate(1.45, phi0, 0.01) for phi0 in phi0s_test] |
- | { | + | plt.plot(phi0s_test, |
- | | + | plt.legend() |
- | return 0; | + | plt.show() |
- | } | + | print(' |
- | </ | + | |
- | wird dargestellt als | + | |
- | <code c [enable_line_numbers=" | + | |
- | #include <stdio.h> | + | |
- | int main() | + | |
- | { | + | |
- | | + | |
- | return | + | |
- | } | + | |
- | </ | + | |
- | ===== Bilder einbinden ===== | ||
- | Ihr Versuchsaufbau sollte so beschrieben sein, dass er für sich stehend verständlich ist - gerne mit einem Foto. | ||
- | Ein Bild laden Sie ins Wiki, indem Sie im Editor | + | def find_sol(length, phi0, prec): |
+ | prev = None | ||
+ | for i in range(10): | ||
+ | sol = evaluate(length, phi0, 0.1 ** i) | ||
+ | if prev and abs(sol - prev) < 0.5 * 0.1 ** prec: | ||
+ | break | ||
+ | prev = sol | ||
+ | return sol | ||
- | Im einfachsten Fall landet ein Bild direkt an der Stelle im Text, an der Sie es eingefügt haben (Siehe [[doku> | ||
- | ===== Tabellen | + | def find_sols(lengths, |
- | Für eine Tabelle mit Ihren Messwerten gibt es im oben im Editfenster des Wikis eine Hilfsfunktion. Sie versteckt sich hinter einem Knopf der so aussieht, wie ein hellblauer Taschenrechner. | + | sols = {} |
+ | for length in lengths: | ||
+ | phi_sols | ||
+ | for phi0 in phi0s: | ||
+ | phi_sols += [find_sol(length, | ||
+ | sols[length] | ||
+ | return sols | ||
+ | </ | ||
+ | |||
+ | ==== 2 Beschleunigung des Endes ==== | ||
- | ===== Syntax und Funktionen im Wiki ===== | + | Die Winkelbeschleunigung beträgt $ \ddot{\phi} |
- | Hier noch Links zu | + | |
- | * den [[doku> | + | |
- | * [[: | + | |
- | * [[: | + | |