Scribd
Upload
25 views

%% Problem 1 - Single Rotor Blade: 'Mode Shape 1'

The document analyzes the natural frequencies and mode shapes of a single rotor blade system and a 4 rotor blade system. For the single rotor blade, it calculates the mass and stiffness matrices, finds the natural frequencies and mode shapes by solving the eigenvalue problem. The first natural frequency is 786.16 Hz and the second is 2794.9 Hz. For the 4 rotor blade system, it similarly calculates the mass and stiffness matrices for the now larger system, finds the 8 natural frequencies and corresponding mode shapes by solving the eigenvalue problem. The first natural frequency is 793.09 Hz and the frequencies range up to 2478 Hz.

Uploaded by

Andrew
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
25 views

%% Problem 1 - Single Rotor Blade: 'Mode Shape 1'

The document analyzes the natural frequencies and mode shapes of a single rotor blade system and a 4 rotor blade system. For the single rotor blade, it calculates the mass and stiffness matrices, finds the natural frequencies and mode shapes by solving the eigenvalue problem. The first natural frequency is 786.16 Hz and the second is 2794.9 Hz. For the 4 rotor blade system, it similarly calculates the mass and stiffness matrices for the now larger system, finds the 8 natural frequencies and corresponding mode shapes by solving the eigenvalue problem. The first natural frequency is 793.09 Hz and the frequencies range up to 2478 Hz.

Uploaded by

Andrew
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 7

%% Problem 1 - Single Rotor Blade

clc
m1=1.81;
m2=4.54;
k1=1.16*10^6;
k2=3.42*10^7;
kc=3*10^7;
M=[m1 0;0 m2];
K=[k1 -k1;-k1 k1+k2];
A=inv(M)*K
[V,D]=eig(A)
w1=sqrt(D(1,1))
w2=sqrt(D(2,2))
modeshape1=[V(1,1);V(2,1)]
modeshape2=[V(1,2);V(2,2)]
subplot(2,1,1)
plotv(modeshape1)
title('Mode Shape 1')
subplot(2,1,2)
plotv(modeshape2)
title('Mode Shape 2')
Results:
A =
1.0e+06 *
0.6409
-0.6409
-0.2555
7.7885
V =
-0.9994
0.0890
-0.0356
-0.9960
D =
1.0e+06 *
0.6180
0
0
7.811
w1 = 786.1599
w2 = 2.7949e+03
modeshape1 =
-0.9994
-0.0356
modeshape2 =
0.0890
-0.9960

%% Problem 2- Four Rotor Blades


clc
NewM=[m1 0 0 0 0 0 0 0; 0 m2 0 0 0 0 0 0; 0 0 m1 0 0 0 0 0; 0 0 0 m2 0 0 0 0;
0 0 0 0 m1 0 0 0; 0 0 0 0 0 m2 0 0; 0 0 0 0 0 0 m1 0; 0 0 0 0 0 0 0 m2];
NewK=[k1 -k1 0 0 0 0 0 0; -k1 k1+k2+kc 0 0 0 0 0 0; 0 0 k1 -k1 0 0 0 0; 0 0
-kc-k1 kc+k1+k2 0 0 0 0; 0 0 0 0 k1 -k1 0 0; 0 0 0 0 -kc-k1 kc+k2+k1 0 0; 0 0
0 0 0 0 k1 -k1; 0 0 0 0 0 0 k2-k1 kc+k1];
B=inv(NewM)*NewK
[V,D]=eig(B)
w1=sqrt(D(1,1))
w2=sqrt(D(2,2))
w3=sqrt(D(3,3))
w4=sqrt(D(4,4))
w5=sqrt(D(5,5))
w6=sqrt(D(6,6))
w7=sqrt(D(7,7))
w8=sqrt(D(8,8))
modeshape1=[V(1,1);V(2,1)]
modeshape2=[V(1,2);V(2,2)]
modeshape3=[V(3,3);V(4,3)]
modeshape4=[V(3,4);V(4,4)]
modeshape5=[V(5,5);V(6,5)]
modeshape6=[V(5,6);V(6,6)]
modeshape7=[V(7,7);V(8,7)]
modeshape8=[V(7,8);V(8,8)]
subplot(4,2,1)
plotv(modeshape1)
title('Mode Shape 1')
subplot(4,2,2)
plotv(modeshape2)
title('Mode Shape 2')
subplot(4,2,3)
plotv(modeshape3)
title('Mode Shape 3')
subplot(4,2,4)
plotv(modeshape4)
title('Mode Shape 4')
subplot(4,2,5)
plotv(modeshape5)
title('Mode Shape 5')
subplot(4,2,6)
plotv(modeshape6)
title('Mode Shape 6')
subplot(4,2,7)
plotv(modeshape7)
title('Mode Shape 7')
subplot(4,2,8)
plotv(modeshape8)
title('Mode Shape 8')

Results:

B=
1.0e+07 *
Columns 1 through 6
0.0641 -0.0641
-0.0256

0.0641 -0.0641

1.4396

0 -0.6863

1.4396

0 -0.6863

1.4396

0.0641 -0.0641

Columns 7 through 8

0.0641 -0.0641
0.7278

0.6863

V=
Columns 1 through 6
-0.9998

0.0465

-0.0186 -0.9989

0 -0.8987

0.0455

0 -0.4385 -0.9990

0 -0.8987

0 -0.4385 -0.9990

Columns 7 through 8
0

0.0455

-0.5924

0.1189

0.8057 -0.9929
D=
1.0e+07 *
Columns 1 through 6
0.0629

1.4408

0.0328

1.4709

0.0328

1.4709

Columns 7 through 8
0

0.1513
0

0.5992

w1 =793.0889
w2 = 3.7958e+03
w3 = 572.9031
w4 = 3.8352e+03
w5 = 572.9031
w6 = 3.8352e+03
w7 = 1.2298e+03
w8 = 2.4478e+03

modeshape1 =
-0.9998
-0.0186
modeshape2 =
0.0465
-0.9989
modeshape3 =
-0.8987
-0.4385
modeshape4 =
0.0455
-0.9990
modeshape5 =
-0.8987
-0.4385
modeshape6 =
0.0455
-0.9990
modeshape7 =
-0.5924
0.8057
modeshape8 =
0.1189
-0.9929

You might also like