ziyouzhendong.asv

内燃机转子仿真

clear all;
clc;
close all;
n=10;
jj=zeros(n,n);
kk=zeros(n,n);

E=[9.8790 2.6700 8.1140 5.9840 8.1140 8.1140 5.9840 8.1140 64.6100 470.8800];        %inertia(kg*m**2)
F=[0.679*1e6 0.2096*1e9 0.3413*1e8 0.3413*1e8 0.3413*1e8 0.3413*1e8 0.3413*1e8 0.66*1e8 0.981*1e8 0]; %stiffness(n*m/rad)

for i=1:n
    jj(i,i)=E(i);
end
%*********************刚度系数矩阵
for i=2:n
    kk(i,i)=F(i-1)+F(i);                  
end
for i=1:n-1
    kk(i,i+1)=-F(i);
    kk(i+1,i)=-F(i);
end
kk(1,1)=F(1);
[fx,fn]=eig(inv(jj)*kk);

for i=1:10

fx(:,i)=fx(:,i)/max(abs(fx(:,i)));

end
fx
nn=(1:10);

figure(1)
subplot(5,2,10)
plot(nn,fx(:,1),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,9)
plot(nn,fx(:,2),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,8)
plot(nn,fx(:,3),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,7)
plot(nn,fx(:,4),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,6)
plot(nn,fx(:,5),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,5)
plot(nn,fx(:,6),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,4)
plot(nn,fx(:,7),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,3)
plot(nn,fx(:,8),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,1)
plot(nn,fx(:,9),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;

subplot(5,2,2)
plot(nn,fx(:,10),'k')
xlim([1,10]);
xlabel('惯量号')
ylabel('归一化振幅');
grid on;