# Propagateur d’orbite à éléments moyens de style ESA PB05 (Brouwer–Lyddane, J2–J6 + luni-solaire)

Propage des orbites proches de la Terre et de la région GNSS avec un modèle semi-analytique d’éléments moyens dans la tradition Brouwer–Lyddane/PB05 : dérive séculaire moyennée sur l’orbite sous le champ zonal EGM96 J2–J6 (moyennage numérique des équations de Gauss — l’analogue numérique de la normalisation de Delaunay avec les formules fermées J2 comme contre-vérification intégrée), reconstruction du premier ordre J2 moyen↔osculant, analyse d’orbite gelée J3, quadrupôle luni-solaire doublement moyenné, contrôles de résonance tessérale de la trace au sol, résolution de l’inclinaison héliosynchrone et validation croisée numérique RK4 dans le rapport.

> Page canonique: https://elysiatools.com/fr/tools/esa-pb05-pce-pck-reddy-porter-keplerian-equation-averaged-element-propagator

- **Catégorie:** Astronomy

- **Mots-clés:** propagateur elements moyens, brouwer lyddane, taux seculaires j2, derivee du noeud, orbite gelee, kozai lidov, perturbation luni-solaire, orbite heliosynchrone, resonance tesseral, moyen vers osculant

## Présentation

Le propagateur intègre les éléments MOYENS par RK4 dont le second membre est la moyenne orbitale des taux variationnels de Gauss sur le champ zonal à la profondeur choisie — l’équivalent numérique de la normalisation de Delaunay de Brouwer–Lyddane/PB05, recoupé avec les formules fermées du premier ordre J2. L’option luni-solaire ajoute le terme de troisième corps quadrupôle à double moyenne (sur l’orbite du satellite et celle du perturbateur), qui reproduit la croissance classique d’inclinaison GEO d’environ 0,75–0,95°/an et porte l’inclinaison mutuelle critique de Kozai–Lidov à 39,23°. Les écarts de courte période J2 sont intégrés à partir des taux instantanés et peuvent être réajoutés pour échantillonner les éléments osculants. Le rapport exécute aussi une propagation numérique cartésienne RK4 indépendante avec le même modèle de forces et affiche leur écart.

## Entrées

- **Demi-grand axe a (km)** (number): 8000
- **Excentricité e** (number): 0.05
- **Inclinaison i (degrés)** (number): 52.5
- **Ascension droite du nœud Ω (degrés)** (number): 30
- **Argument du périgée ω (degrés)** (number): 70
- **Anomalie moyenne M (degrés)** (number): 0
- **Durée de propagation (jours)** (number): 30
- **Profondeur du champ zonal** (select)
- **Perturbations tierces** (select)
- **Ajouter le panneau de reconstruction moyen→osculant** (checkbox)

## Quand l'utiliser

- Évaluer la dérive séculaire du nœud ascendant (RAAN) et de l'argument du périgée sur plusieurs semaines ou mois.
- Analyser la stabilité à long terme d'orbites géostationnaires ou MEO soumises aux effets de marée du Soleil et de la Lune.
- Calculer les conditions de résonance tessérale, d'orbite gelée J3 ou d'inclinaison héliosynchrone pour la planification de mission.

## Fonctionnement

- L'outil intègre numériquement par Runge-Kutta 4 (RK4) les équations variationnelles de Gauss moyennées sur une période orbitale pour les harmoniques zonaux sélectionnés (J2, J2–J4 ou J2–J6).
- Il superpose le terme quadrupolaire luni-solaire à double moyenne pour modéliser les effets séculaires de troisième corps et les dynamiques de type Kozai-Lidov.
- Il reconstruit les éléments osculants à partir des éléments moyens via les corrections périodiques courtes de J2 et valide les résultats par une propagation RK4 cartésienne indépendante.

## Cas d'usage

- Dimensionnement de la dérive de RAAN pour les constellations de satellites d'observation héliosynchrones en orbite basse (LEO).
- Prévision de la dérive en inclinaison induite par la Lune et le Soleil pour le maintien à poste des satellites géostationnaires (GEO).
- Conception d'orbites à périgée figé (orbites gelées) pour minimiser les variations d'altitude au-dessus de zones cibles.

## Questions fréquentes

### Quelle est la différence entre éléments moyens et éléments osculants ?

Les éléments osculants décrivent la trajectoire instantanée réelle avec ses oscillations à court terme, tandis que les éléments moyens filtrent ces variations pour isoler les dérives séculaires et de longue période.

### Quels harmoniques zonaux peuvent être inclus dans la propagation ?

Vous pouvez configurer la profondeur zonale sur J2 seul, J2 à J4, ou le jeu complet EGM96 de J2 à J6.

### Comment les perturbations du Soleil et de la Lune sont-elles modélisées ?

Elles sont intégrées sous forme d'un terme quadrupolaire séculaire à double moyenne, moyenné à la fois sur l'orbite du satellite et sur celle du corps perturbateur.

### Comment l'outil identifie-t-il les résonances de trace au sol ?

Il compare le moyen mouvement moyen du satellite à la vitesse de rotation terrestre pour repérer les commensurabilités rationnelles telles que 1:1 (GEO), 2:1 (Molniya/MEO) ou 15:1 (LEO).

### Le rapport fournit-il une vérification numérique ?

Oui, une propagation numérique cartésienne RK4 directe est exécutée en parallèle pour quantifier l'écart résiduel avec le modèle semi-analytique moyen.

## Outils associés

- [Analyseur de Corrélation](https://elysiatools.com/fr/tools/correlation-analyzer): Outil avancé d'analyse de corrélation qui calcule les coefficients de corrélation entre les variables pour mesurer la force et la direction de leurs relations linéaires.
- [De la palette d’image aux design tokens](https://elysiatools.com/fr/tools/image-to-design-tokens): Extrait les couleurs dominantes d’une image par k-means et les exporte en variables CSS / SCSS / config Tailwind / tokens JSON, avec noms et échelle de teintes
- [Tuteur coût total locataire vs propriétaire (crédit/impôts/assurance/entretien + seuil d'équilibre VAN + TRI)](https://elysiatools.com/fr/tools/rent-vs-buy-total-cost-of-ownership-mortgage-tax-insurance-maintenance-irr-tutor): Simulation mensuelle côte à côte : amortissement du crédit + taxe/assurance/copropriété/entretien contre loyer + assurance locataire ; le locataire place l'apport et l'épargne mensuelle au taux d'opportunité ; sortie : VAN annuelle, année d'équilibre, TRI des flux différentiels, bouclier fiscal sous plafond SALT et diagramme tornado ±20 %.
- [Calculateur de Pertes de Charge Locales (h = K·v²/2g)](https://elysiatools.com/fr/tools/minor-loss-calculator): Calcule les pertes de charge locales dans les raccords, vannes et coudes de tuyauterie : h_m=K·v²/(2g), ΔP_m=K·ρ·v²/2. Inclut les K typiques de 14 pièces courantes (vannes à porte/globe/clapet, coudes 90°/45°, tés, élargissements/contractions brusques, entrées/sorties). Sélection possible du type ou K personnalisé. Supporte n pièces identiques en série ; renvoie la hauteur (m/ft) et la chute de pression (Pa/kPa/bar). Masse volumique en kg/m³, g/cm³ ou lb/ft³.
- [Calculateur de Perte de Charge Darcy-Weisbach (ΔP = f·(L/D)·ρv²/2)](https://elysiatools.com/fr/tools/pipe-pressure-drop-darcy): Calcule la perte de charge principale (par frottement) de Darcy-Weisbach dans une conduite droite : ΔP=f·(L/D)·(ρ·v²/2) Pa ; perte de charge h_f=ΔP/(ρ·g)=f·(L/D)·v²/(2g) m. L'utilisateur saisit le facteur de frottement de Darcy f (pas Fanning). Longueur en m/km/ft, diamètre en m/cm/mm/inch, masse volumique en kg/m³/g/cm³/lb/ft³ (convertie en SI). Renvoie kPa, bar et charge en m/ft.
- [Processeur de Limites de Données](https://elysiatools.com/fr/tools/data-boundary-processor): Outil avancé de traitement des limites qui identifie et gère les valeurs minimales et maximales dans les données numériques. Parfait pour la validation des données, la vérification des plages, l'analyse statistique et le prétraitement des données.
- [Décodeur de factures d'énergie et détecteur de frais cachés (riders CSR/CFL/MHC/HFA · demande heure pleine/creuse · conversion kWh/Btu/therm)](https://elysiatools.com/fr/tools/csr-cfl-cfl-mhc-hfa-btu-ton-hours-utility-bill-line-item-validator-and-rgp-peak-explainer): Collez les lignes de votre facture : classification et explication du fixed customer charge, de la demande, des prix heures pleines/creuses, des riders transport/distribution, fonds renouvelables, découplage, capacity tag ; conversion kWh↔Btu↔therm↔Ccf↔m³↔tonne-heure aux facteurs EIA ; détection des frais fixes cachés et projections d'économies.
- [Calculateur de Densité (ρ = m/V)](https://elysiatools.com/fr/tools/density-calculator): Calcule la densité, la masse ou le volume à partir des deux autres, avec gravité spécifique et flottaison

## Exemples

- [Exemples de Traitement d'Images Web Python](https://elysiatools.com/fr/samples/web-image-processing-python): Exemples de traitement d'images Web Python utilisant PIL/Pillow incluant la lecture, l'enregistrement, le redimensionnement et la conversion de format
- [Exemples de Traitement d'Images Android Java](https://elysiatools.com/fr/samples/android-image-processing-java): Exemples de traitement d'images Android Java incluant lecture/écriture, mise à l'échelle et conversion de format
- [Exemples de Traitement d'Images Android Kotlin](https://elysiatools.com/fr/samples/android-image-processing-kotlin): Exemples de traitement d'images Android Kotlin incluant lecture/écriture, mise à l'échelle et conversion de format
- [Traitement d'Images macOS Objective-C - Exemples](https://elysiatools.com/fr/samples/macos-image-processing-objectivec): Exemples traitement d'images macOS Objective-C incluant lecture/écriture images, redimensionnement images et conversion format