Startseite
Forum
Fragen
Suchen
Formeleditor
Über Uns
Registrieren
Login
FAQ
Suchen
Foren-Übersicht
->
Mechanik
Antwort schreiben
Benutzername
(du bist
nicht
eingeloggt!)
Titel
Nachrichtentext
Smilies
Weitere Smilies ansehen
Schriftfarbe:
Standard
Dunkelrot
Rot
Orange
Braun
Gelb
Grün
Oliv
Cyan
Blau
Dunkelblau
Indigo
Violett
Weiß
Schwarz
Schriftgröße:
Schriftgröße
Winzig
Klein
Normal
Groß
Riesig
Tags schließen
Schreibt eure Formeln hier im Board am besten mit Latex!
So gehts:
Latex-Kurzbeschreibung
|
Formeleditor
Optionen
HTML ist
aus
BBCode
ist
an
Smilies sind
an
BBCode in diesem Beitrag deaktivieren
Smilies in diesem Beitrag deaktivieren
Spamschutz
Text aus Bild eingeben
Alle Zeiten sind GMT + 1 Stunde
Gehe zu:
Forum auswählen
Themenbereiche
----------------
Mechanik
Elektrik
Quantenphysik
Astronomie
Wärmelehre
Optik
Sonstiges
FAQ
Sonstiges
----------------
Off-Topic
Ankündigungen
Thema-Überblick
<table border="0" cellpadding="3" cellspacing="1" width="100%" class="forumline"> <tr> <th class="thCornerL" width="22%" height="26">Autor</th> <th class="thCornerR">Nachricht</th> </tr> <tr> <td width="22%" align="left" valign="top" class="row1"><span class="name"><a name=""></a><b>pleasediaaa</b></span></td> <td class="row1" height="28" valign="top"><table width="95%" border="0" cellspacing="0" cellpadding="0"> <tr> <td width="95%"><img src="templates/subSilver/images/icon_minipost.gif" width="12" height="9" alt="Beitrag" title="Beitrag" border="0" /><span class="postdetails">Verfasst am: 28. Feb 2016 19:44<span class="gen"> </span> Titel: Doppelpendel Differentialgleichung lösen</span></td> </tr> <tr> <td colspan="2"><hr /></td> </tr> <tr> <td colspan="2"><span class="postbody"><span style="font-weight: bold">Meine Frage:</span><br />Guten Tag :)<br /><br />ich bin Physikstudentin im ersten Semester und wir habe in unserem Computerkurs die Hausarbeit bekommen in Python unter anderem die Differentialglichung für ein Doppelpendel zu implementieren und wahlweise mit dem Eulerverfahren oder dem 4 Schritt Runge Kutter Verfahren zu lösen. Das Problem liegt momentan darin, dass wir in Mathe noch keine Differentialgleichungen hatten und uns deshalb 3 Tage lang was angelesen haben, in diesem Punkt jedoch verzweifeln.<br /><br /><span style="font-weight: bold">Meine Ideen:</span><br />Das sind unsere Ideen: Bisher konnten wir mithilfe einiger Beiträge die Bewegungsgleichungen für die beiden Massen des Doppelpendels herleiten und auch verstehen. An Hand eines einfachen Pendels konnten wir auch schon Teile unseres Programmcodes testen und wissen,dass diese funktionieren (Eulerstep, Flowmap). Wie der Name schon sagt soll Eulerstep die differentialgleichung nach dem Eulerverfahren lösen. Flowmap soll den Pointcareschnitt für die Bedingungen, die in der Funktion pointcare aufgelistet sind berechnen.<br />Die Funktion startcondition soll zufällige Anfangswerte für die Geschwindigkeiten generieren. Die gleichung in startcondition ist abhängig von der durch den User eingegeben Systemenergie und der Bedingung, dass x1 und x2 0 sind (also nicht ausgelenkt)<br />x entspricht unserem Winkel<br />y der ersten Ableitung, also der Winkelgeschwindigkeit<br />und z der Winkelbeschleunigung<br />die indizes 1 und 2 jeweils für die Masse 1 oder Masse 2.<br />Unser Problem ist jetzt, dass die Bewegungsgleichungen aus der Funktion Doppelpendel ja voneinander abhängen. Das heißt will ich z1 haben, brauche ich z2 und andersrum. So wie wir es bisher beurteilen können, liegt darin unser Problem. Also eher ein mathematisches. Nichts destotroz, schreibe ich euch unseren Programmcode mit hin.<br /> Im folgenden unser Programmcode:<br /><br />import numpy as np<br />from pylab import*<br />%matplotlib inline<br /><br />def doppelpendel(z1, y1, y2, t, dt, lam, mu, omega, e, x1, x2):<br /> z2 = -lam*omega**2*np.sin(x2)+lam*y1**2*np.sin(x1-x2)-lam*z1*np.cos(x1-x2)<br /> z1 = (omega**2)*np.sin(x1) - (mu/lam)*y2**2*np.sin(x1-x2) - (mu/lam)*z2*np.cos(x1-x2)<br /> return z1, z2<br />def doppelpendelz1(z2, y1, y2, t, dt, lam, mu, omega, e, x1, x2):<br /> z1 = (omega**2)*np.sin(x1) - (mu/lam)*y2**2*np.sin(x1-x2) - (mu/lam)*z2*np.cos(x1-x2)<br /> print("z1 ist {}".format(z1))<br /> return z1<br /><br />def doppelpendelz2(y1, y2, t, dt, lam, mu, omega, e, x1, x2):<br /> z2 = (-lam*omega**2*np.sin(x2) + lam*y1**2*np.sin(x1-x2) - omega**2*lam*np.cos(x1-x2)*np.sin(x1) + mu*y2**2*np.sin(x1-x2)*cos(x1-x2))/(1-mu*(np.cos(x1-x2))**2)<br /> print("z2 ist {}".format(z2))<br /> return z2<br /><br />def startcondition(lam, mu, e, steps):<br /> cachey1 = []<br /> cachey2 = []<br /> for i in range(1, steps+1):<br /> y1 = np.random.uniform(-100,100)<br /> y21 = -lam*y1 + sqrt(lam**2 * y1**2 - (lam**2 * y1**2 - 2*lam*e)/ mu)<br /> y22 = -lam*y1 - sqrt(lam**2 * y1**2 - (lam**2 * y1**2 - 2*lam*e)/ mu)<br /> cachey1.append(y1)<br /> cachey2.append(y21)<br /> cachey1.append(y1)<br /> cachey2.append(y22)<br /> return y1, y21<br /><br />def EulerStep(y1, y2, f1, f2, dt,t, lam, mu, omega, e, x1, x2): <br /> z2 = f1(y1, y2, t, dt, lam, mu, omega, e, x1, x2)<br /> z1 = f2(z2, y1, y2, t, dt, lam, mu, omega, e, x1, x2)<br /> <br /> y2 += dt*z2<br /> x2 += dt*y2<br /> <br /> y1 += dt*z1<br /> x1 += dt*y1<br /> t += dt<br /> return x1, x2, y1, y2, t, dt, z1, z2<br /><br />def FlowMap2(dt,t, steps, f1, f2, Integrator, stopcond, lam, mu, omega, e):<br /> cachey2=[]<br /> cachex2=[]<br /> cachet=[]<br /> cachez1=[]<br /> y1, y21 = startcondition(lam, mu, e, steps)<br /> j=1<br /> y2 = y21<br /> x1=0<br /> x2=0<br /> print( "lam = {} mu = {} omega = {} e = {} y1 = {} y2 = {}".format(lam,mu,omega,e,y1,y2))<br /> for i in range(1, steps+1):<br /> cond =False<br /> while cond==False:<br /> x1, x2, y1, y2, t, dt,z1,z2 = Integrator(y1, y2, f1, f2, dt,t, lam, mu, omega, e, x1, x2)<br /> cond = stopcond(x1,y1)<br /> cachez1.append(z1)<br /> cachet.append(t)<br /> cachex2.append(x2)<br /> cachey2.append(y2)<br /> j +=1<br /> teinschwing = cachet[0]<br /> return plot(cachex2, cachey2, '.')<br /><br />def poincare2(x1, y1):<br /> if fabs(np.sin(x1)) <= 0.01 and y1 > 0:<br /> return True<br /> else:<br /> return <br /><br />Beim ausführen des Codes, haben wir für f1, doppelpendelz2 und für f2 doppelpendelz1 eingesetzt.<br />Es stekkt sich herraus, dass die Werte in die Höhe schießen und irgendwan so groß werden, und zu nan oder inf werden.<br /><br />Wir sind euch sehr dankbar für Anregungen oder Erklärungen wie diese gekoppelten Bewegungsgleichungen miteinander zu handeln sind.<br />Vielen Dank im Vorraus!</span></td> </tr> </table></td> </tr> <tr> <td colspan="2" height="1" class="spaceRow"><img src="templates/subSilver/images/spacer.gif" alt="" width="1" height="1" /></td> </tr> </table>
Registrieren
Login
FAQ
Suchen
