机械原理课程设计—颚式破碎机设计说明书DOC

目 录

一 设计题目……………………………………………………………………1 二 已知条件及设计要求…………………………………………………1

2.1已知条件……………………………………………………………………1 2.2设计要求……………………………………………………………………2

三. 机构的结构分析…………………………………………………………2

3.1六杆铰链式破碎机………………………………………………………2 3.2四杆铰链式破碎机………………………………………………………2

四. 机构的运动分析…………………………………………………………2

4.1六杆铰链式颚式破碎机的运动分析…………………………………2 4.2四杆铰链式颚式破碎机的运动分析…………………………………6

五.机构的动态静力分析……………………………………………………7

5.1六杆铰链式颚式破碎机的静力分析…………………………………7 5.2四杆铰链式颚式破碎机的静力分析…………………………………12

六. 工艺阻力函数及飞轮的转动惯量函数 ……………………17

6.1工艺阻力函数程序 ……………………………………………………17 6.2飞轮的转动惯量函数程序……………………………………………17

七 .对两种机构的综合评价……………………………………………21 八 . 主要的收获和建议…………………………………………………22 九 . 参考文献………………………………………………………………22

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

一.设计题目：铰链式颚式破碎机方案分析 二.已知条件及设计要求 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.5kg?m2, m3 = 200kg, Js3 = 9kg?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。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=9kg?m2，m3 = 900kg, Js3=50kg?m2。曲柄1的质心在O1 点处，2、3构件的质心在各构件的中心。

1

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

2.2 设计要求

试比较两个方案进行综合评价。主要比较以下几方面：

1. 进行运动分析，画出颚板的角位移、角速度、角加速度随曲柄转角的变化曲线。

2. 进行动态静力分析，比较颚板摆动中心运动副反力的大小及方向变化规律，曲柄上的平衡力矩大小及方向变化规律。

3. 飞轮转动惯量的大小。 三.机构的结构分析 3.1六杆铰链式破碎机

六杆铰链式粉碎机拆分为机架和主动件①，②③构件组成的RRR杆组，④⑤构件组成的RRR杆组。

+ +

3.2四杆铰链式破碎机

四杆铰链式破碎机拆分为机架和主动件①，②③构件组成的RRR杆组。

+

四.机构的运动分析

4.1六杆铰链式颚式破碎机的运动分析。 4.1.1 杆件的运动参数。

1) 调用bark函数求主动件①的运动参数。 形式参数 n1 n2 n3 实 值 1 2 0 k 1 r1 r2 game 0.0

t t w w e e p p vp ap vp ap r12 0.0 2）调用rrrk函数求②、③构件组成的RRR杆组进行运动分析。 形式参数 m n1 n2 n3 k1 k2 r1 R2实 值 -1 2 4 3 2 3 3 3）调用rrrk函数对④、⑤构件组成的RRR杆组进行运动分析。 形式参数 m n1 n2 n3 k1 k2 r1 2

r2 t w e p vp ap R34 t w e p vp ap r2 t w e p vp ap 东北大学机械原理课程设计 铰链式颚式破碎机方案分析

实 值 1 3 6 5 4 5 r35 R56 t w e p vp ap 4.1.2 写主程序并运行。

按一定的步长，改变主动件的位置角度，使其在0-360°变化，便可求出机构各点在整个运动循环内的运动参数并打印输出。

（1） 主程序。 #include\#include\main() {

static double p[20][2],vp[20][2],ap[20][2],del; static double t[10],w[10],e[10],pdraw[370], vpdraw[370],apdraw[370]; static int ic;

double r12,r23,r34,r35,r56,r611; double pi,dr;int i; FILE *fp;

r12=0.1;r23=1.250;r34=1.0;r35=1.15;r56=1.96;r611=0.6; pi=4.0*atan(1.0);dr=pi/180.0;

w[1]=-170*2*pi/60;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;p[6][1]=-1; p[6][2]=0.85;

printf(\printf(\

printf(\if((fp=fopen(\{

printf(\exit(0); }

fprintf(fp,\fprintf(fp,\

fprintf(fp,\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);

3

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

bark(2,0,7,2,0.0,r23/2,0.0,t,w,e,p,vp,ap); bark(4,0,8,3,0.0,r34/2,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(6,0,10,5,0.0,r56/2,0.0,t,w,e,p,vp,ap); bark(6,0,11,5,0.0,r611,0.0,t,w,e,p,vp,ap);

printf(\);

fprintf(fp,\e[5]);

pdraw[i]=t[5]; vpdraw[i]=w[5]; apdraw[i]=e[5];

if((i)==0){getch();} }

fclose(fp); getch();

draw1(del,pdraw,vpdraw,apdraw,ic); }

（2）运行结果。 ① 件5的运动参数：

The Kinematic Parameters of Point 5

No THETA1 S5 V5 A5

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.904 10 -135.000 -1.641 -0.345 -3.981 11 -150.000 -1.646 -0.382 -1.008 12 -165.000 -1.652 -0.377 1.519 13 -180.000 -1.657 -0.341 3.297 14 -195.000 -1.662 -0.284 4.237 15 -210.000 -1.666 -0.220 4.436 16 -225.000 -1.668 -0.156 4.121 17 -240.000 -1.670 -0.10 3.584

4

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

18 -255.000 -1.671 -0.051 3.105 19 -270.000 -1.672 -0.007 2.898 20 -285.000 -1.672 0.036 3.063 21 -300.000 -1.671 0.085 3.571 22 -315.000 -1.669 0.142 4.247 23 -330.000 -1.667 0.209 4.791 24 -345.000 -1.663 0.281 4.817 25 -360.000 -1.658 0.346 3.956 ② 动图形：

4.2杆铰链式颚式破碎机的运动分析。 4.2.1运动参数。

1）调用bark函数求主动件①的运动参数。 形式参数 n1 n2 n3 k 实 值 1 2 0 1 r1 r2 game t t w w e e p p vp ap vp ap r12 0.0 0.0 2）调用rrrk函数求②、③构件组成的RRR杆组进行运动分析。 形式参数 m 实 值 1 n1 n2 n3 k1 k2 r1 2 4 3 2 3 r2 t w t w e e p vp ap p vp ap r23 r34 4.2.2 写主程序并运行。

按一定的步长，改变主动件的位置角度，使其在0-360°变化，便可求出机构各点在整个运动循环内的运动参数并打印输出。 (1)主程序。

#include\#include\main()

5

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

{

static double p[20][2],vp[20][2],ap[20][2],del;

static double t[10],w[10],e[10],pdraw[370],vpdraw[370],apdraw[370]; static int ic;

double r12,r23,r34,r47; double pi,dr;int i; FILE *fp;

r12=0.04;r23=1.11;r34=1.96;r47=0.6; pi=4.0*atan(1.0); dr=pi/180.0;

w[1]=-170*2*pi/60;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; printf(\printf(\

printf(\if((fp=fopen(\{

printf(\exit(0); }

fprintf(fp,\fprintf(fp,\

fprintf(fp,\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); bark(2,0,5,2,0.0,r23/2,0.0,t,w,e,p,vp,ap); bark(4,0,6,3,0.0,r34/2,0.0,t,w,e,p,vp,ap); bark(4,0,7,3,0.0,0.6,0.0,t,w,e,p,vp,ap);

printf(\fprintf(fp,\);

pdraw[i]=t[3];vpdraw[i]=w[3];apdraw[i]=e[3]; if((i)==0){getch();} }

fclose(fp); getch();

6

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

draw1(del,pdraw,vpdraw,apdraw,ic); }

(2)运行结果。

①杆件3的运动参数：

The Kinematic Parameters of Point 5

No THETA1 S5 V5 A5

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.848 10 -135.000 -1.666 -0.274 4.339 11 -150.000 -1.669 -0.201 5.544 12 -165.000 -1.671 -0.113 6.358 13 -180.000 -1.672 -0.016 6.703 14 -195.000 -1.672 0.082 6.545 15 -210.000 -1.670 0.174 5.894 16 -225.000 -1.667 0.253 4.807 17 -240.000 -1.663 0.313 3.384 18 -255.000 -1.658 0.351 1.746 19 -270.000 -1.653 0.364 0.030 20 -285.000 -1.647 0.352 -1.639 21 -300.000 -1.642 0.317 -3.149 22 -315.000 -1.638 0.261 -4.415 23 -330.000 -1.635 0.189 -5.375 24 -345.000 -1.632 0.105 -5.988 25 -360.000 -1.632 0.014 -6.232 ②运动图形：

7

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

五.机构的动态静力分析

5.1六杆铰链式颚式破碎机的动态静力分析。

5.1.1 质点7，8，9，10及矿石破碎产生阻力的作用点11的运动参数； ①调用bark函数对质点7进行运动分析： 形式参n1 n2 n3 数 实 值 2 0 7 2 0.0 R23/2 0.0 t w e p vp ap k r1 r2 game t w e p vp ap ②调用bark函数对质点8进行运动分析： 形式参n1 n2 n3 数 实 值 4 0 8 3 0.0 R34/2 0.0 t w e p vp ap k r1 r2 game t w e p vp ap ③调用bark函数对质点9进行运动分析： 形式参n1 n2 n3 数 实 值 3 0 9 4 0.0 R35/2 0.0 t w e p vp ap k r1 r2 game t w e p vp ap ④调用bark函数对质点10进行运动分析： 形式参n1 n2 n3 k r1 r2 game t w e p vp ap 8

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

数 实 值 6 0 10 5 0.0 R56/2 0.0 t w e p vp ap ⑤调用bark函数对质点11进行运动分析： 形式参数 n1 n2 n3 实 值 6 0 11 k 5 r1 r2 game 0.0 t t w w e e p p vp ap vp ap 0.0 r611 （2）调用rrrf函数对④、⑤构件构成的RRR杆组进行动态静力分析： 形式参数 实 值 形式参数 实 值 n1 3 k2 5 n2 6 p p n3 5 vp vp ns1 9 ap ap ns2 10 t t nn1 0 w w nn2 11 e e nexf 11 fr fr k1 4 （3）调用rrrf函数对②、③构件构成的RRR杆组进行动态静力分析： 形式参数 实 值 形式参数 实 值 n1 2 k2 3 n2 4 p p n3 3 vp vp ns1 7 ap ap ns2 8 t t nn1 0 w w

nn2 3 e e nexf 0 fr fr k1 2 （4）调用barf函数对主动件1进行动态静力分析： 形式参数 n1 ns1 nn1 k1 实 值 1 1 2 1 p p ap ap e e fr fr tb

&tb

5.1.2程序并运行。

按一定的步长，改变主动件的位置角度，使其在0-360°变化，便可求出机构各运动副反力及作用在主动件上的平衡力矩。 （1）主程序。

#include \#include\#include\#include\main()

9

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

{

static double p[20][2],vp[20][2],ap[20][2],del;

static double t[10],w[10],e[10],fr[20][2],fe[20][2];

static double sita1[370],fr1draw[370],sita2[370],fr2draw[370], sita3[370],fr3draw[370],tbdraw[370],tb1draw[370];

static double tb,tb1,fr1,bt1,fr4,bt4,fr6,bt6,we1,we2,we3,we4,we5; static int ic;

double r12,r23,r34,r35,r56,r611;double pi,dr;int i; 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; r12=0.1;r23=1.250;r34=1.0;r35=1.15;r56=1.96;r611=0.6; pi=4.0*atan(1.0);dr=pi/180.0;

w[1]=-170*2*pi/60;e[1]=0.0;del=10.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(\printf(\

printf(\

printf(\printf(\

printf(\if((fp=fopen(\{

printf(\exit(0); }

fprintf(fp,\Kineto-static Analysis of a Six-bar Linkase\\n\fprintf(fp,\

fprintf(fp,\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); bark(2,0,7,2,0.0,r23/2,0.0,t,w,e,p,vp,ap); bark(4,0,8,3,0.0,r34/2,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);

10

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

bark(6,0,10,5,0.0,r56/2,0.0,t,w,e,p,vp,ap); bark(6,0,11,5,0.0,r611,0.0,t,w,e,p,vp,ap); rrrf(3,6,5,9,10,0,11,11,4,5,p,vp,ap,t,w,e,fr); rrrf(2,4,3,7,8,3,0,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[7][1]*vp[7][1]+(ap[7][2]+9.81)*vp[7][2])*sm[2]-e[2]*w[2]*sj[2];

we3=-(ap[8][1]*vp[8][1]+(ap[8][2]+9.81)*vp[8][2])*sm[3]-e[3]*w[3]*sj[3];

we4=-(ap[9][1]*vp[9][1]+(ap[9][2]+9.81)*vp[9][2])*sm[4]-e[4]*w[4]*sj[4];

extf(p,vp,ap,t,w,e,11,fe);

we5=-(ap[10][1]*vp[10][1]+(ap[10][2]+9.81)*vp[10][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(\/dr,tb,tb1);

fprintf(fp,\,bt6/dr,tb,tb1);

tbdraw[i]=tb;tb1draw[i]=tb1;fr1draw[i]=fr1;sita1[i]=bt1; fr2draw[i]=fr4;sita2[i]=bt4; fr3draw[i]=fr6;sita3[i]=bt6; if((i)==0){getch();} }

fclose(fp); getch();

draw2(del,tbdraw,tb1draw,ic);

draw3(del,sita1,fr1draw,sita2,fr2draw,sita3,fr3draw,ic); }

#include \

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]; {

11

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

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 Kineto-static Analysis of a Six-bar Linkase NO THETA1 FR6 BT6 TB TB1 (deg.) (N) (deg.) (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 12694.785 130.968 -883.503 -883.503 9 -120.000 15067.304 144.368 -1406.751 -1406.751 10 -135.000 17747.445 153.528 -1623.932 -1623.932 11 -150.000 20694.648 160.499 -1555.124 -1555.124 12 -165.000 24022.489 166.115 -1286.981 -1286.981 13 -180.000 27824.537 170.602 -923.989 -923.989 14 -195.000 32106.137 174.034 -557.444 -557.444 15 -210.000 36785.996 176.536 -253.496 -253.496 16 -225.000 41731.070 178.304 -51.307 -51.307 17 -240.000 46795.855 179.558 36.912 36.912 18 -255.000 51853.137 -179.507 26.634 26.634 19 -270.000 56814.061 -178.745 -39.410 -39.410 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 （3）平衡力矩曲线：

12

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

(4)反力的矢端图曲线：

5.2杆铰链式颚式破碎机的动态静力分析。

5.2.1 质点5，6及矿石破碎产生阻力的作用点7的运动参数； ①调用bark函数对质点5进行运动分析： 形式参数 n1 n2 n3 实 值 2 0 5 k 2 r1 r2 gam t t w w e e p p vp ap vp ap 0.0 r25 0.0 ②调用bark函数对质点6进行运动分析： 形式参数 n1 n2 n3 实 值 4 0 6 k 3 r1 r2 gam t t w w e e p p vp ap vp ap 0.0 r46 0.0 ③调用bark函数对质点7进行运动分析：

13

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

形式参数 n1 n2 n3 实 值 4 0 7 k 3 r1 r2 gam t t w w e e p p vp ap vp ap 0.0 r47 0.0 （2）调用rrrf函数对②，③构件构成的RRR杆组进行动态静力分析： 形式参数 实 值 形式参数 实 值 n1 2 k2 3 n2 4 p p n3 3 vp vp ns1 5 ap ap ns2 6 t t nn1 0 w w

nn2 7 e e nexf 7 fr fr k1 2 （3）调用barf函数对主动件1进行动态静力分析： 形式参数 n1 ns1 nn1 k1 实 值 行。

按一定的步长，改变主动件的位置角度，使其在0-360°变化，便可求出机构各运动副反力及作用在主动件上的平衡力矩。 （1）主程序。

#include \#include \#include \#include \main() {

static double p[20][2],vp[20][2],ap[20][2],del; static double t[10],w[10],e[10];

static double sita1[370],fr1draw[370],sita2[370],fr2draw[370], sita3[370],fr3draw[370],tbdraw[370],tb1draw[370];

static double fr[20][2],fe[20][2];static int ic;double r12,r23,r34,r47; int i;double pi,dr;double fr1,fr4,bt1,bt4,we1,we2,we3,tb,tb1; 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; r12=0.04; r23=1.11; r34=1.96,r47=0.6; pi=4.0*atan(1.0);

w[1]=-170*2*pi/60; e[1]=0.0; del=15.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;

14

p p ap ap e e fr fr tb

1 1 2 1 &tb

5.2．2主程序并运

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

printf(\printf(\

printf(\if((fp=fopen(\{

printf(\exit(0); }

fprintf(fp,\Kineto-static Analysis of a four-bar Linkase\\n\fprintf(fp,\

fprintf(fp,\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); bark(2,0,5,2,0.0,r23/2,0.0,t,w,e,p,vp,ap); bark(4,0,6,3,0.0,r34/2,0.0,t,w,e,p,vp,ap); bark(4,0,7,3,0.0,0.6,0.0,t,w,e,p,vp,ap); rrrf(2,4,3,5,6,0,7,7,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[5][1]*vp[5][1]+(ap[5][2]+9.81)*vp[5][2])*sm[2]-e[2]*w[2]*sj[2];

extf(p,vp,ap,t,w,e,7,fe);

we3=-(ap[6][1]*vp[6][1]+(ap[6][2]+9.81)*vp[6][2])*sm[3]-e[3]*w[3]*sj[3]+fe[7][1]*vp[7][1]+fe[7][2]*vp[7][2]; tb1=-(we1+we2+we3)/w[1];

printf(\\\n\

fprintf(fp,\\\n\tbdraw[i]=tb;tb1draw[i]=tb1; fr1draw[i]=fr1;sita1[i]=bt1;

15

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

fr2draw[i]=fr4;sita2[i]=bt4; fr3draw[i]=fr4;sita3[i]=bt4; if((i)==0){getch();} }

fclose(fp); getch();

draw2(del,tbdraw,tb1draw,ic);

draw3(del,sita1,fr1draw,sita2,fr2draw,sita3,fr3draw,ic); }

#include \

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-3/dr)*(85000.0/181.0)*cos(-t[3]-pi/2); fe[nexf][2]=(-t[1]/dr-3/dr)*(85000.0/181.0)*sin(-t[3]-pi/2); }

else{fe[nexf][1]=0;fe[nexf][2]=0;} }

(2)运行结果：

The Kineto-static Analysis of a four-bar Linkase NO THETA1 FR4 BT4 TB TB1 (deg.) (N) (deg.) (N.m) (N.m)

1 0.000 10261.396 102.685 47.546 47.546 2 -15.000 51448.370 17.786 186.094 186.094 3 -30.000 46996.739 19.370 316.446 316.446 4 -45.000 42725.338 21.193 402.027 402.027 5 -60.000 38603.813 23.247 446.798 446.798 6 -75.000 34590.628 25.540 455.194 455.194 7 -90.000 30637.112 28.117 431.327 431.327 8 -105.000 26695.709 31.108 379.118 379.118 9 -120.000 22733.412 34.813 303.360 303.360 10 -135.000 18753.920 39.876 211.317 211.317 11 -150.000 14843.055 47.694 114.117 114.117 12 -165.000 11289.242 61.310 27.042 27.042 13 -180.000 8894.730 85.799 -32.028 -32.028

16

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

14 -195.000 9702.566 70.652 -75.371 -75.371 15 -210.000 9538.373 72.140 -106.622 -106.622 16 -225.000 9408.540 74.864 -112.809 -112.809 17 -240.000 9336.371 78.581 -90.324 -90.324 18 -255.000 9335.572 82.915 -43.909 -43.909 19 -270.000 9405.058 87.439 15.069 15.069 20 -285.000 9529.828 91.754 72.279 72.279 21 -300.000 9687.086 95.554 114.547 114.547 22 -315.000 9853.708 98.638 133.092 133.092 23 -330.000 10011.499 100.889 125.150 125.150 24 -345.000 10149.104 102.245 93.963 93.963 25 -360.000 10261.396 102.685 47.546 47.546 （3）平衡力矩曲线：

（4）反力的矢端图曲线：

六.工艺阻力函数及飞轮的转动惯量函数 6.1工艺阻力函数程序：

（1）六杆铰链式颚式破碎机的工艺阻力函数程序：

17

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

#include \

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）四杆铰链式颚式破碎机的工艺阻力函数程序 #include \

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-3/dr)*(85000.0/181.0)*cos(-t[3]-pi/2); fe[nexf][2]=(-t[1]/dr-3/dr)*(85000.0/181.0)*sin(-t[3]-pi/2); }

else{fe[nexf][1]=0;fe[nexf][2]=0;} }

6.2 飞轮转动惯量函数程序

6.2.1六杆铰链式颚式破碎机的转动惯量函数程序 #include \#include \#include \#include \

static double sm[10],sj[10]; main() {

static double p[20][2],vp[20][2],ap[20][2],del; static double t[10],w[10],e[10],Tr[370],D[370]; static double fr[20][2],fe[20][2];static int ic;

18

东北大学机械原理课程设计 铰链式颚式破碎机方案分析

double r12,r23,r34,r35,r56,r611,tb,tb1;

double Td,sum1=0.0,E[370],Max=0.0,Min=0.0,Jf=0.0,q;int i,j; double pi,dr,fr1,bt1,fr2,bt2; 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;

r12=0.1; r23=1.25; r34=1.0; r35=1.15;r56=1.96,r611=0.6,del=15.0; pi=4.0*atan(1.0);w[1]=-170*2*pi/60; e[1]=0.0;dr=pi/180.0; q=0.3;

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(\printf(\printf(\if((fp=fopen(\{

printf(\exit(0); }

printf(fp,\printf(fp,\printf(fp,\(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,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(2,0,7,2,0.0,r23/2,0.0,t,w,e,p,vp,ap); bark(4,0,8,3,0.0,r34/2,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(6,0,10,5,0.0,r56/2,0.0,t,w,e,p,vp,ap); bark(6,0,11,5,0.0,r611,0.0,t,w,e,p,vp,ap); rrrf(3,6,5,9,10,0,11,11,4,5,p,vp,ap,t,w,e,fr); rrrf(2,4,3,7,8,3,0,0,2,3,p,vp,ap,t,w,e,fr); barf(1,1,2,1,p,ap,e,fr,&tb); Tr[i]=tb; D[i]=t[1]; }

19

本文来源：https://www.bwwdw.com/article/x7hr.html