Large sparse second order index-1 descriptor systems arise in various disciplines of science and engineering, such as constraint mechanics or multibody dynamics, mechatronics (where mechanical and electrical elements ...
详细信息
ISBN:
(纸本)9781509040919
Large sparse second order index-1 descriptor systems arise in various disciplines of science and engineering, such as constraint mechanics or multibody dynamics, mechatronics (where mechanical and electrical elements are coupled), but also RLC circuit design. Simulation, controller design and design optimization are only some applications of such models. Either of these tasks, just like any other many-query situation becomes unfeasible when the system is high dimensional. This paper discusses an algorithm to obtain a reduced state space model of a large sparse second order index-1 system using an interpolatory projection method based on the iterative rational Krylov algorithm (IRKA). In each iteration of this algorithm, we need to solve a number of linear systems. The main contribution of this paper is to solve these linear systems by exploiting the sparsity of the original model, which reduces the computational cost drastically. The algorithm is applied to a micro-mechanical piezo-actuated structural FEM model of a certain building block of a machine tool. Numerical experiments with a complex 3d model of an adaptive spindle support (a piezo-mechanical multiphysics system) show the effectivity and efficiency of the techniques.
Free boundary and moving boundary problems, that can be used to model crystal growth or the solidification and melting of pure materials, receive growing attention in science and technology. The optimal control of the...
We investigate the numerical solution to a low rank perturbed Lyapunov equation A T X + XA = W via the sign function method (SFM). The sign function method has been proposed to solve Lyapunov equations, see e.g. [1], ...
Differential matrix equations appear in many applications like optimal control of partial differential equations, balanced truncation model order reduction of linear time varying systems and many more. Here, we will f...
We study the numerical solution of control constrained optimal control problems governed by a system of convection diffusion equations with nonlinear reaction terms, arising from chemical processes. control constraint...
详细信息
The spectral abscissa is a fundamental map from the set of complex matrices to the real numbers. Denoted α and defined as the maximum of the real parts of the eigenvalues of a matrix X, it has many applications in st...
详细信息
We propose a new implementation of the sign function based spectral divide-and-conquer method for the generalized non-symmetric eigenvalue problem. The basic idea is to use the generalized matrix sign function to spli...
Linear time-periodic systems arise whenever a nonlinear system is linearized about a periodic trajectory. Stability of the solution may be proven by rigorous bounds on the solution. The key idea of this paper is to de...
The M.E.S.S. software suite for solving large scale matrix equations and related problems is the successor of the obsolete LyaPack MATLAB® toolbox. The software suite consists of a new MATLAB toolbox and a separa...
暂无评论