Weighted model reduction for nonlinear systems and computation method

The very general problem of model reduction of nonlinear systems was made tractable by focusing on the very large subclass consisting of linear subsystems connected by nonlinear interfaces. Such problems constitute a large part of the nonlinear structural problems encountered in addressing the Sandia missions. A synthesis approach to this class of problems was developed consisting of: detailed modeling of the interface mechanics; collapsing the interface simulation results into simple nonlinear interface models; constructing system models by assembling model approximations of the linear subsystems and the nonlinear interface models. These system models, though nonlinear, would have very few degrees of freedom. A paradigm problem, that of machine tool vibration, was selected for application of the reduction approach outlined above. Research results achieved along the way as well as the overall modeling of a specific machine tool have been very encouraging. In order to confirm the interface models resulting from simulation, it was necessary to develop techniques to deduce interface mechanics from experimental data collected from the overall nonlinear structure. A program to develop such techniques was also pursued with good success.

Finite element models (FEMs) are extensively used in the design optimization of utility scale steam turbines. As an example, by simulating multiple startup scenarios of steam power plants, engineers can obtain turbine designs that minimize material utilization, and at the same time, avoid the damaging effects of large thermal stresses or rubs between rotating and stationary parts. Unfortunately, FEMs are computationally expensive and only a limited amount of simulations can be afforded to get the final design. For this reason, numerous model reduction techniques have been developed to reduce the size of the original model without a significant loss of accuracy. When the models are nonlinear, as is the case for steam turbine FEMs, model reduction techniques are relatively scarce and their effectiveness becomes application dependent. Although there is an abundant literature on model reduction for nonlinear systems, many of these techniques become impractical when applied to a realistic industrial problem. This paper focuses on a class of nonlinear FEM characteristic of thermo-elastic problems with large temperature excursions. A brief overview of popular model reduction techniques is presented along with a detailed description of the computational challenges faced when applying them to a realistic problem. The main contribution of this work is a set of modifications to existing methods to increase their computational efficiency. The methodology is demonstrated on a steam turbine model, achieving a model size reduction by four orders of magnitude with only 4% loss of accuracy with respect to the full order FEMs.

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

This paper provides three model reduction methods for linear time-invariant systems in the view of the Riemannian Newton method and the Jacobi-Davidson method. First, the computation of Hankel singular values is converted into the linear eigenproblem by the similarity transformation. The Riemannian Newton method is used to establish the model reduction method. Besides, we introduce the Jacobi-Davidson method with the block version for the linear eigenproblem and present the corresponding model reduction method, which can be seen as an acceleration of the former method. Both the resulting reduced systems can be equivalent to the reduced system originating from a balancing transformation. Then, the computation of Hankel singular values is transformed into the generalized eigenproblem. The Jacobi-Davidson method is employed to establish the model reduction method, which can also lead to the reduced system equivalent to that resulting from a balancing transformation. This method can also be regarded as an acceleration of a Riemannian Newton method. Moreover, the application for model reduction of nonlinear systems with inhomogeneous conditions is also investigated.

This paper extends some methods in model order reduction (MOR) of nonlinear models by trajectory piecewise-linear models based on output weighting (TPWLOW). The devised methods generalize the balanced-truncation method as the core for the proposed algorithms. These algorithms preserve the stability in some mentioned conditions. An error bound of reduction is calculated for these algorithms. In addition, a Krylov–TBR combinational MOR is also casted which benefits from the advantages of both balanced-truncation and Krylov subspace methods. The results show that in balanced-truncation-based algorithms, the efficiency of the first proposed algorithm fairly exceeds the second one. In addition, it is shown that the Krylov–TBR algorithm can reach the same efficiency despite remarkable lower order, provided that the designer can find a stable subspace. The robustness of TPWLOW-reduced models to linearization point movement is also evaluated, which shows the less sensitivity of the proposed methods in comparison with formerly introduced trajectory piecewise-linear models.

The development of structure‐preserving approximation methods, which regard fundamental underlying physical principles, is an active field of research. Especially when the application of model reduction is desirable, systematic and rather generic approaches are of great interest. In this contribution we discuss a structure‐preserving Galerkin approach for a prototypical class of nonlinear partial differential equations on networks. Its derivation is guided by port‐Hamiltonian‐type modeling and appropriate variational principles. Also complexity‐reduction schemes can be integrated in a structure‐preserving way, which becomes crucial in the context of model reduction for nonlinear systems.

Model Reduction for Nonlinear Systems with Distributed Parameters

This article is concerned with the problem of the Hankel-norm model reduction for stochastic discrete-time nonlinear systems in interval type-2 (IT2) Takagi-Sugeno (T-S) fuzzy framework. The IT2 T-S fuzzy model is an efficient model for describing uncertain nonlinear systems, and the model reduction is to simplify the high-order complex systems by reducing the order of the original system. The aim of this article is to reduce the order of the original stochastic discrete-time IT2 fuzzy system into lower order system without ignoring the influence of IT2 membership functions. First, the Hankel-norm performance of the stochastic discrete-time IT2 fuzzy model is analyzed. Then, based on the projection theorem and cone complementary linearization approach, a convex Hankel-norm-based model reduction approach subject to conditions in the form of linear matrix inequalities (LMIs) is obtained. A membership-functions-dependent (MFD) technique is applied to capture the information of IT2 membership functions and further reduce the conservativeness. A numerical example is presented to illustrate the effectiveness of the proposed results.

Interpolatory projection methods for model reduction of nonparametric linear dynamical systems have been successfully extended to nonparametric bilinear dynamical systems. However, this has not yet occurred for parametric bilinear systems. In this work, we aim to close this gap by providing a natural extension of interpolatory projections to model reduction of parametric bilinear dynamical systems. We introduce necessary conditions that the projection subspaces must satisfy to obtain parametric tangential interpolation of each subsystem transfer function. These conditions also guarantee that the parameter sensitivities (Jacobian) of each subsystem transfer function are matched tangentially by those of the corresponding reduced-order model transfer function. Similarly, we obtain conditions for interpolating the parameter Hessian of the transfer function by including additional vectors in the projection subspaces. As in the parametric linear case, the basis construction for two-sided projections does not require computing the Jacobian or the Hessian.

The problem of model reduction of linear systems with certain interconnection structure is considered in this paper. To preserve the interconnection structure between subsystems in the reduction, special care needs to be taken. This problem is important and timely because of the recent focus on complex networked systems in control engineering. Two different model reduction methods are introduced and compared in this paper. Both methods are extensions to the well‐known balanced truncation method. Compared with earlier work in the area these methods use a more general linear fractional transformation framework, and utilize linear matrix inequalities. Furthermore, new approximation error bounds that reduce to classical bounds in special cases are derived. The so‐called structured Hankel singular values are used in the methods, and indicate how important states in the subsystems are with respect to a chosen input–output map for the entire interconnected system. It is shown how these structured Hankel singular values can be used to select an approximation order. Finally, the two methods are applied to a model of a mechanical device. Copyright © 2008 John Wiley & Sons, Ltd.

Time-domain computation of large amplitude 3D ship motions with forward speed

In this paper a new balanced realization method for nonlinear systems is proposed which is based on singular value analysis of Hankel operators. The proposed method balances the relationship between the input-to-state behavior and the state-to-output behavior of nonlinear dynamical systems, whereas the existing results only balance the relationship among the coordinate axes of the state-space. This result is expected to be a basis for new model reduction and system identification of nonlinear systems.

<p style='text-indent:20px;'>We consider model order reduction of a nonlinear cable-mass system modeled by a 1D wave equation with interior damping and dynamic boundary conditions. The system is driven by a time dependent forcing input to a linear mass-spring system at one boundary. The goal of the model reduction is to produce a low order model that produces an accurate approximation to the displacement and velocity of the mass in the nonlinear mass-spring system at the opposite boundary. We first prove that the linearized and nonlinear unforced systems are well-posed and exponentially stable under certain conditions on the damping parameters, and then consider a balanced truncation method to generate the reduced order model (ROM) of the nonlinear input-output system. Little is known about model reduction of nonlinear input-output systems, and so we present detailed numerical experiments concerning the performance of the nonlinear ROM. We find that the ROM is accurate for many different combinations of model parameters.

The settlement calculation model for soft clay foundation is established based on Hooke’s law and the Duncan-Chang model. By introducing the concept of damage ratio, the method for determining the deformation modulus of soil before and after damage under load is presented, and a nonlinear settlement calculation method which considers the structural characteristics of soft clay is proposed. In the end, a practical engineering example is analyzed with the proposed method and some current methods for soft clay settlement calculation. The results indicated that the proposed method is feasible and applicable in practical engineering.

Model order reduction of a dynamical linear time-invariant system appears in many applications from science and engineering. Numerically reliable SVD-based methods for this task require in general $\mathcal{O}(n^3)$ floating-point arithmetic operations, with n being in the range 103−105 for many practical applications. In this paper we investigate the use of graphics processors (GPUs) to accelerate model reduction of large-scale linear systems by off-loading the computationally intensive tasks to this device. Experiments on a hybrid platform consisting of state-of-the-art general-purpose multi-core processors and a GPU illustrate the potential of this approach.