哈工大机械原理大作业二凸轮机构设计(29)
机械原理大作业二-凸轮机构设计..
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机械设计制造及其自动化班级：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 υ===。
哈工大、机械原理大作业、凸轮机构设计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。
机械原理大作业凸轮机构设计
机械原理大作业凸轮机构设计一、凸轮机构概述凸轮机构是一种常见的传动机构，它通过凸轮的旋转运动，带动相应零件做直线或曲线运动。
凸轮机构具有结构简单、运动平稳、传递力矩大等优点，在各种机械设备中得到广泛应用。
二、凸轮基本结构1. 凸轮凸轮是凸起的圆柱体，通常安装在主轴上。
其表面通常为圆弧形或其他曲线形状，以便实现所需的运动规律。
2. 跟随件跟随件是与凸轮配合的零件，它们通过接触面与凸轮相互作用，并沿着规定的路径做直线或曲线运动。
跟随件可以是滑块、滚子、摇臂等。
3. 连杆连杆连接跟随件和被驱动部件，将跟随件的运动转化为被驱动部件所需的运动。
连杆可以是直杆、摇杆等。
三、凸轮机构设计要点1. 几何参数设计设计时需要确定凸轮半径、角度和曲率半径等参数，这些参数的选择将直接影响凸轮机构的运动规律和性能。
2. 运动规律设计根据被驱动部件的运动要求，选择合适的凸轮曲线形状，以实现所需的运动规律。
3. 稳定性设计在设计凸轮机构时，需要考虑其稳定性。
例如，在高速旋转时，可能会发生跟随件脱离凸轮或者产生振动等问题，因此需要采取相应措施提高稳定性。
4. 材料和制造工艺设计在材料和制造工艺方面，需要考虑凸轮机构所承受的载荷和工作环境等因素，选择合适的材料和制造工艺。
四、几种常见凸轮机构及其应用1. 摇臂式凸轮机构摇臂式凸轮机构由摇臂、连杆和被驱动部件组成。
它通常用于实现直线运动或旋转运动，并且具有结构简单、运动平稳等优点。
摇臂式凸轮机构广泛应用于各种机械设备中，如发动机气门控制系统、纺织设备等。
2. 滑块式凸轮机构滑块式凸轮机构由凸轮、滑块、连杆和被驱动部件组成。
它通常用于实现直线运动，并且具有结构简单、运动平稳等优点。
滑块式凸轮机构广泛应用于各种机械设备中，如冲压设备、印刷设备等。
3. 滚子式凸轮机构滚子式凸轮机构由凸轮、滚子、连杆和被驱动部件组成。
它通常用于实现圆弧形运动，并且具有运动平稳、传递力矩大等优点。
滚子式凸轮机构广泛应用于各种机械设备中，如汽车发动机气门控制系统等。
哈工大机械原理大作业直动从动件盘形凸轮机构满分完美版哈尔滨工业大学
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，这时滚子与凸轮间接触应力最小，可提高凸轮寿命。
哈工大机械原理大作业凸轮设计
哈工大机械原理大作业凸轮设计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 )为凸轮轴心位置。
机械原理大作业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。
哈工大机械原理大作业二凸轮机构
void main()
{
float fi=0,fi0=90,h=100,S=0,v=0,a=0,fis=110,fi01=80,fis1=80,T2=0,w1=10,A=0,B=0,v0=0,v01=0,v02=0,f11=30,f31=70,A11=0,A12=0,B11=0,B12=0,A21=0,A22=0,B21=0,C11=0,C21=0,W=0;
{
W=fi*2*PI/360.0;
if(fi<=90)
{
S=h*(W/fi0-(1/(2*PI))*sin((2*PI/fi0)*W));
v=(h*w1/fi0)*(1-cos((2*PI/fi0)*W));
a=(2*PI*h*w1*w1/(fi0*fi0))*sin(2*PI*W/fi0);
}
elseif(fi>=90&&fi<200)
{
A=90;
A=A*2*PI/360.0;
S=h*(A/fi0-(1/(2*PI))*sin((2*PI/fi0)*A));
v=(h*w1/fi0)*(1-cos((2*PI/fi0)*A));
a=(2*PI*h*w1*w1/(fi0*fi0))*sin(2*PI*A/fi0);
}
elseif(fi>=200&&fi<=280)
}
printf("%f\n",a);
}
}
此程序计算推杆位移，速度，加速度线图。
程序二：
#include "stdio.h"
#include "math.h"
哈工大机械原理大作业凸轮 - 黄建青
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 凸轮机构原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角13 45mm 50°余弦加速度35°90°抛-直-抛70°100°120°计算流程框图：2. 凸轮推杆升程，回程运动方程及推杆位移、速度、加速度线图2.1 确定凸轮机构推杆升程、回程运动方程 设定角速度为ω=1 rad/s(1) 升程：0°<φ<50°由公式可得 )]cos(1[20ϕπΦh s -=计算推程、回程的推杆s 、v 、a运动线图及凸轮dss d ϕ-线图 确定凸轮机构基圆半径和偏距计算曲率半径和压力角，确定滚子半径确定凸轮的理论廓线和实际廓线)s i n (2001ϕπωπΦΦh v =)cos(2020212ϕπωπΦΦ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‘——回程运动角。
(完整word版)哈工大机械原理大作业凸轮DOC
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机械原理大作业二课程名称：机械原理设计题目: 凸轮机构设计院系：班级：设计者：学号：指导教师：哈尔滨工业大学一、设计题目如右图所示直动从动件盘形凸轮机构,选择一组凸轮机构的原始参数,据此设计该凸轮机构。
凸轮机构原始参数序号升程（mm）升程运动角升程运动规律升程许用压力角27130150正弦加速度30°回程运动角回程运动规律回程许用压力角远休止角近休止角100°余弦加速度60°30°80°二. 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图凸轮推杆升程运动方程:)]512sin(2156[130s ϕππϕ-= )512sin(4.374)]512cos(1[156v 211ϕπϕπωω=-=a％ t 表示转角,s 表示位移t=0：0.01:5*pi/6;%升程阶段s= [(6＊t)/(5*pi ）- 1/（2＊pi ）*sin(12＊t/5）］＊130; hold on plot(t ，s ）; t= 5＊pi/6:0。
01:pi; ％远休止阶段s=130； hold on plot(t,s ）;t=pi ：0.01:14*pi/9；%回程阶段s=65*［1+cos(9*（t-pi ）/5)]； hold on plot(t ，s ）;t=14*pi/9：0.01:2＊pi ；s=0；hold onplot（t,s）；grid onhold off％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。
哈工大机械原理大作业二凸轮机构设计（29）
哈工大机械原理大作业二凸轮机构设计（29）设计说明书1 设计题目如图所示直动从动件盘形凸轮机构，其原始参数见下表，据此设计该凸轮机构。
2、推杆升程、回程运动方程及位移、速度、加速度线图2.1凸轮运动理论分析推程运动方程：01cos 2h s π=-?? ?Φ?100sin 2h v πωπ=ΦΦ??221200cos 2h a πωπ=ΦΦ??回程运动方程： ()0'01s s h ?-Φ+Φ??=-??Φ??1'0h v ω=-Φ0a =2.2求位移、速度、加速度线图MATLAB 程序 pi=3.1415926; c=pi/180; h=140; f0=120; fs=45; f01=90; fs1=105; %升程 f=0:1:360; for n=0:f0s(n+1)=h/2*(1-cos(pi/f0*f(n+1)));v(n+1)=pi*h/(2*f0*c)*sin(pi/f0*f(n+1));a(n+1)=pi^2*h/(2*f0^2*c^2)*cos(pi/f0*f(n+1));end%远休程for n=f0:f0+fss(n+1)=140;v(n+1)=0;a(n+1)=0;end%回程for n=f0+fs:f0+fs+f01s(n+1)=h*(1-(f(n+1)-(f0+fs))/f01);v(n+1)=-h/(f01*c);a(n+1)=0;end%近休程for n=f0+fs+f01:360;s(n+1)=0;v(n+1)=0;a(n+1)=0;endfigure(1);plot(f,s,'k');xlabel('\phi/\circ');ylabel('s/mm');grid on;title('推杆位移线图') figure(2);plot(f,v,'k');xlabel('\phi/\circ');ylabel('v/(mm/s)');grid on;title('推杆速度线图') figure(3);plot(f,a,'k');xlabel('\phi/\circ');ylabel('a/(mm/s2');grid on;title('推杆加速度线图') 2.3位移、速度、加速度线图3 凸轮机构的dss d ?-线图，确定基圆半径和偏心距 3.1理论分析机构压力角α应按下式计算：ss ed ds +-=0/tan ?α以d s /d φ为横坐标，以s (φ)为纵坐标，可作出d s /d φ-s(φ)曲线如图所示，自D 点作∠BDd '=90?-[α]得直线Dd '，则在Dd '直线或其下方取凸轮轴心。
02-哈工大机械原理大作业任务书-凸轮机构
大作业2 凸轮机构设计
1.设计题目
如图2-1所示直动从动件盘形凸轮机构，其原始参数见表2-1。
从表2-1中选择一组凸轮机构的原始参数，据此设计该凸轮机构。
表2-1 凸轮机构原始参数
2.凸轮机构的设计要求
（1）确定凸轮推杆升程、回程运动方程，并绘制推杆位移、速度、加速度线图；
（2）绘制凸轮机构的s s -ϕ
d d 线图； （3）确定凸轮基圆半径和偏距；
（4）确定滚子半径；
图2-1
（5）绘制凸轮轮廓曲线。
3.要求提交的设计文件
（1）机构运动简图1张，凸轮轮廓图1张（标注凸轮的各部分参数），A3图纸；
（2）设计说明书1份，A4幅面。
4.设计说明书内容要求
设计说明书的封面按机械原理课程网站提供的格式制作。
设计说明书包含以下内容：
（1）设计题目，包括机构运动简图和凸轮机构的原始参数；
（2）凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图；
（3）凸轮机构的s s -ϕ
d d 线图，并依次确定凸轮的基圆半径和偏距； （4）滚子半径的确定及凸轮理论廓线和实际廓线的绘制。
5.凸轮机构的设计流程
（1）计算推杆的位移、速度和加速度线图；
（2）绘制凸轮机构的s s -ϕ
d d 线图； （3）确定凸轮的基圆半径和偏距；
（4）绘制凸轮理论廓线；
（5）确定滚子半径并绘制凸轮实际廓线。
机械原理大作业凸轮结构设计
机械原理大作业(二) 作业名称：机械原理设计题目：凸轮机构设计院系: 机电工程学院班级：设计者:学号:指导教师：丁刚陈明设计时间:哈尔滨工业大学机械设计1、设计题目如图所示直动从动件盘形凸轮机构，根据其原始参数设计该凸轮。
表一:凸轮机构原始参数序号升程(mm) 升程运动角(º)升程运动规律升程许用压力角（º)回程运动角（º)回程运动规律回程许用压力角(º)远休止角(º）近休止角（º)１2 ８0 1５０正弦加速度30 １０0 正弦加速度60 ６0 50２、凸轮推杆运动规律(1)推杆升程运动方程S=ｈ［φ/Φ0-siｎ（２πφ/Φ0)]V＝hω1／Φ０[１-ｃos(2πφ/Φ0）]a=2πhω１2sin(2πφ／Φ0)/Φ02式中:h=1５0，Φ0=5π／6,0<=φ<=Φ０,ω1=1(为方便计算)(2)推杆回程运动方程S=h[1-T/Φ１＋sin(2πＴ/Φ1)/２π]V= -hω１/Φ1[1-coｓ(2πT/Φ1)]ａ＝－２πhω12sin(２πT／Φ1)／Φ1２式中:h=150,Φ１=5π/9，７π/６<＝φ<=3１π／18,T=φ-7π/6３、运动线图及凸轮线图运动线图:用Ｍatlａｂ编程所得源程序如下:t=0:pi/５00:2*pｉ;ｗ1=1;h=150;lｅng=ｌeｎgｔh(ｔ);ｆor m=1:leng;if ｔ(m)＜=5*ｐi/6S(ｍ) = h*(t(m)/(5*pｉ/６)-sin（2*pi*t(m)/(5＊ｐi/6））/(2＊ｐi）);v(m)=h*w1*(1-cos（2*ｐi＊ｔ(m)／(5*pi/6)）)/(５*ｐi/6)；a(m)＝２*ｈ*w1＊w1*ｓin(２*pi*t（m)/(5＊pi/6)）/((5＊pｉ/6)*(5*pｉ/6))；% 求退程位移,速度,加速度elseｉｆt（m)＜=7*pｉ/6Ｓ(m）=h;v（m)=0;a(m)=０；% 求远休止位移，速度，加速度elｓｅif ｔ(ｍ)<=31*pi/１8T(m)＝t(m)-２1＊pi／18;S(m)=ｈ*(1-T(ｍ)／(5*pi/９)＋sｉｎ(2*pi＊T(m)/(5*pｉ／9）)/(2*pi）);v（ｍ)=－h/(5*pi/９）*(１-ｃos(2*pi*T（m)/(5*pi/９）))；a(m)=-２*pi＊h／(5*pi/9）＾2*sin(2＊pi＊Ｔ(m）/（5*pｉ/9));％求回程位移,速度，加速度elｓｅＳ（m)=0;v(m)=0;a(m)=０;% 求近休止位移，速度,加速度endend推杆位移图推杆速度图推杆加速度图4、确定凸轮基圆半径与偏距在凸轮机构得ds/dφ-s线图里再作斜直线Ｄt dｔ与升程得［ｄs/dφ-s(φ)］曲线相切并使与纵坐标夹角为升程许用压力角[α],则D t d t线得右下方为选择凸轮轴心得许用区。
02-机械原理大作业-凸轮机构
序号
升程 （mm） 100
升程运 动规律 正弦加速度
18
升程 许用 压力角 （ ） 40
回程运 动规律 等减等加速
2.确定凸轮推杆升程、回程运动方程
（1）推杆升程运功方程(0 < φ ≤ Φ0 = 150°) 推程为正弦加速度运动规律，因此，运动方程为： ������ 1 2������ s = h[ − sin ( ������)] Φ0 2������ Φ0 v= ℎ������ 2������ ( ������)] [1 − cos Φ0 Φ0 2������ℎ������2 Φ0
如图 2-1 所示直动从动件盘形凸轮机构，其原始参数见表 2-1。从表 2-1 中选择一组凸 轮机构的原始参数，据此设计该凸轮机构。
图 2-1 表 2-1 升程运 动角 （ ） 150 凸轮机构原始参数 回程运 动角 （ ） 100 回程 许用 压力角 （ ） 60 远休 止角 （ ） 70 近休 止角 （ ） 40
′2
Φ0 2
′
′
= 270°时
[������ − (Φ0 + Φ������ )]2
v=−
4ℎ������ Φ0
′2
[������ − (Φ0 + Φ������ )]
a=−
4ℎ������2 Φ0
′2
其中 h=100mm,Φ0 = 150°,Φ������ = 70°,Φ0 = 100°,������ = 10rad/s。 b)在270° = Φ0 + Φ������ +
3
凸轮机构设计 VB 编程
Call S2 Picture2.PSet (i, v), RGB(255, 0, 0) Next i For i = 270 To 320 Step 0.01 φ = i * pa Call S3 Picture2.PSet (i, v), RGB(255, 0, 0) Next i For i = 320 To 360 Step 0.01 φ = i * pa v = 0 Picture2.PSet (i, v), RGB(255, 0, 0) Next i End Sub Private Sub Command3_Click() '加速度 a 曲线 Picture3.Scale (-50, 15000)-(380, -15000) '定义坐标系 Picture3.Line (0, 0)-(380, 0) '画 x 轴 Picture3.Line (0, 15000)-(0, -15000) '画 y 轴 For i = 0 To 360 Step 30 '竖直网格线 Picture3.DrawStyle = 2 Picture3.Line (i, 15000)-(i, -15000) Picture3.CurrentX = i - 14: Picture3.CurrentY = -2 Picture3.Print i Picture3.CurrentX = 365: Picture3.CurrentY = 1800 Picture3.Print "φ" Next i For i = -14000 To 14000 Step 2000 '水平网格线 Picture3.DrawStyle = 2 Picture3.Line (0, i)-(370, i) Picture3.CurrentX = -50: Picture3.CurrentY = i + 500 Picture3.Print i Picture3.CurrentX = 10: Picture3.CurrentY = 15700 Picture3.Print "a" Next i Picture3.DrawWidth = 2 For i = 0 To 150 Step 0.01 φ = i * pa Call S1 Picture3.PSet (i, a), RGB(255, 0, 0) Next i
哈工大机械原理大作业2-凸轮
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮结构设计院系：能源学院指导教师：唐德威赵永强设计时间：2013年6月27日哈尔滨工业大学一、设计题目1、凸轮机构运动简图：2、凸轮机构的原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角21 110 150°3-4-5多项式40°100°3-4-5多项式60°45°65°二、凸轮推杆升程、回程运动方程及推杆位移,速度,加速度线图1、推杆升程，回程运动方程如下：A.推杆升程方程：设11/ rad sω=，由3-4-5多项式可知：当56ϕπ≤≤时，有：3451112211112211112(10156)30(12)60(132)s h T T Th Tv T Tha T T Tωφωφ=-+=-+=-+式中10/ϕφT=H=110 ,056φπ=B.推杆回程方程：当13π/12 ≤φ≤59π/36 时，有：2345222221222'0221222'0[1(10156)]30(12)60(132)s h T T T h v T T T h a T T T ωφωφ=--+=--+=--+ 式中 ()02'0s T ϕφφφ-+=h=110 056φπ= '059φπ= ϕs =5π/92、推杆位移,速度,加速度线图如下（用matlab 编程得）：A 、推杆位移线图clear clcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300); x3=linspace(13*pi/12,59*pi/36,300); x4=linspace(59*pi/36,2*pi,300); t1=x1/(5*pi/6)s1=110*(10*t1.^3-15*t1.^4+6*t1.^5); s2=110;t2=9*x3/(5*pi)-39/20;s3=110*(1-(10*t2.^3-15*t2.^4+6*t2.^5)); s4=0;plot(x1,s1,'k',x2,s2,'k',x3,s3,'k',x4,s4,'k') ; xlabel('角度/rad'); ylabel('位移s/mm'); title('推杆位移线图'); grid;B、推杆速度线图clearclcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300);x3=linspace(13*pi/12,59*pi/36,300);x4=linspace(59*pi/36,2*pi,300);f1=5*pi/6;t1=x1/f1;f2=5*pi/9;t2=9*x3/(5*pi)-39/20;v1=(t1.^2-2*t1.^3+t1.^4)*3300/f1;v2=0;v3=-30*110*(t2.^2-2*t2.^3+t2.^4)/f2;v4=0;plot(x1,v1,'k',x2,v2,'k',x3,v3,'k',x4,v4,'k') xlabel('角度/rad ');ylabel('速度v/(mm/s)');title('推杆速度线图');grid;C、推杆加速度线图clearclcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300);x3=linspace(13*pi/12,59*pi/36,300);x4=linspace(59*pi/36,2*pi,300);f1=5*pi/6;t1=x1/f1;f2=5*pi/9;t2=9*x3/(5*pi)-39/20;a1=60*110*(t1-3*t1.^2+2*t1.^3)/f1^2;a2=0;a3=-60*110*(t2-3*t2.^2+2*t2.^3)/f2^2;a4=0;plot(x1,a1,'k',x2,a2,'k',x3,a3,'k',x4,a4,'k') xlabel('角度/rad');ylabel('加速度a/ ');title('推杆加速度线图');grid;三、凸轮机构的ds/dψ---s线图，并依次确定凸轮的基圆半径和偏距1、凸轮机构的ds/dψ--s线图：clearclcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300);x3=linspace(13*pi/12,59*pi/36,300);x4=linspace(59*pi/36,2*pi,300);f2=5*pi/9;f1=5*pi/6;t1=x1/(5*pi/6)s1=110*(10*t1.^3-15*t1.^4+6*t1.^5);s2=110;t2=9*x3/(5*pi)-39/20;s3=110*(1-(10*t2.^3-15*t2.^4+6*t2.^5));s4=0;v1=(t1.^2-2*t1.^3+t1.^4)*3300/f1;v2=0;v3=-30*110*(t2.^2-2*t2.^3+t2.^4)/f2;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曲线');grid;2、确定凸轮的基圆半径和偏距：clearclcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300);x3=linspace(13*pi/12,59*pi/36,300);x4=linspace(59*pi/36,2*pi,300);f2=5*pi/9;f1=5*pi/6;t1=x1/(5*pi/6)s1=110*(10*t1.^3-15*t1.^4+6*t1.^5);s2=110;t2=9*x3/(5*pi)-39/20;s3=110*(1-(10*t2.^3-15*t2.^4+6*t2.^5));s4=0;v1=(t1.^2-2*t1.^3+t1.^4)*3300/f1;v2=0;v3=-30*110*(t2.^2-2*t2.^3+t2.^4)/f2;v4=0;k1=tan(pi/2-40*pi/180);k2=-tan(pi/6);f=sym('-k1*(2*k/f1^3-6*k^2/f1^4+4*k^3/f1^5)+k^2/f1^3-2*k^3/f1^4+k ^4/f1^5=0');k=solve(f);t01=k/f1;s01=110*(10*t01.^3-15*t01.^4+6*t01.^5);v01=(t01.^2-2*t01.^3+t01.^4)*3300/f1;c=80.5056;d=41.7790;%求出推程切点坐标x=-200:1:200;y5=k1*(x-c)+d;f2=5*pi/9;k2=-tan(pi/6);f=sym('-k2*(-2*(k*9/(5*pi)-39/20)*9/(5*pi)+6*(k*9/(5*pi)-39/20)^2 *9/(5*pi)-4*(k*9/(5*pi)-39/20)^3*9/(5*pi))-(k*9/(5*pi)-39/20)^2+2 *(k*9/(5*pi)-39/20)^3-(k*9/(5*pi)-39/20)^4=0');k=solve(f);t02=k*9/(5*pi)-39/20;s02=110*(1-(10*t02.^3-15*t02.^4+6*t02.^5));v02=-30*110*(t02.^2-2*t02.^3+t02.^4)/f2;o=32.1715;p= -112.4712;%求出回程切点坐标y6=k2*(x-p)+o;y7=x*-k1;plot(v1,s1,v2,s2,v3,s3,v4,s4,x,y5,x,y6,x,y7);xlabel('ds/dψ');ylabel('位移s/mm');title(' ds/dψ—s曲线');grid;所以，由图就可以确定回转中心所在的区域，所以，可取偏距 e=20mm ， 080s = mm ， 所以 222200802082.46r s e =+=+= mm 。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
设计说明书1 设计题目如图所示直动从动件盘形凸轮机构，其原始参数见下表，据此设计该凸轮机构。
2、推杆升程、回程运动方程及位移、速度、加速度线图2.1凸轮运动理论分析 推程运动方程： 01cos 2h s πϕ⎡⎤⎛⎫=-⎢⎥ ⎪Φ⎝⎭⎣⎦100sin 2h v πωπϕ⎛⎫=⎪ΦΦ⎝⎭221200cos 2h a πωπϕ⎛⎫= ⎪ΦΦ⎝⎭回程运动方程： ()0'1s s h ϕ-Φ+Φ⎡⎤=-⎢⎥Φ⎣⎦1'0hv ω=-Φ 0a =2.2求位移、速度、加速度线图MATLAB 程序 pi=3.1415926; c=pi/180; h=140; f0=120; fs=45; f01=90; fs1=105; %升程f=0:1:360; for n=0:f0s(n+1)=h/2*(1-cos(pi/f0*f(n+1)));v(n+1)=pi*h/(2*f0*c)*sin(pi/f0*f(n+1));a(n+1)=pi^2*h/(2*f0^2*c^2)*cos(pi/f0*f(n+1));end%远休程for n=f0:f0+fss(n+1)=140;v(n+1)=0;a(n+1)=0;end%回程for n=f0+fs:f0+fs+f01s(n+1)=h*(1-(f(n+1)-(f0+fs))/f01);v(n+1)=-h/(f01*c);a(n+1)=0;end%近休程for n=f0+fs+f01:360;s(n+1)=0;v(n+1)=0;a(n+1)=0;endfigure(1);plot(f,s,'k');xlabel('\phi/\circ');ylabel('s/mm');grid on;title('推杆位移线图')figure(2);plot(f,v,'k');xlabel('\phi/\circ');ylabel('v/(mm/s)');grid on;title('推杆速度线图')figure(3);plot(f,a,'k');xlabel('\phi/\circ');ylabel('a/(mm/s2');grid on;title('推杆加速度线图')2.3位移、速度、加速度线图3 凸轮机构的dss d ϕ-线图，确定基圆半径和偏心距 3.1理论分析机构压力角α应按下式计算：ss ed ds +-=0/tan ϕα以d s /d φ为横坐标，以s (φ)为纵坐标，可作出d s /d φ-s (φ)曲线如图所示，自D 点作∠BDd '=90˚-[α]得直线Dd '，则在Dd '直线或其下方取凸轮轴心。
同理可做斜直线D t d t 与升程的[d s /d φ-s(φ)]曲线相切并使与纵坐标夹角为升程[α]，则D t d t 线的右下方为选择凸轮轴心的许用区。
作斜直线D t 'd t '与回程的曲线相切，并使与纵坐标夹角为回程的[α](回程的[α]大于升程的[α])，则D t 'd t '线的左下方为选择凸轮轴心的许用区。
考虑到升程开始瞬时机构压力角也不超过许用值，自B 0点作限制线B 0d 0''与纵坐标夹角为升程[α]，则两直线D t d t 和B 0d 0''组成的d t O 1d 0'' 以下区域为选取凸轮中心的许用区，如选O 点作为凸轮回转中心，在推程和回程的任意瞬时，凸轮机构压力角均不会超过许用值，此时凸轮的基圆半径r 0=OB 0，偏距为e 。
3.2绘制dssdϕ-线图和轴心许用区域的MATLAB源程序syms t x%升程s1=h/2*(1-cos(pi/f0/c*t));ds1=diff(s1,'t',1);%远休程s2=h;%回程s3=h*(1-(t-(f0+fs)*c)/f01/c);ds3=diff(s3,'t',1);%近休程s4=0; %升程阶段相切时的角度t01=0.7560;%求出切点坐标x01=subs(ds1,t,t01);y01=subs(s1,t,t01);f1=tan(55*c)*(x-x01)+y01;f2=-tan(55*c)*x;f3=-tan(25*c)*(x+h/f01/c);figure(4)%每个阶段的所对应的转角ezplot(ds1,s1,[0,2*pi/3]);hold on;ezplot(ds3,s3,[165*c,255*c]);hold offgrid onaxis square;title('ds--s线图')figure(5)ezplot(ds1,s1,[0,120*c]);hold on;ezplot(ds3,s3,[165*c,255*c]);hold on;ezplot(f1,[-120,140]);hold on;ezplot(f3,[-120,140]);hold on;ezplot(f2,[-120,140]);hold offgrid ontitle('轴心的许用区域')%相交两直线的方程系数k1=tan(55*c);b1=-tan(55*c)*x01+y01;k2=-tan(25*c)b2=k2*h/f01/c;xo=(b2-b1)/(k1-k2);yo=k1*xo+b1;e=abs(xo) r0=sqrt(xo^2+yo^2)3.2 dssdϕ-线图和轴心许用区域求得：e =28.4666；r0 =61.7834；4 滚子半径的确定及凸轮理论轮廓和实际轮廓的绘制4.1理论分析（1）理论轮廓曲线：根据“反转法”原理，尖顶从动件在反转运动中的位置，即为凸轮理论轮廓对应点的位置。
当凸轮以ω1角速度顺时针转过φ角时，可将从动件从B0点推至B'(x'，y')点；按“反转法”原理，B'(x'，y')点沿-ω1方向(即逆时针方向)转过φ角到达B(x，y)点位置，即为凸轮轮廓上相应点的位置。
这样，凸轮理论轮廓方程可用图中矢量OB'沿逆时针方向旋转φ角到达矢量OB位置表示，故得⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡''yxRyxϕ式中，旋转矩阵⎥⎦⎤-⎢⎣⎡=ϕϕϕϕϕcossinsincosRB'点的坐标(x'，y')：⎥⎦⎤⎢⎣⎡+=⎥⎦⎤⎢⎣⎡essyx0''故得凸轮轮廓上点B(x，y)的坐标：⎥⎦⎤⎢⎣⎡+⎥⎦⎤-⎢⎣⎡=⎥⎦⎤⎢⎣⎡essyx0cossinsincosϕϕϕϕ将上式展开得凸轮理论轮廓方程：ϕϕϕϕcossin)(sincos)(essyessx++=-+=)20(πϕ≤≤式中220∧-∧=ers（2）工作轮廓：以理论轮廓上各点为圆心，以滚子半径为半径的圆族的包络线，即为滚子从动件凸轮的工作轮廓，或称实际轮廓。
即凸轮工作轮廓为单参数滚子圆族的平面曲线族的包络线。
凸轮理论轮廓上各点即滚子圆族圆心的坐标),(yx凸轮工作轮廓上点的坐标),(YX以凸轮转角φ为单参数的平面曲线族的包络线方程为⎭⎬⎫=∂∂=0/),,(0),,(ϕϕϕY X f Y X f当滚子圆半径为rr 时022)(2)(),,(=∧-∧-+∧-=r r y Y x X y x f ϕ0/)(2/)(2/),,(=----=∂∂ϕϕϕϕd dy y Y d dx x X Y X f 解该式得凸轮工作轮廓方程：)20(2)/(2)/(/2)/(2)/(/πϕϕϕϕϕϕϕ≤≤⎪⎪⎭⎪⎪⎬⎫∧+∧±=∧+∧=d dy d dx d dx r y Y d dy d dx d dy r x X rrμ4.2绘制理论轮廓MATLAB 源程序s0=sqrt(r0^2-e^2);x1=(s0+s1)*cos(t)-e*sin(t); y1=(s0+s1)*sin(t)+e*cos(t); figure(6)ezplot(x1,y1,[0,f0*c]) hold onx2=(s0+s2)*cos(t)-e*sin(t); y2=(s0+s2)*sin(t)+e*cos(t); ezplot(x2,y2,[f0*c,(f0+fs)*c]) hold onx3=(s0+s3)*cos(t)-e*sin(t); y3=(s0+s3)*sin(t)+e*cos(t);ezplot(x3,y3,[(f0+fs)*c,(f0+fs+f01)*c]) hold onx4=(s0+s4)*cos(t)-e*sin(t); y4=(s0+s4)*sin(t)+e*cos(t);ezplot(x4,y4,[(f0+fs+f01)*c,(f0+fs+f01+fs1)*c]) grid on4.3确定滚子半径的MATLAB源程序%求解滚子的半径 p表示曲率半径dx1=diff((s0+s1)*cos(t)-e*sin(t));dy1=diff((s0+s1)*sin(t)+e*cos(t));p1=sqrt(dx1^2+dy1^2);figure(7)ezplot(p1,[0,f0*c])hold ondx2=diff((s0+s2)*cos(t)-e*sin(t));dy2=diff((s0+s2)*sin(t)+e*cos(t));p2=sqrt(dx2^2+dy2^2);ezplot(p2,[f0*c,(f0+fs)*c])hold ondx3=diff((s0+s3)*cos(t)-e*sin(t));dy3=diff((s0+s3)*sin(t)+e*cos(t));p3=sqrt(dx3^2+dy3^2);ezplot(p3,[(f0+fs)*c,(f0+fs+f01)*c])hold ondx4=diff((s0+s4)*cos(t)-e*sin(t));dy4=diff((s0+s4)*sin(t)+e*cos(t));p4=sqrt(dx4^2+dy4^2);ezplot(p4,[(f0+fs+f01)*c,(f0+fs+f01+fs1)*c]) hold offgrid onaxis([0 7 0 250]);title('曲率半径曲线')滚子半径Rr=Pmin/2=284.4绘制工作轮廓MATLAB源程序Rr=28X1=x1-Rr*dy1/sqrt(dx1^2+dy1^2);Y1=y1+Rr*dx1/sqrt(dx1^2+dy1^2);figure(6)hold onezplot(X1,Y1,[0,f0*c])X2=x2-Rr*dy2/sqrt(dx2^2+dy2^2);Y2=y2+Rr*dx2/sqrt(dx2^2+dy2^2);hold onezplot(X2,Y2,[f0*c,163*c])X3=x3-Rr*dy3/sqrt(dx3^2+dy3^2);Y3=y3+Rr*dx3/sqrt(dx3^2+dy3^2);hold onezplot(X3,Y3,[166.5*c,(f0+fs+f01)*c])X4=x4-Rr*dy4/sqrt(dx4^2+dy4^2);Y4=y4+Rr*dx4/sqrt(dx4^2+dy4^2);hold onezplot(X4,Y4,[(f0+fs+f01)*c,(f0+fs+f01+fs1)*c]) xm=s0*cos(255*c)-e*sin(255*c);ym=s0*sin(255*c)+e*cos(255*c);Xm=xm+Rr*cos(t);Ym=ym+Rr*sin(t);hold onezplot(Xm,Ym,[103*c,141.5*c])hold offgrid ontitle('凸轮工作轮廓线')5 结果分析凸轮轮廓存在尖点，是由于回程运动为等速运动，推杆由远休程突变为等速运动时有速度的突变及加速度的剧烈变化，凸轮运动在该点附近失真，且易磨损。