机械原理大作业牛头刨床 matlab编程_百度文库

1、子程序：function F=myfun(x,theta1,l1,l3,l4,h,h1,h2F=l4*cosd(x(2+x(1*cosd(x(3-l1*cosd(theta1-h2;l4*sind(x(2+x(1*sind(x(3-l1*sind(theta1-h1;l4*cosd(x(2+l3*cosd(x(3-x(4;l4*sind(x(2+l3*sind(x(3-h;主程序：theta1=0:360;s3=0:360;theta3=0:360;theta4=0:360;sc=0:360;v3=0:360;omega3=0:360;omega4=0:360;vc=0:360;a3=0:36

2、0;alpha3=0:360;alpha4=0:360;ac=0:360;m=0.5;0;70;0.47;l1=160;l3=960;l4=input('Please input Lde 'h=900;h1=460;h2=120;for i=1:1:361x = fsolve(x myfun(x,theta1(i,l1,l3,l4,h,h1,h2,m(1;m(2;m(3;m(4;m=x;s3(i=x(1;theta4(i=x(2;theta3(i=x(3;sc(i=x(4;endfor i=1:1:361V1=cosd(theta3(i -l4*sind(theta4(i -s

3、3(i*sind(theta3(i 0;sind(theta3(i l4*cosd(theta4(i s3(i*cosd(theta3(i 0;0 -l4*sind(theta4(i -l3*sind(theta3(i -1;0 l4*cosd(theta4(i l3*cosd(theta3(i 0;V2=-l1*sind(theta1(i;l1*cosd(theta1(i;0;0;V3=V1V2;v3(i=V3(1;omega4(i=V3(2;omega3(i=V3(3;vc(i=V3(4;A111= omega3(i*sind(theta3(i l4*omega4(i*cosd(theta

4、4(i v3(i*sind(theta3(i+s3(i*omega3(i*cosd(theta3(i 0;-omega3(i*cosd(theta3(i l4*omega4(i*sind(theta4(i -v3(i*cosd(theta3(i+s3(i*omega3(i*sind(theta3(i 0;0 l4*omega4(i*cosd(theta4(i l3*omega3(i*cosd(theta3(i 0;0 l4*omega4(i*sind(theta4(i l3*omega3(i*sind(theta3(i 0;A112=v3(i;omega4(i;omega3(i;vc(i;A1

5、1=A111*A112;A12=-l1*cosd(theta1(i;-l1*sind(theta1(i;0;0;A1=A11+A12;A2=cosd(theta3(i -l4*sind(theta4(i -s3(i*sind(theta3(i 0;sind(theta3(i l4*cosd(theta4(i s3(i*cosd(theta3(i 0;0 -l4*sind(theta4(i -l3*sind(theta3(i -1;0 l4*cosd(theta4(i l3*cosd(theta3(i 0;A3=A2A1;a3(i=A3(1;alpha4(i=A3(2;alpha3(i=A3(3

6、;ac(i=A3(4;endfigure(1;plot(theta1,theta3,'k',theta1,theta4,'k:'grid on;xlabel(theta_1( o ;ylabel(theta_3 theta_4 ( o ;legend('theta_3','theta_4'figure(2plot(theta1,omega3,'k',theta1,omega4,'k:'grid on;xlabel(theta_1( o ;ylabel('omega_3 omega_4 (r

7、a d / s'legend('omega_3','omega_4' figure(3plot(theta1,alpha3,'k',theta1,alpha4,'k:'grid on;xlabel(theta_1( o ;ylabel('alpha_3 alpha_4 (r a d / s2'legend('alpha_3','alpha_4'figure(4 plot(theta1,sc,'k'grid on;xlabel(theta_1( o ;ylabel('s_c (mm'figure(5 plot(theta1,vc,'k'grid on;xlabel(theta_1( o