On the shift-invert Lanczos method for the buckling eigenvalue problem
Abstract
We consider the problem of extracting a few desired eigenpairs of the buckling eigenvalue problem , where is symmetric positive semi-definite, is symmetric indefinite, and the pencil is singular, namely, and share a non-trivial common nullspace. Moreover, in practical buckling analysis of structures, bases for the nullspace of and the common nullspace of and are available. There are two open issues for developing an industrial strength shift-invert Lanczos method: (1) the shift-invert operator does not exist or is extremely ill-conditioned, and (2) the use of the semi-inner product induced by drives the Lanczos vectors rapidly towards the nullspace of , 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 and a regularization of the inner product via a low-rank updating of the semi-positive definiteness of . 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}
}