for i=1:length(phi)-1
s1(i) = s(i);
end
figure(4);
plot(d,s1);
hold on;
k1 = tan(pi/2 - 30*pi/180);
k2 = -tan(pi/2 - 60*pi/180);
y1min = 0;
y2min = 0;
for i=1:length(phi)-1
dm = diff(m)./diff(phi);
dn = diff(n)./diff(phi);
phi1 = zeros(1,length(phi)-1);
机械原理大作业二
题目直动从动件盘形凸轮机构设计
专业
学号
学生
指导教师
1.
设计如图1所示直动从动件盘形凸轮机构。其原始参数见表1。
表1
行程（mm）
升程运动角（）
升程运动规律
升程许用压力角（）
回程运动角（）
回程运动规律
回程许用压力角（）
远休止角（）
近休止角（）
40
90
正弦加速度
30
50
4-5-6-7多项式
end
end
v = diff(s)./diff(t);
t1 = zeros(1,length(t)-1);
for i=1:length(t)-1
t1(i) = t(i);
end
a = diff(v)./diff(t1);
t2 = zeros(1,length(t)-2);
for i=1:length(t)-2
hold on
m = sqrt(e^2+(s0+s).^2).*cos(atan(s0/e)-phi);
n = sqrt(e^2+(s0+s).^2).*sin(atan(s0/e)-phi);
plot(m,n);
axis equal
axis([-80,120,-110,70]);
(5)曲率半径
%曲率半径线图
60
100
120
推程运动规律方程(正弦加速度):
位移方程
,
速度方程
,
加速度方程
,
回程运动规律方程(4-5-6-7多项式):
位移方程
速度方程
加速度方程
式中
3.
(1)从动件位移、速度、加速度线图
%基本参数
phi0 = 90*pi/180;
if(d(i)>0)
if (-k1*d(i)+s1(i) < y1min)
y1min = -k1*d(i)+s1(i);
end
end
if(d(i)<0)
if (-k2*d(i)+s1(i) < y2min)
y2min = -k2*d(i)+s1(i);
end
end
end
y1 = k1*d+y1min;
title('压力角');
(4)理论轮廓线
%理论廓线
m0= sqrt(s0^2+e^2)*cos(phi);
n0= sqrt(s0^2+e^2)*sin(phi);
figure(6);
plot(m0,n0);
hold on
m1 = e*cos(phi);
n1 = e*sin(phi);
plot(m1,n1);
phis = 100*pi/180;
phi00 = 50*pi/180;
phiss = 120*pi/180;
h = 40;
omega = 10;
%从动件位移、速度、加速度线图
t = 0:pi/5000:pi/5;
phi = omega*t;
s = zeros(1,length(phi));
for i=1:length(phi)
t2(i) = t(i);
end
figure(1);
plot(omega*t,s);
title('位移');
set(gca,'XTick',0:pi/2:2*pi);
set(gca,'XTicklabel',{'0','π/2','π','3π/2','2π'});
figure(2);
plot(omega*t1,v);
title('速度');
set(gca,'XTick',0:pi/2:2*pi);
set(gca,'XTicklabel',{'0','π/2','π','3π/2','2π'});
figure(3);
plot(omega*t2,a);
title('加速度');
set(gca,'XTick',0:pi/2:2*pi);
if (0 <= phi(i) && phi(i) <= phi0)
s(i) = h*(phi(i)/phi0-(1/(2*pi))*sin(2*pi*phi(i)/phi0));
elseif (phi0 < phi(i) && phi(i) <= phi0+phis)
s(i) = s(i-1);
elseif (phi0+phis < phi(i) && phi(i) <= phi0+phis+phi00)
(3)压力角线图
%给定基圆半径和偏距
e=20;
s0=60;
%压力角线图
alpha = atan(abs(d-e)./(s1+s0));
figure(5);
plot(omega*t1,alpha*180/pi);
set(gca,'XTick',0:pi/2:2*pi);
set(gca,'XTicklabel',{'0','π/2','π','3π/2','2π'});
y2 = k2*d+y2min;
y0 = -k1*d;
plot(d,y1);
hold on;
plot(d,y2);
hold on;
plot(d,y0);
axis([-120,60,-100,100]);
title('ds/dφ-s线图');
grid on;
三条直线下方的公共区域中，取x0=20,y0=-60,则凸轮基圆半径e=20mm,s0=60mm,r0= 。
set(gca,'XTicklabel',{'0','π/2','π','3π/2','2π'});
(2)ds/dφ-s线图及基圆半径和偏距的确定
%ds/dφ-s线图
ds=diff(s);
dphi=diff(phi);
d = ds./dphi;
s1 = zeros(1,length(phi)-1);
s(i) = h*(1-35*((phi(i)-phi0-phis)/phi00)^4+84*((phi(i)-phi0-phis)/phi00)^5-70*((phi(i)-phi0-phis)/phi00)^6+20*((phi(i)-phi0-phis)/phi00)^7);
else
s(i) = s(i-1);