A hub for civil engineering related designs, analyses, discussion, information, and knowledge...

## Sunday, April 24, 2016

1.1 Introduction
A structure undergoes free vibration when it is brought out of static equilibrium and can then oscillate without any external dynamic excitation. Free vibration of structures occurs with some frequencies which depend only on the parameters of the structures such as boundary conditions, distribution of masses, stiffnesses within the members etc, and not on the reason for the vibration. At each natural frequency of free vibration, the structure vibrates in simple harmonic motion where the displaced shape (mode shape) of the structure is constant but the amplitude of the displacement is varying in a sinusoidal manner with time. The number of natural frequencies in a structure coincides with the number of degrees of freedom in the structure. These frequencies are inherent to the given structure and are often referred to as eigenfrequencies. Each mode shape of vibration shows the form of an elastic curve which corresponds to a specific frequency.

A method of obtaining the natural frequencies and mode shapes of vibration is modal analysis. This is a technique by which the equations of motion, which are originally expressed in physical coordinates, are transformed to modal coordinates using the eigenvalues and eigenvectors gotten by solving the undamped frequency eigenproblem. The transformed equations are called modal equations.
For an undamped free vibration, the equation of motion is;

MÃ¼ + ku = 0 ----------(1)

In the paper downloadable in this post, modal dynamic analysis was carried out on the three storey frame shown above. The results obtained regarding the three assumed degrees of freedom are shown in the pictures below;

(a) Mode 1 vibration parameters and displaced shape

(b) Mode 2 vibration parameters and displaced shape

(c) Mode 3 vibration parameters and displaced shape

The displacement time history response is therefore;

u1(t) = 0.00128448 sin(5.0261t) + 0.00044346 sin(13.7654t) + 0.0002408175 sin(19.0459t)
u2(t) = 0.00107836 sin(5.0261t) – 0.00009 sin(13.7654t) - 0.000313845 sin(19.0459t)
u3(t) = 0.00061299 sin(5.0261t) – 0.00046104 sin(13.7654t) + 0.0002157225 sin(19.0459t)

If velocity (v) = du/dt; The velocity time history response is therefore;
The velocity time history response is obtained by direct differentiation of the displacement time history equations;

v1(t)= 0.0064559 cos(5.0261t) + 0.0061044 cos(13.7654t) + 0.0045865cos(19.0459t)
v2(t)= 0.0054199 cos(5.0261t) – 0.0012388 cos(13.7654t) - 0.00597746 cos(19.0459t)
v3(t)= 0.0030809 cos(5.0261t) – 0.0063464 cos(13.7654t) + 0.0041086 cos(19.0459t)

The acceleration (a) time history response is obtained by direct differentiation of the velocity time history equations;
a1(t) = -0.032448 sin(5.0261t) - 0.0840295 sin(13.7654t) – 0.087354 sin(19.0459t)
a2(t) = -0.02724 sin(5.0261t) + 0.0170525 sin(13.7654t) + 0.113846 sin(19.0459t)
a3(t) = -0.0154849 sin(5.0261t) + 0.0873607 sin(13.7654t) – 0.0782519 sin(19.0459t)