牛头刨床机构方案分析--东北大学车辆工程机械原理课程设计.doc

机械原理课程设计说明书题目： 牛头刨床机构方案分析 班 级 ：姓 名 ： 学 号 ：指导教师 ：成 绩 ： 年 月 日东北大学机械原理课程设计 牛头刨床机构方案分析目 录 一、设计题目11.1已知条件11.2设计要求1二、滑枕初始位置及行程H的确定方法 1三、机构的结构分析2四、机构的运动分析24.1所调用的杆组子程序中虚参与实参对照表24.2源程序 24.3运行结果5五、机构的动态静力分析 65.1所调用的杆组子程序中虚参与实参对照表65.2源程序 75.3运行结果 10六、飞轮的转动惯量126.1源程序126.2运行结果 14七、主要的收获和建议15八、参考文献 15东北大学机械原理课程设计 牛头刨床机构方案分析一、设计题目1.1已知条件（图a,b）所示为两种牛头刨床主机构的运动简图，已知，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。工艺阻力Q如图所示，Q=9000N。主轴1的转速为60r/min(顺时针方向)，许用运转不均匀系数=0.03。1.2设计要求1.2.1进行运动分析，绘制出滑枕5的位移、速度和加速度线图。1.2.2进行动态静力分析，绘制固定铰链处的反力矢端图和平衡力矩Td的线图。1.2.3以主轴为等效构件，计算飞轮转动惯量。l0l1H二、滑枕初始位置及行程H的确定方法构件1和3在左侧垂直的位置为初始位置，此时=+arcsin(l1/l0)行程H等于B点划过的弧所对应的弦长即H=2l1l3/l0b同于a三、机构的结构分析主动件13RPR杆组39247RPP杆组8456四、机构的运动分析4.1所调用的杆组子程序中虚参与实参对照表4.1.1调用bark函数，求点3的运动参数形式参数n1n2n3kr1r2gamtwepvpap实值1301r130.00.0twepvpap4.1.2调用rprk函数求构件的运动参数形式参数mn1n2k1k2r1r2vr2ar2twepvpap实值123230.0&r2&vr2&ar2twepvpap4.1.3调用bark函数求点4的运动参数形式参数n1n2n3kr1r2gamtwepvpap实值2402r240.00.0twepvpap4.1.4调用rppk函数求点5的运动参数形式参数n1n2n3n4k1k2k3r1gam1gam2实值46454560.00.0gam2r2vr2ar2r3vr3ar3twepvpap&r2&vr2&ar2&r3&vr3&ar3twepvpap4.2源程序#include stdio.h#include graphics.h /*图形库*/#include subk.c /*运动分析子程序*/ #include draw.c /*绘图子程序*/ main()static double p202,vp202,ap202,del;static double t10,w10,e10,pdraw370,vpdraw370,apdraw370;static int ic;double r13,r24,r45, r29,r57,r58,gam1,r28; double gam2,r59; double pi,dr,gam7,gam8; double r2,vr2,ar2,r3,vr3,ar3; int i;FILE *fp; /*定义变量*/ char *m=p,vp,ap;r13=0.1;r24=0.75;r45=0.15 ;r29=0.375;gam1=atan(1/4); gam2=-atan(1/3); t6=0.0; w6=0.0; e6=0.0;e1=0.0; del=15.0; /*赋值*/ pi=4.0*atan(1.0); /*求*/ dr=pi/180.0; /*求弧度*/ gam7=gam1*dr; gam8=gam2*dr; /*变角度为弧度*/ w1=-2*pi; p11=0.0; p12=0.0; p21=0.0; p22=-0.4; p61=0.0; p62=0.338; /*赋定铰链值*/ printf( n TheKinematic Parametersof Point 6n);printf(No THETA1 S6 V6 A6n);printf( deg m m/s m/s/sn); /*在屏幕上写表头*/ if(fp=fopen(first file,w)=NULL) printf( Cant open this file.n);exit(0); /*建立并打开文件first file*/ fprintf(fp, n TheKinematic Parametersof Point 6n); fprintf(fp,No THETA1 S6 V6 A6n); fprintf(fp, deg m m/s m/s/sn); /*在文件zouyes k result中写表头*/ic=(int)(360.0/del);for(i=0;i0) fenexf1=-9000.0; else fenexf1=0.0;main()static double p202,vp202,ap202,del;static double t10,w10,e10,tbdraw370,tb1draw370;static double sita1370,fr1draw370,sita2370,fr2draw370,sita3370,fr3draw370,fr3,bt3; static double fr202,fe202,fk202,pk202; static int ic; double r13,r24,r29,r57,r58,r45,gam,gam1,gam2,gam3,gam4,tb; double pi,dr; double r2,vr2,ar2,r3,vr3,ar3; double fr1,bt1,fr2bt2we1,we2,we3,we4,we5,tb1; int i; FILE*fp; char *m=tb,tb1,fr1,fr2; sm1=0.0; sm2=30.0; sm3=0.0; sm4=0.0; sm5=95.0; sj2=0.7; r13=0.1; r24=0.75; r57=0.158;r29=r24/2;r58=0.412;del=15.0; gam=0.0; gam1=0.0; gam2=90.0; gam3=atan(1/4); gam4=-atan(1/3); /*赋值*/ pi=4.0*atan(1.0); /*求*/ dr=pi/180.0; /*求弧度*/ t5=0.0; w5=0.0; e5=0.0; w1=-2*pi; e1=0.0; gam*=dr; gam1*=dr; gam2*=dr; gam3*=dr; gam4*=dr; /*变角度为弧度*/ p11=0.0;p12=0.0; p21=0.0;p22=-0.4; p61=1.0;p62=0.338; /*赋定铰链值*/ printf(n The Kineto-static Analysis of Ntbcn); printf( NO THETA1 FR1 BT1 FR2 BT2 TB TB1n); printf( deg N deg N deg N.m N.mn); /*在屏幕上写表头*/ if(fp=fopen(second file,w)=NULL) printf(Cant open this file./n); exit(0); /*建立并打开文件second file*/ printf(fp,n The Kineto-static Analysis of Ntbcn); printf(fp, NO THETA1 FR1 BT1 FR2 BT2 TB TB1n); printf(fp, deg N deg N deg N.m N.mn); /*在文件中写表头*/ ic=(int)(360.0/del); for(i=0;i0) fenexf1=-9000.0; else fenexf1=0.0;main()static double p202,vp202,ap202,del;static double t10,w10,e10,Tr370;static double fr202,fe202,fk202,pk202;static int ic;int j;double Td,sum1=0.0,D370,E370,Max,Min,Jf,q;double r13,r24,r29,r57,r58,tb,gam2,gam3,gam4;double pi,dr;double r2,vr2,ar2,r3,vr3,ar3;int i;FILE *fp;E0=0.0; Max=E0; Min=E0; q=0.03;sm1=0.0; sm2=30.0; sm3=0.0; sm4=0.0; sm5=95.0; sj2=0.7;r13=0.1; r24=0.75; r29=r24/2; r57=0.158;r58=0.412;del=15.0; gam2=90.0; gam3=atan(1/3); gam4=-atan(1/4); /*赋值*/pi=4.0*atan(1.0); /*求*/dr=pi/180.0; /*求弧度*/t5=0.0;w5=0.0;e5=0.0;w1=-2*pi; e1=0.0; gam2*=dr; gam3*=dr; gam4*=dr; /*变角度为弧度*/ p11=0.0;p12=0.0;p21=0.0;p22=-0.4;p61=0.0;p62=0.338; /*赋定铰链值*/printf(n The Kineto-statia Analysis of Ntbcn);printf(nNo Theta1 Tr E n); printf( deg N.M 1/sn); /*在屏幕上写表头*/ if(fp=fopen(third file,w)=NULL) printf(Cant open this file./n); exit(0); /*建立并打开文件third file*/ fprintf(fp,n The Kineto-statia Analysis of Ntbcn); fprintf(fp,nNo Theta1 Tr E n); fprintf(fp, deg N.M 1/sn); /*在文件中写表头*/ ic=(int)(360.0/del); for(i=0;i=ic;i+) /*建立循环，调用运动分析子程序*/ t1=(double)(i)*del*dr+3.4; bark(1,3,0,1,r13,0.0,0.0,t,w,e,p,vp,ap); rprk(1,2,3,2,3,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); rppk(4,6,4,5,4,5,6,0.0,0.0,gam2,&r2,&vr2,&ar2,&r3,&vr3,&ar3,t,w,e,p,vp,ap); bark(5,0,7,5,0.0,r57,gam3, t, w,e,p,vp,ap); bark(5,0,8,5,0.0,r58,gam4, t, w,e,p,vp,ap); bark(2,0,9,3,0.0,r29,0.0,0.0, t, w,e,p,vp,ap); rppf(4,6,4,5,4,5,0,8,8,4,5,6, p,vp,ap,t,w,e,fr,fk,pk); rprf(2,3,9,3,4,0,0,2,3,p,vp,ap,t,w,e,fr,fk,pk); barf(1,1,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; Jf=(Max-Min)/(w1*w1*q); for(j=0;j=ic;j+) printf(n%3d%13.3f%16.3f%16.3fn,j,Dj,Trj,Ej); fprintf(fp,n%3d%13.3f%16.3f%16.3fn,j,Dj,Trj,Ej); if(j%10=0) getch(); printf(n Jf=%4.3f Td=%4.3f,Jf,Td); /*把运动结果写屏幕上*/ fprintf(fp,n Jf=%4.3f Td=%4.3f,Jf,Td); /*把运动结果写入文件中*/ fclose(fp); getch();6.2运行结果 The Kineto-statia Analysis of NtbcNo Theta1 Tr E deg N.M 1/s 0 3.400 0.000 0.000 1 3.662 47.155 -146.406 2 3.924 109.657 -307.390 3 4.185 155.436 -483.552 4 4.447 124.715 -662.803 5 4.709 1.878 -821.197 6 4.971 -122.522 -944.810 7 5.233 -155.862 -1046.952 8 5.494 -111.280 -1151.663 9 5.756 -48.636 -1