Scribd
Upload
25 views

Matlab Code

Uploaded by

ash17nazurally
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
25 views

Matlab Code

Uploaded by

ash17nazurally
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/ 2

11/1/24 2:54 PM MATLAB Command Window 1 of 2

% Given values
k = 1000; % Spring constant in N/m
m = 10; % Mass in kg
delta = 1e-3; % Initial displacement for m/2 mass in meters

% Define the Mass Matrix (M)


M = [m 0 0;
0 m 0;
0 0 m/2];

% Define the Stiffness Matrix (K)


K = [3*k -2*k 0;
-2*k 3*k -k;
0 -k 3*k];

% (c) Calculate Natural Frequencies and Mode Shapes


% Solve the generalized eigenvalue problem K*Phi = lambda*M*Phi
[Phi, Lambda] = eig(K, M);
omega = sqrt(diag(Lambda)); % Natural frequencies

% Display Natural Frequencies


disp('Natural Frequencies (rad/s):');
disp(omega);

% (d) Plot mode shapes


figure;
plot([1 2 3], Phi(:,1), '-o', 'DisplayName', 'Mode 1'); hold on;
plot([1 2 3], Phi(:,2), '-o', 'DisplayName', 'Mode 2');
plot([1 2 3], Phi(:,3), '-o', 'DisplayName', 'Mode 3');
legend;
xlabel('Mass Index');
ylabel('Amplitude');
title('Mode Shapes');

% (e) Define initial conditions and solve for the time response
x0 = [0; 0; delta]; % Initial displacement (only third mass has initial
displacement)
v0 = [0; 0; 0]; % Initial velocity (all masses start at rest)
initial_conditions = [x0; v0]; % Combine initial displacements and velocities

% Define time span


tspan = [0 1]; % Time from 0 to 1 second

% Solve the system using ode45


[t, X] = ode45(@(t, X) system_dynamics(t, X, M, K), tspan, initial_conditions);

% (f) Plot the response of each mass over time


figure;
plot(t, X(:,1), 'r', 'DisplayName', 'Mass 1'); hold on;
plot(t, X(:,2), 'g', 'DisplayName', 'Mass 2');
plot(t, X(:,3), 'b', 'DisplayName', 'Mass 3');
legend;
xlabel('Time (s)');
11/1/24 2:54 PM MATLAB Command Window 2 of 2

ylabel('Displacement (m)');
title('Response of Each Mass Over Time');

% Function to define the system dynamics for ode45


function dXdt = system_dynamics(~, X, M, K)
n = length(X)/2; % Number of degrees of freedom
dXdt = zeros(2*n, 1); % Initialize derivative vector

% Split the state vector into displacements and velocities


displacements = X(1:n);
velocities = X(n+1:end);

% First half of dXdt is velocities


dXdt(1:n) = velocities;

% Second half of dXdt is accelerations


dXdt(n+1:end) = -M \ (K * displacements);
end

You might also like