转动导杆机构的运动分析

% 转动导杆机构运动分析
% 已知条件
n=200;omega_1=pi*n/30; % 主动件导杆转速和角速度(rad/s) lambda=0.5; % 机构尺度系数
alpha_p=45; % 机构许用压力角
du=180/pi;hd=pi/180; % 角度与弧度转换系数
disp ' ******** 已知条件 ********'
fprintf (' 导杆转速 n = %3.4f r/min \n',n)
fprintf (' 导杆角速度 omega_1 = %3.4f rad/s \n',omega_1)
fprintf (' 机构尺度系数 lambda = %3.4f \n',lambda)
fprintf ('机构许用压力角 alpha_p = %3.4f °\n',alpha_p)
% 1-计算机构运动参数
for phi_1=1:360
phi_3(phi_1)=phi_1-asin(lambda*sin(phi_1*hd))*du; % 从动件曲柄角位移 omega_3_z=lambda*cos(phi_1*hd);
omega_3_m=sqrt(1-lambda^2*sin(phi_1*hd)^2);
omega_3(phi_1)=omega_1*(1-omega_3_z/omega_3_m); % 从动件曲柄角速度 i_31(phi_1)=omega_3(phi_1)/omega_1; % 机构传动比 alpha(phi_1)=asin(lambda*sin(phi_1*hd))*du; % 机构压力角
end
disp ' '
disp ' *** 转动导杆机构的运动参数 ***'
disp ' 导杆转角 曲柄角位移 曲柄角速度'
disp ' (°) (°) (rad/s) '
for phi_1=10:10:360
ydcs=[phi_1 phi_3(phi_1) omega_3(phi_1)];
disp(ydcs)
end
% 2-绘制机构的运动线图
figure(1);
subplot(2,1,1) % 曲柄角位移线图
phi_1=1:360;plot(phi_1,phi_3);
xlabel ('导杆转角 \it \phi_1 / \rm( °)')
ylabel ('\it \phi_3 / \rm(°)')
title('从动件曲柄角位移线图');
subplot(2,1,2) % 曲柄角速度线图
phi_1=1:360;plot(phi_1,omega_3);
xlabel ('导杆转角 \it \phi_1 / \rm( °)')
ylabel ('\it \omega_3 / \rm(rad/s)')
title('从动件曲柄角速度线图');
% 3-绘制机构的压力角线图
figure(2);
phi_1=1:360;plot(phi_1,alpha);
xlabel ('导杆转角 \it \phi_1 / \rm( °)');