SDS Facebook-Gruppe
Moderator: (M) Mod.-Team Allgemein
Ich hab mal die Lagrange-Gleichungen nach s aufgestellt (wenn man nichts vorgibt, ist s ja auch ein Freiheitsgrad!) Ergebnis: DEUTLICH kleinere Terme. Das ganze funktioniert wohl besser, wenn man s verwendet.
Könnte mir jemand mal verraten, wie man die symbolischen Variablen und die s-Funktion überhaupt aufstellt? Ich schreibe meine BA außerhalb und kann daher nicht zur Übung
Könnte mir jemand mal verraten, wie man die symbolischen Variablen und die s-Funktion überhaupt aufstellt? Ich schreibe meine BA außerhalb und kann daher nicht zur Übung
s ist aber keine allgemeine Koordinate im klassischen Sinne. Es kann sich nicht von alleine ändern. Wäre das Seil elastisch, müsste man es tatsächlich wie allgemeine Koordinate behandeln, aber nicht in diesem Fall. Die Gleichungen sind nun mal so lang. Ich habe ein Beispielaufruf von craneODE gemacht und die Ergebnisse mit denen aus der Hilfedatei verglichen, die neulich hochgeladen wurde. Der resultierende dqdt-Vektor stimmte überein, also müssen die Lagrange-Gleichungen auch stimmen. Die Gleichungen selbst will ich aber hier nicht posten, weil die schlicht zu lang sind. Aber länger als 25k Zeichen sind sie auch nicht.
Zuletzt geändert von Supakane am Mi, 12. Jun. 13, 19:19, insgesamt 1-mal geändert.
Zum Beispiel mit interp1(ut, u, t) für u. Mit ut und u als Vektoren, wo Zeitpunkte und die entsprechenden Werte gespeichert sind.slmndr hat geschrieben:Könnte mir jemand mal verraten, wie man die symbolischen Variablen und die s-Funktion überhaupt aufstellt? Ich schreibe meine BA außerhalb und kann daher nicht zur Übung
Kannst du nochmal beschreiben was du bzw. ihr mit "in Richtung alpha" meint? Gibt es auch, außer der Reibkraft noch das Motormoment in beta und gamma richtung?slmndr hat geschrieben:Das ist eine nichtkonservative Kraft "in Richtung" alpha:Dennis Eggers hat geschrieben:Gannz allg: wie rechnet ihr die Widerstandkraft des Motormomentes aus ? Bzw. wir realisiert ihr das in der Gleichung >
Q_alpha = Motormoment + Reibungskräfte in alpha-Richtung
Sehe ich das richtig, dass wir die drei Lagrange-Gleichungen nach alpha, beta und gamma aufstellen müssen?
VG
alpha beta gamma sind die richtungen des Vetorraumes, in dem wir arbeiten.Ayanda hat geschrieben:Kannst du nochmal beschreiben was du bzw. ihr mit "in Richtung alpha" meint? Gibt es auch, außer der Reibkraft noch das Motormoment in beta und gamma richtung?slmndr hat geschrieben:Das ist eine nichtkonservative Kraft "in Richtung" alpha:Dennis Eggers hat geschrieben:Gannz allg: wie rechnet ihr die Widerstandkraft des Motormomentes aus ? Bzw. wir realisiert ihr das in der Gleichung >
Q_alpha = Motormoment + Reibungskräfte in alpha-Richtung
Sehe ich das richtig, dass wir die drei Lagrange-Gleichungen nach alpha, beta und gamma aufstellen müssen?
VG
Die lggleichungen warden automatisch in allen richtungen aufgestellt, da du ja den Vektor q vorgegebn hast.
Ok danke.Supakane hat geschrieben:Zum Beispiel mit interp1(ut, u, t) für u. Mit ut und u als Vektoren, wo Zeitpunkte und die entsprechenden Werte gespeichert sind.slmndr hat geschrieben:Könnte mir jemand mal verraten, wie man die symbolischen Variablen und die s-Funktion überhaupt aufstellt? Ich schreibe meine BA außerhalb und kann daher nicht zur Übung
Wie bekomme ich aus s mein ds für die Geschwindigkeit vom Punkt? Ich habe ja s nur als Vektor. Den kann ich nicht ableiten...?
Und noch eine Verständnisfrage. Bekommt craneODE ut, st als VEKTOREN oder sind das SKALARE für bestimmte Zeiten? Beim lersterem kann man ja innerhalb von craneODE diff() nehmen.
--> Hat sich geklärt. Laut Tutor ist ds nicht einzubeziehen. Das macht auch irgendwie Sinn, weil es in die Energie nicht eingeht, da es eine nicht-holonome Zwangsbedingung ist. s und u sind Skalare.
kann sich evtl. jemand mal meine nr. 1 angucken?
wenn ich die simples
d2alpha = simple(solution.d2alpha)
d2beta = simple(solution.d2beta)
d2gamma = simple(solution.d2gamma)
am ende drinne habe, komm er aus irgendeinem grund zu keinem ergebnis, er hört einfach nicht zu rechnen auf.
hat jemand evtl die zeit sich meinen code anzugucken? würd ihn dann über pm schicken
danke und lg
wenn ich die simples
d2alpha = simple(solution.d2alpha)
d2beta = simple(solution.d2beta)
d2gamma = simple(solution.d2gamma)
am ende drinne habe, komm er aus irgendeinem grund zu keinem ergebnis, er hört einfach nicht zu rechnen auf.
hat jemand evtl die zeit sich meinen code anzugucken? würd ihn dann über pm schicken
danke und lg