搜索
%求解转子系统前三个临界转速和主振型的传递矩阵法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=[l1l1l1l1l1l1l1l1l1l1l1l1l100];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)*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)^k0%求解临界转速k=k+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)*(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=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('振型','弯矩')niwi00.511.5-101第一阶频率的振型和弯矩图轴长不平衡值00.511.5-101第二阶频率的振型和弯矩图轴长不平衡值00.511.5-101第三阶频率的振型和弯矩图轴长不平衡值振型弯矩ni=1.0e+004*0.14680.20650.52541.38372.37592.38323.10363.5473wi=1.0e+003*0.15370.21620.55021.44902.48812.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=rou*pi*D^2/4*h;%轮盘的集质量Jp=m*D^2/8;Jd=Jp/2;J1=Jp-Jd;u1=0.8e-4;%参数的数组形式L=[l1l1l1l1l1l1l1l1l1l1l1l1l1l1];M=[0mmmmmm00000mm];K=[K1000000K1000K100];J=[0J1J1J1J1J1J100000J1J1];Tit=['wi(1)时的振动响应图';'wi(2)时的振动响应图'];U=[0000000000000u1];fori=1:2w=wi(i);n(i)=w*30/piforj=1:14;T(:,:,j)=[1+(L(j)^3)*(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/(6*E*I)L(j)^3/(6*E*I)*U(j)*w^2;(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)^2/(2*E*I)*U(j)*w^2;L(j)*(M(j)*w^2-K(j))J(j)*w^21L(j)L(j)*U(j)*w^2;M(j)*w^2-K(j)001U(j)*w^2;00001];endG=T(:,:,1);forj1=2:14;H=T(:,:,j1)*G;G=H;endD1=H([34],[12]);B=H([34],[52]);B(:,1)=-B(:,1);C=H([34],[15]);C(:,2)=-C(:,2);b=det(B)/det(D1);c=det(C)/det(D1);X(:,1)=([bc001]');forn=2:14;X(:,n)=T(:,:,n-1)*X(:,n-1);%相邻两质点右边的传递关系endy(1)=X(1,1);x(1)=0;forj2=2:13;y(j2)=X(1,j2);x(j2)=x(j2-1)+L(j2-1);endy(14)=X(1,14);x(14)=1.56;xi=0:0.05:1.56;yi=interp1(x,y,xi,'spline');subplot(2,1,i)plot(xi,yi,'b-o','LineWidth',1.5)title(Tit(i,:))xlabel('轴长'),ylabel('不平衡值')gridonend00.20.40.60.811.21.41.6-1012x10-6wi(1)时的振动响应图轴长不平衡值00.20.40.60.811.21.41.6-20246x10-6wi(2)时的振动响应图轴长不平衡值ww=1.0e+003*0.15370.21620.55021.44902.48812.49563.25013.7147n=1.0e+004*0.14680.20650.52541.38372.37602.38313.10363.5473wi=138.3120184.9550n=1.0e+003*0.01501.766200.20.40.60.811.21.41.6-101第一阶频率的振型和弯矩图轴长不平衡值00.20.40.60.811.21.41.6-101第二阶频率的振型和弯矩图轴长不平衡值00.20.40.60.811.21.41.6-101第三阶频率的振型和弯矩图轴长不平衡值振型弯矩00.511.5-101第一阶频率的振型和弯矩图轴长不平衡值00.511.5-101第二阶频率的振型和弯矩图轴长不平衡值00.511.5-101第三阶频率的振型和弯矩图轴长不平衡值振型弯矩00.511.5-101第一阶频率的振型和弯矩图轴长不平衡值00.511.5-101第二阶频率的振型和弯矩图轴长不平衡值00.511.5-101第三阶频率的振型和弯矩图轴长不平衡值振型弯矩00.20.40.60.811.21.41.6-10123x10-6wi(1)时的振动响应图轴长不平衡值00.20.40.60.811.21.41.6-50510x10-6wi(2)时的振动响应图轴长不平衡值00.20.40.60.811.21.41.6-1012x10-6wi(1)时的振动响应图轴长不平衡值00.20.40.60.811.21.41.6-20246x10-6wi(2)时的振动响应图轴长不平衡值