Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
⎡ ⎤ ⎧ ⎫ ⎡<br />
⎤ ⎧ ⎫ ⎧ ⎫<br />
⎣ m 1 0 ⎨x ⎦<br />
′′ ⎬<br />
1<br />
0 m<br />
⎩<br />
2 x ′′ ⎭ + ⎣ k 1 + k 2 −k 2<br />
⎨x ⎦ 1 ⎬ ⎨<br />
2 −k 2 k<br />
⎩<br />
2 x<br />
⎭ = f 1 (t) ⎬<br />
⎩<br />
2 f 2 (t)<br />
⎭<br />
⎡ ⎤ ⎧ ⎫ ⎡ ⎤ ⎧ ⎫ ⎧ ⎫<br />
⎣ 1 0 ⎨x ⎦<br />
′′ ⎬<br />
1<br />
0 3<br />
⎩<br />
x ′′ ⎭ + ⎣ 3 −2 ⎨x ⎦ 1 ⎬ ⎨<br />
−2 2<br />
⎩<br />
x<br />
⎭ = 0 ⎬<br />
⎩<br />
2 sin (5t)<br />
⎭<br />
2<br />
In this example m 11 = 1, m 12 = 0, m 21 = 0, m 22 = 3 and k 11 = 3, k 12 = −2, k 21 = −2, k 22 = 2 and f 1 (t) = 0<br />
and f 2 (t) = sin (5t)<br />
step 2 is now applied which solves the eigenvalue problem in order to find the two natural frequencies<br />
Let ω 2 = λ hence<br />
The solution is λ 1 = 3. 475 and λ 2 = 0.192, therefore<br />
And<br />
det ( [K] − ω 2 [M] ) = 0<br />
⎛⎡<br />
⎤ ⎡ ⎤⎞<br />
det ⎝⎣ 3 −2 ⎦ − ω 2 ⎣ 1 0 ⎦⎠ = 0<br />
−2 2 0 3<br />
⎡<br />
⎤<br />
det ⎣ 3 − ω2 −2<br />
⎦ = 0<br />
−2 2 − 3ω 2<br />
(<br />
3 − ω<br />
2 ) ( 2 − 3ω 2) − (−2) (−2) = 0<br />
3ω 4 − 11ω 2 + 2 = 0<br />
3λ 2 − 11λ + 2 = 0<br />
ω 1 = √ 3.475 = 1.864<br />
ω 2 = √ 0.192 = 0.438<br />
step 3 is now applied which finds the non-normalized eigenvectors. For each natural frequency ω 1 and ω 2 the<br />
corresponding shape function is found by solving the following two sets of equations for the eigen vectors ϕ 1 , ϕ 2<br />
⎛⎡<br />
⎤ ⎡ ⎤⎞<br />
⎧ ⎫ ⎧ ⎫<br />
⎝⎣ 3 −2 ⎦ − ω1<br />
2 ⎣ 1 0 ⎨ϕ ⎦⎠<br />
11 ⎬ ⎨<br />
−2 2 0 3<br />
⎩<br />
ϕ<br />
⎭ = 0⎬<br />
⎩<br />
21 0<br />
⎭<br />
For ω 1 = 1. 864<br />
⎛⎡<br />
⎤ ⎡ ⎤⎞<br />
⎧ ⎫ ⎧ ⎫<br />
⎝⎣ 3 −2 ⎦ − 1.864 2 ⎣ 1 0 ⎨ϕ ⎦⎠<br />
11 ⎬ ⎨<br />
−2 2<br />
0 3<br />
⎩<br />
ϕ<br />
⎭ = 0⎬<br />
⎩<br />
21 0<br />
⎭<br />
⎡<br />
⎤ ⎧ ⎫ ⎧ ⎫<br />
⎣ −0.475 −2 ⎨ 1 ⎬ ⎨<br />
⎦<br />
−2 −8.424<br />
⎩<br />
ϕ<br />
⎭ = 0⎬<br />
⎩<br />
21 0<br />
⎭<br />
This gives one equation to solve for ϕ 21 (the first row equation is only used)<br />
−0.475 − 2ϕ 21 = 0<br />
10