哈工大机械原理大作业凸轮机构32..docxVIP

Harbin Institute of Technology凸轮机构设计说明书课程名称：机械原理设计题目：凸轮机构设计（32）院系：能源科学与工程学院班级：1302403设计者：刘远华学号：1130240333指导教师：焦映厚陈照波设计时间：2015年6月哈尔滨工业大学凸轮机构设计一、设计题目行程/mm升程运动角/。升程运动规律升程许用压力角/。回程运动角/。回程运动规律回程许用压力角/。远休止角/。近休止角/。50603-4-5多项式3080简-直-、凸轮升程、回程运动方程及推杆位移、速度、加速度线图1、运动规律(1)推程规律方程（3-4-5多项式）位移方程S=h[10*(φ/φ0)^3-15*(φ/φ0)^4+6*(φ/φ0)^5]S=50[10*(3φ/π)^3-15*(3φ/π)^4+6*(3φ/π)^5]速度方程V=hω1/φ0[30*(φ/φ0)^2-60*(φ/φ0)^3+30*(φ/φ0)^4]V=150ω1/π[30*(3φ/π)^2-60*(3φ/π)^3+30*(φ/π)^4]加速度方程a=h(ω1)^2/(φ0)^2[60*(φ/φ0) -180*(φ/φ0)^2+120*(φ/φ0)^3]a=450(ω1)^2/π^2[60*(3φ/π) -180*(3φ/π)^2+120*(3φ/π)^3]（2）回程运动规律方程（简-直-简）位移方程S=h-2h/[4+(n-2)π]{1-cos[nπ(φ-φ0-φs速度方程加速度方程2、编程及代码位移、速度、加速度MATLAB代码clear,clc%推程运动规律方程h = 50;omega1 = 10;PHI0 = pi/3;phi1 = 0:pi/1000:PHI0;s1 = h*(10*(phi1/PHI0).^3 - 15*(phi1/PHI0).^4 + 6*(phi1/PHI0).^5);nu1 = ((h*omega1)/PHI0)*(30*(phi1/PHI0).^2 - 60*(phi1/PHI0).^3 + 30*(phi1/PHI0).^4);a1 = ((h*omega1^2)/PHI0^2)*(60*(phi1/PHI0) - 180*(phi1/PHI0).^2 + 120*(phi1/PHI0).^3);%远休止行程PHIs = 5*pi/9;phi2 = [PHI0,PHI0+PHIs];s2 = [s1(end),s1(end)];nu2 = [0,0];a2 = [0,0];%回程运动规律方程（简-直-简）n = 3;PHI00 =4*pi/9;phi3 = (PHI0 + PHIs):pi/1000:(PHI0 + PHIs + PHI00/n); phi4 = (PHI0 + PHIs + PHI00/n):pi/1000:(PHI0 + PHIs + PHI00*(n - 1)/n); phi5 = (PHI0 + PHIs + PHI00*(n - 1)/n):pi/1000:(PHI0 +PHIs + PHI00); s3 = h - (2*h/(4 + (n - 2)*pi))*(1 - cos(n*pi*(phi3 - PHI0 - PHIs)/(2*PHI00)));s4 = h - (h/(4 + (n - 2)*pi))*(n*pi*(phi4 - PHI0 - PHIs)/PHI00 - pi + 2);s5 = (2*h/(4 + (n - 2)*pi))*(1 + cos(n*pi*(phi5 - PHI0 - PHIs)/(2*PHI00) - (n - 2)*pi/2));nu3 = -(h*n*pi*omega1/((4 + (n - 2)*pi)*PHI00))*sin(n*pi*(phi3 - PHI0 - PHIs)/(2*PHI00));nu4 = -h*n*pi*omega1/((4 + (n - 2)*pi)*PHI00) + zeros(size(phi4));nu5 = -(h*n*pi*omega1/((4 + (n - 2)*pi)*PHI00))*sin(n*pi*(phi5 - PHI0 - PHIs)/(2*PHI00) - (n - 2)*pi/2); a3 = -(h*n^2*pi^2*omega1^2/(2*(4 + (n - 2)*pi)*PHI00^2))*cos(n*pi*(phi3 - PHI0 - PHIs)/(2*PHI00));a4 = 0;a5 = -(h*n^2*pi^2*omega1^2/(2*(4 + (n - 2)*pi)*PHI00^2))*cos(n*pi*(phi5 - PHI0 - PHIs)/(2*