CA-AA242B-Ch6

CA-AA242B-Ch6 - AA242B: MECHANICAL VIBRATIONS AA242B:...

Info iconThis preview shows pages 1–7. Sign up to view the full content.

View Full Document Right Arrow Icon
AA242B: MECHANICAL VIBRATIONS AA242B: MECHANICAL VIBRATIONS Solution Methods for the Generalized Eigenvalue Problem These slides are based on the recommended textbook: M. G´ eradin and D. Rixen, “Mechanical Sons, Incorporated, ISBN-13:9780471975465 AA242B: MECHANICAL VIBRATIONS
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
AA242B: MECHANICAL VIBRATIONS Outline 1 Eigenvector Iteration Methods 2 Subspace Construction Methods AA242B: MECHANICAL VIBRATIONS
Background image of page 2
AA242B: MECHANICAL VIBRATIONS Eigenvector Iteration Methods Introduction Generalized eigenvalue problem associated with an n -degree-of-freedom system Kx = ω 2 Mx n eigenvalues 0 ω 2 1 ω 2 2 ≤ ··· ≤ ω 2 n and n associated eigenvectors x 1 , x 2 , ··· , x n AA242B: MECHANICAL VIBRATIONS
Background image of page 3

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
AA242B: MECHANICAL VIBRATIONS Eigenvector Iteration Methods Introduction Dynamic fexibility matrix consider a system that has no rigid-body modes K is non-singular construct the dynamic fexibility matrix deFned as D = K - 1 M the associated eigenvalue problem is Kx = ω 2 Mx Dx = λ x with λ = 1 ω 2 eigenvectors are same as for the generalized eigenproblem Kx = ω 2 Mx eigenvalues are given by λ i = 1 ω 2 i λ 1 λ 2 ≥ ··· ≥ λ n drawbacks o± working with D building cost is O ( n 3 ) non-symmetric matrix and as such must be stored completely AA242B: MECHANICAL VIBRATIONS
Background image of page 4
AA242B: MECHANICAL VIBRATIONS Eigenvector Iteration Methods The Power Algorithm Determination of the fundamental eigenmode an iterative solution to Dx = λ x can be obtained by considering the iteration z p +1 = Dz p proof: expand z 0 in the basis of the eigenmodes of the system z 0 = n X i =1 α i x i Then, at the p -th iteration, z p = D p z 0 = n X i =1 α i D p x i = n X i =1 α i λ p i x i = λ p 1 α 1 x 1 + n X i =2 α i λ i λ 1 « p x i ! Recall that λ 1 λ 2 ≥ ··· ≥ λ n and assume that λ 1 > λ 2 = z p λ p 1 α 1 x 1 k z p +1 k k z p k λ 1 as p → ∞ AA242B: MECHANICAL VIBRATIONS
Background image of page 5

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
AA242B: MECHANICAL VIBRATIONS Eigenvector Iteration Methods The Power Algorithm Convergence analysis p , z p = λ p 1 x 1 + λ p 2 x 2 + ··· = z T p +1 e j z
Background image of page 6
Image of page 7
This is the end of the preview. Sign up to access the rest of the document.

This note was uploaded on 06/17/2010 for the course AA 242B taught by Professor Charbelfarhat during the Spring '10 term at Stanford.

Page1 / 17

CA-AA242B-Ch6 - AA242B: MECHANICAL VIBRATIONS AA242B:...

This preview shows document pages 1 - 7. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online