曲柄滑块机构分析

任务：根据杆长判断机构类型，模拟仿真曲柄滑块机构运动，并绘制滑块速度及加速度图像。 源代码： clear all clc

l1=4,l2=10,e=0,w1=10; if (e==0) & (l1>l2) % load handel % sound(y,Fs)

disp('当e=0时,l1需小于l2。不满足杆长条件，请重新输入l1,l2,e的值') elseif ((e~=0)&(l1+e>l2)) % load handel % sound(y,Fs)

disp('当e~=0时,l1+e需小于l2。不满足杆长条件，请重新输入l1,l2,e的值') else end

fai=linspace(0,2*pi); for i=1:10:length(fai) JGT2(fai(i),l1,l2,e,w1) end

%%%%%%%%%%%计算滑块速度并绘

制图像 figure (2) fai=0:0.01:2*pi;

v=-l1.*sin(fai).*w1-(l2^2-(l1.*sin(fai)-e).^2).^(-0.5)*(l1^2).*(l1.*sin(fai)-e).*sin(fai).*cos(fai).*w1;

plot(fai,v,'-b')

title('滑块角度-速度图像') xlabel('角度\\phi/°') ylabel('速度/m/s') figure (3) A=- l1.*cos(fai) -

(l1^2.*cos(fai).^2)./(l2^2 - (e -

l1.*sin(fai)).^2).^(1/2) - (l1.*sin(fai).*(e - l1.*sin(fai)))./(l2^2 - (e - l1.*sin(fai)).^2).^(1/2) ... - (l1^2.*cos(fai).^2.*(e - l1.*sin(fai)).^2)./(l2^2 - (e - l1.*sin(fai)).^2).^(3/2);

plot(fai,A,'-b')

title('滑块角度-加速度图像')

xlabel('角度\\phi/°') ylabel('速度/m/s')

结果图：

滑块位移图像642滑块位移/m0-2-4-6-20246角度/°8101214

滑块角度-速度图像80604020速度/m/s0-20-40-60-800123角度/°4567

滑块角度-加速度图像3210速度/m/s-1-2-3-4-5-60123角度/°4567