《基于MATLAB的凸轮设计.docx》由会员分享，可在线阅读，更多相关《基于MATLAB的凸轮设计.docx（4页珍藏版）》请在三一文库上搜索。
1、精心整理中国地质大学（武汉）1 .凸轮要求设计一对心直动滚子推杆盘形凸轮机构，滚子半径r二10mm凸轮以等角速度逆时针回转凸轮转角5=0120时，推杆等速上升20mm 8=120180时，推杆远休止；5=180270时，推杆等加速等减速下降 20mm江270360时，推杆近休止。要求推程的最大压力角 二30,试选取合适的基圆半径，并绘制凸轮的廓线。问此凸轮是否有缺陷，应如何补救。2 .列出凸轮运动方程-30I , .0< <2/3S * 6 nV =打1 a = 0|S = 202/3丫<S=- 140 + * 6-320320普6A2TF"-320/X24G0 3
2、60 *njr 480it320-*6A2Ti* 55/4tt <d <q=-320/tT227T3 .由方程写MATLA蹒程序%i.已知参数clear;r0=50; %基圆半径rr=10;%滚子半径h=20; %亍程delta01=120;%推程运动角delta02=60;%远休角delta03=90;%回程运动角hd=pi/180;du=180/pi;n1=delta01+delta02;n2=delta01+delta02+delta03;%2凸轮曲线设计n=360;for i=1:360%计算推杆运动规律if i<=delta01s(i)=30/pi*(i*hd);d
3、s(i)=30/pi;ds=ds(i);elseif i>delta01&&i<=n1;s(i)=h；ds(i)=0;ds=ds(i);elseif i>n1&&i<=(n1+delta03/2)s(i)=-140+320/pi*(i*hd)-160/piA2*(i*hd)A2;ds(i)=320/pi-320/piA2*(i*hd);ds=ds(i);elseif i>(n1+delta03/2)&&i<=n2s(i)=360-480/pi*(i*hd)+160/piA2*(i*hd)A2;ds(i)=-48
4、0/pi+320/piA2*(i*hd);ds=ds(i);elseif i>n2&&i<=ns(i)=0；ds=0;end% 计算凸轮轨迹曲线xx(i)=(r0+s(i)*sin(i*hd);%计算理论轮廓曲线yy(i)=(r0+s(i)*cos(i*hd);dx(i)=ds*sin(i*hd)+(r0+s(i)*cos(i*hd);%算导数dy(i)=ds*cos(i*hd)-(r0+s(i)*sin(i*hd);xp(i)=xx(i)+rr*dy(i)/sqrt(dx(i)A2+dy(i)A2);yp(i)=yy(i)-rr*dx(i)/sqrt(dx(i)A
5、2+dy(i)A2);end%3.输出凸轮轮廓曲线figure(1);hold on ;grid on;axis equal ;axis(-(r0+h-30)(r0+h+10)-(r0+h+10)(r0+rr+10);text(r0+h+3,4,'X');text(3,r0+rr+3,'Y');text(-6,4,'O');title('对心直动滚子推杆盘形凸轮设计)；xlabel( 'x/mm');ylabel( 'y/mm');plot(-(r0+h-40)(r0+h),00,'k')；
6、plot(00,-(r0+h)(r0+rr),'k')；plot(xx,yy,'r-'); %磔凸轮实际轮廓曲线ct=linspace(0,2*pi);plot(r0*cos(ct),r0*sin(ct),'g')；plot(rr*cos(ct),r0+rr*sin(ct),'k'plot(0,r0,'o' ); %滚子圆中心plot(00,r0r0+30,'k')；plot(xp,yp,'b' );%绘凸轮实际轮廓曲线%4.凸轮机构运动仿真% 计算凸轮滚子转角xp0=0;yp0=r
7、0-rr;dss=sqrt(diff(xp).A2+diff(yp).A2);ss(1)=sqrt(xp(1)-xp0)A2+(xp(1)-yp0)A2);for i=1:359%寸轮廓曲线进行差分计算%轮廓曲线第一点长度ss(i+1)=ss(i)+dss(i); end%计算实际廓曲线长度phi=ss/rr;% 运动仿真开始figure(2);%计算滚子转角I 二 3 / .二m=moviein(20);j=0；for i=1:360j=j+1;delta(i)=i*hd;xy=xp',yp'%凸轮转角%凸轮实际轮廓曲线坐标A1=cos(delta(i),sin(delta(
8、i);-sin(delta(i),cos(delta(i);xy=xy*A1;渝转后实际凸轮曲线坐标clf；%绘凸轮plot(xy(:,1),xy(:,2);%凸轮坐标旋转矩阵hold on ;axis equal ;axis(-(120)(470)-(100)(140);plot(-(r0+h-40)(r0+h),0, plot(00,-(r0+h)(r0+rr),'k''k'plot(r0*cos(ct),r0*sin(ct),plot(rr*cos(ct),r0+s(i)+rr*sin(ct),g');%绘凸轮水平轴);%绘凸轮垂直轴)；%绘基圆&
9、#39;k');绘滚子圆plot(0rr*cos(-phi(i),r0+s(i)r0+s(i)+rr*sin(-phi(i),%绘滚子圆标线'k')；plot(00,r0+s(i)r0+s(i)+40,%绘推杆曲线plot(1:360+r0+h,s+r0);plot(r0+h)(r0+h+360),r0rO,'k' ); %绘推杆'k')；'k')；plot(r0+h)(rO+h),r0rO+h,plot(i+r0+h,s(i)+r0,'*');title('对心直动滚子推杆盘形凸轮设计)；xlab
10、el( 'x/mm');ylable( 'y/mm');m(j)=getframe;endmovie(m);4 .运动仿真结果在MATLAB中可以看出轮廓曲线有一处缺口。应用圆弧连接起来。5 .计算结果I由于数据太多，只等间隔取了三十六组数据Xx8.971818.241127.500036.424644.685951.961557.947762.371265.000065.653964.212360.621853.623144.995135.000023.941412.15540.0000-12.0696-23.2658-32.7778-39.9163-44.3
11、549-47.1503-48.8408-49.7267-50.0000-49.2404-46.9846-43.3013-38.3022-32.1394-25.0000-17.1010-8.6824-0.0000Yy50.881750.116947.631443.409237.495930.000021.091210.99770.0000-11.5765-23.3714-35.0000-44.9夕F951-53.6231-60.6218-65.7785-68.9365-70.0000-68.4502-63.9223-56.7728-47.5704-37.2182-27.2222-17.7766-
12、8.7682-0.00008.682417.101 025.000032.139438.302243.301346.984649.240450.0000Dx52.540053.383052.406049.547444.811138.269930.064620.40199.5493-2.1723-14.3980-26.7301-44.99 51-53.6231-60.6218-65.7785-68.9365-70.0000-67.4676-60.0514-48.2845-33.0207-19.8785-12.5201-7.1415-3.19530.00008.682417.1010 25.000
13、032.139438.302243.301346.9846 49.240450.0000Dy0.4324-9.2677-19.2301-29.1094-38.5478-47.1869-54.6817-60.7129-65.0000-67.3121-67.4784-65.3964- 53.6231-44.9951-35.0000-23.9414-12.1554-0.000017.642533.901047.479957.256058.904755.638552.711750.709450.000049.240446.984643 .301338.302232.139425.000017.10108.68240.0000