Dynamics of structures is a special branch of structural analysis, which deals with the behaviour of structures subjected to dynamic loads (loads that vary with time). Such loads develop dynamic reactions, internal forces, and displacements in the structure. These values all change with time, and maximum values often exceed static load reaction values. Dynamic analysis of any structure often begins with the free vibration analysis.

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 truss is loaded as shown below with a lumped mass of 5000kg. We are to obtain the

*eigenfrequencies*and mode shapes by assuming linear behaviour, neglecting damping, and assuming that the stiffness and inertial effects are time independent.

E = 205 KN/mm

^{2}

A = 0.00548 m

^{2}(UB 254 x 146 x 43)

A little consideration will show that the structure has two degrees of freedom which are given by the vertical and horizontal translation at node C. The lumped mass at the node will definitely participate in those movements. The displacements are given by Y

_{1}and Y

_{2}below.

**Geometrical Properties**

Î¸ = tan

^{-1}(3/4) = 36.869°

cos Î¸ = 0.8

sin Î¸ = 0.6

L

_{AC}= √(3

^{2}+ 4

^{2}) = 5m

Î³ = tan

^{-1}(4/4) = 45°

cos Î³ = 0.7071

sin Î³ = 0.7071

L

_{BC}= √(3

^{2}+ 3

^{2}) = 4.243m

**ANALYSIS OF STATE 1**; Y

_{1}= 1.0

**∑F**

_{y}= 0-0.6F

_{AC}- 0.7071F

_{CB}= 1.0 ----------- (1)

**∑F**

_{x}= 0-0.8F

_{AC}+ 0.7071F

_{CB}= 0 ----------- (2)

Solving (1) and (2) simultaneously

F

_{AC}= -0.7143

F

_{CB}= -0.8081

**ANALYSIS OF STATE 2**; Y

_{2}= 1.0

**∑F**

_{y}= 0-0.6F

_{AC}- 0.7071F

_{CB}= 0 ----------- (3)

**∑F**

_{x}= 0-0.8F

_{AC}+ 0.7071F

_{CB}= -1.0 ----------- (4)

Solving (3) and (4) simultaneously

F

_{AC}= 0.7143

F

_{CB}= -0.6061

The frequency equation of a 2 degree of freedom system in terms of displacement is given by;

Where;

M is the lumped mass

A is the amplitudes (Note that the system does not allow us to find the amplitudes but we can find the ratio of the amplitudes)

Ï‰ is the angular velocity

And of course Î´

_{ij }is the axial displacement given by;For the system above, we can easily compute this;

Î´

_{11}= 1/EA[(-0.7143 × -0.7143 × 5) + (-0.8081 × -0.8081 × 4.243)] = 5.3219/EA

Î´

_{22}= 1/EA[(0.7143 × 0.7143 × 5) + (-0.6061 × -0.6061 × 4.243)] = 4.1098/EA

Î´

_{21}= Î´

_{12}= 1/EA[(-0.7143 × 0.7143 × 5) + (-0.6011 × -0.8081 × 4.243)] = -0.4729/EA

To simplify our calculations, let Î´

_{0}= 1/EA

Also let Î» = 1/MÎ´

_{0}Ï‰

^{2}

^{}Therefore;

(5.3219 - Î»)A

_{1}- 0.4729A

_{2}= 0

0.4729A

_{1}+ (4.1098 - Î»)A

_{2}= 0

On expanding;

Î»

^{2}- 9.4317Î» + 21.6484 = 0

On solving the quadratic equation;

Î»

_{1}= 5.4845

Î»

_{2}= 3.9472

You can verify that Ï‰

^{2}= EA/Î»M

Therefore;

Ï‰

_{1}

^{2}= (0.1823 × 205 × 10

^{6}× 0.00548)/5000 = 40.959

Ï‰

_{2}

^{2}= (0.2533 × 205 × 10

^{6}× 0.00548)/5000 = 56.911

Hence;

Ï‰

_{1 }= 6.399 rad/s

Ï‰

_{2}

^{}

_{ }= 7.544 rad/s

To obtain the mode shapes;

For the first mode of vibration; Let A

_{1}= 1.0

(5.3219 - Î»)A

_{1}- 0.4729A

_{2}= 0

(5.3219 - 5.4845) - 0.4729 A

_{2}= 0

Therefore, A

_{2}= -0.1626/0.4729 = -0.344

Î¦

_{1}= [1.000 -0.344]

^{T}

For the second mode of vibration; Let A

_{1}= 1.0

(5.3219 - Î»)A

_{1}- 0.4729A

_{2}= 0

(5.3219 - 3.9472) - 0.4729 A

_{2}= 0

Therefore, A

_{2}= 1.3737/0.4729 = 2.9069

Î¦

_{2}= [1.000 2.907]

^{T}

The modal matrix is then defined as;

Remember Structville....

We have an interesting Facebook interactive page;

Visit us at www.facebook.com/structville

## No comments:

## Post a Comment