Browse > Article
http://dx.doi.org/10.3795/KSME-A.2003.27.4.551

Improvement of Computational Efficiency of the Subspace Iteration Method for Large Finite Element Models  

Joo, Byung-Hyun (한국과학기술원 기계공학과)
Lee, Byung-Chai (한국과학기술원 기계공학과)
Publication Information
Transactions of the Korean Society of Mechanical Engineers A / v.27, no.4, 2003 , pp. 551-558 More about this Journal
Abstract
An efficient and reliable subspace iteration algorithm using the block algorithm is proposed. The block algorithm is the method dividing eigenpairs into several blocks when a lot of eigenpairs are required. One of the key for the faster convergence is carefully selected initial vectors. As the initial vectors, the proposed method uses the modified Ritz vectors for guaranteering all the required eigenpairs and the quasi-static Ritz vectors for accelerating convergency of high frequency eigenvectors. Applying the quasi-static Ritz vectors, a shift is always required, and the proper shift based on the geometric average is proposed. To maximize efficiency, this paper estimates the proper number of blocks based on the theoretical amount of calculation in the subspace iteration. And it also considers the problems generated in the process of combining various algorithms and the solutions to the problems. Several numerical experiments show that the proposed subspace iteration algorithm is very efficient, reliable ,and accurate.
Keywords
Subspace Iteration Method; Large Finite Element Model; Eigenvalue Problem;
Citations & Related Records
연도 인용수 순위
  • Reference
1 Bathe, K. J., 1982, Finite Element Procedures in Engineering Analysis, Prentice-Hall, Englewood Cliffs, New Jersey
2 Bertolini, A. F. and Lam, Y. C., 1998, 'Accelerated Subspace Iteration Using Adaptive Multiple Inverse Iteration,' Computers & Structures, Vol. 66, pp. 45-57   DOI   ScienceOn
3 Qian, Y. and Dhatt, G., 1995, 'An Accelerated Subspace Method for Generalized Eigenproblems,' Computers & Structures, Vol. 54, pp. 1127-1134   DOI   ScienceOn
4 Jennings, A. and Stewart, W. J., 1981, 'A Simultaneous Iteration Algorithm for Real Matrices,' ACM Transactions on Mathematical Software, Vol. 7, pp. 184-198   DOI   ScienceOn
5 Akl, F. A., 1982, 'Acceleration of Subspace Iteration,' International Journal of Numerical Methods in Engineering, Vol. 18, pp. 583-589   DOI   ScienceOn
6 Bathe, K. J. and Ramaswamy, S., 1980, 'An Acceleration Subspace Iteration Method,' Computer Methods in Applied Mechanics and Engineering, Vol. 23, pp. 313-331   DOI   ScienceOn
7 Gu, J., Ma, Z-D. and G. M. Hulbert, 2000, 'A New Load-Dependent Ritz Vector Method for Structural Dynamics Analyses: Quasi-Static Ritz Vectors,' Finite Elements in Analysis and Design, Vol. 36, pp. 261-278   DOI   ScienceOn
8 Jiaxian, X., 1989, 'An Improved Method for Partial Eigensolution of Large Structures,' Computers & Structures, Vol. 32, pp, 1055-1060   DOI   ScienceOn
9 Wilson, E. L. and Itoh, T., 1983, 'An Eigensolution Strategy for Large Systems,' Computers & Structures, Vol. 16, pp259-265   DOI   ScienceOn
10 Trefethen, L. N. and Bau, D., 1997, Numerical Linear Algebra, SIAM, Philadelphia