Application of operator inference to reduced-order modeling of constrained mechanical systems

TL;DR

This paper demonstrates a non-intrusive operator inference method to create reduced-order models of constrained mechanical systems governed by DAEs, enabling efficient simulations without access to system matrices.

Contribution

It introduces a novel application of operator inference to index 2 and 3 DAEs, ensuring stability and interpretability through semidefinite programming.

Findings

Effective reduced models obtained from solution snapshots

Models maintain stability and interpretability

Numerical tests show accurate system representation

Abstract

Constrained mechanical systems occur in many applications, such as modeling of robots and other multibody systems. In this case, the motion is governed by a system of differential-algebraic equations (DAE), often with large and sparse system matrices. The problem dimension strongly influences the effectiveness of simulations for system analysis, optimization, and control, given limited computational resources. Therefore, we aim to obtain a simplified surrogate model with a few degrees of freedom that is able to accurately represent the motion and other important properties of the original high-dimensional DAE model. Classical model reduction methods intrusively exploit the system matrices to construct the projection of the high-fidelity model onto a low-dimensional subspace. In practice, the dynamical equations are frequently an inaccessible part of proprietary software products. In this work, we show an application of the non-intrusive operator inference (OpInf) method to DAE systems of index 2 and 3. Considering the fact that for proper DAEs there exists an ODE realization on the so-called hidden manifold, the OpInf optimization problem directly provides the underlying ODE representation of the given DAE system in the reduced subspace. A significant advantage is that only the DAE solution snapshots in a compressed form are required for identification of the reduced system matrices. Stability and interpretability of the reduced-order model is guaranteed by enforcing the symmetric positive definite structure of the system operators using semidefinite programming. The numerical results demonstrate the implementation of the proposed methodology for different examples of constrained mechanical systems, tested for various loading conditions.