分子立体碰撞模拟.docVIP

分子立体碰撞模拟 作者 hank 文件集合： 1. ball.cpp 2. ball.h 3. Collide_Cube.cpp 4. StdAfx.cpp 5. StdAfx.h 6．vector3f.cpp 7. vector3f.h 效果图： 附录文件： Collide_Cube.cpp： // Collide_Cube.cpp : Defines the entry point for the console application. // // 立体图形.cpp : Defines the entry point for the console application. // #include stdafx.h #include ball.h #include vector3f.h #define M *1 #define DM *0.1 #define CM *0.01 #define MM *0.001 #define L 1.0 M //立方体长度 #define L0 (L*0.5) //#pragma comment(linker, /subsystem:\windows\ /entry:\mainCRTStartup\) //小球的信息 const int n = 5; //每个维度上分布小球的个数，总球数为 n*n*n #define N n*n*n double r = 10 MM; //小球的半径 double speed = 1 MM; //初始速度 ball MyBall[N]; vector3f loc,vel,clr; double m(1.0); double x,y,z; double wall[6]; //小球到立方体壁的距离 vector3f wVtr[6]; //六个壁面的垂直向量 //用于鼠标控制图形变换 float tx(0.0),ty(0.0),tz(0.0); float theta(0),phi(0); int mouse_x,mouse_y; bool mouse_button_pressed = false; //数据存储以及菜单 ofstream fout(data.txt); bool IsWrite = 0; //判断是否写入数据 bool IsW2 = 0;//判断多少毫秒写入数据 int xx; //每多少帧读入一次数据 float yy = 3000; bool Reset = 0; int SpeedState = 0; int w2f,menu; //菜单 int i,j,k; void Init(); void GLScene(); void Reshape(int w,int h); float rndnp(); float rndpi(); void mytimer(int value); void mymenu(int value); void mysubmenu(int value); void motion(int x,int y); void mouse(int button,int state,int x,int y); int main(int argc, char* argv[]) { glutInit(argc,argv); glutInitDisplayMode(GLUT_SINGLE|GLUT_RGB); glutInitWindowSize(600,600); glutInitWindowPosition(100,100); glutCreateWindow(立体图形-Hank); Init(); //副菜单 w2f = glutCreateMenu(mysubmenu); glutAddMenuEntry(Begin,0); glutAddMenuEntry(Stop,1); glutAddMenuEntry(Return,3); //主菜单 menu = glutCreateMenu(mymenu); glutAddMenuEntry(Reset,0); glutAddMenuEntry(Faster,2); glutAddMenuEntry(Slower,3); glutAddMenuEntry(Write2File,4); glutAddMenuEntry(Exit,1); glutSetMenu(menu); glutAttachMenu(GLUT_RIGHT_BUTTON); glutTimerFunc(yy,mytimer,1); glutDisplayFunc(GLScene); glutIdleFunc(GLScene); glutReshapeFu