哈工大机械原理大作业——凸轮——16号.docVIP

机械原理大作业 第二题 （2013年春季学期） 题目 凸轮机构设计 题号：25号 班 级： 班 级： 设 计 者： 设 计 者： 学 号： 学 号： 指导教师： 指导教师： 设计时间： 设计时间： 哈尔滨工业大学 升程 （mm） 升程运动角（度） 升程运动规律 升程许用压力角（度） 回程运动角（度） 回程运动规律 回程许用压力角（度） 远休止角（度） 近休止角（度） 130 90 正弦加速度 30 80 3,4,5多项式 70 130 60 1、凸轮推程方程为： 2、回程方程为： 3、以上是推程和回程的运动方程，在远休止角以及近休止角处，位移不变，速度、加速度都为“0”，根据以上方程组在matlab中编程输出结果，绘制曲线如下： 4、用matlab 编写程序算出凸轮推杆的位移、速度、加速度曲线，并绘制出ds/dfi—s线图 程序如下： syms fi; %定义推杆运动分段函数(0fipi/2) s1=130*(2*fi/pi-(sin(4*fi)/2/pi)); v1=(260/pi)*(1-cos(4*fi)); a1=(8*130/pi)*sin(4*fi); d1=diff(s1,fi); k1=0:pi/90:pi/2; s1=subs(s1,fi,k1); v1=subs(v1,fi,k1); a1=subs(a1,fi,k1); d1=subs(d1,fi,k1); % pi/2= fi 11*pi/9 k2=pi/2:pi/90:11*pi/9; s2=130; v2=0; a2=0; %11*pi/9 =fi 5*pi/3 t2=9*(fi-(pi/2)-(13*pi/18))/4/pi; s3=130*(1-(10*t2^3)+(15*t2^4)-(6*t2^5)); v3=-(30*130*9/4/pi)*t2^2*(1-2*t2+t2^2); a3=-(60*130*81/16/pi/pi)*t2*(1-3*t2+2*t2^2); d3=diff(s3,fi); k3=11*pi/9:pi/90:5*pi/3; s3=subs(s3,fi,k3); v3=subs(v3,fi,k3); a3=subs(a3,fi,k3); d3=subs(d3,fi,k3); % 5*pi/3 =fi2*pi k4=5*pi/3:pi/90:2*pi; s4=0; v4=0; a4=0; %画函数曲线（推杆位移） plot(k1,s1); hold on plot(k2,s2); hold on plot(k3,s3); hold on plot(k4,s4); title(推杆位移曲线); xlabel(x角度fi(rad)),ylabel(y位移s(mm)); hold off pause % 推杆速度 plot(k1,v1); hold on plot(k2,v2); hold on plot(k3,v3); hold on plot(k4,v4); title(推杆速度曲线); xlabel(x角度fi(rad)),ylabel(y速度v(mm/s)); hold off pause %推杆加速度 plot(k1,a1); hold on plot(k2,a2); hold on plot(k3,a3); hold on plot(k4,a4); title(推杆加速度曲线); xlabel(x角度fi(rad)),ylabel(y加速度a(mm/s^2)); hold off pause