A Numerical Method to Solve Nonsymmetric Eigensystems Applied to Dynamics of Turbomachinery

2020 ◽  
Vol 17 (09) ◽  
pp. 1950073
Author(s):  
Alfredo R. de Faria ◽  
Omair Alhatim ◽  
Homero Fonseca Santiago Maciel

In this paper, a canonical transformation is proposed to solve the eigenvalue problem related to the dynamics of rotor-bearing systems. In this problem, all matrices are real, but they may not be symmetric, which leads to the appearance of complex eigenvalues and eigenvectors. The bi-iteration method is selected to solve the original eigenproblem whereas the QR algorithm is adopted to solve the reduced or projected problem. A new canonical transformation of the global eigenproblem which reduces the quadratic eigenproblem to a linear eigenproblem, maintaining numerical stability since all that is required is that the stiffness matrix is well-conditioned, which is always true when it comes to applications in dynamic problems. The proposed technique is good for obtaining dominant eigenvalues and corresponding eigenvectors of real nonsymmetric matrices and it possesses the following properties: (i) the matrix is not transformed, therefore sparsity is maintained, (ii) partial eigensolutions can be obtained and (iii) use may be made of good eigenvectors predictions.

2006 ◽  
Vol 11 (1) ◽  
pp. 13-32 ◽  
Author(s):  
B. Bandyrskii ◽  
I. Lazurchak ◽  
V. Makarov ◽  
M. Sapagovas

The paper deals with numerical methods for eigenvalue problem for the second order ordinary differential operator with variable coefficient subject to nonlocal integral condition. FD-method (functional-discrete method) is derived and analyzed for calculating of eigenvalues, particulary complex eigenvalues. The convergence of FD-method is proved. Finally numerical procedures are suggested and computational results are schown.


1994 ◽  
Vol 77 (5) ◽  
pp. 2481-2495 ◽  
Author(s):  
M. Sammon

A multivariate model is outlined for a distributed respiratory central pattern generator (RCPG) and its afferent control. Oscillatory behavior of the system depends on structure and symmetry of a matrix of phase-switching functions (F omega, phi) that control distribution of central excitation (CE) and inhibition (CI) within the circuit. The matrix diagonal (F omega) controls activation of CI variables as excitatory inputs are altered (e.g., central and afferent contributions to inspiratory off switch); off-diagonal terms (F phi) distribute excitations within the CI system and produce complex eigenvalues at the switching points between inspiration and expiration. For null F phi, phase switchings of saddle equilibria located at end expiration and end inspiration are overdamped all-or-nothing events; graded control of CI is seen for phi > 0. When coupling is significant (phi >> 0), CI dynamics become underdamped, admitting a domain of inputs where chaotic behavior is predictably observed. For the homogeneous RCPG (symmetric F omega, phi), CE oscillations are one-dimensional limit cycles (D = 1) or weakly chaotic (D approximately equal to 1). When perturbations from symmetry are significant, the distributed RCPG becomes partitioned where strongly chaotic oscillations (D > or = 2) and central apnea (D = 0) are seen more frequently. The equations provide means for mapping Silnikov bifurcations that alter the geometry and dimension of the breathing pattern and formalisms for discussing RCPG processing of afferent information.


Author(s):  
Nikta Shayanfar ◽  
Heike Fassbender

The polynomial eigenvalue problem is to find the eigenpair of $(\lambda,x) \in \mathbb{C}\bigcup \{\infty\} \times \mathbb{C}^n \backslash \{0\}$ that satisfies $P(\lambda)x=0$, where $P(\lambda)=\sum_{i=0}^s P_i \lambda ^i$ is an $n\times n$ so-called matrix polynomial of degree $s$, where the coefficients $P_i, i=0,\cdots,s$, are $n\times n$ constant matrices, and $P_s$ is supposed to be nonzero. These eigenvalue problems arise from a variety of physical applications including acoustic structural coupled systems, fluid mechanics, multiple input multiple output systems in control theory, signal processing, and constrained least square problems. Most numerical approaches to solving such eigenvalue problems proceed by linearizing the matrix polynomial into a matrix pencil of larger size. Such methods convert the eigenvalue problem into a well-studied linear eigenvalue problem, and meanwhile, exploit and preserve the structure and properties of the original eigenvalue problem. The linearizations have been extensively studied with respect to the basis that the matrix polynomial is expressed in. If the matrix polynomial is expressed in a special basis, then it is desirable that its linearization be also expressed in the same basis. The reason is due to the fact that changing the given basis ought to be avoided \cite{H1}. The authors in \cite{ACL} have constructed linearization for different bases such as degree-graded ones (including monomial, Newton and Pochhammer basis), Bernstein and Lagrange basis. This contribution is concerned with polynomial eigenvalue problems in which the matrix polynomial is expressed in Hermite basis. In fact, Hermite basis is used for presenting matrix polynomials designed for matching a series of points and function derivatives at the prescribed nodes. In the literature, the linearizations of matrix polynomials of degree $s$, expressed in Hermite basis, consist of matrix pencils with $s+2$ blocks of size $n \times n$. In other words, additional eigenvalues at infinity had to be introduced, see e.g. \cite{CSAG}. In this research, we try to overcome this difficulty by reducing the size of linearization. The reduction scheme presented will gradually reduce the linearization to its minimal size making use of ideas from \cite{VMM1}. More precisely, for $n \times n$ matrix polynomials of degree $s$, we present linearizations of smaller size, consisting of $s+1$ and $s$ blocks of $n \times n$ matrices. The structure of the eigenvectors is also discussed.


2020 ◽  
Vol 10 (01) ◽  
pp. 2150013 ◽  
Author(s):  
Guillaume Dubach

We establish a few properties of eigenvalues and eigenvectors of the quaternionic Ginibre ensemble (QGE), analogous to what is known in the complex Ginibre case (see [7, 11, 14]). We first recover a version of Kostlan’s theorem that was already at the heart of an argument by Rider [1], namely, that the set of the squared radii of the eigenvalues is distributed as a set of independent gamma variables. Our proof technique uses the De Bruijn identity and properties of Pfaffians; it also allows to prove that the high powers of these eigenvalues are independent. These results extend to any potential beyond the Gaussian case, as long as radial symmetry holds; this includes for instance truncations of quaternionic unitary matrices, products of quaternionic Ginibre matrices, and the quaternionic spherical ensemble. We then study the eigenvectors of quaternionic Ginibre matrices. Angles between eigenvectors and the matrix of overlaps both exhibit some specific features that can be compared to the complex case. In particular, we compute the distribution and the limit of the diagonal overlap associated to an eigenvalue that is conditioned to be at the origin. This complements a recent study of overlaps in quaternionic ensembles by Akemann, Förster and Kieburg [1, 2].


2018 ◽  
Vol 931 ◽  
pp. 72-77
Author(s):  
Leonid N. Panasyuk ◽  
Galina M. Kravchenko ◽  
Vakhtang P. Matua

The article considers the modeling of dynamic processes in buildings and structures. A general formulation of the dynamic problem of a massive load motion on a massive structure is considered. The equation of motion is obtained in the form of a finite element method. The equations solving is performed using direct methods of integrating dynamic problems. Absolutely stable schemes of direct integration are constructed, where the system of solving equations is trivial and the matrix of the system is diagonal. Due to this, the complexity at the time step is as low as in explicit schemes. Therefore, the proposed methods can be considered as explicit absolutely stable schemes of direct integration of a dynamical problem with a variable in time mass. These recommendations are for estimating the accuracy of a numerical solution.


Robotica ◽  
2016 ◽  
Vol 35 (6) ◽  
pp. 1310-1326 ◽  
Author(s):  
Guanglei Wu ◽  
Ping Zou

SUMMARYThis paper deals with the stiffness modeling, analysis and comparison of a Biglide parallel grinder with two alternative modular parallelograms. It turns out that the Cartesian stiffness matrix of the manipulator has the property that it can be decoupled into two homogeneous matrices, corresponding to the translational and rotational aspects, through which the principal stiffnesses and the associated directions are identified by means of the eigenvalue problem, allowing the evaluation of the translational and rotational stiffness of the manipulator either at a given pose or the overall workspace. The stiffness comparison of the two alternative Biglide machines reveals the (dis)advantages of the two different spatial modular parallelograms.


1983 ◽  
Vol 66 ◽  
pp. 331-341
Author(s):  
M. Knölker ◽  
M. Stix

AbstractThe differential equations describing stellar oscillations are transformed into an algebraic eigenvalue problem. Frequencies of adiabatic oscillations are obtained as the eigenvalues of a banded real symmetric matrix. We employ the Cowling-approximation, i.e. neglect the Eulerian perturbation of the gravitational potential, and, in order to preserve selfadjointness, require that the Eulerian pressure perturbation vanishes at the outer boundary. For a solar model, comparison of first results with results obtained from a Henyey method shows that the matrix method is convenient, accurate, and fast.


Sign in / Sign up

Export Citation Format

Share Document