哈工大机械原理大作业凸轮机构设计第题
哈工大机械原理大作业凸轮07(可打印修改) (2)
f=linspace(0,2*pi,200); for n=1:length(f)
k1=subs(abs((tx1*yxx1-txx1*yx1)/(tx1^2+yx1^2)^1.5),{x1},{xx1}); v=[v,1/k1]; end
精品内容
s2=70; t2=(s2+s0).*cos(x2)-e*sin(x2); y2=(s0+s2).*sin(x2)+e*cos(x2); tx2=diff(t2,x2); txx2=diff(t2,x2,2); yx2=diff(y2,x2); yxx2=diff(y2,x2,2); for xx2=(pi/3):(pi/100):(175*pi/180);
k3=subs(abs((tx3*yxx3-txx3*yx3)/(tx3^2+yx3^2)^1.5),{x3},{xx3}); v=[v,1/k3]; end
s4=0; t4=(s4+s0).*cos(x4)-e*sin(x4); y4=(s0+s4).*sin(x4)-e*cos(x4); tx4=diff(t4,x4); txx4=diff(t4,x4,2); yx4=diff(y4,x4); yxx4=diff(y4,x4,2); for xx4=(265*pi/180):(pi/100):(2*pi);
elseif f(n)>f01+fs1+f02 && f(n)<=f01+fs1+f02+fs2 s(n)=0; v(n)=0;
哈工大机械原理大作业 第25题 凸轮设计
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院指导教师：林琳刘福利设计时间：2014.06哈尔滨工业大学凸轮机构设计1．设计题目(1)凸轮机构运动简图：2．确定凸轮推杆升程、回程运动方程并绘制推杆位移,速度,加速度线图 ：(1) 凸轮推杆升程，回程运动方程如下：A.推杆升程方程s (ф)=2h (φ∅0)2 = 1040π2ф2,0≤ ф≤π4; s (ф)=130−1040π2(π2−ф)2,π4≤ ф≤π2;v (ф)=4ℎω102= 2080w 2ф, 0≤ ф≤π; v (ф)=2080w 2(π−ф), π≤ ф≤π; a (ф)=12022π2≤ ф≤π4; a (ф)=−2080w 2π2, π4≤ ф≤π2; B.推杆回程方程s (ф)=80-230130]801501[130])(1['00ϕϕϕϕϕϕ⨯=--=+--s h 130(238−9ф4π),5π6≤ ф≤23π18;v (ф)=−5852πw,5π6≤ ф≤23π18; a (ф)=0，5π6≤ ф≤23π18;(2)推杆位移,速度,加速度线图如下：A.推杆位移线图(使用matlab 画图，程序详见附录1)B.推杆速度线图(使用matlab 画图，程序详见附录2)C.推杆加速度线图(使用matlab画图，程序详见附录3)-s线图，并依次确定凸轮的基圆半径和偏距. 3．凸轮机构的dsdф(1) 凸轮机构的ds-s线图：(使用matlab画图，程序详见附录4)dф(2)确定凸轮的基圆半径和偏距：以ds/df-s(f)图为基础，可分别作出三条限制线（推程许用压力角的切界限Dt d t ,回程许用压力角的限制线Dt 'dt'，起始点压力角许用线Bd'')，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。
得图如下： (使用matlab画图，程序详见附录6)得最小基圆对应的坐标位置大约为（55.28，-65.88）经计算取偏距e=55mm，r0=90mm.4．确定滚子半径及绘制凸轮理论轮廓曲线和实际轮廓曲线.为求滚子许用半径，须确定最小曲率半径，以防止凸轮工作轮廓出现尖点或出现相交包络线，确定最小曲率半径数学模型如下：)/)(/()/)(/(])/()/[(22222/322ϕϕϕϕϕϕρd x d d dy d y d d dx d dy d dx -+=其中：ϕϕϕϕcos )(sin ])/[(/0s s e d ds d dx ++-= ϕϕϕϕsin )(cos ])/[(/0s s e d ds d dy +--=ϕϕϕϕϕsin ])/[(cos ])/(2[/02222s s d s d e d ds d x d --+-= ϕϕϕϕϕcos ])/[(sin ])/(2[/02222s s d s d e d ds d y d --+--=利用上式可求的最小曲率半径后可确定实际廓线。
哈工大机械原理大作业-凸轮机构设计
哈工大机械原理大作业-凸轮机构设计(第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ϕϕ⨯=⨯通过编程求其角度。
哈工大机械原理大作业-凸轮结构设计 -05
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院班级：完成者：学号：指导教师：设计时间：设计说明书1 设计题目：如图所示直动从动件盘形凸轮机构，其原始参数见表2-1。
从表2-1中选择一组凸轮机构的原始参数，据此设计该凸轮机构。
注： 5题参数有误，我将远休止角改为45 ，将近休止角改为105 .2、推杆升程、回程运动方程及位移、速度、加速度线图：2.1凸轮运动分析： 推程运动方程; 01cos 2h s πϕ⎡⎤⎛⎫=-⎢⎥ ⎪Φ⎝⎭⎣⎦100sin 2h v πωπϕ⎛⎫=⎪ΦΦ⎝⎭序号 升程 （mm ） 升程运 动角（ ） 升程运 动规律 升程 许用压力角 （） 回程运动角 （） 回程运 动规律 回程许用 压力角 （ ） 远休 止角 （ ） 近休 止角 （ ） 5（改）60120余弦加速度3590等速6545105ω221200cos 2h a πωπϕ⎛⎫=⎪ΦΦ⎝⎭回程运动方程： ()0'01s s h ϕ-Φ+Φ⎡⎤=-⎢⎥Φ⎣⎦1'0h v ω=-Φ 0a =2.2求位移、速度、加速度线图MATLAB 源程序： pi=3.1415926; c=pi/180; h=60; f0=120; fs=45; f01=90; fs1=105;f=0:1:360; %升程阶段 for n=0:f0s(n+1)=h/2*(1-cos(pi/f0*f(n+1)));v(n+1)=pi*h/(2*f0*c)*sin(pi/f0*f(n+1));a(n+1)=pi^2*h/(2*f0^2*c^2)*cos(pi/f0*f(n+1)); end%远休止阶段 for n=f0:f0+fs s(n+1)=60; v(n+1)=0; a(n+1)=0; endfor n=f0+fs:f0+fs+f01 %回程阶段 s(n+1)=h*(1-(f(n+1)-(f0+fs))/f01); v(n+1)=-h/(f01*c); a(n+1)=0; end;for n=f0+fs+f01:360; %近休止阶段 s(n+1)=0; v(n+1)=0; a(n+1)=0; endfigure(1);plot(f,s,'r');xlabel('\phi/\circ');ylabel('s/mm');grid on;title('推杆位移线图')figure(2);plot(f,v,'r');xlabel('\phi/\circ');ylabel('v/\ommiga');grid on;title('推杆速度线图')figure(3);plot(f,a,'r');xlabel('\phi/\circ');ylabel('a/\ommiga');grid on;title('推杆加速度线图')2.3位移、速度、加速度线图：3 凸轮机构的dss d ϕ-线图，确定基圆半径和偏心距： 3.1理论分析：机构压力角α应按下式计算： 0-tan =+ds d es sϕα作出-ds d s ϕ曲线，因凸轮机构压力角0-tan =+ds d e s sϕα，式中，2200=-s r e ，右侧为升程，作与s 轴夹角等于升程许用压力角35的切界线D d t t ，则在直线上或其右下方取凸轮轴心时，可使[]0-=arctan+ds d es sϕαα≤，同理左侧回程，作与s 轴夹角等于回程许用压力角65的切界线D 'd 't t ，则在直线上或其左下方取凸轮轴心时，使[]0-=arctan+ds d es sϕαα≤。
哈工大机械原理大作业 凸轮资料
机械原理大作业课程名称：机械原理设计题目：直动从动件盘形凸轮机构院系：机电学院班级：完成者：学号：指导教师：设计时间：哈尔滨工业大学题目：如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
凸轮运动分为五个阶段 1.升程阶段0~5034500010156s h ϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦ 23410000306030h v ωϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦ 23212000060180120h a ωϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦2.远休止50~1503.回程等加速150~195()20'202sh s h ϕ=--Φ-ΦΦ()10'204s h s ωϕ=--Φ-ΦΦ21'204h v ω=Φ4.回程等减速195~240()2'00'202s h s ϕ=Φ+Φ+Φ-Φ()'100'24s h v ωϕ=-Φ+Φ+Φ-Φ21'204h v ω=-Φ5.近休止240~360其中，Φ 推程运动角sΦ 远休止角 '0Φ 回程运动角使用Matlab实现1.从动件位移%用fi1,fi2,fi3,fi4，fi5代替转角h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;s1=h*(10*(fi1/fio1).^3-15*(fi1/fio1).^4+6*(fi1/fio1).^5); hold onplot(180*fi1/pi,s1);fi2=fio1:0.01:fio1+fis;s2=h;hold onplot(180*fi2/pi,s2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;s3=h-2*h/fio2^2*(fi3-fio1-fis).^2;hold onplot(180*fi3/pi,s3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;s4=2*h/fio2^2*(fio1+fis+fio2-fi4).^2;hold onplot(180*fi4/pi,s4);fi5=fio1+fis+fio2:0.001:2*pi;s5=0;hold onplot(180*fi5/pi,s5);title('位移');xlabel('φ/度'),ylabel('s/mm');grid onhold off2.从动件速度%用fi1,fi2,fi3,fi4，fi5代替转角%设角速度为1rad/sw=1;h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;v1=h*w/fio1*(30*(fi1/fio1).^2-60*(fi1/fio1).^3+30*(fi1/fio1).^4); hold onplot(180*fi1/pi,v1);fi2=fio1:0.01:fio1+fis;v2=0;hold onplot(180*fi2/pi,v2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;v3=-4*h*w/fio2^2*(fi3-fio1-fis);hold onplot(180*fi3/pi,v3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2; v4=-4*h*w/fio2^2*(fio1+fis+fio2-fi4); hold onplot(180*fi4/pi,v4);fi5=fio1+fis+fio2:0.001:2*pi;v5=0;hold onplot(180*fi5/pi,v5);title('速度');xlabel('φ/度'),ylabel('v/(mm/s)');grid onhold off3.从动件加速度%用fi1,fi2,fi3,fi4，fi5代替转角%设角速度为1rad/sw=1;h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;a1=h*w^2/fio1*(60*(fi1/fio1)-180*(fi1/fio1).^2+120*(fi1/fio1).^3) ;hold onplot(180*fi1/pi,a1);fi2=fio1:0.01:fio1+fis;a2=0;hold onplot(180*fi2/pi,a2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;a3=-4*h*w/fio2^2;hold onplot(180*fi3/pi,a3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;a4=4*h*w/fio2^2;hold onplot(180*fi4/pi,a4);fi5=fio1+fis+fio2:0.001:2*pi;a5=0;hold onplot(180*fi5/pi,a5);title('加速度');xlabel('φ/度'),ylabel('v/(mm/s^2)');grid onhold off4.sdsdϕ-线图%用fi1,fi2,fi3,fi4，fi5代替转角h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;s1=h*(10*(fi1/fio1).^3-15*(fi1/fio1).^4+6*(fi1/fio1).^5);ds1=h*(3*10*(fi1/fio1).^2/fio1-4*15*(fi1/fio1).^3/fio1+5*6*(fi1/f io1).^4/fio1);hold onplot(ds1,s1);fi2=fio1:0.01:fio1+fis;s2=h+0*fi2;ds2=0*fi2;hold onplot(ds2,s2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;s3=h-2*h/fio2^2*(fi3-fio1-fis).^2;ds3=-2*h/fio2^2*2*(fi3-fio1-fis);hold onplot(ds3,s3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2; s4=2*h/fio2^2*(fio1+fis+fio2-fi4).^2; ds4=-2*h/fio2^2*2*(fio1+fis+fio2-fi4); hold on plot(ds4,s4);fi5=fio1+fis+fio2:0.001:2*pi; s5=0+0*fi5; ds5=0*fi5; hold on plot(ds5,s5); title('ds/d φ-s');xlabel('ds/d φ(mm/rad)'),ylabel('s(mm)'); grid on hold off5.凸轮轴心位置的确定凸轮压力角的正切值s s e d ds +-=0/tan ϕα，右侧为升程，作与s 轴夹6π角等于升程许用压力角的切界线t t d D ，则在直线上或其下方取凸轮轴心时，可使[]αα≤，同理右侧回程，作与s 轴夹角等于回程许用压力角3π的切界线''t t d D ，则在直线上或其下方取凸轮轴心时，可使[]αα≤。
哈工大—机械原理凸轮大作业
一、题目要求及机构运动简图如图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。
哈工大机械原理大作业凸轮
1 / 17机械原理大作业课程名称：机械原理 设计题目：直动从动件盘形凸轮机构 院 系： 机电学院 班 级： 完 成 者： 学 号： 指导教师： 设计时间：哈尔滨工业大学题目：如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
凸轮运动分为五个阶段1.升程阶段0~5034500010156s h ϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦23410000306030h v ωϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦23212000060180120h a ωϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦2.远休止50~1503.回程等加速150~195()20'202sh s h ϕ=--Φ-ΦΦ()10'204s h s ωϕ=--Φ-ΦΦ21'204h v ω=Φ4.回程等减速195~240()2'00'202s h s ϕ=Φ+Φ+Φ-Φ()'100'204s h v ωϕ=-Φ+Φ+Φ-Φ21'204h v ω=-Φ5.近休止 240~360其中，0Φ 推程运动角 s Φ 远休止角 '0Φ 回程运动角使用Matlab 实现1.从动件位移%用fi1,fi2,fi3,fi4，fi5代替转角h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;s1=h*(10*(fi1/fio1).^3-15*(fi1/fio1).^4+6*(fi1/fio1).^5); hold onplot(180*fi1/pi,s1);fi2=fio1:0.01:fio1+fis;s2=h;hold onplot(180*fi2/pi,s2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;s3=h-2*h/fio2^2*(fi3-fio1-fis).^2;hold onplot(180*fi3/pi,s3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;s4=2*h/fio2^2*(fio1+fis+fio2-fi4).^2;hold onplot(180*fi4/pi,s4);fi5=fio1+fis+fio2:0.001:2*pi;s5=0;hold onplot(180*fi5/pi,s5);title('位移');xlabel('φ/度'),ylabel('s/mm');grid onhold off2.从动件速度%用fi1,fi2,fi3,fi4，fi5代替转角%设角速度为1rad/sw=1;h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;v1=h*w/fio1*(30*(fi1/fio1).^2-60*(fi1/fio1).^3+30*(fi1/fio1).^4); hold onplot(180*fi1/pi,v1);fi2=fio1:0.01:fio1+fis;v2=0;hold onplot(180*fi2/pi,v2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;v3=-4*h*w/fio2^2*(fi3-fio1-fis);hold onplot(180*fi3/pi,v3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;v4=-4*h*w/fio2^2*(fio1+fis+fio2-fi4);hold onplot(180*fi4/pi,v4);fi5=fio1+fis+fio2:0.001:2*pi;v5=0;hold onplot(180*fi5/pi,v5);title('速度');xlabel('φ/度'),ylabel('v/(mm/s)');grid onhold off3.从动件加速度%用fi1,fi2,fi3,fi4，fi5代替转角%设角速度为1rad/sw=1;h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;a1=h*w^2/fio1*(60*(fi1/fio1)-180*(fi1/fio1).^2+120*(fi1/fio1).^3); hold onplot(180*fi1/pi,a1);fi2=fio1:0.01:fio1+fis;a2=0;hold onplot(180*fi2/pi,a2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;a3=-4*h*w/fio2^2;hold onplot(180*fi3/pi,a3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;a4=4*h*w/fio2^2;hold onplot(180*fi4/pi,a4);fi5=fio1+fis+fio2:0.001:2*pi;a5=0;hold onplot(180*fi5/pi,a5);title('加速度');xlabel('φ/度'),ylabel('v/(mm/s^2)');grid onhold off4.sdsdϕ-线图%用fi1,fi2,fi3,fi4，fi5代替转角h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;s1=h*(10*(fi1/fio1).^3-15*(fi1/fio1).^4+6*(fi1/fio1).^5);ds1=h*(3*10*(fi1/fio1).^2/fio1-4*15*(fi1/fio1).^3/fio1+5*6*(fi1/fio1).^4/fio1);hold onplot(ds1,s1);fi2=fio1:0.01:fio1+fis;s2=h+0*fi2;ds2=0*fi2;hold onplot(ds2,s2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;s3=h-2*h/fio2^2*(fi3-fio1-fis).^2; ds3=-2*h/fio2^2*2*(fi3-fio1-fis); hold onplot(ds3,s3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2; s4=2*h/fio2^2*(fio1+fis+fio2-fi4).^2; ds4=-2*h/fio2^2*2*(fio1+fis+fio2-fi4); hold onplot(ds4,s4);fi5=fio1+fis+fio2:0.001:2*pi; s5=0+0*fi5; ds5=0*fi5; hold onplot(ds5,s5);title('ds/d φ-s');xlabel('ds/d φ(mm/rad)'),ylabel('s(mm)'); grid on hold off5.凸轮轴心位置的确定凸轮压力角的正切值s s ed ds +-=0/tan ϕα，右侧为升程，作与s 轴夹6π角等于升程许用压力角的切界线t t dD ，则在直线上或其下方取凸轮轴心时，可使[]αα≤，同理右侧回程，作与s 轴夹角等于回程许用压力角3π的切界线''t t d D ，则在直线上或其下方取凸轮轴心时，可使[]αα≤。
机械原理大作业-凸轮机构
二、凸轮机构一、运动分析凸轮的运动分为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。
哈工大机械原理大作业二凸轮机构
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。
哈工大机械原理大作业第25题凸轮设计
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院指导教师：林琳福利设计时间：2014.06工业大学凸轮机构设计1．设计题目(1)凸轮机构运动简图：(2)凸轮机构的原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角25 130mm 90°等加等减速40°80°等速70°60°130°2．确定凸轮推杆升程、回程运动方程并绘制推杆位移,速度,加速度线图：(1) 凸轮推杆升程，回程运动方程如下：A.推杆升程方程B.推杆回程方程(2)推杆位移,速度,加速度线图如下：A.推杆位移线图(使用matlab画图，程序详见附录1)B.推杆速度线图(使用matlab画图，程序详见附录2)C.推杆加速度线图(使用matlab画图，程序详见附录3)3．凸轮机构的-s线图，并依次确定凸轮的基圆半径和偏距.(1) 凸轮机构的-s线图：(使用matlab画图，程序详见附录4)(2)确定凸轮的基圆半径和偏距：以ds/df-s(f)图为基础，可分别作出三条限制线（推程许用压力角的切界限Dt d t ,回程许用压力角的限制线Dt 'dt'，起始点压力角许用线Bd'')，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。
得图如下：(使用matlab画图，程序详见附录6)得最小基圆对应的坐标位置大约为（55.28，-65.88）经计算取偏距e=55mm，r0=90mm.4．确定滚子半径及绘制凸轮理论轮廓曲线和实际轮廓曲线.为求滚子许用半径，须确定最小曲率半径，以防止凸轮工作轮廓出现尖点或出现相交包络线，确定最小曲率半径数学模型如下：)/)(/()/)(/(])/()/[(22222/322ϕϕϕϕϕϕρd x d d dy d y d d dx d dy d dx -+=其中：ϕϕϕϕcos )(sin ])/[(/0s s e d ds d dx ++-= ϕϕϕϕsin )(cos ])/[(/0s s e d ds d dy +--=ϕϕϕϕϕsin ])/[(cos ])/(2[/02222s s d s d e d ds d x d --+-= ϕϕϕϕϕcos ])/[(sin ])/(2[/02222s s d s d e d ds d y d --+--=利用上式可求的最小曲率半径后可确定实际廓线。
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)。
哈工大机械原理大作业-凸轮机构设计(第题)————————————————————————————————作者：————————————————————————————————日期：机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电学院班级：1208103完成者：xxxxxxx学号：11208103xx指导教师：林琳设计时间：2014.5.2哈尔滨工业大学凸轮机构设计一、设计题目如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
序号 升程（mm ） 升程运动角（°） 升程运动规律 升程许用压力角（°） 回程运动角（°） 回程运动规律 回程许用压力角（°）远休止角（°） 近休止角（°）3 50 150 正弦加速度 30 100 余弦加速度60 30 80二、凸轮推杆升程、回程运动方程及其线图1 、凸轮推杆升程运动方程(650πϕ≤≤) 升程采用正弦加速度运动规律，故将已知条件mm h 50=，650π=Φ带入正弦加速度运动规律的升程段方程式中得：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512sin 215650ϕππϕS ； ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512cos 1601ππωv ； ω⎪⎭⎫ ⎝⎛=512sin 14421ϕπωa ； 2、凸轮推杆推程远休止角运动方程（πϕπ≤≤65） mm h s 50==；0==a v ；3、凸轮推杆回程运动方程（914πϕπ≤≤） 回程采用余弦加速度运动规律，故将已知条件mm h 50=，95'0π=Φ，6s π=Φ带入余弦加速度运动规律的回程段方程式中得：⎥⎦⎤⎢⎣⎡-+=)(59cos 125πϕs ； ()πϕω--=59sin451v ； ()πϕω-=59cos 81-a 21；4、凸轮推杆回程近休止角运动方程（πϕπ2914≤≤） 0===a v s ；5、凸轮推杆位移、速度、加速度线图根据以上所列的运动方程，利用matlab 绘制出位移、速度、加速度线图。
①位移线图编程如下： %用t 代替转角t=0:0.01:5*pi/6;s=50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5));hold onplot(t,s);t=5*pi/6:0.01:pi;s=50;hold onplot(t,s);t=pi:0.01:14*pi/9;s=25*(1+cos(9*(t-pi)/5));hold onplot(t,s);t=14*pi/9:0.001:2*pi;s=0;hold onplot(t,s),xlabel('φ/rad'),ylabel('s/mm'); grid onhold off所得图像为：②速度线图编程如下：%用t代替转角,设凸轮转动角速度为1t=0:0.01:5*pi/6;v=60/pi*(1-cos((12*t)/5));hold onplot(t,v);t=5*pi/6:0.01:pi;v=0;hold onplot(t,v);t=pi:0.01:14*pi/9;v=-45*sin(9*(t-pi)/5);hold onplot(t,v);t=14*pi/9:0.01:2*pi;v=0;hold onplot(t,v),xlabel('φ(rad)'),ylabel('v(mm/s)'); grid onhold off所得图像为：③加速度线图利用matlab编程如下：%用t代替转角,设凸轮转动角速度为1t=0:0.01:5*pi/6;a=144/pi*sin(12*t/5);hold onplot(t,a);t=5*pi/6:0.01:pi;a=0;hold onplot(t,a);t=pi:0.01:14*pi/9;a=-81*cos(9*(t-pi)/5);hold onplot(t,a);t=14*pi/9:0.01:2*pi;a=0;hold onplot(t,a),xlabel('φ(rad)'),ylabel('a(mm/s^2)');grid onhold off 所得图形：三、绘制s d ds -ϕ线图 根据运动方程求得：()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧≤≤≤≤-≤≤≤≤--=πϕππϕππϕπϕππϕπππϕ2914.0914,59sin 4565,0650),512cos 6060(d ds 利用matlab 编程：%用t 代替φ,a 代替ds/d φ,t=0:0.01: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 onplot(a,s);t=5*pi/6:0.01:pi;a=0;s=50;hold onplot(a,s);t=pi:0.01:14*pi/9;a=45*sin(9*(t-pi)/5);s=25*(1+cos(9*(t-pi)/5));hold onplot(a,s);t=14*pi/9:0.01:2*pi;a=0;s=0;hold onplot(a,s),title('ds/dφ-s'),xlabel('ds/d φ(mm/rad)'),ylabel('s(mm)');grid onhold off得s d ds -ϕ图： 凸轮压力角的正切值ss e d ds +-=0/tan ϕα，左侧为升程，作与s 轴夹6π角等于升程许用压力角的切界线t t d D ，则在直线上或其左下方取凸轮轴心时，可使[]αα≤，同理右侧回程，作与s 轴夹角等于回程许用压力角3π的切界线''t t d D ，则在直线上或其右下方取凸轮轴心时，可使[]αα≤。
在升程起始点，有S=0,0/=ϕd ds ，为保证此时[]αα≤，作直线00d B 与纵坐标夹角为6π，凸轮轴心只能在其线上或右下方选取。
三条限制线围成的下方阴影角区域为满足[]αα≤的凸轮轴心的公共许用区域。
编程求得公共许用区域：t0=0:0.001:5*pi/6;t1=5*pi/6:0.001:pi;t2=pi:0.001:14*pi/9;t3=14*pi/9:0.001:2*pi;a0=-(60/pi-60/pi*cos(12*t0/5));a1=0;a2=45*sin(9*(t2-pi)/5);a3=0;s0=50*((6*t0)/(5*pi)-1/(2*pi)*sin(12*t0/5));s1=50;s2=25*(1+cos(9*(t2-pi)/5));s3=0;k1=-tan(pi/3);k2=tan(pi/6);x=-50:0.001:40;b=s0-k1.*a0;m=s2-k2.*a2;g=min(b);h=min(m);y1=k1*x+g;y2=k2*x+h;plot(x,y1,'c',x,y2,'b');hold onx=-50:0.001:0;y3=tan(pi/2-pi/6)*x;plot(x,y3,'m');hold onplot(a0,s0,'b',a1,s1,'g',a2,s2,'r',a3,s3,'c'),title('凸轮轴心位置确定'),xlabel('x(mm)'),ylabel('y(mm)');grid on得到：则可取轴心为x=-10mm,,y=-40mm得：e=10mm ，mm 400=S ，mm S e o 2311.41r 202=+=四、绘制凸轮理论轮廓线以及基圆、偏距圆凸轮理论轮廓方程：ϕϕsin cos )(0e s s x ++-=；ϕϕcos sin )(0e s s y +++=；利用matlab 编程：t=0:0.001:5*pi/6;x=-(40+50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5))).*cos(t)+10*sin(t); y=(40+50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5))).*sin(t)+10*cos(t); hold onplot(x,y);t=5*pi/6:0.001:pi;x=-(40+50).*cos(t)+10*sin(t);y=(40+50).*sin(t)+10*cos(t);hold onplot(x,y);t=pi:0.001:14*pi/9;x=-(40+25*(1+cos(9*(t-pi)/5))).*cos(t)+10*sin(t);y=+(40+25*(1+cos(9*(t-pi)/5))).*sin(t)+10*cos(t);hold onplot(x,y);t=14*pi/9:0.001:2*pi;x=-(40).*cos(t)+10*sin(t);y=+(40).*sin(t)+10*cos(t);hold onplot(x,y);%绘制基园,用绿色线条表示t=0:0.001:2*pi;x=41.2310*cos(t);y=41.2310*sin(t);hold onplot(x,y,'g');%绘制偏距圆，用红色线条表示t=0:0.001:2*pi;x=10*cos(t);y=10*sin(t);hold onplot(x,y,'r'),title('凸轮理论廓线、基圆以及偏距圆'),xlabel('x(mm)'),ylabel('y(mm)');grid onhold off得到：五、确定滚子半径根据曲率半径公式：)/)(/()/)(/(])/()/[(22222/322ϕϕϕϕϕϕρd x d d dy d y d d dx d dy d dx -+= 利用matlab 编程如下：%t 代替转角φ，p 表示曲率半径t=0:0.001:5*pi/6;x=-(40+50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5))).*cos(t)+10*sin(t); y=(40+50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5))).*sin(t)+10*cos(t); dx1=diff(x);dy1=diff(y);dx11=diff(x,2);dy11=diff(y,2);p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11));hold onplot(t,p);t=5*pi/6:0.001:pi;x=-(40+50).*cos(t)+10*sin(t);y=(40+50).*sin(t)+10*cos(t);dx1=diff(x);dy1=diff(y);dx11=diff(x,2);dy11=diff(y,2);p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11));hold onplot(t,p);t=pi:0.001:14*pi/9;x=-(40+25*(1+cos(9*(t-pi)/5))).*cos(t)+10*sin(t);y=+(40+25*(1+cos(9*(t-pi)/5))).*sin(t)+10*cos(t);dx1=diff(x);dy1=diff(y);dx11=diff(x,2);dy11=diff(y,2);p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11));hold onplot(t,p);t=14*pi/9:0.001:2*pi;x=-(40).*cos(t)+10*sin(t);y=+(40).*sin(t)+10*cos(t);dx1=diff(x);dy1=diff(y);dx11=diff(x,2);dy11=diff(y,2);p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11));hold onplot(t,p);hold offtitle('曲率半径ρ'),xlabel('φ(rad)'),ylabel('ρ(mm)');grid on得到：由图可知：m m 2311.41min =ρ，即基圆半径，又因为m m 5~3,min =∆∆-<其中ρr r ，综上，可取滚子半径为10mm 。