1 / 20

Eigensolution Method for Structures Using Accelerated Lanczos Algorithm

2002 한국지진공학회 추계학술발표회. 2002 년 9 월 28 일. Eigensolution Method for Structures Using Accelerated Lanczos Algorithm. Byoung-Wan Kim 1) , Ju-Won Oh 2) and In-Won Lee 3) 1) Graduate Student, Dept. of Civil and Environmental Eng., KAIST

dominique-dudley
Télécharger la présentation

Eigensolution Method for Structures Using Accelerated Lanczos Algorithm

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 2002한국지진공학회 추계학술발표회 2002년 9월 28일 Eigensolution Method for Structures Using Accelerated Lanczos Algorithm Byoung-Wan Kim1), Ju-Won Oh2) and In-Won Lee3) 1) Graduate Student, Dept. of Civil and Environmental Eng., KAIST 2) Professor, Dept. of Civil and Environmental Eng., Hannam Univ. 3) Professor, Dept. of Civil and Environmental Eng., KAIST

  2. Contents • Introduction • Matrix-powered Lanczos method • Numerical examples • Conclusions

  3. Introduction • Background • Dynamic analysis of structures • - Direct integration method • - Mode superposition method  Eigenvalue analysis • Eigenvalue analysis • - Subspace iteration method • - Determinant search method • - Lanczos method • The Lanczos method is very efficient.

  4. Literature review • The Lanczos method was first proposed in 1950. • Erricson and Ruhe (1980): Lanczos algorithm with shifting • Smith et al. (1993): Implicitly restarted Lanczos algorithm • In the fields of quantum physics, Grosso et al. (1993) • modified Lanczos recursion to accelerate convergence

  5. Objective • Application of Lanczos method using the power technique • to the eigenproblem of structures in structural dynamics •  Matrix-powered Lanczos method

  6. Matrix-powered Lanczos method • Eigenproblem of structure

  7. Modified Gram-Schmidt process of Krylov sequence

  8. Modified Lanczos recursion

  9. Reduced tridiagonal standard eigenproblem

  10. Summary of algorithm and operation count n = order of M and K, m = half-bandwidth of M and K q = the number of calculated Lanczos vectors or order of T sj = the number of iterations of jth step in QR iteration

  11. Numerical examples • Structures • Simple spring-mass system (Chen 1993) • Plan framed structure (Bathe and Wilson 1972) • Three-dimensional frame structure (Bathe and Wilson 1972) • Three-dimensional building frame (Kim and Lee 1999) • Physical error norm (Bathe 1996)

  12. Simple spring-mass system (DOFs: 100) • System matrices

  13. Number of operations  Failure in convergence due to numerical instability of high matrix power

  14. Plane framed structure (DOFs: 330) • Geometry and properties A = 0.2787 m2 I = 8.63110-3m4 E = 2.068107Pa  = 5.154102kg/m3

  15. Number of operations  Failure in convergence due to numerical instability of high matrix power

  16. Three-dimensional frame structure (DOFs: 468) • Geometry and properties E = 2.068107Pa  = 5.154102kg/m3 Column in front building : A = 0.2787 m2, I = 8.63110-3m4 Column in rear building : A = 0.3716 m2, I = 10.78910-3m4 All beams into x-direction : A = 0.1858 m2, I = 6.47310-3m4 All beams into y-direction : A = 0.2787 m2, I = 8.63110-3m4

  17. Number of operations

  18. Three-dimensional building frame (DOFs: 1008) • Geometry and properties A = 0.01 m2 I = 8.310-6m4 E = 2.11011Pa  = 7850 kg/m3

  19. Number of operations  Failure in convergence due to numerical instability of high matrix power

  20. Conclusions • Matrix-powered Lanczos method has not only the better convergence but also the less operation count than the conventional Lanczos method. • The suitable power of the dynamic matrix that gives numerically stable solution in the matrix-powered Lanczos method is the second power.

More Related