onizuka_4wyf08e9he33sgce_333.py

🔷 INFORMATIONS GÉNÉRALES

Champ Détail

Nom du script onizuka_4wyf08e9he33sgce_333.py

Version 1.0.0

Date de création 25/03/2026

Auteur / Demandeur refoïa jeremy

Objectif principal Simuler et visualiser la trajectoire 2D d’une fusée prenant en compte la gravité variable, la traînée atmosphérique et la consommation de carburant.

🎯 DESCRIPTION FONCTIONNELLE

📌 Que fait ce script ?

Le script résout un système d’équations différentielles ordinaires (ODE) pour modéliser le vol d’une fusée de type Falcon 9. Il calcule en temps réel la position, la vitesse et la masse de l’engin en fonction de la poussée des moteurs et des forces environnementales (gravité de la Terre et résistance de l’air).

📌 Problème résolu

Il permet de prédire la trajectoire d’un lanceur spatial, d’estimer l’altitude maximale atteinte et de visualiser l’impact du profil de tangage (pitch program) sur la mise en orbite théorique.

⚙️ SPÉCIFICATIONS TECHNIQUES

🐍 Environnement

Élément Valeur

Version Python 3.x

OS cible Tous (Windows / Linux / MacOS)

Mode d’exécution Script d’analyse de données / Simulation

📦 Dépendances / Librairies

numpy : Calculs numériques et gestion des vecteurs.

matplotlib : Génération des graphiques de performance.

scipy (odeint) : Intégration numérique des équations du mouvement.

📥 ENTRÉES (INPUTS)

📂 Paramètres de simulation (Hardcodés)

# Nom Type Description Exemple

1 M_FUEL float Masse initiale de carburant (kg) 400 000

2 THRUST float Poussée constante des moteurs (N) 7.6e6

3 ISP float Impulsion spécifique (efficacité) 300

4 CD float Coefficient de traînée aérodynamique 0.5

📤 SORTIES (OUTPUTS)

📊 Visualisation (Matplotlib)

Le script génère une fenêtre avec 4 graphiques :

Altitude en fonction du temps (km).

Vitesse en fonction du temps (m/s).

Consommation Carburant (Évolution de la masse totale).

Trajectoire (x, z) : Vue en coupe de l’ascension de la fusée.

🧱 STRUCTURE DU SCRIPT

onizuka_…_333.py

│

├── 📌 CONSTANTES PHYSIQUES (G, Masse Terre, Rayon Terre, Atmosphère)

├── 📌 PARAMÈTRES FUSÉE (Masse, Poussée, ISP, Surface)

├── 📌 FONCTION model(y, t) (Cœur de la simulation)

│ ├── Calcul Gravité (loi en 1/r2)

│ ├── Modèle Atmosphérique (isotherme exponentiel)

│ ├── Pilotage (Pitch Program)

│ └── Équations du mouvement (PFD)

└── 📌 MAIN / SIMULATION

├── Intégration avec odeint

└── Génération des graphiques plt.subplot

🔄 LOGIQUE / ALGORITHME

Initialisation : Définition des conditions initiales (x=0,z=0, vitesse nulle, plein de carburant).

Boucle d’Intégration (Toutes les dt) :

Calcul de l’altitude par rapport au centre de la Terre.

Calcul de la densité de l’air ρ=ρ0​⋅e−alt/Hscale​.

Phase de Pitch : Verticale jusqu’à 2km, puis inclinaison exponentielle vers l’horizontale.

Bilan des forces : ∑F=Pousseˊe+Traı^neˊe+Poids=ma.

Mise à jour de la masse (m=m−deˊbit⋅dt).

Résolution : odeint compile ces étapes pour 1000 points temporels.

Rendu : Conversion des vecteurs d’état en graphiques lisibles.

🚨 GESTION DES ERREURS

Masse négative : Le script vérifie if m > M_DRY pour stopper la poussée et la consommation si le réservoir est vide.

Atmosphère : Une condition if alt < 100000 évite les calculs de traînée inutiles ou les erreurs numériques dans le vide spatial (au-delà de 100km). ✅ CONTRAINTES & RÈGLES MÉTIER Règle 1 : La gravité n'est pas constante (9.81); elle diminue avec le carré de la distance au centre de la Terre. Règle 2 : La poussée s'arrête instantanément dès que la masse de carburant est épuisée. Règle 3 : Le profil de tangage est "aveugle" (basé uniquement sur l'altitude et non sur des vecteurs de vitesse). 📝 EXEMPLE D'UTILISATION ▶️ Lancement en ligne de commande Bash python onizuka_4wyf08e9he33sgce_333.py Une fenêtre s'ouvrira affichant les 4 courbes de télémétrie après quelques secondes de calcul.