機(jī)械原理課程設(shè)計(jì)-偏置直動(dòng)滾子推桿盤形凸輪

凸輪大作業(yè)選題：凸輪5-C一：題目及原始數(shù)據(jù)：利用計(jì)算機(jī)輔助設(shè)計(jì)完成偏置直動(dòng)滾子推桿盤形凸輪機(jī)構(gòu)的設(shè)計(jì)，數(shù)據(jù)如下所示，凸輪沿逆時(shí)針?lè)较蜃鰟蛩龠\(yùn)動(dòng)。具體要求如下：1.推程運(yùn)動(dòng)規(guī)律為等加速等減速運(yùn)動(dòng)，回程運(yùn)動(dòng)規(guī)律為五次多項(xiàng)式運(yùn)動(dòng)；2.近休凸輪轉(zhuǎn)角為0°-30°；推程凸輪轉(zhuǎn)角30°-210°；遠(yuǎn)休凸輪轉(zhuǎn)角210°-280°；回程凸輪轉(zhuǎn)角280°-360?！?.初選基圓半徑為22mm；4.偏距為+14mm5.滾子半徑為18mm6.推桿行程為35mm7.許用壓力角為α1=35°，α2=65°。8.最小曲率半徑為0.35rr9.計(jì)算點(diǎn)數(shù)取120.二：推桿運(yùn)動(dòng)規(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/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);x10=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)/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);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);x20=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.*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).*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計(jì)算程序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;x1=(s0+s1).*sin(a1)+e*cos(a1);y1=(s0+s1).*cos(a1)-e*sin(a1);k1=140*a1/pi^2;%對(duì)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=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)/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);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)].*[(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.*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;%推桿近休階段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,'-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('VariableX');%X軸ylabel('VariableY');%Y軸text(-250,-200,'工作廓線')%文字標(biāo)注text(100,-100,'理論廓線');gridon%加網(wǎng)格axisequal%坐標(biāo)相等a1=0:pi/60:pi/2;%壓力角計(jì)算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);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)一次基圓半徑+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);vv22=[0,vv2];p=(1+vv1.^2).^(3/2)./vv22;g4=min(p))ro=sqrt(s0^2+e^2)%求適宜基圓半徑x1y1x11y11x2y2x3y3x4y4x10y10x101y101x20y20x30y30x40y40force1force11force2force3force4c111c3jiaodu11jiaodu3四：計(jì)算結(jié)果及分析：〔1〕工作廓線坐標(biāo)：1;推程加速階段X坐標(biāo)14 24.9479312972938 35.8335795427018 46.6331664128873 57.3230716171326 67.8798591305935 182.841407726263 189.240951783976 205.709540955551 210.235796832989 214.255618712129 217.755304716715 225.009303250807 226.308726472013 227.0328136593412推程加速階段Y坐標(biāo)209.532813659341 208.532370629552 206.998924199561 204.935879972465 202.347567871334 191.488076931747 186.860293712014 181.742572281261 170.077262819810 163.551861181761 156.580882432274 106.201424932735 96.5309400655771 86.5447839816669 76.2628805759777 65.7059313702786 54.8953994088027 43.8534915206956 3推程減速階段X坐標(biāo)227.032813659341 225.331644205755 223.421214513810 220.838411540547 203.884084407342 198.035875547330 191.570525763878 184.505401336072 176.860074548737 140.941425726738 97.8264479001947 49.6409729459048 37.0824807009210 24.3983470898675 11.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.9780327038145 76.9348117792082106.921098432862 115.946980634777 132.644871079527 140.375991577690 147.738568287212 154.762440563530 161.470429890068 167.874092062809 173.969459073374 179.732900602872 9近休X坐標(biāo)10近休Y坐標(biāo)1188.460739555420 193.050970998372 200.632884834880 203.603785729623 209.978359607845 209.532813659341〔2〕理論廓線坐標(biāo);推程加速階段X坐標(biāo)12.8000000000000 22.8728273179084 32.8885297342726 62.3761432336542 71.9473900632018 81.3544007414364 90.5763518517692 99.5926727370081 108.383068927807 116.927546041764 125.206434225271 133.200413203958 140.890538000620 148.258265371979 155.285481007489 161.954527524708 168.248233286539 193.523452696839 197.235553228942 200.468308129584 203.209568133643 208.374821351328 209.044777489697推程加速階段Y坐標(biāo)191.572858202826 190.652383395923 189.241484022951 187.343250393275 184.961604006512 174.967781050798 170.708173911217 165.997083129187 160.843331597179 155.256545793692 149.247151954669 142.826372082411 136.006219708068113.281436600464 104.999590134052 96.3901342229891 87.4697165758390 78.2557315040260 68.7663099834431 59.0203084767176 49.0372963926355 38.8375420595755-3.69999384635310 -14.6561667163081推程