Search

The theory of spectral submanifolds (SSMs) has emerged as a powerful tool for constructing rigorous, low-dimensional reduced-order models (ROMs) of high-dimensional nonlinear mechanical systems. A direct computation of SSMs requires explicit knowledge of nonlinear coefficients in the equations of motion, which limits their applicability to generic finite-element (FE) solvers. Here, we propose a non-intrusive algorithm for the computation of the SSMs and the associated ROMs up to arbitrary polynomial orders. This non-intrusive algorithm only requires system nonlinearity as a black box and hence, enables SSM-based model reduction via generic finite-element software. Our expressions and algorithms are valid for systems with up to cubic-order nonlinearities, including velocity-dependent nonlinear terms, asymmetric damping, and stiffness matrices, and hence work for a large class of mechanics problems. We demonstrate the effectiveness of the proposed non-intrusive approach over a variety of FE examples of increasing complexity, including a micro-resonator FE model containing more than a million degrees of freedom.

Dynamical systems in engineering and physics are often subject to irregular excitations that are best modeled as random. Monte Carlo simulations are routinely performed on such random models to obtain statistics on their long-term response. Such simulations, however, are prohibitively expensive and time consuming for high-dimensional nonlinear systems. Here we propose to decrease this numerical burden significantly by reducing the full system to very low-dimensional, attracting, random invariant manifolds in its phase space and performing the Monte Carlo simulations on that reduced dynamical system. The random spectral submanifolds (SSMs) we construct for this purpose generalize the concept of SSMs from deterministic systems under uniformly bounded random forcing. We illustrate the accuracy and speed of random SSM reduction by computing the SSM-reduced power spectral density of the randomly forced mechanical systems that range from simple oscillator chains to finite-element models of beams and plates., Comment: 26 pages, 15 figures

For mechanical systems subject to periodic excitation, forced response curves (FRCs) depict the relationship between the amplitude of the periodic response and the forcing frequency. For nonlinear systems, this functional relationship is different for different forcing amplitudes. Forced response surfaces (FRSs), which relate the response amplitude to both forcing frequency and forcing amplitude, are then required in such settings. Yet, FRSs have been rarely computed in the literature due to the higher numerical effort they require. Here, we use spectral submanifolds (SSMs) to construct reduced-order models (ROMs) for high-dimensional mechanical systems and then use multidimensional manifold continuation of fixed points of the SSM-based ROMs to efficiently extract the FRSs. Ridges and trenches in an FRS characterize the main features of the forced response. We show how to extract these ridges and trenches directly without computing the FRS via reduced optimization problems on the ROMs. We demonstrate the effectiveness and efficiency of the proposed approach by calculating the FRSs and their ridges and trenches for a plate with a 1:1 internal resonance and for a shallow shell with a 1:2 internal resonance.

We use the recent theory of Spectral Submanifolds (SSM) for model reduction of nonlinear mechanical systems subject to parametric excitations. Specifically, we develop expressions for higher-order nonautonomous terms in the parameterization of SSMs and their reduced dynamics. We provide these results both for general first-order as well as second-order mechanical systems under periodic and quasiperiodic excitation using a multi-index based approach, thereby optimizing memory requirements and the computational procedure. We further provide theoretical results that simplify the SSM parametrization for general second-order dynamical systems. More practically, we show how the reduced dynamics on the SSM can be used to extract the resonance tongues and the forced response around the principal resonances in parametrically excited systems. In the case of two-dimensional SSMs, we formulate explicit expressions for computing the steady-state response as the zero-level set of a two-dimensional function for systems that are subject to external as well as parametric excitation. This allows us to parallelize the computation of the forced response over the range of excitation frequencies. We demonstrate our results on several examples of varying complexity, including finite-element type examples of mechanical systems. Furthermore, we provide an open-source implementation of all these results in the software package SSMTool.

Very high dimensional nonlinear systems arise in many engineering problems due to semi-discretization of the governing partial differential equations, e.g. through finite element methods. The complexity of these systems present computational challenges for direct application to automatic control. While model reduction has seen ubiquitous applications in control, the use of nonlinear model reduction methods in this setting remains difficult. The problem lies in preserving the structure of the nonlinear dynamics in the reduced order model for high-fidelity control. In this work, we leverage recent advances in Spectral Submanifold (SSM) theory to enable model reduction under well-defined assumptions for the purpose of efficiently synthesizing feedback controllers., Comment: 8 pages, 6 figures, conference on decision and control 2022

Dynamical systems are often subject to algebraic constraints in conjunction to their governing ordinary differential equations. In particular, multibody systems are commonly subject to configuration constraints that define kinematic compatibility between the motion of different bodies. A full-scale numerical simulation of such constrained problems is challenging, making reduced-order models (ROMs) of paramount importance. In this work, we show how to use spectral submanifolds (SSMs) to construct rigorous ROMs for mechanical systems with configuration constraints. These SSM-based ROMs enable the direct extraction of backbone curves and forced response curves, and facilitate efficient bifurcation analysis. We demonstrate the effectiveness of this SSM-based reduction procedure on several examples of varying complexity, including nonlinear finite-element models of multi-body systems. We also provide an open-source implementation of the proposed method that also contains all details of our numerical examples.

We show how spectral submanifold theory can be used to construct reduced-order models for harmonically excited mechanical systems with internal resonances. Efficient calculations of periodic and quasi-periodic responses with the reduced-order models are discussed in this paper and its companion, Part II, respectively. The dimension of a reduced-order model is determined by the number of modes involved in the internal resonance, independently of the dimension of the full system. The periodic responses of the full system are obtained as equilibria of the reduced-order model on spectral submanifolds. The forced response curve of periodic orbits then becomes a manifold of equilibria, which can be easily extracted using parameter continuation. To demonstrate the effectiveness and efficiency of the reduction, we compute the forced response curves of several high-dimensional nonlinear mechanical systems, including the finite-element models of a von K\'arm\'an beam and a plate.

Invariant manifolds are important constructs for the quantitative and qualitative understanding of nonlinear phenomena in dynamical systems. In nonlinear damped mechanical systems, for instance, spectral submanifolds have emerged as useful tools for the computation of forced response curves, backbone curves, detached resonance curves (isolas) via exact reduced-order models. For conservative nonlinear mechanical systems, Lyapunov subcenter manifolds and their reduced dynamics provide a way to identify nonlinear amplitude-frequency relationships in the form of conservative backbone curves. Despite these powerful predictions offered by invariant manifolds, their use has largely been limited to low-dimensional academic examples. This is because several challenges render their computation unfeasible for realistic engineering structures described by finite-element models. In this work, we address these computational challenges and develop methods for computing invariant manifolds and their reduced dynamics in very high-dimensional nonlinear systems arising from spatial discretization of the governing partial differential equations. We illustrate our computational algorithms on finite-element models of mechanical structures that range from a simple beam containing tens of degrees of freedom to an aircraft wing containing more than a hundred-thousand degrees of freedom.

With the growing amount of text in health data, there have been rapid advances in large pre-trained models that can be applied to a wide variety of biomedical tasks with minimal task-specific modifications. Emphasizing the cost of these models, which renders technical replication challenging, this paper summarizes experiments conducted in replicating BioBERT and further pre-training and careful fine-tuning in the biomedical domain. We also investigate the effectiveness of domain-specific and domain-agnostic pre-trained models across downstream biomedical NLP tasks. Our finding confirms that pre-trained models can be impactful in some downstream NLP tasks (QA and NER) in the biomedical domain; however, this improvement may not justify the high cost of domain-specific pre-training.

We propose a reformulation for the integral equations approach of Jain, Breunung \& Haller [Nonlinear Dyn. 97, 313--341 (2019)] to steady-state response computation for periodically forced nonlinear mechanical systems. This reformulation results in additional speed-up and better convergence. We show that the solutions of the reformulated equations are in one-to-one correspondence with those of the original integral equations and derive conditions under which a collocation type approximation converges to the exact solution in the reformulated setting. Furthermore, we observe that model reduction using a selected set of vibration modes of the linearized system substantially enhances the computational performance. Finally, we discuss an open-source implementation of this approach and demonstrate the gains in computational performance using three examples that also include nonlinear finite-element models.

Model reduction of large nonlinear systems often involves the projection of the governing equations onto linear subspaces spanned by carefully-selected modes. The criteria to select the modes relevant for reduction are usually problem-specific and heuristic. In this work, we propose a rigorous mode-selection criterion based on the recent theory of Spectral Submanifolds (SSM), which facilitates a reliable projection of the governing nonlinear equations onto modal subspaces. SSMs are exact invariant manifolds in the phase space that act as nonlinear continuations of linear normal modes. Our criterion identifies critical linear normal modes whose associated SSMs have locally the largest curvature. These modes should then be included in any projection-based model reduction as they are the most sensitive to nonlinearities. To make this mode selection automatic, we develop explicit formulas for the scalar curvature of an SSM and provide an open-source numerical implementation of our mode-selection procedure. We illustrate the power of this procedure by accurately reproducing the forced-response curves on three examples of varying complexity, including high-dimensional finite element models.

Distributed word embeddings have yielded state-of-the-art performance in many NLP tasks, mainly due to their success in capturing useful semantic information. These representations assign only a single vector to each word whereas a large number of words are polysemous (i.e., have multiple meanings). In this work, we approach this critical problem in lexical semantics, namely that of representing various senses of polysemous words in vector spaces. We propose a topic modeling based skip-gram approach for learning multi-prototype word embeddings. We also introduce a method to prune the embeddings determined by the probabilistic representation of the word in each topic. We use our embeddings to show that they can capture the context and word similarity strongly and outperform various state-of-the-art implementations., Comment: Accepted at ACL supported conference for Natural Language & Speech Processing. https://www.aclweb.org/anthology/W19-74, Year: 2019

The thermal dynamics in thermo-mechanical systems exhibits a much slower time scale compared to the structural dynamics. In this work, we use the method of multiple scales to reduce the thermo-mechanical structural models with a slowly-varying temperature distribution in a systematic manner. In the process, we construct a reduction basis that adapts according to the instantaneous temperature distribution of the structure, facilitating an efficient reduction in the number of unknown. As a proof of concept, we demonstrate the method on a range of linear and nonlinear beam examples and obtain a consistently better accuracy and reduction in the number of unknowns than standard the Galerkin projection using a constant basis.

We discuss an integral equation approach that enables fast computation of the response of nonlinear multi-degree-of-freedom mechanical systems under periodic and quasi-periodic external excitation. The kernel of this integral equation is a Green's function that we compute explicitly for general mechanical systems. We derive conditions under which the integral equation can be solved by a simple and fast Picard iteration even for non-smooth mechanical systems. The convergence of this iteration cannot be guaranteed for near-resonant forcing, for which we employ a Newton--Raphson iteration instead, obtaining robust convergence. We further show that this integral-equation approach can be appended with standard continuation schemes to achieve an additional, significant performance increase over common approaches to computing steady-state response.

Common trends in model order reduction of large nonlinear finite-element-discretized systems involve the introduction of a linear mapping into a reduced set of unknowns, followed by Galerkin projection of the governing equations onto a constant reduction basis. Though this reduces the number of unknowns in the system, the computational cost for obtaining the solution could still be high due to the prohibitive computational costs involved in the evaluation of nonlinear terms. Hyper-reduction methods are then seen as a fast way of approximating the nonlinearity in the system of equations. In the finite element context, the energy conserving sampling and weighing (ECSW) method has emerged as a stability and structure-preserving method for hyper-reduction. Classical hyper-reduction techniques, however, are applicable only in the context of linear mappings into the reduction subspace. In this work, we extend the concept of hyper-reduction using ECSW to general nonlinear mappings, while retaining its desirable stability and structure-preserving properties. As a proof of concept, the proposed hyper-reduction technique is demonstrated over models of a flat plate and a realistic wing structure, whose dynamics has been shown to evolve over a nonlinear (quadratic) manifold. An online speed-up of over one thousand times relative to the full system has been obtained for the wing structure using the proposed method, which is higher than its linear counterpart using the ECSW.

We apply two recently formulated mathematical techniques, Slow-Fast Decomposition (SFD) and Spectral Submanifold (SSM) reduction, to a von Karman beam with geometric nonlinearities and viscoelastic damping. SFD identifies a global slow manifold in the full system which attracts solutions at rates faster than typical rates within the manifold. An SSM, the smoothest nonlinear continuation of a linear modal subspace, is then used to further reduce the beam equations within the slow manifold. This two-stage, mathematically exact procedure results in a drastic reduction of the finite-element beam model to a one-degree-of freedom nonlinear oscillator. We also introduce the technique of spectral quotient analysis, which gives the number of modes relevant for reduction as output rather than input to the reduction process.

This paper describes the use of a quadratic manifold for the model order reduction of structural dynamics problems featuring geometric nonlinearities. The manifold is tangent to a subspace spanned by the most relevant vibration modes, and its curvature is provided by modal derivatives obtained by sensitivity analysis of the eigenvalue problem, or its static approximation, along the vibration modes. The construction of the quadratic manifold requires minimal computational effort once the vibration modes are known. The reduced order model is then obtained by Galerkin projection, where the configuration-dependent tangent space of the manifold is used to project the discretized equations of motion.