Application of Operator Inference to Reduced‐Order Modeling of Constrained Mechanical Systems

Source abstract

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. Data‐driven identification of the reduced second‐order ODE representations of proper DAE systems, applied for mechanical models with position and velocity constraints; Depending on the constraint type, use displacement and velocity snapshot matrices to eliminate the algebraic conditions; Selection of the number of snapshots based on the impulse response of the system.