English

On the shift-invert Lanczos method for the buckling eigenvalue problem

Numerical Analysis 2019-10-22 v1 Numerical Analysis

Abstract

We consider the problem of extracting a few desired eigenpairs of the buckling eigenvalue problem Kx=λKGxKx = \lambda K_Gx, where KK is symmetric positive semi-definite, KGK_G is symmetric indefinite, and the pencil KλKGK - \lambda K_G is singular, namely, KK and KGK_G share a non-trivial common nullspace. Moreover, in practical buckling analysis of structures, bases for the nullspace of KK and the common nullspace of KK and KGK_G are available. There are two open issues for developing an industrial strength shift-invert Lanczos method: (1) the shift-invert operator (KσKG)1(K - \sigma K_G)^{-1} does not exist or is extremely ill-conditioned, and (2) the use of the semi-inner product induced by KK drives the Lanczos vectors rapidly towards the nullspace of KK, which leads to a rapid growth of the Lanczos vectors in norms and cause permanent loss of information and the failure of the method. In this paper, we address these two issues by proposing a generalized buckling spectral transformation of the singular pencil KλKGK - \lambda K_G and a regularization of the inner product via a low-rank updating of the semi-positive definiteness of KK. The efficacy of our approach is demonstrated by numerical examples, including one from industrial buckling analysis.

Cite

@article{arxiv.1910.08652,
  title  = {On the shift-invert Lanczos method for the buckling eigenvalue problem},
  author = {Chao-Ping Lin and Huiqing Xie and Roger Grimes and Zhaojun Bai},
  journal= {arXiv preprint arXiv:1910.08652},
  year   = {2019}
}
R2 v1 2026-06-23T11:48:18.650Z