哈工大机械原理大作业2-凸轮机构设计-22题()
机械原理大作业2凸轮(18题)
1、运动分析题目如图1所示直动从动件盘形凸轮机构，其原始参数见表1。
图 1表 1序号 升程 （mm ） 升程运动角（） 升程运动规律 升程许用压力角（） 回程运动角（）回程运动规律 回程许用压力角（） 远休止角（） 近休止角（）18 100150正弦加速度 30 100等减等加速 6040702、凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图：2.1从动件运动方程：（1）从动件升程运动方程升程段采用正弦加速度运动规律，运动方程为：()1212112100sin 5/6251001251cos 05/656210012sin 55/6s v a ϕϕππωπϕϕππωϕπ⎫⎡⎤⎛⎫=-⎪⎪⎢⎥⎝⎭⎣⎦⎪⎪⎡⎤⎪⎛⎫⎛⎫=-≤≤⎬ ⎪⎪⎢⎥⎝⎭⎝⎭⎣⎦⎪⎪⨯⎛⎫⎪= ⎪⎝⎭⎪⎭（2）从动件远休止运动方程在远休止s Φ段，即5/619/18πϕπ≤≤时，100s h mm ==,0v =,0a =。
（3）从动件回程运动方程升程段采用等减等加运动规律，运动方程为：()221221220019100518()94001919/184/3518()94005()9s v a πϕπωπϕπϕππωπ⎫⎪⎡⎤=--⎪⎢⎥⎣⎦⎪⎪⎪⎪⎡⎤=--≤≤⎬⎢⎥⎣⎦⎪⎪⎪⎪=-⎪⎪⎭()221221220029()518()94002()4/329/1853()94005()9s v a πϕπωπϕπϕππωπ⎫⎪=-⎪⎪⎪⎪⎪=--≤≤⎬⎪⎪⎪⎪=⎪⎪⎭（4）从动件近休止运动方程在近休止s 'Φ段，即29/182πϕπ≤≤时，0s =,0v =,0a =。
2.2推杆位移、速度、加速度线图：（1）推杆位移线图图 2 推杆位移线图（2）推杆速度线图图 3 推杆速度线图（3）推杆加速度线图图 4 推杆加速度线图3、凸轮机构的ds s d ϕ-线图，并由此确定凸轮的基圆半径和偏距：图 5 凸轮机构的dss d ϕ-线图 4滚子半径的确定及凸轮理论廓线和实际廓线的绘制 4.1凸轮的理论轮廓方程为：00()cos sin (02)()sin cos x s s e y s s e ϕϕϕπϕϕ=+-⎫≤≤⎬=++⎭式中，220031.45s r e mm =-=（1）推程凸轮轮廓方程：11231100sin cos 18sin 5/625(05/6)112(31100sin )sin 18cos 5/625x y ϕϕϕϕππϕπϕϕϕϕππ⎫⎧⎫⎡⎤⎛⎫=+--⎨⎬⎪⎪⎢⎥⎝⎭⎣⎦⎪⎩⎭≤≤⎬⎡⎤⎛⎫⎪=+-+ ⎪⎢⎥⎪⎝⎭⎣⎦⎭（2）远休止凸轮轮廓方程：131cos 18sin (5/619/18)131sin 18cos x y ϕϕπϕπϕϕ=-⎫≤≤⎬=+⎭（3）回程凸轮轮廓方程：222220019131cos 18sin 518()9(19/184/3)20019131sin 18cos 518()9x y πϕϕϕππϕππϕϕϕπ⎫⎧⎫⎪⎪⎪⎡⎤=---⎪⎨⎬⎢⎥⎣⎦⎪⎪⎪⎩⎭⎪≤≤⎬⎧⎫⎪⎪⎪⎡⎤⎪=--+⎨⎬⎢⎥⎪⎣⎦⎪⎪⎪⎩⎭⎭22222002931()cos 18sin 518()9(19/184/3)2002931()sin 18cos 518()9x y πϕϕϕππϕππϕϕϕπ⎫⎧⎫⎪⎪⎪=+--⎪⎨⎬⎪⎪⎪⎩⎭⎪≤≤⎬⎧⎫⎪⎪⎪⎪=+-+⎨⎬⎪⎪⎪⎪⎩⎭⎭（4）近休止凸轮轮廓方程：31cos 18sin (29/182)31sin 18cos x y ϕϕπϕπϕϕ=-⎫≤≤⎬=+⎭4.2凸轮理论轮廓曲线为：图 6 凸轮理论轮廓由上图可编程可求其最小曲率半径为min 10.309110mm ρ=≈，所以滚子半径min 1037r r mm ρ=-∆=-=。
机械原理大作业二-凸轮机构设计..
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机械设计制造及其自动化班级：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。
哈工大机械原理连杆和凸轮大作业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级组。
哈工大机械原理大作业-凸轮机构设计
哈工大机械原理大作业-凸轮机构设计(第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ϕϕ⨯=⨯通过编程求其角度。
哈工大—机械原理凸轮大作业
一、题目要求及机构运动简图如图1所示直动从动件盘形凸轮机构。
其原始参数见表1。
图一凸轮运动简图表一凸轮原始参数二、计算流程框图凸轮机构分析建立数学模型位移方程速度方程加速度方程速度线图位移线图加速线图ds/dΨ-s曲线升程压力角回程压力角确定轴向及基圆半径压力角图确定滚子半径实际轮廓理论轮廓轮廓图结束三、建立数学模型1.从动件运动规律方程首先，由于设计凸轮轮廓与凸轮角速度无关，所以不妨设凸轮运动角速度为w = 1rad/s。
（1）推程运动规律(0 < φ < 90°)s=φ）v=φφ）a=φφ）式中：h=65mm，Φ0=π/2（2）远休程运动规律(90°< φ < 190°)s = 65mmv = 0a = 0（3）回程运动规律(190°< φ < 240°)(190°< φ < 196.25°)(196.25°< φ < 233.75°)(233.75°< φ < 240°) 回程运动中的速度和加速度为位移对时间t的倒数：（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 线的右下方为选择凸轮轴心的许用区。
哈工大机械原理大作业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。
机械原理大作业2凸轮机构大作业..
大作业（二）凸轮机构设计（题号：_________）班级：________________________ 学号：________________________ 姓名：________________________ 同组其他人员：________________________ 完成日期：________________________凸轮机构大作业题目目录1、题目及原始数据；2、推杆的运动规律及凸轮廓线方程；3、计算程序框图4、计算程序；5、计算结果及分析；6、凸轮机构图（包括推杆及凸轮理论和实际廓线，并标出有关尺寸及计算结果7、体会及建议8、参考书利用计算机辅助设计完成下列偏置式直动滚子推杆盘形凸轮机构（推杆的移动副导路位于凸轮盘回转中心的右侧）或摆动滚子推杆盘形凸轮机构的设计，设计已知数据如下表所示，机构中凸轮沿着逆时针方向做匀速转动。
表1 两种凸轮机构的从动件运动规律表2 两种凸轮机构的推杆在近休、推程、远休及回程阶段的凸轮转角表3 偏置直动滚子推杆盘形凸轮机构的已知参数要求：每两人一组，每组中至少打印出一份源程序。
每人都要打印：原始数据；凸轮理论轮廓曲线和实际轮廓曲线的坐标值；推程和回程的最大压力角，以及出现最大压力角时凸轮相应的转角，凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角；凸轮的基圆半径。
整个设计过程所选取的计算点数N=72~120。
利用计算机绘出凸轮的理论轮廓曲线和实际轮廓曲线。
二、推杆运动规律及凸轮轮廓方程推程： 等加速2202/s h δδ=远休:s h =等减速22002()/s h h δδδ=--回程 : 五次多项式整理得；345121212s=h[1-10(()/)+15(()/)6(()/)]o o o δδδδδδδδδδδδ-------理论轮廓廓线方程00()sin cos ()cos sin x s s e y s s e δδδδ=++=+-工作廓线方程cos sin t t x x r y y r δδ=-=-实际廓线方程sin (/)/cos (/)/dx d dy d θδθδ==-'cos 'sin r r x x r y y r θθ=-=-三、计算程序框图四、计算程序#include<stdio.h>#include<math.h>#include<stdlib.h>#define PI 3.141592653double fact[72][2];double theory[72][2];int ang1=180,ang2=225,ang3=315; double h=35, rb=20,b=1;double A1=30*PI/180, A2=75*PI/180; double P=18.2,e=10;double So,r=15;double S(int I){double s;double A;double B;if(I<=ang1/2){A=I*PI/180;B=ang1*PI/180;s=2*h*pow(A/B,2);}else if((I>ang1/2)&&(I<=ang1)){A=I*PI/180;B=ang1*PI/180;s=h-2*h*pow((B-A)/B,2);}else if(I<=ang2)s=h;else if(I<=ang3){A=(I-ang2)*PI/180;B=(ang3-ang2)*PI/180;s=h*(1-10*pow(A/B,3)+15*pow(A/B,4)-6*pow(A/B,5));}else s=0;return(s);}double ds(int Q){double A,B,C;if(Q<=ang1/2){A=Q*PI/180;B=ang1*PI/180;C=4*h*A/(B*B);}else if((Q>ang1/2)&&(Q<=ang1)){A=Q*PI/180;B=ang1*PI/180;C=4*h*(B-A)/(B*B);}else if(Q<=ang2) C=0;else if(Q<=ang3){A=(Q-ang2)*PI/180;B=(ang3-ang2)*PI/180;C=h*(-30*A*A/pow(B,3)+60*pow(A,3)/pow(B,4)-30*pow(A,4)/pow(B,5));}else C=0;return C;}double dss(int B3){double A,B,C;if(B3<=ang1/2){A=B3*PI/180;C=ang1*PI/180;B=4*h/(C*C);}else if(B3>ang1/2&&B3<=ang1){A=B3*PI/180;C=ang1*PI/180;B=-4*h/(C*C);}else if(B3<=ang2)B=0;else if(B3<=ang3){A=(B3-ang2)*PI/180;C=(ang3-ang2)*PI/180;B=h*(-60*A/pow(C,3)+240*A*A/pow(C,4)-120*A*A*A/pow(C,5));}else B=0;return(B);}void xy(int ang){double A,B,C,E,F,dx,dy;A=ang*PI/180;B=S(ang);C=ds(ang);dx=(So+B)*cos(A)+sin(A)*C-e*sin(A);dy=-sin(A)*(So+B)+C*cos(A)-e*cos(A);E=r*dy/sqrt(dx*dx+dy*dy);F=r*dx/sqrt(dx*dx+dy*dy);theory[ang/5][0]=(So+B)*sin(A)+e*cos(A);theory[ang/5][1]=(So+B)*cos(A)-e*sin(A);fact[ang/5][0]=theory[ang/5][0]-E;fact[ang/5][1]=theory[ang/5][1]+F;}double a(int B1)/*****求解压力角****/{double A,B;A=sqrt((ds(B1)-e)*(ds(B1)-e));B=S(B1);return atan(A/(B+So));}double p(int B2){double dx,dy,dxx,dyy;double A,B,C,D,E;A=B2*PI/180;B=ds(B2);C=S(B2);D=dss(B2);dx=(So+C)*cos(A)+sin(A)*B-e*sin(A);dy=-sin(A)*(So+C)+B*cos(A)-e*cos(A);dxx=-(C+So)*sin(A)+cos(A)*B+D*sin(A)-e*cos(A);dyy=-cos(A)*(So+C)-B*sin(A)+D*cos(A)-sin(A)*B+e*sin(A);E=sqrt(pow(dx*dx+dy*dy,3))/sqrt(pow((dx*dyy-dxx*dy),2));return(E);}void main(){ FILE *fp;int i;int k,h,l;double angle1max=0,angle2max=0,pmin=1000;if((fp=fopen("sanying","w"))==NULL){printf("Cann't open this file.\n");exit(0);}fprintf(fp,"\n The Kinematic Parameters of Point 4\n");fprintf(fp," x y x' y' ");//计算数据并写入文件for(;i!=360;){rb=rb+b;So=sqrt(rb*rb-e*e);for(i=0;i<=ang1;i=i+5){if(a(i)>A1||p(i)<P)break;}if(ang1+5-i)continue;for(i=ang1+5;i<=ang2;i=i+5){if(p(i)<P)break;}if(ang2+5-i)continue;for(i=ang2+5;i<=ang3;i=i+5){if(a(i)>A2||p(i)<P)break;}if(ang3+5-i)continue;for(i=ang3+5;i<360;i=i+5){if(p(i)<P)break;}}for(i=0;i<360;i=i+5){xy(i);}for(i=0;i<=ang1;i=i+5){if(angle1max<a(i)){angle1max=a(i);k=i;}if(pmin>p(i)){pmin=p(i);h=i;}}for(i=ang2;i<=ang3;i=i+5){if(angle2max<a(i)){angle2max=a(i);l=i;}if(pmin>p(i)){pmin=p(i);h=i;}}for(i=0;i<72;i++){fprintf(fp,"\n");{fprintf(fp,"%12.3f\t%12.3f\t%12.3f\t%12.3f\t",theory[i][0],theory[i][1],fact[i][0],fact[i][1]);}}fclose(fp);printf(" 理论坐标(x,y) ");printf("实际坐标(x,y)");printf("\n");for(i=0;i<72;i++){printf("%f ",theory[i][0]);printf(" ");printf("%f ",theory[i][1]);printf(" ");printf("%f ",fact[i][0]);printf(" ");printf("%f ",fact[i][1]);printf("\n");}printf("基圆半径是:%f\n",rb);printf("推程最大压力角是:%f\n",angle1max*180/PI);printf("此时角度是是:%d\n",k);printf("回程最大压力角是:%f\n",angle2max*180/PI);printf("此时角度是是:%d\n",l);printf("最小曲率半径是:%f\n",pmin);printf("此时角度是:%d\n",h);}五、计算结果及分析-58.095654100.403660-65.608023113.386892 -49.123827105.084964-55.476046118.673537 -39.778138108.966507-44.921863123.057004 -30.129713112.018750-34.025797126.503933 -20.251983114.218463-22.870774128.988091 -10.220123115.548903-11.541691130.490572 -0.110482115.999947-0.124769130.999941分析基圆半径是:116.000000推程最大压力角是:5.273246 此时角度是是:90回程最大压力角是:21.408137 此时角度是是:275最小曲率半径是:18.384284 此时角度是:315-200-1000100200-200-1000100200理论坐标y 实际坐标y（六）心得体会在解决大作业过程中不仅仅让自己更熟悉课本知识同时使得自己重温C 语言让自己更加熟练与程序的设计，提高了自己的逻辑运用能力，这种对运行的机构的认识，我相信对以后的理论知识求解也有帮助。
哈工大机械原理大作业二凸轮机构
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。
哈工大机械原理大作业2-凸轮机构设计-22题
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) 升程运动角(º) 升程运 动规律升程许用压力角（º） 回程运动角（º）回程运 动规律回程许用压力角（º）远休止角 （º）近休止角 （º） 22 120 90等加等减速 4080等减等加速 70 70120二． 凸轮推杆运动规律1.运动规律（等加速等减速运动） 推程 0450≤≤ϕ2229602ϕπϕ=⎪⎪⎭⎫ ⎝⎛Φ=h sϕπωϕω2219204=Φ=h v2220219204πωω=Φ=h a推程 009045≤≤ϕ ()22220)2(9601202ϕππϕ--=-ΦΦ-=h h s())2(1920422ϕπωπϕω-=-ΦΦ=h v222219204ωπω-=Φ-=h a2.运动规律（等加速等减速运动） 回程 0200160≤≤ϕ()[]22202)98(9601202πϕπϕ--=Φ+Φ-Φ-=Sh h s()[])98(1920-4-202πϕωπϕω-=Φ+Φ-Φ=Sh v222219204ωπω-=Φ-=h a回程 0240200≤≤ϕ()[]222'002)34(9602ϕππϕ-=-Φ+Φ+ΦΦ=S h s()[])34(1920-4-2'002ϕπωπϕω-=-Φ+Φ+ΦΦ=S h v222219204ωπω=Φ=h a三．推杆位移、速度、加速度线图及凸轮sd 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 = 0 Picture1.CurrentY = 0 i = i + 0.1 If i <= 45 Then q = 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 Thenq = is = 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 Thenq = iv = 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) / 6400Picture1.PSet Step(q, -v), vbBlue ElseIf i >= 190 And i <= 230 Thenq = iv = -480 * w * (230 - q) / 6400Picture1.PSet Step(q, -v), vbRedElseIf i >= 230 And i <= 360 Then q = iv = 0Picture1.PSet Step(q, -v), vbBlack ElseEnd 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 / 8100Picture1.PSet Step(q, -a), vbBlackElseIf 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) ^ 2Picture1.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. 求切点转角在图中，右侧曲线为升程阶段的类速度-位移图，作直线D t d t与其相切，且位移轴正方向呈夹角[ 1]=300,则切点处的斜率与直线D t d t的斜率相等，因为k Dtdt=tan600,右侧曲线斜率可以表示为：q;q=tan600继而求出切点坐标（337.272,292.084）。
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
- 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机械原理大作业二课程名称： 机械原理 设计题目： 凸轮机构设计一．设计题目设计直动从动件盘形凸轮机构， 凸轮机构原始参数 序号 升程(mm) 升程运动角(º) 升程运 动规律升程许用压力角（º） 回程运动角（º） 回程运 动规律回程许用压力角（º） 远休止角 （º） 近休止角 （º） 2212090等加等减速4080等减等加速7070120二． 凸轮推杆运动规律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 Thenq = is = 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 Thenq = iv = 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) / 6400Picture1.PSet Step(q, -v), vbBlue ElseIf i >= 190 And i <= 230 Thenq = iv = -480 * w * (230 - q) / 6400Picture1.PSet Step(q, -v), vbRedElseIf i >= 230 And i <= 360 Then q = iv = 0Picture1.PSet Step(q, -v), vbBlack ElseEnd 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 / 8100Picture1.PSet Step(q, -a), vbBlackElseIf 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) ^ 2 Picture1.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）。
同理可求出另一切点（-89.766,16.347）2.确定直线方程y-292.084=1.732*（x-337.272）y-16.347=-0.364*（x+89.766）3.绘图确定基圆半径和偏距程序如下：Private Sub Command5_Click() '开启计时器；画出ds/dq-s图的切线Timer5.Enabled = TruePicture1.Scale (-400, -400)-(400, 400)Picture1.Line (-50, 378.671)-(200, -54.329), vbRedPicture1.Line (-200, -65.683)-(40, -481.363), vbBluePicture1.Line (-200, -209.127)-(200, -63.527), vbBlackPicture1.Line (-200, -56.472)-(200, 89.128), vbGreenEnd SubPrivate Sub Timer5_Timer() 'i为静态变量，控制流程；x代表位移；s代表纵坐标ds/dq；q代表角度；scaley是图形比例系数；t为角度，弧度转化系数Static i As SingleDim x As Single, s As Single, q As Single, scaley As Single, t As Single 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) ^ 2Picture1.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在轴心公共许用区内取轴心位置，能够满足压力角要求，现取直线红线与直线绿线的交点为轴=146.392 偏距e=131.557 s0=64.013心位置，通过解方程组：x=131.557 y=-64.214 r五．绘制凸轮理论廓线和实际廓线1.绘制凸轮理论轮廓Private Sub Command6_Click() '开启计时器；建立坐标系Picture1.Scale (-500, -500)-(500, 500)Timer6.Enabled = TrueEnd SubPrivate Sub Timer6_Timer() 'i为静态变量，控制流程；e为偏距；s为基圆半径平方减去偏距平方再开方所得；q为角度；x为横坐标；y为纵坐标Static i As SingleDim e As Single, q As Single, s As Single, x As Single, y As Single e = 131.557s = 64.213Picture1.CurrentX = 12.764Picture1.CurrentY = -95.125i = i + 0.001If i <= 3.14 / 4 Thenq = ix = (s + 240 * (2 * q / 3.14) ^ 2) * Cos(q) - e * Sin(q)y = (s + 240 * (2 * q / 3.14) ^ 2) * Sin(q) + e * Cos(q)Picture1.PSet Step(x, -y), vbRedElseIf i >= 3.14 / 4 And i <= 3.14 / 2 Thenq = ix = (s + 120 - 97.367 * (3.14 / 2 - q) ^ 2) * Cos(q) - e * Sin(q) y = (s + 120 - 97.367 * (3.14 / 2 - q) ^ 2) * Sin(q) + e * Cos(q) Picture1.PSet Step(x, -y), vbRedElseIf i >= 3.14 / 2 And i <= 5 * 3.14 / 6 Thenq = ix = (s + 120) * Cos(q) - e * Sin(q)y = (s + 120) * Sin(q) + e * Cos(q)Picture1.PSet Step(x, -y), vbRedElseIf i >= 5 * 3.14 / 6 And i <= 19 * 3.14 / 18 Thenq = ix = (s + 120 - 123.23 * (q - 5 / 6 * 3.14) ^ 2) * Cos(q) - e * Sin(q)y = (s + 120 - 123.23 * (q - 5 / 6 * 3.14) ^ 2) * Sin(q) + e * Cos(q)Picture1.PSet Step(x, -y), vbRedElseIf i >= 19 / 18 * 3.14 And i <= 23 / 18 * 3.14 Thenq = ix = (s + 123.23 * (23 / 18 * 3.14 - q) ^ 2) * Cos(q) - e * Sin(q)y = (s + 123.23 * (23 / 18 * 3.14 - q) ^ 2) * Sin(q) + e * Cos(q)Picture1.PSet Step(x, -y), vbRedElseIf i >= 23 / 18 * 3.14 And i <= 2 * 3.14 Thenq = ix = s * Cos(q) - e * Sin(q)y = s * Sin(q) + e * Cos(q)Picture1.PSet Step(x, -y), vbRedElseEnd IfEnd Sub2.确定滚子半径并绘制曲率半径图Private Sub Command7_Click() '开启计时器，建立坐标系Timer7.Enabled = TruePicture1.Scale (-10, -500)-(10, 10)End SubPrivate Sub Timer7_Timer() 'i为静态变量，控制流程；x代表dx/dq;y代表dy/dq；p 代表曲率；e为偏距；s为基圆半径平方减去偏距平方再开方所得Static i As SingleDim s As Single, q As Single, x As Single, y As Single, p As Singles = 64.213e = 131.557Picture1.CurrentX = 0Picture1.CurrentY = -Sqr(e ^ 2 + s ^ 2)i = i + 0.01If i <= 3.14 / 4 Thenq = ix = (480 * 4 * q / 3.14 ^ 2) * Cos(q) + (s + 240 * (2 * q / 3.14) ^ 2) * (-Sin(q)) - e * Cos(q)y = (480 * 4 * q / 3.14 ^ 2) * Sin(q) - e * Sin(q) + (s + 240 * (2 * q / 3.14) ^ 2) * Cos(q)p = Sqr(x ^ 2 + y ^ 2)Picture1.PSet Step(q, -p), vbRedElseIf i >= 3.14 / 4 And i <= 3.14 / 2 Thenq = ix = (2 * 97.367 * (3.14 / 2 - q)) * Cos(q) + (s + 120 - 97.367 * (3.14 / 2 - q) ^ 2) * (-Sin(q)) - e * Cos(q)y = (2 * 97.367 * (3.14 / 2 - q)) * Sin(q) + (s + 120 - 97.367 * (3.14 / 2 - q) ^ 2) * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Picture1.PSet Step(q, -p), vbRedElseIf i >= 3.14 / 2 And i <= 5 / 6 * 3.14 Thenq = ix = (s + 120) * (-Sin(q)) - e * Cos(q)y = (s + 120) * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Picture1.PSet Step(q, -p), vbRedElseIf i >= 5 / 6 * 3.14 And i <= 19 / 18 * 3.14 Thenq = ix = (-2 * 123.23 * (q - 5 / 6 * 3.14)) * Cos(q) + (s + 120 - 123.23 * (q - 5 / 6 * 3.14) ^ 2) * (-Sin(q)) - e * Cos(q)y = (-2 * 123.23 * (q - 5 / 6 * 3.14)) * Sin(q) + (s + 120 - 123.23 * (q - 5 / 6 * 3.14) ^ 2) * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Picture1.PSet Step(q, -p), vbRedElseIf i >= 19 / 18 * 3.14 And i <= 23 / 18 * 3.14 Thenq = ix = (-2 * 123.23 * (23 / 18 * 3.14 - q)) * Cos(q) + (s + 123.23 * (23 / 18 * 3.14 - q) ^ 2) * (-Sin(q)) - e * Cos(q)y = (-2 * 123.23 * (23 / 18 * 3.14 - q)) * Sin(q) + (s + 123.23 * (23 / 18 * 3.14 - q) ^ 2) * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Picture1.PSet Step(q, -p), vbRedElseIf i >= 23 / 18 * 3.14 And i <= 2 * 3.14 Thenq = ix = s * (-Sin(q)) - e * Cos(q)y = s * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Picture1.PSet Step(q, -p), vbRedElseEnd IfEnd Sub其中曲率最小值为126.026mm 可以得到滚子半径为63.013mm3.绘制实际轮廓线Private Sub Command8_Click()Picture1.Scale (-500, -500)-(500, 500) '建立坐标系；开启计时器Timer8.Enabled = TrueEnd SubPrivate Sub Timer8_Timer() 'i为静态变量，控制流程；q代表角度；x代表dx/dy；y 代表dy/dq；a代表理论轮廓的横坐标；b代表理论轮廓的纵坐标；m代表横坐标；n代表纵坐标Static i As SingleDim s As Single, q As Single, x As Single, y As Single, p As Single, rr As Single, a As Single, b As Single, m As Single, n As Singlerr = 63.013s = 64.213e = 131.557Picture1.CurrentX = 12.764Picture1.CurrentY = -95.175i = i + 0.001If i <= 3.14 / 4 Thenq = ix = (480 * 4 * q / 3.14 ^ 2) * Cos(q) + (s + 240 * (2 * q / 3.14) ^ 2) * (-Sin(q)) - e * Cos(q)y = (480 * 4 * q / 3.14 ^ 2) * Sin(q) - e * Sin(q) + (s + 240 * (2 * q / 3.14) ^ 2) * Cos(q)a = (s + 240 * (2 * q / 3.14) ^ 2) * Cos(q) - e * Sin(q)b = (s + 240 * (2 * q / 3.14) ^ 2) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pPicture1.PSet Step(m, -n), vbBlackElseIf i >= 3.14 / 4 And i <= 3.14 / 2 Thenq = ix = (2 * 97.367 * (3.14 / 2 - q)) * Cos(q) + (s + 120 - 97.367 * (3.14 / 2 - q) ^ 2) * (-Sin(q)) - e * Cos(q)y = (2 * 97.367 * (3.14 / 2 - q)) * Sin(q) + (s + 120 - 97.367 * (3.14 / 2 - q) ^ 2) * Cos(q) - e * Sin(q)a = (s + 120 - 97.367 * (3.14 / 2 - q) ^ 2) * Cos(q) - e * Sin(q)b = (s + 120 - 97.367 * (3.14 / 2 - q) ^ 2) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pPicture1.PSet Step(m, -n), vbBlackElseIf i >= 3.14 / 2 And i <= 5 / 6 * 3.14 Thenq = ix = (s + 120) * (-Sin(q)) - e * Cos(q)y = (s + 120) * Cos(q) - e * Sin(q)a = (s + 120) * Cos(q) - e * Sin(q)b = (s + 120) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pPicture1.PSet Step(m, -n), vbBlackElseIf i >= 5 / 6 * 3.14 And i <= 19 / 18 * 3.14 Thenq = ix = (-2 * 123.23 * (q - 5 / 6 * 3.14)) * Cos(q) + (s + 120 - 123.23 * (q - 5 / 6 * 3.14) ^ 2) * (-Sin(q)) - e * Cos(q)y = (-2 * 123.23 * (q - 5 / 6 * 3.14)) * Sin(q) + (s + 120 - 123.23 * (q - 5 / 6 * 3.14) ^ 2) * Cos(q) - e * Sin(q)a = (s + 120 - 123.23 * (q - 5 / 6 * 3.14) ^ 2) * Cos(q) - e * Sin(q)b = (s + 120 - 123.23 * (q - 5 / 6 * 3.14) ^ 2) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pPicture1.PSet Step(m, -n), vbBlackElseIf i >= 19 / 18 * 3.14 And i <= 23 / 18 * 3.14 Thenq = ix = (-2 * 123.23 * (23 / 18 * 3.14 - q)) * Cos(q) + (s + 123.23 * (23 / 18 *3.14 - q) ^ 2) * (-Sin(q)) - e * Cos(q)y = (-2 * 123.23 * (23 / 18 * 3.14 - q)) * Sin(q) + (s + 123.23 * (23 / 18 * 3.14 - q) ^ 2) * Cos(q) - e * Sin(q)a = (s + 123.23 * (23 / 18 * 3.14 - q) ^ 2) * Cos(q) - e * Sin(q)b = (s + 123.23 * (23 / 18 * 3.14 - q) ^ 2) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pPicture1.PSet Step(m, -n), vbBlackElseIf i >= 23 / 18 * 3.14 And i <= 2 * 3.14 Thenq = ix = s * (-Sin(q)) - e * Cos(q)y = s * Cos(q) - e * Sin(q)a = s * Cos(q) - e * Sin(q)b = s * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pPicture1.PSet Step(m, -n), vbBlackElseEnd IfEnd Sub内部轮廓为实际轮廓，外部轮廓为理论轮廓。