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

1、凸轮大作业选题：凸轮5-C一：题目及原始数据：利用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计，已知数据如下所示，凸轮沿逆时针方向做匀速运动。具体要求如下：1.推程运动规律为等加速等减速运动，回程运动规律为五次多项式运动；2.近休凸轮转角为0°-30°；推程凸轮转角30°-210°；远休凸轮转角210°-280°；回程凸轮转角280°-360。°3.初选基圆半径为22mm；4.偏距为+14mm5.滚子半径为18mm6.推杆行程为35mm7.许用压力角为1=35°，2=65°。8.最小曲率

2、半径为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/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.远休阶段：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) %由压力角条件循环求合适基圆半径 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求导 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; %凸轮远休阶段 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;) %文字标注text(100,-100,'理论廓线'); grid on %加网格axis equal %坐标相等 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; %每循环一次基圆半径+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四：计算结果及分析：（1）工作廓线坐标：1;推程加速阶段X坐标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坐标 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坐标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坐标5远休阶段X坐标6远休阶段Y坐标7减速阶段X坐标8减速阶段Y坐标 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坐标10近休Y坐标1 188.460739555420193.050970998372 200.632884834880203.603785729623 209.978359607845209.532813659341（2）理论廓线坐标;推程加速阶段X坐标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坐标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坐标 207.342473187261205.495287697520 203.031928716014 199.955103403567196.269775760348191.9