机械原理大作业头刨床.docVIP

机械原理 大作业 利用解析法对牛头刨床各构件运动分析 姓名： 班级： 学号： 如图所示为一牛头刨床机构简图，各构件的尺寸如图：L1=125 mm、L3=600 mm、L4=150mm,利用解释法求出图中各构件的方位角、角速度和角加速度及E点的位移、速度、加速度析。 解：利用MATLAB编程： 输出为θ3、θ4的值的程序： x1=linspace(0,2*pi,360); l1=125; l3=600; l4=150; l6=275; y=575; w1=1; x3=atan((l6+l1*sin(x1))./(l1*cos(x1))); for i = 1 : 360; if x3(i)0 x3(i)=x3(i)+pi; end end w3=w1*l1*(l1+6*sin(x1))./(l6.^2+l1.^2+2*l1*l6*sin(x1)); A3=(l6.^2-l1.^2)*l6*l1*w1.^2*cos(x1)./(l6.^2+l1.^2+2*l1*l6*sin(x1)).^2; x4= pi-asin((y-l3*sin(x3))./l4); Se=l3*cos(x3)+l4*cos(x4); w4=(-1)*w3.*l3.*cos(x3)./(l4.*cos(x4)); Ve=(-1)*w3.*l3.*sin(x3-x4)./cos(x4); A4=(w3.^2.*l3.*sin(x3)+w4.^2.*l4.*sin(x4)-A3.*l3.*cos(x3))./(l4.*cos(x4)); ae=(-1)*(A3.*l3.*sin(x3-x4)+w3.^2.*l3.*cos(x3-x4)-w4.^2.*l4)./cos(x4); plotyy(180*x1./pi,[180*x3./pi;180*x4./pi],180*x1./pi,Se/1000); legend(x3,x4,Se); xlabel(x1/(°)),ylabel(x3、x4(°)); title(x3、x4、Se曲线); grid on; plotedit on; 运行结果： 运行的图形： 关于W3和W4的输出： x1=linspace(0,2*pi,360); l1=125; l3=600; l4=150; l6=275; y=575; w1=1; x3=atan((l6+l1*sin(x1))./(l1*cos(x1))); for i = 1 : 360; if x3(i)0 x3(i)=x3(i)+pi; end end w3=w1*l1*(l1+6*sin(x1))./(l6.^2+l1.^2+2*l1*l6*sin(x1)); A3=(l6.^2-l1.^2)*l6*l1*w1.^2*cos(x1)./(l6.^2+l1.^2+2*l1*l6*sin(x1)).^2; x4= pi-asin((y-l3*sin(x3))./l4); Se=l3*cos(x3)+l4*cos(x4); w4=(-1)*w3.*l3.*cos(x3)./(l4.*cos(x4)); Ve=(-1)*w3.*l3.*sin(x3-x4)./cos(x4); A4=(w3.^2.*l3.*sin(x3)+w4.^2.*l4.*sin(x4)-A3.*l3.*cos(x3))./(l4.*cos(x4)); ae=(-1)*(A3.*l3.*sin(x3-x4)+w3.^2.*l3.*cos(x3-x4)-w4.^2.*l4)./cos(x4); plotyy(180*x1./pi,[w3;w4],180*x1./pi,Ve/1000); legend(w4,w3,Ve); xlabel(x1/(°)),ylabel(w3、w4(rad/s)); title(w3、w4、Ve曲线); grid on; plotedit on; 在MATLABE下运行： 关于a3和a4的程序： x1=linspace(0,2*pi,360); l1=125; l3=600; l4=150; l6=275; y=575; w1=1; x3=atan((l6+l1*sin(x1))./(l1*cos(x1))); for i = 1 : 360; if x3(i)0