Bonjour,
Dans le cadre de mon master astro, je dois composer un programme en C simulant la dynamique moléculaire des particules soumises a un potentiel de lennard-Jones.
ne sachant pas comment procéder et tournant en rond depuis un moment, je sollicite votre aide.
Merci d'avance si vous pouvez m'aider
Phlegyas
Salut,
La manière la plus simple est de simplement stocker la position de tes particules sous forme de tableau, puis de faire des boucles pour chaque pas de temps où tu calcules pour chaque particule la force qu'elle subit des autres particules et donc son déplacement. Tu mets à jour toutes les positions et c'est reparti pour un tour.
Personnellement j'ai une préférence pour les structures, ça permet d'être un peut plus claire.
Tu peux par exemple faire un truc du genre :
struct Structure_molecule
{
double x, y, z;
double force;
//et tout autre parametres utile
};
Ensuite dans ta fonction principale rien ne t'interdit de crée un tableau de structure
struct Structure_molecule Mes_molecules[1000];
Tu pourra alors récupérer chaque position de la manière suivante :
Mes_molecules[0].x ... ;
Mes_molecules[0].y ... ;
Mes_molecules[0].z ... ;
Mes_molecules[1].x ... ;
Mes_molecules[1].y ... ;
Mes_molecules[1].z ... ;
...
ou alors simplement faire une boucle :
int i = 0;
for(i = 0 ; i < 1000 ; i++)
{
Mes_molecules[i].x ... ;
Mes_molecules[i].y ... ;
Mes_molecules[i].z ... ;
}
De la même manière tu peux récupérer la force appliqué sur chaque molécule, ou tout autre paramètre que tu aura déclaré dans ta structure.
Ok je te suis coincoin... mais pourrais tu me donner un ou 2 tuyau de comment procéder car je suis pas tellement doué en programmation malheureusement....
Merci pour l'explication Jackyzgood... mais comme je l'ai dit, je suis mauvais en programmation donc une explication pas a pas me serais très utile ( sans pour autant me mâcher le boulot bien sur, le but étant de comprendre ce que je fait)
En ce moment j'essaye de crée un petit moteur physique
http://www.youtube.com/watch?v=1CNn3VX9RkA
Il y a encore quelques bugs, mais la base est là. Je ne sais pas exactement ce qu'est le potentiel de lennard-Jones, mais je pense qu'on peut reprendre une bonne partie de mon code source pour faire ce que tu souhaite.
Contacte moi par mp, ou sur msn pour me dire ou je peux t'envoyer les fichiers
Tu peux essayer avec seulement deux particules pour ne pas avoir toutes les difficultés d'un coup.
