机械原理课程设计说明书题目：铰链式颚式破碎机方案分析目录一设计题目 (1)二已知条件及设计要求 (1)2.1已知条件 (1)2.2设计要求 (2)三. 机构的结构分析 (2)3.1六杆铰链式破碎机 (2)3.2四杆铰链式破碎机 (2)四. 机构的运动分析 (2)4.1六杆铰链式颚式破碎机的运动分析 (2)4.2四杆铰链式颚式破碎机的运动分析 (4)五.机构的动态静力分析 (6)5.1六杆铰链式颚式破碎机的静力分析 (6)5.2四杆铰链式颚式破碎机的静力分析 (12)六. 工艺阻力函数及飞轮的转动惯量函数 (16)6.1工艺阻力函数程序 (16)6.2飞轮的转动惯量函数程序 (17)七 .对两种机构的综合评价 (19)八 . 主要的收获和建议 (20)九 . 参考文献 (20)一设计题目铰链式颚式破碎机方案分析二已知条件及设计要求2.1已知条件图(a)所示为六杆铰链式破碎机方案简图。
主轴1的转速为n1 = 170r/min，各部尺寸为：l O1A= 0.1m,l AB = 1.250m, l O3B = 1m, l BC = 1.15m, l O5C = 1.96m, l1=1m, l=0.94m, h1=0.85m, h2=1m。
各构件质量和转动惯量分别为：m2 = 500kg, J s2 = 225.5kg·m2, m3 = 200kg, J s3 = 9kg·m2, m4 = 200kg, J s4 = 9kg·m2, m5=900kg, J=50kg·m2, 构件1的质心位于O1上，其他构件的质心均在各杆的中心处。
D s5为矿石破碎阻力作用点，设L O5D = 0.6m，破碎阻力Q在颚板5的右极限位置到左极限位置间变化，如图(b)所示，Q力垂直于颚板。
图(c)是四杆铰链式颚式破碎机方案简图。
主轴1 的转速n1=170r/min。
l O1A = 0.04m, l AB = 1.11m, l1=0.95m, h1=2m, l O3B=1.96m，破碎阻力Q的变化规律与六杆铰链式破碎机相同，Q力垂直于颚板O3B，Q力作用点为D，且l O3D = 0.6m。
各= 200kg, J s2=9kg·m2，m3 = 900kg, J s3=50kg·m2。
曲杆的质量、转动惯量为m2柄1的质心在O1点处，2、3构件的质心在各构件的中心。
(a) 六杆铰链式破碎机(b) 工艺阻力(c) 四杆铰链式破碎机2.2设计要求试比较两个方案进行综合评价。
主要比较以下几方面：1. 进行运动分析，画出颚板的角位移、角速度、角加速度随曲柄转角的变化曲线。
2. 进行动态静力分析，比较颚板摆动中心运动副反力的大小及方向变化规律，曲柄上的平衡力矩大小及方向变化规律。
3. 飞轮转动惯量的大小。
三. 机构的结构分析3.1六杆铰链式破碎机3.2四杆铰链式破碎机四. 机构的运动分析4.1六杆铰链式颚式破碎机的运动分析（1）调用bark函数求2点的运动参数（2）调用rrrk函数求3点的运动参数（3）调用rrrk函数求5点的运动参数（4）程序:对构件5的运动轨迹分析#include"graphics.h"#include "subk.c"#include "draw.c"main(){static double p[20][2],vp[20][2],ap[20][2],del;static double t[10],w[10],e[10],tdraw[370],wdraw[370],edraw[370];static int ic;double r12,r23,r34,r35,r56;double pi ,dr;double r2,vr2,ar2;int i; FILE*fp;r12=0.1,r23=1.25,r34=1.0,r35=1.15,r56=1.96;w[1]=-17.8;del=15;pi=4.0*atan(1.0);dr=pi/180.0;p[1][1]=0.0;p[1][2]=0.0;p[4][1]=0.94;p[4][2]=-1.0;p[6][1]=-1.0;p[6][2]=0.85;printf("\n The Kinematic Parameters of Point 5\n");printf("No THETA1 t w e\n");printf(" deg rad rad/s rad/s/s\n");if((fp=fopen("file20133098.txt","w"))==NULL){printf("Can't open this file.\n");exit(0);}fprintf(fp,"\n The kinematic parameters of point 10\n");fprintf(fp,"No THETA1 t w e\n");fprintf(fp," deg rad rad/s rad/s/s");ic=(int)(360.0/del);for(i=0;i<=ic;i++){t[1]=(-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);printf("\n%2d %12.3f %12.3f %12.3f %12.3f",i+1,t[1]/dr,t[5],w[5],e[5]);fprintf(fp,"\n%2d%12.3f%12.3f%12.3f%12.3f",i+1,t[1]/dr,t[5],w[5],e[5]);tdraw[i]=t[5];wdraw[i]=w[5];edraw[i]=e[5];if((i%16)==0){getch();};}fclose(fp);getch();draw1(del,tdraw,wdraw,edraw,ic);}（5）数据:随主动件1变化的运动参数The kinematic parameters of point 5No THETA1 t5 w5 e5deg rad rad/s rad/s/s1 0.000 -1.658 0.346 3.9552 -15.000 -1.653 0.392 2.0023 -30.000 -1.647 0.400 -0.9324 -45.000 -1.641 0.362 -4.3545 -60.000 -1.637 0.274 -7.5046 -75.000 -1.633 0.146 -9.6107 -90.000 -1.632 -0.001 -10.1818 -105.000 -1.633 -0.145 -9.1629 -120.000 -1.637 -0.265 -6.90210 -135.000 -1.641 -0.345 -3.98011 -150.000 -1.646 -0.382 -1.00812 -165.000 -1.652 -0.377 1.51813 -180.000 -1.657 -0.341 3.29614 -195.000 -1.662 -0.284 4.23615 -210.000 -1.666 -0.220 4.43516 -225.000 -1.668 -0.156 4.12017 -240.000 -1.670 -0.10 3.58318 -255.000 -1.671 -0.051 3.10519 -270.000 -1.672 -0.007 2.89720 -285.000 -1.672 0.036 3.06321 -300.000 -1.671 0.085 3.57022 -315.000 -1.669 0.142 4.24623 -330.000 -1.667 0.209 4.79024 -345.000 -1.663 0.281 4.81625 -360.000 -1.658 0.346 3.955（6）线图:构件5角位置，角速度，角加速度线图六杆机构颚板角位置、角速度、角加速度随曲柄转角的变化曲线4.2四杆铰链式颚式破碎机的运动分析（1）调用bark函数求2点的运动参数（2）调用rrrk函数求3点的运动参数（3）程序：对构件3的运动轨迹分析#include"graphics.h"#include "subk.c"#include "draw.c"main(){static double p[20][2],vp[20][2],ap[20][2],del;static double t[10],w[10],e[10],tdraw[370],wdraw[370],edraw[370];static int ic;double r12,r23,r34;double pi,dr;double r2,vr2,ar2;int i;FILE*fp; r12=0.04,r23=1.11,r34=1.96;w[1]=-17.8;del=15;pi=4.0*atan(1.0);dr=pi/180.0;p[1][1]=0.0;p[1][2]=0.0;p[4][1]=-0.95;p[4][2]=2.0;printf("\n The Kinematic Parameters of Point 5\n");printf("No THETA1 t3 w3 e3\n");printf(" deg rad rad/s rad/s/s\n");if((fp=fopen("filel20133098.txt","w"))==NULL){printf("Can't open this file.\n");exit(0);}fprintf(fp,"\n The kinematic parameters of point 10\n");fprintf(fp,"No THETA1 t3 w3 e3\n");fprintf(fp," deg rad rad/s rad/s/s");ic=(int)(360.0/del);for(i=0;i<=ic;i++){t[1]=(-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);printf("\n%2d %12.3f %12.3f %12.3f",i+1,t[1]/dr,t[3],w[3],e[3]);fprintf(fp,"\n%2d%12.3f%12.3f%12.3f%12.3f",i+1,t[1]/dr,t[3],w[3],e[3]);wdraw[i]=t[1]/dr;tdraw[i]=t[3];wdraw[i]=w[3];edraw[i]=e[3];if((i%16)==0){getch();};}fclose(fp);getch();draw1(del,tdraw,wdraw,edraw,ic);}（4）数据：随主动件1变化的运动参数The kinematic parameters of point 3No THETA1 t3 w3 e3deg rad rad/s rad/s/s1 0.000 -1.632 0.014 -6.2302 -15.000 -1.632 -0.077 -6.0973 -30.000 -1.634 -0.163 -5.5904 -45.000 -1.637 -0.240 -4.7305 -60.000 -1.641 -0.301 -3.5536 -75.000 -1.646 -0.343 -2.1177 -90.000 -1.651 -0.362 -0.5018 -105.000 -1.656 -0.357 1.1929 -120.000 -1.661 -0.327 2.84710 -135.000 -1.666 -0.274 4.33811 -150.000 -1.669 -0.201 5.54312 -165.000 -1.671 -0.113 6.35613 -180.000 -1.672 -0.016 6.70214 -195.000 -1.672 0.082 6.54315 -210.000 -1.670 0.174 5.89216 -225.000 -1.667 0.253 4.80617 -240.000 -1.663 0.313 3.38318 -255.000 -1.658 0.351 1.74619 -270.000 -1.653 0.364 0.03020 -285.000 -1.647 0.352 -1.63821 -300.000 -1.642 0.317 -3.14822 -315.000 -1.638 0.261 -4.41423 -330.000 -1.635 0.189 -5.37324 -345.000 -1.632 0.105 -5.98625 -360.000 -1.632 0.014 -6.230（6）线图:3点水平位移，速度，加速度线图四杆机构颚板角位置、角速度、角加速度随曲柄转角的变化曲线五.机构的动态静力分析5.1六杆铰链式颚式破碎机的静力分析（1）、（2）、（3）步同运动分析1、2、3（4）调用bark函数求9的运动参数（5）调用bark函数求10的运动参数（6）调用bark函数求8的运动参数（7）调用bark函数求7的运动参数（8）调用bark函数求11的运动参数（9）调用rrrf对4、5杆件组成的rrr杆组进行静力分析（10）调用rrrf对2、3杆组成的rrr杆组进行静力分析（11）调用barf对主动件1进行静力分析（12）程序：对质心的运动分析，对固定铰链的静态动力分析，主动反力偶#include"graphics.h"#include"subk.c"#include"subf.c"#include"draw.c"main(){static double p[20][2],vp[20][2],ap[20][2],del;t[10],w[10],e[10],static double tbdraw[370],tb1draw[370];fr3draw[370],sita1[370],fr1draw[370];static double sita2[370],fr2draw[370],sita3[370],we1,we2,we3,we4,we5;static double fr[20][2],fe[20][2],tb,tb1,fr1,bt1,fr4,bt4,fr6,bt6;static int ic;double pi,dr;int i;FILE *fp;char *m[]={"tb","tb1","fr1","fr4","fr6"};sm[1]=0.0;sm[2]=500;sm[3]=200;sm[4]=200;sm[5]=900;del=15;sj[1]=0.0;sj[2]=25.5;sj[3]=9;sj[4]=9;sj[5]=50;t[1]=0.0;w[1]=-17.8;e[1]=0.0;t[4]=0.0;t[6]=90.0;w[6]=0.0;e[6]=0.0;pi=4.0*atan(1.0);dr=pi/180.0;t[6]=90.0*dr;p[1][1]=0.0;p[1][2]=0.0;p[4][1]=0.94;p[4][2]=-1.0;p[6][1]=-1.0;p[6][2]=0.85;printf("\n The Kinet0-static Analysis of a six-bar Linkase.\n");printf("No HETA1 fr1 sita1 fr4 sita2 fr7 sita3 tb tb1\n");printf(" deg N radian N radian N radian N.m N.m\n");if((fp=fopen("filel20133098","w"))==NULL){printf("Can't open this file.\n");exit(0);}fprintf(fp,"\n The Kinet0-static Analysis of a six-bar Linkase.\n");fprintf(fp,"No HETA1 fr1 sita1 fr4 sita2fr7 sita3 tb tb1\n");fprintf(fp," deg N radian N radianN radian N.m N.m\n");ic=(int)(360.0/del);for(i=0;i<=ic;i++){t[1]=(-i*del)*dr;bark(1,2,0,1,0.1,0.0,0.0,t,w,e,p,vp,ap);rrrk(-1,2,4,3,2,3,1.25,1.0,t,w,e,p,vp,ap);rrrk(1,3,6,5,4,5,1.15,1.96,t,w,e,p,vp,ap);bark(2,0,9,2,0.0,0.625,0.0,t,w,e,p,vp,ap);bark(4,0,10,3,0.0,0.5,0.0,t,w,e,p,vp,ap);bark(3,0,8,4,0.0,0.575,0.0,t,w,e,p,vp,ap);bark(6,0,7,5,0.0,0.98,0.0,t,w,e,p,vp,ap);bark(6,0,11,5,0.0,0.6,0.0,t,w,e,p,vp,ap);rrrf(3,6,5,8,7,0,11,11,4,5,p,vp,ap,t,w,e,fr);rrrf(2,4,3,9,10,0,3,0,2,3,p,vp,ap,t,w,e,fr);barf(1,1,2,1,p,ap,e,fr,&tb);fr1=sqrt(fr[1][1]*fr[1][1]+fr[1][2]*fr[1][2]);bt1=atan2(fr[1][2],fr[1][1]);fr4=sqrt(fr[4][1]*fr[4][1]+fr[4][2]*fr[4][2]);bt4=atan2(fr[4][2],fr[4][1]);fr6=sqrt(fr[6][1]*fr[6][1]+fr[6][2]*fr[6][2]);bt6=atan2(fr[6][2],fr[6][1]);we1=-(ap[1][1]*vp[1][1]+(ap[1][2]+9.81)*vp[1][2])*sm[1]-e[1]*w[1]*sj[1];we2=-(ap[9][1]*vp[9][1]+(ap[9][2]+9.81)*vp[9][2])*sm[2]-e[2]*w[2]*sj[2];we3=-(ap[10][1]*vp[10][1]+(ap[10][2]+9.81)*vp[10][2])*sm[3]-e[3]*w[3]*sj[3];we4=-(ap[8][1]*vp[8][1]+(ap[8][2]+9.81)*vp[8][2])*sm[4]-e[4]*w[4]*sj[4];extf(p,vp,ap,t,w,e,11,fe);we5=-(ap[7][1]*vp[7][1]+(ap[7][2]+9.81)*vp[7][2])*sm[5]-e[5]*w[5]*sj[5]+fe[11][1]*vp[11][1]+fe[11][2]*vp[11][2];tb1=-(we1+we2+we3+we4+we5)/w[1];printf("\n%2d %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f",i+1,t[1]/dr,fr1,bt1/dr,fr4,bt4/dr,fr6/dr,bt6/dr,tb,tb1);fprintf(fp,"\n%2d %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f",i+1,t[1]/dr,fr1,bt1/dr,fr4,bt4/dr,fr6/dr,bt6/dr,tb,tb1);tbdraw[i]=tb;tb1draw[i]=tb1;fr1draw[i]=fr6;sita1[i]=bt6;fr2draw[i]=fr6;sita2[i]=bt6;fr3draw[i]=fr6;sita3[i]=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);}extf(p,vp,ap,t,w,e,nexf,fe)double p[20][2],vp[20][2],ap[20][2],t[10],w[10],e[10],fe[20][2];int nexf;{ double pi,dr;pi=4.0*atan(1.0);dr=pi/180.0;if(w[5]<0){fe[nexf][1]=(-t[1]/dr-90.0)*(85000.0/182.0)*cos(-t[5]-pi/2);fe[nexf][2]=-(-t[1]/dr-90.0)*(85000.0/182.0)*sin(-t[5]-pi/2);}else{fe[nexf][1]=0;fe[nexf][2]=0;}}（13）数据：6点固定铰链力矢；主动件平衡力偶The Kineto-static Analysis of a Six-bar LinkaseNO THETA1 fr6 sita6 tb tb1deg N deg N.m N.m1 0.000 9904.580 77.690 534.273 534.2732 -5.000 10027.035 79.096 703.325 703.3253 -10.000 10141.706 80.761 873.278 873.2784 -15.000 10248.086 82.670 1038.104 1038.1045 -20.000 10346.409 84.798 1191.197 1191.1976 -25.000 10437.556 87.113 1325.633 1325.6337 -30.000 10522.852 89.576 1434.513 1434.5138 -35.000 10603.782 92.138 1511.359 1511.3599 -40.000 10681.663 94.743 1550.551 1550.55110 -45.000 10757.314 97.329 1547.760 1547.76011 -50.000 10830.805 99.833 1500.331 1500.33112 -55.000 10901.317 102.191 1407.585 1407.58513 -60.000 10967.175 104.339 1270.987 1270.98714 -65.000 11026.038 106.222 1094.172 1094.17215 -70.000 11075.232 107.792 882.793 882.79316 -75.000 11112.158 109.009 644.228 644.22817 -80.000 11134.700 109.847 387.143 387.14318 -85.000 11141.565 110.289 120.981 120.98119 -90.000 11132.496 110.330 -144.608 -144.60820 -95.000 11754.579 117.451 -404.323 -404.32321 -100.000 12510.231 123.454 -653.284 -653.28422 -105.000 13356.757 128.454 -883.773 -883.77323 -110.000 14262.754 132.611 -1089.373 -1089.37324 -115.000 15206.996 136.082 -1265.169 -1265.16925 -120.000 16176.401 139.007 -1407.830 -1407.83026 -125.000 17164.019 141.503 -1515.584 -1515.58427 -130.000 18167.323 143.663 -1588.106 -1588.10628 -135.000 19186.867 145.561 -1626.344 -1626.34429 -140.000 20225.248 147.252 -1632.299 -1632.29930 -145.000 21286.309 148.779 -1608.796 -1608.79631 -150.000 22374.519 150.172 -1559.245 -1559.24532 -155.000 23494.484 151.455 -1487.428 -1487.42833 -160.000 24650.573 152.643 -1397.304 -1397.30434 -165.000 25846.624 153.746 -1292.854 -1292.85435 -170.000 27085.723 154.772 -1177.951 -1177.95136 -175.000 28370.060 155.725 -1056.267 -1056.26737 -180.000 29700.830 156.609 -931.213 -931.21338 -185.000 31078.208 157.426 -805.893 -805.89339 -190.000 32501.350 158.178 -683.087 -683.08740 -195.000 33968.451 158.869 -565.243 -565.24341 -200.000 35476.829 159.501 -454.473 -454.47342 -205.000 37023.029 160.078 -352.559 -352.55943 -210.000 38602.952 160.602 -260.954 -260.95444 -215.000 40211.997 161.078 -180.787 -180.78745 -220.000 41845.198 161.511 -112.856 -112.85646 -225.000 43497.372 161.905 -57.628 -57.62847 -230.000 45163.259 162.265 -15.234 -15.23448 -235.000 46837.658 162.595 14.542 14.54249 -240.000 48515.553 162.900 32.265 32.26550 -245.000 50192.226 163.184 38.848 38.84851 -250.000 51863.358 163.450 35.553 35.55352 -255.000 53525.123 163.703 23.978 23.97853 -260.000 55174.261 163.946 6.042 6.04254 -265.000 56808.141 164.181 -16.045 -16.04555 -270.000 58424.813 164.412 -39.809 -39.80956 -275.000 8464.780 79.032 -83.824 -83.82457 -280.000 8468.680 78.880 -146.286 -146.28658 -285.000 8481.823 78.617 -205.306 -205.30659 -290.000 8504.820 78.253 -258.572 -258.57260 -295.000 8538.430 77.804 -303.789 -303.78961 -300.000 8583.465 77.292 -338.729 -338.72962 -305.000 8640.661 76.742 -361.265 -361.26563 -310.000 8710.545 76.186 -369.428 -369.42864 -315.000 8793.293 75.658 -361.459 -361.45965 -320.000 8888.623 75.193 -335.863 -335.86366 -325.000 8995.710 74.829 -291.484 -291.48467 -330.000 9113.158 74.602 -227.576 -227.57668 -335.000 9239.033 74.545 -143.889 -143.88969 -340.000 9370.949 74.689 -40.753 -40.75370 -345.000 9506.210 75.059 80.824 80.82471 -350.000 9642.002 75.674 219.074 219.07472 -355.000 9775.598 76.549 371.389 371.38973 -360.000 9904.580 77.690 534.273 534.273 （14）线图：6点固定铰链力矢；主动件平衡力偶3 1六杆机构曲柄上的平衡力矩的变化规律六杆机构颚板摆动中心运动副反力的大小及方向变化规律5.2四杆铰链式颚式破碎机的静力分析（1）、（2）步同运动分析1、2（3）、调用bark函数对7点进行运动分析（4）、调用bark函数对5点进行运动分析（5）、调用bark函数对6点进行运动分析（6）、调用rrrf函数对2、3杆进行静力分析（7）、调用barf函数对主动件1进行静力分析（8）、程序：对质心的运动分析，对固定铰链的静态动力分析，主动反力偶#include"graphics.h"#include"subk.c"#include"subf.c"#include"draw.c"main(){static double p[20][2],vp[20][2],ap[20][2],t[10],w[10],e[10],del;static double tbdraw[370],tb1draw[370],sita1[370],fr1draw[370];static double sita2[370],fr2draw[370],sita3[370],fr3draw[370];static double fr[20][2],fe[20][2],tb,tb1,fr1,bt1,fr4,bt4,we1,we2,we3;static int ic;double pi,dr;int i;FILE *fp;char *m[]={"tb","tb1","fr1","fr4"};sm[1]=0.0;sm[2]=200.0;sm[3]=900.0;del=15;sj[1]=0.0;sj[2]=9.0;sj[3]=50.0;t[1]=0.0;w[1]=-17.8;e[1]=0.0;pi=4.0*atan(1.0);dr=pi/180.0;t[4]=90.0*dr;p[1][1]=0.0;p[1][2]=0.0;p[4][1]=-0.95;p[4][2]=2.0;printf("\n The Kinet0-static Analysis of a four-bar Linkase.\n");printf("No HETA1 fr1 sita1 fr4 sita2 tb tb1\n");printf(" deg N radian N radian N.m N.m\n");if((fp=fopen("filel20133098.txt","w"))==NULL){printf("Can't open this file.\n");exit(0);}fprintf(fp,"\n The Kinet0-static Analysis of a four-bar Linkase.\n");fprintf(fp,"No HETA1 fr1 sita1 fr4 sita2 tb tb1\n");fprintf(fp," deg N radian N radian N.m N.m\n");ic=(int)(360.0/del);for(i=0;i<=ic;i++){ t[1]=(-i*del)*dr;bark(1,2,0,1,0.04,0.0,0.0,t,w,e,p,vp,ap);rrrk(1,2,4,3,2,3,1.11,1.96,t,w,e,p,vp,ap);bark(2,0,7,2,0.0,0.555,0.0,t,w,e,p,vp,ap);bark(4,0,5,3,0.0,0.98,0.0,t,w,e,p,vp,ap);bark(4,0,6,3,0.0,0.6,0.0,t,w,e,p,vp,ap);rrrf(2,4,3,7,5,0,6,6,2,3,p,vp,ap,t,w,e,fr);barf(1,1,2,1,p,ap,e,fr,&tb);fr1=sqrt(fr[1][1]*fr[1][1]+fr[1][2]*fr[1][2]);bt1=atan2(fr[1][2],fr[1][1]);fr4=sqrt(fr[4][1]*fr[4][1]+fr[4][2]*fr[4][2]);bt4=atan2(fr[4][2],fr[4][1]);we1=-(ap[1][1]*vp[1][1]+(ap[1][2]+9.81)*vp[1][2])*sm[1]-e[1]*w[1]*sj[1];we2=-(ap[7][1]*vp[7][1]+(ap[7][2]+9.81)*vp[7][2])*sm[2]-e[2]*w[2]*sj[2];extf(p,vp,ap,t,w,e,6,fe);we3=-(ap[5][1]*vp[5][1]+(ap[5][2]+9.81)*vp[5][2])*sm[3]-e[3]*w[3]*sj[3]+fe[6][1]*vp[6][1]+fe[6][2]*vp[6][2];tb1=-(we1+we2+we3)/w[1];printf("\n%2d %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f",i+1,t[1]/dr,fr1,bt1/dr,fr4,bt4/dr,tb,tb1);fprintf(fp,"\n%2d %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f %10.3f",i+1,t[1]/dr,fr1,bt1/dr,fr4,bt4/dr,tb,tb1);tbdraw[i]=tb;tb1draw[i]=tb1;fr1draw[i]=fr4;sita1[i]=bt4;fr2draw[i]=fr4;sita2[i]=bt4;fr3draw[i]=fr4;sita3[i]=bt4; if((i%16)==0){getch();}}fclose(fp);getch();draw2(del,tbdraw,tb1draw,ic);draw3(del,sita1,fr1draw,sita2,fr2draw,sita3,fr3draw,ic);}extf(p,vp,ap,t,w,e,nexf,fe)double p[20][2],vp[20][2],ap[20][2],t[10],w[10],e[10],fe[20][2];int nexf;{ double pi,dr;pi=4.0*atan(1.0);dr=pi/180.0; if(w[3]<0){ fe[nexf][1]=(-t[1]/dr-2.0)*(85000.0/181.0)*cos(-t[3]-pi/2);fe[nexf][2]=-(-t[1]/dr-2.0)*(85000.0/181.0)*sin(-t[3]-pi/2);}else{fe[nexf][1]=0;fe[nexf][2]=0;}}}（9）、数据：对质心的运动分析，对固定铰链的静态动力分析，主动反力偶The Kinet0-static Analysis of a four-bar Linkase.No HETA1 fr4 sita2 tb tb1deg N radian N.m N.m1 0.000 10261.190 102.681 47.544 47.5442 -15.000 12216.224 121.676 -19.091 -19.0913 -30.000 15476.149 135.639 -119.576 -119.5764 -45.000 19205.725 144.108 -237.910 -237.9105 -60.000 23110.088 149.501 -356.865 -356.8656 -75.000 27072.663 153.130 -460.673 -460.6737 -90.000 31057.179 155.699 -536.634 -536.6348 -105.000 35068.233 157.604 -575.772 -575.7729 -120.000 39133.494 159.083 -572.442 -572.44210 -135.000 43293.064 160.281 -523.165 -523.16511 -150.000 47590.933 161.289 -425.359 -425.35912 -165.000 52067.066 162.167 -276.765 -276.76513 -180.000 56750.195 162.956 -76.251 -76.25114 -195.000 9702.491 70.657 -75.360 -75.36015 -210.000 9538.359 72.144 -106.600 -106.60016 -225.000 9408.589 74.868 -112.783 -112.78317 -240.000 9336.469 78.583 -90.301 -90.30118 -255.000 9335.692 82.916 -43.895 -43.89519 -270.000 9405.167 87.439 15.070 15.07020 -285.000 9529.897 91.752 72.267 72.26721 -300.000 9687.095 95.552 114.526 114.52622 -315.000 9853.651 98.636 133.067 133.06723 -330.000 10011.380 100.886 125.127 125.12724 -345.000 10148.934 102.242 93.949 93.94925 -360.000 10261.190 102.681 47.544 47.544（3）、线图：对质心的运动分析，对固定铰链的静态动力分析，主动反力偶四杆机构曲柄上的平衡力矩的变化规律四杆机构颚板摆动中心运动副反力的大小及方向变化规律六. 工艺阻力函数及飞轮的转动惯量函数6.1工艺阻力函数程序（1）六杆工艺阻力函数程序#include"math.h"extf(p,vp,ap,t,w,e,nexf,fe)double p[20][2],vp[20][2],ap[20][2],t[10],w[10],e[10],fe[20][2]; int nexf;{double pi,dr;pi=4.0*atan(1.0);dr=pi/180.0;if(w[5]<0){fe[nexf][1]=(-t[1]/dr-90.0)*(85000.0/182.0)*cos(-t[5]-pi/2);fe[nexf][2]=-(-t[1]/dr-90.0)*(85000.0/182.0)*sin(-t[5]-pi/2);} else{fe[nexf][1]=0;fe[nexf][2]=0;}}（2）四杆工艺阻力函数程序#include"math.h"extf(p,vp,ap,t,w,e,nexf,fe)double p[20][2],vp[20][2],ap[20][2],t[10],w[10],e[10],fe[20][2];int nexf; {double pi,dr;pi=4.0*atan(1.0);dr=pi/180.0;if(w[3]<0){fe[nexf][1]=(-t[1]/dr-2.0)*(85000.0/181.0)*cos(-t[3]-pi/2);fe[nexf][2]=-(-t[1]/dr-2.0)*(85000.0/181.0)*sin(-t[3]-pi/2);}else{fe[nexf][1]=0;fe[nexf][2]=0;}}6.2飞轮的转动惯量函数程序（1）程序：六杆铰链式颚式破碎机的转动惯量计算#include "graphics.h"#include "stdio.h"#include "subk.c"#include "subf.c"main(){static double p[20][2],vp[20][2],ap[20][2],del;static double t[10],w[10],e[10];static double tbdraw[370],tb1draw[370],y[370],en[370];static double fr[20][2],fe[20][2];static int ic;int i;double pi,dr,td,tb,emax,emin,x,jf; FILE*fp;sm[1]=0.0;sm[2]=500.0;sm[3]=200.0;sm[4]=200.0;sm[5]=900.0;sj[1]=0.0;sj[2]=25.5;sj[3]=9.0;sj[4]=9.0;sj[5]=50.0;y[0]=0.0;en[0]=0.0;td=0.0;pi=4.0*atan(1.0);dr=pi/180.0;w[1]=-17.8; e[1]=0.0; del=15.0;p[1][1]=0.0;p[1][2]=0.0;p[4][1]=0.94;p[4][2]=-1.0;p[6][1]=-1.0;p[6][2]=0.85;; if((fp=fopen("filejf6","w"))==NULL){printf("Can't open this file./n");exit(0);}fprintf(fp,"\nThe Moment of Inertia Analysis of a Six-bar Linkase:\n"); ic=(int)(360.0/del);for(i=0;i<=ic;i++){ t[1]=(-i)*del*dr;bark(1,2,0,1,0.1,0.0,0.0,t,w,e,p,vp,ap);rrrk(-1,2,4,3,2,3,1.25,1.0,t,w,e,p,vp,ap);rrrk(1,3,6,5,4,5,1.15,1.96,t,w,e,p,vp,ap);bark(2,0,9,2,0.0,0.625,0.0,t,w,e,p,vp,ap);bark(4,0,10,3,0.0,0.5,0.0,t,w,e,p,vp,ap);bark(3,0,8,4,0.0,0.575,0.0,t,w,e,p,vp,ap);bark(6,0,7,5,0.0,0.98,0.0,t,w,e,p,vp,ap);bark(6,0,11,5,0.0,0.6,0.0,t,w,e,p,vp,ap);rrrf(3,6,5,8,7,0,11,11,4,5,p,vp,ap,t,w,e,fr);rrrf(2,4,3,9,10,0,3,0,2,3,p,vp,ap,t,w,e,fr);barf(1,1,2,1,p,ap,e,fr,&tb);tbdraw[i]=tb;td=td+tbdraw[i];} td=td/ic;printf("%10.3f\n",td) ;for(i=0;i<=ic;i++){if(i>0){y[i]=del*dr*td-del*dr/2.0*(tbdraw[i]+tbdraw[i-1]); en[i]=en[i-1]+y[i];} printf("%10.3f%10.3f\n",y[i],en[i]);}for(i=0;i<25;i++)if(en[i]>en[i+1]){x=en[i];en[i]=en[i+1];en[i+1]=x;} emax=x;for(i=0;i<25;i++)if(en[i]<en[i+1]){x=en[i];en[i]=en[i+1];en[i+1]=x;}emin=x;jf=(emax-emin)/(w[1]*w[1]*0.05); printf("%10.3f\n",jf) ;fprintf(fp," jf=%10.3f kg*m*m\n",jf); }#include "math.h"extf(p,vp,ap,t,w,e,nexf,fe)double p[20][2],vp[20][2],ap[20][2],t[10],w[10],e[10],fe[20][2];{double pi,dr;pi=4.0*atan(1.0);dr=pi/180.0;if(w[5]<0){ fe[nexf][1]=(-t[1]/dr-90.0)*(85000.0/182.0)*cos(-t[5]-pi/2);fe[nexf][2]=-(-t[1]/dr-90.0)*(85000.0/182.0)*sin(-t[5]-pi/2);} else{fe[nexf][1]=0;fe[nexf][2]=0;}}（2）六杆铰链式颚式破碎机的转动惯量输出结果The Moment of Inertia Analysis of a Six-bar Linkase:Jf = 129.210 kg*m*m（3）程序：四杆铰链式颚式破碎机的转动惯量计算#include "graphics.h"#include "stdio.h"#include "subk.c"#include "subf.c"main(){static double p[20][2],vp[20][2],ap[20][2],del;static double t[10],w[10],e[10];static double tbdraw[370],tb1draw[370],y[370],en[370];static double fr[20][2],fe[20][2];static int ic;int i;double pi,dr,td,tb,emax,emin,x,jf;FILE*fp;sm[1]=0.0; sm[2]=200.0; sm[3]=900.0;sj[1]=0.0; sj[2]=9.0; sj[3]=50.0;y[0]=0.0;en[0]=0.0;td=0.0;pi=4.0*atan(1.0);dr=pi/180.0;w[1]=-17.8; e[1]=0.0; del=15.0;p[1][1]=0.0;p[1][2]=0.0;p[4][1]=-0.95;p[4][2]=2.0;if((fp=fopen("filejf6","w"))==NULL){printf("Can't open this file./n");exit(0);}fprintf(fp,"\nThe Moment of Inertia Analysis of a Four-bar Linkase:\n"); ic=(int)(360.0/del);for(i=0;i<=ic;i++){t[1]=(-i)*del*dr;bark(1,2,0,1,0.04,0.0,0.0,t,w,e,p,vp,ap);rrrk(1,2,4,3,2,3,1.11,1.96,t,w,e,p,vp,ap);bark(2,0,7,2,0.0,0.555,0.0,t,w,e,p,vp,ap);bark(4,0,5,3,0.0,0.98,0.0,t,w,e,p,vp,ap);bark(4,0,6,3,0.0,0.6,0.0,t,w,e,p,vp,ap);rrrf(2,4,3,7,5,0,6,6,2,3,p,vp,ap,t,w,e,fr);barf(1,1,2,1,p,ap,e,fr,&tb); tbdraw[i]=tb;td=td+tbdraw[i];}td=td/ic;for(i=0;i<=ic;i++){if(i>0){y[i]=del*dr*td-del*dr/2.0*(tbdraw[i]+tbdraw[i-1]); en[i]=en[i-1]+y[i];}printf("%10.3f%10.3f\n",y[i],en[i]);for(i=0;i<25;i++)if(en[i]>en[i+1]){x=en[i];en[i]=en[i+1];en[i+1]=x;}emax=x;for(i=0;i<25;i++)if(en[i]<en[i+1]){x=en[i];en[i]=en[i+1];en[i+1]=x;} emin=x;jf=(emax-emin)/(w[1]*w[1]*0.05); printf("%10.3f\n",jf);fprintf(fp," jf=%10.3f kg*m*m\n",jf);}#include "math.h"extf(p,vp,ap,t,w,e,nexf,fe)double p[20][2],vp[20][2],ap[20][2],t[10],w[10],e[10],fe[20][2];{double pi,dr;pi=4.0*atan(1.0);dr=pi/180.0;if(w[3]<0){ fe[nexf][1]=(-t[1]/dr-2.0)*(85000.0/181.0)*cos(-t[3]-pi/2);fe[nexf][2]=-(-t[1]/dr-2.0)*(85000.0/181.0)*sin(-t[3]-pi/2);;} else{fe[nexf][1]=0;fe[nexf][2]=0;}}（4）四杆铰链式颚式破碎机的转动惯量输出结果The Moment of Inertia Analysis of a Four-bar Linkase:jf= 39.882 kg*m*m七 .对两种机构的综合评价从运动角度分析：颚板的摆动范围：四杆机构Δt=0.040rad，六杆机构Δt=0.040rad；颚板最大角速度：四杆机构w[3]max=0.364306rad/s,六杆机构w[5]max=0.400445rad/s；颚板最大角加速度：四杆机构e[3]max=6.701549rad/s2, 六杆机构e[5]max=10.180768rad/s2。
