机械原理课程设计--铰链式颚式破碎机方案分析.doc

东北大学机械原理课程设计 铰链式颚式破碎机方案分析机械原理课程设计说明书题目： 铰链式颚式破碎机方案分析 班 级 ：姓 名 ：学 号 ：指导教师 ：成 绩 ：2012年 月 日目 录一、设计题目 1二、已知条件及设计要求 12.1已知条件 12.2设计要求 1三、机构的结构分析 23.1六杆铰链式破碎机 23.2四杆铰链式破碎机 2 四、机构的运动分析 24.1六杆铰链式颚式破碎机的运动分析 24.2四杆铰链式颚式破碎机的运动分析 5五、机构的动态静力分析 85.1六杆铰链式颚式破碎机的静力分析 85.2四杆铰链式颚式破碎机的静力分析 12六、工艺阻力函数及飞轮的转动惯量函数 156.1工艺阻力函数程序 156.2 六杆飞轮转动惯量函数程序 166.3 四杆飞轮转动惯量函数程序 17七、对两种机构的综合评价 19八、主要的收获和建议 20九、参考文献 20东北大学机械原理课程设计 铰链式颚式破碎机方案分析一、设计题目及已知条件1.1设计题目：铰链式颚式破碎机方案分析(旧教材题目七) (a) 六杆铰链式破碎机 (b) 工艺阻力 (c) 四杆铰链式破碎机1.2、已知条件：图(a)所示为六杆铰链式破碎机方案简图。主轴1的转速为n1 = 170r/min，各部尺寸为：lO1A = 0.1m, lAB = 1.250m, lO3B = 1m, lBC = 1.15m, lO5C = 1.96m, l1=1m, l2=0.94m, h1=0.85m, h2=1m。各构件质量和转动惯量分别为：m2 = 500kg, Js2 = 25.5kgm2, m3 = 200kg, Js3 = 9kgm2, m4 = 200kg, Js4 = 9kgm2, m5=900kg, Js5=50kgm2, 构件1的质心位于O1上，其他构件的质心均在各杆的中心处。D为矿石破碎阻力作用点，设LO5D = 0.6m，破碎阻力Q在颚板5的右极限位置到左极限位置间变化，如图(b)所示，Q力垂直于颚板。图(c)是四杆铰链式颚式破碎机方案简图。主轴1 的转速n1=170r/min。lO1A = 0.04m, lAB = 1.11m, l1=0.95m, h1=2m, lO3B=1.96m，破碎阻力Q的变化规律与六杆铰链式破碎机相同，Q力垂直于颚板O3B，Q力作用点为D，且lO3D = 0.6m。各杆的质量、转动惯量为m2 = 200kg, Js2=9kgm2，m3 = 900kg, Js3=50kgm2。曲柄1的质心在O1 点处，2、3构件的质心在各构件的中心。二、设计要求试比较两个方案进行综合评价。主要比较以下几方面：1. 进行运动分析，画出颚板的角位移、角速度、角加速度随曲柄转角的变化曲线。2. 进行动态静力分析，比较颚板摆动中心运动副反力的大小及方向变化规律，曲柄上的平衡力矩大小及方向变化规律。3. 飞轮转动惯量的大小。三、机构的结构分析3.1六杆铰链式颚式破碎机六杆铰链式破碎机拆分为机架和主动件构件组成的RRR杆组，构件组成的RRR杆组。 3.2 四杆铰链式破碎机破碎机四杆铰链式破碎机破碎机拆分为机架和主动件，构件组成的RRR杆组。 四、机构的运动分析4.1 六杆铰链式颚式破碎机的运动分析（1）调用bark函数对主动件进行运动分析。见表4.1表4.1形式参数n1n2n3kr1r2gamtwepvpap实值1201r120.00.0twe pvpap（2）调用rrrk函数对由构件组成的RRR杆组运行分析。见表4.2.表4.2形式参数mn1n2n3k1k2r1r2twepvpap实值-126323r23r36twepvpap(3)调用rrrk函数对由构件组成的RRR杆组运行分析。见表4.3.表4.3形式参数mn1n2n3k1k2r1r2twepvpap实值-153454R45R34twepvpap (4)调用bark函数对质点7运动分析。见表4.4.表4.4形式参数n1n2n3kr1r2gamtwepvpap实值50750.0R570.0twepvpap （5）调用bark函数对质点9运动分析。见表4.5.表4.5形式参数n1n2n3kr1r2gamtwepvpap实值20920.0R23/20.0twepvpap （6）调用bark函数对质点10运动分析。见表4.6.表4.6形式参数n1n2n3kr1r2gamtwepvpap实值601030.0r36/20.0twepvpap (7) 调用bark函数对质点11运动分析。见表4.7.表4.7形式参数n1n2n3kr1r2gamtwepvpap实值301140.0R34/20.0twepvpap （8）调用bark函数对质点12运动分析。见表4.8.表4.8形式参数n1n2n3kr1r2gamtwepvpap实值501250.0R45/20.0twepvpap主程序：#includegraphics.h#includesubk.c#includedraw.cmain()static double p202,vp202,ap202,del; static double t10,w10,e10,tbdraw370,tb1draw370,pdraw370,vpdraw370,apdraw370; static int ic; double r12,r23,r34,r45,r36,r57; double pi,dr; double r5,vr5,ar5; int i; FILE *fp; char *m=p,vp,ap; r12=0.1; r23=1.25; r34=1.15;r45=1.96; r36=1.0; r57=0.6; w1=-17.79; e1=0.0; del=1.0; p11=0.0; p12=0.0; p61=0.94; p62=-1.0; p51=-1.0; p52=0.85; pi=4.0*atan(1.0); dr=pi/180.0;printf(n the kinematic parameters of point 5n);printf(NO THETA1 t5 t5 t5n);if(fp=fopen(filel,w)=NULL)printf(Cant open this file.n);exit(0);fprintf(fp,n the kinematic parameters of point 5n);fprintf(fp,NO THETA1 t5 t5 t5n); fprintf(fp, deg m m/s m/s/s);ic = (int)(360.0/del);for (i=0;i=ic;i+) t1=(-i)*del*dr; bark(1,2,0,1,r12,0.0,0.0,t,w,e,p,vp,ap); rrrk(-1,2,6,3,2,3,r23,r36,t,w,e,p,vp,ap); rrrk(-1,5,3,4,5,4,r45,r34,t,w,e,p,vp,ap); bark(5,0,7,5,0.0,r57,0.0,t,w,e,p,vp,ap); bark(2,0,9,2,0.0,r23/2,0.0,t,w,e,p,vp,ap); bark(6,0,10,3,0.0,r36/2,0.0,t,w,e,p,vp,ap); bark(3,0,11,4,0.0,r34/2,0.0,t,w,e,p,vp,ap); bark(5,0,12,5,0.0,r45/2,0.0,t,w,e,p,vp,ap); printf(n%2d%12.3f%12.3f%12.3f%12.3f,i+1,t1/dr,t5,w5,e5);fprintf(fp,n%2d%12.3f%12.3f%12.3f%12.3f,i+1,t1/dr,t5,w5,e5);pdrawi=t5;vpdrawi=w5;apdrawi=e5;if (i%16)=0)getch(); fclose(fp); getch(); draw1(del,pdraw,vpdraw,apdraw,ic,m);运行结果： the kinematic parameters of point 5 NO THETA1 t5 w5 e5 deg m m/s m/s/s 1 0.000 -1.658 0.346 3.950 2 -15.000 -1.653 0.391 1.999 3 -30.000 -1.647 0.400 -0.931 4 -45.000 -1.641 0.362 -4.349 5 -60.000 -1.637 0.274 -7.495 6 -75.000 -1.633 0.146 -9.599 7 -90.000 -1.632 -0.001 -10.169 8 -105.000 -1.633 -0.145 -9.152 9 -120.000 -1.637 -0.265 -6.89410 -135.000 -1.641 -0.345 -3.97511 -150.000 -1.646 -0.381 -1.00612 -165.000 -1.652 -0.377 1.51713 -180.000 -1.657 -0.340 3.29314 -195.000 -1.662 -0.284 4.23215 -210.000 -1.666 -0.220 4.43016 -225.000 -1.668 -0.156 4.11617 -240.000 -1.670 -0.10 3.57918 -255.000 -1.671 -0.051 3.10119 -270.000 -1.672 -0.007 2.89420 -285.000 -1.672 0.036 3.05921 -300.000 -1.671 0.085 3.56622 -315.000 -1.669 0.142 4.24223 -330.000 -1.667 0.209 4.78524 -345.000 -1.663 0.281 4.81125 -360.000 -1.658 0.346 3.9504.2四杆铰链式颚式破碎机的运动分析（9）调用bark函数对主动件进行运动分析。见表4.9表4.9形式参数n1n2n3kr1r2gamtwepvpap实值1201r120.00.0twepvpap（10）调用rrrk函数对由构件组成的RRR杆组运行分析。见表4.10.表4.10形式参数mn1n2n3k1k2r1r2twepvpap实值-142332R34R23twepvpqp（11）调用bark函数对质点6,7和作用点5运动分析。见表4.11.表4.11形式参数n1n2n3kr1r2gamtwepvpap实值40530R450.0twepvpap 表4.12形式参数n1n2n3kr1r2gamtwepvpap实值20620R23/20.0twepvpap 表4.13形式参数n1n2n3kr1r2gamtwepvpap实值40730R34/20.0twepvpap主程序：#includegraphics.h#includesubk.c#includedraw.cmain()static double p202,vp202,ap202,del; static double t10,w10,e10,pdraw370,vpdraw370,apdraw370; static int ic;double r12,r23,r34,r45; double pi,dr;double r5,vr5,ar5;int i;FILE *fp;char *m=p,vp,ap; r12=0.04; r23=1.11; r34=1.96;r45=0.6;w1=-17.79; e1=0.0; del=1.0; pi=4.0*atan(1.0);dr=pi/180.0;p11=0; p12=0;p41=-0.95; p42=2;printf(n the kinematic parameters of point 5n);printf(NO THETA1 t3 w3 e3n);if(fp=fopen(filel,w)=NULL) printf(Cant open this file.n); exit(0);fprintf(fp,n the kinematic parameters of point 5n);fprintf(fp,NO THETA1 t3 w3 e3n); fprintf(fp, deg m m/s m/s/s);ic = (int)(360.0/del);for (i=0;i=ic;i+) t1=(-i)*del*dr;bark(1,2,0,1,r12,0.0,0.0,t,w,e,p,vp,ap);rrrk(-1,4,2,3,3,2,r34,r23,t,w,e,p,vp,ap);bark(4,0,5,3,0.0,r45,0.0,t,w,e,p,vp,ap); bark(2,0,6,2,0.0,r23/2,0.0,t,w,e,p,vp,ap);bark(4,0,7,3,0.0,r34/2,0.0,t,w,e,p,vp,ap);printf(n%2d%12.3f%12.3f%12.3f%12.3f,i+1,t1/dr,t3,w3,e3);fprintf(fp,n%2d%12.3f%12.3f%12.3f%12.3f,i+1,t1/dr,t3,w3,e3);pdrawi=t3;vpdrawi=w3;apdrawi=e3;if (i%16)=0)getch(); fclose(fp); getch();draw1(del,pdraw,vpdraw,apdraw,ic,m);运行结果：the kinematic parameters of point 5NO THETA1 t3 w3 e3 deg m m/s m/s/s 1 0.000 -1.632 0.014 -6.223 2 -15.000 -1.632 -0.077 -6.090 3 -30.000 -1.634 -0.163 -5.584 4 -45.000 -1.637 -0.239 -4.724 5 -60.000 -1.641 -0.301 -3.549 6 -75.000 -1.646 -0.343 -2.114 7 -90.000 -1.651 -0.362 -0.501 8 -105.000 -1.656 -0.357 1.191 9 -120.000 -1.661 -0.327 2.84410 -135.000 -1.666 -0.274 4.33311 -150.000 -1.669 -0.201 5.53612 -165.000 -1.671 -0.113 6.34913 -180.000 -1.672 -0.016 6.69414 -195.000 -1.672 0.082 6.53615 -210.000 -1.670 0.174 5.88516 -225.000 -1.667 0.253 4.80117 -240.000 -1.663 0.313 3.37918 -255.000 -1.658 0.351 1.74419 -270.000 -1.653 0.364 0.03020 -285.000 -1.647 0.352 -1.63721 -300.000 -1.642 0.317 -3.14522 -315.000 -1.638 0.261 -4.40923 -330.000 -1.635 0.188 -5.36724 -345.000 -1.632 0.105 -5.98025 -360.000 -1.632 0.014 -6.223五、机构的动态静力分析5.1六杆铰链式颚式破碎机的动态静力分析（1）调用bark函数对主动件进行运动分析。见表4.1.（2）调用rrrk函数对由构件组成的RRR杆组运行分析。见表4.2.（3）调用rrrk函数对由构件组成的RRR杆组运行分析。见表4.3.（4）求质心8,9,10,11,12和作用点7的运动分析。见表4.4-表4.8.（5）调用rrrf对构件组成的RRR杆组进行静力分析。见表5.1.表5.1形式参数n1n2n3ns1ns2nn1nn2nexfk1k2pvpaptwefr实值534121170754pvpaptwefr（6）调用rrrf对构件组成的RRR杆组进行静力分析。见表5.2.表5.2形式参数n1n2n3ns1ns2nn1nn2nexfk1k2pvpaptwefr实值26391003023pvpaptwefr（5）调用barf函数对主动件进行静力分析。见表5.3表5.3形式参数n1n2nn1k1papefrtb实值1121Papefr&tb主程序：#include graphics.h#include subk.c#include subf.c#include draw.cmain()static double p202,vp202,ap202,del; static double t10,w10,e10;static double sita1370,fr1draw370,sita2370,fr2draw370,sita3370,fr3draw370,tbdraw370,tb1draw370;static double fr202,fe202;static int ic;double r12,r23,r34,r45,r36,r57;int i;double pi,dr;double fr1,fr6,bt1,bt6,fr5,bt5,we1,we2,we3,we4,we5,tb,tb1;FILE*fp; sm1=0.0;sm2=500.0;sm3=200.0;sm4=200.0;sm5=900.0; sj1=0.0;sj2=25.5;sj3=9.0;sj4=9.0;sj5=50.0;r12=0.1; r23=1.25; r34=1.15; r45=1.96;r57=0.6;r36=1.0; pi=4.0*atan(1.0);dr=pi/180.0;w1=-170*2*pi/60; e1=0.0; del=1.0; p11=0.0; p12=0.0; p61=0.94;p62=-1.0; p51=-1.0;p52=0.85; printf(n The Kineto-static Analysis of a Six-bar Linkasen);printf( No THETAL fr1 sita1 fr1 sita1 tb tb1n);printf( (deg.) (N) (deg.) (N.m) (N.m)n);if(fp=fopen(filel,w)=NULL)printf(Cant open this file./n); exit(0);fprintf(fp,n The Kineto-static Analysis of a Six-bar Linkasen);fprintf(fp,No THETAL fr1 sita1 fr5 sita5 tb tb1n );fprintf(fp, eg N radian N.m N.m n );ic=(int)(360.0/del);for(i=0;i=-273.0&t1/dr=-90.0)fenexf1=85000.0/(183.0*dr)*(t1+90*dr)*sin(t5);fenexf2=-85000.0/(183.0*dr)*(t1+90*dr)*cos(t5);else fenexf1=0.0;fenexf2=0.0;运行结果： The Kineto-static Analysis of a Six-bar Linkase No THETAL fr1 sita1 fr5 sita5 tb tb1 eg N radian N.m N.m 1 0.000 9904.580 77.690 534.273 534.273 2 -15.000 10248.086 82.670 1038.104 1038.104 3 -30.000 10522.852 89.576 1434.513 1434.513 4 -45.000 10757.314 97.329 1547.760 1547.760 5 -60.000 10967.175 104.339 1270.987 1270.987 6 -75.000 11112.158 109.009 644.228 644.228 7 -90.000 11132.496 110.330 -144.608 -144.608 8 -105.000 13340.278 128.365 -883.585 -883.585 9 -120.000 16136.071 138.884 -1407.147 -1407.147 10 -135.000 19120.284 145.429 -1625.007 -1625.007 11 -150.000 22280.557 150.044 -1557.275 -1557.275 12 -165.000 25724.638 153.629 -1290.422 -1290.422 13 -180.000 29550.567 156.506 -928.575 -928.575 14 -195.000 33789.984 158.781 -562.676 -562.676 15 -210.000 38396.542 160.527 -258.686 -258.686 16 -225.000 43263.326 161.841 -55.813 -55.813 17 -240.000 48254.151 162.844 33.550 33.550 18 -255.000 53236.591 163.653 24.697 24.697 19 -270.000 58109.329 164.367 -39.702 -39.702 20 -285.000 8481.823 78.617 -205.306 -205.306 21 -300.000 8583.465 77.292 -338.729 -338.729 22 -315.000 8793.293 75.658 -361.459 -361.459 23 -330.000 9113.158 74.602 -227.576 -227.576 24 -345.000 9506.210 75.059 80.824 80.824 25 -360.000 9904.580 77.690 534.273 534.273 5.2四杆铰链式颚式破碎机的动态静力分析（1）调用bark函数对主动件进行运动分析。见表4.9（2）调用rrrk函数对由构件组成的RRR杆组运行分析。见表4.10.（3）调用bark函数对质点6,7和作用点5运动分析。见表4.11.（4）调用rrrf对构件组成的RRR杆组进行静力分析。见表5.4表5.4形式参数n1n2n3ns1ns2nn1nn2nexfk1k2pvpaptwefr实值4237650532pvpaptwef（5）调用barf函数对主动件进行静力分析。见表5.5表5.5形式参数n1n2nn1k1papefrtb实值1121papefr&tb主程序：#include graphics.h#include subk.c#include subf.c#include draw.cmain()static double p202,vp202,ap202,del;static double t10,w10,e10;static double sita1370,fr1draw370,sita2370,fr2draw370,sita3370,fr3draw370,tbdraw370,tb1draw370;static double fr202,fe202;static int ic;double r12,r23,r34,r45;int i;double pi,dr;double fr1,fr4,bt1,bt4,we1,we2,we3,tb,tb1;FILE*fp; sm1=0.0;sm2=200.0;sm3=900.0; sj1=0.0;sj2=9.0;sj3=50.0 ;r12=0.04; r23=1.11; r34=1.96; r45=0.6;pi=4.0*atan(1.0);w1=-170*2*pi/60; e1=0.0; del=1;dr=pi/180.0;p11=0.0;p12=0.0;p41=-0.95;p42=2.0; printf(n The Kineto-static Analysis of a Six-bar Linkasen);printf( NO THETA1 FR1 BT1 TB TB1n);printf( (deg.) (N)