# Trajectoire d'un projectile

En physique, la trajectoire d'un projectile est le chemin que ce projectile empruntera après avoir été lancé.
Par projectile, on entend un objet petit et à masse réduite (par exemple un caillou ou un boulet de canon) par rapport à la taille et à la masse du corps depuis lequel il est propulsé (ce peut être la Terre, la Lune, etc.) à une vitesse et avec un angle tels qu'il ne puisse pas se maintenir en orbite et qu'il retombe. En outre, il faut que sa vitesse soit suffisamment faible par rapport à la vitesse de la lumière — que le facteur de Lorentz qui la caractérise soit suffisamment proche de 1 — pour ne pas devoir quitter les considérations de la physique classique.
Lorsque l’on étudie la trajectoire d’un projectile en négligeant toutes les forces à l'exception du poids, on parle de trajectoire balistique. La résolution des équations de la trajectoire balistique d'un projectile est aisée, mais cette modélisation est éloignée de la réalité lorsque d'autres forces sont considérées, par exemple les forces induites par le fluide dans lequel évolue le projectile (forces de frottement, poussée d'Archimède, etc.).

## Notations

Dans les équations de l'article, les variables suivantes seront utilisées :
* g : l'accélération de la pesanteur en m/s$^2$ (9,81 m/s$^2$ à la surface de la Terre) ;
* θ ou α : l'angle de portée, c'est-à-dire l'angle avec lequel le projectile est lancé, en degrés ;
* v$_0$ ou v$_i$ : la vitesse initiale, à laquelle le projectile est lancé, en m/s ;
* y$_0$ ou h : la hauteur initiale du projectile, en mètres ;
* d : la distance horizontale parcourue par le projectile, du point de lancement au point de tombée au sol, en mètres.
* k : le coefficient de proportionnalité caractéristique de la traînée du fluide.
Aussi, pour exprimer les distances, on pourra utiliser un repère cartésien orthonormé ${\mathcal {R}}=(O;{\vec {i}},{\vec {j}},{\vec {k}})$ où O est le point d'où est lancé le projectile. On appellera l'axe de coordonnées $(Ox)=(O;{\vec {i}})$ l'horizontale, l'axe de coordonnées $(Oy)=(O;{\vec {j}})$ la verticale, et l'axe de coordonnées $(Oz)=(O;{\vec {k}})$ la profondeur.

## Trajectoire balistique

La trajectoire balistique d'un projectile est le chemin que prendra un projectile soumis à la seule action de la gravitation, en négligeant toutes les autres forces.

### Équations horaires et cartésiennes

Théorème (trajectoire balistique) — Dans ${\mathcal {R}}$, si l'on néglige toutes les forces à l'exception du poids, un projectile lancé à t=0 seconde d'une hauteur h (en mètre) à une vitesse v$_0$ (en m/s) en faisant un angle α (en radians) entre l'horizontale et la verticale, et un angle β (en radians) entre la profondeur et l'horizontale, a pour coordonnées (en mètres) au cours du temps (en secondes) :
    $x(t)=\cos(\alpha )v_{0}t$
    $y(t)=-{\frac {1}{2}}gt^{2}+\sin(\alpha )v_{0}t+h$
    $z(t)=\cos(\beta )v_{0}t$.
Équation cartésienne du projectile (trajectoire balistique) — Dans ${\mathcal {R}}$, si l'on néglige toutes les forces à l'exception du poids, un projectile lancé d'une hauteur h (en mètres) à une vitesse v$_0$ (en m/s) en faisant un angle α (en radians) entre l'horizontale et la verticale, et un angle β (en radians) entre la profondeur et l'horizontale, son altitude (en mètres) et sa profondeur (en mètres) en fonction de la distance qu'il a parcouru (en mètres) est :
    $y(x)=-{\frac {1}{2}}{\frac {gx^{2}}{\cos ^{2}(\alpha )v_{0}^{2}}}+\tan(\alpha )x+h$
    $z(x)=x{\frac {\cos(\beta )}{\cos(\alpha )}}$.
On considère un projectile lancé à la surface d'un corps dont le champ de pesanteur est considéré comme uniforme et on se situe dans un référentiel galiléen. Ainsi, la deuxième loi de Newton s'applique à notre système, et celle-ci nous dit que la somme des forces extérieurs qui s'appliquent sur le projectile est égale à la dérivée de sa quantité de mouvement. On a donc
    $\sum _{i}{\vec {F}}_{i}={\frac {\mathrm {d} {\vec {p}}}{\mathrm {d} t}}={\frac {\mathrm {d} \left(m{\vec {v}}\right)}{\mathrm {d} t}}$
où ${\vec {F}}_{i}$ est une force qui s'applique au projectile et ${\vec {p}}=m{\vec {v}}$ la quantité de mouvement du projectile.
Le projectile ne change pas de masse au cours de son lancer, ainsi on peut écrire que
    $\sum _{i}{\vec {F}}_{i}=m{\frac {\mathrm {d} {\vec {v}}}{\mathrm {d} t}}$.
Or, par définition, la dérivée de la vitesse par rapport au temps est l'accélération, donc ${\frac {\mathrm {d} {\vec {v}}}{\mathrm {d} t}}={\vec {a}}$, soit
    $\sum _{i}{\vec {F}}_{i}=m{\vec {a}}$.
D'autre part, toutes les forces sont négligeables par rapport au poids ${\vec {P}}=m{\vec {g}}$. On a donc
    $\sum _{i}{\vec {F}}_{i}={\vec {P}}=m{\vec {g}}=m{\vec {a}}$
d'où
    ${\vec {a}}={\vec {g}}={\begin{pmatrix}0\\-g\\0\end{pmatrix}}$.
On a donc le vecteur accélération du projectile. Pour déterminer son vecteur vitesse, il suffit de trouver la primitive du vecteur accélération par rapport au temps, soit
    ${\vec {v}}={\begin{pmatrix}{C_{x}}_{v}\\-gt+{C_{y}}_{v}\\{C_{z}}_{v}\end{pmatrix}}$ où ${C_{x}}_{v}$, ${C_{y}}_{v}$ et ${C_{z}}_{v}$ sont des constantes.
Or à t=0, on a ${\vec {v}}={\begin{pmatrix}{C_{x}}_{v}\\-g\times 0+{C_{y}}_{v}\\{C_{z}}_{v}\end{pmatrix}}={\begin{pmatrix}{C_{x}}_{v}\\{C_{y}}_{v}\\{C_{z}}_{v}\end{pmatrix}}={\vec {v_{0}}}$.
Et on connaît les valeurs des composantes de ${\vec {v_{0}}}$. En effet, par trigonométrie, on sait que :
    $\cos(\alpha )={\frac {v_{0_{x}}}{v_{0}}}$
    $\sin(\alpha )={\frac {v_{0_{y}}}{v_{0}}}$
    $\cos(\beta )={\frac {v_{0_{z}}}{v_{0}}}$
Ainsi ${\vec {v_{0}}}={\begin{pmatrix}{C_{x}}_{v}\\{C_{y}}_{v}\\{C_{z}}_{v}\end{pmatrix}}={\begin{pmatrix}\cos(\alpha )v_{0}\\\sin(\alpha )v_{0}\\\cos(\beta )v_{0}\end{pmatrix}}$ et donc
    ${\vec {v}}={\begin{pmatrix}\cos(\alpha )v_{0}\\-gt+\sin(\alpha )v_{0}\\\cos(\beta )v_{0}\end{pmatrix}}$
Nous avons donc établi le vecteur vitesse du projectile en fonction du temps. Or, par définition, la dérivée de la position en fonction du temps est la vitesse. Ainsi, en considérant que le projectile occupe le point P de ${\mathcal {R}}$, on définit les coordonnées horaires du projectile en déterminant le vecteur ${\overrightarrow {OP}}$ dont les composantes sont les primitives des composantes du vecteur ${\vec {v}}$.
On a alors
    ${\overrightarrow {OP}}={\begin{pmatrix}\cos(\alpha )v_{0}t+C_{x_{OP}}\\-{\frac {1}{2}}gt^{2}+\sin(\alpha )v_{0}t+C_{y_{OP}}\\\cos(\beta )v_{0}t+C_{z_{OP}}\end{pmatrix}}$ où $C_{x_{OP}}$, $C_{y_{OP}}$ et $C_{z_{OP}}$ sont des constantes.
Or à t=0, on a ${\overrightarrow {OP}}={\begin{pmatrix}C_{x_{OP}}\\C_{y_{OP}}\\C_{z_{OP}}\end{pmatrix}}$ et on connaît la position du projectile. En effet, il part en x=z=0 d'une hauteur h, donc $C_{x_{OP}}=0$, $C_{z_{OP}}=0$ et $C_{y_{OP}}=h$.
D'où
    $x(t)=\cos(\alpha )v_{0}t$
    $y(t)=-{\frac {1}{2}}gt^{2}+\sin(\alpha )v_{0}t+h$
    $z(t)=\cos(\beta )v_{0}t$.
On peut maintenant déterminer l'équation cartésienne du projectile, c'est-à-dire sa hauteur en fonction de la distance qu'il a parcouru.
On sait que
    $x=\cos(\alpha )v_{0}t$
donc, en isolant t, on a
    $t={\frac {x}{\cos(\alpha )v_{0}}}$.
Dès lors, en remplaçant cette expression de t dans l'expression de la composante y et z de l'équation horaire, on obtient l'équation cartésienne du projectile :
    $y(x)=-{\frac {1}{2}}{\frac {gx^{2}}{\cos ^{2}(\alpha )v_{0}^{2}}}+\tan(\alpha )x+h$
    $z(x)=x{\frac {\cos(\beta )}{\cos(\alpha )}}$.

### Distance parcourue

La distance horizontale totale d parcourue jusqu'à ce que le projectile retombe au sol est donnée par :
    $d={\frac {v_{0}}{g}}\cos \alpha \left(v_{0}\sin \alpha +{\sqrt {(v_{0}\sin \alpha )^{2}+2gh}}\right).$
Lorsque la hauteur initiale est égale à zéro, la distance parcourue se calcule ainsi :
    $d={\frac {v_{0}^{2}}{g}}\sin \left(2\alpha \right).$
On cherche ici à résoudre l'équation $y(x)=0$.
Or on sait que
    $y(x)=-{\frac {1}{2}}{\frac {g}{\cos ^{2}(\alpha )v_{0}^{2}}}x^{2}+\tan(\alpha )x+h$.
On peut donc résoudre l'équation. On remarque que l'expression de $y(x)$ est une fonction du second degré en x. On peut donc appliquer une résolution d'équation du second degré.
Donc le discriminant de $y(x)$ est :
    $\Delta =b^{2}-4ac=\tan ^{2}(\alpha )-4\left(-{\frac {1}{2}}{\frac {g}{\cos ^{2}(\alpha )v_{0}^{2}}}\right)h=\tan ^{2}(\alpha )+{\frac {2gh}{\cos ^{2}(\alpha )v_{0}^{2}}}>0$.
Donc $y(x)=0$ a deux solutions :
    ${\begin{aligned}d_{1}&={\frac {-b-{\sqrt {\Delta }}}{2a}}={\frac {-\tan(\alpha )-{\sqrt {\tan ^{2}(\alpha )+{\frac {2gh}{\cos ^{2}(\alpha )v_{0}^{2}}}}}}{-{\frac {g}{\cos ^{2}(\alpha )v_{0}^{2}}}}}\\&={\frac {\tan(\alpha )+{\sqrt {\tan ^{2}(\alpha )+{\frac {2gh}{\cos ^{2}(\alpha )v_{0}^{2}}}}}}{\frac {g}{\cos ^{2}(\alpha )v_{0}^{2}}}}={\frac {\cos ^{2}(\alpha )v_{0}^{2}}{g}}\left(\tan(\alpha )+{\sqrt {\tan ^{2}(\alpha )+{\frac {2gh}{\cos ^{2}(\alpha )v_{0}^{2}}}}}\right)\\&={\frac {\cos(\alpha )v_{0}}{g}}\left(v_{0}\sin(\alpha )+{\sqrt {(v_{0}\sin(\alpha ))^{2}+2gh}}\right)\end{aligned}}$
    $d_{2}={\frac {-b+{\sqrt {\Delta }}}{2a}}={\frac {\cos(\alpha )v_{0}}{g}}\left(v_{0}\sin(\alpha )-{\sqrt {(v_{0}\sin(\alpha ))^{2}+2gh}}\right)$
On ne retiendra que la solution d$_1$, qui est positive (car a est négatif).

### Temps de vol

Le temps de vol t est le temps qu'il faut au projectile pour retomber au sol, soit le temps où on atteint $x(t)=d$:
    $t={\frac {d}{v\cos \theta }}={\frac {v\sin \theta +{\sqrt {(v\sin \theta )^{2}+2gy_{0}}}}{g}}.$

### Hauteur maximale

La hauteur maximale h atteint par le projectile lors de son vol se calcule comme suit :
    $h={\frac {v_{0}^{2}\sin ^{2}(\theta )}{2g}}$.
On cherche ici le maximum d'une fonction, ce qui s'obtient en cherchant les racines de la dérivée de la hauteur, soit résoudre l'équation $y'(x)=0$. Or :
    $y'(x)=-{\frac {g}{\cos ^{2}(\theta )v_{0}^{2}}}x+\tan(\theta ).$
La résolution de cette équation linéaire donne donc :
    $y'(x)=-{\frac {g}{\cos ^{2}(\theta )v_{0}^{2}}}x+\tan(\theta )=0\ \Longleftrightarrow \ x={\frac {\tan(\theta )}{\frac {g}{\cos ^{2}(\theta )v_{0}^{2}}}}={\frac {\cos ^{2}(\theta )v_{0}^{2}}{g}}\tan(\theta )={\frac {v_{0}^{2}\cos(\theta )\sin(\theta )}{g}}={\frac {v_{0}^{2}\sin(2\theta )}{2g}}.$

### Angle de portée pour atteindre une distance d

L'angle de portée est l'angle θ avec lequel un projectile doit être lancé dans le but de parcourir une distance d, compte tenu de la vitesse initiale v_0. La trajectoire étant parabolique, il est plus simple d'utiliser les propriétés de symétrie de la trajectoire et donc que la hauteur maximale est atteinte à la moitié du parcours, donc
    ${\frac {v_{0}^{2}\sin(2\theta )}{2g}}={\frac {d}{2}}$
    $\sin(2\theta )={\frac {gd}{v_{0}^{2}}}$
donc :
    $\theta ={\frac {1}{2}}\arcsin \left({\frac {gd}{v_{0}^{2}}}\right)$
ou
    $\theta ={\frac {\pi }{2}}-{\frac {1}{2}}\arcsin \left({\frac {gd}{v^{2}}}\right)$

### Vitesse à une distance x

La valeur absolue |v| de la vitesse du projectile à une distance x est donnée par :
    $|v|={\sqrt {v^{2}-2gx\tan \theta +\left({\frac {gx}{v\cos \theta }}\right)^{2}}}.$

### Angle θ nécessaire pour atteindre les coordonnées (x ; y)

Pour atteindre une cible à distance x et une altitude y lorsque le projectile est lancé du point de coordonnées (0 ; 0) à une vitesse v, les valeurs de l'angle de portée θ se calculent ainsi :
    $\theta =\arctan {\left({\frac {v^{2}\pm {\sqrt {v^{4}-g(gx^{2}+2yv^{2})}}}{gx}}\right)}.$

## Trajectoire d'un projectile évoluant dans un fluide

Lorsqu'un projectile évolue dans un fluide, il subit une force de traînée qui s'oppose à son mouvement. Cette force de traînée est caractéristique de la forme du projectile, de sa vitesse, et des caractéristiques du fluide.
La traînée est un objet d'étude de la mécanique des fluides, qui ne considère plus le mouvement du mobile dans le fluide, mais qui suppose le mobile fixe et s'intéresse alors à l'écoulement du fluide autour de celui-ci. Cet écoulement est souvent complexe, et est caractérisé par un nombre sans dimension, appelé nombre de Reynolds et noté $\mathrm {R_{e}}$ :
    $\mathrm {R_{e}} ={\frac {v\,L\,\rho }{\eta }}$ avec v la vitesse relative du fluide en mètre par seconde, L la taille caractéristique de l’écoulement en mètres, ρ la masse volumique du fluide en kilogrammes par mètre cube, et η la viscosité dynamique du fluide en pascal-seconde.
On sait alors par l'expérience que :
* lorsque $\mathrm {R_{e}} <\!\!<1$ (il s'agit d'un écoulement de Stokes), la traînée est proportionnelle à la vitesse du mobile : le frottement que subit le projectile est dit linéaire ;
* lorsque $30<\mathrm {R_{e}} <800$ (l'écoulement est laminaire), la traînée est approximativement proportionnelle à la vitesse du mobile à la puissance 1,4 ;
* lorsque $\mathrm {R_{e}} >2000$ (l'écoulement est turbulent), la traînée est proportionnelle au carré de la vitesse du mobile : le frottement que subit le projectile est dit quadratique.
On modélisera la traînée comme étant proportionnelle à la vitesse du projectile élevée à une puissance λ :
    $F_{T}=-k\ v^{\lambda }\quad$ avec F$_T$ la force traînée du fluide en newtons, v la vitesse du projectile, et k le coefficient de proportionnalité caractéristique de la traînée qui s'exprime en kilogrammes par seconde
en prenant en compte la poussée d'Archimède, qui pousse le projectile dans le sens contraire au champ de pesanteur, d'autant plus que son volume est grand et que la masse volumique du fluide et l'intensité de son champ de pesanteur sont importants :
    ${\vec {P_{A}}}=-\rho \,V\,{\vec {g}}\quad$ avec ρ la masse volumique du fluide en kilogrammes par mètre cube et V le volume du projectile en mètres cubes.

### Approximation par analyse numérique

On suppose dans cette section que l’évolution de la vitesse se fait par approximation sur la tangente au point initial (méthode d'Euler). En procédant pas à pas, avec un pas très petit, cette solution est extrêmement satisfaisante, d’autant plus qu’il est très facile de la mettre en œuvre avec un ordinateur.
Avec γ le facteur de précision, c'est-à-dire le facteur qui définit la durée d'une unité de temps, ainsi qu'avec v$_i$ la vitesse initiale (en m/s) du projectile, son angle de portée α (en radians) et sa hauteur initiale h (en mètres), on a :
    $x\left(t+1\right)=x\left(t\right)+v_{x}\left(t\right)\gamma \quad {\text{ avec }}\quad x\left(0\right)=0$
    $y\left(t+1\right)=y\left(t\right)+v_{y}\left(t\right)\gamma \quad {\text{ avec }}\quad y\left(0\right)=h$
    $v_{x}\left(t\right)=v_{x}\left(t-1\right)+a_{x}\left(t-1\right)\gamma \quad {\text{ avec }}\quad v_{x}\left(0\right)=v_{i}\cos \left(\alpha \right)$
    $v_{y}\left(t\right)=v_{y}\left(t-1\right)+a_{y}\left(t-1\right)\gamma \quad {\text{ avec }}\quad v_{y}\left(0\right)=v_{i}\sin \left(\alpha \right)$
    $v\left(t\right)={\sqrt {v_{x}\left(t\right)^{2}+v_{y}\left(t\right)^{2}}}$
    $a_{x}\left(t\right)=-{\frac {1}{m}}k\ v\left(t\right)^{\lambda -1}\ v_{x}\left(t\right)$
    $a_{y}\left(t\right)=-g-{\frac {1}{m}}k\ v\left(t\right)^{\lambda -1}\ v_{y}\left(t\right)+{\frac {1}{m}}\rho \ V\ g$
où λ désigne la puissance caractéristique de la résistance de fluide à laquelle est élevée la vitesse du projectile, et V le volume du projectile rempli du fluide.
On définit $F_{T}$ comme la traînée que subit le projectile
    $F_{T}=-k\ v^{\lambda }$.
Alors ${\vec {F_{T}}}$, le vecteur représentant la traînée, s'écrit
    ${\vec {F_{T}}}=-k\ v^{\lambda -1}\ {\vec {v}}$.
On définit ${\vec {P_{A}}}$ comme le vecteur représentant la poussée d’Archimède qui s’exerce sur le projectile. On a donc
    ${\vec {P_{A}}}=-\rho \ V\ {\vec {g}}$.
On définit ${\vec {P}}$ comme le poids du projectile. On a donc
    ${\vec {P}}=m{\vec {g}}$.
On se place dans un référentiel galiléen. On suppose que la masse du projectile reste constante durant le vol.
Dès lors, d’après la deuxième loi de Newton, on peut écrire que
    $\sum _{i}{\vec {F}}_{i}=m{\vec {a}}={\vec {P}}+{\vec {F_{T}}}+{\vec {P_{A}}}$
soit
    ${\vec {a}}={\vec {g}}+{\frac {\vec {F_{T}}}{m}}+{\frac {\vec {P_{A}}}{m}}$.
Par la méthode d’Euler, nous avons
    ${\vec {v}}\left(t+\Delta t\right)\approx {\vec {v}}\left(t\right)+{\vec {a}}\left(t\right)\Delta t$.
Puis, avec O l’origine du repère et P le point du plan occupé par le projectile, de la même manière la méthode d’Euler dispose que
    ${\vec {OP}}\left(t+\Delta t\right)\approx {\vec {OP}}\left(t\right)+{\vec {v}}\left(t\right)\Delta t$.
On peut alors définir les suites $\left(a_{x_{n}}\right)$ et $\left(a_{y_{n}}\right)$ pour tout entier naturel $n$ par
    $a_{x_{n}}=-{\frac {1}{m}}\ k\ v^{\lambda -1}\ v_{x_{n}}$
    $a_{y_{n}}=-g-{\frac {1}{m}}\ k\ v^{\lambda -1}\ v_{y_{n}}+{\frac {1}{m}}\ \rho \ V\ g$
et les suites $\left(v_{x_{n}}\right)$ et $\left(v_{y_{n}}\right)$ pour tout entier naturel $n$ strictement positif par
    $\left(v_{x_{n}}\right):\left\{{\begin{array}{ll}v_{x_{0}}=v_{i}\cos \left(\alpha \right)\\v_{x_{n}}=v_{x_{n-1}}+a_{x_{n-1}}t\end{array}}\right.$
    $\left(v_{y_{n}}\right):\left\{{\begin{array}{ll}v_{y_{0}}=v_{i}\sin \left(\alpha \right)\\v_{y_{n}}=v_{y_{n-1}}+a_{y_{n-1}}t\end{array}}\right.$
où $v_{i}$ désigne la vitesse initiale et où le réel $t$ est le facteur de précision.
Enfin, on définit les suites $\left(x_{n}\right)$ et $\left(y_{n}\right)$ pour tout entier naturel $n$ par
    $\left(x_{n}\right):\left\{{\begin{array}{ll}x_{0}=0\\x_{n+1}=x_{n}+v_{x_{n}}t\end{array}}\right.$
    $\left(y_{n}\right):\left\{{\begin{array}{ll}y_{0}=h\\y_{n+1}=y_{n}+v_{y_{n}}t\end{array}}\right.$

### Résolution analytique dans le cas quadratique

Lorsque la traînée est proportionnelle au carré de la vitesse du projectile, on fait face à un système d'équations différentielles non-linéaires couplées, et la résolution nécessite une approche plus engagée. Newton lui-même était incapable de résoudre le problème, mais son contemporain Johann Bernoulli réussit après avoir été mis au défi par l'astronome britannique John Keill. Sa solution est connue sous le nom d'équation de l'hodographe, et bien qu'elle soit la plus citée de la littérature physique, elle contient des termes qui ne peuvent être évalués que numériquement, c'est-à-dire que cette solution n'est pas une véritable solution analytique explicite.
Par véritable solution analytique, il faut comprendre une solution qui établisse les équations du mouvement pour toutes les conditions initiales, sans recourir à des approximations ou des simplifications qui ne sont valides que pour certains cas particuliers de conditions initiales. Une solution est explicite si elle dépend explicitement d'une variable indépendante naturelle, ici le temps, et non de plusieurs autres variables auxiliaires, comme le fait la solution de l'hodographe proposée par Johann Bernoulli.
Alors nous avons
    $v_{x}\left(t\right)={\frac {v_{x_{0}}}{\phi \left(t\right)}}\quad {\text{ avec }}\quad v_{x_{0}}=v_{0}\cos \left(\alpha \right)$
et
    $v_{y}\left(t\right)={\frac {v_{y_{0}}-g\ \Phi \left(t\right)}{\phi \left(t\right)}}\quad {\text{ avec }}\quad v_{y_{0}}=v_{0}\sin \left(\alpha \right)$
avec
    $\phi (t)=\exp \left[\omega \ \int _{0}^{t}{\sqrt {v_{x}^{2}\left(s\right)+v_{y}^{2}\left(s\right)}}\ \mathrm {d} s\right]\quad {\text{ avec }}\quad \omega =-{\frac {k}{m}}$
et
    $\Phi ^{\prime }\left(t\right)=\phi \left(t\right)$
de laquelle il existe un développement en série de Taylor qui s’écrit :
    $\Phi \left(t\right)=\sum _{j=0}^{+\infty }{\frac {\Phi ^{(j)}(0)}{j!}}\ t^{j}\quad$
avec $\Phi (0)=0$ et $\Phi ^{\prime }(0)=1$, et pour tout entier naturel j
    $\Phi ^{(j+2)}\left(0\right)=\omega \ {\sqrt {v_{x_{0}}^{2}+\left(v_{y_{0}}-g\ \Phi ^{(j)}\left(0\right)\right)^{2}}}$.

### Détermination de la traînée

On rappelle que $\mathrm {R_{e}} ={\frac {v\ L\ \rho }{\eta }}$ avec v la vitesse relative du fluide en mètre par seconde, L la taille caractéristique de l’écoulement en mètre, ρ la masse volumique du fluide en kilogramme par mètre cube, et η la viscosité dynamique du fluide en pascal-seconde.
Lorsque $\mathrm {R_{e}} <\!\!<1$, l'écoulement suit une loi de Stokes. En particulier, dans le cas d'une sphère, la traînée F$_T$ s'exprime
    $F_{T}=-3\ \pi \ \eta \ D\ v\quad$ avec η la viscosité dynamique du fluide, D le diamètre de la sphère, et v la vitesse du mobile.
Alors le coefficient de proportionnalité k s'écrit $k=-3\pi \,\eta \,D$.
Lorsque $30<\mathrm {R_{e}} <800$, l'écoulement suit une loi de Poiseuille.
Lorsque $\mathrm {R_{e}} >2000$, l'écoulement est turbulent et la traînée F$_T$ s'exprime
    $F_{T}={\frac {1}{2}}\ \rho \ S\ C_{x}\ v^{2}\quad$ avec S la surface de référence du projectile, C$_x$ le coefficient de traînée du projectile.
Alors le coefficient de proportionnalité $k$ s'écrit $k={\frac {1}{2}}\rho \,S\,C_{x}$.

#### Taille caractéristique de l'écoulement

Pour l'étude de la traînée des corps géométriques, lorsque l'objet est profilé, la taille de caractéristique de l'écoulement est la largeur de sa surface frontale. Lorsque l'objet est non-profilé, il s'agit de la longueur de sa surface parallèle à l'écoulement.

#### Masse volumique

La masse volumique de l’air est fonction de la température, de la masse molaire et de la pression. En choisissant pour pression celle de l'atmosphère standard internationale (ISA) au niveau de la mer : P$_0$ = 101 325 Pa = 1 013,25 mbar = 1 013,25 hPa, on peut considérer que
    $\rho _{\mathrm {air} }=1,292\times {\frac {273,15}{T}}\quad$ avec T en kelvins.
Plus généralement, on peut se référer aux tables des masses volumiques de diverses substances.

#### Viscosité dynamique

On peut utiliser la loi semi-empirique de Sutherland pour déterminer la viscosité dynamique dans l'air :
    $\eta _{\mathrm {air} }=B{\frac {\sqrt {T}}{1+{\frac {c}{T}}}}\quad$ avec $B=1,4810^{-6}\;\mathrm {Ns} \ \mathrm {K} ^{-0{,}5}\mathrm {m} ^{2}$ et $c=119,4\;\mathrm {K}$ des constantes déterminées empiriquement, et T la température de l’air en kelvins.

#### Surface de référence

Lorsque l'objet est profilé, la surface de référence du mobile est sa surface frontale, c'est-à-dire la surface occupée par l'objet sur la projection perpendiculaire à l'écoulement. Lorsque l'objet est non-profilé, la surface de référence du mobile est sa surface alaire.

#### Coefficient de traînée

Le coefficient de traînée $C_{x}$ est un nombre sans dimension dont la valeur peut être déterminée expérimentalement en soufflerie.
Bien que le coefficient de traînée soit fonction du nombre de Reynolds, on peut tout de même utiliser en première approximation les coefficients de traînée de l'image ci-contre.
On peut aussi, dans le cas d'une sphère, se référer à ce graphique qui donne le coefficient de traînée en fonction du nombre de Reynolds. Dans certaines plages de ce nombre de Reynolds, les ingénieurs considèrent souvent le $C_{x}$ de la sphère comme constant.