课程设计凸轮程序.docVIP

华南理工大学机械原理课程设计凸轮程序用MATLAB编程，e为偏距，r0为基圆半径，h为升程e=5;r0=55;h=8s0=(r02-e2)0.5;a1=linspace(0,23/18*pi,10000); %近休止角变量s1=0;y1=cos(a1).*(s1+s0)-e.*sin(a1); %近休止理论轮廓方程x1=sin(a1).*(s1+s0)+e.*cos(a1);x11=(s0-5)*sin(a1)+e*cos(a1); %近休止工作轮廓方程y11=(s0-5)*cos(a1)-e*sin(a1);a2=linspace(23/18*pi,29/18*pi,10000); %推程角变量s2=h*(a2-23/18*pi)./(pi/3)-sin(6*(a2-23/18*pi)./(2*pi); %正弦加速度推程运动方程y2=cos(a2).*(s2+s0)-e.*sin(a2); %推程理论轮廓方程x2=sin(a2).*(s2+s0)+e.*cos(a2); %推程理论轮廓方程x20=(h.*(1-cos(2.*pi.*(a2*180/pi-230)/60)/60)-e).*sin(a2)+(s0+s2).*cos(a2); %（x2）对（a2)求导y20=(h.*(1-cos(2.*pi.*(a2*180/pi-230)/60)/60)-e).*cos(a2)-(s0+s2).*sin(a2); %（y2)对（a2）求导x21=x2+5.*y20./sqrt(x20.2+y20.2); %推程的工作轮廓方程y21=y2-5.*x20./sqrt(x20.2+y20.2); %推程的工作轮廓方程a3=linspace(29/18*pi,5/3*pi,10000); %远休止角变量s3=h; y3=cos(a3).*(s3+s0)-e.*sin(a3); %远休止理论轮廓方程x3=sin(a3).*(s3+s0)+e.*cos(a3);x31=(s3+s0-5)*sin(a3)+e*cos(a3); %远休止工作轮廓方程y31=(s3+s0-5)*cos(a3)-e*sin(a3);a4=linspace(5/3*pi,2*pi,10000); %回程角变量s4=h*(1-(3*a4-5*pi)./(pi)+sin(6*a4-10*pi)./(2*pi); %正弦加速度回程运动方程y4=cos(a4).*(s4+s0)-e.*sin(a4); %回程理论轮廓方程x4=sin(a4).*(s4+s0)+e.*cos(a4);x40=(h.*(cos(a4.*180/pi-300)/60*2*pi)-1)/60-e).*sin(a4)+(s0+s4).*cos(a4); %理论轮廓方程对角变量的求导y40=(h.*(cos(a4.*180/pi-300)/60*2*pi)-1)/60-e).*cos(a4)-(s0+s4).*sin(a4);x41=x4+5.*y40./sqrt(x40.2+y40.2); %回程工作轮廓方程y41=y4-5.*x40./sqrt(x40.2+y40.2);a0=linspace(0,2*pi,10000); %基圆角变量x0=s0.*sin(a0)+e.*cos(a0); %基圆方程y0=s0.*cos(a0)-e.*sin(a0);plot(x0,y0,k.);hold onplot(x1,y1,g.);grid onplot(x2,y2,r.);plot(x3,y3,g.);plot(