机械原理 凸轮机构设计

1、机械原理课程设计凸轮机构设计（一）目录·········································1（一）、题目及原始数据····

2、;······················2（二）、推杆运动规律及凸轮廓线方程··············3（三）、计算程序方框图·········

3、;·················5（四）、计算源程序······························6（五）、程序

4、计算结果及分析······················10（六）、凸轮机构图·························

5、·····15（七）、心得体会································16（八）、参考书··········

6、;························16（一）、题目及原始数据试用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计，凸轮以1rad/s的角速度沿逆时针方向转动。要求：（1）、推程运动规律为等加速等减速运动，回程运动规律为五次多项式运动规律；（2）、打印出原始数据；（3）、打印出理论轮廓和实际轮廓的坐标值；（4）、打印出推程和回程的最大压力角

7、，以及出现最大压力角时凸轮的相应转角；（5）、打印出凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角；（6）、打印出凸轮运动的位移；（7）、打印最后所确定的凸轮的基圆半径。原始数据如下：r0=0.015; 初选的基圆半径r0Deltar0=0.0005; 当许用压力角或许用最小曲率半径不满足时，r0以r0为步长增加重新计算rr=0.010; 滚子半径rrh=0.028; 推杆行程he=0.005; 偏距eomega=1; 原动件凸轮运动角速度，逆时针delta1=pi/3; 近休止角1delta2=2*pi/3; 推程运动角2delta3=pi/2; 远休止角3delta4=pi/2; 回程

8、运动角4alpha1=pi/6; 推程许用压力角1alpha2=(70/180)*pi; 回程许用压力角2rho0min=0.3*rr; 许用最小曲率半径 min（二）、推杆运动规律及凸轮廓线方程推杆运动规律：（1） 近休阶段：0o<60 os=0v=0a=0（2） 推程阶段：60o<180 o等加速运动规律：60o<120 os=2h(-60o)2/(120 o)2v=4h(-60o)/(120 o)2a=4h2/(120 o)2等减速运动规律：120o<180 os=h-2h(120o -(-60o)2/(120 o)2 v=4h(120o -(-60o)/(120

9、 o)2 a=-4h2/(120 o)2（3） 远休阶段：180o<270 o s=hv=0a=0（4） 回程阶段：270o360 o 五次多项式运动规律：s=h-(10h(-270o)3/(90 o)3-15h(-270o)4/(90 o)4+6h( v=-(30h(-270o)2/(90 o)3-60h(-270o)3/(90 o4)+30h(-270o)4/(90 o)5)2o o32o2 a=-(60h(-270)/(90)-180h(-270)/(90o4)+120h2(-270o)3/(90 o)5)凸轮廓线方程：（1） 理论廓线方程：s0=sqrt(r02-e2)x=(s0

10、+s)siny=(s0+s)cos+ecos-esin-270o)5/(90 o)5)(2) 实际廓线方程先求x，y的一阶导数x=(v/-e) siny=(v/-e) cos+(s0+s)cos-(s0+s)sin再求sin，cossinx/sqrt(x)+(y) 22cos-y/sqrt(x)+(y) 22最后求实际廓线方程x1=x-rrcosy1=y-rrsin压力角方程：曲率半径计算公式：（四）、计算源程序%凸轮机构大作业 Matlab语言源程序%选题：偏置直动滚子推杆盘形凸轮机构 5A % 推程运动规律：等加速等减速运动% 回程运动规律：五次多项式运动% 作者:WYH 学号：xxxxx

11、xxx 日期：2007.12.26%clear all;%close all;clc;%赋初值r0=0.015;Deltar0=0.0005;rr=0.010;h=0.028;e=0.005;omega=1; %原动件凸轮运动角速度，逆时针delta1=pi/3; %近休止角delta2=2*pi/3; %推程运动角delta3=pi/2; %远休止角delta4=pi/2; %回程运动角alpha1=pi/6; %推程许用压力角alpha2=(70/180)*pi; %回程许用压力角rho0min=0.3*rr; %许用最小曲率半径del1=delta1+delta2;del2=del1+d

12、elta3;temp=0; %判断是否执行r0r0Deltar0的变量while (temp=0)temp=1;s0=sqrt(r02-e2); %求解 s0alpha1max=0;delta1max=0; %定义alpha1的最大值以及对应的delta1值alpha2max=0;delta2max=0; %定义alpha2的最大值以及对应的delta2值rhoamin=r0-rr;deltamin=0; %定义rhoa的最小值以及对应的delta值for I=0:120; %圆周120等分delta=(I*3/180)*pi;if delta>=0&delta<delta

13、1 %近休阶段s=0; %位移v=0; %速度a=0; %加速度elseif delta>=delta1&delta<(delta2/2)+delta1 %等加速推程 s=2*h*(delta-delta1)2/delta22;v=4*h*omega*(delta-delta1)/delta22;a=4*h*omega2/delta22;elseif delta>=(delta2/2)+delta1&delta<del1 %等减速推程 s=h-2*h*(delta2-(delta-delta1)2/delta22;v=4*h*omega*(delta2-

14、(delta-delta1)/delta22;a=-4*h*omega2/delta22;elseif delta>=del1&delta<del2 %远休阶段s=h;v=0;a=0;elseif delta>=del2&delta<=2*pi %五次多项式运动规律回程s=h-(10*h*(delta-del2)3/delta33-15*h*(delta-del2)4/delta34+6*h*(delta-del2)5/delta35);v=-(30*h*omega*(delta-del2)2/delta43-60*h*omega*(delta-del2

15、)3/delta44+30*h*omega*(delta-del2)4/delta45);a=-(60*h*omega2*(delta-del2)/delta43-180*h*omega2*(delta-del2)2/delta44+120*h*omega*(delta-del2)3/delta45);endx=(s0+s)*sin(delta)+e*cos(delta); %理论轮廓方程式y=(s0+s)*cos(delta)-e*sin(delta);x_=(v/omega-e)*sin(delta)+(s0+s)*cos(delta); %理论轮廓对delta求一次导数y_=(v/ome

16、ga-e)*cos(delta)-(s0+s)*sin(delta);x_=(a/omega2-(s0+s)*sin(delta)+(2*v/omega-e)*cos(delta); %理论轮廓对delta求二次导数y_=(a/omega2-(s0+s)*cos(delta)-(2*v/omega-e)*sin(delta);x1=x-rr*(-y_/sqrt(x_2+y_2); %实际轮廓方程式y1=y-rr*(x_/sqrt(x_2+y_2);alpha=atan(v-e)/(sqrt(r02-e2)+s); %求压力角if delta>=del2&delta<=2*p

17、i %判断是否为回程if abs(alpha)>alpha2 %判断是否大于回程许用压力角 r0=r0+Deltar0;temp=0;break;elseif abs(alpha)>alpha2max %满足许用压力角，则找出回程最大压力角alpha2max=abs(alpha);delta2max=delta;endendelseif abs(alpha)>alpha1 %判断是否大于推程许用压力角r0=r0+Deltar0; %不满足许用压力角，则增大基圆半径重新计算temp=0;break;elseif abs(alpha)>alpha1max %满足许用压力角，

18、则找出推程最大压力角alpha1max=abs(alpha);delta1max=delta;endendendrho=(x_2+y_2)(3/2)/(x_*y_-y_*x_); %计算曲率半径if rho<0rhoa=abs(rho)-rr;if rhoa>=rho0min %满足最小曲率半径if rhoa<rhoamin %找出实际轮廓曲线的最小曲率半径及其对应的delta角rhoamin=rhoa;deltamin=delta;endelser0=r0+Deltar0;temp=0;break;endendDelta(I+1)=(delta/pi)*180; %delt

19、a由弧度值转化为角度值 X(I+1)=x*1000;Y(I+1)=y*1000;X1(I+1)=x1*1000;Y1(I+1)=y1*1000;S(I+1)=s;V(I+1)=v;A(I+1)=a;ALPHA(I+1)=(alpha/pi)*180;PHO(I+1)=rho*1000;endenddeltamin=(deltamin/pi)*180;alpha1max=(alpha1max/pi)*180;delta1max=(delta1max/pi)*180;alpha2max=(alpha2max/pi)*180;delta2max=(delta2max/pi)*180;figure(1

20、);axis equal;hold ont=0:0.01:2*pi;xx=r0*cos(t)*1000;yy=r0*sin(t)*1000;xxx=(rr*cos(t)+X(1)/1000)*1000;yyy=(rr*sin(t)+Y(1)/1000)*1000;xxxx=e*cos(t)*1000;yyyy=e*sin(t)*1000;plot(xx,yy,'m-',X,Y,':',X1,Y1,'k',xxx,yyy,'c-',xxxx,yyyy,'y-');%画出理论轮廓及实际轮廓以及基圆legend('

21、;基圆','理论轮廓','实际工作轮廓');plot(0,0,'ko')plot(X(1),Y(1),'ko');title('凸轮轮廓曲线图');xlabel('X/mm');ylabel('Y/mm');figure(2);plot(Delta,S,Delta,V,'r-',Delta,A,'k:'); %画出位移、速度、加速度曲线图title('凸轮运动规律曲线图');xlabel('delta/(o)');ylabel('s/m v/m.s-1 a/m.s-2');legend('位移','速度','加速度');%结果显示：disp(num2str(Delta'),num2str(X&#