• Record: found
  • Abstract: found
  • Poster: found
Is Open Access

Sampling-free parametric model reduction of systems with structured parameter variation

Read this article at

      There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.


      We consider a parametric linear time invariant dynamical systems represented in state-space form as\[E \dot x(t) = A(p) x(t) + Bu(t), \\ y(t) = Cx(t),\]where $E, A(p) \in \mathbb{R}^{n\times n}$, $B\in \mathbb{R}^{n\times m} \(and \)C\in \mathbb{R}^{l\times n}$. Here $x(t)\in \mathbb{R}^{n} \(denotes the state variable, while \)u(t)\in \mathbb{R}^{m}\(and \)y(t)\in \mathbb{R}^{l}\(represent, respectively, the inputs and outputs of the system.We assume that \)A(p)\(depends on \)k\ll n\(parameters \)p=(p_1, p_2, \ldots, p_k)\(such that we may write\[A(p)=A_0+U\,\diag (p_1, p_2, \ldots, p_k)V^T,\]where \)U, V \in \mathbb{R}^{n\times k}\(are given fixed matrices.We propose an approach for approximating the full-order transfer function \)H(s;p)=C(s E -A(p))^{-1}B\(with a reduced-order model that retains the structure of parametric dependence and (typically) offers uniformly high fidelity across the full parameter range. Remarkably, the proposed reduction process removes the need for parameter sampling and thus does not depend on identifying particular parameter values of interest. Our approach is based on the classic Sherman-Morrison-Woodbury formula and allows us to construct a parameterized reduced order model from transfer functions of four subsystems that do not depend on parameters, allowing one to apply well-established model reduction techniques for non-parametric systems. The overall process is well suited for computationally efficient parameter optimization and the study of important system properties.One of the main applications of our approach is for damping optimization: we consider a vibrational system described by\[ \begin{equation}\label{MDK} \begin{array}{rl} M\ddot q(t)+(C_{int} + C_{ext})\dot q(t)+Kq(t)&=E w(t),\\ z(t)&=Hq(t), \end{array} \end{equation} \]where the mass matrix, \)M$, and stiffness matrix, $K$, are real, symmetric positive-definite matrices of order $n$. Here, $q(t)\(is a vector of displacements and rotations, while \) w(t) \(and \)z(t) \(represent, respectively, the inputs (typically viewed as potentially disruptive) and outputs of the system. Damping in the structure is modeled as viscous damping determined by \)C_{int} + C_{ext}\(where \)C_{int}\(and \)C_{ext}\(represent contributions from internal and external damping, respectively. Information regarding damper geometry and positioning as well as the corresponding damping viscosities are encoded in \)C_{ext}= U\diag{(p_1, p_2, \ldots, p_k)} U^T\(where \)U \in \mathbb{R}^{n\times k}\(determines the placement and geometry of the external dampers.The main problem is to determine the best damping matrix that is able to minimize influence of the disturbances, \)w$, on the output of the system $z$. We use a minimization criteria based on the $\mathcal{H}_2$ system norm. In realistic settings, damping optimization is a very demanding problem. We find that the parametric model reduction approach described here offers a new tool with significant advantages for the efficient optimization of damping in such problems.

      Related collections

      Author and article information

      [1 ]Department of Mathematics, Virginia Polytechnic Institute and State University, Blacksburg
      [2 ]Department of Mathematics, University J.J. Strossmayer in Osijek, Osijek, Croatia
      [* ]Correspondence: ztomljan@
      ScienceOpen Posters
      27 April 2018
      Copyright © 2018

      This work has been published open access under Creative Commons Attribution License CC BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Conditions, terms of use and publishing policy can be found at


      Comment on this article