偏置直动滚子推杆盘形凸轮Matlab编程(附录程序).doc

《偏置直动滚子推杆盘形凸轮Matlab编程(附录程序).doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、机械原理大作业学院：机械与电子信息学院授课老师：曾小慧姓名：张京学号：20131004547日期：2015-5-2313机械原理大作业目录1.求轮廓曲线推程阶段远休止阶段回程阶段近休止阶段Matlab程序设计轮廓图形2.求工作廓线推程阶段远休止阶段回程阶段近休止阶段Matlab程序设计轮廓图形３．求解最大压力角　压力角公式　MATLAB程序设计　根据MATLAB程序作图可得出其压力角与角度的关系并分析　失真情况分析４．附录Matlab程序13机械原理大作业凸轮轮廓9-14试设计偏置直动滚子推杆盘形凸轮机构的理论轮廓曲线和工作廓线。已知凸轮轴置于推杆轴线右侧，偏距e=20m

2、m，基圆半径ｒ。=50mm，滚子半径rr=10mm。凸轮以等角速度沿顺时针方向回转，在凸轮转过角d1=120º的过程中，推杆按正弦加速度运动规律上升h=50mm；凸轮继续转过d2=30º时，推杆保持不动；其后，凸轮再回转角度d3=60º时，推杆又按余弦加速度运动规律下降至起始位置；凸轮转过一周的其余角度时，推杆又静止不动。解：1.求理论廓线对于偏置直动滚子推杆盘形凸轮机构，凸轮理论廓线上B点（即滚子中心）的直角坐标为（a）式中①推程阶段（）②远休止阶段13机械原理大作业①回程阶段②近休止阶段Matlab程序设计：a1=linspace(0,2*pi/3);%推程阶段的自变

3、量s1=h*(3*a1/2/pi-sin(3*a1)/2/pi);%推杆产生的相应位移x1=-((s0+s1).*sin(a1)+e*cos(a1));%x函数y1=(s0+s1).*cos(a1)-e*sin(a1);%y函数a2=linspace(0,pi/6);%远休止阶段的自变量s2=50;%推杆位移x2=-((s0+s2).*sin(a2+2*pi/3)+e*cos(a2+2*pi/3));%x函数y2=(s0+s2).*cos(a2+2*pi/3)-e*sin(a2+2*pi/3);%y函数a3=linspace(0,pi/3);%回程阶段的自变量s3=h*(

4、1+cos(3*a3))/2;%推杆位移x3=-((s0+s3).*sin(a3+5*pi/6)+e*cos(a3+5*pi/6));%x函数y3=(s0+s3).*cos(a3+5*pi/6)-e*sin(a3+5*pi/6);%y函数a4=linspace(0,5*pi/6);%近休止阶段的自变量s4=0;%推杆位移x4=-((s0+s4).*sin(a4+7*pi/6)+e*cos(a4+7*pi/6));%x函数y4=(s0+s4).*cos(a4+7*pi/6)-e*sin(a4+7*pi/6);%y函数a0=linspace(0,2*pi);%基圆自变量x5=

5、r0*cos(a0);%x函数y5=r0*sin(a0);%y函数轮廓图形通过Matlab软件，编写程序,将以上各相应值代入式(a)计算理论轮廓线上各点的坐标值。在计算时应注意：在推程阶段取，在远休止阶段取，在回程阶段取，在近休止阶段取13机械原理大作业。画出的图形如下图所示2.求工作廓线（b）其中①推程阶段②远休止阶段13机械原理大作业①回程阶段②近休止阶段Matlab程序设计：%工作廓线m1=-(h*3/2/pi*(1-cos(3*a1))-e).*sin(a1)-(s0+s1).*cos(a1);%中间变量dx/d$n1=(h*3/2/pi*(1-cos(3*a1)

6、)-e).*cos(a1)-(s0+s1).*sin(a1);%中间变量dy/d$p1=-m1./sqrt(m1.^2+n1.^2);%sin&q1=n1./sqrt(m1.^2+n1.^2);%cos&x6=x1-r*q1;%x'函数y6=y1-r*p1;%y'函数m2=-(s0+s2).*cos(a2+2*pi/3)+e*sin(a2+2*pi/3);%中间变量dx/d$n2=-(s0+s2).*sin(a2+2*pi/3)-e*cos(a2+2*pi/3);%中间变量dy/d$p2=-m2./sqrt(m2.^2+n2.^2);%sin&q2=n2./sqrt(m2

7、.^2+n2.^2);%cos&x7=x2-r*q2;%x'函数y7=y2-r*p2;%y'函数m3=(h*3/2*sin(3*a3)+e).*sin(a3+5*pi/6)-(s0+s3).*cos(a3+5*pi/6);%中间变量dx/d$n3=-(h*3/2*sin(3*a3)+e).*cos(a3+5*pi/6)-(s0+s3).*sin(a3+5*pi/6);%中间变量dy/d$p3=-m3./sqrt(m3.^2+n3.^2);%sin&q3=n3./sqrt(m3.^2+n3.^2);%cos&x8=x3-r*q3