机械原理课程设计-偏置直动滚子推杆盘形凸轮

凸轮大作业选题：凸轮5-C一：题目及原始数据：利用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计，数据如下所示，凸轮沿逆时针方向做匀速运动。具体要求如下：1.推程运动规律为等加速等减速运动，回程运动规律为五次多项式运动；2.近休凸轮转角为0°-30°；推程凸轮转角30°-210°；远休凸轮转角210°-280°；回程凸轮转角280°-360。°3.初选基圆半径为22mm；4.偏距为+14mm5.滚子半径为18mm6.推杆行程为35mm7.许用压力角为α1=35°，α2=65°。8.最小曲率半径为0.35rr9.计算点数取120.二：推杆运动规律及凸轮轮廓线方程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.远休阶段：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计算程序e=14;r0=22s0=sqrt(r0^2-e^2);c111=1;c3=1;for(i=1:0.1:200)%由压力角条件循环求适宜基圆半径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;%对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);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;%凸轮远休阶段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,'工作廓线')%文字标注text(100,-100,'理论廓线');gridon%加网格axisequal%坐标相等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);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;%每循环一次基圆半径+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四：计算结果及分析：〔1〕工作廓线坐标：1;推程加速阶段X坐标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坐标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坐标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坐标5远休阶段X坐标6远休阶段Y坐标7减速阶段X坐标8减速阶段Y坐标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坐标10近休Y坐标1188.460739555420 193.050970998372 200.632884834880 203.603785729623 209.978359607845 209.532813659341〔2〕理论廓线坐标;推程加速阶段X坐标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坐标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推程