Structure preserving model order reduction of large sparse second-order index-1 systems and application to a mechatronics model. Issue 6 (1st November 2016)

• Record Type: Journal Article
• Authors: Benner, Peter; Saak, Jens; Uddin, M. Monir
• Is Part Of: Mathematical and computer modelling of dynamical systems. Volume 22:Issue 6(2016)
• Page Start: 509
• Page End: 523
• Publication Date: 2016-11-01
• Subjects: Model-order reduction -- mechatronic system -- balanced truncation -- Lyapunov equations -- large-scale -- differential algebraic system; Engineering -- Mathematical models -- Periodicals; Computer simulation -- Periodicals; 515.39
• DOI: 10.1080/13873954.2016.1218347
• Languages: English
• ISSNs: 1387-3954

Abstract

ABSTRACT: Nowadays, mechanical engineers heavily depend on mathematical models for simulation, optimization and controller design. In either of these tasks, reduced dimensional formulations are obligatory in order to achieve fast and accurate results. Usually, the structural mechanical systems of machine tools are described by systems of second-order differential equations. However, they become descriptor systems when extra constraints are imposed on the systems. This article discusses efficient techniques of Gramian-based model-order reduction for second-order index-1 descriptor systems. Unlike, our previous work, here we mainly focus on a second-order to second-order reduction technique for such systems, where the stability of the system is guaranteed to be preserved in contrast to the previous approaches. We show that a special choice of the first-order reformulation of the system allows us to solve only one Lyapuov equation instead of two. We also discuss improvements of the technique to solve the Lyapunov equation using low-rank alternating direction implicit methods, which further reduces the computational cost as well as memory requirement. The proposed technique is applied to a structural finite element method model of a micro-mechanical piezo-actuators-based adaptive spindle support. Numerical results illustrate the increased efficiency of the adapted method.