Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
FEM_Lecture_Notes.pdf
Скачиваний:
171
Добавлен:
16.05.2015
Размер:
3 Мб
Скачать
☆

Lecture Notes: Introduction to Finite Element Method

Chapter 7. Structural Vibration and Dynamics

II. Free Vibration

Study of the dynamic characteristics of a structure:

•natural frequencies

•normal modes (shapes)

Let f(t) = 0 and C = 0 (ignore damping) in the dynamic equation (8) and obtain

Mu&&+Ku = 0

(12)

Assume that displacements vary harmonically with time, that is,

u (t ) = u sin( ωt ), u& (t ) = ω u cos( ωt ),

u&&(t ) = −ω 2 u sin( ωt ),

where u is the vector of nodal displacement amplitudes. Eq. (12) yields,

[K − ω 2 M ]u = 0

(13)

This is a generalized eigenvalue problem (EVP).

Solutions?

© 1997-2003 Yijun Liu, University of Cincinnati

163

Lecture Notes: Introduction to Finite Element Method

Chapter 7. Structural Vibration and Dynamics

Trivial solution: u = 0 for any values of ω (not interesting).

Nontrivial solutions: u ≠ 0 only if

K −ω2M

 

= 0

(14)

 

This is an n-th order polynomial of ω2, from which we can find n solutions (roots) or eigenvalues ωi.

ωi (i = 1, 2, …, n) are the natural frequencies (or characteristic frequencies) of the structure.

ω1 (the smallest one) is called the fundamental frequency.

For each ωi , Eq. (13) gives one solution (or eigen) vector

[K − ω i 2 M ]u i = 0 .

u i (i=1,2,…,n) are the normal modes (or natural modes, mode shapes, etc.).

Properties of Normal Modes

u iT K u j

= 0 ,

 

u iT M u j = 0 ,

for i ≠ j,

(15)

if ωi ≠ωj . That is, modes are orthogonal (or independent) to each other with respect to K and M matrices.

© 1997-2003 Yijun Liu, University of Cincinnati

164

Lecture Notes: Introduction to Finite Element Method

Chapter 7. Structural Vibration and Dynamics

Normalize the modes:

u iT M u i = 1,

 

u iT K u i = ω i2 .

(16)

Note:

•Magnitudes of displacements (modes) or stresses in normal mode analysis have no physical meaning.

•For normal mode analysis, no support of the structure is necessary.

ωi = 0 there are rigid body motions of the whole or a part of the structure.

apply this to check the FEA model (check for mechanism or free elements in the models).

•Lower modes are more accurate than higher modes in the FE calculations (less spatial variations in the lower modesfewer elements/wave length are needed).

© 1997-2003 Yijun Liu, University of Cincinnati

165

Lecture Notes: Introduction to Finite Element Method

Chapter 7. Structural Vibration and Dynamics

Example:

y

v2

 

 

 

ρ, A, EI

 

θ2

1

2

x

L

 

 

 

 

 

v

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[K −ω2M]

2

 

=

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K =

EI

12

 

−6L

 

 

M =

 

ρAL

156

−22L

 

 

 

 

 

 

 

 

 

2

,

 

 

 

 

 

 

 

 

 

 

 

 

2

.

 

3

 

 

 

 

 

 

 

 

 

 

420

 

 

 

 

 

L

−6L

 

 

4L

 

 

 

 

 

 

 

−22L

4L

 

 

EVP:

 

12−156λ

 

−6L +22Lλ

 

=0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−6L +

22Lλ

2

2

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

4L

−4L

 

 

 

 

 

 

 

 

 

 

 

 

 

in which λ = ω 2 ρAL4 / 420 EI .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Solving the EVP, we obtain,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

 

12

v2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

= 1.38

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω1 = 3.533

 

 

 

 

 

 

 

 

 

#2

 

 

 

 

 

 

 

 

 

 

ρAL

 

 

 

 

 

 

θ2

 

 

 

 

 

 

#3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

1

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω2

=

 

 

 

 

 

 

,

 

 

 

=

 

 

 

.

#1

 

 

 

 

 

 

34.81

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρAL4

 

 

 

 

 

θ

2

7.62

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

L

Exact solutions:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

 

12

 

 

 

 

 

 

 

 

EI

 

 

12

 

 

 

 

 

 

ω1 = 3.516

 

 

 

 

 

 

 

,

 

ω2 = 22.03

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρAL4

 

 

 

 

 

 

 

ρAL4

 

 

 

 

 

 

 

We can see that mode 1 is calculated much more accurately than mode 2, with one beam element.

© 1997-2003 Yijun Liu, University of Cincinnati

166

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]