哈工大机械原理大作业凸轮机构第四题

合集下载

哈工大、机械原理大作业、凸轮机构设计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)连架杆与机架必有最短杆1）.杆件1为曲柄2）.在各杆长度不变的情况下，选取c杆做为机架就可以实现双摇杆机构试以作图法设计一偏置尖底推杆盘形凸轮的轮廓曲线。

已知凸轮以等角速度顺时针回转，正偏距e=10,基园半径r0=30mm.推杆运动规律为：凸轮转角δ=0~150时，推杆00.凸轮转角δ=180~300时推杆等速上升16mm;.凸轮转角δ=150~180时推杆远休；等加速回程16mm;.凸轮转角δ=300~360时推杆近休。

解：解题步骤1）首先绘制位移S与转角δ的关系曲线S-δ曲线。

2）根据S-δ曲线、凸轮基园半径和正偏距，绘制凸轮的轮廓曲线。

000000凸轮仅用了0度，90度，150度，180度，300度几个点绘制轮廓曲线，同学们绘制时英多用些点（一般取12个点，再勾画轮廓曲线）第二篇：机械原理_凸轮机构设计机械原理课程设计——凸轮机构设计（一）目录 (1)＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿（一）、题目及原始数据 (2)（二）、推杆运动规律及凸轮廓线方程 (3)（三）、（四）、（五）、（六）、（七）、（八）、计算程序方框图..........................5 计算源程序..............................6 程序计算结果及分析......................10 凸轮机构图..............................15 心得体会................................16 参考书. (16)（一）、题目及原始数据试用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计，凸轮以1rad/s的角速度沿逆时针方向转动。

要求：（1）、推程运动规律为等加速等减速运动，回程运动规律为五次多项式运动规律；（2）、打印出原始数据；（3）、打印出理论轮廓和实际轮廓的坐标值；（4）、打印出推程和回程的最大压力角，以及出现最大压力角时凸轮的相应转角；（5）、打印出凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角；（6）、打印出凸轮运动的位移；（7）、打印最后所确定的凸轮的基圆半径。

哈工大机械原理连杆和凸轮大作业24题

班级学号机械原理大作业说明书题目 1、连杆机构运动分析2、凸轮机构设计学生姓名1连杆机构运动分析1.设计题目：一、先建立如下坐标系：二、划分杆组如下，进行结构分析：该机构由I级杆组RR（如图1）、II级杆组RPR（如图2、3）和II级杆组RRP（如图4）组成。

（1）（2）（3）（4）三、运动分析数学模型：（1）同一构件上点的运动分析：如右图所示的原动件1,已知杆1的角速度=10/rad s ω，杆长1l =170mm,A y =0，A x =110mm 。

可求得下图中B 点的位置B x 、B y ，速度xB v 、yB v ，加速度xB a 、yB a 。

θcos 1l xB =，θsin 1l yB =θωυsin 1l xB -=，θωυcos 1l yB =，222B2==-cos =-BxB i d x a l x dt ωϕω2222==-sin =-B yB i B d y a l y dtωϕω。

（2）RPRII 级杆组的运动分析：a. 如右图所示是由2个回转副和1个移动副组成的II 级组。

已知两个外运动副C 、B 的位置（B x 、B y 、c x =110mm 、C y =0）、速度（xB υ，yB υ，xC υ=0，yC υ=0）和加速度（0,0,,==yC xC yB xB a a a a ）。

可确定下图中D 点的位置、速度和加速度。

确定构件3的角位移1ϕ、角速度1ω、角加速度1α。

1sin 31..ϕϕl x dt dx C B -= 1s i n 131c o s 13.....2ϕϕϕϕl l x dt x d C B --= 1cos 31..ϕϕl y dt dy C B += 1c o s 131s i n 13.....2ϕϕϕϕl l y dty d C B +-= 根据关系：1111d 122..11.αϕϕωϕϕ====dtd dt ，故可得出： D x =)1cos(4βϕ++l x CD y =)1sin(4βϕ++l y Cb. 如右图所示是由2个回转副和1个移动副组成的II级组。

哈工大机械原理大作业_凸轮机构运动分析12题

Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮结构设计设计时间：2013年6月25日哈尔滨工业大学1、设计题目如图所示直动从动件盘形凸轮机构，原始参数见下表，据此设计该凸轮机构。

2、推杆升程、回程运动方程对于不同运动规律的凸轮结构，其上升与下降的方式不一，但遵循同样的运动顺序：上升、远休止点恒定，下降、近休止点恒定。

因此，在设计时，仅需确定这四个阶段的角度与位置即可。

（1）推杆升程运动方程式中（2）推杆远休程运动方程在远休程段，即时，。

（3）推杆回程运动方程式中（4）推杆近休程运动方程在远休程段，即时，。

3、推杆位移、速度、加速度线图(为方便作图和坐标的度量，取,用Matlab 作图)（1）推杆位移线图升程（mm ）升程运动角（°）升程运动规律升程许用压力角（°）回程运动角（°）回程运动规律回程许用压力角（°）远休止角（°）近休止角（°） 801503-4-5 多项式401003-4-5 多项式604070程序：x1=linspace(0,5*pi/6,1000);x2=linspace(5*pi/6,19*pi/18,1000);x3=linspace(19*pi/18,29*pi/18,1000);x4=linspace(29*pi/18,2*pi,1000);T1=x1/(5*pi/6);T2=(x3-19*pi/18)/(5*pi/9);s1=80*(10*T1.^3-15*T1.^4+6*T1.^5);s2=80;s3=80*(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（2）推杆速度线图程序：x1=linspace(0,5*pi/6,1000);x2=linspace(5*pi/6,19*pi/18,1000);x3=linspace(19*pi/18,29*pi/18,1000);x4=linspace(29*pi/18,2*pi,1000);T1=x1/(5*pi/6);T2=(x3-19*pi/18)/(5*pi/9);v1=30*80*1*T1.^2.*(1-2*T1+T1.^2)/(5*pi/6);v2=0;v3=-30*80*1*T2.^2.*(1-2*T2+T2.^2)/(5*pi/9);v4=0;plot(x1,v1,'k',x2,v2,'k',x3,v3,'k',x4,v4,'k') xlabel('角度ψ/rad');ylabel('速度v/(mm/s)')title('推杆速度线图')grid（3）推杆加速度线图程序：x1=linspace(0,5*pi/6,1000);x2=linspace(5*pi/6,19*pi/18,1000);x3=linspace(19*pi/18,29*pi/18,1000);x4=linspace(29*pi/18,2*pi,1000);T1=x1/(5*pi/6);T2=(x3-19*pi/18)/(5*pi/9);a1=60*80*1*1*T1.*(1-3*T1+2*T1.^2)/((5*pi/6)^2);a2=0;a3=-60*80*1*1*T2.*(1-3*T2+2*T2.^2)/((5*pi/9)^2);a4=0;plot(x1,a1,'k',x2,a2,'k',x3,a3,'k',x4,a4,'k')xlabel('角度ψ/rad');ylabel('加速度a/ ')title('推杆加速度线图')grid4、凸轮机构的线图程序：x1=linspace(0,5*pi/6,1000);x2=linspace(5*pi/6,19*pi/18,1000);x3=linspace(19*pi/18,29*pi/18,1000);x4=linspace(29*pi/18,2*pi,1000);T1=x1/(5*pi/6);T2=(x3-19*pi/18)/(5*pi/9);s1=80*(10*T1.^3-15*T1.^4+6*T1.^5);s2=80;s3=80*(1-(10*T2.^3-15*T2.^4+6*T2.^5));s4=0;v1=30*80*1*T1.^2.*(1-2*T1+T1.^2)/(5*pi/6);v2=0;v3=-30*80*1*T2.^2.*(1-2*T2+T2.^2)/(5*pi/9);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曲线')grid5、凸轮的基圆半径和偏距以凸轮机构的ds/dφ-s线图为基础，可分别作出三条限制线（推程许用压力角的切界线,回程许用压力角的限制线，起始点压力角许用线)，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。

哈工大机械原理大作业凸轮机构第四题教学教材

哈工大机械原理大作业凸轮机构第四题Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计姓名：李清蔚学号： 1140810304班级： 1408103指导教师：林琳一.设计题目设计直动从动件盘形凸轮机构，其原始参数见表1 表一：凸轮机构原始参数升程（mm ) 升程运动角（º）升程运动规律升程许用压力角（º）回程运动角（º）回程运动规律回程许用压力角（º）远休止角（º）近休止角（º）40 90 等加等减速30 50 4-5-6-7多项式60 100 120二.凸轮推杆运动规律（1）推程运动规律（等加速等减速运动）推程F0=90°①位移方程如下：②速度方程如下：③加速度方程如下：（2）回程运动规律（4-5-6-7多项式）回程00240190≤≤ϕ，F0=90°，F s=100°，F0’=50°其中回程过程的位移方程，速度方程，加速度方程如下：三.运动线图及凸轮s d ds -φ线图本题目采用Matlab 编程，写出凸轮每一段的运动方程，运用Matlab 模拟将凸轮的运动曲线以及凸轮形状表现出来。

代码见报告的结尾。

1、程序流程框图开始输入凸轮推程输入凸轮基圆偏距等基本参数输出ds,dv,da 图像输出压力角、曲率半径输出凸轮的构件形状结束2、运动规律ds图像如下：速度规律dv图像如下：加速度da规律如下图：3.凸轮的基圆半径和偏距以ds/dfψ-s图为基础，可分别作出三条限制线（推程许用压力角的切界限D t d t,回程许用压力角的限制线D t'd t'，起始点压力角许用线B0d'')，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。

得图如下：得最小基圆对应的坐标位置O点坐标大约为（13，-50）经计算取偏距e=13mm，r0=51.67mm.2.绘制理论轮廓线上的压力角曲线和曲率半径曲线针对凸轮转向及推杆偏置，令N1=1凸轮逆时针转；N2=1偏距为正。

哈工大机械原理大作业凸轮机构设计

Harbin Institute of Technology大作业设计说明书课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院班级：设计者：学号：指导教师：设计时间：2014年6月6日哈尔滨工业大学一、设计题目(1) 凸轮机构运动简图：二、凸轮推杆升程、回程运动方程及推杆位移,速度,加速度线图（1）推杆升程，回程运动方程如下：A.推杆升程方程：02φϕ≤≤时： 202s h ϕφ⎛⎫= ⎪⎝⎭124h v ωϕφ=2124h a ωφ=02φϕφ≤≤时：()2202hs h φϕφ=--()10204h v ωφϕφ=-21204h a ωφ=-式中 0150,2h πφ==,1ω=1rad/sB.推杆回程方程：'000s s φφϕφφφ+≤≤++()0'01cos 2s h s πϕφφφ⎧⎫=+-+⎡⎤⎨⎬⎣⎦⎭⎩ ()10''00sin 2s h v πωπϕφφφφ=--+⎡⎤⎣⎦ ()2210'2'00cos 2s h a πωπϕφφφφ=--+⎡⎤⎣⎦ 式中：'00150,9080,40s h φφφ====，（2）推杆位移,速度,加速度线图如下（凸轮转角从零到2π变化）：三、凸轮机构的-s 线图，并依次确定凸轮的基圆半径和偏距1)凸轮机构的-s 线图：2)确定凸轮的基圆半径和偏距：在两切线交点右下方为凸轮轴心的许用区域。

将轴心设在点（88,-49.03）处。

于是有偏距e=88mm, S0=49.03,基圆半径=100.72mm。

四、滚子半径的确定及凸轮理论轮廓和实际轮廓的绘制.1）确定凸轮半径可取凸轮半径r=11mm 2）确定凸轮理论廓线及实际轮廓五、计算编程h=150;%推杆上升的最大行程hw1=1;q=90*pi/180;%这里规定推程运动角为90度qs=40*pi/180;%远休止角为40度q1=80*pi/180;%回程运动角为80度for i=1:45%将45度按1度均分，从而得到各个度数上的轮廓坐标qq(i)=i*pi/180;s(i)=2*h*(qq(i)/q)^2;v(i)=4*h*w1*qq(i)/q^2;a(i)=4*h*w1^2/q^2;t(i)=4*h*(qq(i)/q^2);endfor i=45:90qq(i)=i*pi/180;s(i)=h-2*h/q^2*(q-qq(i))^2;v(i)=4*h*w1/q^2*(q-qq(i));a(i)=-4*h*w1^2/q^2;t(i)=4*h/q^2*q-4*h/q^2*qq(i);endfor i=90:130qq(i)=i*pi/180;s(i)=h;v(i)=0;a(i)=0;t(i)=0;endfor i=130:210qq(i)=i*pi/180;s(i)=h/2+h*cos((pi/q1)*(qq(i)-(q+qs)))/2;v(i)=-pi*h*w1*sin((pi/q1)*(qq(i)-(q+qs)))/(2*q1);a(i)=-pi*pi*h*w1*w1*cos((pi/q1)*(qq(i)-(q+qs)))/(2*q1*q1);t(i)=h*cos(pi/q1)/2;endfor i=210:360qq(i)=i*pi/180;s(i)=0;v(i)=0;a(i)=0;t(i)=0;end%推杆运动规律曲线绘制plot(qq,s)plot(qq,v)plot(qq,a)plot(t,s)% 绘制ds/dqq---曲线图hold on;f=@(x)x*(tan(50*pi/180))-154;x=0:195;y=f(x);plot(x,y)hold on;f=@(x)x*(tan(160*pi/180))-17;x=-70:130;y=f(x);plot(x,y)%计算曲率半径sym s,t,i;d=diff(t,s);dd=diff(t,s);R=abs((1+d^2)^(3/2)/dd);%绘制凸轮轮廓线r0=100.72;e=88;S0=49.03;Rr=11;x(i)=(S0-s(i))*cos(qq(i))-e*sin(qq(i));y(i)=(S0+s(i))*sin(qq(i))+e*cos(qq(i));d1(i)=-(S0-s(i))*sin(qq(i))-e*cos(qq(i)); d2(i)= (S0+s(i))*cos(qq(i))-e*sin(qq(i)); X=x-Rr*d1(i)/((d1(i))^2+(d2(i))^2)^(1/2); Y=y+Rr*d2(i)/((d1(i))^2+(d2(i))^2)^(1/2); plot(X,Y,x,y)。

哈工大机械原理大作业凸轮设计

哈工大机械原理大作业凸轮设计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 )为凸轮轴心位置。

哈工大机械原理大作业凸轮机构设计题

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机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计一．设计题目设计直动从动件盘形凸轮机构，1.运动规律（等加速等减速运动）推程 0450≤≤ϕ 推程 009045≤≤ϕ2.运动规律（等加速等减速运动）回程 00200160≤≤ϕ 回程 00240200≤≤ϕ三．推杆位移、速度、加速度线图及凸轮s d ds-φ线图采用VB 编程，其源程序及图像如下： 1.位移：Private Sub Command1_Click() Timer1.Enabled = True '开启计时器 End SubPrivate Sub Timer1_Timer() Static i As SingleDim s As Single, q As Single 'i作为静态变量，控制流程；s代表位移；q代表角度Picture1.CurrentX = 0Picture1.CurrentY = 0i = i + 0.1If i <= 45 Thenq = is = 240 * (q / 90) ^ 2Picture1.PSet Step(q, -s), vbRedElseIf i >= 45 And i <= 90 Thenq = is = 120 - 240 * ((90 - q) ^ 2) / (90 ^ 2)Picture1.PSet Step(q, -s), vbGreenElseIf i >= 90 And i <= 150 Thenq = is = 120Picture1.PSet Step(q, -s), vbBlackElseIf i >= 150 And i <= 190 Thenq = is = 120 - 240 * (q - 150) ^ 2 / 6400Picture1.PSet Step(q, -s), vbBlueElseIf i >= 190 And i <= 230 Thens = 240 * (230 - q) ^ 2 / 6400Picture1.PSet Step(q, -s), vbRedElseIf i >= 230 And i <= 360 Thenq = is = 0Picture1.PSet Step(q, -s), vbBlackElseEnd IfEnd Sub2.速度Private Sub Command2_Click()Timer2.Enabled = True '开启计时器End SubPrivate Sub Timer2_Timer()Static i As SingleDim v As Single, q As Single, w As Single 'i为静态变量，控制流程；q代表角度；w代表角速度，此处被赋予50Picture1.CurrentX = 0Picture1.CurrentY = 0w = 50i = i + 0.1If i <= 45 Thenv = 480 * w * q / 8100Picture1.PSet Step(q, -v), vbRedElseIf i >= 45 And i <= 90 Thenq = iv = 480 * w * (90 - q) / 8100Picture1.PSet Step(q, -v), vbBlack ElseIf i >= 90 And i <= 150 Thenq = iv = 0Picture1.PSet Step(q, -v), vbGreen ElseIf i >= 150 And i <= 190 Then q = iv = -480 * w * (q - 150) / 6400 Picture1.PSet Step(q, -v), vbBlue ElseIf i >= 190 And i <= 230 Then q = iv = -480 * w * (230 - q) / 6400Picture1.PSet Step(q, -v), vbRed ElseIf i >= 230 And i <= 360 Then q = iv = 0ElseEnd IfEnd Sub3.加速度Private Sub Command3_Click()Timer3.Enabled = True '开启计时器End SubPrivate Sub Timer3_Timer()Static i As SingleDim a As Single, w As Single, q As Single 'i为静态变量，控制流程；a代表加速度；q代表角度；w代表角速度w = 50Picture1.CurrentX = 0Picture1.CurrentY = 0i = i + 0.1If i <= 45 Thenq = ia = 480 * w ^ 2 / 8100Picture1.PSet Step(q, -a), vbRedElseIf i >= 45 And i <= 90 Thenq = ia = -480 * w ^ 2 / 8100ElseIf i >= 90 And i <= 150 Thenq = ia = 0Picture1.PSet Step(q, -a), vbGreenElseIf i >= 150 And i <= 190 Thenq = ia = -480 * w ^ 2 / 6400Picture1.PSet Step(q, -a), vbBlueElseIf i >= 190 And i <= 230 Thenq = ia = 480 * w ^ 2 / 6400Picture1.PSet Step(q, -a), vbRedElseIf i >= 230 And i <= 360 Thenq = ia = 0Picture1.PSet Step(q, -a), vbBlackElseEnd IfEnd Sub4.ds/dq---dsPrivate Sub Command4_Click()Timer4.Enabled = True '开启计时器；建立坐标系Picture1.Scale (-400, -400)-(400, 400)End SubPrivate Sub Timer4_Timer()Static i As SingleDim x As Single, s As Single, q As Single, scaley As Single, t As Single 'i为静态变量，控制流程；x代表位移；s代表纵坐标ds/dq；q代表角度Picture1.CurrentX = 0Picture1.CurrentY = 0scaley = 1t = 3.14 / 180i = i + 0.1If i <= 45 Thenq = i * tx = 194.734 * qs = 240 * (2 * q / 3.14) ^ 2Picture1.PSet Step(x, -s), vbRedElseIf i >= 45 And i <= 90 Thenq = i * tx = 194.734 * (3.14 / 2 - q)s = 120 - 97.367 * (3.14 / 2 - q) ^ 2Picture1.PSet Step(x, -s), vbRedElseIf i >= 90 And i <= 150 Thenq = i * tx = 0s = 120 * scaleyPicture1.PSet Step(x, -s), vbRedElseIf i >= 150 And i <= 190 Thenq = i * tx = -246.46 * (q - 5 * 3.14 / 6)s = 120 - 123.23 * (q - 5 * 3.14 / 6) ^ 2 Picture1.PSet Step(x, -s), vbRedElseIf i >= 190 And i <= 230 Thenq = i * tx = -246.46 * (23 * 3.14 / 18 - q)s = 123.23 * (23 * 3.14 / 18 - q) ^ 2Picture1.PSet Step(x, -s), vbRedElseIf i >= 230 And i <= 360 Thenq = i * tx = 0s = 0Picture1.PSet Step(x, -s), vbRedElseEnd IfEnd Sub四．确定凸轮基圆半径和偏距1. 求切点转角在图中，右侧曲线为升程阶段的类速度-位移图，作直线Dt dt与其相切，且位移轴正方向呈夹角[ 1]=300,则切点处的斜率与直线D t d t的斜率相等，因为k Dtdt=tan600,右侧曲线斜率可以表示为：q;q=tan600继而求出切点坐标（337.272,292.084）。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计姓名：李清蔚学号：1140810304 班级：1408103指导教师：林琳一.设计题目设计直动从动件盘形凸轮机构，其原始参数见表1表一：凸轮机构原始参数升程（mm )升程运动角（º）升程运动规律升程许用压力角（º）回程运动角（º）回程运动规律回程许用压力角（º）远休止角（º）近休止角（º） 4090等加等减速30 50 4-5-6-7多项式60 100 120二.凸轮推杆运动规律（1）推程运动规律（等加速等减速运动）推程F0=90°①位移方程如下：②速度方程如下：③加速度方程如下：（2）回程运动规律（4-5-6-7多项式）回程00240190≤≤ϕ，F0=90°，F s=100°，F0’=50°其中回程过程的位移方程，速度方程，加速度方程如下：三.运动线图及凸轮s d ds -φ线图本题目采用Matlab 编程，写出凸轮每一段的运动方程，运用Matlab 模拟将凸轮的运动曲线以及凸轮形状表现出来。

代码见报告的结尾。

1、程序流程框图开始输入凸轮推程回程的运动方程输入凸轮基圆偏距等基本参数输出ds,dv,da 图像输出压力角、曲率半径图像输出凸轮的构件形状结束2、运动规律ds图像如下：速度规律dv图像如下：加速度da规律如下图：3.凸轮的基圆半径和偏距以ds/dfψ-s图为基础，可分别作出三条限制线（推程许用压力角的切界限D t d t,回程许用压力角的限制线D t'd t'，起始点压力角许用线B0d'')，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。

压力角数学模型：曲率半径数学模型：)/)(/()/)(/(])/()/[(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+-=++=其中R0为确定的滚子半径，令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),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)-deg2 rad(100))/deg2rad(50)).^6+20*((phi4-deg2rad(90)-deg2rad(100))/deg2rad (50)).^7);plot(phi4,s4)hold onphi5=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*(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-deg2rad(90)-deg2 rad(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*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))/d eg2rad(50)).^2-1680*((phi4-deg2rad(90)-deg2rad(100))/deg2rad(50)).^3+ 2100*((phi4-deg2rad(90)-deg2rad(100))/deg2rad(50)).^4-840*((phi4-deg2 rad(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 onplot(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),1000);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)-deg 2rad(100))/deg2rad(50)).^6+20*((phi4-deg2rad(90)-deg2rad(100))/deg2ra d(50)).^7);phi5=linspace(deg2rad(240),deg2rad(360),1000)s05=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-de g2rad(90)-deg2rad(100))/deg2rad(50)).^5-140*((phi4-deg2rad(90)-deg2ra d(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(line([0,x3],[0,y3]))plot(line([x3,x3],[0,y3]))plot(line([0,x3],[y3,y3]))plot(x1,d1)hold onplot(x2,d2)hold ongrid onaxis([-110,60,-100,50])x4=13;y4=-50;plot(line([0,x4],[0,y4]))plot(line([x4,x4],[0,y4]),'k')plot(line([0,x4],[y4,y4]),'k')⑤压力角曲线，曲率半径曲线xo=13;yo=50;%题目2铰链O点取值e=xo;r0=sqrt(xo^2+yo^2);%基圆半径与偏距取值N1=1;N2=2;R0=15;%小滚子半径n=2;omiga=2*pi*n/60;w=omiga;phi1=linspace(0,90/2/180*pi,1000);phi2=linspace(90/2/180*pi,90/180*pi,1000);phi3=linspace(deg2rad(90),deg2rad(190),1000);phi4=linspace(deg2rad(190),deg2rad(240),1000);phi5=linspace(deg2rad(240),deg2rad(360),1000);phi0=90/180*pi;h=40;%推程s1=2*h*(phi1/phi0).^2+0*phi1;v1=4*h*w/(deg2rad(90)).^2*phi1+0*phi1;a1=4*h*w.^2/(phi0).^2+0*phi1;s2=h-2*h*(phi0-phi2).^2/(phi0).^2+phi2*0;v2=4*h*w*(phi0-phi2)/(deg2rad(90)).^2+phi2*0;a2=-4*h*w.^2/(phi0).^2+phi2*0;%远休程s3=h+phi3*0;v3=0+phi3*0;a3=0+phi3*0;%回程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)-deg2 rad(100))/deg2rad(50)).^6+20*((phi4-deg2rad(90)-deg2rad(100))/deg2rad (50)).^7)+phi4*0;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-deg2rad(90)-deg2 rad(100))/deg2rad(50)).^6]+phi4*0;a4=-(h*w.^2/(deg2rad(50)).^2)*[420*((phi4-deg2rad(90)-deg2rad(100))/d eg2rad(50)).^2-1680*((phi4-deg2rad(90)-deg2rad(100))/deg2rad(50)).^3+ 2100*((phi4-deg2rad(90)-deg2rad(100))/deg2rad(50)).^4-840*((phi4-deg2 rad(90)-deg2rad(100))/deg2rad(50)).^5]+phi4*0;s5=0+phi5*0;v5=0+phi5*0;a5=0+phi5*0;%全部运动过程s=[s1,s2,s3,s4,s5];v=[v1,v2,v3,v4,v5];a=[a1,a2,a3,a4,a5];phi=[phi1,phi2,phi3,phi4,phi5];%初始推程s0,起点坐标x0,y0%理论廓线s0=sqrt(r0^2-e^2);x0=(s0+s).*sin(N1*phi)+N2*e.*cos(N1*phi);y0=(s0+s).*cos(N1*phi)-N2*e.*sin(N1*phi);LXphi=(v./omiga-N1*N2*e).*sin(N1*phi.*pi/180)+N1*(s0+s).*cos(N1*phi.* pi/180);LYphi=(v./omiga-N1*N2*e).*cos(N1*phi.*pi/180)-N1*(s0+s).*sin(N1*phi.* pi/180);%计算压力角%推程压力角r1=s0+s1;ang1=abs(atan((v1/omiga-N1*N2*e)./r1))*180/pi;r2=s0+s2;ang2=abs(atan((v2/omiga-N1*N2*e)./r2))*180/pi;%回程压力角r4=s0+s4;ang4=abs(atan((v4/omiga-N1*N2*e)./r4))*180/pi;%远休程压力角r3=(s0+h)*ones(1,size(s3,5));ang3=abs(atan((-N1*N2*e)./r3))*180/pi+0*phi3;%近休程压力角r5=s0*ones(1,size(s5,5));ang5=abs(atan((-N1*N2*e)./r5))*180/pi+0*phi5;%凸轮的曲率半径LLXphi=(a./omiga^2.-s0-s).*sin(N1*phi.*pi/180)+(2*v./omiga-N1*N2*e).* N1.*cos(N1*phi.*pi/180);LLYphi=(a./omiga^2.-s0-s).*cos(N1*phi.*pi/180)-(2*v./omiga-N1*N2*e).* N1.*sin(N1*phi.*pi/180);A=(LXphi.^2+LYphi.^2).^1.5;B=-LXphi.*LLYphi+LYphi.*LLXphi;rou=A./B;M=-1;rou0=abs(rou+M*R0);plot(phi,rou0)axis([0,6,10,100])xlabel('凸轮转角/rad')ylabel('长度/mm')hold onangle=[ang1,ang2,ang3,ang4,ang5];plot(phi,angle)axis([0,6,0,100])4、理论轮廓线和实际轮廓线xo=13;yo=-50;%题目2铰链O点取值e=xo;r0=sqrt(xo^2+yo^2);%基圆半径与偏距取值N1=1;N2=1;R0=15;%小滚子半径n=2;omiga=2*pi*n/60;w=omiga;phi1=linspace(0,90/2/180*pi,1000);phi2=linspace(90/2/180*pi,90/180*pi,1000);phi3=linspace(deg2rad(90),deg2rad(190),1000);phi4=linspace(deg2rad(190),deg2rad(240),1000);phi5=linspace(deg2rad(240),deg2rad(360),1000);phi0=90/180*pi;h=40;%推程s1=2*h*(phi1/phi0).^2+0*phi1;v1=4*h*w/(deg2rad(90)).^2*phi1+0*phi1;a1=4*h*w.^2/(phi0).^2+0*phi1;s2=h-2*h*(phi0-phi2).^2/(phi0).^2+phi2*0;v2=4*h*w*(phi0-phi2)/(deg2rad(90)).^2+phi2*0;a2=-4*h*w.^2/(phi0).^2+phi2*0;%远休程s3=h+phi3*0;v3=0+phi3*0;a3=0+phi3*0;%回程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)-deg2 rad(100))/deg2rad(50)).^6+20*((phi4-deg2rad(90)-deg2rad(100))/deg2rad (50)).^7)+phi4*0;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-deg2rad(90)-deg2 rad(100))/deg2rad(50)).^6]+phi4*0;a4=-(h*w.^2/(deg2rad(50)).^2)*[420*((phi4-deg2rad(90)-deg2rad(100))/d eg2rad(50)).^2-1680*((phi4-deg2rad(90)-deg2rad(100))/deg2rad(50)).^3+ 2100*((phi4-deg2rad(90)-deg2rad(100))/deg2rad(50)).^4-840*((phi4-deg2 rad(90)-deg2rad(100))/deg2rad(50)).^5]+phi4*0;%近休程s5=0+phi5*0;v5=0+phi5*0;a5=0+phi5*0;%全部运动过程s=[s1,s2,s3,s4,s5];v=[v1,v2,v3,v4,v5];a=[a1,a2,a3,a4,a5];phi=[phi1,phi2,phi3,phi4,phi5];%初始推程s0,起点坐标x0,y0%理论廓线s0=sqrt(r0^2-e^2);x0=(s0+s).*sin(N1*phi)+N2*e.*cos(N1*phi);y0=(s0+s).*cos(N1*phi)-N2*e.*sin(N1*phi);LX=(v./omiga-N1*N2*e).*sin(N1*phi)+N1*(s0+s).*cos(N1*phi);LY=(v./omiga-N1*N2*e).*cos(N1*phi)-N1*(s0+s).*sin(N1*phi);plot(x0,y0,'R','linewidth',2)axis equalhold on%实际廓线M=-1;if R0==0Sx=x0;Sy=y0;elseH=sqrt(LX.^2+LY.^2);Sx=x0-N1*M*R0*LY./H;Sy=y0+N1*M*R0*LX./H;endplot(Sx,Sy)hold onaxis equal%画出基圆r1=r0;k=0:0.005:2*pi;x2=r1*sin(k);y2=r1*cos(k);plot(x2,y2,'-.')axis equalxlabel('凸轮转角/rad')ylabel('长度/mm')%画出偏距圆r2=e;k1=0:0.001:2*pi;x3=r2*sin(k1);y3=r2*cos(k1);plot(x3,y3,'k')r3=R0;k2=0:0.001:2*pi;x4=40.85+r3*sin(k2);y4=44.7+r3*cos(k2);plot(x4,y4,'k')r4=R0/5;k3=0:0.001:2*pi;x5=40.85+r4*sin(k3);y5=44.7+r4*cos(k3);plot(x5,y5,'k')hold ony6=44.7+50;plot(line([40.85,40.85],[44.7,y6]),'k') axis equalgrid onlegend('理论廓线','实际廓线','基圆','偏距圆');。