東北大學機械原理課程設計(牛頭刨床)

1、精選優(yōu)質文檔-傾情為你奉上精選優(yōu)質文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質文檔-傾情為你奉上專心-專注-專業(yè)機械原理課程設計說明書題目： 牛頭刨床機構方案分析 班 級 ：姓 名 ：學 號 ：指導教師 ：成 績 ：2012年 月 日目錄 TOC o 1-2 h z u 1題目說明（圖a）所示為牛頭刨床主機構的運動簡圖，已知，l1=0.1m,l0=0.4m,l3=0.75m,l4=0.15m,ly=0.738m,l3=0.375m,a=0.05m,b=0.15,c=0.4m,d=0.1m。只計構件3、5的質量，其余略去不計，m3=30kg,Js3=0.7kgm2,m5=95kg。

2、工藝阻力Q如圖所示，Q=9000N。主軸1的轉速為60r/min(順時針方向)，許用運轉不均勻系數=0.03。2內容要求與作法1.進行運動分析，繪制滑枕5的位移、速度和加速度線圖。2.進行動態(tài)靜力分析，繪制固定鉸鏈處的反力矢端圖金額平衡力矩Tb的線圖。3.以主軸為等效構件，計算飛輪的轉動慣量。上機前認真讀懂所用子程序，自編主程序，初始位置取滑枕5的左極限位置。3.滑枕初始位置及行程H的確定方法滑枕初始位置為左極限位置,此時構件1和構件3在O1O2左側垂直；滑枕行程H由運動分析的結果得知?；淼淖髽O限位置為-0.337，右極限為0.038，所以行程為0.375；4機構的結構分析213213 =

3、2 * GB3 84 = 3 * GB3 7 = 4 * GB3 主動件RPR桿組RRP桿組36955-1桿組運動分析子程序函數及形參與真值對照表1)調用bark函數，求3點的運動參數。形式參數n1n2n3kr1r2gamtwepvpap實 值1301r130.00.0twepvpap 2)調用rprk函數，求構件3的運動參數。形式參數mn1n2k1k2r1r2vr2ar2twepvpap實 值123320.0&r2&vr2&ar2twepvpap3)調用bark函數，求4點的運動參數。形式參數n1n2n3kr1r2gamtwepvpap實 值2403r240.00.0twepvpap4)調用

4、rrpk函數，求5點的位置和5構件的運動參數。形式參數mn1n2n3k1k2k3r1r2vr2ar2twepvpap實 值-1495456r45&r2&vr2&ar2twepvpap 5-2主程序、運行結果、運行圖像主程序/* Note:Your choice is C IDE */#include ”graphics.h” /*圖形庫*/#include subk.c /*運動分析子程序*/#include draw.c /*繪圖子程序*/main() static double p202,vp202,ap202,del; static double t10,w10,e20,pdraw370

5、,vpdraw370,apdraw370; static int ic; /*定義靜態(tài)變量*/ double r13,r24,r45,r12; double pi,dr; double r2,vr2,ar2; int i; /*定義局部變量*/ FILE *fp; /*定義文件指針變量*/char *m=p,vp,ap; /*定義字符串，改變字符串會改變繪圖結果的標注*/ r13=0.1;r24=0.75;r45=0.15;r12=0.4; pi=4.0*atan(1.0); /*求*/ dr=pi/180.0 /*求弧度*/; t6=0.0;w6=0.0;e6=0.0; w1=-2*pi;e

6、1=0.0; del=15.0; p21=0.0; p22=0.0; p11=0.0; p12=0.4; p91=-100; p92=0.738; /*變量賦值*/ t6=0.0*dr; /*變角度為弧度*/ printf(n The Kinematic Parameters of Point5n); printf(No THETA1 S5 V5 A5n); printf(n deg m m/s m/s/sn); /*在屏幕上寫表頭*/ if(fp=fopen(file1,w)=NULL) printf(Cant open this file.n); exit(0); /*建立并打開文件fil

7、e1*/ fprintf(fp,n The Kinematic Parameters of Point5n); fprintf(fp,No THETA1 S5 V5 A5n); fprintf(fp, deg m m/s m/s/sn); /*在文件file1中寫表頭*/ ic=(int)(360.0/del); for(i=0;i=ic;i+) /*建立循環(huán)，調用運動分析子程序*/ t1=(-(pi-asin(r13/r12)+(-i)*del*dr; bark(1,3,0,1,r13,0.0,0.0,t,w,e,p,vp,ap); rprk(1,2,3,3,2,0.0,&r2,&vr2,&

8、ar2,t,w,e,p,vp,ap); bark(2,4,0,3,r24,0.0,0.0,t,w,e,p,vp,ap); rrpk(-1,4,9,5,4,5,6,r45,&r2,&vr2,&ar2,t,w,e,p,vp,ap); /*r2,vr2,ar2在子程序中已定義為指針變量，所以其前面要加&符號*/ printf(n%2d%12.3f%12.3f%12.3f%12.3f,i+1,t1/dr,p51,vp51,ap51); /*把運算結果寫屏幕上*/ fprintf(fp,n%2d%12.3f%12.3f%12.3f%12.3f,i+1,t1/dr,p51,vp51,ap51); /*把運

9、算結果寫在文件中*/ pdrawi=p51; vpdrawi=vp51; apdrawi=ap51; /*將運算結果傳給pdrawi 、vpdrawi 、apdrawi以備繪圖用*/ if(i%16)=0)getch(); /*屏幕滿16行停頓*/ fclose(fp); /*關閉文件*/ getch(); draw1(del,pdraw,vpdraw,apdraw,ic,m); /*調用繪圖程序*/ /*主程序結束*/運行結果The Kinematic Parameters of Point5No THETA1 S5 V5 A5 deg m m/s m/s/s1 -165.522 -0.33

10、7 -0.000 7.251 2 -180.522 -0.331 0.273 5.878 3 -195.522 -0.315 0.492 4.688 4 -210.522 -0.291 0.665 3.611 5 -225.522 -0.260 0.794 2.593 6 -240.522 -0.225 0.882 1.623 7 -255.522 -0.187 0.930 0.715 8 -270.522 -0.148 0.942 -0.123 9 -285.522 -0.109 0.921 -0.90310 -300.522 -0.072 0.867 -1.66811 -315.522 -

11、0.037 0.781 -2.49512 -330.522 -0.007 0.657 -3.47013 -345.522 0.017 0.489 -4.67114 -360.522 0.033 0.265 -6.10815 -375.522 0.038 -0.022 -7.66016 -390.522 0.030 -0.371 -9.00017 -405.522 0.006 -0.762 -9.57818 -420.522 -0.033 -1.148 -8.64619 -435.522 -0.088 -1.450 -5.45920 -450.522 -0.152 -1.571 -0.05421

12、 -465.522 -0.216 -1.450 5.69922 -480.522 -0.270 -1.128 9.21023 -495.522 -0.308 -0.723 9.82024 -510.522 -0.330 -0.333 8.74625 -525.522 -0.337 0.000 7.251 運行圖像 圖5-1圖5-16.牛頭刨床的靜力分析6-1桿組運動分析子程序函數及形參實值對照表1)調用bark函數，求3點的運動參數。形式參數n1n2n3kr1r2gamtwepvpap實 值1301r130.00.0twepvpap2）調用rprk函數，求構件3的運動參數。形式參數mn1n2k

13、1k2r1r2vr2ar2twepvpap實 值123320.0&r2&vr2&ar2twepvpap3)調用bark函數，求4點的運動參數。形式參數n1n2n3kr1r2gamtwepvpap實 值24030.0r240.0twepvpap4)調用rrpk函數，求5點的位置和5構件的運動參數。形式參數mn1n2n3k1k2k3r1r2vr2ar2twepvpap實 值-1495456r45&r2&vr2&ar2twepvpap5)調用bark函數，求構件2的質心6點的運動參數。形式參數n1n2n3kr1r2gamtwepvpap實 值2603r260.00.0twepvpap6)調用bark

14、函數，求構件5的質心7點的運動參數。形式參數n1n2n3kr1r2gamtwepvpap實 值50750.0r57-161.57*drtwepvpap7)調用bark函數，求構件5的工藝阻力作用點8點的運動參數。形式參數n1n2n3kr1r2gamtwepvpap實 值50850.0r58165.96*drtwepvpap8) 調用rrpf函數，求4、5點的反作用力。形式參數n1n2n3ns1ns2nn1nn2nexfk1k2pvpaptwefr實 值4950708845pvpaptwefr9) 調用rprf函數求3點的反作用力。形式參數n1n2ns1ns2nn1nn2nexfk1k2pvpa

15、ptwefrfkpk實 值236040032pvpaptwefrfkpk10）調用barf函數，求1點運動副反力及平衡力矩。形式參數n1ns1nn1k1papefrtb實 值1031papefr&tb 6-2主程序、運行結果、運行圖像主程序/* Note:Your choice is C IDE */#include ”graphics.h” /*圖形庫*/#include subk.c /*運動分析子程序*/ #include subf.c /*動態(tài)靜力分析子程序*/#include draw.c /*繪圖子程序*/main() static double p202,vp202,ap202,

16、del; static double t10,w10,e10,tbdraw370,tb1draw370; static double sita1370,fr1draw370,sita2370,fr2draw370,sita3370,fr3draw370; static double fr202,fe202,fk202,pk202, tb,tb1,fr1,bt1,fr2,bt2,we1,we2,we3,we4,we5; static int ic; double r12,r13,r24,r45,r26,r58,r57; double pi,dr; double r2,vr2,ar2; int i

17、; FILE *fp; char *m=tb,tb1,fr1,fr2;/*定義字符串，前兩個字符串是為畫平衡力矩曲線提供的；后三個字符串是為畫反力矢端圖提供的；改變字符串會改變繪圖結果的標注*/ sm1=0.0;sm2=0.0;sm3=30.0;sm4=0.0;sm5=95.0; sj1=0.0;sj2=0.0;sj3=0.7;sj4=0.0; r13=0.1;r24=0.75;r45=0.15;r26=0.375; r58=sqrt(0.1*0.1+0.25*0.25); r57=sqrt(0.15*0.15+0.05*0.05);r12=0.4; pi=4.0*atan(1.0); dr=

18、pi/180.0; t6=0.0;w6=0.0;e6=0.0; w1=-2*pi;e1=0.0;del=15.0; t6=0.0*dr; p21=0.0; p22=0.0; p11=0.0; p12=0.4; p91=-100; p92=0.738; printf(n The Kineto-static Analysis of a six-bar Linkasen); printf(No HETAL fr1 sita1 fr2 sita2 tb tb1n); printf( deg N radian N radian N.m N.m ); if(fp=fopen(file1,w)=NULL)

19、printf(Cant open this file.n); exit(0); fprintf(fp,n The Kineto-static Analysis of a six-bar Linkasen); fprintf(fp,No HETAL fr1 sita1 fr2 sita2 tb tb1n); fprintf(fp, deg N radian N radian N.m N.m ); ic=(int)(360.0/del); for(i=0;i=(-0.31825)&p510)fenexf1=-9000.0;elsefenexf1=0.0;運行結果The Kineto-static

20、Analysis of a six-bar LinkaseNo HETAL fr1 sita1 fr2 sita2 tb tb1 deg N radian N radian N.m N.m 1 -165.522 1465.961 14.478 627.211 -171.015 0.000 0.000 2 -180.522 1127.912 14.005 446.916 179.556 -28.293 -28.293 3 -195.522 15767.371 12.723 6867.477 -148.658 -746.184 -746.184 4 -210.522 15065.899 10.81

21、8 6002.447 -154.166 -995.161 -995.161 5 -225.522 14516.560 8.455 5324.070 -162.398 -1174.074 -1174.074 6 -240.522 14074.069 5.769 4865.571 -172.533 -1288.627 -1288.627 7 -255.522 13725.674 2.881 4635.881 176.812 -1344.548 -1344.548 8 -270.522 13479.810 -0.104 4606.793 167.177 -1347.945 -1347.945 9 -

22、285.522 13355.544 -3.086 4730.512 159.624 -1304.217 -1304.21710 -300.522 13373.618 -5.964 4969.198 154.494 -1216.382 -1216.38211 -315.522 13549.209 -8.631 5310.503 151.601 -1083.634 -1083.63412 -330.522 13886.197 -10.968 5763.750 150.487 -900.827 -900.82713 -345.522 14373.388 -12.834 6345.168 150.61

23、5 -659.496 -659.49614 -360.522 1208.101 165.934 504.947 7.117 28.293 28.29315 -375.522 1602.230 165.525 711.386 -1.585 -2.929 -2.92916 -390.522 1999.234 166.143 957.246 -5.892 -57.335 -57.33517 -405.522 2257.789 167.966 1159.948 -6.456 -124.576 -124.57618 -420.522 2148.950 171.064 1170.305 -3.036 -1

24、68.380 -168.38019 -435.522 1409.272 175.286 812.249 7.743 -133.095 -133.09520 -450.522 7.272 -179.826 244.641 88.754 -0.727 -0.72721 -465.522 1484.336 5.037 866.590 169.198 138.980 138.98022 -480.522 2280.194 9.192 1242.566 -177.203 175.402 175.40223 -495.522 2270.202 12.203 1157.293 -170.237 121.22

25、4 121.22424 -510.522 1883.416 13.937 884.112 -167.926 50.459 50.45925 -525.522 1465.961 14.478 627.211 -171.015 -0.000 -0.000運行圖像圖6-1 圖6-1圖6-2圖6-27.飛輪的轉動慣量計算7-1飛輪轉動慣量的計算方法 1、一周期內驅動力矩功等于阻力功，所以：t*Td=1/2(Tr0+Tr1)+1/2(Tr1+Tr2)+.+1/2(Trn-1+Trn)因為Trn=Tr0，所以由上式可得：Td=/tTri=1/nTri2、間隔i-1、i內的盈虧功變化量E:E=Td-1/2*

26、（Tri+Tri-1）3、計算各點的盈虧功Ei：Ei=Ei-1+Ei4、找出最大和最小盈虧功：Emax,Emin5、計算飛輪的轉動慣量Jf: Jf=(Emax-Emin)/W， 7-2飛輪轉動慣量的程序及運行結果主程序/* Note:Your choice is C IDE */#include”stdio.h”#include”graphics,h” /*圖形庫*/#includesubk.c /*運動分析子程序*/ #includesubf.c /*動態(tài)靜力分析子程序*/ main() static double p202,vp202,ap202,del; static double t1

27、0,w10,e10,Tr370; static double fr202,fk202,pk202,fe202; static int ic; int j; double Td,sum1=0.0,D370,E370,Max,Min,Jf,q; double r12,r13,r24,r45,r26,r58,r57,tb; double pi,dr; double r2,vr2,ar2; int i; FILE *fp; E0=0.0; Max=E0; Min=E0; q=0.03; sm1=0.0;sm2=0.0;sm3=30.0;sm4=0.0;sm5=95.0; sj1=0.0;sj2=0.0

28、;sj3=0.7;sj4=0.0; r12=0.4;r13=0.1;r24=0.75;r45=0.15;r26=0.375; r58=sqrt(0.1*0.1+0.25*0.25); r57=sqrt(0.15*0.15+0.05*0.05); pi=4.0*atan(1.0); dr=pi/180.0; t6=0.0;w6=0.0;e6=0.0; w1=-2*pi;e1=0.0;del=15.0; t6=0.0*dr; p21=0.0; p22=0.0; p11=0.0; p12=0.4; p91=-100; p92=0.738; printf(n The Kineto-statia Ana

29、lysis of Ntbcn ); printf(nNo Theta1 Tr E n ); printf( deg. N.M 1/sn ); /*在屏幕上寫表頭*/ if(fp=fopen(file1,w)=NULL) printf(Cant open this file.n); exit(0); /*建立并打開文件file1*/ fprintf(fp,n The Kineto-statia Analysis of Ntbcn); fprintf(fp,nNo Theta1 Ter En); fprintf(fp, deg. N.M Jn); /*在文件中寫表頭*/ ic=(int)(360.

30、0/del); for(i=0;i=ic;i+) /*建立循環(huán)，調用運動分析子程序*/ t1=(-(pi-asin(r13/r12)+(-i)*del*dr; bark(1,3,0,1,r13,0.0,0.0,t,w,e,p,vp,ap); rprk(1,2,3,3,2,0.0,&r2,&vr2,&ar2,t,w,e,p,vp,ap); bark(2,4,0,3,r24,0.0,0.0,t,w,e,p,vp,ap); rrpk(-1,4,9,5,4,5,6,r45,&r2,&vr2,&ar2,t,w,e,p,vp,ap); bark(2,6,0,3,r26,0.0,0.0,t,w,e,p,vp

31、,ap); bark(5,0,7,5,0.0,r58,-161.57*dr,t,w,e,p,vp,ap); bark(5,0,8,5,0.0,r57,165.96*dr,t,w,e,p,vp,ap); rrpf(4,9,5,0,7,0,8,8,4,5,p,vp,ap,t,w,e,fr); rprf(2,3,6,0,4,0,0,3,2,p,vp,ap,t,w,e,fr,fk,pk); barf(1,0,3,1,p,ap,e,fr,&tb); /*r2,vr2,ar2,tb在子程序中已定義為指針變量，所以前面要加&符號*/ Tri=tb; Di=t1; for(j=1;j=ic;j+) sum1=sum1+Trj-1; Td=sum1/ic; for(j=1;j=ic;j+) Ej=Ej-1+del*dr*(Td-0.5*(Trj+Trj-1); for(j=1;j=ic;j+) if(Max=Ej) Min=Ej; J