機(jī)械原理大作業(yè)牛頭刨床

機(jī)械原理大作業(yè)利用解析法對牛頭刨床各構(gòu)件運(yùn)動分析姓名：班級：學(xué)號：如圖所示為一牛頭刨床機(jī)構(gòu)簡圖，各構(gòu)件的尺寸如圖：L1=125mm、L3=600mm、L4=150mm,利用解釋法求出圖中各構(gòu)件的方位角、角速度和角加速度及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)));fori=1:360;ifx3(i)<0x3(i)=x3(i)+pi;endendw3=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曲線');gridon;plotediton;運(yùn)行結(jié)果：運(yùn)行的圖形：關(guān)于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)));fori=1:360;ifx3(i)<0x3(i)=x3(i)+pi;endendw3=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曲線');gridon;plotediton;在MATLABE下運(yùn)行：關(guān)于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)));fori=1:360;ifx3(i)<0x3(i)=x3(i)+pi;endendw3=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,[A3;A4],180*x1./pi,ae/1000);legend('a3','a4','ae');xlabel('x1/(°)