機(jī)械原理大作業(yè)2-凸輪機(jī)構(gòu)設(shè)計(jì)

1、Harbin Institute of Technology機(jī)械原理大作業(yè)二課程名稱： 機(jī)械原理 設(shè)計(jì)題目： 凸輪機(jī)構(gòu)設(shè)計(jì) 院 系： 機(jī)電工程學(xué)院 姓 名： 學(xué) 號(hào)： 班 級(jí)： 指導(dǎo)教師： 1.設(shè)計(jì)題目設(shè)計(jì)直動(dòng)從動(dòng)件盤形凸輪機(jī)構(gòu)，其原始參數(shù)見(jiàn)表1表一： 凸輪機(jī)構(gòu)原始參數(shù)升程（mm)升程運(yùn)動(dòng)角（º）升程運(yùn)動(dòng)規(guī)律升程許用壓力角（º）回程運(yùn)動(dòng)角（º）回程運(yùn)動(dòng)規(guī)律回程許用壓力角（º）遠(yuǎn)休止角（º）近休止角（º）3070等加等減速30170正弦加速度601001202.凸輪推桿運(yùn)動(dòng)規(guī)律（1） 推程運(yùn)動(dòng)規(guī)律（等加速等減速運(yùn)動(dòng)） 推程 s=2h02

2、=607362=4h102=1201(718)2a=4h1202=12012(718)2 推程 s=h-2h020-2=30-60(718)20-2=4h1020-=12017182718-a=-4h1202=-12012(718)2（2） 回程運(yùn)動(dòng)規(guī)律（正弦加速度） 回程 s=h1-0-s0'+12sin20'-0-s=-h10'1-cos20'-0-sa=-2h120'2sin20'-0-s 0推程運(yùn)動(dòng)角；s遠(yuǎn)休止角；0'回程運(yùn)動(dòng)角開(kāi)始輸入初始參數(shù)運(yùn)行各部分程序從動(dòng)件位移、速度、加速度曲線輸出1 輸出2ds/d-s曲線，確定基圓半徑和

3、偏距理論輪廓線上的壓力角和曲率半徑圖輸出3輸出4繪制理論輪廓線和實(shí)際輪廓線結(jié)束3.運(yùn)動(dòng)線圖及凸輪線圖采用Matlab編程，其所有源程序見(jiàn)附頁(yè)：令1=2*160 rad/s可得運(yùn)動(dòng)規(guī)律圖如下：1. 凸輪的基圓半徑和偏距 以ds/df-s圖為基礎(chǔ)，可分別作出三條限制線（推程許用壓力角的切界限D(zhuǎn)tdt,回程許用壓力角的限制線Dt'dt'，起始點(diǎn)壓力角許用線B0d'')，以這三條線可確定最小基圓半徑及所對(duì)應(yīng)的偏距e,在其下方選擇一合適點(diǎn)，即可滿足壓力角的限制條件。得圖如下：得最小基圓對(duì)應(yīng)的坐標(biāo)位置大約為（20，-35）經(jīng)計(jì)算取偏距e=20mm，r0=40.3mm.2.繪

4、制理論輪廓線上的壓力角曲線和曲率半徑曲線針對(duì)凸輪轉(zhuǎn)向及推桿偏置，令N1=1凸輪逆時(shí)針轉(zhuǎn)；N2=1偏距為正。壓力角數(shù)學(xué)模型：=tan-1|dsd-N1N2e|s+r02-e2曲率半徑數(shù)學(xué)模型：其中：3.凸輪理論廓線和實(shí)際廓線理論廓線數(shù)學(xué)模型： 凸輪實(shí)際廓線坐標(biāo)方程式： 其中rt為確定的滾子半徑。 根據(jù)上面公式，利用matlab編程求解，其代碼如下：令rt=10mm得圖如下： *附頁(yè)相關(guān)程序代碼：1>位移曲線phi1=linspace(0,70/2/180*pi); %推程phi2=linspace(70/2/180*pi,70/180*pi);phi0=70/180*pi; h=30;s1

5、=2*h*(phi1/phi0).2;s2=h-2*h*(phi0-phi2).2/(phi0).2;plot(phi1,s1)hold onplot(phi2,s2)hold onphi3=linspace(deg2rad(70),deg2rad(170),1000); %遠(yuǎn)休程s3=h;plot(phi3,s3,'-b')hold onphi4=linspace(deg2rad(170),deg2rad(240); %回程s4=h*(1-(phi4-phi0-deg2rad(100)/deg2rad(70)+1/2/pi*sin(2*pi*(phi4-phi0-deg2ra

6、d(100)/deg2rad(70);plot(phi4,s4)hold onphi5=linspace(deg2rad(240),deg2rad(360),1000); %近休程s5=0;plot(phi5,s5,'-b')axis(0,4/3*deg2rad(360),0,0+3/2*h);grid onylabel('位移/mm');xlabel('凸輪轉(zhuǎn)角/rad');2>速度曲線phi0=deg2rad(70); %推程phi1=linspace(0,deg2rad(70/2);h=30; n=1; w=2*pi*n/60; v1

7、=4*h*w/(deg2rad(70)2*phi1;plot(phi1,v1)hold onphi2=linspace(deg2rad(70/2),deg2rad(70);v2=4*h*w*(phi0-phi2)/(deg2rad(70)2;plot(phi2,v2)hold onphi3=linspace(deg2rad(70),deg2rad(170),1000); %遠(yuǎn)休程v3=0;plot(phi3,v3,'-b')hold onphi4=linspace(deg2rad(170),deg2rad(240); %回程v4=-h*w/deg2rad(70)*(1-cos(

8、2*pi/deg2rad(70)*(phi4-deg2rad(70)-deg2rad(100);plot(phi4,v4)hold onphi5=linspace(deg2rad(240),deg2rad(360),1000); %近休程v5=0;plot(phi5,v5,'-b')grid onylabel('速度/mm/s');xlabel('凸輪轉(zhuǎn)角/rad');3>加速度曲線phi0=deg2rad(70); phi1=linspace(0,deg2rad(70/2),103);h=30; n=1; w=2*pi*n/60; a1=

9、4*h*w2/(phi0)2;plot(phi1,a1,'-b')hold onphi2=linspace(deg2rad(70/2),deg2rad(70),103); a2=-4*h*w2/(phi0)2;plot(phi2,a2,'-b')hold onphi3=linspace(deg2rad(70),deg2rad(170),1000); a3=0;plot(phi3,a3,'-b')hold onphi4=linspace(deg2rad(170),deg2rad(240); a4=-2*pi*h*w2/(deg2rad(70)2*s

10、in(2*pi/deg2rad(70)*(phi4-deg2rad(70)-deg2rad(100);plot(phi4,a4)hold onphi5=linspace(deg2rad(240),deg2rad(360),1000); a5=0;plot(phi5,a5,'-b')hold onphi6=deg2rad(35),deg2rad(35);a6=4*h*w2/(phi0)2,-4*h*w2/(phi0)2;plot(phi6,a6)grid onylabel('加速度/mm/s2');xlabel('凸輪轉(zhuǎn)角/rad');4>d

11、sd-s曲線phi1=linspace(0,70/2/180*pi,1000); phi2=linspace(70/2/180*pi,70/180*pi,1000);phi0=70/180*pi; h=30;s01=2*h*(phi1/phi0).2;s02=h-2*h*(phi0-phi2).2/(phi0).2;phi3=linspace(deg2rad(70),deg2rad(170),1000); s03=h;phi4=linspace(deg2rad(170),deg2rad(240); s04=h*(1-(phi4-phi0-deg2rad(100)/deg2rad(70)+1/2

12、/pi*sin(2*pi*(phi4-phi0-deg2rad(100)/deg2rad(70);phi5=linspace(deg2rad(240),deg2rad(360),1000) s05=0; %ds/fs1=4*h*phi1/(phi0)2;s2=4*h*(phi0-phi2)/(phi0).2;s3=0;s4=h*(-1/deg2rad(70)+1/deg2rad(70)*cos(2*pi*(phi4-phi0-deg2rad(100)/deg2rad(70);s5=0;x=s1,s2,s3,s4,s5;y=s01,s02,s03,s04,s05;plot(x,y)hold on

13、 %確定凸輪半徑和偏距k1=tan(pi/2-30/180*pi);k2=-tan(pi/2-60/180*pi);y1=-k1*x+y;y2=-k2*x+y;y1s=min(y1);y2s=min(y2);x=linspace(-50,50,1000);d1=k1*x+y1s;d2=k2*x+y2s;x0=linspace(0,60,200);d0=-k1*x0;plot(x,d1)hold onplot(x,d2)hold onplot(x0,d0)grid onaxis(-55,55,-80,40)axis equal%由圖像得r0=40.3,e=20xlabel('ds/d&#

14、166;×')ylabel('s/mm')5>壓力角曲線h=30;r0=40.3;e=20; N1=1;% 凸輪逆時(shí)針轉(zhuǎn) N2=1;% 偏距為正 n=1;omiga=2*pi*n/60;s0=sqrt(r02-e2);phi1=linspace(0,70/2/180*pi); phi2=linspace(70/2/180*pi,70/180*pi);phi0=70/180*pi;s1=2*h*(phi1/phi0).2;s20=h-2*h*(phi0-phi2).2/(phi0).2;v1=4*h*omiga/(deg2rad(70)2*phi1;v2=

15、4*h*omiga*(phi0-phi2)/(deg2rad(70)2;phi3=linspace(deg2rad(70),deg2rad(170),1000); s3=h;phi4=linspace(deg2rad(170),deg2rad(240); s40=h*(1-(phi4-phi0-deg2rad(100)/deg2rad(70)+1/2/pi*sin(2*pi*(phi4-phi0-deg2rad(100)/deg2rad(70);v4=-h*omiga/deg2rad(70)*(1-cos(2*pi/deg2rad(70)*(phi4-deg2rad(70)-deg2rad(1

16、00);phi5=linspace(deg2rad(240),deg2rad(360),1000); s5=0;%計(jì)算壓力角DsDphi1=v1./omiga;rs1=s0+s1;ang1=abs(atan(DsDphi1-N1*N2*e)./rs1)*180/pi;% ang1：推程壓力角DsDphi2=v2./omiga;rs3=s0+s20;ang2=abs(atan(DsDphi2-N1*N2*e)./rs3)*180/pi; %ang3：推程壓力角DsDphi3=v4./omiga;rs5=s0+s40ang4=abs(atan(DsDphi3-N1*N2*e)./rs5)*180/

17、pi; %ang5：回程壓力角rs2=(s0+h)*ones(1,size(s3,2);rs4=s0*ones(1,size(s5,2);ang3=abs(atan(-N1*N2*e)./rs2)*180/pi+0*phi3; %ang2：遠(yuǎn)休止壓力角ang5=abs(atan(-N1*N2*e)./rs4)*180/pi+0*phi5; %ang4：近休止壓力角ang=ang1,ang2,ang3,ang4,ang5;phi=phi1,phi2,phi3,phi4,phi5;plot(phi,ang)grid onylabel('壓力角/°');xlabel(

18、9;凸輪轉(zhuǎn)角/rad');title('壓力角曲線')6>曲率半徑曲線h=30; r0=40.3; e=20;N1=1;% 凸輪逆時(shí)針轉(zhuǎn) N2=1;% 偏距為正 rt=10;% 滾子半徑 n=1; M=-1;omiga=2*pi*n/60;s0=sqrt(r02-e2);phi1=linspace(0,70/2/180*pi); phi2=linspace(70/2/180*pi,70/180*pi);phi0=70/180*pi;s1=2*h*(phi1/phi0).2;s2=h-2*h*(phi0-phi2).2/(phi0).2;v1=4*h*omiga/(

19、deg2rad(70)2*phi1;a1=4*h*omiga2/(phi0)2+0*phi1;v2=4*h*omiga*(phi0-phi2)/(deg2rad(70)2;a2=-4*h*omiga2/(phi0)2+0*phi2;phi3=linspace(deg2rad(70),deg2rad(170); s3=h+0*phi3;v3=0+0*phi3;a3=0+0*phi3;phi4=linspace(deg2rad(170),deg2rad(240); s4=h*(1-(phi4-phi0-deg2rad(100)/deg2rad(70)+1/2/pi*sin(2*pi*(phi4-p

20、hi0-deg2rad(100)/deg2rad(70);v4=-h*omiga/deg2rad(70)*(1-cos(2*pi/deg2rad(70)*(phi4-deg2rad(70)-deg2rad(100);a4=-2*pi*h*omiga2/(deg2rad(70)2*sin(2*pi/deg2rad(70)*(phi4-deg2rad(70)-deg2rad(100);phi5=linspace(deg2rad(240),deg2rad(360); %½üÐݳÌs5=0+0*phi5;v5=0+0*phi5;a5=0+0

21、*phi5;phi=phi1,phi2,phi3,phi4,phi5;s=s1,s2,s3,s4,s5;v=v1,v2,v3,v4,v5;a=a1,a2,a3,a4,a5;Dx0Dphi=(v./omiga-N2*N1*e).*sin(N1*phi.*pi/180)+N1*(s0+s).*cos(N1*phi.*pi/180);Dy0Dphi=(v./omiga-N2*N1*e).*cos(N1*phi.*pi/180)-N1*(s0+s).*sin(N1*phi.*pi/180);phi=phi1,phi2,phi3,phi4,phi5;s=s1,s2,s3,s4,s5;v=v1,v2,v3

22、,v4,v5;a=a1,a2,a3,a4,a5;DsDphi=v./omiga;DDsDphi=(a./omiga2);DDxDphi=(DDsDphi-s0-s).*sin(N1*phi.*pi/180)+(2*DsDphi-N1*N2*e).*N1.*cos(N1*phi.*pi/180);DDyDphi=(DDsDphi-s0-s).*cos(N1*phi.*pi/180)-(2*DsDphi-N1*N2*e).*N1.*sin(N1*phi.*pi/180);A=(Dx0Dphi.2+Dy0Dphi.2).1.5;B=Dx0Dphi.*DDyDphi-Dy0Dphi.*DDxDphi;

23、rou=A./B; %理論輪廓曲率半徑rou0=abs(rou+M*rt);% 實(shí)際輪廓曲率半徑plot(phi,rou0)axis(1,7,0,100)xlabel('凸輪轉(zhuǎn)角/rad')ylabel('曲率半徑/mm')grid ontitle('曲率半徑曲線')7>輪廓線h=30; r0=40.3; e=20; N1=1;% 凸輪逆時(shí)針轉(zhuǎn) N2=1;%偏距為正rt=10;%滾子半徑 n=1; M=-1;omiga=2*pi*n/60;s0=sqrt(r02-e2);phi1=linspace(0,70/2/180*pi); %推程ph

24、i2=linspace(70/2/180*pi,70/180*pi);phi0=70/180*pi;s1=2*h*(phi1/phi0).2;s2=h-2*h*(phi0-phi2).2/(phi0).2;v1=4*h*omiga/(deg2rad(70)2*phi1;a1=4*h*omiga2/(phi0)2+0*phi1;v2=4*h*omiga*(phi0-phi2)/(deg2rad(70)2;a2=-4*h*omiga2/(phi0)2+0*phi2;phi3=linspace(deg2rad(70),deg2rad(170); %遠(yuǎn)休程s3=h+0*phi3;v3=0+0*phi3

25、;a3=0+0*phi3;phi4=linspace(deg2rad(170),deg2rad(240); %回程s4=h*(1-(phi4-phi0-deg2rad(100)/deg2rad(70)+1/2/pi*sin(2*pi*(phi4-phi0-deg2rad(100)/deg2rad(70);v4=-h*omiga/deg2rad(70)*(1-cos(2*pi/deg2rad(70)*(phi4-deg2rad(70)-deg2rad(100);a4=-2*pi*h*omiga2/(deg2rad(70)2*sin(2*pi/deg2rad(70)*(phi4-deg2rad(70)-deg2rad(100); phi5=linspace(deg2rad(240),deg2rad(360); %近休程s5=0+0*phi5;v5=0+0*phi5;a5=0+0*phi5;phi=phi1,phi2,phi3,ph