铰链式颚式破碎机机械原理设计

1、机械原理课程设计说明书题目：铰链式颚式破碎机方案分析班级：姓名：学号：指导教师： 成绩：2011 年 9月26 日目 录一 设计题目 1二已知条件及设计要求32.1已知条件32.2设计要求 3三. 机构的结构分析43.1六杆铰链式破碎机43.2四杆铰链式破碎机 4四. 机构的运动分析44.1六杆铰链式颚式破碎机的运动分析44.2四杆铰链式颚式破碎机的运动分析 7五. 机构的动态静力分析 105.1六杆铰链式颚式破碎机的静力分析 105.2四杆铰链式颚式破碎机的静力分析 16六. 工艺阻力函数及飞轮的转动惯量函数 126.1工艺阻力函数程序 216.2飞轮的转动惯量函数程序21七. 对两种机构的

2、综合评价 21八. 主要的收获和建议 22九. 参考文献22一设计题目：铰链式颚式破碎机方案分析二已知条件及设计要求2.1已知条件J6 IB曲柄骋箱图1.1六杆铰链式破碎机图1.2工艺阻力/dI A图1.3四杆铰链式破碎机图(a)所示为六杆铰链式破碎机方案简图。主轴1的转速为n1 = 170r/min ,各部尺寸为：I01A = 0.1m, lAB = 1.250m, IO3B = 1m, 1BC = 1.15m, IO5C = 1.96m, I1=1m, l2=0.94m, h1=0.85m, h2=1m 。各构件质量和转动惯量分别为： m2 = 500kg, Js2 = 25.5kg?m2

3、, m3 = 200kg, Js3 = 9 kg ?m2, m4 = 200kg, Js4=9kg?m2, m5=900kg, Js5=50kg?m2,构件1的质心位于 O1上，其他构件的质心 均在各杆的中心处。D为矿石破碎阻力作用点，设 LO5D = 0.6m，破碎阻力Q在 颚板5的右极限位置到左极限位置间变化，如图(b)所示，Q力垂直于颚板。图(c)是四杆铰链式颚式破碎机方案简图。主轴1的转速n1=170r/min IO1A =0.04m, IAB = 1.11m, I1=0.95m, h1=2m, IO3B=1.96m ，破碎阻力 Q的变化规 律与六杆铰链式破碎机相同，Q力垂直于颚板O3

4、B Q力作用点为D,且IO3D = 0.6m。各杆的质量、转动惯量为 m2= 200kg, Js2=9kg?m2 m3= 900kg, Js3=50kg ?m2曲柄1的质心在O1点处，2、3构件的质心在各构件的中心。2.2设计要求试比较两个方案进行综合评价。主要比较以下几方面：1. 进行运动分析，画出颚板的角位移、角速度、角加速度随曲柄转角的变化曲线。2. 进行动态静力分析，比较颚板摆动中心运动副反力的大小及方向变化规 律，曲柄上的平衡力矩大小及方向变化规律。3. 飞轮转动惯量的大小。三机构的结构分析3.1六杆铰链式破碎机六杆铰链式破碎机拆分为机架和主动件,构件组成的RRF杆组。23.2四杆铰

5、链式破碎机构件组成的RRR干组，3RRR杆组四杆铰链式破碎机拆分为机架和主动件，构件组成的O3四机构的运动分析4.1六杆铰链式颚式破碎机的运动分析(1)调用bark函数求主动件中2点的运动参数。见表4.1表4.1形式参数n1n2n3kr1r2gamtwePvpap实值1201r120.00.0twePvpap(2)调用rrrk函数对由构件组成的RRR杆组进行运动分析。见表4.2表4.2形式参数mn1n2n3k1k2r1r2twepvpap实值-124323r23r34twepvpap(3)调用rrrk函数对由构件组成的RRF杆组进行运动分析。见表4.3表4.3形式参数mn1n2n3k1k2r1

6、r2twepvpap实值136545r35,r56twepvpap(4)六杆运动程序：#in elude graphics.h #i nclude subk.c#i nclude draw.cmai n()static double p202,vp202,ap202;static double t10,w10,e10,del;static double draw370,tdraw370,wdraw370,edraw370; static int ic;double r12,r23,r34,r35,r56,r67,gam1;double pi,dr;double r2,vr2,ar2;int i

7、;FILE *fp;r12=0.1; r34=1.0;r23=1.250;r56=1.96;r35=1.15; r67=0.6;gam1=0.0;pi=4.0*ata n(1.0);w1=-170.0*2*pi/60.0;e1=0.0; del=15.0; dr=pi/180.0;p11=0.0;/*try aga in */p12=0.0;p41=0.94;p42=-1.0;p61=-1.0;p6 2=0.85;rad/sprin tf(n The Kin ematic Parameters of Poi nt6n); prin tf(No THETA1 S7 V7 A7n);printf(

8、”degradrad/s/sn);ic=(i nt)(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);drawi=t1/dr;tdrawi=t5;wdrawi=w5;edrawi=e5;if(fp=fope n( file_ww.txt,w)=NULL) _prin tf(Ca nt open this file.

9、/n);exit(0);99609 寸 eb899 L- 000.09?gCXI卜s寸sCXIbC99V000.9 寸?寸CXIL6、寸 60CXI.0卜 99V000.000- eCXI卜寸CXI寸CXI寸 Lb699V000 9GCXICXI=90Soo =9 000.0 000.0CXI卜 9 868CXI卜 00.9CXI卜 9 goLoLgo.9 =9 寸 8900L9 o卜 9 LCXIL寸 99L.0- 899 9C寸寸 OCXICXIO 999in ooo ocm coiioCXI0 9 0 L CXI 寸 CXI CXI CXIIIIin m cxii6LLJ卜6 CXI c

10、o卜 co CN G卜MdpCXI9CDT16 a) pgluod 40 sa)lueedolelu i|j_-MU-(0Mep9 Mep/vV MeprOJp) l Mep01100)6Ad4sosONolpa)6(OHH(8L&一)七uTe p -=TePM-=Tep5Tep-L-=uai& & & & Pcxl&urd&ld4二二 Mep sMeP2 =TePMTep-L-=uvecxlL&scxlL&scxlL&scxlL& pCXI&ur)匕 d)(+oHvo上)04Snax J0.0Z593O Uhsm 0.4B049OAwax 10 P183466图4.1六杆机构颚板角位置、角速

11、度、角加速度随曲柄转角的变化曲线4.2四杆铰链式颚式破碎机的运动分析(1)调用bark函数求主动件中2点的运动参数。见表4.4表4.4形式参数n1n2n3kr1r2gamtwePvpap实值1201r120.00.0twePvpap(2)调用rrrk函数对由构件组成的RRF杆组进行运动分析。见表4.5表4.5形式参数mn1n2n3k1k2r1r2twepvpap实值124323r23r34twepvpap（3）四杆运动程序：#in elude graphics.h #i nclude subk.c#i nclude draw.cmai n()static double p202,vp202,a

12、p202;static double t10,w10,e10,del;static double draw370,tdraw370,wdraw370,edraw370; static int ic;double r12,r34,r23,r45,gam1;double pi,dr;double r2,vr2,ar2;int i;FILE *fp;r12=0.04;r34=1.96;r23=1.11;r45=0.6;gam1=0.0;pi=4.0*ata n(1.0);w1=-170.0*2*pi/60.0; e1=0.0; del=15.0; dr=pi/180.0;p11=0.0;/*try

13、aga in */p12=0.0;p41=-0.95;p42=2.0;rad/sprin tf(n The Kin ematic Parameters of Poi nt6n); prin tf(No THETA1 S3 V3 A3n);printf(”degradrad/s/sn);ic=(i nt)(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);drawi=t1/dr;tdrawi=t3;w

14、drawi=w3;edrawi=e3;if(fp=fope n( file_ww.txt,w)=NULL) prin tf(Ca nt open this file./n);exit(0);COEsECXI co CXICD1卜T cxi1CXI6TT00 m卜 o 00 寸3CO co9龙T80CD6CO9T0TCO0COTTC)1龙T oco m CXICDcoTco CDTIDCOCD寸9COCDCXI m coCDT9 CXICDco sE9 m CDT109T卜99TCO99T00 m9TCO m 9T卜STCXIsT6 a) po00TCN1o0in CXI CXI1o0CD 特i

15、o0in m CXI1o00CNo0in00CXI1o0CD810NTCXIco寸m9卜006oTTCXIco寸TmT9T卜T0060CXICXI CXIeluod 40 sa)lueedolelu i|j_-MU宀-(0Mep9 Mep/vV Mepra)p) l Mep01100)6_(d4sos宀olpa)6(OHH(8L&一)七uTe p -=TePM-=Tep5Tep-L-=uai& & & & Pcxl&urd&ld4二二 Mep sMeP2 =TePMTep-L-=uvecxlL&scxlL&scxlL&scxlL& pCXI&ur)匕 d)(+oHvo上)040,048B140

16、.3643S4如a G,?033Z5图4.1四杆机构颚板角位置、角速度、角加速度随曲柄转角的变化曲线五.机构的动态静力分析5.1六杆铰链式颚式破碎机的静力分析（1）调用bark函数对主动件进行运动分析。见表 4.1。（2）调用rrrk函数对由构件组成的RRR干组进行运动分析。见表4.2 ：（3）调用rrrk函数对由构件组成的RRR干组进行运动分析。见表4.3 ：（4）求各构件的质心11、8、9、10点及矿石破碎阻力作用点7点的运动参 数。见表5.1表5.5。表5.1 7点运动参数形式参数n1n2n3kr1r2gamtwePvpap实值60750.0r670.0twePvpap表5.2 8 点运

17、动参数形式参数n1n2n3kr1r2gamtwePvpap实值20820.0r23/20.0twePvpap表5.3 9点运动参数形式参数n1n2n3kr1r2gamtwepvpap实值40930.0r34/20.0twepvpap表5.4 10 点运动参数形式参数n1n2n3kr1r2gamtwepvpap实值301040.0r35/20.0twepvpap表5.5 11点运动参数形式参数n1n2n3kr1r2gamtwepvpap实值601150.0r56/20.0twepvpap(5)调用rrrf对由杆组成的RRR干组进行静力分析。见表5.6表5.6形式参数n1n2n3ns1ns2nn1

18、nn2n exfk1k2pvpaptwefr实值365101107745pvpaptwefr(6)调用rrrf对由杆组成的RRF杆组进行静力分析。见表5.7表5.7形式参数n1n2n3ns1ns2nn1nn2n exfk1k2pvpaptwefr实值2438903023pvpaptwefr(7)调用barf对主动件进行静力分析。见表5.8表5.8形式参数n1ns1nn1k1papefrtb实值1121papefr&tb六杆受力程序#in elude graphics.h#inelude subk.c#i nclude subf.c#include draw.cmai n()static dou

19、ble p202,vp202,ap202,del;staticdoublet10,w10,e10,tbdraw370,tb1draw370,fr1draw370,sita1370,fr 2draw370,sita2370,fr3draw370,sita3370;static double fr202,fe202;static int ic;double r12,r34,r23,r35,r47,r56,r67;double gam1,gam2,tb;int i;double pi,dr,fr6,bt6,we1,we2,we3,we4,we5,tb1;FILE*fp;char *m=tb,tb1,

20、fr1,fr2;sm1=0.0;sm2=500.0;sm3=200.0;sm4=200.0;sm5=900.0;sj1=0.0;sj2=25.2;sj3=9.0;sj4=9;sj5=50.0;r12=0.1; r34=1.0;r23=1.250;r56=1.96;r35=1.15; r67=0.6;pi=4.0*ata n(1.0);dr=pi/180.0;w1=-170.0*2*pi/60.0;e1=0.0;del=15.0;p11=0.0；p12=0.0;p41=0.94;p42=-1.0;p61=-1.0;p62=0.85;prin tf(n The Kin eto-static An

21、alysis of a Six-bar Lin kasen);prin tf( NO THETA1 FR1 BT TB TB1n);printf( (deg.) (N) (deg.) (N.m)(N.m)n);if(fp=fope n(file.txt,w)=NULL)prin tf(Ca nt open this file./n);exit(0);Lin kasen);(deg.)fprintf(fp,n TheKineto-staticAnalysis of a Six-barfprin tf(fp,NO THETA1 FR1 FR2 BT2 TB TB1n);fprin tf(fp,(d

22、eg.)(N)(deg.)(N)(N.m)(N.m)n);ic=(i nt)(360.0/del);for(i=0;i=ic;i+)t1=(double)(-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,r67,0.0,t,w,e,p,vp,ap); bark(2,0,8,2,0.0,r23/2,0.0,t,w,e,p,vp,ap); b

23、ark(4,0,9,3,0.0,r34/2,0.0,t,w,e,p,vp,ap); bark(3,0,10,4,0.0,r35/2,0.0,t,w,e,p,vp,ap); bark(6,0,11,5,0.0,r56/2,0.0,t,w,e,p,vp,ap); rrrf(3,6,5,10,11,0,7,7,4,5,p,vp,ap,t,w,e,fr); rrrf(2,4,3,8,9,0,3,0,2,3,p,vp,ap,t,w,e,fr);barf(1,1,2,1,p, ap,e,fr, &tb);fr6=sqrt(fr 1*fr 1+fr2*fr 2);bt6=ata n2(fr62,fr 1);

24、we1=-(ap11*vp11+(ap12+9.81)*vp12)*sm1-e1*w1*sj1 ;we2=-(ap81*vp81+( ap82+9.81)*vp82)*sm2-e2*w2*sj2 ；we3=-(ap91*vp91+( ap92+9.81)*vp92)*sm3-e3*w3*sj3 ；we4=-(ap101*vp101+( ap102+9.81)*vp102)*sm4-e4*w4 *sj4;extf(p,vp,ap,t,w,e,7,fe);we5=-(ap111*vp111+( ap112+9.81)*vp112)*s m 5-e5*w 5 *sj5+fe71*vp71+fe72*

25、vp72;%10.3ftb1=-(we1+we2+we3+we4+we5)/w1;prin tf(%2d%10.3f%10.3f%10.3f%10.3fn,i+1,t1/dr,fr6,bt6/dr,tb,tb1);fprin tf(fp,%2d %10.3f %10.3f %10.3f %10.3f %10.3f,i+1,t1/dr,fr6,b t6/dr,tb,tb1);tbdrawi=tb;tb1drawi=tb1;fr1drawi=fr6;sita1i=bt6;fr2drawi=fr6;sita2i=bt6;fr3drawi=fr6;sita3i=bt6;if(i%16)=0)getch

26、();fclose(fp);getch();draw2(del,tbdraw,tb1draw,ic,m); draw3(del,sita1,fr1draw,sita2,fr2draw,sita3,fr3draw,ic,m);extf(p,vp,ap,t,w,e, nexf,fe)double p202,vp20, ap202,t10,w10,e10,fe202;int n exf;double pi,dr;pi=4.0*ata n(1.0);dr=pi/180.0;if(t1/dr=-273.0)fen exf1=-85000.0/(-183.0*dr)*(t1+90*dr)*si n(t5)

27、;fen exf2=85000.0/(-183.0*dr)*(t1+90*dr)*cos(t 5);elsefen exf1=0.0;fe nexf2=0.0;运行结果:The Ki neto-static An alysis of a Six-bar Lin kaseNOTHETA1FR4BT2 TBTB1(deg.)(N)(deg.) (N.m)(N.m)10.000 9904.58077.690534.454534.4542-15.00010248.08682.6701038.4481038.4483-30.00010522.85289.5761434.9051434.9054-45.0

28、0010757.31497.3291548.0671548.0675-60.00010967.175104.3391271.1131271.1136-75.00011112.158109.009644.146644.1467-90.00011132.496110.330-144.853-144.8538-105.00013340.278128.365-883.912-883.9129-120.00016136.071138.884-1407.4681 -1407.46810-135.00019120.284145.429-1625.253 -1625.25311-150.00022280.55

29、7150.044-1557.401 -1557.40112-165.00025724.638153.629-1290.413 -1290.41313-180.00029550.567156.506-928.440-928.44014-195.00033789.984158.781-562.443-562.44315-210.00038396.542160.527-258.404-258.40416-225.00043263.326161.841-55.542-55.54217-240.00048254.151162.84433.74733.74718-255.00053236.591163.6

30、5324.76824.76819-270.00058109.329164.367-39.786-39.78620-285.0008481.82378.617-205.535-205.53521-300.0008583.46577.292-339.048-339.04822-315.0008793.29375.658-361.780-361.78023-330.0009113.15874.602-227.795-227.79524-345.0009506.21075.05980.792(30.79225-360.0009904.58077.690534.454534.454图5.2六杆机构颚板摆

31、动中心运动副反力的大小及方向变化规律图5.1六杆机构曲柄上的平衡力矩的变化规律5.2四杆铰链式颚式破碎机的静力分析(1) 调用bark函数对主动件进行运动分析。见表 4.4。(2) 调用rrrk函数对由构件组成的RRR干组进行运动分析。见表4.5(3) 求各构件的质心6、7点及矿石破碎阻力作用点5点的运动参数。 见表5.9表5.11。表5.9 5点运动参数形式参数n1n2n3kr1r2gamtwepvpap实值40530.0r450.0twepvpap表5.10 6 点运动参数形式参数n1n2n3kr1r2gamtwepvpap实值20620.0r23/20.0twepvpap表5.11 7点

32、运动参数形式参数n1n2n3kr1r2gamtwepvpap实值40730.0r34/20.0twepvpap(4)调用rrrf对由杆组成的RRF杆组进行静力分析。见表5.12表 5.12形式参数n1n2n3ns1ns2nn1nn2n exfk1k2pvpaptwefr实值2436705523pvpaptwefr(5)调用barf对主动件进行静力分析。见表5.13表 5.13形式参数n1ns1nn1k1papefrtb实值1121papefr&tb四杆力程序：#in elude graphics.h#i nclude subk.c#in clude subf.c#include draw.cm

33、ai n()static double p202,vp202,ap202,del;staticdoublet10,w10,e10,tbdraw370,tb1draw370,fr1draw370,sita1370,fr 2draw370,sita2370,fr3draw370,sita3370;static double fr202,fe202;static int ic;double r12,r34,r23,r45,r47,r36,f;double gam1,gam2,tb;int i;double pi,dr,fr4,bt4,we1,we2,we3,tb1;FILE*fp;char *m=t

34、b,tb1,fr1,fr2; sm1=0.0;sm2=200.0;sm3=900.0;sj1=0.0;sj2=9.0;sj3=50.0;r12=0.04;r34=1.96;r23=1.11;r45=0.6;gam仁 0.0; gam2=0.0;pi=4.0*ata n(1.0);dr=pi/180.0;w1=-170.0*2*pi/60.0;e1=0.0;del=15.0;p11=0.0;/*try aga in */p12=0.0;p41=-0.95;p42=2.0;printf(nThe Kineto-static Analysis of a Six-bar Linkasen);print

35、f( NO THETA1 FR1 BT1 TB TB1f n);prin tf(deg.) (N) (deg.) (N.m)(N.m)(N) n);if(fp=fope n(file.txt,w)=NULL)prin tf(Ca nt open this file./n);exit(0);fprintf(fp,nThe Kineto-staticAnalysisof a Six-barLin kasen);fprin tf(fp,NO THETA1 FR1 FR2 BT2 TB TB1n);fprin tf(fp,(deg.)(N) (deg.) (N) (deg.) (N.m)(N.m)n)

36、;ic=(i nt)(360.0/del);for(i=0;i=-183.0&t1/dr=-2.0)fen exf1=-85000.0/(-181.0*dr)*(t1+2*dr)*si n( t3);fen exf2=85000.0/(-181.0*dr)*(t1+2*dr)*cos(t3);elsefen exf1=0.0;fe nexf2=0.0; 运行结果:The Ki neto-static An alysis of a Six-bar Lin kaseNO THETA1fr1sita1radia n Nfr4 sita4 radia n N.tb tb1 m N.mdegN10.00

37、010261.396102.68547.54647.5462-15.00012216.646121.678-19.102-19.1023-30.00015476.665135.641-119.595-119.5954-45.00019206.225144.108-237.935-237.9355-60.00023110.501149.501-356.889-356.8896-75.00027072.942153.130-460.689-460.6897-90.00031057.291155.699-536.638-536.6388-105.00035068.160157.604-575.763

38、-575.7639-120.00039133.233159.082-572.422-572.42210-135.00043292.628160.280-523.140-523.14011-150.00047590.352161.289-425.335-425.33512-165.00052066.382162.167-276.749-276.74913-180.00056749.461162.956-76.249-76.24914-195.0009702.56670.652-75.371-75.37115-210.0009538.37372.140-106.622-106.62216-225.0009408.54074.864-112.809-112.80917-240.0009336.37178.581-90.324-90.32418-255.0009335.57282.915-43.909-43.90919-270.0009405.05887.43915.06915.06920-285.0009529.82891.75472.27972.27921-300.0009687.