Dimensionality reduction for data sets representing high-dimensional nonlinear mechanical systems is of crucial importance in science and technology. Low-dimensional models are expected to reduce computational cost and capture the essential physics of a high-dimensional system from data. Given the growing interest, for example, in light mechanical structures and MEMS devices, there is the need for truly nonlinear models, capturing amplitude-dependent properties and com-peting steady states solutions, which are increasingly important to identify, as highlighted in experiments of nonlinear mechanical vibrations [1, 2, 3, 4, 5, 6, 7, 8]. Predicting coexisting stable and unstable forced responses for a broad range of forcing amplitudes and frequencies is paramount in structural dynamics. However, a generally applicable technique returning such reliable low-dimensional models of nonlinear mechanical vibrations has not emerged yet.
The most common approaches to data-driven reduced-order modeling are the Proper Orthogonal Decomposition (POD) followed by a Galerkin projection [9, 10, 11, 12] and the Dynamic Mode Decomposition (DMD) [13, 14]. The former approach requires the knowledge of the governing equations of motion and, once a relevant number of modes is identified from data, projects these equations onto those modes to construct a reduced-order model. DMD and its improved versions [15, 16, 17, 18, 19], supported by Koopman operator theory [20, 21], seek a low-rank approximation to the dynamics of observable data without reliance on the governing equations of motion. With this approach, DMD and Koopman mode expansions are able to linearize the observed dynamics around attracting fixed points on domains that cannot include additional fixed points or limit cycles [22, 23, 24, 25]. Therefore, while truly powerful for globally linearizable dynamics , these linear techniques cannot capture essentially nonlinear dynamical systems (or non-linearizable systems) with multiple coexisting steady states.
Other approaches treat the dimensionality reduction and the data-driven dynamical modeling as separate problems. Typically, the data is first processed via a dimensionality reduction algorithm, which ranges from POD or Principal Component Analysis (PCA), its kernelized version , manifold learning techniques [29, 30, 31]
or autoencoders[32, 33]. Afterwards, the dynamics are identified in the reduced coordinates using classic regression techniques (least-squares , LASSO , SINDy ), Bayesian learning techniques 
or neural networks in different architectures (fully connected, convolutional, recurrent)[37, 32, 38, 39, 40]. Some of these techniques return complex, black-box models (which may be non-physical ), while others offer sparse models (LASSO, SINDy, Bayesian learning), which allow for easy interpretation and analysis . The resulting dynamics, however, are intrinsically determined by the representation offered by dimensionality reduction algorithms, unless penalized in the optimization 
. Indeed, the advocated simplicity of sparse models depends critically on the reduction method, as even a linear coordinate change will dramatically destroy the sparsity of a model. In addition, those methods feature a high number of hyperparameters that need to be tuned extensively for good performance. Most importantly, the eventual lack of predictive capabilities often makes the models unattractive for practical use. Indeed, the insertion of parameter variations, disturbances or external forcing into these models is generally heuristic, and hence returns questionable conclusions.
Our objective here is to discuss a new data-driven reduced-order modeling approach in the context of mechanical vibrations. Based on the recent theory of spectral submanifolds (SSMs) 
, this approach identifies very low dimensional, sparse models over different time scales by restricting the full system dynamics to a nested family of attractors. The SSMs forming this family are the smoothest nonlinear continuations of the eigenspaces of the linear part of the dynamical system. When transformed to a normal form, the reduced dynamics on each SSM is low-dimensional, sparse and relevant for all trajectories in the domain of attraction of the SSM. Importantly, each SSM may contain multiple coexisting steady states and hence capture non-linearizable dynamics.
and an open-sourceMatlab® implementation, SSMTool, for an arbitrary, finite-dimensional dynamical system is available in . Another concept closely linked to SSMs is that of invariant foliations , which provides a rigorous nonlinear extension of classic linear modal analysis. Our present discussion of data-driven SSM-based models follows the terminology and notation of the more technical exposition in . The open-source Matlab® code, SSMLearn, used in our examples is available in .
The remainder of this paper is organized as follows. Section 2 introduces SSMs and discusses their relevance for data-driven model reduction, also depending on the type of experiments that generate the data. We also discuss how our method is complementary to (non-parametric) signal processing techniques in nonlinear system identification , ranging from the Hilbert transform and its variants [56, 57] to wavelet decompositions . In section 2, we summarize our data-driven identification of SSMs and the resulting explicit models on SSMs. We demonstrate the method in section 3 on a set of examples, which are all analyzed via the Matlab® implementation of our approach SSMLearn . In particular, after a preliminary numerical example, we examine two experimental datasets of nonlinear mechanical systems, one of which regards an internally resonant structure. The data sets in these examples come from diverse sources, from non-contact measurement systems (e.g., digital image correlation, laser scanner vibrometry) to classic accelerometers.
2 Spectral submanifolds and data-driven models on them
In this paper, we consider
-degree-of-freedom mechanical systems of the form
is a generalized coordinate vector,is a positive definite, symmetric mass matrix. The forcing vector contains all conservative and non-conservative forces, including linear and nonlinear ones. The matrix , its inverse and are of class with (smooth functions) or (analytic functions).
The equivalent first-order form of Eq. (1), with and , reads
where . We assume that is an asymptotically stable equilibrium and that is a semi-simple matrix and has
complex conjugate pairs of eigenvalues with negative real parts. We order these eigenvalueswith decreasing real parts, and we denote by the corresponding two-dimensional eigenspaces (or modal subspaces).
We denote by the direct sum of of these modal subspaces, i.e., . The -dimensional, spectral subspace is invariant for the linearization of system (2). Its reduced dynamics is governed by the eigenvalues , which, along with the conjugate ones, form the set . The spectral submanifold (SSM), , is the smoothest nonlinear continuation of the linear subspace , as can be deduced from the more abstract invariant manifold results of [59, 60, 61, 62]. Specifically, is the unique -dimensional, class invariant manifold of system (2) tangent to the spectral subspace at the origin. The existence of is guaranteed whenever the eigenvalues with are not in resonance with those in , i.e., for ,
as discussed in . From a numerical perspective, the nonresonance condition in Eq. (3) is violated if the absolute value of the left-hand side of the inequality is below a certain tolerance. In that case, one needs to add the resonant modal subspace to , resulting in the SSM of the form . This larger SSM can be used to capture nonlinear modal interactions, e.g., in weakly damped systems with rationally dependent linearized frequencies. The dynamics restricted to SSMs gives exact nonlinear reduced-order models for system (2) [42, 46, 50].
The most important SSMs from a data-driven perspective are slow SSMs, which are constructed over the spectral subspace spanned by the slowest modes: . Slow SSMs are attracting normally hyperbolic invariant manifolds  to which nearby trajectories converge exponentially fast , as illustrated in Fig. 1(a). Therefore, generic experiments on mechanical systems in the form of (2) will yield trajectories converging exponentially fast to slow SSMs, which in turn capture the asymptotic dynamics near the equilibrium. Faster time-scales of the dynamics can be extracted from trajectory data by model reduction to higher dimensional members of the nested slow SSM family . As an illustration, trajectories on slow SSMs are shown in Fig. 1
(b), along with their spectrogram (or short time Fourier transform) in Fig.1(c). The trajectory on the top is in , while the middle one is initialized close to and hence converges to , as seen from the disappearance of higher frequencies. Finally, the trajectory at the bottom in Fig. 1(b,c) belongs to , where two modal contributions can be clearly identified.
For trajectories with generic initial conditions—as those generated by hammer impacts—discarding the initial part of the measured signal yields trajectory data close to a slow -dimensional SSM, where is the number of dominant frequencies in the signal. This number is inferred from a preliminary time-frequency analysis (e.g. spectrograms, wavelet transforms ), as those of Fig. 1(c). The larger , the more data is needed to properly explore the SSM. In contrast, targeted experiments can focus on specific SSMs, which is the case for two-dimensional SSMs in resonance decay experiments [1, 4]. In that setting, a near-resonant oscillation is first isolated using a shaker, which is then turned off. This generates a system trajectory that decays towards the equilibrium along the targeted two-dimensional SSM, provided that no internal resonance occurs. Using the shaker, we can typically excite higher amplitudes in comparison to hammer impact testing, where the energy is spread among multiple modes.
2.1 Learning SSMs from data
To learn SSMs from data, we use the methodology presented in , which is implemented in the open-source Matlab® package, SSMLearn . In what follows, we sketch the main ideas of this method before going into the details of the data-driven reduced-order models that SSMLearn can identify.
Measuring all phase space variables of a mechanical system is generally unrealistic. Typically, only a limited set of observed quantities is available, so that we need to embed the SSM, , into a lower dimensional space of observables. According to the prevalence version of Whitney’s embedding theorem , almost all sets of independent and simultaneous measurements form an embedding space for -dimensional SSMs if . This is the case, for example, when displacements and velocities of at least material points of a mechanical system are available. Practical experiments, however, generally only provide the displacement, velocity or acceleration of a single material point, denoted as , recorded at time intervals. To this end, we exploit Takens delay embedding theorem , which, in its prevalence version , states that
forms an embedding space with probability one ifunder generic nondegeneracy conditions on the sampling time . Further spaces may also qualify in practice, e.g., featuring or constructed from multiple measurements augmented by delays, but one needs to examine on a case-by-case basis whether these are embedding spaces or not.
We denote by the embedded SSM, for which now we need to construct a reduced-order model in the embedding space. We assume that the equilibrium is at and that does not fold over its tangent space at the origin , so that we can construct a data-driven graph-style parametrization for over . We let be the matrix whose orthonormal columns span and we define the SSM parametrization, , as
where we assume that is a multivariate polynomial from order to . The matrix , as well as the coefficients of the polynomial
, can be found via constrained maximum likelihood estimation of (4), as discussed in .
Once trajectories in the projection coordinates are known, we can identify the SSM-reduced dynamics. Here, the idea is to find the extended normal form of the vector field governing the dynamics in the projection coordinate (or reduced) domain , motivated by classic studies of bifurcations [66, 67]. Specifically, we need to find an invertible change of coordinates (and its inverse) that brings the SSM-reduced dynamics to its simplest possible complex polynomial form with . The linear part of is the diagonal matrix of the eigenvalues related to the SSM, with denoting complex modal coordinates for the linearized system. The maps , and are multivariate polynomials with their coefficients determined from an extended normal form approach used in classic unfoldings of bifurcations [66, 67]. In this approach, the classic Poincaré  normal form construct is relaxed in that not only resonant but also near-resonant terms are kept in the normal form (see  for more details). This normalization renders a sparse vector field extracting the fundamental physics, as we discuss in the next section. We determine resonant coefficients from an initial estimate of the linearized dynamics, and we identify from data the maps , and by minimizing the conjugacy error, as explained in detail in . For example, the structure of the cubic normal form for a two-dimensional SSM is
which resembles the classic Hopf normal form . These normal form models are particularly simple to handle in polar coordinates , defined as for .
2.2 Interpretability and extrapolation from SSM-reduced models
The most general normal form on a -dimensional SSM is
Some explicit examples are presented in the examples of section 3, including cubic polar normal forms of two-dimensional and four dimensional SSMs, the latter appearing both for non-resonant eigenvalues and for a resonance. If the linearized frequencies are non-resonant, then and only depend on the amplitudes . The normal form (6) then decouples the amplitude dynamics from the phase dynamics. This enables us to distinguish different modal contributions, perform a slow-fast decomposition, detect modal interactions and analyze the uncoupled oscillator limit. The zero-amplitude limit of the functions and converges to the linearized damping and frequency of mode , i.e.,
Hence, and are the nonlinear continuations of these linear quantities, characterizing how dissipation and frequency change with respect to the amplitudes (and phases for internally resonant systems). For a two-dimensional SSM, the parametrized curves and are the backbones of transient oscillations [1, 44, 47, 8], representing the instantaneous damping and frequency as nonlinear functions of the normal form amplitude . Normal form amplitudes do not, however, have any direct physical meaning. For physical insights, we need to express any amplitude of interest via the SSM parametrization and the normal form transformation . For instance, for two-dimensional SSMs and for a scalar quantity defined on the observable space , the amplitude of the oscillations can be defined as [44, 46]
Then, backbone curves can be expressed as parametric curves and .
SSMs are robust features of the dynamics, because they survive under small autonomous perturbations and even under some non-autonomous perturbations of the vector field (2) . The most important class of these perturbations in our context is that of small external time-periodic forcing appearing on the right-hand side of Eq. (1). In that case, the autonomous SSM will serve as the leading order approximation for a non-autonomous, time-periodic SSM that carries reduced, time-periodic dynamics [42, 47, 48, 50]. With the addition of such forcing, the normal form (6) becomes 
where is the forcing frequency and the forcing amplitudes for each mode. Generally, numerical continuation is necessary for studying periodic responses and eventual bifurcations of (9) depending on forcing frequencies and amplitudes. In the simplest case of , however, we can introduce the phase shift to obtain from (9) the forced normal form
which yields closed form predictions for amplitudes and phases of the forced periodic solutions:
known as frequency response curves (FRCs), parametrized by the amplitude . Predictions of these curves from unforced data, however, have generally been unavailable. Physical amplitudes can be derived from the predictions of (11) using Eq. (8) and the stability of the predicted forced response can be derived from the Jacobian of the vector field (10) . We note find from Eq. (11) that the forced backbone curve (the location of maximal amplitude responses of FRCs under varying ) coincides with that of decaying oscillations, given by . Specifically, maximal amplitude responses occur at amplitudes satisfying , and phase-lag quadrature, i.e., . These maximal amplitude responses can be used to calibrate the normal form forcing amplitude with experimentally exerted forcing levels.
Equations (9,10) have accuracy [47, 53], but higher-order approximations can improve this accuracy further . We expect, for example, that forced backbone curves depart from those of decaying oscillations at large motion and/or large forcing amplitude values . From a data-driven perspective, once the autonomous core of Eqs. (9,10) is identified, we only need to calibrate the forcing amplitudes for predicting FRCs. In matching experimental results, one calibration point is sufficient if the forcing amplitude is kept constant during experimental frequency sweeps. The SSM-parametrization under forcing is also modulated by a small time-periodic component [49, 53]. This may be disregarded since the leading-order non-linearizable dynamics is already captured by forcing in the reduced dynamics (9,10).
We now discuss some examples that illustrate the power of the SSM-based, data-driven model reduction method we have discussed. Our first example is a chain of lumped oscillators, while the other two involve data from laboratory experiments. Additional details and further examples can also be found in the Matlab® live-scripts in .
To express trajectory reconstruction errors, we use the normalized mean trajectory error NMTE, which, for a dataset of instances of observable points and their reconstruction , is defined as
Here is a relevant normalization vector, which is usually taken to be the data point with the maximum norm in the dataset. To validate the reduced dynamics on a test trajectory, we integrate the reduced-order model from the same initial condition and compare the results.
3.1 Identification of SSMs in a chain of oscillators
We consider the chain of oscillators sketched in Fig. 2(a), where we set the first mass as [kg] and the others as [kg]. We also assume all spring-dampers to be linear with unitary stiffness, except for the leftmost one that exerts a nonlinear force on the first mass. The linear damping matrix for the system is proportional to the mass and stiffness matrices with constants and , with the resulting eigenvalues at the trivial equilibrium reported in Fig. 2(a).
We start with the study of the slow two-dimensional SSM of the oscillator chain. We compute this SSM via SSMTool , from which we initialize the two decaying trajectories shown in Fig. 2(b); one of these trajectories is used for testing the constructed model. We identify reduced-order models from two different observables. The first observable set is the set of all phase space variables, while the second is a set of delayed samples of the (scalar) displacement of the rightmost mass . We select the delay embedding of minimal dimension (five) required by the Takens theorem. The embedding coordinates are, therefore, , where the sampling time is [s]. A cubic-order parametrization for the phase space embedding and a parametrization for the delay embedding show good accuracy. The SSM and its embedding in the delay space are shown in Fig. 2(b,c). The flat appearance of the manifold in 2(d) in the delay space is a general phenomenon, as shown mathematically in . The cubic polar normal form on the phase-space-embedded SSM is found by SSMLearn to be
A similar model is identified for the delay embedding. Both reduced-order models capture well the dynamics of the testing trajectories, as seen in Fig. 2(e-h), with less than NMTE error. The instantaneous damping and frequency are shown in 3(b,c), displaying only a minimal disagreement. We note that this identification is robust against perturbations of the initial condition. Indeed, if we initialize trajectories slightly off the SSM as shown in Fig. 3(a), then SSMLearn still finds a good approximation for the reduced dynamics, as demonstrated by the curves in Fig. 3(c). If these perturbations are not small enough for the dynamics to be described by a two-dimensional SSMs, then we need to increase the SSM dimension. For instance, we computed trajectories, shown in Fig. 4(a), decaying along the slow four-dimensional SSM, . Five of these trajectories are used for training and one is left for testing our reduced-order model.
The cubic normalized, SSM-reduced dynamics identified by SSMLearn has a NMTE error and is of the form
Eventual differences between the dynamics of the slowest mode in (13) with respect to those in (14) are due to different amplitude scalings. Prediction of a test trajectory based on the model (14) is shown in Fig. 4(b). The instantaneous frequencies for the two modes are shown in Fig. 4(c,d). These are surfaces since both frequencies depend (either weakly or strongly) on both modal amplitudes.
For additional validation, we show FRCs of the models (13)-(14) for different forcing amplitudes in Fig. 4(e,f) around the first two eigenfrequencies. While on the two-dimensional SSM we have the closed-form solution (11), FRCs on the four-dimensional SSM are computed using the periodic orbit toolbox of the numerical continuation core coco . These plots in Fig. 4(e,f) are completed with backbones curves and forced responses obtained via numerical integration of the full model. The forcing only acts in the direction of the first two modes, with amplitudes and [mN]. Our data-driven predictions, which are based only on unforced data and a simple calibration procedure for the normal form forcing amplitudes, are in close agreement with the responses from the full system.
3.2 Resonance decay in the Brake-Reuss beam
The Brake-Reuss beam (BRB) is a benchmark system in the study of jointed structures [72, 7, 8]. In our study, it consists of two 304 stainless steel beams assembled with a lap joint, as shown in Fig. 5(a). While full models for these structures may not be smooth, we find that trajectory data can be fitted well to smooth models, thereby justifying an SSM-based approach. The data considered here arises from a single resonance decay test, available from , targeting the slowest structural mode. One observable is the measurement from an accelerometer mounted, as shown in Fig. 5(a), on the shaker with time history illustrated in Fig. 5(c). Another observable is the displacement field of the bottom side of the beam, measured using digital image correlation (DIC). The latter dataset, consisting of points over [cm] of beam length, has a limited time range due to limitations in camera memory. The initial evolution of the measured displacement field is depicted in Fig. 5(b).
Our goal in this example is to construct a nonlinear reduced-order model using displacement data and validate it on the acceleration measurement. We truncate the time signals after shaker release to eliminate the influence of disturbances from non-perfect detachment. Nevertheless, high frequency contributions decay rapidly and the transient settles along the slowest SSM. To diversify the data, we augment the displacement with four delayed measurements, so that the observable phase space has dimension . The SSM is approximately a plane in this space, but the reduced dynamics is highly nonlinear. For adequate accuracy, the normal form indeed needs terms up to to capture the dynamics:
The model can be used to approximate the beam kinetic energy as
where is the number of DIC measurement locations and [kg] is the beam mass. As discussed in , the kinetic energy amplitude is a good proxy for the instantaneous decay properties, i.e., the instantaneous damping ratio and frequency, shown in Fig. 5(d,e), respectively. The instantaneous damping ratio is defined from the normal form dynamics (15) as , expressed in percentage. The damping exhibits a strong variation from its linear limit, while the frequency here shows a peculiar softening-hardening trend.
We validate our displacement-based SSM-reduced model on the data from the accelerometer located at [mm] from the left end of the beam, as shown in Fig. 5
(a). This signal is reconstructed from the model by interpolating in the grid to obtain the accelerometer location and differentiating in time. The NMTE error is onlyand the model practically denoises the signal, as seen in Fig. 5(f,g). A further validation in Fig. 5(h,i) compares the instantaneous decay properties of the data-driven model (15) to those extracted using the Peak Finding and Fitting [57, 8] signal processing technique. There is close agreement among these curves, especially in the strongly nonlinear oscillation regime.
3.3 Impacts on an internally resonant tester structure
Our final example is the resonant tester shown in Fig. 6(a,b). It consists of two beam-like parts made of aluminum 6061-T6, where the external beam is C-shaped and clamped to the ground on one side, while the internal beam is jointed to the external one via bolts three bolts: two side bolts are torqued to [Nm] for structural integrity whereas the middle bolt is torqued to [Nm] for enhanced frictional slip and associated nonlinearities. Additionally, a linear spring (Model 1NCH2, Grainger, Inc.) connects the tip of the external beam to a fixed rigid frame in the direction of the z-axis. The system possesses an internal resonance between its slowest transverse bending modes, whose frequencies indeed clock at [Hz] and [Hz]. We consider transverse vibrations in the out of plane direction—the z-axis in Fig. 6(a). The available observable is the velocity of the inner beam tip, measured via laser scanner vibrometry (PSV400, Polytec Inc.). Transient vibrations are recorded for 3 seconds at a sampling rate of [Hz].
A modally-tuned impulse hammer (PCB 086C01, PCB Piezotronics, Inc.) is used to excite transverse vibrations from three different impact locations in Fig. 6(c), so that our dataset features trajectories (four per impact location), shown in Fig. 6(c). We label these trajectories as where refers to the location and to the test number. Time-frequency analyses of the velocity signals, two of which are reported in Fig. 6(d,e), show that only the two slowest frequencies are present in the signal, so that the time responses can be well approximated by the slowest four-dimensional SSM of the system. The impact locations, the hammer tip and the forcing amounts were selected to achieve a sufficient trajectory diversity in the dataset without exciting further structural modes. For constructing an SSM-reduced model, we truncate the velocity signals after the hammer impact, use ten trajectories for training and leave two trajectories for testing.
The minimal embedding dimensions (nine for a four-dimensional manifold) fails to produce accurate reduced-order models (the NMTE error amounts to more than ). We therefore augment the delay embedding space so that each embedding vector captures approximately two cycles of the slowest oscillation. This procedure yields a -dimensional delay embedding space. The result of our identification remains robust if we consider more embedding dimensions. A linear approximation to the embedded SSM has a good accuracy and our automated normal form algorithm, after estimating linearized eigenvalues, identifies a resonance among them. Defining , we obtain from SSMLearn the cubic SSM-reduced polar normal form
This data-driven model reconstructs both test trajectories with an average NMTE error, cf. Fig. 7(a,b). The decay of the slow mode amplitude and that of the fast one are shown in Fig. 7(c). Due to modal interactions, these decays are not monotone. From the plot, we notice a great diversity of decays depending on the impact location, and location three (the closest to the inner beam tip) is characterized by the highest amplitudes variability.
Figure 8(a) shows energy repartition among the modes for the third impact location. This repartition is defined as the instantaneous ratio between the amplitude of a mode and the amplitude sum. Clearly, the slow mode tends to accumulate energy over time, while the fast mode dissipates it quickly. These trends are not monotonic, showing simultaneous and opposite changes in growth/decay rates, which implies that the faster mode is absorbing energy from the slower one. This behavior can also be inferred by the instantaneous properties illustrated in 8(b,c). The uncoupled limit of the oscillators suggests that the modes admit frequency softening and damping intensification when the oscillation amplitude increase. This is consistent with typical observations of jointed structures . In particular, the fast mode is coupled to the slow one and its damping undergoes consistent variation, becoming also negative for some times . Note that nonlinearity and coupling can be reduced at higher bolt torques, which, in the beam assembly used here, corresponds to tightening of the middle bolt. Coupling revealed by SSMLearn suggests that nearly decoupled modal oscillator models employed elsewhere [74, 75] are only valid for high bolt torques and small frictional slip, i.e., weak contact nonlinearities.
In addition to measuring decaying vibrations, we also perform some forced testing near the linearized frequency of the slow mode. We trigger forced responses in near-resonance with the slow (first bending) mode by using the Brüel Kjær 4810 shaker shown in Fig. 6(b), mounted on one end of the linear spring, and acquire velocity response from the tip of the inner beam by laser vibrometry. We also monitor the amplitude of shaker tip velocities and keep them constant while sweeping the frequencies around the first bending mode. In that sense, the response we obtain can be seen as transmissibilty rather than a classic frequency response curve, with forcing amplitudes kept constant throughout frequency sweeps. Starting from forced velocity time histories, we estimate the normal form forcing to be added to the vector field (17) as in (9). The resulting predictions are in very good agreement with experimental measurements in this weakly nonlinear regime, as shown in Fig. 8(d). Deeper analyses on forced responses are currently under investigation. Thanks to feedback loops used to track forcing, forced response curves can be extracted with improved accuracy, especially at non-linearizable amplitudes.
We have reviewed a general methodology for constructing sparse reduced-order models for potentially high-dimensional, nonlinear mechanical systems from data. Our approach constructs normal forms on attracting spectral submanifolds (SSMs), which are the smoothest nonlinear continuation of spectral subspaces of the linearized dynamics. Implemented in the publicly available Matlab® code SSMLearn , our algorithm takes generic observable data as input to identify robust and predictive nonlinear models that also capture for non-linearizable dynamics. SSM theory offers a systematic basis for model reduction and allows a simplification of the reduced dynamics via normal forms, which are particularly insightful for mechanical systems. Indeed, SSM-reduced models can handle multi-modal interactions, identify amplitude-dependent damping and frequency, and predict the forced structural response.
We have illustrated SSM-reduced modeling in numerical and experimental case studies, featuring different types of observables, nonlinearities and SSM dimensions. Specifically, we have discussed different dynamical regimes and the relevance of slow SSMs in a chain of oscillators, derived a reduced-order model from digital image correlation measurements of the Brake-Reuss beam, and unfolded the internally resonant dynamics of a tester structure, also predicting forced responses. These examples were analyzed using the open-source Matlab® package SSMLearn  that performs data-driven, SSM-based model reduction starting from vibrations data. This algorithm only requires a minimal number of input parameters: the SSM dimension, the polynomial order for SSM parametrization and the polynomial order of the reduced dynamics. The SSM dimension is either known a priori from targeted experiments (e.g. resonance decay) or can be estimated via time-frequency signal processing analysis of the input data. This makes our method a parametric complement to non-parametric identification tools. Polynomial orders can be adjusted to improve accuracy, noting that excessive orders may lead to overfitting. With the help of the numerical continuation core coco  included in SSMLearn , users can compute forced response curves or design nonlinear control strategies from the identified nonlinear models.
Further examples, both numerical and experimental, with detailed code are available in the SSMLearn repository , which is also suitable for high-dimensional fluid flows and fluid structure interaction problems . Current limitations of the present approach include weaker performance for large forcing amplitudes. These appear, for example, in the Brake-Reuss beam experiments of [7, 8], which we expect to capture only with a more refined forced-reduced dynamics and improved calibration procedures. The same requirement holds for more complicated forcing types (e.g., quasi-periodic or random), which are relevant in structural dynamics. We are addressing these challenges in ongoing work.
-  M. Peeters, G. Kerschen, and J.C. Golinval. Dynamic testing of nonlinear vibrating structures using nonlinear normal modes. Journal of Sound and Vibration, 330(3):486–509, 2011.
-  M. Peeters, G. Kerschen, and J.C. Golinval. Modal testing of nonlinear vibrating structures based on nonlinear normal modes: experimental demonstration. Mechanical Systems and Signal Processing, 25(4):1227–1247, 2011.
-  M. Kurt, M. Eriten, D.M. McFarlan, L.A. Bergman, and A.F. Vakakis. Strongly nonlinear beats in the dynamics of an elastic system with a strong local stiffness nonlinearity: Analysis and identification. Journal of Sound and Vibration, 333(7):2054–2072, 2014.
-  D.A. Ehrhardt and M.S. Allen. Measurement of nonlinear normal modes using multi-harmonic stepped force appropriation and free decay. Mechanical Systems and Signal Processing, 76-77:612–633, 2016.
-  L. Renson, A. Gonzalez-Buelga, D.A.W. Barton, and S.A. Neild. Robust identification of backbone curves using control-based continuation. Journal of Sound and Vibration, 367:145–158, 2016.
-  K.J. Moore, M. Kurt, M. Eriten, D.M. McFarlan, L.A. Bergman, and A.F. Vakakis. Direct detection of nonlinear modal interactions from time series measurements. Mechanical Systems and Signal Processing, 125:311–329, 2019. Exploring nonlinear benefits in engineering.
-  W. Chen, D. Jana, A. Singh, M. Jin, M. Cenedese, G. Kosova, M.W.R. Brake, C.W. Schwingshackl, S. Nagarajaiah, K. Moore, and J.P. Noël. Measurement and identification of the nonlinear dynamics of a jointed structure using full-field data; Part I - Measurement of nonlinear dynamics. Mechanical Systems and Signal Processing, 166:108401, 2022.
-  M. Jin, G. Kosova, M. Cenedese, W. Chen, D. Jana, A. Singh, M.W.R. Brake, C.W. Schwingshackl, S. Nagarajaiah, K. Moore, and J.P. Noël. Measurement and identification of the nonlinear dynamics of a jointed structure using full-field data; Part II - Nonlinear system identification. Mechanical Systems and Signal Processing, 166:108402, 2022.
-  J. Awrejcewicz, V.A. Krys’ko, and A.F. Vakakis. Order Reduction by Proper Orthogonal Decomposition (POD) Analysis, pages 279–320. Springer, Berlin, Heidelberg, 2004.
-  P. Holmes, J.L. Lumley, G. Berkooz, and C.W. Rowley. Turbulence, Coherent Structures, Dynamical Systems and Symmetry. Cambridge Monographs on Mechanics. Cambridge University Press, 2 edition, 2012.
-  K. Lu, Y. Jin, Y. Chen, Y. Yang, L. Hou, Z. Zhang, Z. Li, and C. Fu. Review for order reduction based on proper orthogonal decomposition and outlooks of applications in mechanical systems. Mechanical Systems and Signal Processing, 123:264–297, 2019.
-  S. Hijazi, G. Stabile, A. Mola, and G. Rozza. Data-driven POD-Galerkin reduced order model for turbulent flows. Journal of Computational Physics, 416:109513, 2020.
-  P.J. Schmid. Dynamic mode decomposition of numerical and experimental data. Journal of Fluid Mechanics, 656:5–28, 2010.
-  J.N. Kutz, S.L. Brunton, B.W. Brunton, and J.L. Proctor. Dynamic Mode Decomposition. Society for Industrial and Applied Mathematics, Philadelphia, PA, 2016.
-  K.K. Chen, J.H. Tu, and C.W. Rowley. Variants of dynamic mode decomposition: Boundary condition, Koopman, and Fourier analyses. Annual Review of Fluid Mechanics, 22(6):887–915, 2012.
-  M.O. Williams, C.W. Rowley, and I.G. Kevrekidis. A kernel-based method for data-driven koopman spectral analysis. Journal of Computational Dynamics, 2(2):247–265, 2015.
-  C.W. Rowley and S.T.M. Dawson. Model reduction for flow analysis and control. Annual Review of Fluid Mechanics, 49(1):387–417, 2017.
-  A. Alla and J.N. Kutz. Nonlinear model order reduction via dynamic mode decomposition. SIAM Journal on Scientific Computing, 39(5):B778–B796, 2017.
B. Lusch, J.N. Kutz, and S.L. Brunton. Deep learning for universal linear embeddings of nonlinear dynamics.Nature Communications, 9:4950, 2018.
-  M. Budišić, R. Mohr, and I. Mezić. Applied koopmanism. Chaos: An Interdisciplinary Journal of Nonlinear Science, 22(4):047510, 2012.
-  I. Mezić. Analysis of fluid flows via spectral properties of the koopman operator. Annual Review of Fluid Mechanics, 45(1):357–378, 2013.
-  S. Bagheri. Koopman-mode decomposition of the cylinder wake. Journal of Fluid Mechanics, 726:596–623, 2013.
-  S.L. Brunton, B.W. Brunton, J.L. Proctor, and J.N. Kutz. Koopman invariant subspaces and finite linear representations of nonlinear dynamical systems for control. PLOS ONE, 11:1–19, 2016.
-  J. Page and R.R. Kerswell. Koopman mode expansions between simple invariant solutions. Journal of Fluid Mechanics, 879:1–27, 2019.
E. Kaiser, J.N. Kutz, and S.L. Brunton. Data-driven discovery of koopman eigenfunctions for control.Machine Learning: Science and Technology, 2(3):035023, 2021.
-  J. Page and R.R. Kerswell. Koopman analysis of burgers equation. Phys. Rev. Fluids, 3:071901, Jul 2018.
-  K.F.R.S. Pearson. LIII. On lines and planes of closest fit to systems of points in space. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 2(11):559–572, 1901.
-  C.M. Bishop. Pattern Recognition and Machine Learning. Information Science and Statistics. Springer-Verlag New York, 2006.
-  S.T. Roweis and L.K. Saul. Nonlinear dimensionality reduction by locally linear embedding. Science, 290(5500):2323–2326, 2000.
-  Y. Bengio, A. Courville, and P. Vincent. Representation learning: a review and new perspectives. IEEE Transactions on Pattern Analysis and Machine Intelligence, 35(8):1798–1828, 2013.
-  J.-C. Loiseau, S.L. Brunton, and B.R. Noack. From the POD-Galerkin method to sparse manifold models. In P. Benner, editor, Handbook of Model Order Reduction, Volume II: Applications, pages 1–47. De Gruyter GmbH, 2019.
-  I. Goodfellow, Y. Bengio, and A. Courville. Deep Learning. MIT Press, 2016.
-  K. Champion, B. Lusch, J.N. Kutz, and S.L. Brunton. Data-driven discovery of coordinates and governing equations. Proceedings of the National Academy of Sciences, 116(45):22445–22451, 2019.
-  R. Tibshirani. Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society. Series B (Methodological), 58(1):267–288, 1996.
-  S.L. Brunton, J.L. Proctor, and J.N. Kutz. Discovering governing equations from data by sparse identification of nonlinear dynamical systems. Proceedings of the National Academy of Sciences, 113(15):3932–3937, 2016.
-  A. Ben Abdessalem, N. Dervilis, D.J. Wagg, and K. Worden. Model selection and parameter estimation of dynamical systems using a novel variant of approximate bayesian computation. Mechanical Systems and Signal Processing, 122:364–386, 2019.
-  S.A. Billings. Nonlinear System Identification: NARMAX Methods in the Time, Frequency, and Spatio-Temporal Domains. Springer, 2013.
-  D. Hartman and L.K. Mestha. A deep learning framework for model reduction of dynamical systems. In 2017 IEEE Conference on Control Technology and Applications (CCTA), pages 1917–1922, 2017.
-  H.F.S. Lui and W.R. Wolf. Construction of reduced-order models for fluid flows using deep feedforward neural networks. Journal of Fluid Mechanics, 872:963–994, 2019.
M. Raissi, P. Perdikaris, and G.E. Karniadakis. Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations.Journal of Computational Physics, 378:686–707, 2019.
-  G.E. Karniadakis, I.G. Kevrekidis, L. Lu, P. Perdikaris, S. Wang, and L. Yang. Physics-informed machine learning. Nature Reviews Physics, 3:422–440, 2021.
-  G. Haller and S. Ponsioen. Nonlinear normal modes and spectral submanifolds: existence, uniqueness and use in model reduction. Nonlinear Dynamics, 86(3):1493–1534, 2016.
-  G. Haller and S. Ponsioen. Exact model reduction by a slow-fast decomposition of nonlinear mechanical systems. Nonlinear Dynamics, 90:617–647, 2017.
-  R. Szalai, D. Ehrhardt, and G. Haller. Nonlinear model identification and spectral submanifolds for multi-degree-of-freedom mechanical vibrations. Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 473:20160759, 2017.
-  S. Jain, P. Tiso, and G. Haller. Exact nonlinear model reduction for a von Kármán beam: slow-fast decomposition and spectral submanifolds. Journal of Sound and Vibration, 423:195–211, 2018.
-  S. Ponsioen, T. Pedergnana, and G. Haller. Automated computation of autonomous spectral submanifolds for nonlinear modal analysis. Journal of Sound and Vibration, 420:269–295, 2018.
-  T. Breunung and G. Haller. Explicit backbone curves from spectral submanifolds of forced-damped nonlinear mechanical systems. Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 474:20180083, 2018.
-  S. Ponsioen, T. Pedergnana, and G. Haller. Analytic prediction of isolated forced response curves from spectral submanifolds. Nonlinear Dynamics, 98:2755–2773, 2019.
-  S. Ponsioen, S. Jain, and G. Haller. Model reduction to spectral submanifolds and forced-response calculation in high-dimensional mechanical systems. Journal of Sound and Vibration, 488:115640, 2020.
-  S. Jain and G. Haller. How to compute invariant manifolds and their reduced dynamics in high-dimensional finite-element models? Nonlinear Dynamics, in press, 2021.
-  S. Jain, T. Thurnher, M. Li, and G. Haller. SSMTool-2.0: Computation of invariant manifolds & their reduced dynamics in high-dimensional mechanics problems, www.georgehaller.com, 2021.
-  R. Szalai. Invariant spectral foliations with applications to model order reduction and synthesis. Nonlinear Dynamics, 101:2645–2669, 2020.
-  M. Cenedese, J. Axås, B. Bäuerlein, K. Avila, and G. Haller. Data-driven modeling and prediction of non-linearizable dynamics via spectral submanifolds. Submitted, 2021.
-  M. Cenedese, J. Axås, and G. Haller. SSMLearn: Data-driven reduced order models for nonlinear dynamical systems. To be released, 2021.
-  J.P. Noël and G. Kerschen. Nonlinear system identification in structural dynamics: 10 more years of progress. Mechanical Systems and Signal Processing, 83:2–35, 2017.
-  M. Feldman. Hilbert transform in vibration analysis. Mechanical Systems and Signal Processing, 25(3):735–802, 2011.
-  M. Jin, W. Chen, M.R.W. Brake, and H. Song. Identification of instantaneous frequency and damping from transient decay data. Journal of Vibration and Acoustics, 142(5):051111, 2020.
-  K.J. Moore, M. Kurt, M. Eriten, D.M. McFarland, L.A. Bergman, and A.F. Vakakis. Wavelet-bounded empirical mode decomposition for measured time series analysis. Mechanical Systems and Signal Processing, 99:14–29, 2018.
-  X. Cabré, E. Fontich, and R. De la lLave. The parameterization method for invariant manifolds I: manifolds associated to non-resonant subspaces. Indiana University Mathematics Journal, 52(2):283–328, 2003.
-  X. Cabré, E. Fontich, and R. De la lLave. The parameterization method for invariant manifolds II: regularity with respect to parameters. Indiana University Mathematics Journal, 52(2):329–360, 2003.
-  X. Cabré, E. Fontich, and R. De la lLave. The parameterization method for invariant manifolds III: overview and applications. Journal of Differential Equations, 218(2):444–515, 2005.
-  A. Haro, M. Canadell, J. Figueras, and J.M. Luque, A. Mondelo. The Parameterization Method for Invariant Manifolds, volume 195 of Applied Mathematical Sciences. Springer International Publishing, 1st edition, 2016.
-  N. Fenichel. Asymptotic stability with rate conditions. Indiana University Mathematics Journal, 23(12):1109–1137, 1974.
-  T. Sauer, J.A. Yorke, and M. Casdagli. Embedology. Journal of Statistical Physics, 65:579–616, 1997.
-  F. Takens. Detecting strange attractors in turbulence. In D. Rand and L. Young, editors, Dynamical Systems and Turbulence, Warwick 1980, pages 366–381. Springer Berlin Heidelberg, Berlin, Heidelberg, 1981.
-  J. Guckenheimer and P.J. Holmes. Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, volume 42 of Applied Mathematical Sciences. Springer-Verlag New York, 1983.
-  J. Murdock. Normal Forms and Unfoldings for Local Dynamical Systems. Springer Monographs in Mathematics. Springer-Verlag New York, 2003.
-  H. Poincaré. Les Méthodes Nouvelles de la Mécanique Céleste. Gauthier-Villars et Fils, Paris, 1892.
-  J.E. Marsden and M. McCracken. The Hopf Bifurcation and Its Applications, volume 19 of Applied Mathematical Sciences. Springer-Verlag New York, 1976.
-  M. Cenedese and G. Haller. How do conservative backbone curves perturb into forced responses? A Melnikov function analysis. Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 476:20190494, 2020.
-  H. Dankowicz and F. Schilder. Recipes for Continuation. Society for Industrial and Applied Mathematics, 2013.
-  M.R.W. Brake. The mechanics of jointed structures: recent research and open challenges for developing predictive models for structural dynamics. Springer International Publishing, 2018.
-  T.P. Sapsis, D.D. Quinn, A.F. Vakakis, and L.A. Bergman. Effective stiffening and damping enhancement of structures with strongly nonlinear local attachments. Journal of Vibration and Acoustics, 134(1):011016, 2012.
-  M. Eriten, M. Kurt, G. Luo, D.M. D. McFarland, L.A. Bergman, and A.F. Vakakis. Nonlinear system identification of frictional effects in a beam with a bolted joint connection. Mechanical Systems and Signal Processing, 39(1):245–264, 2013.
-  D.J. Segalman, M.S. Allen, M. Eriten, and K. Hoppman. Experimental assessment of joint-like modal models for structures. Proceedings of the ASME 2015 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 8: 27th Conference on Mechanical Vibration and Noise. Boston, Massachusetts, USA. August 2–5, 2015. V008T13A025. ASME.