颚式破碎机设计项目说明指导书最终版.doc

机械原理课程设计阐明书题目： 铰链式颚式破碎机方案分析 班 级 ：机械工程1405班姓 名 ：施坤军学 号 ：3104指引教师 ：杨强成 绩 ： 9月 29日目 录 一 设计题目1二 已知条件及设计规定12.1已知条件12.2设计规定1三. 机构构造分析23.1六杆铰链式破碎机23.2四杆铰链式破碎机2四. 机构运动分析24.1六杆铰链式颚式破碎机运动分析24.2编写主程序并运营24.3四杆铰链式颚式破碎机运动分析54.4编写主程序并运营5五.机构动态静力分析85.1六杆铰链式颚式破碎机静力分析85.2编写主程序并运营95.3四杆铰链式颚式破碎机静力分析135.4编写主程序并运营14六. 工艺阻力函数及飞轮转动惯量函数 176.1工艺阻力函数程序 176.2飞轮转动惯量函数程序17七 .对两种机构综合评价21八 . 重要收获和建议21九 . 参照文献21一.设计题目：铰链式颚式破碎机方案分析二.已知条件及设计规定2.1 已知条件 图1.1 六杆铰链式破碎机 图1.2 工艺阻力 图1.3 四杆铰链式破碎机图(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构件质心在各构件中心。2.2 设计规定试比较两个方案进行综合评价。重要比较如下几方面：1. 进行运动分析，画出颚板角位移、角速度、角加速度随曲柄转角变化曲线。2. 进行动态静力分析，比较颚板摆动中心运动副反力大小及方向变化规律，曲柄上平衡力矩大小及方向变化规律。3. 飞轮转动惯量大小。三.机构构造分析3.1六杆铰链式破碎机 六杆铰链式粉碎机拆分为机架和积极件，构件构成RRR杆组，构件构成RRR杆组。21 1042311 8796533.2四杆铰链式破碎机四杆铰链式破碎机拆分为机架和积极件，构件构成RRR杆组。21 753642四.机构运动分析4.1六杆铰链式颚式破碎机运动分析。1) 调用bark函数求积极件运动参数。形式参数n1n2n3kr1r2gamtwepvpap实 值1201r120.00.0twepvpap2）调用rrrk函数求、构件构成RRR杆组进行运动分析。形式参数mn1n2n3k1k2r1r2t wepvpap实 值-124323r23r34twepvpap3） 调用rrrk函数对、构件构成RRR杆组进行运动分析。形式参数mn1n2n3k1k2r1r2t wepvpap实 值136545r35r56twepvpap4.2编写主程序并运营。按一定步长，变化积极件位置角度，使其在0-360变化，便可求出机构各点在整个运动循环内运动参数并打印输出。(1) 主程序。#includegraphics.h#includesubk.c#includedraw.cmain() static double p202,vp202,ap202; static double t10,w10,e10,del; static double pdraw370,vpdraw370,apdraw370; double r12,r23,r34,r35,r56,r68; double pi,dr; int i,ic; FILE *fp; char *m=p,vp,ap; pi=4.0*atan(1.0); dr=pi/180.0; r12=0.1;r23=1.25;r34=1.0;r35=1.15;r56=1.96;r68=0.6; w1=(-pi*17.0/3.0);e1=0.0;del=1.0; p11=0.0;p12=0.0; p41=0.94;p42=-1.0; p61=-1.0;p62=0.85; printf(nThe Kinematic Parameters of Member5n); printf(No THETA1 T5 W5 E5n); printf( deg rad rad/s rad/s/sn); if(fp=fopen(file_ww.txt,w)=NULL) printf(Cant open this file.n); exit(0); fprintf(fp,nThe Kinematic Parameters of Member5n); fprintf(fp,No THETA1 T5 W5 E5n); fprintf(fp, deg rad rad/s rad/s/sn); 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,4,3,2,3,r23,r34,t,w,e,p,vp,ap); rrrk(1,3,6,5,4,5,r35,r56,t,w,e,p,vp,ap); bark(6,0,7,5,0.0,r56/2,0.0,t,w,e,p,vp,ap); bark(6,0,8,5,0.0,r68,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); (2)运营成果。杆件5运动参数： The Kinematic Parameters of Pole5No THETA1 T5 W5 E5 deg rad rad/s rad/s/s 1 0.000 -1.658 0.346 3.956 2 -15.000 -1.653 0.392 2.002 3 -30.000 -1.647 0.400 -0.932 4 -45.000 -1.641 0.362 -4.355 5 -60.000 -1.637 0.274 -7.506 6 -75.000 -1.633 0.146 -9.612 7 -90.000 -1.632 -0.001 -10.183 8 -105.000 -1.633 -0.145 -9.165 9 -120.000 -1.637 -0.265 -6.90410 -135.000 -1.641 -0.345 -3.98111 -150.000 -1.646 -0.382 -1.00812 -165.000 -1.652 -0.377 1.51913 -180.000 -1.657 -0.341 3.29714 -195.000 -1.662 -0.284 4.23715 -210.000 -1.666 -0.220 4.43616 -225.000 -1.668 -0.156 4.12117 -240.000 -1.670 -0.10 3.58418 -255.000 -1.671 -0.051 3.10519 -270.000 -1.672 -0.007 2.89820 -285.000 -1.672 0.036 3.06321 -300.000 -1.671 0.085 3.57122 -315.000 -1.669 0.142 4.24723 -330.000 -1.667 0.209 4.79124 -345.000 -1.663 0.281 4.81725 -360.000 -1.658 0.346 3.956运动图形：4.3四杆铰链式颚式破碎机运动分析。1）调用bark函数求积极件运动参数。形式参数n1n2n3kr1r2gametwepvpap实 值1201r120.00.0twepvpap2）调用rrrk函数求、构件构成RRR杆组进行运动分析。形式参数mn1n2n3k1k2r1r2t wepvpap实 值126423r24r46twepvpap4.4编写主程序并运营。按一定步长，变化积极件位置角度，使其在0-360变化，便可求出机构各点在整个运动循环内运动参数并打印输出。(1)主程序。#includegraphics.h#includesubk.c#includedraw.cmain() static double p202,vp202,ap202; static double t10,w10,e10,del; static double pdraw370,vpdraw370,apdraw370; static int ic;double r12,r24,r46;double pi,dr;int i; FILE *fp; char *m=p,vp,ap; pi=4.0*atan(1.0); dr=pi/180.0; r12=0.04;r24=1.11;r46=1.96; w1=(-(17.0/3.0)*pi);e1=0.0;del=1.0; p11=0.0; p12=0.0; p61=-0.95; p62=2.0; printf(nThe Kinematic Parameters of Member3n); printf(No THETA1 T3 W3 E3n); printf( deg rad rad/s rad/s/sn); if(fp=fopen(file_ww.txt,w)=NULL) printf(Cant open this file.n); exit(0); fprintf(fp,nThe Kinematic Parameters of Member3n); fprintf(fp,No THETA1 T3 W3 E3n); fprintf(fp, deg rad rad/s rad/s/sn); 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,4,2,3,r24,r46,t,w,e,p,vp,ap); bark(2,0,3,2,0.0,r24/2,0.0,t,w,e,p,vp,ap); bark(6,0,5,3,0.0,r46/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/dr,w3,e3)； fprintf(fp,n%2d %12.3f %12.3f %12.3f %12.3f,i+1,t1/dr,t3/dr,w3,e3); pdrawi=t3; vpdrawi=w3; apdrawi=e3; if(i%16)=0)getch(); fclose(fp); getch(); draw1(del,pdraw,vpdraw,apdraw,ic,m); (2)运营成果。杆件3运动参数： The Kinematic Parameters of Pole3No THETA1 T3 W3 E3 deg rad rad/s rad/s/s 1 0.000 -1.632 0.014 -6.232 2 -15.000 -1.632 -0.077 -6.098 3 -30.000 -1.634 -0.163 -5.591 4 -45.000 -1.637 -0.240 -4.731 5 -60.000 -1.641 -0.301 -3.553 6 -75.000 -1.646 -0.343 -2.117 7 -90.000 -1.651 -0.362 -0.501 8 -105.000 -1.656 -0.357 1.192 9 -120.000 -1.661 -0.327 2.84810 -135.000 -1.666 -0.274 4.33911 -150.000 -1.669 -0.201 5.54412 -165.000 -1.671 -0.113 6.35813 -180.000 -1.672 -0.016 6.70314 -195.000 -1.672 0.082 6.54515 -210.000 -1.670 0.174 5.89416 -225.000 -1.667 0.253 4.80717 -240.000 -1.663 0.313 3.38418 -255.000 -1.658 0.351 1.74619 -270.000 -1.653 0.364 0.03020 -285.000 -1.647 0.352 -1.63921 -300.000 -1.642 0.317 -3.14922 -315.000 -1.638 0.261 -4.41523 -330.000 -1.635 0.189 -5.37524 -345.000 -1.632 0.105 -5.98825 -360.000 -1.632 0.014 -6.232运动图形：五.机构动态静力分析5.1六杆铰链式颚式破碎机动态静力分析。（1）求质点7，9，10,11及矿石破碎产生阻力作用点8运动参数；调用bark函数对质点7进行运动分析： 形式参数n1n2n3kr1r2gametwepvpap实 值60750.0r56/20.0twepvpap调用bark函数对质点8进行运动分析： 形式参数n1n2n3kr1r2gametwepvpap实 值60850.0r680.0twepvpap调用bark函数对质点9进行运动分析： 形式参数n1n2n3kr1r2gametwepvpap实 值30940.0r35/20.0twepvpap调用bark函数对质点10进行运动分析： 形式参数n1n2n3kr1r2gametwepvpap实 值401030.0r34/20.0twepvpap调用bark函数对质点11进行运动分析： 形式参数n1n2n3kr1r2gametwepvpap实 值201120.0r23/20.0twepvpap（2）调用rrrf函数对、构件构成RRR杆组进行动态静力分析：形式参数n1n2n3ns1ns2nn1nn2nexfk1实 值635798085形式参数k2pvpaptwefr实 值4pvpaptwefr（3）调用rrrf函数对、构件构成RRR杆组进行动态静力分析：形式参数n1n2n3ns1ns2nn1nn2nexfk1实 值42310113003形式参数k2pvpaptwefr实 值2pvpaptwefr（4）调用barf函数对积极件1进行动态静力分析：形式参数n1ns1nn1k1papefrtb 实 值1121papefr&tb5.2编写主程序并运营。按一定步长，变化积极件位置角度，使其在0-360变化，便可求出机构各运动副反力及作用在积极件上平衡力矩。（1） 主程序。#includegraphics.h#includesubk.c#includeextf.c#includesubf.c#includedraw.cmain()static double p202,vp202,ap202,fr202; static double t10,w10,e10,del,sita1370; static double tb,tb1,fr1,bt1,fr4,bt4,fr6,bt6,we1,we2,we3,we4,we5; static double fe202,tbdraw370,tb1draw370,fr1draw370; static double fr2draw370,sita2370,fr3draw370,sita3370; static int ic; double r12,r23,r34,r35,r56,r68; double pi,dr; int i; FILE *fp; char *m=tb,tb1,fr1,fr2; pi=4.0*atan(1.0); dr=pi/180.0; r12=0.1;r23=1.25;r34=1.0;r35=1.15;r56=1.96;r68=0.6; w1=(-pi*17.0/3.0);e1=0.0;del=1.0; p11=0.0;p12=0.0; p41=0.94;p42=-1.0; p61=-1.0;p62=0.85; 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; printf(nThe Kineto-static Analysis of six-bar Linkbasen); printf(No THETA1 fr6 sita6 tb tb1n); printf( deg N radian N.m N.m ); if(fp=fopen(file_hhhh.txt,w)=NULL) printf(Cant open this file.n); exit(0); fprintf(fp,nThe Kineto-static Analysis of six-bar Linkbasen); fprintf(fp,No THETA1 fr6 sita6 tb tb1n); fprintf(fp, deg N radian N.m N.m ); 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,4,3,2,3,r23,r34,t,w,e,p,vp,ap); rrrk(1,3,6,5,4,5,r35,r56,t,w,e,p,vp,ap); bark(6,0,7,5,0.0,r56/2,0.0,t,w,e,p,vp,ap); bark(6,0,8,5,0.0,r68,0.0,t,w,e,p,vp,ap); bark(3,0,9,4,0.0,r35/2,0.0,t,w,e,p,vp,ap); bark(4,0,10,3,0.0,r34/2,0.0,t,w,e,p,vp,ap); bark(2,0,11,2,0.0,r23/2,0.0,t,w,e,p,vp,ap); rrrf(6,3,5,7,9,8,0,8,5,4,p,vp,ap,t,w,e,fr); rrrf(4,2,3,10,11,3,0,0,3,2,p,vp,ap,t,w,e,fr); barf(1,1,2,1,p,ap,e,fr,&tb); fr1=sqrt(fr11*fr11+fr12*fr12); bt1=atan2(fr12,fr11); fr4=sqrt(fr41*fr41+fr42*fr42); bt4=atan2(fr42,fr41); fr6=sqrt(fr61*fr61+fr62*fr62); bt6=atan2(fr62,fr61);we1=-(ap11*vp11+(ap12+9.81)*vp12)*sm1-e1*w1*sj1； we2=-(ap111*vp111+(ap112+9.81)*vp112)*sm2-e2*w2*sj2； we3=-(ap101*vp101+(ap102+9.81)*vp102)*sm3-e3*w3*sj3； we4=-(ap91*vp91+(ap92+9.81)*vp92)*sm4-e4*w4*sj4; extf(p,vp,ap,t,w,e,8,fe)； we5=-(ap71*vp71+(ap72+9.81)*vp72)*sm5-e5*w5*sj5+fe81*vp81+fe82*vp82; tb1=-(we1+we2+we3+we4+we5)/w1； printf(n%2d%10.3f%10.3f%10.3f%10.3f%10.3f,i+1,t1/dr,fr6,bt6/dr,tb,tb1)； fprintf(fp,n%2d%10.3f%10.3f%10.3%10.3f%10.3f,i+1,t1/dr,fr6,bt6/dr,tb,tb1); tbdrawi=tb；tb1drawi=tb1; fr1drawi=fr1；sita1i=bt1; fr2drawi=fr4;sita2i=bt4; fr3drawi=fr6；sita3i=bt6; if(i%16)=0)getch(); fclose(fp); getch(); draw2(del,tbdraw,tb1draw,ic,m); draw3(del,sita1,fr1draw,sita2,fr2draw,sita3,fr3draw,ic,m);运营成果： The Kineto-static Analysis of six-bar LinkbaseNo THETA1 fr6 sita6 tb tb1 deg 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 12680.958 130.862 -883.317 -883.317 9 -120.000 150