机械原理课程设计编程说明书

1、机械原理课程设计编写说明书设计题目：牛头刨床机构的设计及运动分析（1） 指导教师：王春华 何凡设计者：冯帅学号：01 班级：机械07-2班2009年7月10日辽宁工程技术大学目 录1设计任务及要求2数学模型的建立3程序框图4程序清单及运行结果5设计总结6参考文献一、设计任务及要求已知：曲柄每分钟转数n2，各构件尺寸及重心位置，且刨头导路X-X位于导杆端点B所作圆弧的平分线上， 数据见下表：设计内容导 杆 机 构 的 运 动 分 析单位r/minmm符号n2lo2o4lo2Alo4BlBClo4s4参数603801105400.25lo4B 0.5lo4B要求：（1）作机构的运动简图（2）用C语

2、言编写主程序调用子程序，对机构进行运动分析，动态显示，并打印程序及运算结果。（3）画出导轨4的角位移，角速度，角加速度。二、 数学模型如图四个向量组成封闭四边形，于是有 (1)按复数式可以写成由于3=90，上式可化简为 (2)根据（2）式中实部、虚部分别相等得 (3) (4)(3) (4)联立解得 (5)b=a2+d2+2adsin (6)将（2）对时间求一阶导数得 (7) (8)将（2）对时间求二阶导数得 (9) (10)即滑块沿杆方向的加速度，通常曲柄可近似看作均角速转动，则。三、程序框图输入a,d,输入N作循环，For(i=0;i=N;i+)依次计算,b,ac将计算的,b,ac结果分别存

3、入数组或文件中按格式输出所有计算结果初始化图形系统绘制直角坐标系直角坐标系下分别绘出角位移、角速度、角加速度图象四、程序清单及运算结果(1)数据#includemath.h#includeconio.h#include#define PI 3.#define M 0. /* M=PI/180 */main() int P1,i=0,j=0,P_171;float P_371,P_3171,P_31171;float P3,P31,P311,w3;float s3,s31,s311;FILE *f1;if(f1=fopen(fdata.txt,w)=NULL)printf(fdata.txt f

4、ile cannot open!n);exit(0);clrscr();printf(jiao du wei yi su du jiao jia su du n);for(P1=0;P1=0&P1270&P190&P1270)P3=PI+atan(380+110*sin(P1*M)/(110*cos(P1*M);P3/=M;if(P3!=90&P3!=270)s3=110*cos(P1*M)/cos(P3*M);else if(P3=90)s3=440;elses3=240;s31=-110*6.3*sin(P1-P3)*M);P31=110*6.3*cos(P1-P3)*M)/s3;w3=P

5、31;P311=(6.3*6.3*110*sin(P3-P1)*M)-2*w3*s31;P_1j=P1;P_3j=P3;P_31j=P31;P_311j=P311;printf(%d %6.4f %6.4f %6.4fn,P_1j,P_3j,P_31j,P_311j);fprintf(f1,%d %6.4f %6.4f %6.4fn,P_1j,P_3j,P_31j,P_311j);j=j+1;fclose(f1);程序运行结果Jiao du wei yi su du jiao jia su du0 73.8569 0.4871 3545.24105 74.2912 0.6056 3298.72

6、3610 74.8152 0.7131 3056.405315 75.4205 0.8105 2819.880620 76.0993 0.8984 2590.126025 76.8443 0.9775 2367.655030 77.6487 1.0484 2152.623035 78.5064 1.1116 1944.927140 79.4114 1.1677 1744.271745 80.3582 1.2172 1550.224150 81.3417 1.2603 1362.259460 83.4000 1.3293 1002.154565 84.4658 1.3557 828.710470

7、 85.5506 1.3771 658.767075 86.6503 1.3934 491.628080 87.7612 1.4051 326.597885 88.8793 1.4120 162.972490 90.0000 1.5750 0.000095 91.1230 1.4120 -162.8737100 92.2412 1.4051 -326.4991105 93.3520 1.3934 -491.5290110 94.4518 1.3771 -658.6653Jiao du wei yi su du jiao jia su du115 95.5365 1.3558 -828.6074

8、120 96.6024 1.3294 -1002.0498125 97.6452 1.2977 -1179.6742130 98.6607 1.2604 -1362.1470135 99.6443 1.2172 -1550.1101140 100.5911 1.1678 -1744.1528145 101.4961 1.1117 -1944.8036150 102.3538 1.0484 -2152.4961155 103.1583 0.9776 -2367.5242160 103.9033 0.8985 -2589.9907165 104.5821 0.8106 -2819.7402170

9、105.1875 0.7132 -3056.2615180 106.1459 0.4872 -3545.0913185 106.4817 0.3571 -3793.4160190 106.7095 0.2149 -4040.1492195 106.8195 0.0600 -4280.6279200 106.8013 -0.1081 -4508.6533205 106.6445 -0.2892 -4716.2061210 106.3388 -0.4832 -4893.1987215 105.8743 -0.6891 -5027.3159220 105.2423 -0.9052 -5104.034

10、7225 104.4354 -1.1291 -5106.9404230 103.4489 -1.3572 -5018.4829Jiao du wei yi su du jiao jia su du235 102.2813 -1.5845 -4821.2866240 100.9355 -1.8052 -4500.0879245 99.4196 -2.0119 -4044.2471250 97.7478 -2.1967 -3450.5139255 95.9407 -2.3516 2725.5784260 94.0251 -2.4687 -1887.6504265 92.0338 -2.5418 -

11、966.4163270 90.0000 -1.5750 -0.3448275 87.9729 -2.5419 964.6959280 85.9814 -2.4689 1886.0481285 84.0657 -2.3518 2724.1545290 82.2584 -2.1970 3449.3188265 92.0338 -2.5418 -966.4163270 90.0000 -1.5750 -0.3448275 87.9729 -2.5419 964.6959280 85.9814 -2.4689 1886.0481285 84.0657 -2.3518 2724.1545290 82.2

12、584 -2.1970 3449.3188Jiao du wei yi sudu jiao jia su du295 80.5863 -2.0122 4043.2969300 79.0701 -1.8056 4499.3887305 77.7240 -1.5850 4820.8208310 76.5561 -1.3576 5018.2305315 75.5693 -1.1296 5106.8682320 74.7621 -0.9056 5104.1118325 74.1297 -0.6895 5027.5117330 73.6649 -0.4836 4893.4844335 73.3590 -

13、0.2896 4716.5571340 73.2019 -0.1084 4509.0503345 73.1834 0.0597 4281.0542350 73.2932 0.2147 4040.5911Jiao du wei yi su du jiao jia su du355 73.5208 0.3569 3793.865(2)导杆机构运动模拟程序#include#include#includevoid pist(float x0,float y0,float l,float h,float theta);void pirot(float x0,float y0,float l);void

14、sgd();main() int gd=DETECT,gmode,n; initgraph(&gd,&gmode,); for(n=1;n=60;n+) setbkcolor(0); sgd(); getch(); /* 摇块子程序 */ void pist(float x0,float y0,float l,float h,float theta) float x,y; x=x0-cos(theta)*1/2+h/2*sin(theta); y=y0+1/2*sin(theta)+h/2*cos(theta); moveto(x,y); linerel(l*cos(theta),-l*sin

15、(theta); linerel(-h*sin(theta),-h*cos(theta); linerel(-l*cos(theta),l*sin(theta); lineto(x,y); /* 支点子程序 */ void pirot(float x0,float y0,float l) float pi=3.; float x,y; int i,n; int h=4; setcolor(3); circle(x0,y0,3); setcolor(150); x=x0-1/2; y=y0+sin(pi/3)*l; moveto(x0,y0); lineto(x,y); linerel(20,0

16、); lineto(x0,y0); n=1/4; for(i=0;i=n;i+) moveto(x+i*h,y+h); lineto(x+(i+1)*h,y); /* 连杆，摇块，导轨程序*/ void sgd() int i; int x04,y04,x02,y02; int l04b=580; int l02a=110; float xa,ya,xb,yb,l; float ss,theta; for(i=0;i=59;i+) /*循环动画*/ ss=i*0.1; x04=380; x02=380; /*坐标计算*/ y04=320; y02=135; xa=x02+l02a*cos(ss

17、); ya=y02+l02a*sin(ss); l=sqrt(xa-x04)*(xa-x04)+(ya-y04)*(ya-y04); xb=x04+(xa-x04)*l04b/l; yb=y04+(ya-y04)*l04b/l; theta=-atan(ya-y04)/(xa-x04); setcolor(14); cleardevice(); /*清屏*/ circle(xa+4,ya-8,3); setcolor(2); line(x02,y02,xa+4,ya-8); line(xb,yb,x04,y04); setcolor(9); pirot(x02,y02,16); /*画支点Xo

18、4,Xo2*/ pirot(x04,y04,16); pist(xa,ya,20,10,theta); delay(3000); 程序运行结果图如下：(3)导杆3的角位移的，角速度，角加速度曲线图的设计程序#include#include#include#define pi 3.main() float w1=2*pi*60/60,M1=110,M6=380;float P3,P1,s3,w3,s31,s311;int gd=DETECT,gmode,i;initgraph(&gd,&gmode,);clrscr();for(P1=0;P1=0&P1pi*1.5&P1=2*pi)P3=atan

19、(M6+M1*sin(P1)/(M1*cos(P1);else if(P1=pi/2|P1=1.5*pi)P3=pi/2;else P3=pi+atan(M6+M1*sin(P1)/(M1*cos(P1);if(P1!=pi/2&P1!=1.5*pi) s3=M1+M6; else s3=M6-M1;s31=-M1*w1*sin(P1-P3);w3=M1*w1*cos(P1-P3)/s3;s3=sqrt(M1*M1+M6*M6+2*M1*M6*sin(P1);putpixel(100+P1*180/pi,60+s3*sin(P3)/2.5,5);/* 位移 */line(100,200,500,200); /* 坐标系 */line(495,205,500,200);line(495,195,500,200);line(100,10,100,3