西工大机械原理大作业2凸轮机构作业
凸轮机构大作业 西工大机械原理
大作业（二）凸轮机构设计（题号：4-A）（一）题目及原始数据···············（二）推杆运动规律及凸轮廓线方程·········（三）程序框图·········（四）计算程序·················（五）程序计算结果及分析·············（六）凸轮机构图·················（七）心得体会··················（八）参考书···················一题目及原始数据试用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计（1）推程运动规律为五次多项式运动规律，回程运动规律为余弦加速度运动规律；（2）打印出原始数据；（3）打印出理论轮廓和实际轮廓的坐标值；（4）打印出推程和回程的最大压力角，以及出现最大压力角时凸轮的相应转角；（5）打印出凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角；（6）打印最后所确定的凸轮的基圆半径。
机械原理大作业凸轮
机械原理大作业二题目：凸轮机构设计19班号： 1408301学号： 1140830118姓名：高奎教师：焦映厚完成时间： 2016.6.11.从动件位移，速度，加速度图syms fai1fai2fai3=pi/3:0.01:8/9*pi;fai4=4/3*pi:0.01:2*pi;omiga=1;h=25;fai1=0:0.01:pi/3;fai2=160/180*pi:0.01:240/180*pi;s1=h*(fai1*3/pi-1/(2*pi)*sin(360/60*fai1));s2=h/2*(1+cos(180/80*(fai2-160/180*pi)));figure(1);subplot(3,1,1);subs(s1,'fai1',fai1);subs(s2,'fai2',fai2);plot(fai1,s1);hold on;plot(fai2,s2);plot(fai3,25);plot(fai4,0);xlabel('凸轮转角（rad）');ylabel('位移（mm）');title('位移图');v1=h*omiga*3/pi*(1-cos(360/60*fai1));v2=-h*omiga*180/160*sin(180/80*(fai2-160/180*pi));subplot(3,1,2);plot(fai1,v1,'g');hold on;plot(fai2,v2,'g');plot(fai3,0,'g');plot(fai4,0,'g');xlabel('凸轮转角(rad)');ylabel('速度(mm/s)');title('速度图');a1=2*pi*h*omiga^2/((60/180*pi)^2)*sin(360/60*fai1);a2=-(180/80)^2/2*h*omiga^2*cos(180/80*(fai2-160/180*pi)); subplot(3,1,3);plot(fai1,a1,'r');hold on;plot(fai2,a2,'r');plot(fai3,0,'r');plot(fai4,0,'r');xlabel('凸轮转角(rad)');ylabel('加速度(mm/s^2)');title('加速度图');2.类速度-位移图x1=60;t1=100;x2=80;t2=120;h=25;x1=x1*pi/180;x2=x2*pi/180;t1=t1*pi/180;t2=t2*pi/180; x= 0:0.001:60*pi/180;%升程s = h*(x/x1-sin(2*pi*x/x1)/(2*pi));k =-h*(1-cos(2*pi*x/x1))/x1;plot(k,s,'r'),hold on;x=160*pi/180:0.001:240*pi/180;%回程s = h*(1+cos(pi*(x-(x1+t1))/x2))/2;k = pi*h*sin(pi*(x-(x1+t1))/x2)/(2*x2);plot(k,s,'r'),hold on;%回程切线for i=-3.9:1:-3.9;f=@(k)k*tan(20/180*pi)+i;k =-50:0.1:50;s=f(k);plot(k,s),hold on;end%升程切线for i=-57:0.2:-57;f=@(k)-k*tan(55*pi/180)+i;k =-50:0.1:50;s=f(k);plot(k,s),hold on;endgrid onf=@(k)k*tan(55*pi/180);k=-50:0.1:0;s=f(k);plot(k,s);hold on;xlabel('ds/dψ');ylabel('s');title('类位移-速度图');3.压力角和曲率半径图figure(3);e=20;r0=102;s0=sqrt(r0.^2-e.^2);rs1=s0+s1;rs2=s0+s2;ang1=abs(atan((v1/omiga-e)./rs1))*180/pi; ang2=abs(atan((v2/omiga-e)./rs2))*180/pi; plot(fai1,ang1);hold on;plot(fai2,ang2);hold on;plot(fai3,9.0789);hold on;plot(fai4,11.5257);title('压力角图');h=25;t0=pi*60/180;t01=pi*80/180;ts=pi*100/180;ts1=pi*120/180;e=20;s0=100;t=0:0.001:60*pi/180;s=h*(t/t0-sin(2*pi*t/t0)/(2*pi));dx1 =(h/t0-h*cos(2*pi*t/t0)).*cos(t)-(s0+s).*sin(t)- e*cos(t); dy1=(h/t0-h*cos(2*pi*t/t0)).*sin(t)+(s0+s).*cos(t)- e*sin(t); p=sqrt(dx1.^2+dy1.^2);hold onplot(t,p);t=60*pi/180:pi/200:160*pi/180;s=h;dx2 =- sin(t).*(s + s0) - e*cos(t);dy2 =cos(t).*(s + s0) - e*sin(t);p=sqrt(dx2.^2+dy2.^2);hold onplot(t,p);t=160*pi/180:pi/200:240*pi/180;s=0.5*h*(1+cos(pi*(t-(t0+ts))/t01));dx3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*cos(t)-sin(t).*(s + s0) - e*cos(t);dy3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*sin(t)+ cos(t).*(s + s0) - e*sin(t);p=sqrt(dx3.^2+dy3.^2);hold onplot(t,p);t=240*pi/180:pi/200:2*pi;s=0;dx4 =- sin(t).*(s + s0) - e*cos(t);dy4 =cos(t).*(s + s0) - e*sin(t);p=sqrt(dx4.^2+dy4.^2);hold on;plot(t,p);hold off;title('曲率半径');grid on;4.凸轮理论轮廓和实际轮廓的绘制fai=0:0.01:2*pi;x1=60;t1=100;x2=80;t2=120;h=25;x1=x1.*pi./180;x2=x2.*pi./180;t1=t1.*pi./180;t2=t2.*pi./180; e=20;r0=102;s0=100;rr=20;%滚子半径x=0:pi/200:60.*pi/180;s = h.*(x./x1-sin(2.*pi.*x./x1)./(2.*pi));X1=(s0+s).*cos(x)-e.*sin(x);Y1=(s0+s).*sin(x)+e.*cos(x);X11=X1-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y11=Y1-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X1,Y1,'r',X11,Y11,'k'),hold on;plot(e*cos(fai),e*sin(fai));plot(r0*cos(fai),r0*sin(fai),'--g');x=60.*pi/180:pi/200:160.*pi/180;s=25;X2=(s0+s).*cos(x)-e.*sin(x);Y2=(s0+s).*sin(x)+e.*cos(x);X22=X2-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y22=Y2-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X2,Y2,'r',X22,Y22,'k'),hold on;x=160.*pi/180:pi/200:240.*pi/180;s=h.*(1+cos(pi.*(x-(x1+t1))./x2))./2;X3=(s0+s).*cos(x)-e.*sin(x);Y3=(s0+s).*sin(x)+e.*cos(x);X33=X3-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y33=Y3-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X3,Y3,'r',X33,Y33,'k'),hold on;x=240*pi/180:pi/200:2*pi;s=0;X4=(s0+s).*cos(x)-e.*sin(x);Y4=(s0+s).*sin(x)+e.*cos(x);X44=X4-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y44=Y4-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X4,Y4,'r',X44,Y44,'k'),hold on;x=240:pi/200:2*pi;X4=(s0+s).*cos(x)-e.*sin(x);Y4=(s0+s).*sin(x)+e.*cos(x);plot(X4,Y4,'b');legend('凸轮实际轮廓','凸轮理论轮廓','偏距圆','基圆');grid on;axis equal;。
凸轮机构大作业___西工大机械原理要点
大作业（二）凸轮机构设计（题号：4-A）（一）题目及原始数据···············（二）推杆运动规律及凸轮廓线方程·········（三）程序框图·········（四）计算程序·················（五）程序计算结果及分析·············（六）凸轮机构图·················（七）心得体会··················（八）参考书···················一题目及原始数据试用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计（1）推程运动规律为五次多项式运动规律，回程运动规律为余弦加速度运动规律；（2）打印出原始数据；（3）打印出理论轮廓和实际轮廓的坐标值；（4）打印出推程和回程的最大压力角，以及出现最大压力角时凸轮的相应转角；（5）打印出凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角；（6）打印最后所确定的凸轮的基圆半径。
机械原理大作业二凸轮(完整版) 20
Harbin Institute of Technology机械原理大作业二作业名称：凸轮机构设计设计题目： 20 院系：机电工程学院班级：设计者：学号：指导教师：**设计时间： 2014年5月哈尔滨工业大学一、设计题目如图1所示直动从动件盘形凸轮机构，其原始参数见表1。
表1 凸轮机构参数升程/h mm 升程运动角 0/ϕ 升程运动规律 升程许用压力角[]/α 回程运动角 0/ϕ'回程运动规律 回城许用压力角[]/α'远休止角/s ϕ 近休止角/s ϕ' 110 120 正弦加速度 3590 正弦加速度 65 50 100二、运动方程式及运动线图本实验假设凸轮逆时针旋转。
1.确定凸轮机构推杆升程、回程运动方程（ϕ为凸轮转角，ω为凸轮角速度）。
推程(余弦加速度)（203ϕπ≤≤）： 图1远休止段：当21738πϕπ≤≤时，110s h==，0v=，0a=回程(正弦加速度)（171389πϕπ≤≤）：近休止段：当1329πϕπ≤≤时，0s=，0v=，0a=2.绘制推杆位移、速度、加速度线图图2三、绘制凸轮机构的ds s d ϕ-线图，并由此确定凸轮的基圆半径和偏距图3图4如图3、图4所示，在ds d ϕ轴正侧（对应于推程），以tan(90[])α-为斜率做ds s d ϕ-曲线的切线1L ，在ds d ϕ轴负侧（对应于回程），以tan(90[])α'+为斜率做ds s d ϕ-曲线的切线2L ，再过点(0,0)O 做斜率为tan(90[])α+的直线3L ，则直线1L 、2L 、3L 与s 轴的夹角分别为[]α、[]α'、[]α。
显然，1L 、2L 、3L 三条直线下方的公共部分即为满足推程压力角不超过[]α和回程压力角不超过[]α'时，凸轮回转中心的可取区域。
记直线1L 与2L 的交点为12P ，直线1L 与3L 的交点为13P ，则最小基圆半径013r OP ==45.79mm ，对应的偏距e =26.27mm （点13P 到s 轴的距离），2200s r e =-=37.51mm 。
机械原理大作业二-凸轮机构设计..
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机械设计制造及其自动化班级：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 υ===。
机械原理大作业凸轮结构设计
机械原理大作业(二) 作业名称：机械原理设计题目：凸轮机构设计院系: 机电工程学院班级：设计者:学号:指导教师：丁刚陈明设计时间:哈尔滨工业大学机械设计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线得右下方为选择凸轮轴心得许用区。
机械原理大作业——凸轮.docx
大作业（二）凸轮机构设计题号：6班级：姓名：学号：同组者：成绩：完成时间：目录一凸轮机构题目要求 (1)二摆杆的运动规律及凸轮轮廓线方程 (2)三计算程序 (3)四运算结果及凸轮机构图 (9)4.1 第一组（A组）机构图及计算结果 (9)4.2 第二组（B组）机构图及计算结果 (14)4.3 第三组（C组）机构图及计算结果 (19)五心得体会 (24)第一组（A组） (24)第二组（B组） (24)第三组（C组） (24)六参考资料 (25)附录程序框图 (26)一凸轮机构题目要求（摆动滚子推杆盘形凸轮机构）题目要求：试用计算机辅助设计完成下列偏置直动滚子推杆盘形凸轮机构或摆动滚子推杆盘形凸轮机构的设计，已知数据如下各表所示。
凸轮沿逆时针方向作匀速转动。
表一摆动滚子推杆盘形凸轮机构的已知参数题号初选的基圆半径R0/mm机架长度Loa/mm摆杆长度Lab/mm滚子半径Rr/mm推杆摆角φ许用压力角许用最小曲率半径[ρamin][α1] [α2]A 15 60 55 10 24°35°70°0.3RrB 20 70 65 14 26°40°70°0.3RrC 22 72 68 18 28°45°65°0.35Rr 要求：1）凸轮理论轮廓和实际轮廓的坐标值2）推程和回程的最大压力角，及凸轮对应的转角3）凸轮实际轮廓曲线的最小曲率4）半径及相应凸轮转角5）基圆半径6）绘制凸轮理论廓线和实际廓线7）计算点数：N：72～120推杆运动规律：1)推程运动规律：等加速等减速运动2)回程运动规律：余弦加速度运动二摆杆的运动规律及凸轮轮廓线方程1）推程：1，运动规律：等加速等减速运动；2，轮廓线方程：A：等加速推程段设定推程加速段边界条件为：在始点处δ=0，s=0，v=0。
在终点处。
整理得：（注意：δ的变化范围为0~δ0/2。
机械原理凸轮大作业程序、图与设计说明说
一、 设计题目：27题如图所示直动从动件盘形凸轮机构，其原始参数见下表，据此设计该凸轮机构。
二、 凸轮推杆升程、回程运动方程和位移、速度、加速度线图对于不同运动规律的凸轮结构，其上升与下降的方式不一，但遵循同样的运动顺序：上升、远休止点恒定、下降、近休止点恒定。
因此设计它仅需确定这四个阶段的角度与位置即可。
1、对于凸轮的升程的运动方程如下所示：12[sin()]2s h ϕπϕφπφ=-12[1cos()]h v ωπϕφφ=- 21222sin()h a πωπϕφφ=升程 （mm ） 升程运动角（°） 升程运动规律 升程许用压力角（°） 回程运动角（°） 回程运动规律 回程许用压力角（°） 远休止角（°） 近休止角（°） 130 150正弦加速度30100余弦加速度6030802、对于凸轮的回程运动方程如下所示：01{1cos[()]}2s h s πϕφφφ=+-+1011sin[()]s h v πωπϕφφφφ=--+2210121cos[()]2s h a πωπϕφφφφ=--+3、推杆位移线图（横坐标凸轮转过角度，单位：°；纵坐标推杆位移，单位：mm ）4、推杆速度线图（横坐标图轮转过角度，单位°；纵坐标为推杆速度，单位mm/s)5、推杆加速度线图（横坐标为凸轮转过角度，单位：°；纵坐标为推杆加速度，单位：mm/s^2）三、 凸轮机构的ds sd ϕ-线图，凸轮的基圆半径和偏距1、根据ds s d ϕ-线图得无偏心距时最小基圆半径0116r =。
具体作图过程见下页。
2、查阅相关文献，知直动从动件盘形凸轮最小基圆半径求解公式为：00022(arctan)tan tan h r πππφαφα=-经计算得：0116.7835r =四、 滚子半径的确定及凸轮理论廓线和实际廓线的绘制查阅有关文献，得知在有偏置机构的凸轮最小基圆半径求解公式为 最小基圆半径：02cos p r r α=其对应的偏心距：0tan 2r e α=计算得：67.55p r =；33.78e =。
机械原理大作业凸轮机构题DOC
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：连杆机构运动分析院系：机械设计制造及其自动化班级：设计者：学号：指导教师：设计时间：一．设计题目设计直动从动件盘形凸轮机构，凸轮机构原始参数 序号 升程(mm) 升程运动角(º) 升程运 动规律升程许用压力角（º） 回程运动角（º）回程运 动规律回程许用压力角（º）远休止角 （º）近休止角 （º） 22 120 90等加等减速 4080等减等加速 70 70120二． 凸轮推杆运动规律1.运动规律（等加速等减速运动） 推程 0450≤≤ϕ2229602ϕπϕ=⎪⎪⎭⎫ ⎝⎛Φ=h s ϕπωϕω2219204=Φ=h v2220219204πωω=Φ=h a 推程 009045≤≤ϕ()222020)2(9601202ϕππϕ--=-ΦΦ-=hh s())2(1920422ϕπωπϕω-=-ΦΦ=h v222219204ωπω-=Φ-=h a2.运动规律（等加速等减速运动） 回程 00200160≤≤ϕ ()[]2222)98(9601202πϕπϕ--=Φ+Φ-Φ-=S h h s ()[])98(1920-4-22πϕωπϕω-=Φ+Φ-Φ=S h v 222219204ωπω-=Φ-=h a回程 00240200≤≤ϕ ()[]222'002)34(9602ϕππϕ-=-Φ+Φ+ΦΦ=S h s ()[])34(1920-4-2'002ϕπωπϕω-=-Φ+Φ+ΦΦ=S h v222219204ωπω=Φ=h a三．推杆位移、速度、加速度线图及凸轮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 = 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) ^ 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的斜率相等，因为kDtdt=tan600,右侧曲线斜率可以表示为：q;q=tan600继而求出切点坐标（337.272,292.084）。
哈工大机械原理大作业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。
机械原理大作业二-丁烨
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
机械原理大作业二
课程名称：机械原理
设计题目：凸轮机构设计
院系：机电工程学院
班级：1008301
设计者：丁烨
学号：1100801020
指导教师：丁刚
设计时间：2010年6月4日
（一）、设计题目
1、机构运动简图：
2、凸轮机构的原始参数：
选择第32组，数据如下
(二)、设计过程
（1）、计算推杆的位移、速度和加速度线图
位移-转角曲线
速度-转角曲线
加速度-转角曲线
(2)、绘制凸轮机构的
s
d
ds
-
ϕ线图
s
d ds -ϕ
线图
（3）、确定凸轮的基圆半径和偏距
基圆半径为r 0 = 142mm ，偏距e = 20mm 。
（4）、绘制凸轮理论廓线
（5）、确定滚子半径并绘制凸轮实际轮廓线最小曲率半径 = 90.00mm
滚子半径= 10mm
得到的内包络轮廓线为：
取内包络轮廓线为凸轮的工作轮廓曲线。
机械原理大作业-凸轮机构
二、凸轮机构一、运动分析凸轮的运动分为4个阶段：推程运动、远休程、回程运动、近休程。
该凸轮机构4个阶段的运动角分别为推程运动角90˚、远休止角100 ˚、回程运动角50 ˚、近休止角120 ˚。
推程运动阶段的运动规律为正弦加速度运动，回程运动的运动规律为4-5-6-7多项式运动。
凸轮的简图如图1所示。
图1对该机构进行简单的运动分析：1.升程阶段采用正弦加速度的运动规律，从动件的位移、速度、加速度、压力角的计算公式如下：计算时将相应的量带入公式即可得到。
类速度可以直接将位移方程对凸轮转角ϕ求导得到。
2.远休程阶段的位移不变，与凸轮升程阶段最后的位移相等，速度、加速度则变为0。
3.回程阶段位移、速度、加速度可通过代入4-5-6-7多项式的方程求得。
4.近休程阶段的位移与回程阶段最后的位移相等，且为0，速度、加速度均变为0.二、流程框图图2三、运用VC编程#include<stdio.h>#include<math.h>#define pi 3.141592654 //定义全局变量int main() //主函数{int i,j,k,l;double s; //定义位移量double v; //定义速度量double a; //定义加速度量double r; //定义弧度制角度量double d,o,m,t=40,x1,x2,y1,y2,d1,d2; //定义中间变量double p; //定义角度制角度量double w=1; //定义并角速度量赋值double R=50; //定义基圆半径double e=30; //定义偏距double n; //定义压力角double u; //定义曲率半径double Rr=17; //定义滚子半径并赋值double x,y,X,Y; //定义实际与理论廓线上点的坐标r=0;for(i=0;i<20;i++){s=20/pi*(4*r-sin(4*r));x=-(t+s)*sin(r)-e*cos(r);y=(t+s)*cos(r)-e*sin(r);d1=-(s+t)*cos(r)+e*sin(r);d2=-(s+t)*sin(r)-e*cos(r);X=x-Rr*d2/pow(d1*d1+d2*d2,0.5);Y=y+Rr*d1/pow(d1*d1+d2*d2,0.5);d=80/pi*(1-cos(4*r));v=80/pi*(1-cos(4*r));a=320/pi*sin(4*r);m=atan(fabs(d-e)/(s+t));n=180*m/pi;x1=(t+s)*cos(r)+v/w*sin(r)-e*sin(r);y1=-(t+s)*sin(r)+v/w*cos(r)-e*cos(r);x2=-(t+s)*sin(r)+v/w*cos(r)+a/(w*w)*sin(r)+v/w*cos(r)-e*cos(r);y2=-(t+s)*cos(r)-v/w*sin(r)+a/(w*w)*cos(r)-v/w*sin(r)+e*sin(r);u=pow(x1*x1+y1*y1,1.5)/fabs(x1*y2-y1*x2);r=r+pi/40;p=180/pi*r;printf("%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf\n", p,s,v,a,d,n,u,x,y,X,Y);}r=pi/2;for(j=0;j<5;j++){s=s;x=-(t+s)*sin(r)-e*cos(r);y=(t+s)*cos(r)-e*sin(r);d1=-(s+t)*cos(r)+e*sin(r);d2=-(s+t)*sin(r)-e*cos(r);X=x-Rr*d2/pow(d1*d1+d2*d2,0.5);Y=y+Rr*d1/pow(d1*d1+d2*d2,0.5);d=0;v=0;a=0;m=atan(fabs(d-e)/(s+t));n=180*m/pi;x1=(t+s)*cos(r)+v/w*sin(r)-e*sin(r);y1=-(t+s)*sin(r)+v/w*cos(r)-e*cos(r);x2=-(t+s)*sin(r)+v/w*cos(r)+a/(w*w)*sin(r)+v/w*cos(r)-e*cos(r);y2=-(t+s)*cos(r)-v/w*sin(r)+a/(w*w)*cos(r)-v/w*sin(r)+e*sin(r);u=pow(x1*x1+y1*y1,1.5)/fabs(x1*y2-y1*x2);r=r+pi/9;p=180/pi*r;printf("%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf\n", p,s,v,a,d,n,u,x,y,X,Y);}r=(19*pi)/18;for(k=0;k<20;k++){o=(18*r-19*pi)/(5*pi);s=40*(1-35*pow(o,4)+84*pow(o,5)-70*pow(o,6)+20*pow(o,7));x=-(t+s)*sin(r)-e*cos(r);y=(t+s)*cos(r)-e*sin(r);d1=-(s+t)*cos(r)+e*sin(r);d2=-(s+t)*sin(r)-e*cos(r);X=x-Rr*d2/pow(d1*d1+d2*d2,0.5);Y=y+Rr*d1/pow(d1*d1+d2*d2,0.5);d=18*40/5/pi*(-35*4*pow(o,3)+84*5*pow(o,4)-70*6*pow(o,5)+20*7*pow(o,6));v=-80/pi*(140*pow(o,3)-420*pow(o,4)+420*pow(o,5)-140*pow(o,6));a=-160/pi*(420*pow(o,2)-1680*pow(o,3)+2100*pow(o,4)-840*pow(o,5));m=atan(fabs(d-e)/(s+t));n=180*m/pi;x1=(t+s)*cos(r)+v/w*sin(r)-e*sin(r);y1=-(t+s)*sin(r)+v/w*cos(r)-e*cos(r);x2=-(t+s)*sin(r)+v/w*cos(r)+a/(w*w)*sin(r)+v/w*cos(r)-e*cos(r);y2=-(t+s)*cos(r)-v/w*sin(r)+a/(w*w)*cos(r)-v/w*sin(r)+e*sin(r);u=pow(x1*x1+y1*y1,1.5)/fabs(x1*y2-y1*x2);r=r+pi/72;p=180/pi*r;printf("%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf\n", p,s,v,a,d,n,u,x,y,X,Y);}r=(4*pi)/3;for(l=0;l<5;l++){s=s;x=-(t+s)*sin(r)-e*cos(r);y=(t+s)*cos(r)-e*sin(r);d1=-(s+t)*cos(r)+e*sin(r);d2=-(s+t)*sin(r)-e*cos(r);X=x-Rr*d2/pow(d1*d1+d2*d2,0.5);Y=y+Rr*d1/pow(d1*d1+d2*d2,0.5);d=0;v=0;a=0;m=atan(fabs(d-e)/(s+t));n=180*m/pi;x1=(t+s)*cos(r)+v/w*sin(r)-e*sin(r);y1=-(t+s)*sin(r)+v/w*cos(r)-e*cos(r);x2=-(t+s)*sin(r)+v/w*cos(r)+a/(w*w)*sin(r)+v/w*cos(r)-e*cos(r);y2=-(t+s)*cos(r)-v/w*sin(r)+a/(w*w)*cos(r)-v/w*sin(r)+e*sin(r);u=pow(x1*x1+y1*y1,1.5)/fabs(x1*y2-y1*x2);r=r+2*pi/15;p=180/pi*r;printf("%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf\n", p,s,v,a,d,n,u,x,y,X,Y);}return 0;}四、计算结果处理1.输出数据位移s、速度v、加速度a、类速度ds/dϕ、压力角α、曲率半径ρ（其中曲率半径缺失的数据为太大而不合题意的数据，已将其舍去）:表1凸轮轮廓：理论廓线坐标、实际廓线坐标：表22.根据输出数据做出图像：图2图3图4图5图6图7图8。
机械原理大作业——凸轮
欢迎共阅大作业（二）凸轮机构设计题号： 6班级：姓名：学号：同组者：成绩：完成时间：目录一凸轮机构题目要求 (1)二摆杆的运动规律及凸轮轮廓线方程 (2)三计算程序 (3)四运算结果及凸轮机构图 (9)4.1 第一组（A组）机构图及计算结果 (9)4.2 第二组（B组）机构图及计算结果 (14)4.3 第三组（C组）机构图及计算结果 (19)五心得体会 (24)第一组（A组） (24)第二组（B组） (24)第三组（C组） (24)六参考资料 (25)附录程序框图 (26)一凸轮机构题目要求（摆动滚子推杆盘形凸轮机构）题目要求：试用计算机辅助设计完成下列偏置直动滚子推杆盘形凸轮机构或摆动滚子推杆盘形凸轮机构的设计，已知数据如下各表所示。
凸轮沿逆时针方向作匀速转动。
表一摆动滚子推杆盘形凸轮机构的已知参数题号初选的基圆半径R0/mm机架长度Loa/mm摆杆长度Lab/mm滚子半径Rr/mm推杆摆角φ许用压力角许用最小曲率半径[ρamin][α1] [α2]A 15 60 55 10 24°35°70°0.3RrB 20 70 65 14 26°40°70°0.3RrC 22 72 68 18 28°45°65°0.35Rr 要求：1）凸轮理论轮廓和实际轮廓的坐标值2）推程和回程的最大压力角，及凸轮对应的转角3）凸轮实际轮廓曲线的最小曲率4）半径及相应凸轮转角5）基圆半径6）绘制凸轮理论廓线和实际廓线7）计算点数：N：72～120推杆运动规律：1)推程运动规律：等加速等减速运动2)回程运动规律：余弦加速度运动二摆杆的运动规律及凸轮轮廓线方程1）推程：1，运动规律：等加速等减速运动； 2，轮廓线方程： A ：等加速推程段设定推程加速段边界条件为： 在始点处 δ=0，s=0，v=0。
在终点处 h /2 s ,2/==δοδ。
机械原理课程设计凸轮机构
Part Three
机械原理课程设计 凸轮机构方案
设计目的和要求
设计目的：掌握凸轮机构的基本原 理和设计方法
设计内容：包括凸轮机构的设计、 制造、装配和调试
添加标题
添加标题
添加标题
添加标题
设计要求：满足凸轮机构的运动要 求，如速度、加速度、行程等
设计步骤：明确设计任务、选择设 计方案、进行设计计算、绘制设计 图纸、制作模型、进行实验验证等
凸轮轮廓曲线的设计方法包括解析法、图 解法和计算机辅助设计等。
凸轮轮廓曲线的设计需要满足凸轮机构 的运动规律、负载、速度、加速度等要 求，同时需要考虑到凸轮的制造工艺和 成本等因素。
凸轮机构压力角计算
压力角定义：凸轮与从动件接触点 处法线与凸轮轮廓线之间的夹角
压力角影响因素：凸轮轮廓线形状、 从动件形状、凸轮半径、从动件半 径
凸轮机构工作原理
凸轮机构通过凸轮与从动件 的接触，实现从动件的位移 和运动
凸轮机构由凸轮、从动件和 机架组成
凸轮机构的工作原理是利用 凸轮的轮廓曲线，使从动件
产生预定的运动
凸轮机构的应用广泛，如汽 车、机床、机器人等领域
凸轮机构分类
按照凸轮运动规律分类：等 速运动凸轮、等加速运动凸 轮、等减速运动凸轮等
Part Six
凸轮机构运动仿真 与优化
运动仿真模型的建立
确定凸轮机构的类型和参数 建立凸轮机构的三维模型 设定运动仿真的初始条件和边界条件 设定运动仿真的时间步长和仿真时间 设定运动仿真的输出变量和观察点 运行运动仿真，观察仿真结果，并进行优化
运动仿真结果分析
凸轮机构运动仿 真结果：包括位 移、速度、加速 度等参数
凸轮从动件的类 型：滚子从动件、 滑块从动件、圆 柱从动件等
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
机械原理大作业凸轮..
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮结构设计院系：机电工程学院班级：1308108设计者：仲星光学号：1130810816指导教师：林琳设计时间：2015年6月7日一、设计题目如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
行程（mm)升程运动角（°）升程运动规律升程许用压力角（°）回程运动角（°）回程运动规律回程许用压力角（°）远休止角（°）近休止角（°）6080余弦加速度30 60 摆抛摆60 100 120二．数学计算方法（设计书中所给出的回程运动方程有误）1.理论轮廓和工作轮廓廓线方程正偏置平面凸轮，反转法，小滚子中心在反转运动中的轨迹即为凸轮的理论轮廓，推杆所在直线一直与偏心圆相切，由几何关系可得盘形凸轮理论轮廓线方程为：滚子从动件盘形凸轮的实际轮廓线是以理论轮廓上各点为圆心，以小滚子半径为半径的圆族的包络线，理论廓线B点相对应的实际廓线B'点的方程为：2.压力角3.凸轮轮廓曲率半径由数学分析知，凸轮理论轮廓曲线上任一点的曲率半径的计算公式为：（以上方程来自《基于MatLab语言的机构设计与分析》上海科学技术出版社）三．MatLab程序程序设计流程1.主程序：推杆位clc;clear;h=60;%行程phi01=80;%推程运动角phis1=100;%远休止角phi02=60;%回程运动角phis2=120;%近休止角alpha1=pi/6;%升程许用压力角alpha2=pi/3;%回程许用压力角omiga=1;%凸轮角速度设为1%计算凸轮转角、从动件位移，速度和加速度%绘制从动件位移、速度、加速度线图[psi,s,v,a]=CanShuJiSuan(h,phi01,phis1,phi02,phis2,... omiga,alpha1,alpha2);e=22;r0=57;%由ds/dpsi-s图像确定基圆半径和偏距%计算凸轮轮廓线曲率半径及压力角%绘制凸轮理论轮廓线上的压力角线图和曲率半径图[ang,rou,DxDpsi,DyDpsi]=YaLiJiao_QuLvBanJin( r0,psi,s, v,a,e,omiga );Rr=9;%由曲率半径最小值确定小滚子半径%计算凸轮轮廓曲线%绘制理论和实际轮廓曲线LunKuoXian( r0,psi,s,e,Rr,DxDpsi,DyDpsi );2.余弦加速运动规律子函数function [ s1,v1,a1,psi1 ] = Yuxian( phi01,h,omiga )%计算余弦加速度运动规律psi1=linspace(0,phi01,round(phi01));s1=(h/2).*[1 - cos(pi.*psi1./phi01)];v1=(pi*h*omiga./(2.*(phi01*pi/180))) .*sin(pi.*psi1./phi01);a1=pi^2*h*omiga^2/(2*(phi01*pi/180)^2)*cos(pi.*psi1./p hi01)end3.远休止程子函数function [ s2,v2,a2,psi2 ] = YuanXiu( phi01,phis1,h ) psi2=linspace(phi01+1,phi01+phis1,round(phis1));s2=h*psi2./psi2;v2=0*psi2;a2=0*psi2;end4.摆抛摆子函数function [ s3,v3,a3,psi3 ] =BaiPaoBai( phi01,phis1,phi02,h,omiga )%计算回程摆线-抛物线-摆线运动规律%计算当phi01+phis1<psi<=phi01+phis1+phi02/8时，从动件的位移、速度、加速度psi3_1=linspace(phi01+phis1+phi02/8/100,phi01+phis1+ph i02/8,100);s3_1=h-h*(2*(psi3_1-phi01-phis1)/phi02-sin(4*pi*(psi3_ 1-phi01-phis1)/...phi02)/(2*pi))/(2+pi);v3_1=-2*h*omiga*(1-cos(4*pi*(psi3_1-phi01-phis1)/phi02 ))/((2+pi)*...phi02*pi/180);a3_1=-8*h*pi*omiga^2*sin(4*pi*(psi3_1-phi01-phis1)/phi 02)/((2+pi)*...(phi02*pi/180)^2);%计算当phi01+phis1<psi<=phi01+phis1+3*phi02/8时，从动件的位移、速度、加速度psi3_2=linspace(phi01+phis1+phi02/8+phi02/4/100,phi01+ phis1+3*phi02/8,100);s3_2=h-h*(4*pi*((psi3_2-phi01-phis1).^2)/(phi02^2)-(pi -2)*(psi3_2-phi01...-phis1)/phi02+pi/16-1/(2*pi))/(2+pi);v3_2=-h*omiga*(8*pi*(psi3_2-phi01-phis1)/phi02-pi+2)/( (2+pi)*phi02*pi/180);a3_2=-8*h*pi*omiga^2/((2+pi)*(phi02*pi/180)^2)*psi3_2. /psi3_2;%计算当phi01+phis1<psi<=phi01+phis1+5*phi02/8时，从动件的位移、速度、加速度psi3_3=linspace(phi01+phis1+3*phi02/8+phi02/4/100,phi0 1+phis1+5*phi02/8,100);s3_3=h-h*(2*(pi+1)*(psi3_3-phi01-phis1)/phi02-pi/2-sin (4*pi*(psi3_3-phi01...-phis1)/phi02-pi)/(2*pi))/(2+pi);v3_3=-2*h*omiga*(pi+1-cos(4*pi*(psi3_3-phi01-phis1)/ph i02-pi))/((2+pi)...*phi02*pi/180);a3_3=-8*h*pi*omiga^2*sin(4*pi*(psi3_3-phi01-phis1)/phi 02-pi)/((2+pi)*...(phi02*pi/180)^2);%计算当phi01+phis1<psi<=phi01+phis1+7*phi02/8时，从动件的位移、速度、加速度psi3_4=linspace(phi01+phis1+5*phi02/8+phi02/4/100,phi0 1+phis1+7*phi02/8,100);s3_4=h-h*(-4*pi*((psi3_4-phi01-phis1).^2)/(phi02^2)+(7 *pi+2)*(psi3_4-...phi01-phis1)/phi02-33*pi/16+1/(2*pi))/(2+pi);v3_4=-h*omiga*(-8*pi*(psi3_4-phi01-phis1)/phi02+7*pi+2 )/((2+pi)*phi02*pi/180);a3_4=8*h*pi*omiga^2/((2+pi)*(phi02*pi/180)^2)*psi3_4./ psi3_4;%计算当phi01+phis1<psi<=phi01+phis1+phi02时，从动件的位移、速度、加速度psi3_5=linspace(phi01+phis1+7*phi02/8+phi02/8/100,phi0 1+phis1+phi02,100);s3_5=h-h*(2*(psi3_5-phi01-phis1)/phi02+pi-sin(4*pi*(ps i3_5-phi01-phis1).../phi02-2*pi)/(2*pi))/(2+pi);v3_5=-2*h*omiga*(1-cos(4*pi*(psi3_5-phi01-phis1)/phi02 -2*pi))/((2+pi)...*phi02*pi/180);a3_5=-8*h*pi*omiga^2*sin(4*pi*(psi3_5-phi01-phis1)/phi 02-2*pi)/((2+pi)...*(phi02*pi/180)^2);%输出回程阶段凸轮的转角、从动件的位移速度加速度的数组s3=[s3_1,s3_2,s3_3,s3_4,s3_5];v3=[v3_1,v3_2,v3_3,v3_4,v3_5];a3=[a3_1,a3_2,a3_3,a3_4,a3_5];psi3=[psi3_1,psi3_2,psi3_3,psi3_4,psi3_5];end5.近休止程子函数function [ s4,v4,a4,psi4 ] =JinXiu( phi01,phis1,phi02,phis2 )psi4=linspace(phi01+phis1+phi02+1,phi01+phis1+phis2+ph i02,round(phis2));s4=0*psi4;v4=0*psi4;a4=0*psi4;End6.参数计算及图像绘制[s1,v1,a1,psi1]=YuXian(phi01,h,omiga);%推程[s2,v2,a2,psi2]=YuanXiu(phi01,phis1,h);%远休程[s3,v3,a3,psi3]=BaiPaoBai(phi01,phis1,phi02,h,omiga);%回程[s4,v4,a4,psi4]=JinXiu(phi01,phis1,phi02,phis2);%近休程psi=[psi1,psi2,psi3,psi4];%凸轮转角s=[s1,s2,s3,s4];%从动件位移v=[v1,v2,v3,v4];%从动件速度a=[a1,a2,a3,a4];%从动件加速度%绘制从动件位移、速度、加速度线图figure(1)subplot(3,1,1)plot(psi,s);grid onxlabel('凸轮转角(度)');ylabel('位移（mm）');subplot(3,1,2);plot(psi,v);grid onxlabel('凸轮转角(度)');ylabel('速度（mm/s）');subplot(3,1,3)plot(psi,a);grid onxlabel('凸轮转角(度)');ylabel('加速度(mm/s^2)');%绘制ds/dphi-s线图，并确定凸轮基圆半径和偏距DsDphi=v/omiga;figure(2)plot(DsDphi,s);hold onx1=-150:100;y1=tan(pi/2-alpha1)*(x1-65.26)+22.33; plot(x1,y1);hold ony2=tan(pi/2+alpha2)*(x1+107.6)+22.36; plot(x1,y2);7.轮廓线绘制if Rr==0x=x0;y=y0;elseA=sqrt(DxDpsi.^2+DyDpsi.^2);x=x0+Rr*DyDpsi./A;y=y0-Rr*DxDpsi./A;endfigure(4)plot(r0.*cos(psi.*pi/180),r0.*sin(psi.*pi/180),'-.',.. .x0,y0,'--',x,y,e*cos(psi.*pi/180),e*sin(psi.*pi/180)); grid on;legend('基圆','凸轮理论轮廓','凸轮实际轮廓','偏距圆');axis equalend8.压力角曲线及曲率半径function[ang,rou,DxDpsi,DyDpsi]=YaLiJiao_QuLvBanJin( r0,psi,s, v,a,e,omiga )s0=sqrt(r0.^2-e.^2);rs1=s0+s;ang=abs(atan((v/omiga-e)./rs1))*180/pi;DxDpsi=(v./omiga-e).*sin(psi.*pi/180)+...(s0+s).*cos(psi.*pi/180);DyDpsi=(v./omiga-e).*cos(psi.*pi/180)-...(s0+s).*sin(psi.*pi/180);DDxDpsi=(a./(omiga^2)-(s0+s)).*sin(psi.*pi/180)... +(2*v./omiga-e).*cos(psi.*pi/180);DDyDpsi=(a./(omiga^2)-(s0+s)).*cos(psi.*pi/180)... +(2*v./omiga-e).*sin(psi.*pi/180);A=(DxDpsi.^2+DyDpsi.^2).^1.5;B=abs(DxDpsi.*(DDyDpsi)-DyDpsi.*(DDxDpsi));rou=A./B;%凸轮理论轮廓曲率半径%绘图figure(3)axis onplotyy(psi,rou,psi,ang);%绘制曲率半径图[AX]=plotyy(psi,rou,psi,ang);set(AX(1),'yTick',[0:20:200]);set(AX(2),'yTick',[0:10:100]);end四．程序运行结果与生成图像。
机械原理大作业 凸轮机构程序
theta_1=input('请输入凸轮近休止角 theta1=');theta_2=input('请输入凸轮推程止角 theta2=');theta_3=input('请输入凸轮远修止角 theta3=');r_0=input('请输入初选基圆半径 r0=');l_oa=input('请输入机架长度 l_OA=');l_ab=input('请输入摆杆长度 l_AB=');r_r=input('请输入滚子半径 r_r=');phi=input('请输入推杆摆角 phi=');alpha_1=input('请输入许用压力角 aloha_1=');alpha_2=input('请输入许用压力角 aloha_2=');rho_amin=input('请输入许用最小曲率半径半径/r_r=')*r_r;delta_0=0.005;omega=1;delta_1=theta_1;delta_2=theta_2-theta_1;delta_3=theta_3-theta_2;delta_4=360-theta_3;temp=0;while(temp==0)phi_0=acosd((l_ab^2+l_oa^2-r_0^2)/(2*l_ab*l_oa));temp=1;%最大压力角及相应凸轮转角初值alpha_1max=0;alpha_2max=0;delta_1max=0;delta_2max=0;%最小曲率半径及相应凸轮转角初值rho_0min=r_0-r_r;delta_min=0;for i=0:120delta=i*3;if delta>=0&delta<theta_1 %近休s=0;v=0;a=0;elseif delta>=theta_1&delta<theta_2 %推程s=phi*((delta-theta_1)/delta_2-sin(2*pi*(delta-theta_1)/delta_2)/(2*pi)); v=phi*omega*(1-cos(2*pi*(delta-theta_1)/delta_2))/delta_2;a=2*pi*phi*omega^2*sin(2*pi*(delta-theta_1)/delta_2)/delta_2^2;elseif delta>=theta_2&delta<theta_3 %远休s=phi;v=0;elseif delta>=theta_3&delta<(theta_3+delta_4/2) %回程加速s=phi-2*phi*(delta-theta_3)^2/delta_4^2;v=-4*phi*omega*(delta-theta_3)/delta_4^2;a=-4*phi*omega^2/delta_4^2;elseif delta>=(theta_3+delta_4/2)&delta<=360 %回程减速s=2*phi*(delta_4-(delta-theta_3))^2/delta_4^2;v=-4*phi*omega*(delta_4-(delta-theta_3))/delta_4^2;a=4*phi*omega^2/delta_4^2;endx=l_oa*sind(delta)-l_ab*sind(delta+s+phi_0); %理论轮廓方程y=l_oa*cosd(delta)-l_ab*cosd(delta+s+phi_0);d_x=l_oa*cosd(delta)-l_ab*cosd(delta+s+phi_0); %理论方程一阶导数d_y=l_oa*(-sind(delta))+l_ab*sind(delta+s+phi_0);dd_x=-l_oa*sind(delta)+l_ab*sind(delta+s+phi_0); %理论方程二阶导数dd_y=-l_oa*cosd(delta)+l_ab*cosd(delta+s+phi_0);x_0=x-r_r*(-d_y/sqrt(d_x^2+d_y^2)); %实际轮廓方程y_0=y-r_r*(d_x/sqrt(d_x^2+d_y^2));if delta>=theta_2&delta<theta_3 %推程压力角alpha_11=atand((l_ab*abs(v/omega)-(l_oa*cosd(phi_0+s)-l_ab))/(l_oa*sind(phi_0+s) ));%推程压力角方程if alpha_11<0alpha_11=-alpha_11;endALPHA_11(i+1)=alpha_11;if alpha_11>alpha_1r_0=r_0+delta_0;temp=0;break;else if alpha_11>alpha_1maxalpha_1max=alpha_11;delta_1max=delta;endendelse if delta>=theta_3&delta<=360 %回程压力角alpha_22=atand((l_ab*abs(v/omega)+(l_oa*cosd(phi_0+s)-l_ab))/(l_oa*sind(phi_0+s) ));%回程压力角方程if alpha_22<0alpha_22=-alpha_22;endALPHA_22(i+1)=alpha_22;if alpha_22>alpha_2r_0=r_0+delta_0;temp=0;break;else if alpha_22>alpha_2maxalpha_2max=alpha_22;delta_2max=delta;endendendendrho=1/((abs(d_x*dd_y-dd_x*d_y))/((d_x^2+d_y^2)^(3/2))); %曲率半径方程if rho>=rho_aminif rho<rho_0minrho_0min=rho;delta_min=delta;endelser_0=r_0+delta_0;temp=0;breakendDELTA(i+1)=delta;X(i+1)=x;Y(i+1)=y;X_0(i+1)=x_0;Y_0(i+1)=y_0;S(i+1)=s;V(i+1)=v;A(i+1)=a;RHO(i+1)=rho;endendfigure(1);axis equalhold ont=0:360;X_1=r_0*cosd(t); %画基圆Y_1=r_0*sind(t);X_2=r_r*cosd(t)+X(10); %画滚子Y_2=r_r*sind(t)+Y(10);plot(X_1,Y_1,'m--',X,Y,':',X_0,Y_0,'k',X_2,Y_2,'c--'); legend('基圆','理论轮廓','实际轮廓');plot(0,0,'ko'); %固定凸轮点plot(X(10),Y(10),'ko'); %固定滚子点title('凸轮轮廓曲线图');xlabel('x/mm');ylabel('y/mm');hold offfigure(2);plot(DELTA,S);title('凸轮运动规律曲线图——摆角');xlabel('{\delta}/(^o)');ylabel('{\phi}/(^o)');grid onfigure(3);plot(DELTA,V);title('凸轮运动规律曲线图——角速度');xlabel('{\delta}/(^o)');ylabel('{\omega}/(rad/s)');grid onfigure(4);plot(DELTA,A);title('凸轮运动规律曲线图——角加速度');xlabel('{\delta}/(^o)');ylabel('{\alpha}/(m/s^2)');grid ondisp('推程最大压力角：');alpha_1maxdisp('推程最大压力角相应转角：');delta_1maxdisp('回程最大压力角：');alpha_2maxdisp('回程最大压力角相应转角：');delta_2maxdisp('凸轮最小曲率半径：');rho_0mindisp('凸轮最小曲率半径相应转角：');delta_mindisp('最后确定的凸轮基圆半径：');r_0。