资源描述:
《 转子动力学作业(传递矩阵法).doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、%求解转子系统前三个临界转速和主振型的传递矩阵法clcclear%等截面轴参数l1=0.12;d=0.04;A=pi*d*d/4;%轮盘参数D=0.5;h=0.025;%盘轴材料参数(忽略轴的质量)a=1;u=0.3;rou=7800;E=2.0e11;G=E/(2*(1+u));I=pi*(d^4)/64;K1=2.0e7;v1=6*E*I/(a*G*A*l1*l1);mi=rou*pi*D^2/4*h;%轮盘的集质量Jp=mi*D^2/8;Jd=Jp/2;Ji=Jp-Jd;%参数的数组形式L=[l1l1l1l1l1l1l1l1l1l1l
2、1l1l100];M=[0mimimimimimi00000mimi0];K=[K1000000K1000K1000];v=[v1v1v1v1v1v1v1v1v1v1v1v1v100];J=[0JiJiJiJiJiJi00000JiJi0];k=0;Tit=['第一阶频率的振型和弯矩图';'第二阶频率的振型和弯矩图';'第三阶频率的振型和弯矩图'];forw=0:0.01:4000;fori=1:15;T(:,:,i)=[1+(L(i)^3)*(1-v(i))*(M(i)*w^2-K(i))/(6*E*I)L(i)+L(i)^2*J(i)*
3、w^2/(2*E*I)L(i)^2/(2*E*I)L(i)^3*(1-v(i))/(6*E*I);(L(i)^2)*(M(i)*w^2-K(i))/(2*E*I)1+L(i)*J(i)*w^2/(E*I)L(i)/(E*I)L(i)^2/(2*E*I);L(i)*(M(i)*w^2-K(i))J(i)*w^21L(i);M(i)*w^2-K(i)001];endH=T(:,:,1);fori2=2:15;H=T(:,:,i2)*H;endF=H(3,1)*H(4,2)-H(3,2)*H(4,1);ifF*(-1)^k<0%求解临界转速k=k
4、+1;wi(k)=w;w=wi(k)ni(k)=wi(k)*30/pi;endendfori1=1:3;w=wi(i1);forj=1:14;T(:,:,j)=[1+(L(j)^3)*(1-v(j))*(M(j)*w^2-K(j))/(6*E*I)L(j)+L(j)^2*J(j)*w^2/(2*E*I)L(j)^2/(2*E*I)L(j)^3*(1-v(j))/(6*E*I);(L(j)^2)*(M(j)*w^2-K(j))/(2*E*I)1+L(j)*J(j)*w^2/(E*I)L(j)/(E*I)L(j)^2/(2*E*I);L(j)*
5、(M(j)*w^2-K(j))J(j)*w^21L(j);M(j)*w^2-K(j)001];endH=T(:,:,1);forj=2:15;H=T(:,:,j)*H;endb=-H(4,1)/H(4,2);X(:,1)=([1b00]');forn=2:16;X(:,n)=T(:,:,n-1)*X(:,n-1);%相邻两质点右边的传递关系endforj1=1:15;y(j1)=X(1,j1);z(j1)=X(3,j1);x(j1)=(j1-1)*l1;endy(16)=X(1,16);x(16)=1.56;z(16)=X(3,16);y=
6、y/max(abs(y));%归一化z=z/max(abs(z));subplot(3,1,i1)plot(x,y,'b-',x,z,'r:')title(Tit(i1,:))xlabel('轴长'),ylabel('不平衡值')axis([0,1.56,-1.2,1.2])gridonz;endlegend('振型','弯矩')niwini=1.0e+004*0.14680.20650.52541.38372.37592.38323.10363.5473wi=1.0e+003*0.15370.21620.55021.44902.48812
7、.49563.25013.7147%转子系统的不平衡响应clcclearww=[153.68216.23550.2214492488.12495.63250.13714.7]%前8阶固有频率n=ww*30/pi%前8阶转频wi=[0.9*ww(1)(ww(1)+ww(2))/2]%0.9w(1)和(w(1)+w(2))/2)%等截面轴参数l1=0.12;d=0.04;A=pi*d*d/4;%轮盘参数D=0.5;h=0.025;%盘轴材料参数(忽略轴的质量)rou=7800;E=2.0e11;I=pi*(d^4)/64;K1=2.0e7;m=
8、rou*pi*D^2/4*h;%轮盘的集质量Jp=m*D^2/8;Jd=Jp/2;J1=Jp-Jd;u1=0.8e-4;%参数的数组形式L=[l1l1l1l1l1l1l1l1l1