哈工大机械原理大作业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 υ===。
哈工大机械设计大作业二
哈工大机械设计大作业二————————————————————————————————作者：————————————————————————————————日期：大作业计算说明书题目:盘形凸轮轮廓的图解法设计学院:英才学院班号:班学号:6１2１820510姓名:林海奇日期：2014年9月27日哈尔滨工业大学大作业任务书题目:盘形凸轮轮廓的图解法设计设计原始数据：用图解法设计偏置滚子直动从动盘形凸轮轮廓。
原始数据如表格1(推杆的偏置方向及推杆推程和回程运动规律代号见表下方的注）。
表格错误!不能识别的开关参数。
设计原始数据凸轮角速度ω方向基圆半径r（ｍm)偏距e(mm)滚子半径rr(mｍ)推杆运动规律推程回程升程h（ｍm)推程角φ远休止角sφ回程角'φ近休止角'sφ逆时针50 12 １2错误!错误!4５1３°50°130°50°注：（１)推杆的偏置方向应使机构推程压力角较小。
(2）推杆的运动规律(推程、回程)①——等速运动规律;②——等加速度等减速运动规律；③——余弦加速度运动规律；——正弦加速度运动规律。
设计要求:1. 用A3图纸,按１:1比例绘图。
2. 凸轮理论轮廓线用点划线，实际轮廓线用粗实线。
3. 用虚线画出机架和从动件。
4. 作图过程中用到的线用细实线画。
5. 不校验压力角。
目录1. 设计过程…………………………………………………………………………………1 （1)取比例尺并作基圆(2)作反转运动，量取''00s s φφφφ、、、 ,并等分'00φφ、(３)计算推杆的预期位移 （4）确定理论轮廓线上的点 （5)绘制理论轮廓线 (６）绘制实际轮廓线2. 参考文献 (2)1.设计过程(1）取比例尺并作基圆,比例尺选为１:1，实际基圆半径为基圆半径0r 与滚子半径r r 之和，即６2ｍm 如图1所示。
(2)作反转运动，量取''00s s φφφφ、、、 ,并等分'00φφ、。
哈工大、机械原理大作业、凸轮机构设计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。
机械原理大作业——凸轮
大作业（二）凸轮机构设计题号： 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。
在终点处 h /2 s ,2/==δοδ。
整理得：⎪⎩⎪⎨⎧===^2^2/*h *4a ^2/**h *4v ^2^2/*h *2s δοωδοδωδοδ（ 注意：δ的变化范围为0~δ0/2。
机械原理大作业——凸轮.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。
哈工大机械原理大作业凸轮机构设计第题
哈工大机械原理大作业-凸轮机构设计(第题)————————————————————————————————作者：————————————————————————————————日期：机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电学院班级：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 绘制出位移、速度、加速度线图。
哈工大机械原理大作业直动从动件盘形凸轮机构满分完美版哈尔滨工业大学
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]);
机械原理大作业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 为确定的滚子半径。
- 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机械原理大作业说明书课程名称：机械原理设计题目：凸轮机构设计院系：能源科学与工程班级：1302403设计者：闫昭琦学号：1130240312指导教师：陈照波焦映厚设计时间：2015年6月哈尔滨工业大学一、设计题目 如右图所示直动从动件盘形凸轮机构，选择一组凸轮机构的原始参数，据此设计该凸轮机构。
二. 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图凸轮推杆升程运动方程：)]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);。