Scribd
Upload
35 views3 pages

Question 7

This document analyzes the natural frequencies and damping ratios of a multi-degree of freedom system. It calculates the system matrices, finds the eigenvalues and eigenvectors, and uses the eigenvalues to determine the natural frequencies and damping ratios through solving equations. It provides the results of the analysis.

Uploaded by

omar mohsen
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
35 views3 pages

Question 7

This document analyzes the natural frequencies and damping ratios of a multi-degree of freedom system. It calculates the system matrices, finds the eigenvalues and eigenvectors, and uses the eigenvalues to determine the natural frequencies and damping ratios through solving equations. It provides the results of the analysis.

Uploaded by

omar mohsen
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Question 7 M File

% Mohamed Yasser
% Problem Set 1 Question (7)
close all
clear
clc

global s

M = [10 , 0 , 0 ; 0 , 20 , 0 ; 0 , 0 , 30];
C = [15 , -5 , 0 ; -5 , 20 , -15 ; 0 , -15 , 15];
K = [140 , -60 , 0 ; -60 , 100 , -40 ; 0 , -40 , 40];
Q = zeros(length(M),length(M));

A = [M , Q ; C , M];
B = [Q , -M ; K , Q];
% A&B = System Augmented Matrices

[U,L] = eig(-A\B);
% U = eigenvectors/modes of vibrations
% L = eigenvalues

w_positive = zeros(1,length(L));
z_positive = zeros(1,length(L));

for i = 1:length(L)
s = L(i,i);
fun = @solve_natfreq_dampratio_positive;
X0 = [5 , 1/sqrt(2)];
X = fsolve(fun,X0);
w_positive(i) = X(1);
z_positive(i) = X(2);
end

w_negative = zeros(1,length(L));
z_negative = zeros(1,length(L));

for i = 1:length(L)
s = L(i,i);
fun = @solve_natfreq_dampratio_negative;
X0 = [5 , 1/sqrt(2)];
X = fsolve(fun,X0);
w_negative(i) = X(1);
z_negative(i) = X(2);
end

for i = 1:length(L)
w_angle = angle(w_positive(i));
z_angle = angle(z_positive(i));

if w_angle ~= 0
w_positive(i) = 0;
end

if z_angle ~= 0
z_positive(i) = 0;
end
end

for i = 1:length(L)
w_angle = angle(w_negative(i));
z_angle = angle(z_negative(i));

if w_angle ~= 0
w_negative(i) = 0;
end

if z_angle ~= 0
z_negative(i) = 0;
end
end

w = w_positive + w_negative;
% w = undamped natural frequency [rad/s]
z = z_positive + z_negative;
% z = damping ratio

w_d = w.*sqrt(1-z.^2);
% w_d = damped natural frequency [rad/s]

Question 7 User Defined Function 1


% Mohamed Yasser
% This function calculates the undamped natural frequency and
damping ratio from the eigenvalues
% x(1) = undamped natural frequency
% x(2) = damping ratio
% s = -x(2)*x(1) + i*x(1)*sqrt(1-(x(2))^2)
function F = solve_natfreq_dampratio_positive(x)
global s
s_r = real(s);
s_i = imag(s);
F(1) = s_r + x(2)*x(1);
F(2) = s_i - x(1)*sqrt(1-(x(2))^2);
end
Question 7 User Defined Function 2
% Mohamed Yasser
% This function calculates the undamped natural frequency and
damping ratio from the eigenvalues
% x(1) = undamped natural frequency
% x(2) = damping ratio
% s = -x(2)*x(1) - i*x(1)*sqrt(1-(x(2))^2)
function F = solve_natfreq_dampratio_negative(x)
global s
s_r = real(s);
s_i = imag(s);
F(1) = s_r + x(2)*x(1);
F(2) = s_i + x(1)*sqrt(1-(x(2))^2);
end

Answers

A = B =
10 0 0 0 0 0 0 0 0 -10 0 0
0 20 0 0 0 0 0 0 0 0 -20 0
0 0 30 0 0 0 0 0 0 0 0 -30
15 -5 0 10 0 0 140 -60 0 0 0 0
-5 20 -15 0 20 0 -60 100 -40 0 0 0
0 -15 15 0 0 30 0 -40 40 0 0 0

w_d =
3.8653 3.8653 1.9329 1.9329 0.7039 0.7039

U =
-0.0503 -0.0503 -0.0396 -0.0396 0.1898 0.1898
-0.2302i +0.2302i -0.2023i +0.2023i +0.0120i -0.0120i
0.0062 0.0062 -0.1048 -0.1048 0.4266 0.4266
+0.0661i -0.0661i -0.3313i +0.3313i +0.0334i -0.0334i
0.0065 0.0065 0.0163 0.0163 0.6697 0.6697
-0.0077i +0.0077i +0.1794i -0.1794i +0.0000i +0.0000i
0.9324 0.9324 0.4153 0.4153 -0.0169 -0.0169
+0.0000i +0.0000i +0.0471i -0.0471i +0.1330i -0.1330i
-0.2607 -0.2607 0.7044 0.7044 -0.0424 -0.0424
-0.0318i +0.0318i +0.0000i +0.0000i +0.2988i -0.2988i
0.0242 0.0242 -0.3568 -0.3568 -0.0296 -0.0296
+0.0317i -0.0317i -0.0782i +0.0782i +0.4714i -0.4714i

You might also like