機械原理課程設(shè)計 偏置直動滾子推桿盤形凸輪

1、凸輪大作業(yè)選題：凸輪5-C一：題目及原始數(shù)據(jù)：利用計算機輔助設(shè)計完成偏置直動滾子推桿盤形凸輪機構(gòu)的設(shè)計，已知數(shù)據(jù)如下所示，凸輪沿逆時針方向做勻速運動。具體要求如下：1.推程運動規(guī)律為等加速等減速運動，回程運動規(guī)律為五次多項式運動；2.近休凸輪轉(zhuǎn)角為0°-30°；推程凸輪轉(zhuǎn)角30°-210°；遠(yuǎn)休凸輪轉(zhuǎn)角210°-280°；回程凸輪轉(zhuǎn)角280°-360。°3.初選基圓半徑為22mm；4.偏距為+14mm5.滾子半徑為18mm6.推桿行程為35mm7.許用壓力角為1=35°，2=65°。8.最小曲率

2、半徑為0.35rr9.計算點數(shù)取120.二：推桿運動規(guī)律及凸輪輪廓線方程1.推程加速階段：s1=70.*a1.*a1/pi/pi; x1=(s0+s1).*sin(a1)+e*cos(a1); y1=(s0+s1).*cos(a1)-e*sin(a1); k1=140*a1/pi2; 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); x1

3、0=x1-18*j1; y10=y1-18*i1;2.推程減速階段：s11=35-70.*(pi-a11).*(pi-a11)/pi/pi; x11=(s0+s11).*sin(a11)+e*cos(a11); y11=(s0+s11).*cos(a11)-e*sin(a11); k11=140.*(pi-a1)/pi2; 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).

4、*(k11-e)+(s0+s11).*(s0+s11).(-1/2); x101=x11-18*j11; y101=y11-18*i11;3.遠(yuǎn)休階段：x2=(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); x

5、20=x2-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); x3=(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.*

6、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); x30=x3-18*j3; y30=y3-18*i3;5.近休階段：s4=0;x4=(s0+s4).*sin(a4)+e*cos(a4);y4=(s0+s4).*cos(a4)-e*sin(a4);k4=0;i4=(k4-e).*sin(a4)+(s0+s4)

7、.*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); x40=x4-18*j4; y40=y4-18*i4;三：matlab計算程序e=14;r0=22s0=sqrt(r02-e2);c111=1;c3=1; for(i=1:0.1:200) %由壓力角條件循環(huán)求合適基圓半徑 if (c111<35/180)&&(c3<65/180) %判斷條件 else a1=0

8、:pi/60:pi/2; % 推程加速階段 s1=70.*a1.*a1/pi/pi; x1=(s0+s1).*sin(a1)+e*cos(a1); y1=(s0+s1).*cos(a1)-e*sin(a1); k1=140*a1/pi2; % 對s1求導(dǎo) 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); x10=x1-18*j1; y10

9、=y1-18*i1; a11=pi/2:pi/60:pi; % 推程減速階段 s11=35-70.*(pi-a11).*(pi-a11)/pi/pi; x11=(s0+s11).*sin(a11)+e*cos(a11); y11=(s0+s11).*cos(a11)-e*sin(a11); k11=140.*(pi-a1)/pi2; 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).*(k1

10、1-e).*(k11-e)+(s0+s11).*(s0+s11).(-1/2); x101=x11-18*j11; y101=y11-18*i11; a2=pi:pi/60:25*pi/18; %凸輪遠(yuǎn)休階段 s2=35; %推桿行程 x2=(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).

11、*(k2-e).*(k2-e)+(s0+s2).*(s0+s2).(-1/2); x20=x2-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); x3=(s0+s3).*sin(a3)+e*cos(a3); y3=(s0+s3).*cos(a3)-e*sin(a3); k3=-1050

12、.*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); x30=x3-18*j3; y30=y3-18*i3; a4=33*pi/18:pi/60:2*pi; %推桿近休階

13、段 s4=0; x4=(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); x40=x4-18*j4; y40=y4-18*i4; plot(x10,y10,'-g*',x101,y101,'

14、;-g+',x20,y20,'-r*',x30,y30,'-b*',x40,y40,'-k*',x1,y1,'-g*',x11,y11,'-g+',x2,y2,'-r*',x3,y3,'-b*',x4,y4,'-k*') %凸輪輪廓曲線繪制title(凸輪輪廓曲線繪制');xlabel('Variable X'); %X軸ylabel('Variable Y'); %Y軸text(-250,-200,'工作廓線&

15、#39;) %文字標(biāo)注text(100,-100,'理論廓線'); grid on %加網(wǎng)格axis equal %坐標(biāo)相等 a1=0:pi/60:pi/2; %壓力角計算force1=abs(atan(k1).*(s0+s1).-1); c1=max(force1); a11=pi/2:pi/60:pi; force11=abs(atan(k11).*(s0+s11).-1);c11=max(force11);c111=max(c1,c11);a2=pi:pi/60:25*pi/18;force2=abs(atan(k2).*(s0+s2).-1);c2=max(force2

16、); a3=25*pi/18:pi/60:33*pi/18; force3=abs(atan(k3).*(s0+s3).-1);c3=max(force3); a4=33*pi/18:pi/60:2*pi;force4=abs(atan(k4).*(s0+s4).-1);c4=max(force4); r0=r0+1; %每循環(huán)一次基圓半徑+1 s0=sqrt(r02-e2); end end %求最大壓力角位置 c111 c3 m1,n1=sort(force1); bend1=n1(end-1+1:end) jiaodu1=(bend1(end)-1)*3 m11,n11=sort(for

17、ce11); 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+250 vv=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

18、=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);vv22=0,vv2;p=(1+vv1.2).(3/2)./vv22;g4=min(p) ) ro=sqrt(s02+e2) %求合適基圓半

19、徑x1y1x11y11x2y2x3y3x4y4x10y10x101y101x20y20x30y30x40y40force1force11force2force3force4c111c3jiaodu11jiaodu3四：計算結(jié)果及分析：（1）工作廓線坐標(biāo)：1;推程加速階段X坐標(biāo)1424.9479312972938 35.833579542701846.6331664128873 57.323071617132667.8798591305935 182.841407726263189.240951783976 205.709540955551 210.235796832989 214.255618

20、712129217.755304716715 225.009303250807226.308726472013 227.0328136593412推程加速階段Y坐標(biāo) 209.532813659341 208.532370629552206.998924199561 204.935879972465202.347567871334 191.488076931747186.860293712014 181.742572281261 170.077262819810163.551861181761 156.580882432274 106.20142493273596.5309400655771 8

21、6.5447839816669 76.2628805759777 65.705931370278654.8953994088027 43.8534915206956 3推程減速階段X坐標(biāo)227.032813659341 225.331644205755223.421214513810 220.838411540547 203.884084407342 198.035875547330191.570525763878184.505401336072 176.860074548737 140.941425726738 97.8264479001947 49.6409729459048 37.082

22、4807009210 24.398347089867511.6292027035831 4推程減速階段Y坐標(biāo)5遠(yuǎn)休階段X坐標(biāo)6遠(yuǎn)休階段Y坐標(biāo)7減速階段X坐標(biāo)8減速階段Y坐標(biāo) 42.7749931638176 54.5773610652433 65.978032703814576.9348117792082 106.921098432862115.946980634777 132.644871079527 140.375991577690147.738568287212 154.762440563530 161.470429890068 167.874092062809173.969459073

23、374 179.732900602872 9近休X坐標(biāo)10近休Y坐標(biāo)1 188.460739555420193.050970998372 200.632884834880203.603785729623 209.978359607845209.532813659341（2）理論廓線坐標(biāo);推程加速階段X坐標(biāo)12.8000000000000 22.872827317908432.8885297342726 62.3761432336542 71.9473900632018 81.354400741436490.5763518517692 99.5926727370081 108.383068927

24、807 116.927546041764125.206434225271 133.200413203958140.890538000620 148.258265371979155.285481007489 161.954527524708168.248233286539 193.523452696839197.235553228942 200.468308129584 203.209568133643 208.374821351328 209.044777489697推程加速階段Y坐標(biāo)191.572858202826 190.652383395923189.241484022951 187.3

25、43250393275184.961604006512 174.967781050798170.708173911217 165.997083129187160.843331597179 155.256545793692149.247151954669 142.826372082411136.006219708068 113.281436600464104.999590134052 96.390134222989187.4697165758390 78.2557315040260 68.766309983443159.020308476717649.0372963926355 38.8375420595755 -3.69999384635310 -14.6561667163081推程減速階段X坐標(biāo) 207.342473187261205.495287697520 203.031928716014 199.955103403567196.269775760348191.9