东北大学机械原理课程设计(牛头刨床)

1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业机械原理课程设计说明书题目： 牛头刨床机构方案分析 班 级 ：姓 名 ：学 号 ：指导教师 ：成 绩 ：2012年 月 日目录 TOC o 1-2 h z u 1题目说明（图a）所示为牛头刨床主机构的运动简图，已知，l1=0.1m,l0=0.4m,l3=0.75m,l4=0.15m,ly=0.738m,l3=0.375m,a=0.05m,b=0.15,c=0.4m,d=0.1m。只计构件3、5的质量，其余略去不计，m3=30kg,Js3=0.7kgm2,m5=95kg。

2、工艺阻力Q如图所示，Q=9000N。主轴1的转速为60r/min(顺时针方向)，许用运转不均匀系数=0.03。2内容要求与作法1.进行运动分析，绘制滑枕5的位移、速度和加速度线图。2.进行动态静力分析，绘制固定铰链处的反力矢端图金额平衡力矩Tb的线图。3.以主轴为等效构件，计算飞轮的转动惯量。上机前认真读懂所用子程序，自编主程序，初始位置取滑枕5的左极限位置。3.滑枕初始位置及行程H的确定方法滑枕初始位置为左极限位置,此时构件1和构件3在O1O2左侧垂直；滑枕行程H由运动分析的结果得知。滑枕的左极限位置为-0.337，右极限为0.038，所以行程为0.375；4机构的结构分析213213 =

3、2 * GB3 84 = 3 * GB3 7 = 4 * GB3 主动件RPR杆组RRP杆组36955-1杆组运动分析子程序函数及形参与真值对照表1)调用bark函数，求3点的运动参数。形式参数n1n2n3kr1r2gamtwepvpap实 值1301r130.00.0twepvpap 2)调用rprk函数，求构件3的运动参数。形式参数mn1n2k1k2r1r2vr2ar2twepvpap实 值123320.0&r2&vr2&ar2twepvpap3)调用bark函数，求4点的运动参数。形式参数n1n2n3kr1r2gamtwepvpap实 值2403r240.00.0twepvpap4)调用

4、rrpk函数，求5点的位置和5构件的运动参数。形式参数mn1n2n3k1k2k3r1r2vr2ar2twepvpap实 值-1495456r45&r2&vr2&ar2twepvpap 5-2主程序、运行结果、运行图像主程序/* Note:Your choice is C IDE */#include ”graphics.h” /*图形库*/#include subk.c /*运动分析子程序*/#include draw.c /*绘图子程序*/main() static double p202,vp202,ap202,del; static double t10,w10,e20,pdraw370

5、,vpdraw370,apdraw370; static int ic; /*定义静态变量*/ double r13,r24,r45,r12; double pi,dr; double r2,vr2,ar2; int i; /*定义局部变量*/ FILE *fp; /*定义文件指针变量*/char *m=p,vp,ap; /*定义字符串，改变字符串会改变绘图结果的标注*/ r13=0.1;r24=0.75;r45=0.15;r12=0.4; pi=4.0*atan(1.0); /*求*/ dr=pi/180.0 /*求弧度*/; t6=0.0;w6=0.0;e6=0.0; w1=-2*pi;e

6、1=0.0; del=15.0; p21=0.0; p22=0.0; p11=0.0; p12=0.4; p91=-100; p92=0.738; /*变量赋值*/ t6=0.0*dr; /*变角度为弧度*/ printf(n The Kinematic Parameters of Point5n); printf(No THETA1 S5 V5 A5n); printf(n deg m m/s m/s/sn); /*在屏幕上写表头*/ if(fp=fopen(file1,w)=NULL) printf(Cant open this file.n); exit(0); /*建立并打开文件fil

7、e1*/ fprintf(fp,n The Kinematic Parameters of Point5n); fprintf(fp,No THETA1 S5 V5 A5n); fprintf(fp, deg m m/s m/s/sn); /*在文件file1中写表头*/ ic=(int)(360.0/del); for(i=0;i=ic;i+) /*建立循环，调用运动分析子程序*/ t1=(-(pi-asin(r13/r12)+(-i)*del*dr; bark(1,3,0,1,r13,0.0,0.0,t,w,e,p,vp,ap); rprk(1,2,3,3,2,0.0,&r2,&vr2,&

8、ar2,t,w,e,p,vp,ap); bark(2,4,0,3,r24,0.0,0.0,t,w,e,p,vp,ap); rrpk(-1,4,9,5,4,5,6,r45,&r2,&vr2,&ar2,t,w,e,p,vp,ap); /*r2,vr2,ar2在子程序中已定义为指针变量，所以其前面要加&符号*/ printf(n%2d%12.3f%12.3f%12.3f%12.3f,i+1,t1/dr,p51,vp51,ap51); /*把运算结果写屏幕上*/ fprintf(fp,n%2d%12.3f%12.3f%12.3f%12.3f,i+1,t1/dr,p51,vp51,ap51); /*把运

9、算结果写在文件中*/ pdrawi=p51; vpdrawi=vp51; apdrawi=ap51; /*将运算结果传给pdrawi 、vpdrawi 、apdrawi以备绘图用*/ if(i%16)=0)getch(); /*屏幕满16行停顿*/ fclose(fp); /*关闭文件*/ getch(); draw1(del,pdraw,vpdraw,apdraw,ic,m); /*调用绘图程序*/ /*主程序结束*/运行结果The Kinematic Parameters of Point5No THETA1 S5 V5 A5 deg m m/s m/s/s1 -165.522 -0.33

10、7 -0.000 7.251 2 -180.522 -0.331 0.273 5.878 3 -195.522 -0.315 0.492 4.688 4 -210.522 -0.291 0.665 3.611 5 -225.522 -0.260 0.794 2.593 6 -240.522 -0.225 0.882 1.623 7 -255.522 -0.187 0.930 0.715 8 -270.522 -0.148 0.942 -0.123 9 -285.522 -0.109 0.921 -0.90310 -300.522 -0.072 0.867 -1.66811 -315.522 -

11、0.037 0.781 -2.49512 -330.522 -0.007 0.657 -3.47013 -345.522 0.017 0.489 -4.67114 -360.522 0.033 0.265 -6.10815 -375.522 0.038 -0.022 -7.66016 -390.522 0.030 -0.371 -9.00017 -405.522 0.006 -0.762 -9.57818 -420.522 -0.033 -1.148 -8.64619 -435.522 -0.088 -1.450 -5.45920 -450.522 -0.152 -1.571 -0.05421

12、 -465.522 -0.216 -1.450 5.69922 -480.522 -0.270 -1.128 9.21023 -495.522 -0.308 -0.723 9.82024 -510.522 -0.330 -0.333 8.74625 -525.522 -0.337 0.000 7.251 运行图像 图5-1图5-16.牛头刨床的静力分析6-1杆组运动分析子程序函数及形参实值对照表1)调用bark函数，求3点的运动参数。形式参数n1n2n3kr1r2gamtwepvpap实 值1301r130.00.0twepvpap2）调用rprk函数，求构件3的运动参数。形式参数mn1n2k

13、1k2r1r2vr2ar2twepvpap实 值123320.0&r2&vr2&ar2twepvpap3)调用bark函数，求4点的运动参数。形式参数n1n2n3kr1r2gamtwepvpap实 值24030.0r240.0twepvpap4)调用rrpk函数，求5点的位置和5构件的运动参数。形式参数mn1n2n3k1k2k3r1r2vr2ar2twepvpap实 值-1495456r45&r2&vr2&ar2twepvpap5)调用bark函数，求构件2的质心6点的运动参数。形式参数n1n2n3kr1r2gamtwepvpap实 值2603r260.00.0twepvpap6)调用bark

14、函数，求构件5的质心7点的运动参数。形式参数n1n2n3kr1r2gamtwepvpap实 值50750.0r57-161.57*drtwepvpap7)调用bark函数，求构件5的工艺阻力作用点8点的运动参数。形式参数n1n2n3kr1r2gamtwepvpap实 值50850.0r58165.96*drtwepvpap8) 调用rrpf函数，求4、5点的反作用力。形式参数n1n2n3ns1ns2nn1nn2nexfk1k2pvpaptwefr实 值4950708845pvpaptwefr9) 调用rprf函数求3点的反作用力。形式参数n1n2ns1ns2nn1nn2nexfk1k2pvpa

15、ptwefrfkpk实 值236040032pvpaptwefrfkpk10）调用barf函数，求1点运动副反力及平衡力矩。形式参数n1ns1nn1k1papefrtb实 值1031papefr&tb 6-2主程序、运行结果、运行图像主程序/* Note:Your choice is C IDE */#include ”graphics.h” /*图形库*/#include subk.c /*运动分析子程序*/ #include subf.c /*动态静力分析子程序*/#include draw.c /*绘图子程序*/main() static double p202,vp202,ap202,

16、del; static double t10,w10,e10,tbdraw370,tb1draw370; static double sita1370,fr1draw370,sita2370,fr2draw370,sita3370,fr3draw370; static double fr202,fe202,fk202,pk202, tb,tb1,fr1,bt1,fr2,bt2,we1,we2,we3,we4,we5; static int ic; double r12,r13,r24,r45,r26,r58,r57; double pi,dr; double r2,vr2,ar2; int i

17、; FILE *fp; char *m=tb,tb1,fr1,fr2;/*定义字符串，前两个字符串是为画平衡力矩曲线提供的；后三个字符串是为画反力矢端图提供的；改变字符串会改变绘图结果的标注*/ sm1=0.0;sm2=0.0;sm3=30.0;sm4=0.0;sm5=95.0; sj1=0.0;sj2=0.0;sj3=0.7;sj4=0.0; r13=0.1;r24=0.75;r45=0.15;r26=0.375; r58=sqrt(0.1*0.1+0.25*0.25); r57=sqrt(0.15*0.15+0.05*0.05);r12=0.4; pi=4.0*atan(1.0); dr=

18、pi/180.0; t6=0.0;w6=0.0;e6=0.0; w1=-2*pi;e1=0.0;del=15.0; t6=0.0*dr; p21=0.0; p22=0.0; p11=0.0; p12=0.4; p91=-100; p92=0.738; printf(n The Kineto-static Analysis of a six-bar Linkasen); printf(No HETAL fr1 sita1 fr2 sita2 tb tb1n); printf( deg N radian N radian N.m N.m ); if(fp=fopen(file1,w)=NULL)

19、printf(Cant open this file.n); exit(0); fprintf(fp,n The Kineto-static Analysis of a six-bar Linkasen); fprintf(fp,No HETAL fr1 sita1 fr2 sita2 tb tb1n); fprintf(fp, deg N radian N radian N.m N.m ); ic=(int)(360.0/del); for(i=0;i=(-0.31825)&p510)fenexf1=-9000.0;elsefenexf1=0.0;运行结果The Kineto-static

20、Analysis of a six-bar LinkaseNo HETAL fr1 sita1 fr2 sita2 tb tb1 deg N radian N radian N.m N.m 1 -165.522 1465.961 14.478 627.211 -171.015 0.000 0.000 2 -180.522 1127.912 14.005 446.916 179.556 -28.293 -28.293 3 -195.522 15767.371 12.723 6867.477 -148.658 -746.184 -746.184 4 -210.522 15065.899 10.81

21、8 6002.447 -154.166 -995.161 -995.161 5 -225.522 14516.560 8.455 5324.070 -162.398 -1174.074 -1174.074 6 -240.522 14074.069 5.769 4865.571 -172.533 -1288.627 -1288.627 7 -255.522 13725.674 2.881 4635.881 176.812 -1344.548 -1344.548 8 -270.522 13479.810 -0.104 4606.793 167.177 -1347.945 -1347.945 9 -

22、285.522 13355.544 -3.086 4730.512 159.624 -1304.217 -1304.21710 -300.522 13373.618 -5.964 4969.198 154.494 -1216.382 -1216.38211 -315.522 13549.209 -8.631 5310.503 151.601 -1083.634 -1083.63412 -330.522 13886.197 -10.968 5763.750 150.487 -900.827 -900.82713 -345.522 14373.388 -12.834 6345.168 150.61

23、5 -659.496 -659.49614 -360.522 1208.101 165.934 504.947 7.117 28.293 28.29315 -375.522 1602.230 165.525 711.386 -1.585 -2.929 -2.92916 -390.522 1999.234 166.143 957.246 -5.892 -57.335 -57.33517 -405.522 2257.789 167.966 1159.948 -6.456 -124.576 -124.57618 -420.522 2148.950 171.064 1170.305 -3.036 -1

24、68.380 -168.38019 -435.522 1409.272 175.286 812.249 7.743 -133.095 -133.09520 -450.522 7.272 -179.826 244.641 88.754 -0.727 -0.72721 -465.522 1484.336 5.037 866.590 169.198 138.980 138.98022 -480.522 2280.194 9.192 1242.566 -177.203 175.402 175.40223 -495.522 2270.202 12.203 1157.293 -170.237 121.22

25、4 121.22424 -510.522 1883.416 13.937 884.112 -167.926 50.459 50.45925 -525.522 1465.961 14.478 627.211 -171.015 -0.000 -0.000运行图像图6-1 图6-1图6-2图6-27.飞轮的转动惯量计算7-1飞轮转动惯量的计算方法 1、一周期内驱动力矩功等于阻力功，所以：t*Td=1/2(Tr0+Tr1)+1/2(Tr1+Tr2)+.+1/2(Trn-1+Trn)因为Trn=Tr0，所以由上式可得：Td=/tTri=1/nTri2、间隔i-1、i内的盈亏功变化量E:E=Td-1/2*

26、（Tri+Tri-1）3、计算各点的盈亏功Ei：Ei=Ei-1+Ei4、找出最大和最小盈亏功：Emax,Emin5、计算飞轮的转动惯量Jf: Jf=(Emax-Emin)/W， 7-2飞轮转动惯量的程序及运行结果主程序/* Note:Your choice is C IDE */#include”stdio.h”#include”graphics,h” /*图形库*/#includesubk.c /*运动分析子程序*/ #includesubf.c /*动态静力分析子程序*/ main() static double p202,vp202,ap202,del; static double t1

27、0,w10,e10,Tr370; static double fr202,fk202,pk202,fe202; static int ic; int j; double Td,sum1=0.0,D370,E370,Max,Min,Jf,q; double r12,r13,r24,r45,r26,r58,r57,tb; double pi,dr; double r2,vr2,ar2; int i; FILE *fp; E0=0.0; Max=E0; Min=E0; q=0.03; sm1=0.0;sm2=0.0;sm3=30.0;sm4=0.0;sm5=95.0; sj1=0.0;sj2=0.0

28、;sj3=0.7;sj4=0.0; r12=0.4;r13=0.1;r24=0.75;r45=0.15;r26=0.375; r58=sqrt(0.1*0.1+0.25*0.25); r57=sqrt(0.15*0.15+0.05*0.05); pi=4.0*atan(1.0); dr=pi/180.0; t6=0.0;w6=0.0;e6=0.0; w1=-2*pi;e1=0.0;del=15.0; t6=0.0*dr; p21=0.0; p22=0.0; p11=0.0; p12=0.4; p91=-100; p92=0.738; printf(n The Kineto-statia Ana

29、lysis of Ntbcn ); printf(nNo Theta1 Tr E n ); printf( deg. N.M 1/sn ); /*在屏幕上写表头*/ if(fp=fopen(file1,w)=NULL) printf(Cant open this file.n); exit(0); /*建立并打开文件file1*/ fprintf(fp,n The Kineto-statia Analysis of Ntbcn); fprintf(fp,nNo Theta1 Ter En); fprintf(fp, deg. N.M Jn); /*在文件中写表头*/ ic=(int)(360.

30、0/del); for(i=0;i=ic;i+) /*建立循环，调用运动分析子程序*/ t1=(-(pi-asin(r13/r12)+(-i)*del*dr; bark(1,3,0,1,r13,0.0,0.0,t,w,e,p,vp,ap); rprk(1,2,3,3,2,0.0,&r2,&vr2,&ar2,t,w,e,p,vp,ap); bark(2,4,0,3,r24,0.0,0.0,t,w,e,p,vp,ap); rrpk(-1,4,9,5,4,5,6,r45,&r2,&vr2,&ar2,t,w,e,p,vp,ap); bark(2,6,0,3,r26,0.0,0.0,t,w,e,p,vp

31、,ap); bark(5,0,7,5,0.0,r58,-161.57*dr,t,w,e,p,vp,ap); bark(5,0,8,5,0.0,r57,165.96*dr,t,w,e,p,vp,ap); rrpf(4,9,5,0,7,0,8,8,4,5,p,vp,ap,t,w,e,fr); rprf(2,3,6,0,4,0,0,3,2,p,vp,ap,t,w,e,fr,fk,pk); barf(1,0,3,1,p,ap,e,fr,&tb); /*r2,vr2,ar2,tb在子程序中已定义为指针变量，所以前面要加&符号*/ Tri=tb; Di=t1; for(j=1;j=ic;j+) sum1=sum1+Trj-1; Td=sum1/ic; for(j=1;j=ic;j+) Ej=Ej-1+del*dr*(Td-0.5*(Trj+Trj-1); for(j=1;j=ic;j+) if(Max=Ej) Min=Ej; J