哈工大机械原理大作业凸轮 - 黄建青
哈工大、机械原理大作业、凸轮机构设计20题
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：能源科学与工程学院班级：1102301设计者：刘平成学号：1110200724指导教师：唐德威设计时间：2013年6月7日凸轮机构设计1．设计题目(1) 凸轮机构运动简图：(2)凸轮机构的原始参数表2-1.凸轮机构原始参数 序号 升程（mm ）升程运动角 升程运动规律升程许用压力角20 110 120° 正弦加速度35°回程运动角回程运动规律 回程许用压力角 远休止角近休止角 90°正弦加速度 65°90°60°（二）凸轮运动方程及相关图像、程序凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图： ○1 凸轮推杆升程、回程方程 πϕπϕϕs)650(πϕ≤≤140)(2=ϕs 511()69πφπ≤≤pi))*5708)/(23.2289)/1.-(sin(2+57083.2289)/1.-(-140(1)(3ϕπϕϕ=s1116()99πφπ≤≤)2914(πϕπ≤≤ 0)(4=ϕs ○2速度方程/2.0944;/2.09440))cos(2-140(1)(1πϕϕ=v 16(2)9πφπ≤≤ 0)(2=ϕv 511()69πφπ≤≤ 708;5708))/1.53.2289)/1.-(cos(2-140(1)(3ϕπϕ=v 1116()99πφπ≤≤ 0)(4=ϕv 16(2)9πφπ≤≤○3加速度方程 .0944^2;/2.0944)/2sin(2280)(1πϕπϕ=a )650(πϕ≤≤0)(2=ϕa 511()69πφπ≤≤08^25708)/1.573.2289)/1.-(sin(2280)(3ϕππϕ=a 1116()99πφπ≤≤ 0)(4=ϕa 16(2)9πφπ≤≤推杆位移、速度、加速度线图matlab编程clear,clcpu=0*pi/180:0.0001:120*pi/180; %升程运动角范围pf=120*pi/180:0.0001:210*pi/180; %远休止角范围pd=210*pi/180:0.0001:300*pi/180; %回程运动角范围pn=300*pi/180:0.0001:2*pi; %近休止角范围h=110e-3; %升程w=10; %凸轮角速度p0=120*pi/180; %升程运动角p01=90*pi/180; %回程运动角ps=90*pi/180; %远休止角%----------推程-----------------------------------------su=h.*(pu./p0-sin(2.*pi.*pu./p0)/(2*pi)); %推杆位移vu=h*w/p0*(1-cos(2*pi*pu./p0)); %推程速度au=2*pi*h*w^2/p0^2*sin(2*pi*pu./p0); %推程加速度%------------远休止----------------------------nf=size(pf);sf=h*ones(nf); %推杆位移vf=zeros(nf); %推程速度af=zeros(nf); %推程加速度%---------------回程------------------------------T=pd-(p0+ps);sd=h/2*(1+cos(pi/p01*T)); %回程位移vd=-pi*h*w/(2*p01)*sin(pi/p01*T); %回程速度ad=-pi^2*h*w^2/(2*p01^2)*cos(pi/p01*T); %回程加速度%--------------------近休止---------------------------------nn=size(pn);sn=zeros(nn); %推杆位移vn=zeros(nn); %推程速度an=zeros(nn); % 推程加速度%------画出推杆位移、速度、加速度线图---------------p=[pu,pf,pd,pn];s=[su,sf,sd,sn];subplot(2,3,1),hold onplot(p,s*1e3,'linewidth',2),xlabel('\phi/rad'),ylabel('s/mm'),grid on,title('推杆位移'),axis([0,2*pi,1.1*min(s)*1e3,1.1*max(s)*1e3]) subplot(2,3,2) v=[vu,vf,vd,vn];plot(p,v,'linewidth',2),xlabel('\phi/rad'),ylabel('v/m/s'),grid on,title('推杆速度'),axis([0,2*pi,1.1*min(v),1.1*max(v)]) subplot(2,3,3) a=[au,af,ad,an];plot(p,a,'linewidth',2),xlabel('\phi/rad'),ylabel('a/m/s^2'),grid on,title('推杆加速度'),axis([0,2*pi,1.1*min(a),1.1*max(a)]) hold off（三）凸轮机构s d ds-ϕ图像及程序代码 %--------------------求ds/d_phi-------------------subplot(2,3,4),plot(v/w*1e3,s*1e3,'linewidth',2),xlabel('ds/d\phi/mm'),ylabel('s/mm'),axis equal,grid on,title('ds/d\phi —s')%---------------------凸轮轴心许用区域--------------------------- alpha_up=35*pi/180; %升程许用压力角 alpha_down=65*pi/180; %回程许用压力角 p1=pi/2-alpha_up; %推程斜率角 p2=alpha_down-pi/2; %回程斜率角 ku=tan(p1); %推程切线斜率 kd=tan(p2); %回程切线斜率 R2=[cos(-p2),-sin(-p2);sin(-p2),cos(-p2)];%推程旋转矩阵 R1=[cos(-p1),-sin(-p1);sin(-p1),cos(-p1)];%推程旋转矩阵 nu=size(pu); for i=1:nu(2)Temp=R1*[vu(i)/w;su(i)];vut(i)=Temp(1); %旋转推程ds/dp-s 曲线 sut(i)=Temp(2); endnd=size(pd); for i=1:nd(2)Temp=R2*[vd(i)/w;sd(i)];vdt(i)=Temp(1); %旋转回程ds/dp-s 曲线 sdt(i)=Temp(2); endfor j=1:nu(2)if sut(j)==min(sut)temu=j; %旋转推程ds/dp-s 曲线后求最低点 end endfor j=1:nd(2)if sdt(j)==min(sdt)temd=j; %旋转回程ds/dp-s曲线后求最低点endendt1=1.2*min(vd/w):0.01:1.2*max(vu/w); %切线定义域t2=min(vd/w)/6:0.01:1.2*max(vu/w);t3=0:0.01:1.2*max(vu/w);s1=ku*(t2-vu(temu)/w)+su(temu); %推程切线s2=kd*(t1-vd(temd)/w)+sd(temd); %回程切线s3=tan(-p1)*t3; %推程起点压力角限制线subplot(2,3,5) %画图hold on,axis equal,grid onplot(v/w*1e3,s*1e3,'linewidth',2)plot(t2*1e3,s1*1e3,'linewidth',1,'color','r')plot(t1*1e3,s2*1e3,'linewidth',1,'color','r')plot(t3*1e3,s3*1e3,'linewidth',1,'color','r')xlabel('ds/d\phi/mm'),ylabel('s/mm'),hold off,title('ds/d\phi—s，轴向许用范围')（四）确定凸轮的基圆半径和偏距、绘制凸轮机圆、偏距圆、理论轮廓曲线---------------画理论廓线图-------------------------e=36e-3;s0=52e-3;r0=sqrt(s0^2+e^2);x=(s0+s).*cos(p)-e.*sin(p);y=(s0+s).*sin(p)+e.*cos(p);x1=r0*cos(p);y1=r0*sin(p);subplot(2,3,6)plot(x*1e3,y*1e3,'linewidth',1),axis equal,grid on,hold on,title('廓线图')plot(x1*1e3,y1*1e3,'linewidth',1,'color','r')%-------------求最小曲率半径-----------------------nx=size(x);nx1=nx(2)-2;dydp=diff(y)./diff(p);%求微分dxdp=diff(x)./diff(p);d2ydp2=diff(dydp)./diff(p(1:nx1+1));d2xdp2=diff(dxdp)./diff(p(1:nx1+1));rho=(dxdp(1:nx1).^2+dydp(1:nx1).^2).^1.5./abs((dxdp(1:nx1).*d2ydp2(1:nx1)-dydp(1:nx1).*d2x dp2(1:nx1)));%理论廓线曲率半径rhomin=min(rho);%最小曲率半径rr=rhomin-3e-3;%----------------实际廓线图----------------X=x(1:nx(2)-1)-rr*dydp./(dxdp.^2+dydp.^2).^0.5;%求实际廓线坐标Y=y(1:nx(2)-1)+rr*dxdp./(dxdp.^2+dydp.^2).^0.5;plot(X*1e3,Y*1e3,'linewidth',2,'color','k')%画实际廓线图Legend('理论廓线','基圆','实际廓线'),axis([1.1*min(x)*1e3,1.1*max(x)*1e3,1.1*min(y)*1e3,1.1*max(y)*1e3])得到基圆半径311mm、偏距36mm。
哈工大机械原理大作业凸轮机构设计第题
哈工大机械原理大作业-凸轮机构设计(第题)————————————————————————————————作者：————————————————————————————————日期：机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电学院班级：1208103完成者：xxxxxxx学号：11208103xx指导教师：林琳设计时间：2014.5.2哈尔滨工业大学凸轮机构设计一、设计题目如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
序号 升程（mm ） 升程运动角（°） 升程运动规律 升程许用压力角（°） 回程运动角（°） 回程运动规律 回程许用压力角（°）远休止角（°） 近休止角（°）3 50 150 正弦加速度 30 100 余弦加速度60 30 80二、凸轮推杆升程、回程运动方程及其线图1 、凸轮推杆升程运动方程(650πϕ≤≤) 升程采用正弦加速度运动规律，故将已知条件mm h 50=，650π=Φ带入正弦加速度运动规律的升程段方程式中得：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512sin 215650ϕππϕS ； ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512cos 1601ππωv ； ω⎪⎭⎫ ⎝⎛=512sin 14421ϕπωa ； 2、凸轮推杆推程远休止角运动方程（πϕπ≤≤65） mm h s 50==；0==a v ；3、凸轮推杆回程运动方程（914πϕπ≤≤） 回程采用余弦加速度运动规律，故将已知条件mm h 50=，95'0π=Φ，6s π=Φ带入余弦加速度运动规律的回程段方程式中得：⎥⎦⎤⎢⎣⎡-+=)(59cos 125πϕs ； ()πϕω--=59sin451v ； ()πϕω-=59cos 81-a 21；4、凸轮推杆回程近休止角运动方程（πϕπ2914≤≤） 0===a v s ；5、凸轮推杆位移、速度、加速度线图根据以上所列的运动方程，利用matlab 绘制出位移、速度、加速度线图。
哈工大机械原理大作业凸轮结构设计3
仅供个人参考Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮结构设计院系：机电工程学院班级：设计者：学号：指导教师：设计时间：哈尔滨工业大学1、设计题目2、凸轮机构推杆升程、回程运动方程，推杆位移、速度、加速度线图。
（1）推杆各行程运动方程（设定角速度为s s rad /2/1πω==）①从动件推程运动方程(650πϕ≤≤) ⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫ ⎝⎛-=ϕφπφω0012cos 1h v ; 代入数据，可得：②从动件远休程运动方程(πϕπ≤≤65) ③从动件回程运动方程(914πϕπ≤≤) 代入数据，可得:④从动件进休程运动方程(πϕπ2914≤≤) （2）、推杆位移、速度、加速度线图①推杆的位移线图如下②推杆的速度线图如下③推杆的加速度线图如下 3、凸轮机构的dss d ϕ-线图，并依次确定凸轮的基圆半径和偏距凸轮机构的s d ds -ϕ线图如下图所示 由图中范围选定点（-10，-50）为凸轮转轴O 点，则mm r 99.501050220=+=取基圆半径为r 0 =51mm ，偏距e = 10mm 。
4、滚子半径的确定及凸轮理论廓线和实际廓线的绘制由程序计算得凸轮理论轮廓线最小曲率半径mm r 51min = .由滚子半径选择范围∆-<min ρr r ,mm 5~3=∆得到滚子半径mm r r 46≤.又因为凸轮整体尺寸较小，此范围明显过大，故适当减小滚子半径,这里取半径为mm r r 12= .得到图线为：附录1.求位移、速度、加速度的程序（matlab ）function f = tulunh=50;x1=150;t1=30;x2=100;t2=80;w=2*pi;x1=x1*pi/180;x2=x2*pi/180;t1=t1*pi/180;t2=t2*pi/180;%升程x=0:0.001:x1;s = h*(x/x1-sin(2*pi*x/x1)/(2*pi));v = h*w*(1-cos(2*pi*x/x1))/x1;a = 2*pi*h*w*w*sin(2*pi*x/x1)/(x1*x1);subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%远休x = x1:0.001:x1+t1;s = h;v=0;a=0;subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%回程x= x1+t1:0.001:x1+t1+x2;s = h*(1+cos(pi*(x-(x1+t1))/x2))/2;v = -pi*h*w*sin(pi*(x-(x1+t1))/x2)/(2*x2);a = -pi*pi*h*w*w*cos(pi*(x-(x1+t1))/x2)/(2*x2*x2);subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%近休x=x1+t1+x2:0.001:x1+x2+t1+t2;s = 0;v = 0;a = 0;subplot(3,1,1),plot(x,s),xlabel('φ/rad'),ylabel('S/mm'),title('位移-转角图线'),hold onsubplot(3,1,2),plot(x,v),xlabel('φ/rad'),ylabel('v/(mm/s)'),title('速度-转角图线'),hold onsubplot(3,1,3),plot(x,a),xlabel(φ/rad'),ylabel('a/(mm/s^2)'),title('加速度-转角图线'),hold on2.绘制凸轮机构s d ds -ϕ线图 function f= jiyuan;x1=150;t1=30;x2=100;t2=80;h=50;x1=x1*pi/180;x2=x2*pi/180;t1=t1*pi/180;t2=t2*pi/180;x= 0:0.001:150*pi/180;%升程 v/ws = h*(x/x1-sin(2*pi*x/x1)/(2*pi));k =-h*(1-cos(2*pi*x/x1))/x1;plot(k,s,'r'),hold on ;x=180*pi/180:0.001:280*pi/180;%回程 v/ws = h*(1+cos(pi*(x-(x1+t1))/x2))/2;k = pi*h*sin(pi*(x-(x1+t1))/x2)/(2*x2);plot(k,s,'g'),hold on ;%回程切线for i=-11:1:-11;f=@(k)k*tan(pi/6)+i;k =-40:0.1:50;s=f(k);plot(k,s),hold on ;end%升程切线for i=-45:0.2:-45;f=@(k)-k*tan(60*pi/180)+i;k =-40:0.1:50;s=f(k);plot(k,s),hold on ;endgrid onf=@(k)k*tan(50*pi/180);k=-50:0.1:0;s=f(k);plot(k,s),hold onxlabel('ds/d φ');ylabel('s(φ)');title('类速度-位移图线 ');plot(-10,-50,’o ’);3.绘制凸轮轮廓曲线function f= lunkuo;h=50;x1=150;t1=30;x2=100;t2=80;x1=x1.*pi./180;x2=x2.*pi./180;t1=t1.*pi./180;t2=t2.*pi./180;s0=51;e=10;rr=12;%升程x=0:pi/200:150.*pi/180;s = h.*(x./x1-sin(2.*pi.*x./x1)./(2.*pi));X1=(s0+s).*cos(x)-e.*sin(x);Y1=(s0+s).*sin(x)+e.*cos(x);%实际轮廓X11=X1-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y11=Y1-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X1,Y1,'r',X11,Y11,'r'),hold on;%远休x=150.*pi/180:pi/180:180.*pi/180;s=50;X2=(s0+s).*cos(x)-e.*sin(x);Y2=(s0+s).*sin(x)+e.*cos(x);X22=X2-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y22=Y2-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X2,Y2,'g',X22,Y22,'g'),hold on;%回程x=180.*pi/180:pi/180:280.*pi/180;s = h.*(1+cos(pi.*(x-(x1+t1))./x2))./2;X3=(s0+s).*cos(x)-e.*sin(x);Y3=(s0+s).*sin(x)+e.*cos(x);X33=X3-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y33=Y3-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X3,Y3,'k',X33,Y33,'k'),hold on;%近休x=280*pi/180:pi/180:2*pi;s=0;X4=(s0+s).*cos(x)-e.*sin(x);Y4=(s0+s).*sin(x)+e.*cos(x);X44=X4-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y44=Y4-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X4,Y4,'b',X44,Y44,'b'),hold on;x=0:pi/200:2*pi;X4=(s0+s).*cos(x)-e.*sin(x); Y4=(s0+s).*sin(x)+e.*cos(x); plot(X4,Y4,'b');text(-40,90,'理论轮廓线');text(-40,72,'实际轮廓线');text(-5,55,'基圆');grid on;axis equal仅供个人用于学习、研究；不得用于商业用途。
哈工大机械原理大作业——凸轮——5号
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院班级：设计者：学号：指导教师：设计时间：凸轮机构设计一． 设计题目表1 凸轮机构原始参数二． 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图1. 推杆升程运动方程（ϕ为转角）（203ϕπ≤≤） 令1/rad s ω= 位移：001212040[sin()]sin(3)2s h ϕπϕϕϕφπφππ=-=- 速度：1002120120[1cos()]cos(3)h v ωπϕϕφφππ=-=- 加速度：2120022360[sin()]sin(3)h a πωπϕϕφφπ==2. 推杆回程运动方程（ϕ为转角）7563πϕπ≤≤ 位移：''00128001604014[1sin()]sin(4)]233T s h T πϕϕπφπφππ=-+=-+- 速度：1''00216016014[1cos()]cos(4)]3h v T ωπϕϕπφφππ=--=-+- 加速度：21'2'002264014sin()sin(4)]3h a T πωπϕπφφπ=-=-- 式中：027()()326s T πππϕϕϕϕϕ=-+=-+=-由MATLAB编程得到线位移图像：线速度图像：线加速度图像：三．绘制机构的ds sdϕ-线图由dsds vdtdd wtϕϕ==可知线图即vsw-线图，由MATLAB编程后得图四．确定基圆半径和偏距1.画切线图如上页图，得其切线图2.求基圆半径和偏距从图中读取e=50mmr0=111.80mm五．画理论及实际轮廓附：MATLAB编程1.推杆位移线图代码%升程阶段t=0:0.001:2*pi/3;s=120.*t./pi-40./pi.*sin(3.*t);hold onplot(t,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;s=80;hold onplot(t,s);%回程阶段t=7*pi/6:0.001:5*pi/3;s=800/3-160.*t./pi+40./pi.*sin(4.*t-14*pi/3) hold onplot(t,s);%近休止阶段t=5*pi/3:0.001:2*pi;s=0;hold onplot(t,s);grid onhold off2.推杆速度线图代码%升程阶段t=0:0.001:2*pi/3;v=120/pi-120/pi*cos(3.*t);hold onplot(t,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;v=0;hold onplot(t,s);%回程阶段t=7*pi/6:0.001:5*pi/3;v=-160/pi+160/pi*cos(4.*t-14/3*pi) hold onplot(t,s);%近休止阶段t=5*pi/3:0.001:2*pi;v=0;hold onplot(t,s);grid onhold off3.推杆加速度线图代码%升程阶段t=0:0.001:2*pi/3;a=360/pi.*sin(3.*t)hold onplot(t,a);%远休止阶段t=2*pi/3:0.001:7*pi/6;a=0;hold onplot(t,a);%回程阶段t=7*pi/6:0.001:5*pi/3;a=-640/pi.*sin(4.*t-14/3*pi) hold onplot(t,a);%近休止阶段t=5*pi/3:0.001:2*pi;a=0;hold onplot(t,a);grid onhold off4.dssdϕ-线图代码%升程阶段t=0:0.001:2*pi/3;s=120.*t./pi-40./pi.*sin(3.*t);v=120/pi-120/pi*cos(3.*t);hold onplot(v,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;s=80;v=0;hold onplot(v,s);%回程阶段t=7*pi/6:0.001:5*pi/3;s=800/3-160.*t./pi+40./pi.*sin(4.*t-14*pi/3);v=-160/pi+160/pi*cos(4.*t-14/3*pi)hold onplot(v,s);5.最终轮廓线图代码h=80;w=1;e=50;rr=20;s0=100;m=120*pi/180;ms=(120+90)*pi/180;m1=(120+90+90)*pi/180; for i=1:1:120mm(i)=i*pi/180.0;s1=h.*( mm(i)./m -(1/(2*pi)).*sin(2.*pi.*mm(i)./m));v1=(w*h./m)*(1-cos(2*pi*mm(i)/m));x(i)=(s0+s1)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s1)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s1)*cos(mm(i))-e*sin(mm(i))+v1/w*sin(mm(i)); b(i)=-(s0+s1)*sin(mm(i))-e*cos(mm(i))+v1/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=121:1:210mm(i)=i*pi/180;s2=h;v2=0;x(i)=(s0+s2)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s2)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s2)*cos(mm(i))-e*sin(mm(i))+v2/w*sin(mm(i));b(i)=-(s0+s2)*sin(mm(i))-e*cos(mm(i))+v2/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=211:1:300mm(i)=i*pi/180;s3=800/3-160/pi.*mm(i)+40/pi.*sin(4.*mm(i)-14/3*pi);v3=-160/pi+160/pi.*cos(4*mm(i)-14/3*pi)x(i)=(s0+s3)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s3)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s3)*cos(mm(i))-e*sin(mm(i))+v3/w*sin(mm(i));b(i)=-(s0+s3)*sin(mm(i))-e*cos(mm(i))+v3/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=301:1:360mm(i)=i*pi/180;x(i)=(s0+0)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+0)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+0)*cos(mm(i))-e*sin(mm(i))+v3/w*sin(mm(i));b(i)=-(s0+0)*sin(mm(i))-e*cos(mm(i))+v3/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endplot(x,y,'r',xx,yy,'g')text(0,20,'实际轮廓线')text(120,100,'理论轮廓线')hold on。
哈工大机械原理大作业-凸轮机构设计
哈工大机械原理大作业-凸轮机构设计(第3题)(共15页)-本页仅作为预览文档封面，使用时请删除本页-机械原理大作业二课程名称：机械原理设计题目：凸轮设计院系：机电学院班级： 1208103完成者： xxxxxxx学号： xx指导教师：林琳设计时间：工业大学凸轮设计一、设计题目如图所示直动从动件盘形凸轮，其原始参数见表，据此设计该凸轮。
二、凸轮推杆升程、回程运动方程及其线图1 、凸轮推杆升程运动方程(650πϕ≤≤) 升程采用正弦加速度运动规律，故将已知条件mm h 50=，650π=Φ带入正弦加速度运动规律的升程段方程式中得：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512sin 215650ϕππϕS ；⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512cos 1601ππωv ； ⎪⎭⎫ ⎝⎛=512sin 14421ϕπωa ； 2、凸轮推杆推程远休止角运动方程（πϕπ≤≤65） mm h s 50==； 0==a v ；3、凸轮推杆回程运动方程（914πϕπ≤≤）回程采用余弦加速度运动规律，故将已知条件mm h 50=，95'0π=Φ，6s π=Φ带入余弦加速度运动规律的回程段方程式中得：⎥⎦⎤⎢⎣⎡-+=)(59cos 125πϕs ；()πϕω--=59sin451v ； ()πϕω-=59cos 81-a 21；4、凸轮推杆回程近休止角运动方程（πϕπ2914≤≤） 0===a v s ；5、凸轮推杆位移、速度、加速度线图根据以上所列的运动方程，利用matlab 绘制出位移、速度、加速度线图。
①位移线图 编程如下： %用t 代替转角 t=0::5*pi/6;s=50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)); hold on plot(t,s); t=5*pi/6::pi; s=50; hold on plot(t,s); t=pi::14*pi/9;s=25*(1+cos(9*(t-pi)/5));hold onplot(t,s);t=14*pi/9::2*pi;s=0;hold onplot(t,s),xlabel('φ/rad'),ylabel('s/mm'); grid onhold off所得图像为：②速度线图编程如下：%用t代替转角,设凸轮转动角速度为1t=0::5*pi/6;v=60/pi*(1-cos((12*t)/5));hold onplot(t,v);t=5*pi/6::pi;v=0;hold onplot(t,v);t=pi::14*pi/9;v=-45*sin(9*(t-pi)/5);hold onplot(t,v);t=14*pi/9::2*pi;v=0;hold onplot(t,v),xlabel('φ(rad)'),ylabel('v(mm/s)'); grid onhold off所得图像为：③加速度线图利用matlab编程如下：%用t代替转角,设凸轮转动角速度为1t=0::5*pi/6;a=144/pi*sin(12*t/5);hold onplot(t,a);t=5*pi/6::pi;a=0;hold onplot(t,a);t=pi::14*pi/9;a=-81*cos(9*(t-pi)/5);hold onplot(t,a);t=14*pi/9::2*pi; a=0; hold onplot(t,a),xlabel('φ(rad)'),ylabel('a(mm/s^2)'); grid on hold off所得图形：三、绘制s d ds -ϕ线图根据运动方程求得：()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧≤≤≤≤-≤≤≤≤--=πϕππϕππϕπϕππϕπππϕ2914.0914,59sin 4565,0650),512cos 6060(d ds 利用matlab 编程：%用t 代替φ,a 代替ds/d φ, t=0::5*pi/6;a=-(60/pi-60/pi*cos(12*t/5));s=50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)); hold on plot(a,s); t=5*pi/6::pi; a=0; s=50; hold on plot(a,s); t=pi::14*pi/9;a=45*sin(9*(t-pi)/5); s=25*(1+cos(9*(t-pi)/5)); hold on plot(a,s);t=14*pi/9::2*pi; a=0; s=0; hold onplot(a,s),title('ds/d φ-s'),xlabel('ds/d φ(mm/rad)'),ylabel('s(mm)'); grid on hold off 得s d ds-ϕ图：凸轮压力角的正切值s s e d ds +-=0/tan ϕα，左侧为升程，作与s 轴夹6π角等于升程许用压力角的切界线t t d D ，则在直线上或其左下方取凸轮轴心时，可使[]αα≤，同理右侧回程，作与s 轴夹角等于回程许用压力角3π的切界线''t t d D ，则在直线上或其右下方取凸轮轴心时，可使[]αα≤。
哈工大机械原理大作业凸轮设计
哈工大机械原理大作业凸轮设计Harb inIn stituteofTech no logy大作业设计说明书课程名称：设计题目：院班学级：机械原理凸轮机构设计1208103系：机械设计制造及其自动化设计指导教师：设计时间：林琳2019425哈尔滨工业大学一、运动分析题目如图所示直动从动件盘形凸轮机构，其原始数据参数见表2-1,。
从表2-1中选择一组凸轮机构原始参数，据此设计该凸轮机构。
二、凸轮运动规律升程运动角（°）90升程运动规律生程许回程运用压力动角角（°）等加等4080减速回程运动规律回程许远休用压力止角角（°）余弦加7040速度近休止角（°）150升程（mm）1501 、升程运动规律（0 /4）位移s=2h（速度v2 /2 4*150*w( /2)A24*150*w A2( /2)人2加速度a2 、升程运动规律(/4 /2)位移s 1502*150( /2 )A2(/2)A2速度v4*60*w( /2 )( /2)A24*60*wA2( /2)A2加速度a3 、回程运动规律(/2 2/2 2 ) 93 位移s 75*{1 cos[ (/2 2 )]}949速度vhw*sin*[ ( /2 2 )]92*4 499加速度aA2hw A2cos[ ( /2 2 )]94 2*(4 )A299根据运动规律做出的曲线以及源代码如图所示位移线图速度线图加速度线图位移线图源代码fl=pi/180;x0=0:fl:pi/4;x1=pi/4:fl:pi/2;x2=pi/2:fl:13*pi/18;x3=13*pi/18:fl:7*pi/6;x4=7*pi/6:fl:2*pi;s0=300*(2*x0/pi)A 2；s1=150-1200*(pi/2-x1).*(pi/2-x1)/(pi.*pi);s2=150+x2*0;s3=75*(1+cos(9/4*(x3-13*pi/18)));s4=x4*0;Plot(x0,s0,x1,s1,'b',x2,s2,'b',x3,s3,'b',x4,s4,'b')axis([070200])title('杆位移线图')xlabel(' 0 (rad)')ylabel('V(mm⑸')gridon速度源代码fl=pi/180;x0=0:fl:pi/4;x1=pi/4:fl:pi/2;x2=pi/2:fl:13*pi/18;x3=13*pi/18: fl:7*pi/6;x4=7*pi/6:fl:2*pi;w=30;v0=600.*w.*x0/(pi/2)A2;v1=600.*w.*(pi/2-x1)/(pi/2)A2;v2=0*x2;v3=-150*30*pi/(2*4*pi/9).*si n(9/4*(x3-13*pi/18));v4=0*x4;Plot(x0,v0,'b',x1,v1,'b',x2,v2,'b',x3,v3,'b',x4,v4,'b')title('推杆速度')xlabel(' 0 (rad)')ylabel('v(mm/s')gridon加速度源代码fl=pi/180;x0=0:fl:pi/4;x1=pi/4:fl:pi/2;x2=pi/2:fl:13*pi/18;x3=13*pi/18: fl:7*pi/6;x4=7*pi/6:fl:2*pi; w=30;a0=600*w.A2/(pi/2).A2+xO*0;a 仁-600*w.A2/(pi/2)A2+x1*0;a2=x2*0; a3=-pi*pi*150*30*30/(2*4*pi/9)A2.*cos(9/4*(x3-13*pi/18));a4=x4*0;Plot(x0,a0,'b',x1,a1,'b',x2,a2,'b',x3,a3,'b',x4,a4,'b')title('推杆加速度')xlabel(' 0 (rad)')ylabel('a(mm/sA2')gridon三、凸轮机构的dss曲线绘制d由凸轮机构位移公式可知4h(/2)A2(0 /4)ds 4*60 ( /2 )( /4 /2)d ( /2)A275*4*s in 9( 13 )( /2 2 /2 2 4 ) 9418999 则其曲线如图所示其源代码如下clcfl=pi/180;x0=0:fl:pi/4;x1=pi/4:fl:pi/2;x2=pi/2:fl:13*pi/18;x3=13*pi/18: fl:7*pi/6;x4=7*pi/6:fl:2*pi;w=30;d0=-600.*x0/(pi/2)A2;d1=-600.*(pi/2-x1)/(pi/2)A2;d2=0*x2;d3=75*4/9.*si n(9/4*(x3-13*pi/18));d4=0*x4;s0=300*(2*x0/pi)A2;s1=150-1200*(pi/2-x1).*(pi/2-x1)/(pi.*pi);s2=150+x2*0;s3=75*(1+cos(9/4*(x3-13*pi/18)));s4=x4*0;Plot(d0,s0,'b',d1,s1,'b',d2,s2,'b',d3,s3,'b',d4,s4,'b')title(' 类速度-位移曲线')xlabel(' 类速度(mm/rad)')ylabel(' 位移(mmm)')gridon四、确定凸轮的基圆半径和偏距以ds/df-s图为基础，可分别作出二条限制线，以这二条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件图像如图所示由图像可知，设置点(50，-100 )为凸轮轴心位置。
哈工大—机械原理凸轮大作业
一、题目要求及机构运动简图如图1所示直动从动件盘形凸轮机构。
其原始参数见表1。
图一凸轮运动简图表一凸轮原始参数二、计算流程框图凸轮机构分析建立数学模型位移方程速度方程加速度方程速度线图位移线图加速线图ds/dΨ-s曲线升程压力角回程压力角确定轴向及基圆半径压力角图确定滚子半径实际轮廓理论轮廓轮廓图结束三、建立数学模型1.从动件运动规律方程首先，由于设计凸轮轮廓与凸轮角速度无关，所以不妨设凸轮运动角速度为w = 1rad/s。
（1）推程运动规律(0 < φ < 90°)s=φ）v=φφ）a=φφ）式中：h=65mm，Φ0=π/2（2）远休程运动规律(90°< φ < 190°)s = 65mmv = 0a = 0（3）回程运动规律(190°< φ < 240°)(190°< φ < 196.25°)(196.25°< φ < 233.75°)(233.75°< φ < 240°) 回程运动中的速度和加速度为位移对时间t的倒数：（4）近休程运动规律(240°< φ < 360°)s = 0v = 0a = 02.从动件位移、速度、加速度线图（1）位移线图（2）速度线图（3）加速度线图（4）位移、速度、加速度线图MATLAB源程序%% 已知条件h = 65; %mmphi_0 = 90./180*pi; %radalpha_up_al = 35./180*pi; %升程许用压力角phi_00 = 50./180*pi;alpha_down_al = 70./180*pi; %回程许用压力角phi_s = 100./180*pi;phi_ss = 120./180*pi;w = 1;%% 绘制从动件位移、速度、加速度线图% 推程阶段t_up = 0 : 0.5 : 90;t_up1 = t_up./180*pi;syms t_up1phi_up s_up v_up a_upphi_up = w.*t_up1;s_up = h./2.*(1 - cos(pi.*phi_up./phi_0));v_up = diff(s_up,t_up1);a_up = diff(v_up,t_up1);s_up1 = double(subs(s_up,t_up./180*pi));v_up1 = double(subs(v_up,t_up./180*pi));a_up1 = double(subs(a_up,t_up./180*pi));% 远休程t_s = 90 : 0.5 : (90+100);t_s1 = t_up./180*pi;s_s(1:201) = h;v_s(1:201) = 0;a_s(1:201) = 0;% 回程阶段1t_down1 = (90+100) : 0.5 : (90+100+50/8);t_down11 = t_down1./180*pi;syms t_down11phi_down1s_down1v_down1a_down1phi_down1 = w.*t_down11;s_down1 = h - h./(4+pi).*(pi.*(phi_down1 - phi_0 - phi_s)./phi_00 - ...sin(4.*pi.*(phi_down1 - phi_0 - phi_s)./phi_00)./4);v_down1 = diff(s_down1,t_down11);a_down1 = diff(v_down1,t_down11);s_down11 = double(subs(s_down1,t_down1./180*pi));v_down11 = double(subs(v_down1,t_down1./180*pi));a_down11 = double(subs(a_down1,t_down1./180*pi));% 回程阶段2t_down2 = (90+100+50/8) : 0.5 : (90+100+7*50/8);t_down22 = t_down2./180*pi;syms t_down22phi_down2s_down2v_down2a_down2phi_down2 = w.*t_down22;s_down2 = h - h./(4+pi).*(2+pi.*(phi_down2 - phi_0 - phi_s)./phi_00 - 9.*sin(pi./3 + 4.*pi.*(phi_down2 - phi_0 - phi_s)./(3.*phi_00))./4);v_down2 = diff(s_down2,t_down22);a_down2 = diff(v_down2,t_down22);s_down22 = double(subs(s_down2,t_down2./180*pi));v_down22 = double(subs(v_down2,t_down2./180*pi));a_down22 = double(subs(a_down2,t_down2./180*pi));% 回程阶段3t_down3 = (90+100+7*50/8) : 0.5 : (90+100+50);t_down33 = t_down3./180*pi;syms t_down33phi_down3s_down3v_down3a_down3phi_down3 = w.*t_down33;s_down3 = h - h./(4+pi).*(4+pi.*(phi_down3 - phi_0 - phi_s)./phi_00 - …sin(4.*pi.*(phi_down3 - phi_0 - phi_s)./phi_00)./4);v_down3 = diff(s_down3,t_down33);a_down3 = diff(v_down3,t_down33);s_down33 = double(subs(s_down3,t_down3./180*pi));v_down33 = double(subs(v_down3,t_down3./180*pi));a_down33 = double(subs(a_down3,t_down3./180*pi));% 近休程t_ss = (90+100+50) : 0.5 : 360;s_ss(1:241) = 0;v_ss(1:241) = 0;a_ss(1:241) = 0;% 绘图位移t = [t_up t_s t_down1 t_down2 t_down3 t_ss];phi = w .* t ./ 180 .*pi;s = [s_up1 s_s s_down11 s_down22 s_down33 s_ss];v = [v_up1 v_s v_down11 v_down22 v_down33 v_ss];a = [a_up1 a_s a_down11 a_down22 a_down33 a_ss];figure('Name','从动件位移-时间线图');plot(t,s,'k','linewidth',1.0);grid on;title('从动件位移-时间线图');xlabel('转角\phi / 度');ylabel('位移h/mm');% 绘图速度figure('Name','从动件速度-时间线图');plot(t,v,'k','linewidth',1.0);grid on;title('从动件速度-时间线图');xlabel('转角\phi / 度');ylabel('速度v/mm*s^{-1}');% 绘图加速度figure('Name','从动件加速度-时间线图');plot(t,a,'k','linewidth',1.0);grid on;title('从动件加速度-时间线图');xlabel('转角\phi / 度');ylabel('加速度a/mm*s^{-2}');3.绘制ds/dΦ线图并确定基圆半径和偏距（1）绘制ds/dΦ线图及源程序①MATLAB源程序：%% 绘制ds/dphi-s线图，确定基圆半径和偏距ds_dphi = v ./ w;figure('Name','凸轮ds/dphi - s线图');plot(ds_dphi,s,'k','linewidth',1.5);hold on;axis([-150 150 -70 70]);grid on;title('凸轮ds/dphi - s线图');xlabel('ds/dphi / (mm*s^{-2})');ylabel('s/mm');% 三条临界线x = linspace(-150,150,301);k_up = tan(pi/2 - alpha_up_al);y_up = k_up.*x - 66;plot(x,y_up,'linewidth',1.5);k_down = - tan(pi/2 - alpha_down_al);y_down = k_down.*x - 24.7;plot(x,y_down,'linewidth',1.5);x0 = linspace(0,150,151);k0 = - tan(alpha_up_al);y0 = k0.*x0;plot(x0,y0,'--');% 由图像选取凸轮基圆半径为r0 = sqrt(23^2 + 34^2) = 41 mm,偏距e = 23mm plot(23,-34,'or');r0 = 41;e = 23;plot(linspace(0,23,10),linspace(0,-34,10),'r',linspace(0,23,10),linspace(-34,-34,10),'r',li nspace(23,23,10),linspace(0,-34,10),'r','linewidth',1.0);（2）确定基圆半径和偏距在凸轮机构的ds/d φ-s 线图里再作斜直线Dt-dt 与升程的[ds/d φ-s]曲线相切并使与纵坐标夹角为升程许用压力角[α]，则Dt-dt 线的右下方为选择凸轮轴心的许用区。
哈工大机械原理课程设计—产品生产包装线设计(方案)-黄建青
H a r b i n I n s t i t u t e o f T e c h n o l o g y课程设计说明书（论文）课程名称：机械原理课程设计设计题目：产品包装生产线(方案7)院系：能源科学与工程学院班级： 1302402设计者：黄建青学号： 1130240222指导教师：陈明设计时间： 2015.07.6-2015.07.11哈尔滨工业大学产品包装生产线(方案7)1.设计课题概述如图1所示，输送线1上为小包装产品，其尺寸为长*宽*高=600*200*200，采取步进式输送方式，送第一包产品至托盘A 上（托盘A 上平面与输送线1的上平面同高）后，托盘A 下降200mm ，第二包产品送到后，托盘A 上升205mm 、顺时针旋转90°，把产品推入输送线2，托盘A 逆时针回转90°、下降5mm 。
原动机转速为1430rpm ，产品输送量分三档可调，每分钟向输送线2分别输送10、16、22件小包装产品。
图1功能简图2.设计课题工艺分析（1）工艺方法分析推动产品在输送线1上运动的是执行机构1，在A 处使产品上升、转位的是执行构件2，在A 处把产品推到下一个工位的是执行构件3，三个执行构件的运动协调关系如图所示。
执行构件1、2、3的工作周期分别为T 1、T 2、T 3，构件3的动作周期是T 3’。
三个执行构件的工作周期关系为：2T 1= T 2= T 3,T 3’=1/20T 3。
图2 运动循环图3.运动功能分析及运动功能系统图根据前面的分析可知，驱动执行构件1工作的执行机构应该具有运动功能如图3所示。
该运动功能把一个连续的单向转动转换为连续的往复移动，主动件每转动一周，从动件（执行构件1）往复运动一次，主动件的转速分别为10、16、22 rpm 。
10、16、22 rpm图3 执行机构1的运动功能 由于电动机转速为1430rpm ，为了在执行机构1的主动件上分别得到10、16、22 rpm 的转速，则由电动机到执行机构1之间的传动比i z 有3种,分别为：123143014310143089.3751614306522z z z i i i ====== 总传动比由定传动比i c 与变传动比i v 组成，满足以下关系式：112233z c v z c v z c v i i i i i i i i i =⋅=⋅=⋅三种传动比中i z1最大，i z3最小。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
H a r b i n I n s t i t u t e o f T e c h n o l o g y机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：能源学院班级： 1302402 设计者：黄建青学号： 1130240222 指导教师：焦映厚陈照波设计时间： 2015年06月23日凸轮机构设计说明书1. 设计题目设计直动从动件盘形凸轮机构，机构运动简图如图1，机构的原始参数如表1所示。
图1 机构运动简图表1 凸轮机构原始参数计算流程框图：2. 凸轮推杆升程，回程运动方程及推杆位移、速度、加速度线图2.1 确定凸轮机构推杆升程、回程运动方程 设定角速度为ω=1 rad/s(1) 升程：0°<φ<50°由公式可得 )]cos(1[20ϕπΦh s -=)sin(201ϕπωπΦΦh v =)cos(202212ϕπωπΦΦh a =(2) 远休止：50°<φ<150° 由公式可得s = 45 v = 0a = 0(3) 回程：150°<φ<240° 由公式得：()()2200002000000022000000,2(1)(1)1,12(1)(1),2(1)s s s s s s s s s Φhn s h ΦΦΦΦΦΦn ΦnΦΦn h n s h ΦΦΦΦΦΦn Φn n ΦΦΦn hn s ΦΦΦΦΦn Φn ϕϕϕϕϕϕ'⎧=---+<≤++⎪'-⎪⎪⎡⎤''-⎪=----++<≤++⎨⎢⎥'-⎣⎦⎪⎪'---⎪'=-++<≤++'-⎪⎩2010000001000200100000n (),(1)(1)n ,(1)(1)n (1),(1)s s s s s s s s Φh v ΦΦΦΦΦΦn Φn ΦΦn h v ΦΦΦΦn Φn n ΦΦΦn h v ΦΦΦΦΦn ΦΦn ωϕϕωϕϕωϕ⎧'=---+<≤++⎪'-⎪⎪''-⎪=-++<≤++⎨'-⎪⎪'---'⎪=--++<≤++''-⎪⎩220100200000220100020n ,(1)(1)0,(1)n ,(1)s s s s s s Φh a ΦΦΦΦn ΦnΦΦn a ΦΦΦΦn n Φn h a ΦΦΦΦΦn Φn ωϕϕωϕ'⎧=-+<≤++⎪'-⎪⎪''-=++<≤++⎨⎪⎪'-'=-++<≤++⎪'-⎩式中 Φ0——推程运动角； Φs ——远休止角；Φ0‘——回程运动角。
(4) 近休止：240°<φ<360° 由公式可得： s = 0 v = 0a = 02.2 位移、速度、加速度曲线源代码n=3;p0=50*pi/180; p1=150*pi/180;p2=150*pi/180+pi/(2*n); p3=150*pi/180+pi*(n-1)/(2*n); p4=240*pi/180; x1=0:(pi/300):p0;s1=(45/2)*(1-cos(pi*x1/p0)); v1=pi*45/(2*p0)*sin(pi*x1/p0); a1=(pi^2)*45/(2*p0^2)*cos(pi*x1/p0); x2=p0:(pi/300):p1; s2=45; v2=0; a2=0;x3=p1:(pi/300):p2;s3=45-9*45*(x3-5*pi/6).^2/(pi^2);v3=-45*(n^2)*(x3-p1)/((n-1)*((pi/2)^2)); a3=-45*(n^2)/((n-1)*(pi/2)^2); x4=p2:(pi/300):p3;s4=45-45*(n*(x4-p1)/(pi/2)-0.5)/(n-1); v4=-45*n/((n-1)*(pi/2)); a4=0;x5=p3:(pi/300):p4;s5=45*n^2*(1-(x5-p1)/(pi/2)).^2/(2*(n-1)); v5=-45*n^2*(1-(x5-p1)/(pi/2))/((pi/2)*(n-1)); a5=45*n^2/((n-1)*(pi/2)^2);x6=p4:(pi/300):(2*pi);s6=0;v6=0;a6=0;figure(1);plot(x1,s1,'b',x2,s2,'b',x3,s3,'b',x4,s4,'b',x5,s5,'b',x6,s6,'b'); title('推杆线位移图 ');xlabel('φ(rad)');ylabel('S(mm)');figure(2);plot(x1,v1,'b',x2,v2,'b',x3,v3,'b',x4,v4,'b',x5,v5,'b',x6,v6,'b'); title('推杆速度线图');xlabel('φ(rad)');ylabel('V(mm/s)');figure(3);plot(x1,a1,'b',x2,a2,'b',x3,a3,'b',x4,a4,'b',x5,a5,'b',x6,a6,'b'); title('推杆加速度线图');xlabel('φ(rad)');ylabel('a(mm/s^2)');图1 推杆位移曲线图2 推杆速度曲线图3 推杆加速度曲线3. 运动线图及凸轮sdds-φ线图n=3;p0=50*pi/180;p1=150*pi/180;p2=150*pi/180+pi/(2*n);p3=150*pi/180+pi*(n-1)/(2*n); p4=240*pi/180;x1=0:(pi/300):p0;s1=(45/2)*(1-cos(pi*x1/p0));f1=45*pi*sin(pi*x1/p0)/(2*p0); x2=p0:(pi/300):p1;s2=45;f2=0;x3=p1:(pi/300):p2;s3=45-9*45*(x3-5*pi/6).^2/(pi^2); f3=-2*9*45*(x3-5*pi/6)/(pi^2); x4=p2:(pi/300):p3;s4=45-45*(n*(x4-p1)/(pi/2)-0.5)/(n-1); f4=-90*n/(pi*(n-1)); x5=p3:(pi/300):p4;s5=45*n^2*(1-(x5-p1)/(pi/2)).^2/(2*(n-1)); f5=-90*n^2*(1-(x5-p1)/(pi/2))/(pi*(n-1)); x6=p4:(pi/300):(2*pi); s6=0; f6=0;plot(f1,s1,f2,s2,f3,s3,f4,s4,f5,s5,f6,s6); title('凸轮机构ds/d ψ-s 线图'); gird on ;图4s d ds-φ线图 4.凸轮机构基圆半径和偏距的确定以图4为基础，可分别作出三条限制线（推程许用压力角的切界限D t d t ，回程许用压力角的限制线D t 'd t '，起始点压力角许用线B 0d '')，以这三条线可确定最小基圆半径及所对应的偏距e ，在其下方选择一合适点，即可满足压力角的限制条件。
利用MATLAB 作图，其代码如下，得出图如图5所示。
n=3;p0=50*pi/180; p1=150*pi/180;p2=150*pi/180+pi/(2*n); p3=150*pi/180+pi*(n-1)/(2*n); p4=240*pi/180;x1=0:(pi/100):p0;s1=(45/2)*(1-cos(pi*x1/p0));v1=pi*45/(2*p0)*sin(pi*x1/p0);a1=(pi^2)*45/(2*p0^2)*cos(pi*x1/p0);x2=p0:(pi/100):p1;s2=45;v2=0;a2=0;x3=p1:(pi/100):p2;s3=45-9*45*(x3-5*pi/6).^2/(pi^2);v3=-45*(n^2)*(x3-p1)/((n-1)*((pi/2)^2));a3=-45*(n^2)/((n-1)*(pi/2)^2);x4=p2:(pi/100):p3;s4=45-45*(n*(x4-p1)/(pi/2)-0.5)/(n-1);v4=-45*n/((n-1)*(pi/2));a4=0;x5=p3:(pi/100):p4;s5=45*n^2*(1-(x5-p1)/(pi/2)).^2/(2*(n-1));v5=-45*n^2*(1-(x5-p1)/(pi/2))/((pi/2)*(n-1)); a5=45*n^2/((n-1)*(pi/2)^2);x6=p4:(pi/100):(2*pi);s6=0;v6=0;a6=0;N=-100:(pi/100):85;k1=tan(pi/2-35*pi/180);K1=v1./a1;t=2;while abs(K1(t)-k1)>0.5t=t+1;endX1=v1(t);Y1=s1(t);m1=k1*(N-X1)+Y1;k5=-tan(pi/2-70*pi/180);K5=v5./a5;t=2;while abs(K5(t)-k5)>0.1t=t+1;endX5=v5(t);Y5=s5(t);m5=k5*(N-X5)+Y5;m2=tan(pi/2+35*pi/180)*N;plot(v1,s1,v2,s2,v3,s3,v4,s4,v5,s5,v6,s6,N,m5,N,m2,N,m1); title('凸轮机构ds/d\phi-s 线图'); grid on ;图5由图5得：可取x 0=20mm ，y 0=-70mm则e =20mm ，2200072.82mm r x y =+= 5.滚子半径的确定为求滚子许用半径，须确定最小曲率半径，以防止凸轮工作轮廓出现尖点或出现相交包络线，确定最小曲率半径数学模型如下：)/)(/()/)(/(])/()/[(22222/322ϕϕϕϕϕϕρd x d d dy d y d d dx d dy d dx -+=其中 ϕϕϕϕcos )(sin ])/[(/0s s e d ds d dx ++-=ϕϕϕϕsin )(cos ])/[(/0s s e d ds d dy +--=ϕϕϕϕϕsin ])/[(cos ])/(2[/02222s s d s d e d ds d x d --+-= ϕϕϕϕϕcos ])/[(sin ])/(2[/02222s s d s d e d ds d y d --+--=理论廓线数学模型：ϕϕϕϕsin cos )(cos sin )(00e s s y e s s x -+=++=凸轮实际廓线坐标方程式：22'22')/()/()/()/()/()/(ϕϕϕϕϕϕd dy d dx d dy r y y d dy d dx d dx r x x rr+-=++= 其中r r 为确定的滚子半径。