Voici les valeurs des différents paramêtres du vol:

gv=10                    (accélération du champ de pesanteur terrestre)

M0= 0.570            (masse initiale de la fusée avec moteurs)

Mpoud= 0.029      (masse de poudre par moteur)

Tcomb=1.75          (temps de combustion des moteurs)

       (débit massique de poudre de propulsion, supposé constant pendant la propulsion ???)

rau= 1.22              (masse volumique de l'air)

d= 0.059               (diamêtre de la fusée)

Ca= 0.7                 (coefficient de pénétration dans l'air, sans doute surévalué)

               (le 1.1 permet de tenir compte de l'épaisseur des ailettes de stabilisation)  

L'algorithme de résolution de l'Equation différentielle du Mouvement de la fusée correspond à la méthode d'Euler et a été programmé dans MathCad 6.0. Il permet d'étudier l'évolution des T premières secondes du vol (propulsé puis balistique) lorsque l'engin est propulsé par nmoteur moteurs(1,2 ou 3). L'intervalle de temps pendant lequel gamma est supposé constant vaut 0,01s.

L'algorithme renvoie un tableau (nommé move) de n lignes et 5 colonnes qui permet de tracer les diagrammes d'évolution de l'altitude, de la vitesse, de l'accélération et de la masse en fonction du temps. 

Le source du document MathCad de calcul et d'affichage des caractéristiques cinématiques du vol est disponible au format ZIP. Vous pouvez le télécharger en cliquant ici puis l'exécuter grâce à un module (WorkModel.exe) proposé gracieusement par MathSoft et teléchargeable sur le serveur de MathSoft.