哈工大机械原理大作业二上传版及答案.doc

哈尔滨工业大学 机械原理大作业 二 利用Matlab进行直动从动件盘形凸轮机构设计说明书 题目：（13） 课程名称：机械原理 学院：外国语学院 姓名：XX 班号：XXX 学号：XXX 一：设计题目 二：建立凸轮轮廓设计计算数学模型 （1）对于绘制从动件位移，速度，加速度线图 凸轮推程（余弦加速度）（w1:50)： s1=h/2*(1-cos(pi/W0*w1)) v1=pi*h*wv/(2*W0)*sin(pi/W0*w1); a1=pi.^2*h*wv.^2/(2*W0.^2)*cos(pi/W0*w1); 凸轮远休止(w2:50-150)： s2=h/2*(1-cos(pi/W0*W0))+(1-w2*eps)*eps; v2=w2*0; a2=w2*0; 凸轮回程（抛-直-抛）(w3,w4,w5:150-240) s3=h-h*n.^2/((2*(n-1))*W01.^2)*(w3-W0-WS).^2; v3=-h*n.^2*wv/((n-1)*W01.^2)*(w3-W0-WS); a3=-h*n.^2*wv.^2/((n-1)*W01.^2)+w3*eps*eps; s4=h-h/(n-1)*(n/W01*(w4-W0-WS)-1/2); v4=-h*n*wv/((n-1)*W01)+w4*eps*eps; a4=w4*0 s5=h*(n^2)/(2*(n-1))*((1-(w5-W0-WS)/W01).^2); v5=-h*n.^2*wv/((n-1)*W01)*(1-(w5-W0-WS)/W01); a5=h*n.^2*wv.^2/((n-1)*W01.^2)+w5*eps*eps; 凸轮近休止：(w6:240-360) w6=W0+WS+W01:0.01:2*pi s6=0*w6; v6=0*w6; a6=0*w6; （2）对于绘制ds/dw-s线图 即为绘制V-S线图。 其中，找到符合范围的轴心，需要作至少三条切界线，它们分别是： y1=tan(55/180*pi)*x1-k1; y2=tan(160/180*pi)*x1-k2; y3=-tan(55/180*pi)*x1 找到直线y1,y2与理论轮廓相切的时候的K的值，此时位于辅助线下方的地方可以取轴心的位置，在这里，我们取e=5.s0=15 （3）对于绘制理论轮廓上的压力角与曲率半径 r0=sqrt(e.^2+s0.^2); 定义a即为压力角 a=arctan(tan(abs(V-e)./(s0+S))) P为曲率半径，则有 P=abs(((dx1.^2+dy1.^2).^(3/2))./(dx1.*ddy1-dy1.*ddx1) 其中： dx1=(s0+S).*cos(W)-e*sin(W)+V.*sin(W); dy1=-(s0+S).*sin(W)-e*cos(W)+V.*cos(W); ddx1=-(s0*sin(W)-V.*cos(W)+S.*sin(W)+e*cos(W)-A.*sin(W)-V.*cos(W)); ddy1=-s0*cos(W)-V.*sin(W)-S.*cos(W)+e*sin(W)+A.*sin(W)+V.*cos(W); （4）对于确定滚子半径，绘制理论以及实际轮廓 滚子半径根据推荐公式为：rr=min(P)-5mm 理论轮廓应当满足方程： X1=(s0+S).*sin(W)+e*cos(W); Y1=(s0+S).*cos(W)-e*sin(W); 实际轮廓应当满足方程： X2=X1-rr*dy1./((dx1.^2+dy1.^2).^(1/2)) Y2=Y1+rr*dx1./((dx1.^2+dy1.^2).^(1/2)) 三：绘制计算流程框图 四：计算程序清单 %位移 wv=0.1 %凸轮角速度 C=180/pi; h=45; W0=50/180*pi; WS=100/180*pi; W01=90/180*pi; n=3 w1=0:0.01:W0 s1=h/2*(1-cos(pi/W0*w1)) w2=W0:0.01:(WS+W0) s2=h/2*(1-cos(pi/W0*W0))+(1-w2*eps)*eps; w3=(WS+W0):0.01:(W0+WS+W01/n); s3=h-h*n.^2/((2*(n-1))*W01.^2)*(w3-W0-WS).^2; w4=(WS+W0+W01/n):0.01:W0+WS+W01*(n-1)/n; s4=h-h/(n-1)*(n/W01*(w4-W0-WS)-1/2); w5=W0+WS+W01*(n-1)/n:0.01:W0+W