哈工大机械原理大作业凸轮27..docVIP

机械原理大作业三 课程名称：机械原理 设计题目：直动从动件盘形凸轮机构设计 院 系：机电工程学院 班 级：1308301 完 成 者：于华清 学 号：1130830104 指导教师：赵永强 唐德威 设计时间：2015.6.16 哈尔滨工业大学 设计题目 凸轮机构原始参数 行程（mm) 升程运动角（o） 升程运动规律 升程许用压力角（o） 回程运动角（o） 回程运动规律 回程许用压力角（o） 远休止角（o） 近休止角（o） 20 50 3-4-5多项式 30 90 等加等减速 60 100 120 凸轮推杆运动规律 推程 回程 三、运动线图及凸轮线图 采用Matlab编程，其源代码如下： function [psi,s,v,a,f,F]=tl(h,phi01,phis01,phi02,phis02,n) w=2*pi*n/60; psi1=linspace(0,phi01,5*round(phi01)+1); s1=h/2*(1-cos(pi*psi1/phi01)); f1=h/2*(pi/(phi01*pi/180))*sin(pi*psi1/phi01); F1=h*pi^2/2*sin(pi*psi1/phi01)/((phi01*pi/180)^2); v1=pi*h*w*sin(pi*psi1/phi01)/2/(phi01*pi/180); a1=pi^2*h*w^2*cos(pi*psi1/phi01)/2/((phi01*pi/180)^2); psi2=linspace(phi01+1,phis01+phi01,5*round(phis01)); s2=h*ones(1,length(psi2)); f2=zeros(1,length(psi2));F2=f2; v2=zeros(1,length(psi2)); a2=zeros(1,length(psi2)); psi3=linspace(phi01+phis01,phi01+phis01+phi02,5*round(phi02)); s3=h-10*h*((psi3-phi01-phis01)/phi02).^3+15*h*((psi3-phi01-phis01)/phi02).^4-6*h*((psi3-phi01-phis01)/phi02).^5; f3=-30*h/(phi02*pi/180)*((psi3-phi01-phis01)/phi02).^2+60*h/(phi02*pi/180)*((psi3-phi01-phis01)/phi02).^3-30*h/(phi02*pi/180)*((psi3-phi01-phis01)/phi02).^4; F3=-60*h/(phi02*pi/180)^2*((psi3-phi01-phis01)/phi02)+180*h/(phi02*pi/180).^2*((psi3-phi01-phis01)/phi02).^2-120*h/(phi02*pi/180)^2*((psi3-phi01-phis01)/phi02).^3; v3=f3*w; a3=F3*w^2; psi4=linspace(phi01+phis01+phi02+1,360,5*round(phis02)); s4=zeros(1,length(psi4)); f4=zeros(1,length(psi4));F4=f4; v4=zeros(1,length(psi4)); a4=zeros(1,length(psi4)); psi=[psi1,psi2,psi3,psi4]; s=[s1,s2,s3,s4]; f=[f1,f2,f3,f4];F=[F1,F2,F3,F4]; v=[v1,v2,v3,v4]; a=[a1,a2,a3,a4]; end 绘图程序： function huatu1(psi,s,v,a,f) figure(1) subplot(3,1,1); plot(psi,s),grid on xlabel(凸轮转角(^o) );ylabel(位移（mm）); title(位移（mm）) subplot(3,1,2); plot(psi,v),grid on xlabel(凸轮转角(^o));ylabel(速度（mm/s）); subplot(3,1,3); plot(psi,a),grid on xlabel(凸轮转角(^o