Gelöste Aufgaben/Kerb: Unterschied zwischen den Versionen
Zur Navigation springen
Zur Suche springen
Phasen-Diagramme
Phasengeschwindigkeit
Weg-Zeit-Diagramm
Zeile 21: | Zeile 21: | ||
==tmp== | |||
<!-------------------------------------------------------------------------------->{{MyCodeBlock|title=Header | |||
|text=Ein Satellit im Orbit ist einer Erdbeschleunigung g ausgesetzt, die nichtlinear von seinem Abstand zum Erdmittelpunkt abhängt. Dieses nichtlineare Anfangswertproblem lösen wir hier als [[Anfangswertprobleme|Anfangswertproblem]] numerisch. | |||
|code= | |||
<syntaxhighlight lang="lisp" line start=1> | |||
/*********************************************************/ | |||
/* MAXIMA script */ | |||
/* version: wxMaxima 15.08.2 */ | |||
/* author: Andreas Baumgart */ | |||
/* last updated: 2018-03-21 */ | |||
/* ref: Kerb (TM-C, Labor 6) */ | |||
/* description: finds persiod solution for */ | |||
/* the tracectory of a satellite */ | |||
/*********************************************************/ | |||
</syntaxhighlight> | |||
}} | |||
==tmp== | |||
<!-------------------------------------------------------------------------------->{{MyCodeBlock|title=Declarations | |||
|text= | |||
Wir brauchen die Systemparameter | |||
{| class="wikitable" | |||
! Symbol !! Bedeutung | |||
|- | |||
|R = 6378,137 km||Erdradius am Äquator | |||
|- | |||
|g<sub>0</sub> = 9.81 m/s2||Erdbeschleunigung an der Erdoberfläche | |||
|- | |||
|T = 1 d||Periodendauer einer Erdumdrehung (d .. day) | |||
|O<sub>r</sub> = 35786 km||Höhe eines Satelliten über der Erdoberfläche in einem geostationären Orbit. | |||
|- | |||
|} | |||
|code= | |||
<syntaxhighlight lang="lisp" line start=1> | |||
/*********************************************************/ | |||
/* declarations */ | |||
/* declare variational variables - see 6.3 Identifiers */ | |||
declare("δW", alphabetic); | |||
declare("δr", alphabetic); | |||
declare("δu", alphabetic); | |||
declare("δφ", alphabetic); | |||
assume(R>0, g>0, rho>0, mu>0, T>0); | |||
/* parameter */ | |||
params: [R = 6378137*m, /* earth radius at equator */ | |||
g = 9.81*m/s^2, /* gravitational constant */ | |||
T = 24*60*60*s]; /* 24 hours */ | |||
verify: [O=R+35786000*m]; /* geostat. orbit */ | |||
dimles: [g=mu*R/T^2]; /* dimensionless parameter */ | |||
</syntaxhighlight> | |||
}} | |||
==tmp== | |||
<!-------------------------------------------------------------------------------->{{MyCodeBlock|title=Kinematics | |||
|text=.. | |||
|code= | |||
<syntaxhighlight lang="lisp" line start=1> | |||
XX | |||
</syntaxhighlight> | |||
}} | |||
==tmp== | |||
<!-------------------------------------------------------------------------------->{{MyCodeBlock|title=Equilibrium Conditions | |||
|text=.. | |||
|code= | |||
<syntaxhighlight lang="lisp" line start=1> | |||
XX | |||
</syntaxhighlight> | |||
}} | |||
==tmp== | |||
<!-------------------------------------------------------------------------------->{{MyCodeBlock|title=Stationary Solution | |||
|text=.. | |||
|code= | |||
<syntaxhighlight lang="lisp" line start=1> | |||
XX | |||
</syntaxhighlight> | |||
}} | |||
==tmp== | |||
<!-------------------------------------------------------------------------------->{{MyCodeBlock|title=Solving | |||
|text=.. | |||
|code= | |||
<syntaxhighlight lang="lisp" line start=1> | |||
XX | |||
</syntaxhighlight> | |||
}} | |||
<!-------------------------------------------------------------------------------->{{MyCodeBlock|title= | ==tmp== | ||
|text= | <!-------------------------------------------------------------------------------->{{MyCodeBlock|title=Post-Processing | ||
|text=.. | |||
|code= | |code= | ||
<syntaxhighlight lang="lisp" line start=1> | <syntaxhighlight lang="lisp" line start=1> | ||
XX | |||
</syntaxhighlight> | </syntaxhighlight> | ||
}} | }} | ||
<hr/> | |||
'''Links''' | '''Links''' | ||
* ... | * ... |
Version vom 25. März 2021, 08:43 Uhr
Aufgabenstellung
Ein Satellit soll eine stabile Umlaufbahn um die Erde beschrieben.
Berechnen Sie mögliche Lösungen.
Lösung mit Maxima
Bei unserer Suche geht es um kartesische- und polare-Koordinatensysteme, Gleichgewichtsbedingungen formuliert mit dem Prinzip der virtuellen Verrückungen und er Stabilität von Lösungen.
tmp
Header
Ein Satellit im Orbit ist einer Erdbeschleunigung g ausgesetzt, die nichtlinear von seinem Abstand zum Erdmittelpunkt abhängt. Dieses nichtlineare Anfangswertproblem lösen wir hier als Anfangswertproblem numerisch.
/*********************************************************/
/* MAXIMA script */
/* version: wxMaxima 15.08.2 */
/* author: Andreas Baumgart */
/* last updated: 2018-03-21 */
/* ref: Kerb (TM-C, Labor 6) */
/* description: finds persiod solution for */
/* the tracectory of a satellite */
/*********************************************************/
tmp
Declarations
Wir brauchen die Systemparameter {
/*********************************************************/
/* declarations */
/* declare variational variables - see 6.3 Identifiers */
declare("δW", alphabetic);
declare("δr", alphabetic);
declare("δu", alphabetic);
declare("δφ", alphabetic);
assume(R>0, g>0, rho>0, mu>0, T>0);
/* parameter */
params: [R = 6378137*m, /* earth radius at equator */
g = 9.81*m/s^2, /* gravitational constant */
T = 24*60*60*s]; /* 24 hours */
verify: [O=R+35786000*m]; /* geostat. orbit */
dimles: [g=mu*R/T^2]; /* dimensionless parameter */
tmp
Kinematics
..
XX
tmp
Equilibrium Conditions
..
XX
tmp
Stationary Solution
..
XX
tmp
Solving
..
XX
tmp
Post-Processing
..
XX
Links
- ...
Literature
- ...