機械原理課程設計報告 偏置直動滾子推桿盤形凸輪

-.z.凸輪大作業(yè)選題：凸輪5-C一：題目及原始數(shù)據(jù)：利用計算機輔助設計完成偏置直動滾子推桿盤形凸輪機構的設計，數(shù)據(jù)如下所示，凸輪沿逆時針方向做勻速運動。具體要求如下：1.推程運動規(guī)律為等加速等減速運動，回程運動規(guī)律為五次多項式運動；2.近休凸輪轉角為0°-30°；推程凸輪轉角30°-210°；遠休凸輪轉角210°-280°；回程凸輪轉角280°-360。°3.初選基圓半徑為22mm；4.偏距為+14mm5.滾子半徑為18mm6.推桿行程為35mm7.許用壓力角為α1=35°，α2=65°。8.最小曲率半徑為0.35rr9.計算點數(shù)取120.二：推桿運動規(guī)律及凸輪輪廓線方程1.推程加速階段：s1=70.*a1.*a1/pi/pi;*1=(s0+s1).*sin(a1)+e*cos(a1);y1=(s0+s1).*cos(a1)-e*sin(a1);k1=140*a1/pi^2;i1=[(k1-e).*sin(a1)+(s0+s1).*cos(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);j1=[-(k1-e).*cos(a1)+(s0+s1).*sin(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);*10=*1-18*j1;y10=y1-18*i1;2.推程減速階段：s11=35-70.*(pi-a11).*(pi-a11)/pi/pi;*11=(s0+s11).*sin(a11)+e*cos(a11);y11=(s0+s11).*cos(a11)-e*sin(a11);k11=140.*(pi-a1)/pi^2;i11=[(k11-e).*sin(a11)+(s0+s11).*cos(a11)].*[(k11-e).*(k11-e)+(s0+s11).*(s0+s11)].^(-1/2);j11=[-(k11-e).*cos(a11)+(s0+s11).*sin(a11)].*[(k11-e).*(k11-e)+(s0+s11).*(s0+s11)].^(-1/2);*101=*11-18*j11;y101=y11-18*i11;3.遠休階段：*2=(s0+s2).*sin(a2)+e*cos(a2);y2=(s0+s2).*cos(a2)-e*sin(a2);k2=0;i2=[(k2-e).*sin(a2)+(s0+s2).*cos(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);j2=[-(k2-e).*cos(a2)+(s0+s2).*sin(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);*20=*2-18*j2;y20=y2-18*i2;4.回程階段：a30=33*pi/18-a3;s3=(350.*a30.*a30.*a30/(4*pi/9)^3-525.*a30.*a30.*a30.*a30/(4*pi/9)^4+210.*a30.*a30.*a30.*a30.*a30/(4*pi/9)^5);*3=(s0+s3).*sin(a3)+e*cos(a3);y3=(s0+s3).*cos(a3)-e*sin(a3);k3=-1050.*a30.*a30/(4*pi/9)^3+2100.*a30.*a30.*a30/(4*pi/9)^4-1050.*a30.*a30.*a30.*a30/(4*pi/9)^5;i3=[(k3-e).*sin(a3)+(s0+s3).*cos(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);j3=[-(k3-e).*cos(a3)+(s0+s3).*sin(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);*30=*3-18*j3;y30=y3-18*i3;5.近休階段：s4=0;*4=(s0+s4).*sin(a4)+e*cos(a4);y4=(s0+s4).*cos(a4)-e*sin(a4);k4=0;i4=[(k4-e).*sin(a4)+(s0+s4).*cos(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);j4=[-(k4-e).*cos(a4)+(s0+s4).*sin(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);*40=*4-18*j4;y40=y4-18*i4;三：matlab計算程序e=14;r0=22s0=sqrt(r0^2-e^2);c111=1;c3=1;for(i=1:0.1:200)%由壓力角條件循環(huán)求適宜基圓半徑if(c111<35/180)&&(c3<65/180)%判斷條件elsea1=0:pi/60:pi/2;%推程加速階段s1=70.*a1.*a1/pi/pi;*1=(s0+s1).*sin(a1)+e*cos(a1);y1=(s0+s1).*cos(a1)-e*sin(a1);k1=140*a1/pi^2;%對s1求導i1=[(k1-e).*sin(a1)+(s0+s1).*cos(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);j1=[-(k1-e).*cos(a1)+(s0+s1).*sin(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);*10=*1-18*j1;y10=y1-18*i1;a11=pi/2:pi/60:pi;%推程減速階段s11=35-70.*(pi-a11).*(pi-a11)/pi/pi;*11=(s0+s11).*sin(a11)+e*cos(a11);y11=(s0+s11).*cos(a11)-e*sin(a11);k11=140.*(pi-a1)/pi^2;i11=[(k11-e).*sin(a11)+(s0+s11).*cos(a11)].*[(k11-e).*(k11-e)+(s0+s11).*(s0+s11)].^(-1/2);j11=[-(k11-e).*cos(a11)+(s0+s11).*sin(a11)].*[(k11-e).*(k11-e)+(s0+s11).*(s0+s11)].^(-1/2);*101=*11-18*j11;y101=y11-18*i11;a2=pi:pi/60:25*pi/18;%凸輪遠休階段s2=35;%推桿行程*2=(s0+s2).*sin(a2)+e*cos(a2);y2=(s0+s2).*cos(a2)-e*sin(a2);k2=0;i2=[(k2-e).*sin(a2)+(s0+s2).*cos(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);j2=[-(k2-e).*cos(a2)+(s0+s2).*sin(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);*20=*2-18*j2;y20=y2-18*i2;a3=25*pi/18:pi/60:33*pi/18;%推桿回程階段a30=33*pi/18-a3;s3=(350.*a30.*a30.*a30/(4*pi/9)^3-525.*a30.*a30.*a30.*a30/(4*pi/9)^4+210.*a30.*a30.*a30.*a30.*a30/(4*pi/9)^5);*3=(s0+s3).*sin(a3)+e*cos(a3);y3=(s0+s3).*cos(a3)-e*sin(a3);k3=-1050.*a30.*a30/(4*pi/9)^3+2100.*a30.*a30.*a30/(4*pi/9)^4-1050.*a30.*a30.*a30.*a30/(4*pi/9)^5;i3=[(k3-e).*sin(a3)+(s0+s3).*cos(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);j3=[-(k3-e).*cos(a3)+(s0+s3).*sin(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);*30=*3-18*j3;y30=y3-18*i3;a4=33*pi/18:pi/60:2*pi;%推桿近休階段s4=0;*4=(s0+s4).*sin(a4)+e*cos(a4);y4=(s0+s4).*cos(a4)-e*sin(a4);k4=0;i4=[(k4-e).*sin(a4)+(s0+s4).*cos(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);j4=[-(k4-e).*cos(a4)+(s0+s4).*sin(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);*40=*4-18*j4;y40=y4-18*i4;plot(*10,y10,'-g*',*101,y101,'-g+',*20,y20,'-r*',*30,y30,'-b*',*40,y40,'-k*',*1,y1,'-g*',*11,y11,'-g+',*2,y2,'-r*',*3,y3,'-b*',*4,y4,'-k*')%凸輪輪廓曲線繪制title(‘凸輪輪廓曲線繪制');*label('Variable*');%*軸ylabel('VariableY');%Y軸te*t(-250,-200,'工作廓線')%文字標注te*t(100,-100,'理論廓線');gridon%加網(wǎng)格a*isequal%坐標相等a1=0:pi/60:pi/2;%壓力角計算force1=abs(atan((k1).*(s0+s1).^-1));c1=ma*(force1);a11=pi/2:pi/60:pi;force11=abs(atan((k11).*(s0+s11).^-1));c11=ma*(force11);c111=ma*(c1,c11);a2=pi:pi/60:25*pi/18;force2=abs(atan((k2).*(s0+s2).^-1));c2=ma*(force2);a3=25*pi/18:pi/60:33*pi/18;force3=abs(atan((k3).*(s0+s3).^-1));c3=ma*(force3);a4=33*pi/18:pi/60:2*pi;force4=abs(atan((k4).*(s0+s4).^-1));c4=ma*(force4);r0=r0+1;%每循環(huán)一次基圓半徑+1s0=sqrt(r0^2-e^2);endend%求最大壓力角位置c111c3[m1,n1]=sort(force1);bend1=n1(end-1+1:end)jiaodu1=(bend1(end)-1)*3[m11,n11]=sort(force11);bend11=n11(end-1+1:end)jiaodu11=(bend11(end)-1)*3+90[m3,n3]=sort(force3);bend3=n3(end-1+1:end)jiaodu=(bend3(end)-1)*3+250vv=i1./j1;vv1=diff(vv);vv2=diff(vv,2);vv22=[0,vv2];p=(1+vv1.^2).^(3/2)./vv22;g=min(p)vv1=i11./j11;vv11=diff(vv1);vv21=diff(vv1,2);vv221=[0,vv21];p1=(1+vv11.^2).^(3/2)./vv221;g1=min(p1)vv2=i2./j2;vv222=diff(vv2);vv223=diff(vv2,2);vv2211=[0,vv223];p2=(1+vv222.^2).^(3/2)./vv2211;g2=min(p2)vv=i3./j3;vv1=diff(vv);vv2=diff(vv,2);vv22=[0,vv2];p=(1+vv1.^2).^(3/2)./vv22;g3=min(p)vv=i4./j4;vv1=diff(vv);vv2=diff(vv,2);v