《哈工大机械原理大作业凸轮机构第四题.doc》由会员分享，可在线阅读，更多相关《哈工大机械原理大作业凸轮机构第四题.doc（18页珍藏版）》请在金锄头文库上搜索。
1、Harbin Institute of Technology机械原理大作业二课程名称： 机械原理 设计题目： 凸轮机构设计 姓 名： 李清蔚 学 号： 1140810304 班 级： 1408103 指导教师： 林 琳 一.设计题目设计直动从动件盘形凸轮机构，其原始参数见表1表一： 凸轮机构原始参数升程（mm)升程运动角（）升程运动规律升程许用压力角（）回程运动角（）回程运动规律回程许用压力角（）远休止角（）近休止角（）4090等加等减速30504-5-6-7多项式60100120二.凸轮推杆运动规律（1） 推程运动规律（等加速等减速运动） 推程 F0=90 位移方程如下： 速度方程如下： 加
2、速度方程如下：（2） 回程运动规律（4-5-6-7多项式） 回程 ，F0=90，Fs=100，F0=50其中回程过程的位移方程，速度方程，加速度方程如下：三.运动线图及凸轮线图本题目采用Matlab编程，写出凸轮每一段的运动方程，运用Matlab模拟将凸轮的运动曲线以及凸轮形状表现出来。代码见报告的结尾。1、程序流程框图开始输入凸轮推程回程的运动方程输入凸轮基圆偏距等基本参数输出压力角、曲率半径图像输出ds,dv,da图像结束输出凸轮的构件形状2、运动规律ds图像如下：速度规律dv图像如下：加速度da规律如下图：3.凸轮的基圆半径和偏距 以ds/df-s图为基础，可分别作出三条限制线（推程许用
3、压力角的切界限Dtdt,回程许用压力角的限制线Dtdt，起始点压力角许用线B0d)，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。得图如下：得最小基圆对应的坐标位置O点坐标大约为（13，-50）经计算取偏距e=13mm，r0=51.67mm.2.绘制理论轮廓线上的压力角曲线和曲率半径曲线针对凸轮转向及推杆偏置，令N1=1凸轮逆时针转；N2=1偏距为正。压力角数学模型：曲率半径数学模型：其中：曲率半径以及压力角的图像如下图：压力角图像：曲率半径图像：3.凸轮理论廓线和实际廓线理论廓线数学模型： 凸轮实际廓线坐标方程式： 其中R0为确定的滚子半径，令
4、R0=15mm利用Matlab程序写出凸轮的理论廓线和实际廓线以及基圆、偏距圆的曲线方程，并且进行模拟 ，可以得出如下凸轮图像： *附页相关程序代码： 位移曲线phi1=linspace(0,90/2/180*pi); %phi2=linspace(90/2/180*pi,90/180*pi);phi0=90/180*pi;h=40;s1=2*h*(phi1/phi0).2;s2=h-2*h*(phi0-phi2).2/(phi0).2;plot(phi1,s1)hold onplot(phi2,s2)hold onphi3=linspace(deg2rad(90),deg2rad(190),
5、1000); %s3=h;plot(phi3,s3,-b)hold onphi4=linspace(deg2rad(190),deg2rad(240); %s4=h*(1-35*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).4+84*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).5-70*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).6+20*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).7);plot(phi4,s4)hold onp
6、hi5=linspace(deg2rad(240),deg2rad(360),1000); %s5=0;plot(phi5,s5,-b)axis(0,0+11/10*deg2rad(360),0,0+11/10*h);grid onylabel(位移/mm);xlabel(凸轮转角 /rad); 速度曲线phi0=deg2rad(90); phi1=linspace(0,deg2rad(90/2);h=40;n=1;w=2*pi*n/60;v1=4*h*w/(deg2rad(90).2*phi1;phi2=linspace(deg2rad(90/2),deg2rad(90);v2=4*h*w*
7、(phi0-phi2)/(deg2rad(90).2;phi3=linspace(deg2rad(90),deg2rad(190),1000); v3=0;phi4=linspace(deg2rad(190),deg2rad(240); v4=-h*w/deg2rad(50)*140*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).3-420*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).4+420*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).5-140*(phi4-deg2
8、rad(90)-deg2rad(100)/deg2rad(50).6;phi5=linspace(deg2rad(240),deg2rad(360),1000); v5=0;plot(phi1,v1)hold onplot(phi2,v2)hold onplot(phi3,v3,-b)hold onplot(phi4,v4)hold onplot(phi5,v5,-b)grid onylabel(速度 /mm);xlabel(凸轮转角/rad)加速度曲线phi0=deg2rad(90); %phi1=linspace(0,deg2rad(90/2),10.3);h=40;n=1;w=2*pi*
9、n/60;a1=4*h*w.2/(phi0).2;phi2=linspace(deg2rad(90/2),deg2rad(90),10.3); a2=-4*h*w.2/(phi0).2;phi3=linspace(deg2rad(90),deg2rad(190),1000); %a3=0;phi4=linspace(deg2rad(190),deg2rad(240); %a4=-(h*w.2/(deg2rad(50).2)*420*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).2-1680*(phi4-deg2rad(90)-deg2rad(100)/d
10、eg2rad(50).3+2100*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).4-840*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).5;phi5=linspace(deg2rad(240),deg2rad(360),1000); %a5=0;phi6=deg2rad(45),deg2rad(45);a6=4*h*w.2/(phi0).2,-4*h*w.2/(phi0).2;plot(phi1,a1,-b)hold onplot(phi2,a2,-b)hold onplot(phi3,a3,-b)hold o
11、nplot(phi4,a4)hold onplot(phi5,a5)hold onplot(phi6,a6)grid onylabel(加速度/mm);xlabel(凸轮转角/rad);3、 曲线 %ds/f-sphi1=linspace(0,90/2/180*pi,1000); phi2=linspace(90/2/180*pi,90/180*pi,1000);phi0=90/180*pi;h=40;s01=2*h*(phi1/phi0).2;s02=h-2*h*(phi0-phi2).2/(phi0).2;phi3=linspace(deg2rad(90),deg2rad(190),100
12、0); s03=h;phi4=linspace(deg2rad(190),deg2rad(240);s04=h*(1-35*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).4+84*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).5-70*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).6+20*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).7);phi5=linspace(deg2rad(240),deg2rad(360),1000) s0
13、5=0;%ds/fs1=4*h/(deg2rad(90)2*phi1;s2=4*h*(phi0-phi2)/(deg2rad(90)2;s3=0;s4=-h/deg2rad(50)*140*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).3-420*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).4+420*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).5-140*(phi4-deg2rad(90)-deg2rad(100)/deg2rad(50).6;s5=0; x=s1,s2,s3,s4,s5;y=s01,s02,s03,s04,s05;plot(x,y)hold on %确定凸轮半径和偏距k1=tan(pi/2-30/180*pi);k2=-tan(pi/2-60/180*pi);y1=-k1*x+y;y2=-k2*x+y;y11=min(y1);y22=min(y2);x1=linspace(-5,50,10000);x2=linspace(-95,55,10000);d1=k1*x1+y11;d2=k2*x2+y22;x3=linspace(0,45,500);x3=12;y3=-48;plot(l
