机械原理 牛头刨床附录MATLAB编程部分.docVIP

机械原理课程设计 课题 ：牛头刨床机构设计分析 指导老师 ：王玉丹 附录： MATLAB编程分析 源程序： clear all;clc; %初始条件 theta1=linspace(-20.77,339.23,100);%单位度 theta1=theta1*pi/180;%转换为弧度制 W1=80*pi/30;%角速度 单位rad/s H=0.5;%行程 单位m L1=0.1605;%O2A的长度 单位m L3=0.6914;%O3B的长度 单位m L4=0.2074;%BF的长度 单位m L6=0.430;%O2O3的长度 单位m L6u=0.6572;%O3D的长度 单位m Z=pi/180;%角度与弧度之间的转换 dT=(theta1(3)-theta1(2))/W1;%时间间隔 for j=1:100 t(j)=dT*(j-1);%时间因素 end %求解S3、Theta3、Theta4和SE四个变量 S3=((L6)^2+(L1)^2-2*L6*L1*cos(theta1+pi/2)).^0.5;%求出O3A的值 for i=1:100%求解角度theta3、Theta4和SE的长度 theta3(i)=acos(L1*cos(theta1(i))/S3(i)); theta4(i)=asin((L6u-L3*sin(theta3(i)))/L4); SE(i)=L3*cos(theta3(i))+L4*cos(theta4(i)); end%求解完成 %求解完成 %求解VS3、W3、W4和VE四个变量 for i=1:100 J= inv([cos(theta3(i)),-S3(i)*sin(theta3(i)),0,0; sin(theta3(i)),S3(i)*cos(theta3(i)),0,0; 0,-L3*sin(theta3(i)),-L4*sin(theta4(i)),-1; 0,L3*cos(theta3(i)),L4*cos(theta4(i)),0]); K=J*W1*[-L1*sin(theta1(i));L1*cos(theta1(i));0;0]; VS3(i)=K(1); W3(i)=K(2); W4(i)=K(3); VE(i)=K(4); end%求解完成 %求解aS3、a3、a4、aE四个变量 for i=1:100 J= inv([cos(theta3(i)),-S3(i)*sin(theta3(i)),0,0; sin(theta3(i)),S3(i)*cos(theta3(i)),0,0; 0,-L3*sin(theta3(i)),-L4*sin(theta4(i)),-1; 0,L3*cos(theta3(i)),L4*cos(theta4(i)),0]); P=W1*W1*[-L1*cos(theta1(i));-L1*sin(theta1(i));0;0]; M=[-W3(i)*sin(theta3(i)),-VS3(i)*sin(theta3(i))-S3(i)*W3(i)*cos(theta3(i)),0,0; W3(i)*cos(theta3(i)),VS3(i)*cos(theta3(i))-S3(i)*W3(i)*sin(theta3(i)),0,0; 0,-L3*W3(i)*cos(theta3(i)),-L4*W4(i)*cos(theta4(i)),0; 0,-L3*W3(i)*sin(theta3(i)),-L4*W4(i)*sin(theta4(i)),0]; N=[VS3(i);W3(i);W4(i);VE(i)]; K=J*(-M*N+P); aS3(i)=K(1); a3(i)=K(2); a4(i)=K(3); aE(i)=K(4); end%求解完成 %动态静力分析 %初始条件 M4=20; M5=3; M6=62; Js4=1.2; Js5=0.0255; Fc=1500; Ls4=0.5*L3; Ls5=0.5*L4; %给切削阻力赋值 for i=1:100 if((abs(SE(1)-SE(i))0.05*Habs(SE(1)-SE(i))0.