Scribd
Upload
18 views5 pages

Assignment 4 Magnification Factor

The document discusses a student's submission of MATLAB code for calculating the magnification factor and phase angle of a vibrating system. The code analyzes how damping ratio and forcing frequency affect maximum vibration amplitude. It also plots phase angle versus frequency ratio for different damping values.

Uploaded by

Muhammad Tayyab
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
Download as pdf or txt
18 views5 pages

Assignment 4 Magnification Factor

The document discusses a student's submission of MATLAB code for calculating the magnification factor and phase angle of a vibrating system. The code analyzes how damping ratio and forcing frequency affect maximum vibration amplitude. It also plots phase angle versus frequency ratio for different damping values.

Uploaded by

Muhammad Tayyab
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
Download as pdf or txt
Download as pdf or txt
You are on page 1/ 5

Mechanical Vibrations

Assignment # 4

Magnification Factor

Submission Date: 20-11-2023

Name: Muhammad Tayyab

Roll No 035

Batch: BSME 2020-24

Submitted To Dr. Zahid Iqbal


MATLAB Code for Magnification Factor
% System parameters
m = 1; % mass
%c = 0.1; % damping coefficient
k = 1; % stiffness
F0 = 1; % amplitude of the force
delta=F0/k;
% Frequency range
omega_range = 0: 0.1: 10; % adjust the range as needed
Cc=2*sqrt(k*m);
% Damping ratios to be considered
zeta_values = [0.1, 0.2, 0.3, 0.4, 0.5,1, 1.2, 1.5];

% Initialize arrays to store results


max_amplitudes = zeros(length(zeta_values), length(omega_range));

% Loop over each damping ratio


for i = 1:length(zeta_values)
zeta = zeta_values(i);
c=zeta*Cc
% Loop over each forcing frequency
for j = 1:length(omega_range)
omega = omega_range(j);

% Solve the differential equation


tspan = [0, 20];
initial_conditions = [0; 0]; % initial position and velocity
[t, x] = ode45(@(t, x) [(x(2)); (F0 * sin(omega * t) - c * x(2) - k *
x(1)) / m - 2 * zeta * sqrt(k/m) * x(2)], tspan, initial_conditions);

% Find maximum amplitude


max_amplitudes(i, j) = max(x(:, 1));
end
end

% Plot results
figure;
plot(omega_range, max_amplitudes/delta);
xlabel('Frequency Ratio (\omega/\omega_n)');
ylabel('Maximum Amplitude');
title('Force Vibration with Viscous Damping');
legend(strcat('\zeta = ', cellstr(num2str(zeta_values'))));
grid on;
Results
MATLAB Code for Phase Angle
% Define the range of r values
r = linspace(0, 5, 100);

% Define different values for zeta


zeta_values = [0.1,0.2,0.3,0.4 0.5, 1, 2];

% Plot phi vs r for each zeta value


figure;
hold on;

for i = 1:length(zeta_values)
zeta = zeta_values(i);
phi_rad =atan((2 * zeta * r) ./ (1 - r.^2));

% Convert phi from radians to degrees


phi_deg = rad2deg(phi_rad);

% Ensure phi is in the range [0, 180]


phi_deg(phi_deg < 0) = phi_deg(phi_deg < 0) + 180;

plot(r, phi_deg, 'DisplayName', ['\zeta = ' num2str(zeta)]);


end

hold off;

% Add labels and legend


xlabel('r');
ylabel('\phi (degrees)');
title('Plot of \phi vs r for different \zeta values');
legend('show');
Results

You might also like