REFINING COMPLEX FREQUENCIES OF VIBRATING STRUCTURES
DOI:
https://doi.org/10.26512/ripe.v2i20.15003Abstract
In structural analysis, and related, the computation of natural frequencies of vibrating systems given by its mass, damping and stiffness matrices, when at least one of this matrices is dense and the matrix size is not prohibitively large, is usually done through reduction to a generalized first order system. The QZ iteration method is then applied to achieve the Generalized Schur form, from where all the natural frequencies, and associated mode shape vectors, can be
taken from, as eigenvalues and eigenvectors of some mathematical problem. As a drawback, this reduction not only increases by a factor of 2 the sizes of the working matrices, but also allows important properties like positive definiteness (or semi-definiteness), and even symmetry, to be lost. Even further, it must be considered the fact that the condition number of the new problem (as an upper bound for the ratio between changes in the solution and changes in the data, by
means of matrix and vector norms) usually increases, sometimes by a large amount. This work investigates how the final computation of complex-valued system eigenvalues (and then complex system frequencies) can take advantage of an iterative refinement technique, based in a shifted and inverted Krylov-type strategy that uses the original linear second-order matrices, in the same computational environment. The proposed strategy takes advantage that, in its original second-order form, the eigenvalue problem seams to be far better conditioned and well suited for high performance computations than in the generalized first-order form. The proposed strategy is applied to some test matrices from the Harwell-Boeing Collection of the Matrix Market
website. Numerical examples using professional software are also provided.
Keywords: second-order, complex frequency, vibration, structure
Downloads
References
E. Anderson, et al., 1999. LAPACK User’s Guide. Third Edition, SIAM.
Z. Bai and Y. Su, 2005. SOAR: A second-order Arnoldi method for the solution of the quadratic eigenvalue problem, SIAM J. Matrix Anal. Appl., vol. 26, n.3, 640”“659 , DOI:10.1137/S0895479803438523.
J. Carvalho and J. Claeyssen, 2011. Computation of extreme clustered natural frequencies of damped second-order linear systems, Proceedings of Conferˆencia Brasileira de Dinˆamica, Controle e Aplicac¸ ˜oes, vol. 1, 454”“457.
J. Carvalho and J. Claeyssen, 2013. Refining natural frequencies of damped linear vibrating structures, Proceedings of 34th Ibero-Latin American Congress on Computational Methods in Engineering.
J. Claeyssen, 1990. On predicting the response of nonconservative linear vibrating systems by using dynamical matrix solutions, J. Sound and Vibration, vol. 140, 73”“84.
S. Eisenstat and I. C. Ipsen, 1998. Three absolute perturbation bounds for matrix eigenvalues imply relative bounds. SIAM J. on Matrix Analysis and Appl., vol. 20, n. 1, 149”“158, DOI:10.1137/S0895479897323282.
G. Golub and C. Van Loan, 1996. Matrix computations, Johns Hopkins University Press, Baltimore, MD, USA.
D. Higham and N. Higham, 1999. Structured Backward Error and condition of generalized eigenvalue problems, SIAM J. Matrix Anal. Appl, vol. 20, n.2, 493”“512, DOI:10.1137/S0895479896313188.
G. Sleijpen and H. VanderVorst, 1996. A Jacobi-Davidson iteration method for linear eigenvalue problems, SIAM J. Matrix Anal. Appl., vol. 17, 401”“425, DOI:10.1137/S0036144599363084.
A. Sondipon, 1999. Rates of change of eigenvalues and eigenvectors in damped dynamic system, AIAA Journal, vol. 37, n.11, 1452”“1458, DOI: 10.2514/2.622.
G. W. Stewart, 1972. On the Sensitivity of the Eigenvalue Problem Ax = Bx. SIAM J. on Numerical Analysis vol. 9, n.4, 669”“686, DOI:10.1137/0709056.
A. Varga, 1990. Computation of irreducible generalized state-space realizations, Kybernetika, vol. 26, 89”“106, DOI: 10.1.1.539.2629.
Y. Zhou, 2006. Studies on Jacobi-Davidson, Rayleigh quotient iteration, inverse iteration generalized Davidson and Newton updates, Numerical Linear Algebra with Applications, vol.13, 621”“642, DOI: 10.1002/nla.490.
The Matrix Market directory on the web. Online reference and repository at http://math.nist.gov/MatrixMarket.
Downloads
Published
How to Cite
Issue
Section
License
Given the public access policy of the journal, the use of the published texts is free, with the obligation of recognizing the original authorship and the first publication in this journal. The authors of the published contributions are entirely and exclusively responsible for their contents.
1. The authors authorize the publication of the article in this journal.
2. The authors guarantee that the contribution is original, and take full responsibility for its content in case of impugnation by third parties.
3. The authors guarantee that the contribution is not under evaluation in another journal.
4. The authors keep the copyright and convey to the journal the right of first publication, the work being licensed under a Creative Commons Attribution License-BY.
5. The authors are allowed and stimulated to publicize and distribute their work on-line after the publication in the journal.
6. The authors of the approved works authorize the journal to distribute their content, after publication, for reproduction in content indexes, virtual libraries and similars.
7. The editors reserve the right to make adjustments to the text and to adequate the article to the editorial rules of the journal.
