We present a formal procedure for structure-preserving model reduction of linear second-order and Hamiltonian control problems that appear in a variety of physical contexts, e.g., vibromechanical systems or electrical circuit design. Typical balanced truncation methods that project onto the subspace of the largest Hankel singular values fail to preserve the problem's physical structure and may suffer from lack of stability. In this paper, we adopt the framework of generalized Hamiltonian systems that covers the class of relevant problems and that allows for a generalization of balanced truncation to second-order problems. It turns out that the Hamiltonian structure, stability, and passivity are preserved if the truncation is done by imposing a holonomic constraint on the system rather than standard Galerkin projection. Abstract. We present a formal procedure for structure-preserving model reduction of linear second-order and Hamiltonian control problems that appear in a variety of physical contexts, e.g., vibromechanical systems or electrical circuit design. Typical balanced truncation methods that project onto the subspace of the largest Hankel singular values fail to preserve the problem's physical structure and may suffer from lack of stability. In this paper, we adopt the framework of generalized Hamiltonian systems that covers the class of relevant problems and that allows for a generalization of balanced truncation to second-order problems. It turns out that the Hamiltonian structure, stability, and passivity are preserved if the truncation is done by imposing a holonomic constraint on the system rather than standard Galerkin projection.Key words. structure-preserving model reduction, generalized Hamiltonian systems, balanced truncation, strong confinement, invariant manifolds, Hankel norm approximation AMS subject classifications. 70Q05, 70J50, 93B11, 93C05, 93C70.
DOI. 10.1137/0807327171. Introduction. Model reduction is a major issue for control, optimization, and simulation of large-scale systems. In particular, for linear time-invariant systems, balanced truncation is a well-established tool for deriving reduced (i.e., low-dimensional) models that have an input-output behavior similar to the original model [1, 2]; see also [3] and the references therein. The general idea of balanced truncation is to restrict the system onto the subspace of easily controllable and observable states which can be determined by the computing Hankel singular values associated with the system. Moreover the method is known to preserve certain properties of the original system such as stability or passivity and gives an error bound that is easily computable. One drawback of balanced truncation is that there is no straightforward generalization to second-order systems; see, e.g., [4] or the recent articles [5,6,7,8] for a discussion of various possible strategies. Second-order equations occur in modeling and control of many physical systems, e.g., electrical circuits, structural mechanics, or vibroacoustic models (see...