哈工大机械原理大作业 凸轮机构 32 (3)

ds4 = diff(s4,phi4); ds5 = diff(s5,phi5); ds6 = diff(s6,phi6);

alpha1 = atan(abs(ds1 - e)/(s0 + s1)); alpha2 = atan(abs(ds2 - e)/(s0 + s2)); alpha3 = atan(abs(ds3 - e)/(s0 + s3)); alpha4 = atan(abs(ds4 - e)/(s0 + s4)); alpha5 = atan(abs(ds5 - e)/(s0 + s5)); alpha6 = atan(abs(ds6 - e)/(s0 + s6));

dx11 = diff(dx1,phi1); dy11 = diff(dy1,phi1);

dx22 = diff(dx2,phi2); dy22 = diff(dy2,phi2); dx33 = diff(dx3,phi3); dy33 = diff(dy3,phi3); dx44 = diff(dx4,phi4); dy44 = diff(dy4,phi4); dx55 = diff(dx5,phi5); dy55 = diff(dy5,phi5); dx66 = diff(dx6,phi6); dy66 = diff(dy6,phi6);

rou1 = abs(((dx1^2 + dy1^2)^1.5)/(dx1*dy11 - dx11*dy1)); rou2 = abs(((dx2^2 + dy2^2)^1.5)/(dx2*dy22 - dx22*dy2)); rou3 = abs(((dx3^2 + dy3^2)^1.5)/(dx3*dy33 - dx33*dy3)); rou4 = abs(((dx4^2 + dy4^2)^1.5)/(dx4*dy44 - dx44*dy4)); rou5 = abs(((dx5^2 + dy5^2)^1.5)/(dx5*dy55 - dx55*dy5)); rou6 = abs(((dx6^2 + dy6^2)^1.5)/(dx6*dy66 - dx66*dy6));

phi1 = 0:pi/1000:PHI0;

phi2 = PHI0:pi/1000:PHI0+PHIs;

phi3 = (PHI0 + PHIs):pi/1000:(PHI0 + PHIs + PHI00/n);

phi4 = (PHI0 + PHIs + PHI00/n):pi/1000:(PHI0 + PHIs + PHI00*(n - 1)/n); phi5 = (PHI0 + PHIs + PHI00*(n - 1)/n):pi/1000:(PHI0 +PHIs + PHI00); phi6 = (PHI0 + PHIs + PHI00):pi/1000:2*pi;

plot(phi1,subs(alpha1),

phi2,subs(alpha2),phi3,subs(alpha3),phi4,subs(alpha4),phi5,subs(alpha5),phi6,subs(alpha6))

plot(phi1,subs(rou1),

phi2,subs(rou2),phi3,subs(rou3),phi4,subs(rou4),phi5,subs(rou5),phi6,subs(rou6))

压力角线图

曲率半径线图

三、确定滚子半径，绘制凸轮理论轮廓和实际轮廓

凸轮逆时针转动，从动件导路偏置于凸轮轴心右侧，凸轮理论轮廓线方程为：??x=(s0+s)sin?+ecos?，(0???2?)

y=(s+s)cos?-esin?0?逐段代入位移表达式则可求凸轮理论轮廓线方程，用Matlab编写程序绘制如图所示。

逐段求出最小值得，，取滚子半径，以理论

轮廓上各点为圆心、以滚子半径为半径作出的滚子圆族包络线即为滚子从动件盘形凸轮的工作轮廓，其方程为：

dyd??X=x+rr?22(dxd?)+(dyd?)?，0???2? ?dxd??Y=y-rr?(dxd?)2+(dyd?)2?MATLAB代码

clear,clc r0 = 90; e = 40;

s0 = sqrt(r0^2 - e^2); h = 50;

PHI0 = pi/3; n = 3;

PHI00 =4*pi/9; PHIs = 5*pi/9; rr = 10;

syms phi1phi2phi3phi4phi5phi6

s1 = h*(10*(phi1/PHI0).^3 - 15*(phi1/PHI0).^4 + 6*(phi1/PHI0).^5);

s2 = h*(10 - 15 + 6);

s3 = h - (2*h/(4 + (n - 2)*pi))*(1 - cos(n*pi*(phi3 - PHI0 - PHIs)/(2*PHI00)));

s4 = h - (h/(4 + (n - 2)*pi))*(n*pi*(phi4 - PHI0 - PHIs)/PHI00 - pi + 2);

s5 = (2*h/(4 + (n - 2)*pi))*(1 + cos(n*pi*(phi5 - PHI0 - PHIs)/(2*PHI00) - (n - 2)*pi/2));

s6 = 0;

x1 = -(s0 + s1)*sin(phi1) - e*cos(phi1); %àí????àa y1 = (s0 + s1)*cos(phi1) - e*sin(phi1);

x2 = -(s0 + s2)*sin(phi2) - e*cos(phi2); %àí????àa y2 = (s0 + s2)*cos(phi2) - e*sin(phi2);

x3 = -(s0 + s3)*sin(phi3) - e*cos(phi3); %àí????àa y3 = (s0 + s3)*cos(phi3) - e*sin(phi3);

x4 = -(s0 + s4)*sin(phi4) - e*cos(phi4); %àí????àa y4 = (s0 + s4)*cos(phi4) - e*sin(phi4);

x5 = -(s0 + s5)*sin(phi5) - e*cos(phi5); %àí????àa y5 = (s0 + s5)*cos(phi5) - e*sin(phi5);

x6 = -(s0 + s6)*sin(phi6) - e*cos(phi6); %àí????àa y6 = (s0 + s6)*cos(phi6) - e*sin(phi6);

dy1 = diff(y1,phi1); dx1 = diff(x1,phi1);

dy2 = diff(y2,phi2); dx2 = diff(x2,phi2);

dy3 = diff(y3,phi3); dx3 = diff(x3,phi3);

dy4 = diff(y4,phi4); dx4 = diff(x4,phi4);

dy5 = diff(y5,phi5); dx5 = diff(x5,phi5);

dy6 = diff(y6,phi6); dx6 = diff(x6,phi6);

X1 = x1 - rr*dy1/sqrt(dx1^2 + dy1^2); %êμ?ê??àa Y1 = y1 + rr*dx1/sqrt(dx1^2 + dy1^2);

X2 = x2 - rr*dy2/sqrt(dx2^2 + dy2^2); %êμ?ê??àa Y2 = y2 + rr*dx2/sqrt(dx2^2 + dy2^2);

X3 = x3 - rr*dy3/sqrt(dx3^2 + dy3^2); %êμ?ê??àa Y3 = y3 + rr*dx3/sqrt(dx3^2 + dy3^2);

X4 = x4 - rr*dy4/sqrt(dx4^2 + dy4^2); %êμ?ê??àa

Y4 = y4 + rr*dx4/sqrt(dx4^2 + dy4^2);

X5 = x5 - rr*dy5/sqrt(dx5^2 + dy5^2); %êμ?ê??àa Y5 = y5 + rr*dx5/sqrt(dx5^2 + dy5^2);

X6 = x6 - rr*dy6/sqrt(dx6^2 + dy6^2); Y6 = y6 + rr*dx6/sqrt(dx6^2 + dy6^2);

phi1 = 0:pi/1000:PHI0;

phi2 = PHI0:pi/1000:PHI0+PHIs;

phi3 = (PHI0 + PHIs):pi/1000:(PHI0 + PHIs + PHI00/n);

phi4 = (PHI0 + PHIs + PHI00/n):pi/1000:(PHI0 + PHIs + PHI00*(n - 1)/n); phi5 = (PHI0 + PHIs + PHI00*(n - 1)/n):pi/1000:(PHI0 +PHIs + PHI00); phi6 = (PHI0 + PHIs + PHI00):pi/1000:2*pi;

plot(subs(X1),subs(Y1),subs(X2),subs(Y2),subs(X3),subs(Y3),subs(X4),subs(Y4),subs(X5),subs(Y5),subs(X6),subs(Y6)) hold on

plot(subs(x1),subs(y1),subs(x2),subs(y2),subs(x3),subs(y3),subs(x4),subs(y4),subs(x5),subs(y5),subs(x6),subs(y6))hold off axis square

凸轮轮廓图

搜索“diyifanwen.net”或“第一范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，第一范文网，提供最新经管营销哈工大机械原理大作业 凸轮机构 32 (3)全文阅读和word下载服务。