哈工大机械原理大作业2-凸轮27
机械原理大作业二-凸轮机构设计..
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机械设计制造及其自动化班级：1208104完成者：郑鹏伟学号：**********指导教师：林琳刘福利设计时间：2014年6月4日哈尔滨工业大学一、 设计题目：凸轮的机构运动简图如下图所示:序 号 升程 （mm ） 升程运动角（°） 升程运 动规律 升程 许用压力角（°）回程运动角 （°） 回程运动规律回程许用压力角（°） 远休止角（°） 近休止角（°） 14 90120余弦 加速 度3590等减等加速657575二、 凸轮推杆升程、回程运动方程及推杆位移，速度加速度线图：（1）凸轮推杆升程运动方程： 根据题意知：00120759075s s Φ=Φ='Φ='Φ=（1）从动件升程运动方程（设为1rad sω=）122212s [1cos()]2sin()2cos()2h h h a πϕπωπυϕπωπϕ=-Φ=ΦΦ=ΦΦ（2）从动件远休止运动方程在远休止s Φ段，即213312πϕπ≤≤时，90,0,0s h mm a υ====。
（3）从动件回程运动方程升程段采用等减等加运动规律，运动方程为：①当回程0002s s ϕ'ΦΦ+Φ≤≤Φ+Φ+134()123πϕπ≤≤时： 20s 201022122[-+]4[()]4s hs h h h a ϕωυϕω=-ΦΦ'Φ=--Φ+Φ'Φ=-'Φ（）②当回程0002s s ϕ'Φ'Φ+Φ+≤≤Φ+Φ+Φ419()312πϕπ≤≤时：20020100202122[)]4[)]4s s h s h h a ϕωυϕω'=Φ+Φ+Φ-'Φ'=-Φ+Φ+Φ-'Φ='Φ（（（4）从动件近休止运动方程在近休止s 'Φ段，即19212πϕπ≤≤时，s 0,0,0a υ===。
机械原理大作业——凸轮.docx
大作业（二）凸轮机构设计题号：6班级：姓名：学号：同组者：成绩：完成时间：目录一凸轮机构题目要求 (1)二摆杆的运动规律及凸轮轮廓线方程 (2)三计算程序 (3)四运算结果及凸轮机构图 (9)4.1 第一组（A组）机构图及计算结果 (9)4.2 第二组（B组）机构图及计算结果 (14)4.3 第三组（C组）机构图及计算结果 (19)五心得体会 (24)第一组（A组） (24)第二组（B组） (24)第三组（C组） (24)六参考资料 (25)附录程序框图 (26)一凸轮机构题目要求（摆动滚子推杆盘形凸轮机构）题目要求：试用计算机辅助设计完成下列偏置直动滚子推杆盘形凸轮机构或摆动滚子推杆盘形凸轮机构的设计，已知数据如下各表所示。
凸轮沿逆时针方向作匀速转动。
表一摆动滚子推杆盘形凸轮机构的已知参数题号初选的基圆半径R0/mm机架长度Loa/mm摆杆长度Lab/mm滚子半径Rr/mm推杆摆角φ许用压力角许用最小曲率半径[ρamin][α1] [α2]A 15 60 55 10 24°35°70°0.3RrB 20 70 65 14 26°40°70°0.3RrC 22 72 68 18 28°45°65°0.35Rr 要求：1）凸轮理论轮廓和实际轮廓的坐标值2）推程和回程的最大压力角，及凸轮对应的转角3）凸轮实际轮廓曲线的最小曲率4）半径及相应凸轮转角5）基圆半径6）绘制凸轮理论廓线和实际廓线7）计算点数：N：72～120推杆运动规律：1)推程运动规律：等加速等减速运动2)回程运动规律：余弦加速度运动二摆杆的运动规律及凸轮轮廓线方程1）推程：1，运动规律：等加速等减速运动；2，轮廓线方程：A：等加速推程段设定推程加速段边界条件为：在始点处δ=0，s=0，v=0。
在终点处。
整理得：（注意：δ的变化范围为0~δ0/2。
哈工大机械原理大作业直动从动件盘形凸轮机构满分完美版哈尔滨工业大学
end
%曲率半径
dx=diff(x);
dx(36001)=0;
dy=diff(y);
dy(36001)=0;
dydx=dy./dx;
ddy=diff(dydx);
ddy(36001)=0;
ddy=ddy./dx;
~7~
xlabel('凸轮转角φ/°');
ylabel('从动件加速度 a/(mm/s^2)');
title('加速度');
subplot(2,3,4);
plot(dsdPhi,s);
axis equal
axis([-30 50 -50 30]);
hold on
plot(x1,f1,'r');
hold on
axis equal
hold on
plot(x,y,'k')
legend('理论廓线',-1);
hold on
%轮廓图
%实际廓线
%理论廓线
~8~
哈尔滨工业大学
直动从动件盘型凸轮机构设计说明书
plot(x0,y0,':')
hold on
%基圆
plot(xe,ye,'k:')
%偏距圆
legend('实际廓线','理论廓线','基圆','偏距圆',-1);
subplot(2,3,1);
plot(Phi,s);
grid on
axis([0 360 0 27]);
哈工大机械原理大作业-凸轮机构设计
哈工大机械原理大作业-凸轮机构设计(第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 ，则在直线上或其右下方取凸轮轴心时，可使[]αα≤。
哈工大机械原理大作业——凸轮——2号
哈工大机械原理大作业——凸轮——2号————————————————————————————————作者：————————————————————————————————日期：Harbin Institute of Technology机械原理大作业课程名称：机械原理设计题目：凸轮机构设计一、设计题目(1)凸轮机构运动简图：(2)凸轮机构的原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角14 90°120°余弦加速度35°90°3-4-5多项式65°80°70°(1) 推杆升程、回程运动方程如下：A.推杆升程方程：设为1rad sω=升程位移为：()()1cos451cos1.52hsπψψψ⎡⎤⎛⎫=-=-⎢⎥⎪Φ⎝⎭⎣⎦23ψπ≤≤升程速度为：()()1100sin67.5sin1.52hvπωπψψωψ⎛⎫==⎪ΦΦ⎝⎭23ψπ≤≤升程加速度为：()()2221100cos101.25cos1.52haπωπψψωψ⎛⎫==⎪ΦΦ⎝⎭23ψπ≤≤B.推杆回程方程：回程位移为：()()345111110156s h T T T ψ⎡⎤=--+⎣⎦1029918ψπ≤≤ 回程速度为：()()22111103012h v T T T ωψ=--+'Φ 1029918ψπ≤≤ 回程加速度为：()()221111260132h a T T T ωψ=--+'Φ 1029918ψπ≤≤其中：()010s T ψ-Φ+Φ='Φ1029918ψπ≤≤ (2) 利用Matlab 绘制推杆位移、速度、加速度线图 A. 推杆位移线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300); x3=linspace(10*pi/9,29*pi/18,300); x4=linspace(29*pi/18,2*pi,300); T1=(x3-10*pi/9)/(pi/2); s1=45*(1-cos(1.5*x1)) s2=90;s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5)); s4=0;plot(x1,s1,'r',x2,s2,'r',x3,s3,'r',x4,s4,'r') xlabel('角度ψ/rad'); ylabel('位移s/mm') title('推杆位移线图') gridaxis([0,7,-10,100]) 得到推杆位移线图：B.推杆速度线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);v1=67.5*1*sin(1.5*x1);v2=0;v3=-30*90*1*T1.^2/(pi/2).*(1-2*T1+T1.^2);v4=0;plot(x1,v1,'r',x2,v2,'r',x3,v3,'r',x4,v4,'r') xlabel('角度ψ/rad');ylabel('速度v/(mm/s)')title('推杆速度线图')Grid得到推杆速度线图：C.推杆加速度线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);a1=101.25*1.^2.*cos(1.5*x1);a2=0;a3=-60.*90.*T1./((pi/2).^2).*(1-3*T1+2*T1.^2); a4=0;plot(x1,a1,'r',x2,a2,'r',x3,a3,'r',x4,a4,'r') xlabel('角度ψ/rad');ylabel('加速度a/')title('推杆加速度线图')Grid得到推杆加速度线图：三、凸轮机构的ds/dψ-s线图，并依次确定凸轮的基圆半径和偏距.1、凸轮机构的ds/dψ--s线图：x1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);s1=45*(1-cos(1.5*x1))s2=90;s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5));s4=0;v1=67.5*1*sin(1.5*x1);v2=0;v3=-30*90*1*T1.^2/(pi/2).*(1-2*T1+T1.^2);v4=0;plot(v1,s1,'r',v2,s2,'r',v3,s3,'r',v4,s4,'r')xlabel('ds/dψ');ylabel('(位移s/mm)')title('ds/dψ—s曲线')gridaxis([-120,80,-10,100])得到ds/dψ—s曲线：2、确定凸轮的基圆半径和偏距：在dssdϕ-线图中，右侧曲线为升程阶段的类速度-位移图，作直线D t d t与其相切，且与位移轴正方向呈夹角[α1]=350, 故该直线斜率：32sin2=tan5533cos2okϕϕ⨯=⨯通过编程求其角度。
(完整word版)机械原理大作业2-1120810417-凸轮
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院班级：1208104完成者：学号：1120810417指导教师：林琳刘福利设计时间：2014年6月2日哈尔滨工业大学一、设计题目如下图所示为直动从动件盘形凸轮机构，据此设计该凸轮机构：二、原始参数 序号升程升程运动角 升程运动规律 升程许用压力角 回程运动角 回程运动规律 回程许用压力角 远休止角 近休止角 15 90mm150°正弦加速度30°100°余弦加速度60°55°55°三、推杆升程方程和推杆回程方程： 在这里取ω=1rad/s. （1）推杆升程方程：650,)512sin(215690)(πφφππφφ≤≤⎥⎦⎤⎢⎣⎡-=s650),512cos(108)(πφφφπφν≤≤-=650,512sin 2.259)(πφφπφ≤≤=a（2）推杆回程方程：36613641,)05.059cos(145)(πφππφφ≤≤⎥⎦⎤⎢⎣⎡-+=sω36613641,)05.059sin(181)(πφππφφν≤≤⎥⎦⎤⎢⎣⎡---= 36613641),05.059cos(8.145)(≤≤--=φππφφa四、matlab 程序及曲线图像注：每一段都为完整程序，可直接运行。
1.推杆位移曲线clear allp1=0:pi/360:(5*pi/6-pi/360); w=1;s1=90*(6*p1/(5*pi)-1/(2*pi)*sin(12*p1/5)); p2=5*pi/6:pi/360:(41*pi/36-pi/360); s2=90*ones(1,length(p2));p3=41*pi/36:pi/360:(61*pi/36-pi/360); s3=45*(1+cos(9*p3/5-1*pi/20)); p4=61*pi/36:pi/360:2*pi; s4=0*p4;p=[p1,p2,p3,p4]; s=[s1,s2,s3,s4];plot(p,s)xlabel('Φ(角度)');ylabel('S(位移)'); title('推杆位移曲线');2.推杆速度曲线clear allp1=0:pi/360:(5*pi/6-pi/360);w=1;v1=108*w/pi*(1-cos(12*p1/5));p2=5*pi/6:pi/360:(41*pi/36-pi/360);v2=0*p2;p3=41*pi/36:pi/360:(61*pi/36-pi/360);v3=-81*w*sin(9*p3/5-1*pi/20);p4=61*pi/36:pi/360:2*pi;v4=0*p4;p=[p1,p2,p3,p4];v=[v1,v2,v3,v4];plot(p,v)xlabel('Φ(角度)');ylabel('V(速度)'); title('推杆速度曲线');3.推杆加速度曲线clear allp1=0:pi/360:(5*pi/6-pi/360);w=1;a1=36*36*w^2/5/pi*sin(12*p1/5);p2=5*pi/6:pi/360:(41*pi/36-pi/360);a2=0*p2p3=41*pi/36:pi/360:(61*pi/36-pi/360);a3=-18*81*w^2/10*cos(9*p3/5-1*pi/20);p4=61*pi/36:pi/360:2*pi;a4=0*p4;p=[p1,p2,p3,p4];a=[a1,a2,a3,a4];plot(p,a)xlabel('Φ(角度)');ylabel('a(加速度)'); title('推杆加速度曲线');4.凸轮机构的ds/dφ-s线图clear allp1=0:pi/360:(5*pi/6-pi/360);w=1;s1=90*(6*p1/(5*pi)-1/(2*pi)*sin(12*p1/5)); p2=5*pi/6:pi/360:(41*pi/36-pi/360);s2=90*ones(1,length(p2));p3=41*pi/36:pi/360:(61*pi/36-pi/360);s3=45*(1+cos(9*p3/5-1*pi/20));p4=61*pi/36:pi/360:2*pi;s4=0*p4;p=[p1,p2,p3,p4];s=[s1,s2,s3,s4];p1=0:pi/360:(5*pi/6-pi/360);w=1;v1=108*w/pi*(1-cos(12*p1/5));p2=5*pi/6:pi/360:(41*pi/36-pi/360);v2=0*p2;p3=41*pi/36:pi/360:(61*pi/36-pi/360);v3=-81*w*sin(9*p3/5-1*pi/20);p4=61*pi/36:pi/360:2*pi;v4=0*p4;p=[p1,p2,p3,p4]; v=[v1,v2,v3,v4]; vx=-v; hold on plot(vx,s)%直线Dtdty=-100:0.01:100; x=-69; hold onplot(x,y,'-r'); % 直线Dt’dt’ x=-100:0.01:100; y=-0; hold onplot(x,y,'-r'); grid on hold offtitle('ds/d φ-s 曲线');曲线为升程阶段的类速度-位移图，根据升程压力角与回城压力角做直线与其相切，, 其直线斜率分别为：K 1=)30150tan(+=0 K 2=)60150tan(-为∞；两直线方程为：}{0,69=-=y x进而确定凸轮偏距和基圆半径：在轴心公共许用区内取轴心位置，能够满足压力角要求，由图可得:取s0=200mm ，e=30；r0=（2002 ＋502）1/2=206.2mmclear allp1=0:pi/360:(5*pi/6-pi/360);w=1;s1=90*(6*p1/(5*pi)-1/(2*pi)*sin(12*p1/5));p2=5*pi/6:pi/360:(41*pi/36-pi/360);s2=90*ones(1,length(p2));p3=41*pi/36:pi/360:(61*pi/36-pi/360);s3=45*(1+cos(9*p3/5-1*pi/20));p4=61*pi/36:pi/360:2*pi;s4=0*p4;p=[p1,p2,p3,p4];s=[s1,s2,s3,s4];s0=200;e=30;x=(s0+s).*cos(p)-e*sin(p);y=(s0+s).*sin(p)+e*cos(p);plot(x,y)title('凸轮理论轮廓');6.凸轮实际轮廓工作轮廓曲率半径ρ、理论轮廓曲率半径ρ与滚子半径r三者存在如下关系aρa=ρ+r，不妨最终设定滚子半径为30mm，这时滚子与凸轮间接触应力最小，可提高凸轮寿命。
机械原理大作业2-凸轮机构设计
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院姓名：学号：班级：指导教师：1.设计题目设计直动从动件盘形凸轮机构，其原始参数见表1 表一：凸轮机构原始参数升程（mm ) 升程运动角（º）升程运动规律升程许用压力角（º）回程运动角（º）回程运动规律回程许用压力角（º）远休止角（º）近休止角（º）30 70 等加等减速30 170 正弦加速度60 100 1202.凸轮推杆运动规律（1）推程运动规律（等加速等减速运动）推程035≤ϕ0≤推程0070≤ϕ35≤（2）回程运动规律（正弦加速度）回程00240≤ϕ170≤开始输入初始参数运行各部分程序输出1 输出2 输出3 输出4 结束从动件位移、速度、加速度曲线ds/dψ-s曲线，确定基圆半径和偏距理论轮廓线上的压力角和曲率半径图绘制理论轮廓线和实际轮廓线3.运动线图及凸轮s d ds -φ线图采用Matlab 编程，其所有源程序见附页：令可得运动规律图如下：1.凸轮的基圆半径和偏距以ds/dfψ-s图为基础，可分别作出三条限制线（推程许用压力角的切界限D t d t,回程许用压力角的限制线D t'd t'，起始点压力角许用线B0d'')，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。
得图如下：得最小基圆对应的坐标位置大约为（20，-35）经计算取偏距e=20mm，r0=40.3mm.2.绘制理论轮廓线上的压力角曲线和曲率半径曲线针对凸轮转向及推杆偏置，令N1=1凸轮逆时针转；N2=1偏距为正。
压力角数学模型：曲率半径数学模型：)/)(/()/)(/(])/()/[(22222/322ϕϕϕϕϕϕρd x d d dy d y d d dx d dy d dx -+= 其中：)sin(])/[()cos(])/(2[/102212122ϕϕϕϕϕN s s d s d N e N N d ds d x d --+-=)cos(])/[(）（sin ])/(2[/202212122ϕϕϕϕϕN s s d s d N e N N d ds d y d --+--=3.凸轮理论廓线和实际廓线理论廓线数学模型：ϕϕϕϕ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 tt+-=++=其中rt 为确定的滚子半径。
哈工大机械原理大作业2凸轮机构设计
机械原理大作业（二）作业名称：凸轮机构设计设计题目：23题院系：班级：设计者：学号：指导教师：设计时刻：哈尔滨工业大学机械设计1.运动分析题目：设计直动从动件盘形凸轮机构，其原始参数见下表2.确信凸轮机构推杆升程、回程运动方程（设定角速度为ω=10 rad/s）升程：0°< Φ < 120°由公式可得：s=60-60*cos(3*Φ/2)；v=90*ω*sin(3*Φ/2)；a=135*ω2 *cos(3*Φ/2)；远停止：120°< Φ < 200°由公式可得：s=120；v=0；a=0；回程：200°< Φ < 290°由公式可得：s=h[1-(10T23-15T24+6T25)]v=(-30hω1/Φ0')T22(1–2T2+T22)a=(-60hω12/Φ0'2)T2(1–3T2+2T22)式中：T2=(Φ-Φ0-Φs)/ Φ0'近停止： 290°< Φ < 360°由公式可得：s=0；v=0；a=0；3.绘制推杆位移、速度、加速度线图(设ω=10rad/s）1) 推拉位移曲线代码：%推杆位移曲线;x=0:(pi/1000):(2*pi/3);s1=60-60*cos(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);s2=120;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;s3=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));m=(29*pi/18):(pi/1000):(2*pi);s4=0;plot(x,s1,'b',y,s2,'b',z,s3,'b',m,s4,'b'); xlabel('角度(rad)');ylabel('行程(mm)');title('推杆位移曲线');grid;2）推杆速度曲线代码：%推杆速度曲线;w=10;x=0:(pi/1000):(2*pi/3);v1=90*w*sin(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);v2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;v3=(-30*120*w/(pi/2))*T2.^2.*(1-2*T2.^2+T2.^2); % v3=-120*w*sin(2*z-20*pi/9);m=(29*pi/18):(pi/1000):(2*pi);v4=0;plot(x,v1,'r',y,v2,'r',z,v3,'r',m,v4,'r'); xlabel('角度(rad)');ylabel('速度(mm/s)');title('推杆速度曲线(w=10rad/s)');grid;3）凸轮推杆加速度曲线代码：%凸轮推杆加速度曲线;w=10;x=0:(pi/1000):(2*pi/3);a1=135*w^2*cos(3*x/2);y=(2*pi/3):(pi/1000):(10*pi/9);a2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;a3=(-60*120*w^2/(pi/2)^2)*T2.*(1-3*T2.^2+2*T2.^2); m=(29*pi/18):(pi/1000):(2*pi);a4=0;plot(x,a1,'m',y,a2,'m',z,a3,'m',m,a4,'m');xlabel('角度(rad)');ylabel('加速度(mm/s^2)');title('凸轮推杆加速度曲线(w=10rad/s)');grid;4）绘制凸轮机构的dd/dd−d线图，并依次确信凸轮的基圆半径和偏距代码：%dd/dd−d线图，确信e，s0;x=0:(pi/1000):(2*pi/3);s1=60-60*cos(1.5*x);ns1=90*sin(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);s2=120;ns2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;s3n=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));ns3=-120*10*3*T2.^2+120*15*4*T2.^3-120*6*5*T2.^4 ;m=(29*pi/18):(pi/1000):(2*pi);s4=0;ns4=0;x1=0:pi/36000:pi/2;s1n=60-60*cos(1.5*x1);v1=90*sin(1.5*x1);m1=diff(s1n);%求切线1n1=diff(v1);z=m1./n1;for i=1:length(z);if abs(z(i)+tan(-55*pi/180))<0.001;breakendendb11=s1n(i)-z(i)*v1(i);x1=-300:200;y01=z(i)*x1+b11;%切线1k1=z(i);plot(x1,y01)x3=10*pi/9:pi/36000:14*pi/9;%求切线2s3n=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));v3=-120*10*3*T2.^2+120*15*4*T2.^3-120*6*5*T2.^4 ;m3=diff(s3n);n3=diff(v3);p=m3./n3;for o=1:length(p);if abs(p(o)-tan(-25*pi/180))<0.01;breakendendo;b33=s3n(o)-p(o)*v3(o);x3=-300:700;y03=p(o)*x3+b33;%切线2plot(x3,y03);sym uv[u,v]=solve('u= 1.4281*v-81.7665','u=-0.4663*v-59.6715');%v=11.66332347972972972972972972973 x%u=-65.110107738597972972972972972973 yplot(ns1,s1,'m',ns2,s2,'b',ns3,s3n,'b',ns4,s4,'b',x1,y01,'g',x3,y03,'g',v,u,'*'); xlabel('ds/d¦µ');ylabel('S');axis([-300,200,-300,300]);title('s0,e 的确信');grid;确信凸轮基圆半径与偏距：偏距e=90mm，d020mm；基圆半径为d0=150mm。
哈工大—机械原理凸轮大作业
一、题目要求及机构运动简图如图1所示直动从动件盘形凸轮机构。
其原始参数见表1。
图一凸轮运动简图表一凸轮原始参数二、计算流程框图凸轮机构分析建立数学模型位移方程速度方程加速度方程速度线图位移线图加速线图ds/dΨ-s曲线升程压力角回程压力角确定轴向及基圆半径压力角图确定滚子半径实际轮廓理论轮廓轮廓图结束三、建立数学模型1.从动件运动规律方程首先，由于设计凸轮轮廓与凸轮角速度无关，所以不妨设凸轮运动角速度为w = 1rad/s。
（1）推程运动规律(0 < φ < 90°)s=h2×[1−cos(πφ0×φ）]v=πhw2φ0×sin(πφ0×φ）a=π2ℎw22φ02×cos(πφ0×φ）式中：h=65mm，Φ0=π/2（2）远休程运动规律(90°< φ < 190°)s = 65mmv = 0a = 0（3）回程运动规律(190°< φ < 240°)s1=h−h4+π∗(π∗(φ−φ0−φs)φ0′−sin(4∗π∗φ−φ0−φsφ0′)4)(190°< φ < 196.25°)s2=h−h4+π∗(2+π∗(φ−φ0−φs)φ0′−9∗sin(π3+4∗π∗φ−φ0−φs3∗φ0′)4)(196.25°< φ < 233.75°)s3=h−h4+π∗(4+π∗(φ−φ0−φs)φ0′−sin(4∗π∗φ−φ0−φsφ0′)4)(233.75°< φ < 240°)回程运动中的速度和加速度为位移对时间t的倒数：v=ds dta=dv dt（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 线的右下方为选择凸轮轴心的许用区。
哈工大机械原理大作业凸轮
1 / 17机械原理大作业课程名称：机械原理 设计题目：直动从动件盘形凸轮机构 院 系： 机电学院 班 级： 完 成 者： 学 号： 指导教师： 设计时间：哈尔滨工业大学题目：如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
凸轮运动分为五个阶段1.升程阶段0~5034500010156s h ϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦23410000306030h v ωϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦23212000060180120h a ωϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦2.远休止50~1503.回程等加速150~195()20'202sh s h ϕ=--Φ-ΦΦ()10'204s h s ωϕ=--Φ-ΦΦ21'204h v ω=Φ4.回程等减速195~240()2'00'202s h s ϕ=Φ+Φ+Φ-Φ()'100'204s h v ωϕ=-Φ+Φ+Φ-Φ21'204h v ω=-Φ5.近休止 240~360其中，0Φ 推程运动角 s Φ 远休止角 '0Φ 回程运动角使用Matlab 实现1.从动件位移%用fi1,fi2,fi3,fi4，fi5代替转角h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;s1=h*(10*(fi1/fio1).^3-15*(fi1/fio1).^4+6*(fi1/fio1).^5); hold onplot(180*fi1/pi,s1);fi2=fio1:0.01:fio1+fis;s2=h;hold onplot(180*fi2/pi,s2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;s3=h-2*h/fio2^2*(fi3-fio1-fis).^2;hold onplot(180*fi3/pi,s3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;s4=2*h/fio2^2*(fio1+fis+fio2-fi4).^2;hold onplot(180*fi4/pi,s4);fi5=fio1+fis+fio2:0.001:2*pi;s5=0;hold onplot(180*fi5/pi,s5);title('位移');xlabel('φ/度'),ylabel('s/mm');grid onhold off2.从动件速度%用fi1,fi2,fi3,fi4，fi5代替转角%设角速度为1rad/sw=1;h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;v1=h*w/fio1*(30*(fi1/fio1).^2-60*(fi1/fio1).^3+30*(fi1/fio1).^4); hold onplot(180*fi1/pi,v1);fi2=fio1:0.01:fio1+fis;v2=0;hold onplot(180*fi2/pi,v2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;v3=-4*h*w/fio2^2*(fi3-fio1-fis);hold onplot(180*fi3/pi,v3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;v4=-4*h*w/fio2^2*(fio1+fis+fio2-fi4);hold onplot(180*fi4/pi,v4);fi5=fio1+fis+fio2:0.001:2*pi;v5=0;hold onplot(180*fi5/pi,v5);title('速度');xlabel('φ/度'),ylabel('v/(mm/s)');grid onhold off3.从动件加速度%用fi1,fi2,fi3,fi4，fi5代替转角%设角速度为1rad/sw=1;h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;a1=h*w^2/fio1*(60*(fi1/fio1)-180*(fi1/fio1).^2+120*(fi1/fio1).^3); hold onplot(180*fi1/pi,a1);fi2=fio1:0.01:fio1+fis;a2=0;hold onplot(180*fi2/pi,a2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;a3=-4*h*w/fio2^2;hold onplot(180*fi3/pi,a3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;a4=4*h*w/fio2^2;hold onplot(180*fi4/pi,a4);fi5=fio1+fis+fio2:0.001:2*pi;a5=0;hold onplot(180*fi5/pi,a5);title('加速度');xlabel('φ/度'),ylabel('v/(mm/s^2)');grid onhold off4.sdsdϕ-线图%用fi1,fi2,fi3,fi4，fi5代替转角h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;s1=h*(10*(fi1/fio1).^3-15*(fi1/fio1).^4+6*(fi1/fio1).^5);ds1=h*(3*10*(fi1/fio1).^2/fio1-4*15*(fi1/fio1).^3/fio1+5*6*(fi1/fio1).^4/fio1);hold onplot(ds1,s1);fi2=fio1:0.01:fio1+fis;s2=h+0*fi2;ds2=0*fi2;hold onplot(ds2,s2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;s3=h-2*h/fio2^2*(fi3-fio1-fis).^2; ds3=-2*h/fio2^2*2*(fi3-fio1-fis); hold onplot(ds3,s3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2; s4=2*h/fio2^2*(fio1+fis+fio2-fi4).^2; ds4=-2*h/fio2^2*2*(fio1+fis+fio2-fi4); hold onplot(ds4,s4);fi5=fio1+fis+fio2:0.001:2*pi; s5=0+0*fi5; ds5=0*fi5; hold onplot(ds5,s5);title('ds/d φ-s');xlabel('ds/d φ(mm/rad)'),ylabel('s(mm)'); grid on hold off5.凸轮轴心位置的确定凸轮压力角的正切值s s ed ds +-=0/tan ϕα，右侧为升程，作与s 轴夹6π角等于升程许用压力角的切界线t t dD ，则在直线上或其下方取凸轮轴心时，可使[]αα≤，同理右侧回程，作与s 轴夹角等于回程许用压力角3π的切界线''t t d D ，则在直线上或其下方取凸轮轴心时，可使[]αα≤。
- 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机械原理大作业说明书课程名称：机械原理设计题目：凸轮机构设计院系：班级：设计者：学号：指导教师：设计时间：哈尔滨工业大学一、设计题目 如右图所示直动从动件盘形凸轮机构，选择一组凸轮机构的原始参数，据此设计该凸轮机构。
二. 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图凸轮推杆升程运动方程：)]512sin(2156[130s ϕππϕ-=)512sin(4.374)]512cos(1[156v 211ϕπϕπωω=-=aclear;clc;hu=pi/180; %输入题目要求的参数 h=130; %从动件最大线位移 w1=10; %假设的凸轮 fi0=150*hu; %推程运动角 fis=100*hu; %远休止角 fi02=70*hu; %回程运动角 fis2=40*hu; %近休止角 %开始计算 %推程计算xfi1=0:0.01:fi0; %凸轮转角 T1=xfi1/fi0;s1=h*(10*T1.^3-15*T1.^4+6*T1.^5);v1=30*h*w1*T1.^2.*(1-2*T1+T1.^2)/fi0;a1=60*h*w1^2*T1.*(1-3*T1+2*T1.^2)/fi0^2; xfi2=(fi0+fis):0.01:(fi0+fis+fi02);%远休程计算xfi3=fi0:0.01:(fi0+fis);s3=h*ones(1,length(xfi3));v3=zeros(1,length(xfi3));a3=zeros(1,length(xfi3));%回程计算T=xfi2-fi0-fis;s2=h*(1-T/fi02+1/(2*pi)*sin(2*pi*T/fi02)); v2=-h*w1*(1-cos(2*pi*T/fi02))/fi02;a2=-2^pi*h*w1^2*sin(2*pi*T/fi02)/fi02^2;%近休程计算xfi4=(2*pi-fis2):0.01:2*pi;s4=zeros(1,length(xfi4));v4=zeros(1,length(xfi4));a4=zeros(1,length(xfi4));% t表示转角，令ω1=1t=0:0.01:5*pi/6;%升程阶段v=156*1*[1-cos(12*t/5)]/pi hold onplot(t,v);t= 5*pi/6:0.01:pi;%远休止阶段v=0hold onplot(t,v);t=pi:0.01:14*pi/9;%回程阶段v=-117*1*sin(9*(t-pi)/5) hold onplot(t,v);t=14*pi/9:0.01:2*pi;%近休止阶段v=0hold ont=0:0.001:5*pi/6;a=374.4*sin(12*t/5)/pi;hold onplot(t,a);t=5*pi/6:0.01:pi;a=0;hold onplot(t,a);t=pi:0.001:14*pi/9;a=-210.6*cos(9*(t-pi)/5); hold onplot(t,a);t=14*pi/9:0.001:2*pi;a=0;hold on三. 绘制凸轮机构的sdds-ϕ线图% t表示转角，x(横坐标)表示速度ds/dφ，y（纵坐标）表示位移s t=0:0.001:5*pi/6;% 升程阶段x= 156*1*(1-cos(12*t/5))/pi;y= 130*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5));hold onplot(x,y,'-r');t= 5*pi/6:0.01:pi;%远休止阶段x=0;y=130;hold onplot(x,y,'-r');t=pi:0.001:14*pi/9;% 回程阶段x=-117*1*sin(9*(t-pi)/5); y=65*(1+cos(9*(t-pi)/5)); hold onplot(x,y,'-r');t=14*pi/9:0.01:2*pi;%近休止阶段x=0;y=0;hold onplot(x,y,'-r');grid onhold off四.按许用压力角确定凸轮基圆半径和偏距1. 求切点转角（1）在图-4中，右侧曲线为升程阶段的类速度-位移图，作直线D t d t与其相切，且位移轴正方向呈夹角[α1]=300,则切点处的斜率与直线D t d t的斜率相等，因为k Dtdt=tan300,右侧曲线斜率可以表示为错误！未找到引用源。
,所以，错误！未找到引用源。
，通过编程求其角度。
编码：求得转角t =1.1123进而求的切点坐标（x,y）=(93.8817, 45.8243)（2）在图-4中，左侧曲线为回程阶段的类速度-位移图，作直线D’t d’t与其相切，且位移轴正方向呈夹角[α1]=600,则切点处的斜率与直线D’t d’t的斜率相等，因为k Dtdt=tan300同理求得切点坐标（x,y）=( -110.0654, 42.3144)2. 确定直线方程直线D t d t：y =tan(pi/3)(x-93.8817)=45.8243;直线D t’d t’: y =-tan(pi/3)(x+84.3144)+110.0654;3. 绘图确定基圆半径和偏距% 直线Dtdtx=-125:1:150;y= tan(pi/3)*(x-93.8798)+45.8243;hold onplot(x,y);% 直线Dt’dt’x=-125:1:150;y=-tan(pi/6)*(x+110.0654)+34.3144;hold onplot(x,y);%直线Ddx=0:1:150;y=tan(2*pi/3)*x;hold onplot(x,y);t=0:0.001:5*pi/6;% 升程阶段x= 156*1*[1-cos(12*t/5)]/pi;y= 130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]; hold onplot(x,y,'-r');t= 5*pi/6:0.01:pi;%远休止阶段x=0;y=130;hold onplot(x,y,'-r');t=pi:0.001:14*pi/9;% 回程阶段x=-117*1*sin(9*(t-pi)/5);y=65*[1+cos(9*(t-pi)/5)];hold onplot(x,y,'-r');t=14*pi/9:0.01:2*pi;%近休止阶段x=0;y=0;hold onplot(x,y,'-r');grid onhold off如图，在这三条直线所围成的公共许用区域，只要在公共许用区域内选定凸轮轴心O的位置，凸轮基圆半径r0和偏距e就可以确定了。
现取轴心位置为x=20，y=-125，则可得偏距e=20，基圆半径=127五.绘制凸轮理论轮廓线编码：%凸轮的理论轮廓，t表示转角，x表示横坐标，y表示纵坐标t=0:0.0001:5*pi/6;x=(125+130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]).*cos(t)-20*sin(t);y=(125+130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]).*sin(t)+20*cos(t);hold onplot(x,y);t= 5*pi/6:0.0001:pi;x=(125+130).*cos(t)- 20*sin(t);y=(125+130).*sin(t)+ 20*cos(t);hold onplot(x,y);t=pi:0.0001:14*pi/9;x=(125+65*[1+cos(9*(t-pi)/5)]).*cos(t)- 20*sin(t);y=(125+65*[1+cos(9*(t-pi)/5)]).*sin(t)+ 20*cos(t);hold onplot(x,y);t= 14*pi/9:0.0001:2*pi;x=(125).*cos(t)- 20*sin(t);y=(125).*sin(t)+ 20*cos(t);hold onplot(x,y);%基圆t=0:0.001:2*pi;x=20.1074*cos(t);y=20.1074*sin(t);hold onplot(x,y);% 偏心圆t=0:0.001:2*pi;x=13.3509*cos(t);y=13.3509*sin(t);hold onplot(x,y);grid onhold off六、在理论廓线上分别绘出基圆与偏距圆：编码：%凸轮的理论轮廓，t表示转角，x表示横坐标，y表示纵坐标t=0:0.0001:5*pi/6;x=(125+130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]).*cos(t)-20*sin(t);y=(125+130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]).*sin(t)+20*cos(t);hold onplot(x,y);t= 5*pi/6:0.0001:pi;x=(125+130).*cos(t)- 20*sin(t);y=(125+130).*sin(t)+ 20*cos(t);hold onplot(x,y);t=pi:0.0001:14*pi/9;x=(125+65*[1+cos(9*(t-pi)/5)]).*cos(t)- 20*sin(t); y=(125+65*[1+cos(9*(t-pi)/5)]).*sin(t)+ 20*cos(t); hold onplot(x,y);t= 14*pi/9:0.0001:2*pi;x=(125).*cos(t)- 20*sin(t);y=(125).*sin(t)+ 20*cos(t);hold onplot(x,y);%基圆t=0:0.001:2*pi;x=127*cos(t);y=127*sin(t);hold onplot(x,y);% 偏心圆t=0:0.001:2*pi;x=20*cos(t);y=20*sin(t);hold onplot(x,y);grid onhold off七.确定滚子半径1. 绘制曲率半径图% 凸轮理论轮廓半径，t表示转角，p表示曲率半径,%dxi表示dx/dφ, dyi表示dy/dφ,i=1,2,3,4h=130; %升程t0=pi*5/6; % 升程角t01=pi*5/9; % 回程角ts=pi/6; %远休止角ts1=pi*4/9; %近休止角e=20; %偏距s0=125;% 升程阶段t=linspace(0,pi*5/6,1000);s=h*(t/t0-sin(2*pi*t/t0)/(2*pi));dx1 =(h/t0-h*cos(2*pi*t/t0)).*cos(t)-(s0+s).*sin(t)- e*cos(t); dy1=(h/t0-h*cos(2*pi*t/t0)).*sin(t)+(s0+s).*cos(t)- e*sin(t); p=sqrt(dx1.^2+dy1.^2);hold onplot(t,p);% 远休止阶段t=linspace(pi*5/6,pi,1000);s=h;dx2 =- sin(t).*(s + s0) - e*cos(t);dy2 =cos(t).*(s + s0) - e*sin(t);p=sqrt(dx2.^2+dy2.^2);hold onplot(t,p);% 回程阶段t=linspace(pi,pi*14/9,1000);s=0.5*h*(1+cos(pi*(t-(t0+ts))/t01));dx3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*cos(t)- sin(t).*(s + s0) - e*cos(t); dy3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*sin(t)+ cos(t).*(s + s0) - e*sin(t); p=sqrt(dx3.^2+dy3.^2);hold onplot(t,p);%近休止阶段t=linspace(pi*14/9,pi*2,1000);s=0;dx4 =- sin(t).*(s + s0) - e*cos(t);dy4 =cos(t).*(s + s0) - e*sin(t);p=sqrt(dx4.^2+dy4.^2);hold onplot(t,p);hold offtitle('曲率半径ρ','FontSize',20);grid on八. 绘制实际轮廓线% 凸轮理论轮廓半径，t表示转角，p表示曲率半径,%dxi表示dx/dφ, dyi表示dy/dφ,i=1,2,3,4h=130; %升程t0=pi*5/6; % 升程角t01=pi*5/9; % 回程角ts=pi/6; %远休止角ts1=pi*4/9; %近休止角e=20; %偏距s0=125;rr=10; %滚子半径% 升程阶段t=linspace(0,pi*5/6,1000);s=h*(t/t0-sin(2*pi*t/t0)/(2*pi));x1=(s0+s).*cos(t)-e*sin(t);y1=(s0+s).*sin(t)+e*cos(t);dx1 =(h/t0-h*cos(2*pi*t/t0)).*cos(t)-(s0+s).*sin(t)- e*cos(t); dy1=(h/t0-h*cos(2*pi*t/t0)).*sin(t)+(s0+s).*cos(t)- e*sin(t);X1=x1-rr*dy1./(sqrt(dx1.^2+dy1.^2));Y1=y1+rr*dx1./(sqrt(dx1.^2+dy1.^2));hold onplot(x1,y1);plot(X1,Y1);% 远休止阶段t=linspace(pi*5/6,pi,1000);s=h;x2=(s+s0).*cos(t)-e*sin(t);y2=(s+s0).*sin(t)+e*cos(t);dx2 =- sin(t).*(s + s0) - e*cos(t);dy2 =cos(t).*(s + s0) - e*sin(t);X2=x2-rr*dy2./(sqrt(dx2.^2+dy2.^2));Y2=y2+rr*dx2./(sqrt(dx2.^2+dy2.^2));hold onplot(x2,y2);plot(X2,Y2);% 回程阶段t=linspace(pi,pi*14/9,1000);s=0.5*h*(1+cos(pi*(t-(t0+ts))/t01));x3=(s+s0).*cos(t)-e*sin(t);y3=(s+s0).*sin(t)+e*cos(t);dx3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*cos(t)- sin(t).*(s + s0) - e*cos(t); dy3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*sin(t)+ cos(t).*(s + s0) - e*sin(t);X3=x3-rr*dy3./(sqrt(dx3.^2+dy3.^2));Y3=y3+rr*dx3./(sqrt(dx3.^2+dy3.^2));hold onplot(x3,y3);plot(X3,Y3);%近休止阶段t=linspace(pi*14/9,pi*2,1000);s=0;x4=(s+s0).*cos(t)-e*sin(t);y4=(s+s0).*sin(t)+e*cos(t);dx4 =- sin(t).*(s + s0) - e*cos(t);dy4 =cos(t).*(s + s0) - e*sin(t);X4=x4-rr*dy4./(sqrt(dx4.^2+dy4.^2));Y4=y4+rr*dx4./(sqrt(dx4.^2+dy4.^2));hold onplot(x4,y4);plot(X4,Y4);hold offgrid ontitle('凸轮实际轮廓线','FontSize',20);。