# ESA-PB05-artiger Bahnpropagator mit mittleren Elementen (Brouwer–Lyddane, J2–J6 + lunisolare)

Propagiert erdnahe und GNSS-Bahnen mit einem semianalytischen Mittelerelement-Modell in der Tradition Brouwer–Lyddane/PB05: bahngemittelte säkulare Drift unter dem EGM96-Zonalfeld J2–J6 (numerische Gauss-Mittelung — das numerische Analogon der Delaunay-Normalisierung mit den geschlossenen J2-Raten als eingebauter Gegenprüfung), Rekonstruktion erster Ordnung J2 kurzperiodisch mittel↔oskulierend, J3-Frozen-Orbit-Analyse, doppelt gemittelter lunisolarer Quadrupol, Resonanzprüfungen der Bodenspur, heliosynchrone Bahnneigungs-Lösung und RK4-Numerik-Kreuzvalidierung im Bericht.

> Kanonische Seite: https://elysiatools.com/de/tools/esa-pb05-pce-pck-reddy-porter-keplerian-equation-averaged-element-propagator

- **Kategorie:** Astronomy

- **Schlagwörter:** mittelerelemente-propagator, brouwer lyddane, j2-sakulare raten, knoten-drift, frozen orbit, kozai lidov, lunisolare storung, sonnensynchrone bahn, tesseral-resonanz, mittel zu oskulierend

## Überblick

Der Propagator integriert die MITTLEREN Elemente per RK4, deren rechte Seite das Bahnmittel der Gaußschen Variationsraten über das Zonalfeld der gewählten Tiefe ist — das numerische Äquivalent der Delaunay-Normalisierung nach Brouwer–Lyddane/PB05, abgeglichen mit den geschlossenen J2-Raten erster Ordnung. Die lunisolare Option addiert den quadrupolartigen Drittkörperterm mit Doppelmittelung (über Satelliten- und Störerbahn), der das klassische GEO-Bahnneigungswachstum von rund 0,75–0,95°/Jahr reproduziert und die kritische gegenseitige Neigung von Kozai–Lidov 39,23° enthält. Kurzperiodische J2-Abweichungen werden aus den momentanen Raten integriert und können zur Abtastung der oskulierenden Elemente zurückaddiert werden.

## Eingaben

- **Große Halbachse a (km)** (number): 8000
- **Exzentrizität e** (number): 0.05
- **Bahnneigung i (Grad)** (number): 52.5
- **Rektaszension des Knotens Ω (Grad)** (number): 30
- **Perigäumsargument ω (Grad)** (number): 70
- **Mittlere Anomalie M (Grad)** (number): 0
- **Propagationsdauer (Tage)** (number): 30
- **Zonalfeld-Tiefe** (select)
- **Drittkörper-Störungen** (select)
- **Panel Mittel→oskulierend hinzufügen** (checkbox)

## Wann verwenden

- Bei der Missionsplanung zur Abschätzung langfristiger säkularer Knoten- und Perigäumsdriften von LEO-, MEO- oder GEO-Satelliten.
- Zur Auslegung sonnensynchroner Orbits (SSO) oder geodätischer Frozen Orbits unter Berücksichtigung höherer Zonalharmonischer (J2–J6).
- Zur Analyse von Bahnneigungsdriften geostationärer Bahnen durch kombinierte lunisolare Gravitationsstörungen über Monate oder Jahre.

## Funktionsweise

- Sie definieren die anfänglichen Bahnelemente (große Halbachse, Exzentrizität, Inklination, RAAN, Argument des Perigäums, mittlere Anomalie) sowie die Propagationsdauer in Tagen.
- Wählen Sie die gewünschte Zonalfeld-Tiefe (J2, J2–J4 oder J2–J6) und aktivieren Sie bei Bedarf lunisolare Drittkörper-Störungen (Sonne, Mond oder beide) über ein doppelt gemitteltes Quadrupolmodell.
- Das System integriert die mittleren Elemente mittels RK4 über die Gaußschen Bewegungsgleichungen und gleicht die Ergebnisse mit analytischen J2-Raten sowie einer numerischen kartesischen RK4-Referenzbahn ab.
- Der generierte HTML-Bericht visualisiert den zeitlichen Verlauf der Bahnelemente, listet Bodenspurresonanzen auf und stellt auf Wunsch die Rücktransformation von mittleren zu oskulierenden Elementen dar.

## Anwendungsfälle

- Langzeit-Driftanalyse der Knotengeraden (RAAN) für Erdbeobachtungs- und SAR-Satelliten im sonnensynchronen Orbit.
- Abschätzung des jährlichen Korrekturbedarfs für die Bahnneigung (North-South Stationkeeping) von GEO-Nachrichtensatelliten.
- Prüfung von Resonanzbedingungen zwischen Satellitenumlaufzeit und Erdrotation bei GNSS-Konstellationen im MEO-Bereich.

## Häufig gestellte Fragen

### Was unterscheidet mittlere Bahnelemente von oskulierenden Elementen?

Oskulierende Elemente beschreiben die exakte Momentanellipse inklusive kurzperiodischer Oszillationen, während mittlere Elemente diese Schwingungen herausmitteln und die reine säkulare sowie langperiodische Drift abbilden.

### Welche zonalen Schwerefeldkoeffizienten werden berücksichtigt?

Je nach Auswahl werden J2 (erste Ordnung), J2 bis J4 oder der vollständige Satz zonaler Harmonischer von J2 bis J6 aus dem EGM96-Schwerefeldmodell ausgewertet.

### Wie wird die lunisolare Störung modelliert?

Die Störung durch Sonne und Mond wird über ein doppelt gemitteltes Quadrupol-Potential berechnet, das sowohl über die Satellitenumlaufbahn als auch über die Bahn des Störkörpers gemittelt ist.

### Kann das Tool Bedingungen für sonnensynchrone Bahnen ermitteln?

Ja, der Bericht berechnet automatisch die für die gegebene Bahnhöhe erforderliche heliosynchrone Bahnneigung und gleicht die berechnete RAAN-Driftrate damit ab.

### Wozu dient die integrierte RK4-Numerik-Kreuzvalidierung?

Sie führt eine unabhängige numerische Integration der Bewegungsgleichungen durch, um die Abweichung zwischen dem semianalytischen Mittelwertmodell und der Punktmassen- bzw. Vollfelddynamik transparent auszuweisen.

## Ähnliche Tools

- [Korrelationsanalyse](https://elysiatools.com/de/tools/correlation-analyzer): Fortschrittetes Korrelationsanalyse-Werkzeug zur Berechnung von Korrelationskoeffizienten zwischen Variablen zur Messung der Stärke und Richtung ihrer linearen Beziehungen.
- [Bild-Palette zu Design-Tokens](https://elysiatools.com/de/tools/image-to-design-tokens): Extrahiert dominante Farben per k-Means und exportiert sie als CSS-/SCSS-Variablen, Tailwind-Config oder JSON-Tokens mit Namen und Farbskala
- [Berater Mieten vs. Kaufen — Gesamtbetriebskosten (Hypothek/Steuern/Versicherung/Instandhaltung + NPV-Break-even + IRR)](https://elysiatools.com/de/tools/rent-vs-buy-total-cost-of-ownership-mortgage-tax-insurance-maintenance-irr-tutor): Monatliche Nebeneinander-Simulation: Hypothekentilgung + Grundsteuer/Versicherung/HOA/Instandhaltung gegen Miete + Mieterversicherung; der Mieter investiert Anzahlung und monatliche Ersparnis zum Opportunitätssatz; Ausgabe: jährlicher Kapitalwert, Break-even-Jahr, IRR der Differenzflüsse, SALT-Deckel-Steuerschild und ±20-%-Tornado-Diagramm.
- [Örtlicher Druckverlustrechner (h = K·v²/2g)](https://elysiatools.com/de/tools/minor-loss-calculator): Berechnet die örtlichen Druckverluste in Formstücken, Ventilen und Bögen der Rohrleitung: h_m=K·v²/(2g), ΔP_m=K·ρ·v²/2. Enthält typische K-Werte für 14 gängige Bauteile (Schieber/Durchgangs-/Rückschlagventile, 90°/45°-Knie, T-Stücke, plötzliche Erweiterungen/Verengungen, Ein-/Austritte). Bauteiltyp wählbar oder K frei eingebbar. Unterstützt n gleiche Bauteile in Reihe; ausgabeseitig werden Höhe (m/ft) und Druckabfall (Pa/kPa/bar) ausgegeben. Dichte in kg/m³, g/cm³ oder lb/ft³.
- [Rohr-Druckverlustrechner nach Darcy-Weisbach (ΔP = f·(L/D)·ρv²/2)](https://elysiatools.com/de/tools/pipe-pressure-drop-darcy): Berechnet den Druckverlust durch Reibung nach Darcy-Weisbach in einem geraden Rohr: ΔP=f·(L/D)·(ρ·v²/2) Pa; Druckhöhenverlust h_f=ΔP/(ρ·g)=f·(L/D)·v²/(2g) m. Der Benutzer gibt den Darcy-Reibungsfaktor f ein (nicht Fanning). Länge in m/km/ft, Durchmesser in m/cm/mm/inch, Dichte in kg/m³/g/cm³/lb/ft³ (auf SI umgerechnet). Ausgabe in kPa, bar sowie Höhe in m/ft.
- [Daten-Grenzverarbeiter](https://elysiatools.com/de/tools/data-boundary-processor): Erweitertes Grenzwertverarbeitungswerkzeug zur Identifizierung und Handhabung von Minimum- und Maximumwerten in numerischen Daten. Perfekt für Datenvalidierung, Bereichsprüfung, statistische Analyse und Datenvorverarbeitung.
- [Strom-/Gas-/Wasserrechnungs-Decoder und Detektor versteckter Gebühren (CSR/CFL/MHC/HFA-Rider · Spitzenlast · kWh/Btu/therm-Umrechnung)](https://elysiatools.com/de/tools/csr-cfl-cfl-mhc-hfa-btu-ton-hours-utility-bill-line-item-validator-and-rgp-peak-explainer): Fügen Sie Rechnungspositionen ein: Klassifizierung und Erklärung von Grundpreis, Leistungsentgelt, Hoch-/Niedertarif, Transport-/Verteilungs-Ridern, Erneuerbaren-Fonds, Decoupling, Capacity Tag; Umrechnung kWh↔Btu↔therm↔Ccf↔m³↔Tonnenstunden nach EIA-Faktoren; versteckte Fixkosten und Einsparprojektionen.
- [Dichterechner (ρ = m/V)](https://elysiatools.com/de/tools/density-calculator): Berechnet Dichte, Masse oder Volumen aus den anderen beiden, mit spezifischer Dichte und Schwimmverhalten

## Beispiele

- [Web Python Bildverarbeitung Beispiele](https://elysiatools.com/de/samples/web-image-processing-python): Web Python Bildverarbeitungsbeispiele mit PIL/Pillow einschließlich Lesen, Speichern, Skalieren und Formatkonvertierung
- [Android Java Bildverarbeitungsbeispiele](https://elysiatools.com/de/samples/android-image-processing-java): Android Java Bildverarbeitungsbeispiele einschließlich Lesen/Schreiben, Skalierung und Formatkonvertierung
- [Android Kotlin Bildverarbeitungsbeispiele](https://elysiatools.com/de/samples/android-image-processing-kotlin): Android Kotlin Bildverarbeitungsbeispiele einschließlich Lesen/Schreiben, Skalierung und Formatkonvertierung
- [Bildverarbeitung macOS Objective-C - Beispiele](https://elysiatools.com/de/samples/macos-image-processing-objectivec): macOS Objective-C Bildverarbeitungsbeispiele einschließlich Bild-Lesen/Schreiben, Bild-Skalierung und Formatkonvertierung