Calcul numérique en Physique en C (mouvement d'un obus)
Description du problème
Soit un obus soumis à son poids et à des forces de frottements (proportionelles à la vitesse ou au carré de la vitesse).
Nous cherchons à décrire la trajectoire de cet obus.
Réponse
Téléchargez le code source∞ C répondant à ce problème.
On compile en tapant :
g main.cc -o main
On exécute en tapant :
./main
Afin de visualiser le résultat (dans
main.dat∞),
GNUPlot est nécessaire pour interpréter le fichier
main.gp∞. Vous pouvez aussi télécharger le
Makefile∞ du projet.
Remarques
On peut observer la trajectoire de l'obus dans les trois cas suivant :
- Cas 1 : absence de frottements
sans_frott.png∞
- Cas 2 : frottements proportionnels à la vitesse
frott_visq.png∞
- Cas 3 : frottements proportionnels au carré de la vitesse
frott_fluide.png∞
On remarque que :
- Cas 1 : la trajectoire est bien parabolique.
- Cas 2 : la trajectoire présente une asymptote oblique caractéristique d'une force de frottement proportionnelle à la vitesse.
- Cas 3 : la trajectoire présente une asymptote verticale caractéristique d'une force de frottement proportionnelle au carré de la vitesse.
Obtenir le code source
Accès web
http://svn.berlios.de/viewcvs/openphysic/cpp/physique/meca/rfd_obus Accès anonyme
svn checkout svn://svn.berlios.de/openphysic/cpp/physique/meca/rfd_obus Accès développeur
svn checkout svn+ssh://scls19fr@svn.berlios.de/svnroot/repos/openphysic/cpp/physique/meca/rfd_obus
There are no comments on this page. [Add comment]