1、哈工大机械原理大作业凸轮11Harbin Institute of Technology机械原理大作业二课程名称： 机械原理 设计题目： 凸轮结构设计 院 系： 机电工程学院 班 级： 设 计 者： 学 号： 指导教师： 设计时间： 一、设计题目如图2-1所示直动从动件盘形凸轮机构，其原始参数见表2-1。从表2-1中选择一组凸轮机构的原始参数，据此设计该凸轮机构。 序号升程（mm）升程运动角（度）升程运动规律升程许用压力角（度）回程运动角（度）回程运动规律回程许用压力角（度）远休止角（度）近休止角（度）2080120正弦加速度3590正弦加速度659060二、凸轮推杆运动规律（1）、凸轮推杆升
2、程运动方程 （2）、凸轮推杆回程运动方程推杆位移图for n=1:1:360 hd=pi/180; x=n*hd; if n120&n210&n=300 s(n)=80*(1-(x-7*pi/6)/(pi/2)+1/(2*pi)*sin(4*x-14*pi/3); v(n)=-80/(pi/2)*(1-cos(4*x-14*pi/3); a(n)=-2*pi*80/(pi*pi/4)*sin(4*x-14*pi/3); else s(n)=0; v(n)=0; a(n)=0; endendn=1:360;plot(n,s)grid on推杆速度图for n=1:1:360 hd=pi/180;
3、 x=n*hd; if n120&n210&n=300 s(n)=80*(1-(x-7*pi/6)/(pi/2)+1/(2*pi)*sin(4*x-14*pi/3); v(n)=-80/(pi/2)*(1-cos(4*x-14*pi/3); a(n)=-2*pi*80/(pi*pi/4)*sin(4*x-14*pi/3); else s(n)=0; v(n)=0; a(n)=0; endendn=1:360;plot(n,v)grid on推杆加速度图：for n=1:1:360 hd=pi/180; x=n*hd; if n120&n210&n=300 s(n)=80*(1-(x-7*pi/
4、6)/(pi/2)+1/(2*pi)*sin(4*x-14*pi/3); v(n)=-80/(pi/2)*(1-cos(4*x-14*pi/3); a(n)=-2*pi*80/(pi*pi/4)*sin(4*x-14*pi/3); else s(n)=0; v(n)=0; a(n)=0; endendn=1:360;plot(n,a)grid on绘制凸轮机构的曲线for n=1:1:360 hd=pi/180; x=n*hd; if n120&n210&n=300 s(n)=80*(1-(x-7*pi/6)/(pi/2)+1/(2*pi)*sin(4*x-14*pi/3); v(n)=-80
5、/(pi/2)*(1-cos(4*x-14*pi/3); a(n)=-2*pi*80/(pi*pi/4)*sin(4*x-14*pi/3); else s(n)=0; v(n)=0; a(n)=0; endendn=1:360;plot(v,s)grid on确定基圆的半径和偏距升程图切线斜率：，为升程许用压力角回程图切线斜率：，为回程许用压力角在许用范围中，任意取一个点（0，-100）作为凸轮轴心。=100.4988，e=0，s0=100。绘制理论廓线：理论廓线方程四、滚子半径的确定由图可得最小曲率半径为基圆半径，p=77.43，故下面的计算中，取rr=15五、实际廓线的输出根据上述的公式绘
6、制出凸轮实际廓线附编程：for n=1:1:360 hd=pi/180; x=n*hd; if n120&n210&n=300 s(n)=80*(1-(x-7*pi/6)/(pi/2)+1/(2*pi)*sin(4*x-14*pi/3); v(n)=-80/(pi/2)*(1-cos(4*x-14*pi/3); a(n)=-2*pi*80/(pi*pi/4)*sin(4*x-14*pi/3); else s(n)=0; v(n)=0; a(n)=0; endendn=1:360;%plot(v,s,r)grid onhold ona1=cot(35*hd);aa1=-cot(65*hd);x=
7、-100:100;xx=-150:150;b1=-75.1 y=a1*x+b1;%plot(x,y) hold onbb1=-20 yy=aa1*xx+bb1;% plot(xx,yy) hold on %plot(0,-100,*)e=0;s0=100;r=sqrt(e2+s02);for n=1:360 x=n*hd; w(n)=r*cos(x); z(n)=r*sin(x); xxx(n)=(s0+s(n)*cos(x)-e*sin(x); yyy(n)=(s0+s(n)*sin(x)+e*cos(x);endplot(w,z,-r,linewidth,3);hold on plot(x
8、xx,yyy,linewidth,3)hold onplot(0,0,o)for n=1:360 x=n*hd; x1(n)=-(s0+s(n)*sin(x)-e*cos(x)+v(n)*cos(x); y1(n)=(s0+s(n)*cos(x)-e*sin(x)+v(n)*sin(x); x2(n)=-(s0+s(n)*cos(x)+e*sin(x)+v(n)*sin(x)-v(n)*sin(x)+a(n)*cos(x); y2(n)=-(s0+s(n)*sin(x)-e*cos(x)+v(n)*cos(x)+v(n)*cos(x)+a(n)*sin(x); p(n)=abs(x1(n)2+y1(n)2)1.5/(x1(n)*y2(n)-y1(n)*x2(n);endn=1:360;% plot(n,p)for n=1:360; x=n*hd; xr(n)=xxx(n)-15*y1(n)/(x1(n)2+y1(n)2)0.5; yr(n)=yyy(n)+15*x1(n)/(x1(n)2+y1(n)2)0.5;endn=1:360 plot(xr,yr,r,linewidth,3)% plot(n,p)grid onaxis autoaxis equal
