Block-partitioned Rayleigh-Ritz method for efficient eigenpair reanalysis of large-scale finite element models

被引:0
作者
Jeong, Yeon-Ho [1 ]
Boo, Seung-Hwan [1 ]
Yim, Solomon C. [2 ]
机构
[1] Korea Maritime & Ocean Univ, Div Naval Architecture & Ocean Syst Engn, 727 Taejong Ro, Busan 49112, South Korea
[2] Oregon State Univ, Sch Civil & Construct Engn, Corvallis, OR 97331 USA
基金
新加坡国家研究基金会;
关键词
structural analysis; eigenvalue problem; Rayleigh-Ritz method; approximate reanalysis; finite element analysis; block-partitioning; COMBINED APPROXIMATIONS; DYNAMIC REANALYSIS; ALGORITHM; PERTURBATIONS;
D O I
10.1093/jcde/qwad030
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In this manuscript, we propose a new effective method for eigenpair reanalysis of large-scale finite element (FE) models. Our method utilizes the matrix block-partitioning algorithm in the Rayleigh-Ritz approach and expresses the Ritz basis matrix using thousands of block matrices of very small size. To avoid significant computational costs from the projection procedure, we derive a new formulation that uses tiny block computations instead of global matrix computations. Additionally, we present an algorithm that recognizes which blocks are changed in the modified FE model to achieve computational cost savings when computing new eigenpairs. Through selective updating for the recognized blocks, we can effectively construct the new Ritz basis matrix and the new reduced mass and stiffness matrices corresponding to the modified FE model. To demonstrate the performance of our proposed method, we solve several practical engineering problems and compare the results with those of the combined approximation method, the most well-known eigenpair reanalysis method, and ARPACK, an eigenvalue solver embedded in many numerical programs.
引用
收藏
页码:959 / 978
页数:20
相关论文
共 56 条
[1]   Density-based shape optimization for fail-safe design [J].
Ambrozkiewicz, Olaf ;
Kriegesmann, Benedikt .
JOURNAL OF COMPUTATIONAL DESIGN AND ENGINEERING, 2020, 7 (05) :615-629
[2]   AN ACCELERATED SUBSPACE ITERATION METHOD [J].
BATHE, KJ ;
RAMASWAMY, S .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 1980, 23 (03) :313-331
[3]  
Bathe KJ, 2006, Finite Element Procedures
[4]   An estimator for the diagonal of a matrix [J].
Bekas, C. ;
Kokiopoulou, E. ;
Saad, Y. .
APPLIED NUMERICAL MATHEMATICS, 2007, 57 (11-12) :1214-1229
[5]   Computing charge densities with partially reorthogonalized Lanczos [J].
Bekas, C ;
Saad, Y ;
Tiago, ML ;
Chelikowsky, JR .
COMPUTER PHYSICS COMMUNICATIONS, 2005, 171 (03) :175-186
[6]   An automated multilevel substructuring method for eigenspace computation in linear elastodynamics [J].
Bennighof, JK ;
Lehoucq, RB .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2004, 25 (06) :2084-2106
[8]  
Buluç A, 2009, SPAA'09: PROCEEDINGS OF THE TWENTY-FIRST ANNUAL SYMPOSIUM ON PARALLELISM IN ALGORITHMS AND ARCHITECTURES, P233
[9]  
Bunch J. R., 2014, Sparse Matrix Computations
[10]  
Calvetti D., 1994, Electron. Trans. Numer. Anal., V2, P21