## Experimental continuation in nonlinear dynamics: recent advances and future challenges

Ghislain Raze, Ga¨ etan Abeloos and Ga¨tan Kerschen e

* Aerospace and Mechanical Engineering Department, University of Li`ge, All´e de la e e D´ ecouverte, 9, Li`ge, 4000, Belgium. e

*Corresponding author(s). E-mail(s): g.raze@uliege.be;

## Abstract

Experimental continuation encompasses a set of methods that combine control and continuation to obtain the full bifurcation diagram of a nonlinear system experimentally, including responses that would be unstable in the system without feedback control. Such control-based methods thus allow the experimenter to directly and exhaustively explore the dynamics of the system without the need for a good mathematical model. The objective of this paper is twofold, namely (i) to review and present the state-of-the-art methods in a unified manner and (ii) to introduce a novel experimental derivative-free arclength continuation procedure, termed arclength control-based continuation. These methods are also demonstrated and compared on an electronic Duffing oscillator and a clamped thin plate featuring geometrical nonlinearity. Finally, the current state of the art is reflected upon, and the challenges lying ahead for this growing field are discussed.

Keywords: Experimental nonlinear dynamics, Control-based methods, Experimental continuation, Control-based continuation, Phase-locked loop, Response-controlled testing

## 1 Introduction

The theory of nonlinear dynamical systems was developed based on Poincar´ e's seminal work. It witnessed extraordinary advances during the 20th century thanks to Arnold, Holmes, Lorenz, Ruelle and Smale to name a few. A couple of decades ago, the theory spread across engineering thanks to a series of reference monographs thoroughly discussing the different phenomena and attractors nonlinear systems can exhibit. Nayfeh and Mook [1] applied perturbation methods to study nonlinear phenomena in single- and multi-degree-of-freedom systems. Guckenheimer and Holmes [2] adopted a different, geometrical viewpoint, appearing as an ideal companion to the perturbation approach. Kuznetsov [3] published a complete treatise on bifurcations.

Since the 1970s, impressive progress was also realized in computational nonlinear dynamics with the development of nonlinear finite element methods [4] and of numerical continuation (NC) [5-8]. Numerical continuation amounts to calculating a bifurcation diagram, a key concept in nonlinear dynamics [9] representing the different branches of responses of a nonlinear system together with their stability and bifurcations as one or more parameters are varied. Today, the calculation of bifurcation diagrams of complex numerical models is within reach [10-13]. Noteworthy are also the advances in the area of nonlinear model reduction [14-17]. We also

note that there now exists the possibility to perform direct, data-driven reduction onto invariant manifolds [18].

Although not all challenges have been overcome yet, the theoretical understanding of nonlinear dynamical phenomena and their prediction using numerical models have reached a high level of maturity. Conversely, experimental nonlinear dynamics [19] received comparatively far less attention and is the main thrust of this paper. In this context, nonlinear system identification (NSI) which focuses on the development of mathematical models directly from experimental data is the most popular approach to date in the literature. There exists a wide variety of techniques see, e.g., [20, 21]. Once a good model of the system has been obtained thanks to NSI, numerical continuation techniques can then be exploited to calculate a bifurcation diagram, as carried out, e.g., in [11, 22].

Despite the recent advances in the field, NSI remains a particularly difficult exercise and imposes a number of challenges and constraints to the experimenter. An elegant, model-less alternative to NSI is experimental continuation which leverages feedback control for the identification of a bifurcation diagram. This novel experimental strategy represents a true paradigm shift in nonlinear dynamics. There are two key differences with numerical continuation, namely (i) bifurcation diagrams have to be obtained in a model-less manner and (ii) feedback control must be applied to stabilize unstable periodic orbits which otherwise cannot be measured experimentally. The first experimental continuation method, control-based continuation (CBC), pioneered by Sieber et al [23, 24], can be seen as an experimental equivalent to pseudo-arclength continuation. The fact that CBC requires to calculate derivatives using finite differences from experimental data, a computationallyintensive and highly sensitive-to-noise process, motivated researchers to develop methods which do not rely on experimental derivatives. The three most effective derivative-free experimental continuation methods to date, namely the simplified CBC (SCBC) method [25], phase-locked loops (PLLs) [26] and response-controlled testing (RCT) [27], resort to natural parameter continuation, which simply increases the bifurcation parameter using the previous solution as an initial guess for the next solution. Thanks to these advances, the frequency responses and backbone curves of nonlinear systems around their primary resonances can now be identified in a robust and effective manner.

The objective of this feature article is to offer the first extensive review of the developments in the EC area since its invention, to present the three derivative-free methods in an unified manner and to introduce a novel experimental continuation method, termed arclength control-based continuation (ACBC) [28], the first derivative-free arclength continuation scheme.

The paper is organized as follows. Bifurcation diagrams and the means to compute them are introduced in Section 2. Section 3 provides a brief introduction to basic concepts in feedback control. Section 4 presents the core ideas of EC and details its first forms. Section 5 then reviews the most recent methods and presents them under a unified framework. The novel ACBC approach is presented in Section 6. All the methods are compared experimentally with an electronic Duffing oscillator and a clamped plate featuring geometrical nonlinearity in Sections 7 and 8, respectively. Finally, the conclusions of the present study and the challenges lying ahead of us are summarized in Section 9.

## 2 Bifurcation diagrams: a central tool in nonlinear dynamics

Nonlinear dynamical systems usually behave in a complex manner, requiring adequate tools to study their behavior. Bifurcation diagrams are ubiquitously used in nonlinear dynamics for that purpose. These diagrams typically present how one quantity characterizing the response of the system under study is influenced by a parameter, often called a bifurcation parameter . In this work, we focus on the periodic solutions of nonlinear systems and their stability (but fixed points can be seen as a special subcase).

In this section, the main features of analytical and numerical approaches to draw out bifurcation diagrams are summarized, and used to illustrate the challenges faced in experimental nonlinear dynamics.

## 2.1 Periodic solutions

We consider a generic, non-autonomous system governed by the set of first-order ordinary differential equations (ODEs)

$$\dot { \mathbf x } ( t ) = \mathbf f ( \mathbf x ( t ), t ), \quad \quad ( 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where x is a vector gathering the states of the system under consideration, t is the time, and f is a state evolution law. We note that this framework encompasses the case of a set of second-order ODEs generally used in vibration problems

$$\begin{array} { c } M \ddot { q } ( t ) + C \dot { q } ( t ) + K \dot { q } ( t ) + \mathbf f _ { n l } ( \mathbf q ( t ), \dot { q } ( t ) ) = \mathbf f _ { \text{ext} } ( t ), \\ \cdot \cdot \cdot \cdot \cdot \end{array} \quad \begin{array} { c } \ r e s \mathbb { F } \\ \cdot \cdot \end{array} \quad \begin{array} { c } \ r e s \mathbb { F } \\ \cdot \cdot \cdot \end{array} \quad \begin{array} { c } \ r e s \mathbb { F } \end{array}$$

where M C , and K are linear structural mass, damping and stiffness matrices, respectively, f nl is the vector of nonlinear forces and f ext is the vector of external forces. Indeed, this set of second-order ODEs can be recast into a set of first-order ODEs (provided that the mass matrix is invertible) using

$$\begin{array} { c c c } ( \text{province} \, \text{and} \, \text{me mass matrix is inverse} ) \, \text{using} \\ \dot { x } = \begin{bmatrix} \dot { q } \\ \ddot { q } \end{bmatrix} = & \begin{bmatrix} 0 & I \\ - M ^ { - 1 } K & - M ^ { - 1 } C \end{bmatrix} \begin{bmatrix} q \\ \dot { q } \end{bmatrix} & \begin{bmatrix} \text{$\text{$sn}$} \\ \text{var} \\ = f ( x, t ). \end{bmatrix} & ( 3 ) & ( \text{to} \\ \text{and} \\ \text{lat} \end{bmatrix}$$

When the right-hand side of Equation (1) is periodic in time, periodic solutions x generally exist, satisfying

$$\mathbf x ( t + T ) = \mathbf x ( t ), \quad \forall t \in \mathbb { R } \quad \ \ ( 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where T is a period of the solution, being related to its angular frequency ω by ωT = 2 π .

## 2.2 Analytical solutions

Equation (1) with given initial conditions has a unique solution provided that f satisfies mild smoothness conditions given by Picard's existence theorem. In some cases, analytical insight into these solutions can be gained.

## 2.2.1 Exact solutions

Due to the very nature of nonlinear systems, exact analytical solutions seldom exist. A few special or nongeneric cases can be treated, such as the unforced, undamped pendulum or Duffing oscillator with the use of Jacobi elliptic functions [29].

For illustration, we consider the free response q of an undamped Duffing oscillator, a typical academic example qualitatively representative of geometrically nonlinear structures, which is governed by the following ODE:

$$m \ddot { q } ( t ) + k q ( t ) + k _ { 3 } q ^ { 3 } ( t ) = 0, \quad \ ( 5 )$$

where m k , and k 3 are the mass, linear stiffness and cubic stiffness coefficients, respectively, of the oscillator and an overdot denotes a derivation with respect to time t . Dividing Equation (5) by m ,

$$\ddot { q } ( t ) + \omega _ { 0 } ^ { 2 } q ( t ) + \alpha _ { 3 } q ^ { 3 } ( t ) = 0 \text{ \quad \ \ } ( 6 )$$

with ω 2 0 = k/m and α 3 = k /m 3 . The exact solution of Equation (6) is given by

$$q ( t ) = A \text{sn} ( \Omega t + b | \kappa ), \text{ \quad \ \ } ( 7 )$$

where A is the amplitude of the free oscillations, sn is the elliptic sine function, b is a phase-like variable to be adjusted to the initial conditions (together with A ), Ω is a frequency-like variable and κ is the elliptic parameter [29]. These two latter parameters are related to A by

$$\begin{array} { c c c c } \cdot & & \cdot \\ \Omega = \sqrt { \omega _ { 0 } ^ { 2 } + \frac { \alpha _ { 3 } A ^ { 2 } } { 2 } }, & & \kappa = - \frac { \alpha _ { 3 } A ^ { 2 } } { 2 \Omega ^ { 2 } }. & & ( 8 ) \\ - & & \cdot & & \cdot & \cdot \end{array}$$

## 2.2.2 Perturbation methods

Perturbation methods constitute a means to obtain approximate analytical solutions to nonlinear problems. They rely on a small parameter, usually denoted ϵ . In most cases, the external forcing, response amplitude, nonlinear coefficients and/or damping coefficients are assumed to scale with ϵ . When ϵ = 0, a (generally known) solution forms the basis around which a power series expansion in terms of ϵ is performed.

Perturbation methods are a powerful tool to study nonlinear systems analytically. The most famous ones are the method of multiple scales [1] and the averaging method [2]. These methods can provide insight into the system's behavior, but are inherently restricted by the limitations associated with power series expansions, and generally break down when one or several assumed small parameters become too large.

## 2.3 Numerical procedures

If analytical solutions are unavailable, intractable or inaccurate, periodic solutions for a fixed ω can be found by numerically integrating Equation (1) until the transients die out. This approach is generally inefficient, and more sophisticated approaches can be used, such as the shooting [30], harmonic balance [11, 31] and orthogonal collocation [32] methods. These methods cast the ODEs into a set of algebraic equations that can be solved by well-known numerical routines such as the Newton-Raphson method. This usually yields accurate results.

To gain a broader understanding of the system's behavior, a bifurcation diagram can be built by varying a bifurcation parameter. To find the family of periodic solutions that form a complete bifurcation diagram, numerical continuation can be used, as will be presented next.

## 2.4 Numerical continuation

Numerical continuation (also known as branch tracing or path-following method) denotes a class of numerical methods used to solve a system of algebraic equations with one (or more) free parameter. The solutions of these equations generally form a set whose dimension is equal to the number of bifurcation parameters. Since the most common case considers a single parameter, this set is generically a curve. We treat this case without loss of generality, because the case with multiple bifurcation parameters can be treated, e.g., by sampling the values of all but one parameter, and using a one-parameter continuation for all samples, thereby effectively 'slicing' the solution set.

Numerical continuation is an incremental procedure. The main idea behind it is to build a branch of solutions point by point (or piece by piece), using a constructive approach where the information at the previous point (or piece) of the curve is known to be well estimated and is used to find a suitable initial guess for the next one. This defines the prediction step. This initial guess is then refined to ensure that it lies on the sought curve within a prescribed tolerance, which defines the correction step. Figure 1 schematically depicts the main ideas behind this predictor-corrector scheme, where x is one of the unknowns, and λ is a bifurcation parameter.

Fig. 1 : Schematics of a general continuation method.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-571549cc-6914-4b1c-82b1-14095f27d535.jpg)

Numerical continuation is a mature tool that is now treated in a variety of reference books and exploited in a large number of bifurcation analysis software, such as AUTO [33], MATCONT [7], COCO [12] and MANLAB [34]. In this work, we will focus on the two most popular variants of numerical continuation, namely natural parameter and arclength continuations. There also are numerous books discussing this topic, such as [6, 8, 30, 32, 35].

## 2.4.1 Natural parameter continuation

Natural parameter continuation (also called sequential continuation) is the simplest form of continuation. It uses the bifurcation parameter of the problem to parametrize the solution curve. The implicit function theorem then gives conditions under which this choice is guaranteed to be sensible, i.e., when the solution curve is indeed (locally) parametrizable by the bifurcation parameter. When the solution is known for a given value of the parameter, the prediction step simply consists in increasing the parameter value, using the previous solution as initial guess for the next point. The correction phase generally consists in using an iterative solver such as the NewtonRaphson method. Figure 2 represents a general schematics of a natural parameter continuation method.

Fig. 2 : Schematics of the natural parameter continuation method.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-8a9ded54-2514-4b7f-adbf-2003f757dc52.jpg)

Natural parameter continuation fails at points where the implicit function theorem does not hold when one parametrizes a curve with the continuation parameter. Two notable cases are limit points and branch points, represented in Figures 3a and 3b, respectively. In the presence of the former, natural parameter continuation fails most of the time because there is no nearby solution at the imposed parameter value, as shown in Figure 3a. With branch points, it is generally hard to predict and control which path the method will follow close to the bifurcation; this depends on the step length and on the problem at hand (cf. Figure 3b). To address these issues, more elaborate continuation approaches were developed.

Fig. 3 : Typical scenarii where the natural parameter continuation fails or has an unpredictable outcome: folds (a) and branch points (b).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-702cff18-f9b9-4a46-b197-e56fd03d744c.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-a68a5b38-4b65-4630-9090-c8ec3294c81c.jpg)

## 2.4.2 Arclength continuation

The first ideas behind arclength continuation can be traced back to the seminal work of Haselgrove in the early sixties [36], and were extensively developed and popularized by Keller [5]. It was (seemingly independently) popularized nearly twenty years later in the field of static mechanics with the famous works of Riks [37] and Crisfield [38], who respectively proposed what is now known as the pseudo-arclength and arclength continuation procedures. Adaptations to dynamics came ten years later using the harmonic balance [39] and shooting [40] methods. The underlying idea of these methods is to parametrize the curve with an (approximation of its) arclength.

These methods both use a more refined prediction step with either a tangent or secant predictor. The correction is then performed either on a hyperplane whose normal is the prediction vector (Figure 4a), or on a hypersphere centered around the previous point and whose radius is equal to the prediction step length (Figure 4b). Mathematically, λ is let free, and the correction locus adds another constraint to the problem to make it wellposed. Clearly, these methods offer better robustness than natural parameter continuation when the curve folds. In addition, the tangent predictor generally makes the method more likely to stay on the same branch when a branching behavior is encountered. These are some of the reasons for their current popularity and ubiquitous presence in numerical continuation packages.

## 2.5 Stability

Equipped with the possibility to draw bifurcation diagrams, we now have to wonder whether the obtained family of solutions will be observable in the nonlinear system under study. By contrast with linear time-invariant (LTI) systems for which solutions are either all stable or all unstable, nonlinear systems can possess solutions of different stability. In addition, while unstable solutions of LTI systems grow unbounded with time, trajectories repelled by a nearby unstable nonlinear solution have the additional possibility to converge to other stable solutions.

In an experimental context, the unstable solutions of a system are impossible to observe in steady-state conditions without additional control. As for numerical aspects, unstable solutions can be computed, but an important question is to know whether this solution is stable or not. To answer this question, we look at the evolution of infinitesimally small perturbations around the equilibrium, which gives us results about the

Fig. 4 : Schematics of the pseudo-arclength (a) and arclength (b) methods.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-5b6c7b19-8668-42a4-a72b-17cee45da0f4.jpg)

asymptotic stability of a solution. For more elaborate discussions on stability, we refer to, e.g., [2, 41].

We now assume that we know a periodic solution of Equation (1) (obtained analytically or numerically) and denote it by x ∗ . A small perturbation s is added to the equilibrium x ∗ , i.e.,

$$\mathbf x ( t ) = \mathbf x _ { * } ( t ) + \mathbf s ( t ), \quad | \mathbf s ( t ) | < < | \mathbf x _ { * } ( t ) | \quad ( 9 ) \quad 2. \mathbf t$$

Inserting this perturbed solution into Equation (1) and expanding the result in Taylor series around s ( ) = 0, we obtain t

$$\sigma _ { \vee \nu } - \vartheta, \, w c \, \varnothing \text{vam\ } \\ \dot { x } _ { * } ( t ) + \dot { s } ( t ) & = \mathbf f ( \mathbf x _ { * }, t ) + \frac { \partial } { \partial \mathbf x } \mathbf f ( \mathbf x, t ) \Big | _ { \mathbf x = \mathbf x _ { * } } \mathbf s ( t ) \quad \text{can} \\ + O ( | \mathbf s ( t ) | ^ { 2 } ). & ( 1 0 ) \\ \text{Neglecting terms of orders higher than linear in } | \mathbf s | \quad \text{.}$$

Neglecting terms of orders higher than linear in | s | and accounting for the fact that x ∗ is a solution of

Equation (1), one gets the first-order variational problem

$$\tau \pi \pi \\ \dot { s } ( t ) = \frac { \partial } { \partial x } \mathbf f ( \mathbf x, t ) \Big | _ { \mathbf x = \mathbf x _ { * } ( t ) } \, \mathbf s ( t ). \quad ( 1 1 ) \\ \tilde { s } \, \text{is a linear, time-varying equation featuring} \\ \text{Jacobian matrix} \, \mathbf f / \partial \mathbf x \, \text{and governing the evo-}$$

This is a linear, time-varying equation featuring the Jacobian matrix ∂ /∂ f x and governing the evolution of the (small) perturbations. The solution of such an equation cannot be known in closedform in general, and requires numerical integration. Nevertheless, the time-varying part of this equation is periodic, and this periodicity entails special consequences which are described by the Floquet theory [9]. The variational equation can be solved (using numerical time integration) with its principal fundamental matrix Φ ( ) determined t by

$$\dot { \Phi } ( t ) = \frac { \partial } { \partial x } \mathfrak { f } ( \mathbf x, t ) \Big | _ { \mathbf x = \mathbf x _ { * } ( t ) } \Phi ( t ), \ \Phi ( 0 ) = \mathbf I. \ ( 1 2 ) \\ \text{Due to the linearity and periodic nature of the} \\ \text{first-order variational nroblem. it can be shown}$$

Due to the linearity and periodic nature of the first-order variational problem, it can be shown that the growth or decay of any small perturbation is related to the eigenvalues of this principal fundamental matrix evaluated at t = T (also known as the monodromy matrix) Φ ( T ), called the Floquet (or characteristic) multipliers. Specifically, if any Floquet multiplier is greater than one in modulus, the solution is asymptotically unstable, and it is stable otherwise.

The monodromy matrix is thus a determining factor for stability. Unfortunately, there generally exists no explicit expression of this matrix from ∂ /∂ f x , but the key point to remember is that this Jacobian matrix plays a pivotal role for the asymptotic stability of a solution .

## 2.6 Illustration with a Duffing oscillator

For illustration, we will consider the damped, forced Duffing oscillator. This example is wellknown for the nonlinear dynamicist. However, it can exhibit rather complex dynamics associated with a wide variety of secondary resonances [42, 43].

The response q of the considered oscillator to an external excitation f is governed by a

generalization of Equation (5):

$$m \ddot { q } ( t ) + c \dot { q } ( t ) + k q ( t ) + k _ { 2 } q ^ { 2 } ( t ) + k _ { 3 } q ^ { 3 } ( t ) = f ( t ), \quad \text{ or } u \\. \. \. \. \. \. \. \ \text{ } \ \text{ } ( 1 3 ) \quad \text{ trial }$$

̸

where c is the damping coefficient of the oscillator, and we added a quadratic stiffness force with a coefficient k 2 that can be used to represent asymmetries in the system. When k 2 = 0, the oscillator characterized by Equation (13) is often called the Helmholtz-Duffing oscillator. Upon using the normalizations ¯ t = ω t 0 and ¯ = q √ k /kq 3 , Equation (13) can be put into the dimensionless form

$$\bar { q } ^ { \prime \prime } ( \bar { t } ) + 2 \zeta _ { 0 } \bar { q } ^ { \prime } ( \bar { t } ) + \bar { q } ( \bar { t } ) + \beta _ { 2 } \bar { q } ^ { 2 } ( t ) + \bar { q } ^ { 3 } ( \bar { t } ) = \bar { f } ( \bar { t } ), \ ( 1 4 ) \quad o b t.$$

with ζ 0 = c/ (2 √ km ), β 2 = k / 2 √ kk 3 , ¯ f = √ k /k f 3 3 , and a prime denotes derivation with respect to the dimensionless time ¯ . For the examt ples shown hereafter, we chose ζ 0 = 0 05 . and β 2 = 0 (hence we consider a Duffing oscillator), unless stated otherwise.

Figure 5 features a bifurcation diagram obtained by coupling a harmonic balance procedure with numerical continuation [11, 31]. This diagram represents the response amplitude of the Duffing oscillator under a harmonic forcing of amplitude ¯ ( ¯ ) f t = sin(¯ ¯ ) ωt (with a dimensionless frequency such that ω = ω ω 0 ¯), taking ¯ ω as bifurcation parameter. This bifurcation diagram is also often called the nonlinear frequency response (NFR). In addition to the well-known primary (1:1) resonance that bends rightward due to the hardening nonlinearity of the oscillator, a series of peaks appear throughout the bifurcation diagram and are associated to non-primary resonances [44]. Specifically, odd superharmonic ( m :1, m ∈ N 0 ) resonances appear on the main branch of the NFR at frequencies lower than that of the primary resonance, whereas the even superharmonic resonances bifurcate from it through symmetry-breaking bifurcations.

Figure 6 represents the NFR of the Duffing oscillator under a higher-amplitude harmonic forcing, namely ¯ ( ¯ ) f t = 3sin(¯ ¯ ). ωt In addition to the primary and superharmonic resonances, subharmonic (1: , l l ∈ N 0 ) and ultrasubharmonic ( ml : , m,l ∈ N 0 , gcd( m,l ) = 1) resonances appear, and they are isolated from the main NFR. We note that these isolated resonances were not obtained directly from numerical continuation of the main branch; stochastic procedures were used to find a set of initial conditions resulting in a subharmonic or ultrasubharmonic response. Once a suitable initial point was found, numerical continuation was used to trace out the complete isola. Thus, even in this simple oscillator, complex nonlinear features can appear.

## 2.7 Consequences for open-loop testing

We now illustrate how the complexity featured by nonlinear systems as deceptively simple as the Duffing oscillator impacts the experimental results obtained while testing a system.

During a test, a drive signal is generated by the experimenter and fed to the system. In vibration testing, this drive signal is fed to an actuator (generally through an amplifier) whose purpose is to excite the structure under test, either via a force or via base excitation. In both cases, the drive is seldom proportional to the forcing signal because the amplifier and the actuator have their dynamics [45], hence we distinguish the drive from the actual forcing. For a harmonic excitation, the two parameters that the user can control are the drive amplitude and frequency. The response of a nonlinear system is generally characterized for multiple values of these parameters. For instance, the drive amplitude may be fixed and the frequency varied. This variation can be continuous or discrete, representing swept-sine (SWS) or stepped-sine (STS) excitations, respectively. The set of points obtained from different frequencies can be interpreted as the NFR. Alternatively, one can fix the frequency and vary the drive amplitude and obtain a so-called S-curve [25].

We apply a SWS with logarithmically varying frequency to the Duffing oscillator, as shown in Figure 7. As the frequency increases, a smooth change in amplitude is observed, except in the vicinity of resonances. The most noticeable feature is the occurrence of a sudden jump-down at a frequency ¯ ω ≈ 7 2. . Conversely, by progressively decreasing the frequency, a jump-up occurs around ¯ ω ≈ 2 1. . The jump-up and jumpdown phenomena occur where the NFR folds , as seen in Figure 7. Indeed, SWS can be seen as an experimental equivalent to natural parameter continuation.

Fig. 5 : Bifurcation diagram of the Duffing oscillator with a dimensionless forcing amplitude of 1: main branch of the NFR ( ) and bifurcating even superharmonic resonances ( ). Full and dashed lines represent stable and unstable solutions, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-2c379eba-fbda-4d10-954c-ce93f5649850.jpg)

Fig. 6 : Bifurcation diagram of the Duffing oscillator with a dimensionless forcing amplitude of 3: main branch of the NFR ( ) bifurcating even superharmonic resonances ( ), isolated subharmonic resonances ( ) and isolated ultrasubharmonic resonances ( ). Full and dashed lines represent stable and unstable solutions, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-1f3857d2-0336-4c7c-af3a-9beba15b9180.jpg)

Closing up on the low-frequency region, we can observe most odd superharmonic resonances, whereas the even ones are sometimes partially missed due to the transient nature of the excitation. Sub- and ultrasubharmonic resonances are completely missed by the (simulated) testing procedure because of their isolated nature.

## 2.8 Challenges for experimental nonlinear dynamics

To summarize this section, we have seen that a combination of theoretical and numerical techniques can be used to obtain bifurcation diagrams that provide a great deal of information about the response of nonlinear systems. Thanks to the

Fig. 7 : Bifurcation diagram of the Duffing oscillator with a dimensionless forcing amplitude of 3: sweptup ( ) and swept-down ( ) excitation frequency. The results from numerical continuation are given for reference ( ), where full and dashed lines represent stable and unstable solutions, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-bec8c51e-5017-4675-a252-7a15b187039e.jpg)

maturity of these approaches, and in particular the numerical continuation procedure, stable and unstable parts of a bifurcation diagram can be obtained.

test. These methods are thus paradigm-shifting for nonlinear vibration testing and allow the experimenter to perform experimental continuation.

As for experiments, the classical open-loop testing procedures cannot offer the same picture. There are in particular specific challenges associated with nonlinear dynamics:

- 1. Unstable solutions are by nature impossible to follow with open-loop approaches.
- 2. Jumps are thus susceptible to occur when a stability change is encountered through bifurcations.
- 3. Several solutions are missed during a test, resulting in an incomplete bifurcation diagram.
- 4. Isolated solutions are challenging to reach with classical continuation approaches. In an open-loop setting, they can be obtained in some cases by using a stochastic approach [46], or leveraging jumps [47].

A novel class of testing methods, herein called control-based methods, has emerged slightly before the beginning of the last decade [23, 48]. The salient capabilities of these methods are that, through suitable automatic adjustments of the drive signal, they ( ) allow the experimenter i to prescribe given parameters of the system to extract complete bifurcation diagrams and ( ii ) can stabilize the equilibria of the system under

## 3 Feedback control: a brief introduction

Before delving into control-based methods, we review a few essential aspects about feedback control, which is a key enabler for experimental continuation.

## 3.1 Linear feedback control

Humans have always wanted to control their environment, i.e., to stir the systems that surround them toward desired states of operation. A natural way to perform this is to use feedback control, wherein the output of a system to be controlled (usually called plant) is sensed and used to adapt its input, eventually (and hopefully) reaching a desired output through the action of that input. When feedback is devoid of human interaction, one speaks of automatic control. We shall refer a system with and without feedback to as closedand open-loop systems, respectively.

Feedback is naturally used by most animals, but it can be surprising to realize that this concept was purposefully exploited rather late and formalized even later in human history [49]. Starting

from the twentieth century, the study of feedback control has known a great evolution, mostly motivated by technological needs and advances, leading to the ubiquitous presence of controlled systems around us today. For a historical perspective on this subject, we refer, e.g., to [50, 51].

Figure 8 presents the typical architecture of a linear control system. The plant, excited by the drive (or control signal) u and unknown disturbances d , is monitored through its response y (typically an element of x or a linear combination of them). This response is compared to a desired reference y ∗ (in the sequel, a subscript with a star indicates a target quantity) and fed to a controller, that automatically adapts the drive u . The controller can be designed to force y to follow the reference y ∗ as closely as possible, to reject the effect of the external disturbance d , or to do both these things at the same time.

In this linear time-invariant setting, the plant response can be expressed as the superposition of the effect of the drive u and the disturbance d . We use a Laplace transform formalism to express this as

$$Y ( p ) = H _ { y u } ( p ) U ( p ) + H _ { y d } ( p ) D ( p ), \quad ( 1 5 ) \quad \underset { 1 \colon n } { \longrightarrow }$$

where Y , U and D are the Laplace transforms of y , u and d , respectively, H yu and H yd are the transfer functions of the plant, and p is Laplace's variable. The drive is the controller output given by

$$U ( p ) = C ( p ) ( Y _ { * } ( p ) - Y ( p ) ), \quad ( 1 6 ) \quad \text{a1} _ { \text{st} }$$

where C is the controller's transfer function, and Y ∗ is the Laplace transform of y ∗ . Injecting Equation (16) into Equation (15) yields the closed-loop relation

$$\begin{smallmatrix} ( 1 + H _ { y u } ( p ) C ( p ) ) Y ( p ) \\ = H _ { y u } ( p ) C ( p ) Y _ { * } ( p ) + H _ { y d } ( p ) D ( p ) \end{smallmatrix} \quad ( 1 7 ) \quad \text{vat}$$

In general, the controller transfer function C is chosen to be large compared to the other transfer functions. Hence,

$$\lim _ { \substack { | C ( p ) | \to \infty \\ = \lim _ { | C ( p ) | \to \infty } } H _ { y u } ( p ) C ( p ) Y ( p ) \quad \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$$$

Hence, Y ( p ) ≈ Y ∗ ( p ), and y t ( ) ≈ y ∗ ( ) t if the controller has a large enough gain. It is thus theoretically possible to achieve perfect tracking between y ∗ and y . In practice, finite gains are chosen for several reasons. First, closed-loop stability may require finite gains (and in general, unknown, unmodeled dynamics will impose a finite gain). Second, sensor noise can pollute the sensed signal y and transmit fully to (and even be amplified by) the closed-loop system. Third, a large controller gain can require a large drive u , exceeding the actuator capabilities. There may also be other practical reasons depending on the application at hand.

## 3.2 Nonlinear control

Nonlinear control encompasses feedback control of systems where either the plant or the controller is nonlinear, or both. Naturally, its development was closely related to linear feedback control and nonlinear dynamics theories. The inherent difficulties associated with nonlinear dynamical systems unsurprisingly arise (and are arguably even more problematic) in nonlinear control systems, making them a contemporary active research field. For historical perspectives, we refer to [52-56]. There exists a variety of reference books on nonlinear control [41, 57, 58], with one specifically aimed toward nonlinear vibration control [59].

One of the challenges associated with nonlinear (and linear) control (except for adaptive and optimal nonlinear control, to a certain extent) stems from the need for a representative model of the plant for the controller design. Although substantial progress has been made in NSI, there still remains inherent challenges associated with this endeavor [20, 21], and the complexity associated with modeling a nonlinear system generally reverberates into nonlinear control. This motivated the use of model-less control approaches, which emerged in the context of chaos control.

In the early nineties, a seminal article written by Ott, Grebogi and Yorke proposed a new modelless method to control the behavior of a chaotic system [60]. This OGY algorithm sparked a lot of interest in the late twentieth and early twenty-first centuries; reviews on chaos control methods can be found in [61-63]. A notable contribution by Pyragas later proposed a simplified continuous delayed feedback law to achieve the same goal [64]. By

Fig. 8 : Schematics of a typical control system.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-d5335f48-d4f5-4acb-96f1-8b781770a8ef.jpg)

feeding back the difference between an observable and its time-delayed self, this control law was also shown to be able to stabilize the targeted periodic orbit. We note that choosing this delay to be equal to the period of the target makes the control noninvasive on that orbit, an important notion that we shall come back to and define later. Chaos control can be considered as the ancestor of the first control-based methods discussed in this work.

- 1. Continuation: unstable branches are connecting various stable branches, which might not (or not fully) be reached with open-loop system testing methods. For instance, in a HelmholtzDuffing oscillator, SWS tests may systematically miss the highest amplitude of the NFR, as illustrated in Figure 9. There can also be stable responses that can only be accessed through unstable branches (so-called islands of stability [67]).

## 4 Experimental continuation: the pioneers and first efforts

Control-based methods combine control and continuation to trace out a complete bifurcation diagram of a nonlinear system. This idea was first proposed for chaotic systems (or chaotic at least in some region of the investigated parameter range) by combining the OGY algorithm (or its variants) to a natural parameter continuation and was coined under the term 'tracking' [65, 66]. Methods combining continuous feedback and continuation to trace out bifurcation diagrams made up by the collection of periodic orbits of nonchaotic systems were introduced following the seminal works of Sieber and Krauskopf [23, 48] and are the focal point of the present work.

## 4.1 Control-based approaches: generalities

## 4.1.1 Why are closed-loop tests necessary?

Besides the academic thrill of tracing out complete experimental bifurcation diagrams, one may legitimately wonder what interest there is in attempting to follow unstable solutions in a general setting, given that they will be practically impossible to observe in the bare system. There are multiple reasons for this, and some of them are given below, the most important arguably being the first one.

- 2. Global dynamics: unstable equilibria of saddle type carry important information about the global dynamics. They lie on the separatrices of the system and provide valuable information about the basins of attraction of the system, and especially their boundaries. Loosely speaking, if a stable and an unstable equilibria are 'close', one would expect the basin of attraction of the former to be small. On a related note, due to the hyperbolic nature of the dynamics around such sets, saddles also play an important role in transient dynamics.
- 3. Qualitative changes through bifurcations: branching on a stable path through bifurcations might lead the system to a completely different type of behavior, whereas following an unstable path may yield a system response that is similar to its pre-bifurcation behavior, easing the analysis and interpretation of the results.
- 4. Data diversity: in an NSI or model updating framework, measuring unstable solutions offers more data points for building an accurate model of the system under test.
- 5. System interconnections: finally, when one considers nonlinear substructuring or nonlinear control, due to the feedback action enacted upon the system by its surroundings, unstable equilibria might become stable, or play a prominent role in the appearance of new stable equilibria, and thus become obviously practically relevant.

Fig. 9 : Bifurcation diagram of the HelmholtzDuffing oscillator with ζ 0 = 0 05, . β 2 = 1 9 and . a dimensionless forcing amplitude of 0.05: sweptup ( ) and swept-down ( ) excitation frequency. The results from numerical continuation are given for reference ( ), where full and dashed lines represent stable and unstable solutions, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-fbee3d96-f2e2-403b-bb7f-86ee6c862a30.jpg)

## 4.1.2 Common architecture and non-invasiveness

Figure 10 presents the general architecture of a control-based testing method (note the similarity with Figure 8). As in open-loop testing, the plant (e.g., in the vibration testing case, a structure under test, and a set of actuators and sensors) is excited by a drive signal u . The response of the plant is then analyzed through any of its observables (e.g., in the vibration testing case, a displacement q or a force f ) to obtain a representative signal y which is fed back to a controller. This controller compares this representative signal to a reference y ∗ , and adapts the drive signal to enforce as much as possible an equality between the reference and actual signals.

The target (or a set of targets) z ∗ is predetermined by the experimenter depending on the goals of the tests to be performed, and the selected experimental continuation method. In some methods, the continuation procedure will adapt y ∗ through a continuation procedure. In other methods (discussed in Section 5), one simply sets y ∗ = z ∗ , and the continuation block is trivial.

A feedback controller can and generally will completely modify the dynamics of the system.

Since the goal of experimental continuation is to test the plant, such a modification should be avoided. A requirement put on the control (in addition to its two main goals) is thus that it should be non-invasive , i.e., the equilibria of the plant in the closed-loop system should be identical to those of the open-loop system. A sufficient condition to enforce this is to subject the closedloop system to the same excitation as in open-loop conditions. Hence, the total input to the plant in closed-loop can be seen as the superposition of an open-loop excitation signal used to investigate the dynamics of the plant, and a stabilizing signal which vanishes once the system reaches the targeted open-loop equilibrium.

## 4.1.3 Stabilizing effect

As outlined in Section 4.1.2 and particularly in Figure 10, the main difference between an openloop test and a test with a control-based method is the presence of a feedback signal. We thus reconsider the results of Section 2.5 with this difference, staying at an abstract level for generality. A way to account for this feedback is to modify Equation (1) as

$$\dot { \mathbf x } ( t ) = \tilde { \mathbf f } ( \mathbf x ( t ), \mathbf u ( t ; \mathbf x ( t ) ), t ), \quad \ ( 1 9 )$$

where the right-hand side ˜ f now depends on u , the control signal. Among other possibilities, u may either be determined from a prescribed continuous-time evolution

$$\dot { \mathbf u } ( t ; \mathbf x ( t ) ) = g _ { 1 } ( \mathbf x ( t ), \mathbf u ( t ; \mathbf x ( t ) ), t ), \quad ( 2 0 )$$

be updated by a discrete-time law, or be a direct function of the states and time

$$\mathbf u ( t ; \mathbf x ( t ) ) = \mathbf g _ { 2 } ( \mathbf x ( t ), t ). \quad \quad ( 2 1 )$$

In our setting, the control is non-invasive if, at an equilibrium of the open-loop system x ∗ ( ), t

$$\mathbf u ( t ; \mathbf x _ { * } ( t ) ) = 0 \text{ \quad \ \ } ( 2 2 )$$

and

$$\tilde { f } ( x ( t ), 0, t ) = f ( x ( t ), t ), \quad \forall x, t. \quad ( 2 3 )$$

If these two conditions hold, we thus have

$$\mathfrak { n }. \quad \tilde { \mathfrak { f } } ( \mathbf x _ { * } ( t ), \mathbf u ( t ; \mathbf x _ { * } ( t ) ), t ) = \mathfrak { f } ( \mathbf x _ { * } ( t ), t ) = \dot { \mathbf x } _ { * } ( t ), \ ( 2 4 )$$

Fig. 10 : General schematics of a control-based method.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-ea2f8bae-9b96-4070-a496-d41a838db97b.jpg)

i.e., x ∗ is an equilibrium of both the open- and closed-loop systems. We note that in experiments, being exactly on the equilibrium x ∗ is not possible due to, e.g., noise. Equation (22) can therefore be seen as purely theoretical and can be replaced by the condition that the control signal should be smaller and smaller as the system approaches its equilibrium state. Instead of Equation (22), noninvasiveness can therefore be satisfied if the control is required to satisfy the Lipschitz continuity condition the different nature of the control signals used by the different approaches.

$$| | \mathbf u ( t, \mathbf x ( t ) ) | | \leq L | | \mathbf x ( t ) - \mathbf x _ { * } ( t ) | | \quad ( 2 5 ) \quad \text{th}$$

around x ∗ ( ), t were L is a finite positive constant and || · || is the Euclidean norm.

Following the same developments as in Section 2.5, we may study the stability of this equilibrium by introducing small perturbations, and linearizing the dynamics to reveal the Jacobian of the controlled system

$$\frac { \text{d} \tilde { f } } { \text{dx} } \Big | _ { x = x _ { * } ( t ) } = \Big ( \frac { \tilde { \partial } \tilde { f } } { \partial x } + \frac { \tilde { \partial } \tilde { f } } { \partial \bar { u } } \frac { \text{du} } { \partial \bar { x } } \Big ) \Big | _ { x = x _ { * } ( t ) } ( 2 6 ) \quad \text{res} \\ ( \text{where the total derivative is taken with fixed $t$} ), \quad \text{flic} \\ \text{and. using Equations} ( 2 2 ) \text{ and} ( 2 3 ). \text{ we find}$$

(where the total derivative is taken with fixed t ), and, using Equations (22) and (23), we find

$$\frac { d \tilde { f } } { d x } \Big | _ { x = x _ { * } ( t ) } = \Big ( \frac { \partial f } { \partial x } + \frac { \partial \tilde { f } } { \partial u } \frac { \partial u } { d x } \Big ) \Big | _ { x = x _ { * } ( t ) } \cdot \quad ( 2 7 ) \quad \text{$wh$} \\ \text{We thus observe that the Jacobian of the original} \quad \text{is $t$} \\ \text{system is modified by the effect of the control sig} \quad \text{the}$$

We thus observe that the Jacobian of the original system is modified by the effect of the control signal u . Consequently, according to the discussion in Section 2.5, the stability of x ∗ may be altered by the action of the controller. Unfortunately, going beyond this discussion while staying in a general framework may prove complicated, first because of the non-trivial relation between the Jacobian and the Floquet multipliers, but also because of

We may nevertheless illustrate this section with a specific example. Figure 11 presents the counterpart of Figure 6 when the Duffing oscillator is stabilized with a differential controller with a fixed reference signal ¯ q ∗ correctly ensuring the non-invasiveness condition, i.e., ¯ = ¯ f k d (¯ q ′ ∗ -¯ ), q ′ where ¯ f is sinusoidal and ¯ q ′ corresponds to an open-loop (stable or unstable) equilibrium solution to that forcing (as done in the controlbased continuation method, CBC, introduced in the next section). With the differential gain ¯ k d , the open-loop unstable branches can be partially (Figure 11a) or totally (Figure 11b) stabilized. When the gain is not large enough, the remaining unstable solution co-exists with other stable attractors for which the control is invasive (keeping ¯ q ∗ fixed to make the control non-invasive on the unstable solution). Interpretations of the stabilizing effects of the CBC in terms of folding and closed-loop NFR can be found in [28, 68] and [69], respectively.

The stabilizing and non-invasive properties of the control signal may intuitively seem like conflicting requirements to the newcomer: how can this signal stabilize anything if it is supposed to vanish? The key element of the answer lies in understanding that while it does indeed vanish when the system is exactly at the equilibrium, it does not when the system is near this equilibrium. To illustrate this, a simple and evocative example is that of an inverted pendulum on a cart, where the cart motion is controlled to ensure that the pendulum stays in its upright position. This control therefore stabilizes what would normally be an unstable equilibrium and acts against any perturbation. It is also non-invasive, because the cart stays still when the pendulum is upright in the

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-70ca2c60-5bfe-407d-b61c-da02ac97d232.jpg)

10

8

6

4

2

0

Amplitude (-)

0

1

2

3

4

Frequency (-)

(b)

Fig. 11 : Bifurcation diagram of the Duffing oscillator with a derivative controller with ¯ k d = 0 2 (a) . and ¯ k d = 2 (b), with a dimensionless forcing amplitude of 3: main branch of the NFR ( ) bifurcating even superharmonic resonances ( ), isolated subharmonic resonances ( ) and isolated ultrasubharmonic resonances ( ). Full and dashed lines represent (closed-loop) stable and unstable solutions, respectively.

(theoretical) absence of perturbations. Another illustrative example is given in [70, 71].

## 4.2 Control-based continuation

We now introduce the first experimental continuation method, namely CBC. Inspired by previous works [60, 64, 72], it was proposed by Sieber and Krauskopf [23, 48] in a very general setting by seeking an experimental equivalent to pseudoarclength continuation. The two key innovations of this work were ( i ) to use a PID-type controller 1 to stabilize the unstable responses of the system under test and ( ii ) to formulate a zero problem with measurable functions that indicate

1 We note that other types of controller can be used, as long as they satisfy the non-invasiveness requirement.

5

6

7

8

Fig. 12 : Schematics of the CBC-FD method; the red dashed frame indicates a controller block.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-5b8af13d-6514-4739-8d42-70c3bc25f66c.jpg)

a non-invasive control. Figure 12 schematizes this original version of the CBC method.

and Equation (28) becomes

This method is organized around the reference signal x ∗ , which is set so as to make the control non-invasive by prescribing the force applied on the structure f to a preset profile f ∗ (alternatively, there can be versions where the drive u is specified instead of f ), i.e.,

$$f ( t ; \mathbf w _ { * } ) - f _ { * } ( t ) \approx 0. \quad \quad ( 2 8 ) \quad \ n _ { b }$$

This non-invasiveness requirement is dealt with by using truncated Fourier series of the signals of interest,

$$x _ { * } ( t ) & \approx \mathbf h ( t ) \mathbf w _ { * }, \quad ( 2 9 ) \quad \mathbf i \\ \omega \mathbf m n \text{ vector with a rather than for prior}$$

where the column vector w ∗ gathers the Fourier coefficients of x ∗ , and h ( ) is the associated vector t of 2 h +1 harmonic functions

$$\text{w} \circ \, \vec { w } \,, \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \vec { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \, \hat { w } \\ \text{h} ^ { T } ( t ) = \begin{bmatrix} 1 \\ \sin ( \omega t ) \\ \vdots \\ \sin ( h \omega t ) \\ \cos ( \omega t ) \\ \vdots \\ \cos ( h \omega t ) \end{bmatrix}, \quad ( 3 0 ) \quad \hat { h } ^ { \ } a t _ { \ } t _ { \ } c _ { \ } c _ { \ } a \\ \text{upscript $T$ denotes the transposition} \quad n \\. Bv adjusting the Fourier coefficients w _ { \ } a$$

where superscript T denotes the transposition operator. By adjusting the Fourier coefficients w ∗ of the reference signal fed to the controller, the difference between the excitation f and a desired excitation profile f ∗ is sought to be as close to zero as possible. A well-posed problem can be obtained if the force profile is also discretized with Fourier series, i.e.,

$$f ( t ; w _ { * } ) \approx h ( t ) b ( w _ { * }, \omega ), \quad \ \ ( 3 1 )$$

$$b ( w _ { * }, \omega ) - b _ { * } = 0, \text{ \quad \ \ } ( 3 2 )$$

with b ∗ the desired force profile in a Fourier basis. For instance, if f ∗ ( ) t = sin( ωt ), b ∗ = e 2 , where e i is the i th canonical basis vector of R 2 h +1 . Equation (32) defines a zero problem with a nonlinear function b w ( ∗ , ω ), where the number of unknowns ( w ∗ and ω ) exceeds the number of equations by one. To solve this problem, a numerical pseudo-arclength continuation algorithm can be used. Since this algorithm relies on Newton-type methods to find the zeros of a residual function, it requires derivatives, which are not known in an experimental setup, but can be approximated using finite differences. Hence, in this work, this original version of the CBC is called CBC-FD. From the output of this continuation method, x ∗ ( ) can be reconstructed from t w ∗ using Equations (29) and (30).

In a vibration testing context, as the CBCFD method prescribes the forcing profile, so is its amplitude, and in general this method provides the NFR of a system, although variants can be conceptualized to provide other types of response curves. Figure 13 presents the main resonance of an energy harvester from [73], showing that this method is indeed able to obtain a full branch of an NFR.

## 4.2.1 CBC-FD applications

The first experimental validation of CBC-FD was with a vertically forced pendulum [24, 74]. An alternative way to compute the residual function using a time delay (similarly to [64]) to

## Figure subject to copyright. Please see the original publication.

Fig. 13 : Primary resonance of a nonlinear energy harvester obtained experimentally with the CBC-FD method [73]. Reprinted from [73], Copyright (2012), with permission from SAGE.

handle non-fundamental harmonics was also proposed in [48], and later used on energy harvesting systems [73, 75]. This extended time delay approach is simpler than the CBC-FD algorithm presented herein, but has less stabilization capabilities. The CBC-FD method was validated on a system with impacts [70] and on a geometrically nonlinear beam [76]. Interestingly, similar ideas have recently been introduced in quasistatic tests to avoid snapthrough and snapback phenomena [67, 77, 78], in real-time hybrid testing [79], and in testing procedures for microwave circuits [80].

## 5 Derivative-free experimental natural parameter continuation

Strategies to determine the stability of the found periodic orbit relying on temporarily switching the control off were investigated in [81], and the effect of noise was thoroughly studied in [82]. To the authors knowledge, Bureau et al [81] were also the first to measure a complete isola experimentally, as shown in Figure 14. The possibility to use B-splines (instead of Fourier decomposition) was recently explored by [83] for slow-fast systems. An in-depth dissertation on CBC-FD can be found in [71], and a software package, CONTINEX, is available in COCO [84].

## 4.2.2 Limitations

The CBC-FD method leverages pseudo-arclength continuation. Folding is almost guaranteed to be absent, which makes this method inherently more robust (continuation-wise) than traditional SWS tests. Unfortunately, this robustness comes at the expense of substantial complexity and longer tests due to the operations required by the method, such as the computation of a Jacobian through finite differences. This motivates the derivativefree methods that are treated in the next sections.

We now treat the set of state-of-the-art methods that can be classified as 'derivative-free', in the sense that they do not rely (explicitly) on a Jacobian of the residual function, and thus do not require finite differences and the likes (e.g. Broyden updates). We first present them in a unified framework, showing that they can be seen as derivative-free experimental natural parameter continuation methods (cf. Section 2.4.1), before reviewing them.

## 5.1 Continuation

Similarly to its numerical counterpart presented in Section 2.4, experimental continuation can be performed in two main variants. The first one, arclength-type continuation, is used in the CBCFD method as presented in Section 4.2. The second one is natural parameter continuation, for which a series of relevant parameters can be chosen, generally based on how representative they are of the excitation or response of the system.

We now specialize the discussion to nonlinear vibrations, given that most of the methods treated herein were developed in this area. The discussion can however be translated to other disciplines by considering relevant variables. To represent the state of a vibrating structure, four distinct variables are generally used:

- · f , an amplitude of the external forcing applied to the system,
- · a , the amplitude of some harmonic of the response at some point of the structure (typically, the same harmonic as the forcing of the collocated degree of freedom) and
- · ω , the frequency of this forcing,

## Figure subject to copyright. Please see the original publication.

Fig. 14 : First complete isola (1:3 subharmonic resonance of an impact oscillator) obtained experimentally thanks to two runs of the CBC-FD method (full lines). Data from SWS-up (+) and SWS-down ( ) ◦ is indicated for reference. Stable and unstable responses are indicated in black and gray, respectively. Reprinted from [81], Copyright (2014), with permission from Elsevier.

## · θ , the phase of this harmonic.

There are of course more variables than this, but these four are of particular interest in control-based methods. In the most basic set-ups described in Section 2.7, the operator controls f and ω , and a and θ are naturally determined by the evolution of the system (or, in a mathematical setting, they are the unknowns of the problem, whereas f and ω are given parameters). For linear systems, this evolution is often represented with a fixed f and a varying ω in Bode diagrams. Due to the invariance of the linear properties with respect to f , the results are usually normalized with the latter. For nonlinear systems however, variations of f retain their full importance, resulting in more than a mere scaling of the response.

'slices' of that manifold). We note that the terminology 'backbone curve' refers herein to a special case of a constantθ curve where θ = -π/ 2, which is often close to the amplitude maximum of the NFR [85].

Control-based methods change the paradigm of Section 2.7, in the sense that f and/or ω can be automatically adjusted to prescribe a and/or θ as parameters set by the experimenter. But in all cases, there are always two parameters among the four variables that can be set by the operator, while the two remaining variables adapt according to the dynamics of the system and can be used to represent its response.

In addition, one of the two chosen parameters is usually kept constant during a test, while the other is varied. The fixed parameter determines the type of response, whereas the changing parameter is most of the time used as the abscissa when plotting the response of the system. For illustration, Figure 15 represents all possible responses of a Duffing oscillator under tonal excitation as a two-dimensional manifold in the three-dimensional ( ω , f , a ) and ( ω , f , θ ) spaces, together with response curves obtained when one of these parameters is fixed (which can be seen as

With this view in mind, we can classify most of the experimental continuation methods (that shall be presented hereafter in this section) following Table 1. In this table, we can distinguish methods that require no, one or two control procedures. By control procedures, we mean, in a broad sense, any procedure that will adjust the input to the systems, using continuous feedback control (for SCBC and PLL) or discrete iterative methods (for SCBC and RCT). Methods simultaneously imposing f and ω require no control procedure, and are thus always open-loop. Most of the existing control-based methods require just one control procedure, with the exception of the SCBC method applied to the search of specific phase resonance [86], SCBC combined with PLL (CBPLL) [28], or PLL with amplitude control (PLL-A) [87].

Fixing a parameter while varying another one to parametrize the response makes the methods listed in Table 1 inherently equivalent to natural parameter continuation. As such, they are all susceptible to folding issues. This is well-known for the open-loop procedures, but should also be kept in mind when using the control-based methods listed in Table 1.

We can now present in more detail the main derivative-free experimental natural parameter continuation methods, namely, SCBC, PLL and RCT.

Fig. 15 : Response manifold of a Duffing oscillator under tonal excitation in the ( ω, f, a ) (a) and ( ω, f, θ ) (b) spaces, with constantf (NFR, , (c)) constantω (S-curve, , (d)), constantθ (backbone curve, , (e)) and constanta (harmonic force curve, , (f)) curves.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-0dc82282-cc8a-4296-a8e8-484ffd1cc7bd.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-c5bb92c0-e43d-4850-9a02-fefb60d99155.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-a7b59bb4-bfce-4569-a443-a65351765b1a.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-78806515-8958-4c69-ab5c-999deb4c6cbe.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-93961c7e-b4eb-443a-bade-fbdfce571d98.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-a7d5ef59-5cc4-4c09-a28b-5c7a97b95e7f.jpg)

## 5.2 Simplified control-based continuation

A substantial simplification of the original CBC algorithm was proposed by Barton et al [25] and is 18

Table 1 : Classification of the existing experimental natural parameter continuation methods based on the imposed (fixed and parametrization) parameters. Cells in white, light and dark gray require zero, one and two control procedures, respectively, and cells in black are impossible. The brackets indicate which type of curve is obtained.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-b0040db6-cbd6-479a-9715-7bb192abd727.jpg)

Fig. 16 : Schematics of the SCBC method; the red dashed frame indicates a controller block.

schematized in Figure 16. They noted that, since the PID controller is linear and time-invariant, a single-harmonic input x ∗ -x to the PID yields a single-harmonic drive u . They thus proposed to set the fundamental amplitude of the reference signal to a target value a ∗ ,

We refer to this simplified implementation as simplified CBC (SCBC).

$$x _ { * } ( t ) = a _ { * } \sin ( \omega t ) + \mathbf h _ { \mathbf n f } ( t ) \mathbf w _ { *, \mathbf n f } \quad ( 3 3 ) \quad c \mathbf u$$

and to set its other harmonic coefficients, w ∗ , nf , equal to those of x , thus satisfying the noninvasiveness requirement. Thus, if w are the Fourier coefficients of x (and w nf its nonfundamental Fourier coefficients), Equation (32) is changed into

$$w _ { *, n f } - w _ { n f } ( w _ { * }, \omega ) = 0. \quad \ \ ( 3 4 ) \quad \ \ a$$

They also simplified the continuation algorithm by noting that Equation (34) can be solved with a simple Picard iteration scheme, thereby removing the need for experimental finite differences.

The main reason for using a ∗ as a bifurcation parameter instead of ω is that it removes the aforementioned folding issue in most cases of interest, provided the controller's gain are properly chosen. Indeed, it can be observed in Figure 15d that Scurves can be parametrized with the amplitude a . It should be noted that the actual amplitude a will not be exactly equal to the target amplitude a ∗ , but will tend to it as the controller gains are increased. We finally note that, thanks to its broad theoretical framework, the SCBC can be performed with another bifurcation parameter than a ∗ , but the latter is preferred in the vast majority of cases.

A SCBC test is performed at a constant frequency ω , by progressively varying the value of the target amplitude a ∗ , obtaining so-called S-curves. Collecting these curves for different target amplitudes gives a surface (in the ( ω, f, a ∗ ) space) that can be sliced at a constant force level to obtain an NFR. This procedure essentially amounts to

fitting experimental data and performing interpolation on the fit. Figure 17, reproduced from [25], illustrates such a process.

## 5.2.1 SCBC applications

The SCBC method was applied to systems with friction [88] and beams [69], and advanced versions were developed by Renson et al, enabling the tracking of backbone curves [89], fold bifurcations [86] and internal resonances [90]. A simplified implementation using adaptive filters and circumventing the use of Picard iterations was proposed in [91], thereby speeding up the process. Interestingly, the SCBC approach can also be used in a purely numerical setting to solve challenging problems [92].

Noise can be handled using local Gaussian process regression [93]. From a NSI perspective, noise effects were shown to be smaller with CBC than with traditional testing methods [94], and successes using machine learning and/or Gaussian regression processes using SCBC results were demonstrated in [95, 96]. Procedures to determine the open-loop stability were developed based on data-driven linearization [97, 98]. A systematic method to choose the gains of the PD controller to guarantee the absence of folding was also developed for a Duffing oscillator [68]. Finally, we note that the SCBC found applications beyond structural mechanics, such as aeroelastic systems [99], biochemical processes [100], microelectromechanical systems [101] and pedestrian flow control [102, 103]. An in-depth dissertation on the SCBC method can be found in [28].

## 5.3 Phase-locked loops

Resonance conditions are oftentimes complicated to reach in nonlinear systems due to their closeness to unstable regions. This is what motivated Babitsky [104] to introduce velocity feedback, creating auto-resonant systems that automatically maintain themselves at resonance. This idea was later extended to continuation by Sokolov and Babitsky [105] who proposed the concept of phasecontrolled systems: the phase shift between the excitation and the response θ can be used as a parameter for the NFR. An experimental implementation of this concept was realized with a delayed acceleration feedback by Mojrzisch, Wallaschek and Bremer [106]. The unfolding of the

NFR also comes with feedback stabilization, offering the possibility to obtain its unstable parts.

Delayed displacement, velocity or acceleration feedback is very simple to realize in practice. It can be used, e.g., for (linear) experimental modal analysis [107]. However, it is generally invasive and stabilization can be tricky to achieve. The latter point has recently been discussed in [108] and treated by measuring the structural response at multiple locations. An alternative to this type of approach is to use phase-locked loops (PLLs).

Invented nearly one hundred years ago, PLLs are extensively used in electronics for signal synchronization purposes and are actually thought to be one of the largest engineering application of feedback control in the world [109]. Given a reference signal, PLLs are able to adjust the frequency and phase of another signal so as to lock the phase between these two signals to a prescribed value. Coupling PLLs with (electro-) mechanical systems was achieved even prior to the works of Babitsky for structural control as a solution to modal spillover [110, 111] and to track resonance in fatigue testing [112]. Phase parametrization of frequency response functions (FRFs) using PLLs was proposed in [113, 114], and the first time PLLs were used in an experimental continuation framework for nonlinear structural systems was the relatively recent work of Peter and Leine [26]. The main advantage that PLLs have over delayed feedback is their natural non-invasiveness.

A PLL test set-up consists of a structure under test, a phase detector, a PID controller, and the equivalent of a voltage-controlled oscillator, as depicted in Figure 18. The difference between a reference θ ∗ and θ , the phase shift between the response x and either the drive u or the force f , is fed to the PID controller, which automatically adapts the frequency of excitation ω to reach the desired phase. Indeed, if θ ∗ is constant, since ω is the output of the PID,

$$\omega ( t ) & = \omega ( 0 ) + k _ { \text{p} } ( \theta ( t ) - \theta _ { * } ) \\ & \quad + k _ { \text{i} } \int _ { 0 } ^ { t } ( \theta ( \tau ) - \theta _ { * } ) d \tau + k _ { \text{d} } \dot { \theta } ( t ), \\ \cdot \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

where k p , k i and k d are the proportional, integral and derivative gains of the controller, respectively. Differentiating Equation (35) with respect to time,

$$\dot { \omega } ( t ) = k _ { \text{p} } \dot { \theta } ( t ) + k _ { \text{i} } ( \theta ( t ) - \theta _ { \ast } ) + k _ { \text{d} } \ddot { \theta } ( t ), \quad ( 3 6 )$$

## Figure subject to copyright. Please see the original publication.

Fig. 17 : Manifold of the response of a nonlinear energy harvester obtained with the SCBC method. Black and gray curves represent S-curves and NFRs, respectively. The dark gray part of the manifold is openloop unstable. Reprinted figure with permission from [25]. Copyright (2013) by the American Physical Society.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-025611fb-a750-4628-835b-561a265bf701.jpg)

Fig. 18 : Schematics of the PLL method; the red dashed frame indicates a controller block.

̸

and assuming steady-state conditions ( ˙ ω = ˙ = 0) θ shows that, at steady state, θ t ( ) = θ ∗ as long as k i = 0. Once ω reaches a steady value, the control is naturally non-invasive and can stabilize unstable responses. By fixing the forcing amplitude f ∗ and varying the phase θ ∗ , one can obtain an NFR. Alternatively, by fixing the phase θ ∗ (typically to obtain phase resonance) and varying the force f ∗ , one can obtain a constant-phase curve (typically, a phase-resonance backbone curve).

an indication of technological maturity, and some of the associated results are plotted in Figure 19. In addition, an adaptation for multi-input-multioutput test case was sketched in [126] based on a centralized modal-based control scheme, and a framework for phase resonance testing of baseexcited systems was proposed in [127].

## 5.3.1 PLL applications

The main thrust for PLL application was (and still is) nonlinear modal analysis, given how easily it can be performed with this control-based method. Early works provided a theoretical study based on normal forms that was developed and validated in [115], and an NFR synthesis method based on the measured phase resonance curve was proposed in [116]. Since then, the PLL has been applied to a large variety of structures, including beams with friction [117, 118], impacts [87, 119], repelling magnets [120] or featuring extremely large vibration amplitude [121], a T-structure [122], plates with internal resonances [123, 124] and a turbine blade with friction [125]. This latter example is

Finally, we note that PLL was compared to velocity feedback in [108] and to CBC in [128]. All methods were shown to be capable of providing relevant results to the experimenter. In-depth dissertations on PLL can be found in [28, 129, 130].

## 5.4 Response-controlled stepped-sine testing

The most recent control-based method is the RCT approach, whose basic principles were first drafted in [131] and then elaborated in [132] together with a theoretical justification. Similarly to the SCBC and PLL, its underlying motivation came from the observation that the frequency-force curves are generally unfolded when the experiment is performed with a constant amplitude response a ∗ and varied excitation frequency ω . It can thus be seen as a dual of the SCBC method, as illustrated

## Figure subject to copyright. Please see the original publication.

Fig. 19 : NFRs (colored lines, featuring repeatability) and phase resonance backbone (black dashed line) obtained with a PLL of a turbine blade component with root friction. Reprinted figure with permission from [125]. Copyright (2020) by the American Society of Mechanical Engineers.

by Table 1. The general schematics of the RCT method are given in Figure 20.

SCBC, NFRs can then be extracted using interpolation of this HFS.

A key strength of RCT is that, due to numerous cases of modal tests where the structure is base-excited [133], conventional linear vibration testing software generally implement a possibility to perform the test at a constant response amplitude. The control principle used in Siemens Testlab (the commercial software used for all RCT demonstrations) is conceptually explained in [134], although a detailed formulation of the control scheme does not seem to be available. Simply put, it is an automatic trial-and-error procedure that adapts the drive amplitude at a given frequency to match the desired response amplitude. A test at a given target amplitude is started at a given frequency with a trial drive amplitude maintained constant over several periods of excitation, and the harmonic response of the structure is estimated, e.g., with a digital tracking filter. A frequency response function between the response and drive is estimated (or updated) based on the ratio of these amplitudes, and the next drive amplitude is calculated based on this update. The process is repeated until the response amplitude reaches its target (within tolerance bounds) or until a maximum number of corrections have been carried out. The test then carries on to the next frequency. With this method, the control is naturally non-invasive, because the drive is always tonal.

For a given target amplitude, the test yields a harmonic force curve (HFC). Collecting them for multiple target amplitudes, one obtains a harmonic force surface (HFS), and, similarly to

Another interesting feature of the RCT is the individual linear-like look of the frequency-force curves, making it possible to identify them by standard linear system identification tools providing modal characteristics. This quasi-linearity is supported by the single-nonlinear mode theory [135], which states that if a structure approximately vibrates in a tonal fashion according to one of its mode and if the amplitude of one degree of freedom is maintained constant, the amplitude at the other degrees of freedom will remain constant. This effectively 'freezes' the effect of geometrical nonlinearities throughout the structure, giving an intuitive understanding for the observed quasi-linear behavior. This approach also offers an alternative to harmonic force surface interpolation for obtaining NFRs by synthesizing them using identified modal properties.

The RCT method thus provides significant advantages, such as a simplified post-processing procedure leveraging quasi-linearity, and the possibility to use the method directly with commercial vibration measure and control apparatus. When the assumptions for quasi-linearity do not hold, the post-processing might become more complicated but the method still offers a valid control-based framework.

## 5.4.1 RCT applications

In spite of its relative youth, the RCT has successfully been applied to a variety of examples, including beam setups [27, 136, 137], a bolted missile [132], a stack piezoelectric actuator [138]

Fig. 20 : Schematics of the RCT method.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-a2eba578-b078-4647-8c2a-65429c606675.jpg)

and a highly-damped control fin actuation mechanism [139]. We refer to [140] for an in-depth dissertation on RCT. The RCT has not been experimentally compared to any other controlbased method yet, and such a comparison is carried out in this work.

based on two criteria. First, its non-fundamental harmonics are set by the non-invasiveness requirements similarly to the SCBC method. Second, its fundamental harmonic is used together with the excitation frequency to perform continuation.

## 6 A derivative-free arclength control-based continuation procedure

According to the discussion in Section 5.1, current control-based methods require a choice from the experimenters. They can use simple, derivativefree methods such as those presented in Table 1, but the test is then exposed to potential folding issues and the associated possible occurrence of a jump. Some examples will be given later in Section 7.4.1. Alternatively, they can opt for the more general CBC-FD method, which is however more complex, sensitive to noise and makes the test longer. The main reason for these two downsides is the need for experimental derivatives evaluated through finite differences or Broyden updates.

The purpose of this section is to propose a new method which is almost as generic as the CBC-FD method in the sense that it does not assume any parametrization a priori, but which also circumvents the need for experimental derivatives. This procedure is called arclength CBC (ACBC) herein, was originally presented in [28], is schematically depicted in Figure 21 and explained hereafter.

## 6.1 Non-invasive control

As in the CBC, a PID-type controller with a reference signal is used to stabilize the structure under test. This reference signal is determined

Since the method works along the principles of CBC, non-invasiveness is guaranteed as long as the reference and actual signals have identical non-fundamental harmonics. As in [91], the noninvasive character of the control signal is ensured by the use of adaptive filters. These filters are able to automatically perform a truncated Fourier decomposition of a signal x with h harmonics, i.e.,

$$x ( t ) \approx \mathbf h ( t ) \mathbf w, \text{ \quad \ \ } ( 3 7 )$$

where w approximate the 2 h + 1 Fourier coefficients of x .

The Widrow-Hoff least mean squares (LMS) algorithm can be used to update the Fourier coefficients w of a signal x through the simple discrete-time law

$$\begin{array} { c } \mathbf w ( ( n + 1 ) t _ { s } ) - \mathbf w ( n t _ { s } ) \\ = \mu t _ { s } ( x ( n t _ { s } ) - \mathbf h ( n t _ { s } ) \mathbf w ( n t _ { s } ) ) \mathbf h ^ { T } ( n t _ { s } ), \end{array} \quad ( 3 8 )$$

where µ is the filter gain and t s is the sampling time.

Assuming the Fourier coefficients converge, the control can thus be made non-invasive by defining the target signal as

$$x _ { * } ( t ) & \, = \, a _ { * } \sin ( \omega t ) + b _ { * } \cos ( \omega t ) \\ & \quad + \text{h} ( t ) \left ( \text{I} - e _ { 2 } e _ { 2 } ^ { T } - e _ { h + 2 } e _ { h + 2 } ^ { T } \right ) \text{w} \\ & \, = \, a _ { * } \sin ( \omega t ) + b _ { * } \cos ( \omega t ) + \text{h} ( t ) \text{P} _ { \text{nf} } \text{w}. \\ \text{P} _ { - e } \, \text{can be interpreted as a matrix whose action}$$

P nf can be interpreted as a matrix whose action retains only the non-fundamental harmonics so that, if Equations (37) and (39) hold, the signal x ∗ -x is mono-harmonic.

Fig. 21 : Schematics of the ACBC method; the red dashed frame indicates a controller block.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-43525321-e9fc-4d54-9e15-f50d090f7d94.jpg)

## 6.2 Simplified continuation problem

Having secured the non-invasiveness of the method, three parameters remain to be set by the method, namely a ∗ , b ∗ and ω . Without loss of generality, one of the Fourier coefficients can be set to zero, say b ∗ = 0. The two remaining unknowns are set by enforcing the fundamental forcing amplitude f to be equal to the target one f ∗ .

anomaly α of a point on the ellipse. The procedure can then be repeated using ellipses centered around each next point, thereby obtaining a set of points on the NFR in the same fashion as Crisfield's arclength continuation [38].

The method can probably be best understood by looking at the force imposed on the structure as a function of the frequency and target amplitude. On this representation of the HFS depicted in Figure 22a, the NFR can be seen as a level set of this surface, where the forcing amplitude f corresponds to the target f ∗ . In this example, we consider for simplicity infinitely large control gains, such that a = a ∗ . However, we note that the foregoing discussion also holds for finite gains provided they are large enough.

The new approach for experimental continuation leverages this, and is inspired from arclength continuation [37]. Because the NFR is onedimensional in the ( ω, a ∗ ) space, a sufficiently small ellipse with semi-axes ∆ ω and ∆ a ∗ and centered on a point of the branch intersects it twice, as depicted in Figure 22a. The intersections are characterized by f -f ∗ = 0, whereas every other point of the ellipse does not (in general) satisfy this equality. The continuation procedure can thus be formulated as follows: given a current point ( ω , a n n ) and an ellipse centered around it in the ( ω, a ∗ ) space, the next point of the NFR ( ω n +1 , a n +1 ) is found as the intersection of the ellipse and the level set where f = f ∗ , as seen in Figure 22b. To do so, the procedure sweeps along the ellipse by automatically varying the eccentric

A key advantage of this approach is that the harmonic forcing amplitude can easily be measured online, and there are several strategies to find the eccentric anomaly that solves the zero problem f -f ∗ = 0 without need for experimental derivatives:

- 1. A first simple possibility consists in sweeping the angle α t ( ) at a constant rate r α until f = f ∗ .
- 2. A continuous integral control law ˙ α = -k α ( f -f ∗ ) can be used. An advantage of this law is that it automatically makes the algorithm converge toward a point where f = f ∗ , and orients the curve so that the method does not turn back. Indeed, this law makes the point where f = f ∗ and ∂f/∂α &gt; 0 a stable equilibrium, whereas the other solution is unstable. As a consequence, this law converges to intersections which pierce the NFR outward as α is increased. However, since the sweep rate depends on the error f -f ∗ , the sweeping rate of the ellipse depends on the amplitude of the profile f -f ∗ , and is typically low and high away and close to resonances, respectively.
- 3. A discontinuous integral control law ˙ ( α t ) = -r α sign( f -f ∗ ) combines the advantages of the sweeping and integral control methods. However, the discontinuity of this control law can make the convergence toward a solution where f -f ∗ = 0 difficult.

The sweeping law can be chosen freely by the user depending on their preference, and is denoted by

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-03a6aa70-2a58-4ce2-8466-0acd376dbd7b.jpg)

## 6.3 The algorithm

f (-) Algorithm 1 summarizes the proposed derivativefree ACBC procedure. In this implementation, the prediction is made with the previous converged value of α , which amounts to implementing a secant predictor. ρ represents a relative tolerance, and σ &lt; 1 is a factor that makes the algorithm stop sweeping in more strict conditions than the convergence check in order to anticipate transient effects. If this preemptive procedure is not enough to meet the convergence conditions, one can adapt the sweeping law by making it slower.

One can also decide to measure the NFR in one direction or another (with globally increasing or globally decreasing forcing frequencies) by selecting a positive or negative value for k α or r α .

Fig. 22 : (a): HFS (levels of gray) and definition of the ellipse ( ) built around a given point ( · ) used for the ACBC method to obtain points ( ♦ ■ , ) on the NFR ( ) at ¯ = f ¯ f ∗ = 0 0875. (b): Profile of . the function ¯ f -¯ f ∗ along the ellipse, whose roots ( ♦ ■ , ) correspond to the intersection of the ellipse with the NFR.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-cfe5713e-ef22-410b-bef2-c248aec63bbd.jpg)

˙ α = s f ( -f ∗ ) in the sequel. The sweeping rate should be high for a fast continuation, but is nevertheless limited by the system's transients, as well as the adaptive filters speed to guarantee a non-invasive control.

Algorithm 1 Derivative-free arclength continuation for CBC experiment.

- 1: Input: ω , a 0 0 , ∆ ω, ∆ a ∗ selected by the user.

2:

n

0

3:

←

4:

5:

6:

7:

8:

9:

10:

11:

12:

13:

14:

15:

16:

17:

18:

19:

20:

21:

22:

23:

24:

k

if

α

&gt;

α

else

←

α

end if

←

loop

Prediction

(

ω, a

∗

)

(

ω

Wait

←

t

steady

Correction while

f

˙

|

-

α

= (

s f

(

f

r

α

&gt;

0

then

n

+∆ cos(

ω

α , a

)

n

+∆

a

∗

seconds for steady-state

∗

&gt; σρf

∗

do

-

ω, a

∗

)

|

f

∗

)

(

end while

←

Wait

t

steady seconds for steady-state

Convergence check if

f

f

∗

ρf

∗

then

(

|

-

|

≤

ω

n

+1

n

, a

←

n

else

Adapt end if

25:

end loop

The proposed algorithm is equivalent to the CBC-FD in terms of generality of the continuation approach, as it does not make any assumption on the unfolding property of a parameter. A

n

+1

+1

←

s

(optional)

)

(

ω, a

∗

ω

n

+∆

ω

cos(

α , a

)

n

+∆

a

∗

0

π

0

or

)

sin(

α

sin(

α

))

))

key advantage of this approach is its derivativefree nature, brought by the combined use of the adaptive filters and the arclength continuation procedure in a two-dimensional space.

## 6.4 Illustration with a Helmholtz-Duffing oscillator

The Helmholtz-Duffing example in Figure 9 is revisited with a numerical simulation of the ACBC algorithm in Simulink with ¯ k d = 1, k α = 1, h = 15, µ = 0 0025, . ∆¯ a ∗ = ∆¯ = 0 05 and ω . ¯ t steady = 50. Clearly, the ACBC algorithm is able to capture the full NFR in this example, including its open-loop unstable parts.

## 7 Illustration with an electronic Duffing oscillator

Having reviewed the different control-based approaches for experimental continuation, we illustrate them with two experimental examples.

A first demonstration is made with an analog electronic circuit that replicates the behavior of a Duffing oscillator [141, 142]. This circuit is described in more details in the supplementary materials of this article and in [143]. The main advantages of using such an experimental demonstrator are its near perfect repeatability and actuation (i.e., the drive and the force are almost exactly equal), which allows us to test the inherent features of control-based methods. Yet, this setup is subject to practical experimental challenges, such as noise, unknown non-ideal behavior, and real-time constraints.

## 7.1 Experimental setup

A picture of the experimental apparatus, comprising the electronic Duffing oscillator and an acquisition and control system, is given in Figure 24a and schematized in Figure 24b. The oscillator takes a forcing signal V in as input, and outputs a displacement signal V x and a velocity-proportional signal V v , with V v ≈ -10 -3 ˙ V x . Compared to Figure 10, the acquisition and control system encompasses the continuation, controller and signal analyzer blocks. For all methods except the RCT, this system was a MicroLabBox from dSPACE (as in Figure 24a), whereas the RCT was performed with a Scadas Mobile from Siemens.

Fig. 23 : Bifurcation diagram of the HelmholtzDuffing oscillator with ζ 0 = 0 05, . β 2 = 1 9 and a . dimensionless forcing amplitude of 0.05 obtained with the ACBC: input parameter space (a) and NFR (b). The full ( ω, a ∗ ) curve is shown ( ) together with the results at convergence ( ◦ ). The results from numerical continuation are also given for reference ( ), where full and dashed lines represent open-loop stable and unstable solutions, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-e7a90bec-7169-4640-b698-6fc60ed2abef.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-774ca954-8eea-48f1-a372-9010f0f85e60.jpg)

The parameters c , k and k 3 (cf. Equation (13)) of the electronic Duffing oscillator can be tuned with knobs (potentiometers), and their theoretical

Fig. 24 : Picture (a) and schematics (b) of the experimental setup with an electronic Duffing oscillator.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-d103e472-435a-4771-9f3b-95477447204a.jpg)

values can be computed based on the electrical parameters of the board (see supplementary materials). They are reported in Table 2. They were compared to an identification with the nonlinear frequency subspace identification (FNSI) method [144] using data from a swept-up sine test. A close agreement between the theoretical and identified values can be observed, with a discrepancy of up to 8% for the damping coefficient. These parameters correspond to a natural frequency ω 0 = 131 5 . rad/s (20.9 Hz) and a damping ratio ζ 0 = 0 005. . This damping ratio can be considered rather small although it is still representative for a number of metallic structures. It was initially chosen to result in large amplitudes for moderate forcing, and strong nonlinearity, but it also revealed some limitations of various methods that seemed interesting to report herein.

## 7.2 Implementation details

All methods (except for the RCT) were programmed in a fully online way with dSPACE

Table 2 : Theoretical and identified parameters of the electronic Duffing oscillator.

| Parameter     | Theoretical value   | Identified value   |
|---------------|---------------------|--------------------|
| m ( s 2 )     | 10 - 4              | 1 . 0461 × 10 - 4  |
| c ( s )       | 1.3 × 10 - 4        | 1.4231 × 10 - 4    |
| k (-)         | 1.923               | 1.8099             |
| k 3 ( V - 2 ) | 0.9887              | 1.0077             |

RTI for Simulink. Almost all methods require a Fourier decomposition. This decomposition was performed online by adaptive filters, which were shown to be quick, effective and accurate [28, 87, 91]. For a tonal signal at frequency ω , an expression for an optimal gain (minimizing the settling time of the Fourier decomposition) was derived in [28] 2 . Here, based on a linearization of this expression with respect to the (small) sampling time, a frequency-dependent gain was used

$$\mu = \omega \bar { \mu }, \text{ \quad \ \ } \text{(4)}$$

where the optimal dimensionless gain based on this linearization is ¯ = 2. For a signal with multiµ ple harmonics, this value was observed to lead to a slightly too aggressive filter, and a smaller value was usually selected. Another procedure to tune µ was recently proposed in [145].

Some iterative methods propose specific procedures to detect if the system is at steady-state. In this work, we simply waited a certain time t steady before subjecting the system to another set of parameters.

All CBC methods were implemented with the same differential controller. It should be noted that the velocity signal was used instead of the displacement one, thereby avoiding the use of an approximate derivation with a high-pass filter (as can commonly be done in practical PID implementations).

Stepped methods (CBC-FD, SCBC, ACBC and RCT) require the definition of steps ∆ ω and ∆ a ∗ , which are based on the desired resolution of the NFR. For CBC methods, since the feedback is made with a velocity-type signal, the steps were defined in terms of this signal, ∆ | V v 1 | for simplicity, where | · | k denotes the k th harmonic amplitude of a signal. For the RCT, a (nearly) equivalent step in terms of displacement

2 Note that the symbol µ used herein corresponds to µt s in [28].

∆ | V x 1 | was defined instead. The two are related by ∆ | V v 1 | = 10 -3 ω ∆ | V x 1 | .

Most of the parameters were tuned on a trial and error basis. A fine-tuning toward optimal performance was not sought given the very large number of parameters and their variety. The results given therein are thus not to be interpreted as the best each method can yield, but rather as an order of idea of how they compare.

In the sequel, most of the results will be presented in terms of the total amplitude of the response over a period, being represented by the operator | · | (without subscript), i.e.,

$$| x | = \max _ { t \in [ 0, 2 \pi / \omega ] } \L a b s ( x ( t ) ), \quad \ \ ( 4 1 )$$

where abs is the absolute value operator.

## 7.3 Primary resonance

The primary resonance of the oscillator was first tested. The time for the two methods where parameters are swept, namely SWS and PLL, were chosen based on the fastest stepped method (in this case, ACBC). All parameters are given in the supplementary materials. We focus on measuring an NFR at a specific excitation f = 0 03 V . between 50 and 200 rad/s.

## 7.3.1 Comparison of the NFR obtained with different measurement approaches

To start this analysis off, Figure 25 shows the velocity signal NFR directly obtained from methods that ran without issues, namely, the SWS, PLL and ACBC approaches. The agreement between the two latter is excellent, even at the resonance peak. Some discrepancies are observable with the SWS data: as usual, transient effects are the cause for a slightly lower peak amplitude. We also note that the SWS data includes higher harmonics as well. The PLL test experiences large noise at low frequency, which is due to the 2:1 superharmonic resonance of the oscillator. While it is too small to be noticeable in the amplitude of the NFR, Figure 25b shows that the phase of the response evolves in a non-monotonous way at low frequencies, due to this superharmonic resonance.

Once the phase becomes unique again (in particular near the primary resonance), this is no longer an issue.

Fig. 25 : Velocity NFR of the electronic Duffing oscillator obtained with the SWS-up ( ) SWSdown ( ), PLL ( ) and ACBC ( ) methods ((a): amplitude, (b): phase).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-e02ed9ad-0be9-4f92-877a-2169cc9d1427.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-3055477e-99e1-4af5-affb-2525d7bd890e.jpg)

## 7.3.2 Challenges for the pseudo-arclength CBC method

The results obtained with the CBC-FD method are compared to those of the ACBC method in

Figure 26. Convergence issues were encountered with the former method in the vicinity of the resonance peak, requiring to restart the method from the high-frequency limit (200 rad/s) to approach the peak from both sides.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-58a10ac3-173d-4a43-b115-49a75505d40a.jpg)

Fig. 26 : Velocity NFR of the electronic Duffing oscillator obtained with the ACBC ( ) and CBCFD ( , □ : starting from 50 rad/s, △ : starting from 200 rad/s) methods.

The difficulties encountered with CBC-FD are most likely due to two main causes. First, owing to the small damping of the oscillator, the NFR curve experiences a sharp turn near the fold bifurcation. As a consequence, the pseudo-arclength method can struggle to converge to the curve during a correction phase. This issue could be solved by using an arclength continuation scheme, or taking smaller steps. The second issue is linked to the use of finite differences, and again, is mostly problematic in the vicinity of the resonance peak. Indeed, close to this peak, the residual function varies extremely rapidly and has a large curvature (as will also be illustrated by Figure 27a). Similar issues were encountered in [75, 78, 79]. Smaller finite difference steps could have been taken, but this would shift the issues to zones away from resonance, where the residual function is much less sensitive to variations in the input, and where large enough finite difference steps need to be taken to overcome noise. This issue could be solved using an adaptive finite difference step by, e.g., prescribing an ideal change in the residual function (suitably larger than the noise level).

## 7.3.3 Resolution requirements for the methods using interpolation

As explained in the respective works where the methods were proposed, both SCBC and RCT obtain a collection of response curves that can be combined into a response surface. NFR curves can be extracted from this surface by slicing it with a plane characterized by | f | = f ∗ , which essentially amounts to performing interpolation on experimental data. An off-the-shelf solution used in this work was to use the griddedInterpolant function in Matlab to fit the measured force and amplitude as 2D spline functions of the frequency and target amplitude, as ˆ ( f ω, a ∗ ) and ˆ( a ω, a ∗ ), respectively (where a hat denotes an interpolated function). The curves satisfying ˆ ( f ω, a ∗ ) -f ∗ = 0 were then found using the fimplicit function. Figure 27a shows such curves obtained for the SCBC. Points of these curves were finally fed to the fitted function ˆ( a ω, a ∗ ) to obtain the NFR curve.

The results of the SCBC and RCT obtained with this method are compared to those of the ACBC in Figure 27b. In this case, we analyze the fundamental harmonic amplitude rather than the total (multi-harmonic) one, because the latter is more complicated to extract with the results from TestLab with the RCT. These two amplitudes do not differ much anyway. While the interpolated NFRs are close to the one from ACBC away from the resonance peak, there is a clear discrepancy between the two types of results in the vicinity of the peak. This issue is mainly due to the interpolation, and to an insufficient resolution of the ( ω, a ∗ ) test grid: near resonance, the function ˆ ( f ω, a ∗ ) varies rapidly and requires a very fine resolution for the interpolation to be correct. Similar issues were encountered in [69] for the detection of fold bifurcations from the response surface. This can be solved using either a finer grid (but this requires a longer test time), or using an adaptive sampling strategy instead of a regular grid. Figure 28 supports this hypothesis by showing the results obtained by narrowing the test frequencies from 120 to 180 rad/s with steps of 2 rad/s, and amplitude steps of 0.05 V, effectively doubling the resolution in both directions. Clearly, the

Fig. 27 : Contour plot of the measured force for the SCBC test with the contour ˆ f = 0 03 . V highlighted ( ) (a), and velocity NFR of the electronic Duffing oscillator obtained with the ACBC method ( ) and estimated with the SCBC ( ) and RCT ( ) methods (b).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-0df461b7-5268-446f-8e74-23d7e13e1c35.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-50a78e9b-0ea3-41df-9fc4-f9047a6f6ec7.jpg)

agreement between both methods is improved, at the expense of a longer test time (increasing by a factor 4). The interpolation with the SCBC still inadequately represents a range of about 10 rad/s near the resonance peak, though.

In addition to this issue, we noted that the SCBC method could encounter convergence struggles, as will be discussed in Section 7.3.4. As for

Fig. 28 : Velocity NFR of the electronic Duffing oscillator obtained with the ACBC method ( ) and estimated with the SCBC method with higher frequency and amplitude resolution ( ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-f9cf822a-2fbb-483f-96a9-3ed7e31153d6.jpg)

the RCT, another issue related to stability was encountered, which is why the curve in Figure 27b appears incomplete. This issue is discussed in Section 7.3.5.

## 7.3.4 Invasiveness of the SCBC in the vicinity of superharmonic resonances

The SCBC method was implemented with both the original scheme based on Picard iterations [25], as well as a scheme that exploits the output of an adaptive filter [91] to make the control non-invasive. The former was observed to be much slower than the latter at some measurement points, forcing a raise of the tolerance to 5% for the Picard scheme (instead of 1% for the adaptive scheme) to allow for convergence. As shall be shown hereafter, this occurs in the vicinity of strong superharmonic resonances. Away from them, both methods performed similarly. The difference between the two methods in terms of results presented in Figure 27b were negligible, which is why only the results of the method based on adaptive filters were shown.

Figure 29a compares the results of the two methods in terms of the response surface. Differences are observable at low frequencies, while the agreement is almost perfect for the rest of the surface (not shown here). Figure 29b shows

the relative error, given by the invasiveness (the norm of the left-hand side of Equation (34) divided by a ∗ ) made by the Picard scheme, featuring the largest error in the zones of largest discrepancies. The distribution of these zones clearly follows a trend, making superharmonic resonances a potential suspect for this effect. The invasiveness of the SCBC with adaptive filters was also checked a posteriori, and it was also observed that the control shows the largest degree of invasiveness in the same regions, in spite of the expected action of the adaptive filter. This was not detected during the test, because the check for non-invasiveness was performed based on the error made by the signal synthesized by the adaptive filter, x t ( ) -h ( ) t w ( ). t

To ascertain that superharmonic resonances can be an issue, the 2:1 and 3:1 resonances were tracked with the PLL, by locking the phase of the second and third harmonic to -π and -π/ 2, respectively 3 , and varying the excitation amplitude (implementation details are given in the supplementary materials). Overlaying the obtained results 4 with those of the SCBC test with Picard iterations in Figure 29b clearly highlights that the issue indeed occurs close to the superharmonic resonances. The non-fundamental harmonic coefficients are set with a Picard iteration scheme, and a failure of this scheme in the vicinity of a resonance is coherent with numerical experience [146].

## 7.3.5 Instabilities with the RCT approach

We now turn to the incomplete aspect of the NFR curve obtained with the RCT, which is due to stability issues during the test. Figure 30a plots the actual amplitude of the displacement as a function of the target amplitude and forcing frequency. Zones with large errors are clearly apparent. The points with a relative error larger than 20% were rejected when performing the interpolation, explaining the incomplete aspect of the NFR curve obtained with the RCT in Figure 27b. These zones also appear to organize into a particular pattern reminiscent of the zone of instability

3 The resonant phase lag for the 2:1 resonance differs from that observed in a perfect Duffing oscillator due to the inevitable presence of asymmetries, making it more similar to a Helmholtz-Duffing oscillator in this regard.

4 The forcing amplitude was converted to an equivalent imposed amplitude for the reference signal via a ∗ = f/k d + . a

Fig. 29 : Comparison of the response surfaces obtained with the SCBC method with a Picard scheme ( ) and adaptive filters ( ) (a), and contour plot of the invasiveness of the SCBC method with a Picard scheme with the locus of 2:1 ( ) and 3:1 ( ) superharmonic resonances (b).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-33f65651-84a8-4934-8ead-28428779193b.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-cfe50d88-74b7-4da3-b9df-2d681d4cf29a.jpg)

of the open-loop oscillator, delineated by the loci of fold bifurcations.

To verify this hypothesis, the experimental fold bifurcation tracking procedure proposed by [86] was implemented (implementation details are given in the supplementary materials), and the resulting loci of fold bifurcations are overlaid in Figure 30b. The near complete depletion

Fig. 30 : Contour plot of the actual harmonic (displacement) amplitude of the response of the RCT method (where × indicate points with a relative error larger than 20%) (a) and scatter plot of the frequency vs actual (displacement) amplitude (color scales represent the relative error between the target and actual amplitude) with the locus of fold bifurcations ( ) (b).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-27b3b099-ba7f-404f-9f54-130994d3a02b.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-b58a6be7-5411-4e4f-971d-b9548afe294b.jpg)

of points between the two fold bifurcation curves confirms that the issue is very likely to come from the instability of the oscillator. It thus appears that the RCT was not able to stabilize unstable solutions in this case, and the oscillatory look of the amplitude observed in the large-error region in Figure 30a corresponds to alternating convergence of the response to the upper and lower stable branches. The presence of many unstable points makes the controller struggle and systematically reach its maximum number of iterations, explaining why the RCT took about an order of magnitude longer than the other methods (cf. Table 3).

Some points lay inside the unstable zone delineated in Figure 30b, and an analysis of the time series revealed that they are associated with transient phenomena. The resulting amplitude thus cannot be considered as a steady-state quantity. For illustration, Figure 31 depicts the time series associated with the point ( ω, a ∗ ) = (170 0 9) for , . which the amplitude indicated by the TestLab software results is | V x 1 | = 0 9135, hence seemingly . correct enough. The time series of the amplitude estimated by TestLab was unfortunately unavailable and was retrieved through time-frequency analysis (using the short-time Fourier transform from the spectrum function in Matlab). Its estimation is overlaid together with the target amplitude. It appears that the amplitude of the considered point, indicated by the middle step in Figure 31, is calculated during transients (right before the end of that step) and is thus not representative of the steady-state response of the oscillator. Assessing this in a systematic way is unfortunately difficult given the little amount of information that can be extracted from TestLab results.

The accumulation of points at the left of the leftmost branch, and the lack of points at the right of the rightmost branch in Figure 30b, which are both stable regions, was first thought to be due to the stepped-up frequency during the test. However, a stepped-down frequency surprisingly showed the same trend. The reason for this puzzling trend was not investigated further.

Finally, we note that there are points in the upper left corner of Figure 30a which also feature a large error. Similarly to the case exposed in Figure 29b, this is most likely due to superharmonic resonances, where folding is also expected to occur at high amplitudes.

Because the RCT works with a black-box commercial software, there is no potential implementation issue. Different parameters were tried out without success, as a zone of instability always

Fig. 31 : Amplitude of the displacement ( ), estimated harmonic amplitude ( ) and target amplitude ( ) during the RCT at a ∗ = 0 9 V and . around ω = 170 rad/s.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-8da9ca07-30b7-4a1b-998c-9a400a515a85.jpg)

seemed to appear in the same range (more details are given in the supplementary materials). The authors believe that the main difference between this work and the original [27, 132] is the presence of an electrodynamic shaker. This is known to stabilize unstable responses in some cases [147], which might be the case therein. This issue would deserve further investigation but was deemed beyond the scope of this work. Nevertheless, if this explanation is correct, the extraction of the folded part of the NFR of the structure alone is more related to continuation aspects than stabilization ones, since the RCT itself performs continuation on curves that do not fold but does not appear to inherently change the stability of the system under test (composed of the structure and the shaker).

## 7.3.6 Testing times

We conclude this first series of tests by providing the times needed to perform each experiment in Table 3. In terms of testing times, the fact that we focused on one NFR puts direct NFR measurement techniques, namely SWS, PLL, CBC-FD and ACBC, at an advantage compared to the other methods. It should thus be kept in mind when looking at the total time taken by each methods that SCBC and RCT allow for the extraction of multiple NFRs in one measurement run, and a more meaningful measure of the time efficiency of a method can be given by the average time taken per measurement point.

## 7.4 Superharmonic resonance

We now turn to the 3:1 superharmonic resonance of the Duffing oscillator. This type of resonance is seldom studied in control-based methods, but actually challenges all of them.

## 7.4.1 Numerical illustration

Before moving forward with the experiments, a numerical study has been made around the 3:1 superharmonic resonance of the Duffing oscillator. Variants of a harmonic balance method coupled with numerical continuation [11] have been used to compute constant forcing amplitude, forcing frequency and response amplitude curves with the theoretical parameters in Table 2. They are displayed in Figure 32.

Clearly, all methods based on natural parameter continuation, namely SWS, SCBC, PLL and RCT will fail, because folding is present in all cases. Similar observations were made in [28]. As for the ACBC approach, success is a priori uncertain given that the first harmonic amplitude curve features a loop that does not form a contour level.

A solution nevertheless exists for the PLL, as exploited in [148], which consists in replacing the phase lag of the first harmonic by that of the third harmonic. As Figure 33b evidences, such a parametrization does not exhibit any fold in the vicinity of the 3:1 resonance, making it a viable approach.

Finally, we noted experimental difficulties with the CBC-FD method as well. In the vicinity of the 3:1 resonance, non-invasiveness was found to be crucial for the accuracy of the results and for the success of the method itself. To investigate the impact of invasiveness, a series of numerical simulations were performed by imposing a non-tonal forcing to the oscillator. To the nominal fundamental harmonic forcing was added a forcing on the third harmonic of 1% relative amplitude and fixed phase. Performing the continuation for multiple values of the phase, we obtained an idea of locus of responses if a 1% relative tolerance is allowed on the residual, as depicted in Figure 34. In spite of the small tolerance, the loop associated with the resonance can vary by up to about 50%

Table 3 : Times taken by each test for the primary resonance.

| Method            | Total time (s)   | Time per point (s)   |
|-------------------|------------------|----------------------|
| SWS (up and down) | 240              | /                    |
| CBC-FD            | 396              | 4.96                 |
| SCBC (Picard)     | 453              | 0.62                 |
| SCBC (adaptive)   | 397              | 0.55                 |
| ACBC              | 240              | 3.19                 |
| PLL               | 240              | /                    |
| RCT               | 7 847            | 10.13                |

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-df77ef82-1f63-4d7d-86fc-17a1708470ad.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-c5bea966-2d37-48a6-91b3-aac4a4a4db60.jpg)

0.05

Fig. 32 : Numerical simulation of the response of the first harmonic of the Duffing oscillator with f = 1 V ((a): frequency parametrization, (b): phase parametrization), with ω = 55 5 rad/s (c), and . | V x 1 | = 0 5 . V (d).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-7b77885d-4045-4815-9f39-8f728927f05d.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-ad1de8ed-1836-4345-aa37-d586e841e45a.jpg)

Fig. 33 : Numerical simulation of the response of the third harmonic of the Duffing oscillator with f = 1 V ((a): amplitude, (b): phase).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-da288251-3223-4c11-99d4-320e02275e85.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-239ab9a6-1fc4-454e-af14-bde202fee57f.jpg)

in size, stressing the importance of a non-invasive control.

Fig. 34 : Numerical simulation of the response of the first harmonic of the Duffing oscillator with nominal mono-harmonic forcing with f = 1 V ( ) and locus of responses with a 1% additive forcing amplitude and arbitrary phase on the third harmonic ( ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-c3b7c311-819a-4fd5-a4a9-bc2b9342a2da.jpg)

Numerical simulations of the Duffing oscillator tested with the CBC-FD uncovered a bottleneck linked to this issue. If too large a relative tolerance was chosen (such as 1%), the continuation procedure was progressing in an erratic fashion inside a locus similar to that depicted in Figure 34. The resulting irregular curve posed challenges for the continuation by making inappropriate predictions which could not be corrected. Alternatively, choosing a small tolerance made the curve less erratic. However, inaccuracies on the residual (coming from the Fourier decomposition) were observed to increase the vicinity of the resonance, eventually exceeding the tolerance threshold, rendering the convergence of the method impossible.

## 7.4.2 Experimental testing

In view of the issues outlined in Sections 7.3.4, 7.3.5 and 7.4.1, the SCBC and RCT methods were not deemed suitable for the 3:1 superharmonic resonance. A SWS was performed in spite of its deficiencies to obtain a reference method, together with the CBC-FD, ACBC and

PLL (with the phase of the third harmonic), with a forcing amplitude of 1 V. All parameters are given in the supplementary materials.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-da9279d4-c0e6-4863-9d45-511ba212e845.jpg)

(a)

Fig. 35 : (a): Velocity NFR around the 3:1 superharmonic resonance of the electronic Duffing oscillator obtained with the PLL ( ), SWS-up ( ) and SWS-down ( ). (b): fundamental harmonic of the velocity signal around the 3:1 superharmonic resonance of the electronic Duffing oscillator obtained with the PLL ( ), ACBC ( ) and CBC-FD ( , □ : starting from 45 rad/s, △ : starting from 55 rad/s) methods.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-2974eab2-bed5-44ee-a908-b4ddaa60bea0.jpg)

Figure 35a shows the results of the most successful method, namely the PLL, for the 3:1 resonance, together with swept-sine results for reference. The PLL is able to robustly go through the unstable branches and trace the complete NFR, and it was the only method able to do so thanks to the use of the phase of the third harmonic.

The NFR follows a loop, which is easier to see in Figure 35b where only the first harmonic amplitude is depicted. The results of the CBCFD and ACBC are also shown therein. For the former, convergence of the Newton-Raphson procedure did not occur at some point. This lack of convergence could be due to the increased noise level observed around the 3:1 resonance, making the finite-difference estimation of the Jacobian inaccurate. As for the ACBC, the control was seen to be invasive, signalling that the adaptive filters were failing to fulfill their purpose, as in Section 7.3.4. The estimated NFR consequently significantly diverges from what is estimated to be its true shape, similarly to the issue raised with Figure 34. In addition, the fact that the NFR loops around the 3:1 resonance (Figure 32a) indicates that the HFS cannot be parametrized solely by ω and a ∗ in this case, making the proposed ellipse sweeping procedure uncertain.

The success of the PLL based on the phase lag of the third harmonic could lead to think that the SCBC and RCT would potentially benefit from a similar unfolding provided they work with the amplitude of the third harmonic. Unfortunately, as Figure 33a reveals, such an unfolding does not necessarily occur for this harmonic amplitude, making it uncertain whether these methods would be able to handle superharmonic resonances.

This series of tests on a superharmonic resonance showed that the assumptions of all methods based on natural parameter continuation can fail even with a single-degree-of-freedom system. Such resonances are nevertheless usually less important and/or relevant than primary resonances, that remain the main interest in nonlinear vibration testing.

## 8 Illustration with a clamped thin plate

We now move to a real structure with challenging nonlinear dynamics, namely multiple nonlinear

modes with a distributed geometrical nonlinearity. This set-up does not feature the same exceptional repeatability as the electronic Duffing oscillator and therefore gives a more realistic idea of what can be expected from a test with a control-based method.

## 8.1 Experimental setup

A thin rectangular steel plate of dimensions 500 mm × 300 mm × 0.5 mm and clamped on its shortest edges was considered as the next experimental setup. The plate was excited with an electrodynamic shaker (TIRA TV 51075) in current mode, and the response of the plate was monitored with an impedance head (DYTRAN 5860B) as well as a laser vibrometer (Polytec NLV2500-5). In addition to these measurements from the structure, the voltage across the electrodes of the shaker was also monitored.

Figure 36 shows pictures of the experimental clamped plate. The clamping device was realized with two cantilever arms fixed on an optical table. The plate was then fixed to the arms with bolts, nuts and wide washers to distribute the load as evenly as possible on the clamped side. The cantilever arms were observed to be rigid enough to provide an approximate clamping condition for the bending dynamics of the plate. However, they were not rigid enough with regards to the plate axial stiffness, and a bar was added on top of the plate to stiffen the assembly, allowing for the observation of clear nonlinear midplane stretching effects. Details about the dimensions of the setup are given in the supplementary materials.

Compared with the electronic Duffing oscillator, several new challenges arise. The most problematic one is the occurrence of shaker-structure interactions. These interactions have two main consequences when a harmonic drive signal is fed to the amplifier of the shaker: ( ) the fundameni tal forcing harmonic amplitude is frequency- and amplitude-dependent and ( ii ) harmonics appear in the force. These issues were overlooked in this work by considering the signal fed to the amplifier, V in (proportional to the current in the shaker), as the input to the system (as opposed to the force measured by the impedance head). The price to pay for this convenience is that shaker-structure interactions can then substantially dampen the structural modes if V in is considered as the input.

The shaker was thus placed following a tradeoff between conflicting requirements: on the one hand, it should be placed on a responsive location of the targeted modes to excite them well, but on the other hand it should be placed close to the clamping to minimize shaker-structure interactions.

## 8.2 Linear modal analysis

A low-level pseudorandom excitation signal was applied to the plate to compute the FRFs when the system behaves linearly. The PolyMAX modal identification method [149] was then used to extract the linear modal characteristics of the first nine modes (below 100 Hz). The characteristics of the first five modes, which will be tested in nonlinear regimes of motion, are gathered in Table 4.

Table 4 : Modal characteristics of the first five modes of the plate.

|   Mode |   Frequency (Hz) |   Damping ratio (%) |
|--------|------------------|---------------------|
|      1 |            15.51 |                0.21 |
|      2 |            25.86 |                0.89 |
|      3 |            41.11 |                0.91 |
|      4 |            49.02 |                0.33 |
|      5 |            52.9  |                0.33 |

Figures 37 and 38 show the two FRFs related to the NFRs of interest in the sequel, namely the shaker voltage and laser velocity.

Figure 39 shows the FRF of the force applied to the structure, as measured by the impedance head. Clear shaker-structure interactions are visible, as the FRF is far from being constant. This is due to the high flexibility of the plate combined with the relatively high mass of the impedance head [150]. Nevertheless, the shaker does not appear to dampen the modes of interest excessively. Controlling the force would prove at best extremely challenging anyway, because this FRF has a nonminimum phase zero (a pair of complex conjugate zeros with positive real part) around 137 Hz, which is very likely to destabilize any attempt of linear feedback control [151].

The dynamics of the plate in higher-frequency regions feature a high modal density. While this might not be a limitation for several methods, preliminary SWS tests showed that, in a nonlinear regime of motion, extremely complex phenomena

Fig. 36 : Pictures of the experimental setup: (a): overview of the clamped plate and (b): view of the backside of the clamping device.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-71ad0c18-35dd-4f17-a85d-c94d1a690857.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-cb5e2817-0b04-47b9-aa6f-0417f6db8845.jpg)

Fig. 37 : Bode plot of the shaker voltage FRF.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-219adcda-76f7-4922-9640-8f857f97c95b.jpg)

took place therein, with the appearance of internal and combination resonances. Both these phenomena are generally associated with the appearance of quasiperiodic oscillations, that none of the state-of-the-art control-based method can handle. This high frequency range was thus not investigated further in this work.

## 8.3 Implementation details

The CBC requires feedback control similar to direct velocity feedback. As is well known in the field of linear control of flexible structures, collocation of the applied forcing and the feedback signal is almost mandatory to hope for a stable closedloop system [151]. Loosely speaking, collocation allows for the feedback controller to work with power conjugate quantities, thereby controlling the power flow out of the structure. In particular, the control law can be designed to avoid power input to the structure which could lead to closedloop instabilities. In most works about the CBC, one-degree-of-freedom-like structures were considered, thereby dismissing this question. For the plate, flexibility effects bear a strong importance. To allow for collocated control, one possibility is to measure the motion of the moving part of the

Fig. 38 : Bode plot of the laser velocity FRF.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-0aa427e0-7d88-4eac-bc79-534de346fbc0.jpg)

Fig. 39 : Bode plot of the impedance head force FRF.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-d9df2c6c-a5bd-46e5-985e-2e226ac53767.jpg)

shaker [91, 152]. Alternatively, the strategy used in this work consists in sensing the voltage across the electrodes of the shaker (which is measured by the shaker amplifier). Since the current and voltage are power conjugates, a truly collocated control is achieved. More specifically, this is confirmed in Figure 37 , where the phase bounded between ± 90 ◦ ensures an interlacing pole-zero pattern, typical of a collocated sensor-actuator pair [151]. A first-order, unit-gain low-pass filter with a cutoff frequency at 500 Hz was nevertheless used on the voltage signal to avoid the excitation of highfrequency dynamics and the undesirable effects associated with digital control.

For the PLL and RCT, such a collocation is not necessary. Since the phase evolution around the first considered mode shows better monotonicity for the velocity signal (cf. Figure 38), the latter was chosen as the feedback signal for PLL.

Finally, the CBC-FD and RCT methods were unsuccessful. Concerning the former, this is probably due to the stronger noise in this setup. This lead the method to diverge quite far away from

the NFR at the start of the test at low amplitudes. Since this behavior was worrisome near resonance, the CBC-FD was not tested in this case. Better results can be hoped for if noise is handled properly [82], but this is beyond the scope of the present work. The RCT method also had stability issues similar to those encountered for the electronic Duffing. However, since these systematically lead to an overload, the RCT could not be performed fully to guarantee the structural integrity of the plate.

The tested methods, namely SCBC, PLL and ACBC, were compared once again with a SWS for reference. All parameters are given in the supplementary materials.

## 8.4 Broadband NFR

We start off the nonlinear analysis with the aim to obtain an NFR encompassing multiple resonant modes by investigating the frequency range between 10 and 60 Hz. In this case, several of the control-based methods are impractical or even inapplicable. The PLL is inherently hampered by the fact that the phase is non-monotonous (cf. Figures 37 and 38). In addition, the phase sensitivity away from resonance is rather small, indicating that phase parametrization is inadequate, and phase noise may cause large frequency noise. As for SCBC and RCT, their main drawback is that they require a certain response amplitude at all frequencies. Near antiresonances, reaching such an amplitude would require unreasonably large forcing levels. Furthermore, since the frequency band is quite large, the structure would have to be excited at maximum response amplitude for a long time, which is obviously undesirable. Hence, while these three methods are suitable to measure the response of a structure around its resonance, they are unsuited for broadband NFR testing. This could nevertheless be different for the two latter if an adapted sampling procedure was chosen instead of a regular grid in the ( ω, a ∗ ) space.

Using the ACBC for different levels of excitation, one can obtain a global picture of the plate dynamics. Figures 40 and 41 feature the voltage and velocity NFRs, respectively, for five different drive amplitudes. Although the test started at 10 Hz, the frequency range shown therein starts at 20 Hz because the first mode is barely visible and barely excited by the shaker (as seen in Figure 37), thus remaining linear and being rather uninteresting. Each of these NFRs were obtained under 8 minutes of testing.

The ACBC method was tested at the highest level of excitation ( V in = 0.05 V) against SWS for reference. Figure 42 presents the voltage NFR. The second and third modes are hardening, while the fifth mode is initially softening, but seems to start its hardening phase. The softening trend may come from a small but inevitable initial curvature of the plate. The results agree almost perfectly with SWS, except for the inevitable hysteresis, small transients near the jump phenomena, and for the fourth mode. The ACBC sometimes had trouble converging in the vicinity of the fourth mode. The issue is thought to come from its poor observability and controllability (this mode being almost absent from the considered driving-point NFR at a linear level, cf. Figure 37). Superharmonic resonances are appearing throughout the test. They are strong enough to be noticed as peaks in the NFR, but not enough to form loops; hence the ACBC method was able to characterize them.

Figure 43 shows the velocity NFR at the laser measurement point. Mode 5 is more visible therein, as well as the superharmonic resonances.

All modes feature an interesting and strongly nonlinear behavior at the highest considered forcing amplitude. To close this analysis, we investigate the dynamics of the second mode more in depth. It should be noted that since the frequency of mode 5 is about twice that of mode 2, the structure is in near-internal resonance conditions, resulting in an unconventional shape of the resonance peak.

## 8.5 Second mode

We now focus on a single resonance, the second flexible mode of the assembly. Tests that were performed for the Duffing oscillator are repeated for the considered methods, namely SWS, SCBC, ACBC and PLL. All tests were performed on the same day within a time span of 3 h 30 min, and Table 5 provides the time taken by each method (NB: as in the Duffing case, the times for SWS and PLL were chosen based on the fastest method). Figure 44a presents the results for the voltage NFR. All methods agree relatively well away from resonance, and in particular along the unstable

Fig. 40 : Shaker voltage NFR of the clamped plate for different drive signal amplitudes: 0 01 V ( . ), 0 02 . V ( ), 0 03 V ( . ), 0 04 V ( . ) and 0 05 V ( . ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-142d934f-b181-43b1-9768-c4caa34b9c40.jpg)

Fig. 41 : Laser velocity NFR of the clamped plate for different drive signal amplitudes: 0 01 V ( . ), 0 02 . V ( ), 0 03 V ( . ), 0 04 V ( . ) and 0 05 V ( . ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-d87a5d50-75b9-4fc9-803e-f2012161a5d1.jpg)

branch (except of course for SWS), but small discrepancies are visible near resonance. These discrepancies are more marked for the velocity NFR, as illustrated by Figure 44b, where two peaks are generally visible.

Table 5 : Times taken by each test for the primary resonance of the second mode of the plate.

Adaptive filters provide means to analyze the reason for the presence of a double peak in the NFR. As Figure 45 reveals, the second harmonic becomes more prominent around 26.5 Hz, which

| Method            |   Total time (s) | Time per point (s)   |
|-------------------|------------------|----------------------|
| SWS (up and down) |              190 | /                    |
| SCBC (adaptive)   |              397 | 0.58                 |
| ACBC              |              190 | 1.13                 |
| PLL               |              190 | /                    |

is slightly more than half the linear resonance frequency of mode 5. A 2:1 superharmonic resonance

Fig. 42 : Shaker voltage NFR of the clamped plate for V in = 0 05 V: SWS-up ( . ), SWS-down ( ) and ACBC ( ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-e491d6b5-95ca-4360-a266-ac5427d28cde.jpg)

Fig. 43 : Edge velocity NFR of the clamped plate for V in = 0 05 V: SWS-up ( . ), SWS-down ( ) and ACBC ( ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-e69684a2-12e7-4e0a-a042-7bbb740bbb5e.jpg)

of mode 5 thus occurs right before the primary resonance of mode 2, explaining the appearance of the second peak.

In addition, Figure 46 shows the evolution of the phase of the first harmonic for the voltage and velocity NFRs obtained with the ACBC. Similarly to the FRFs (Figures 37 and 38), the voltage NFR does not feature a monotonous phase decrease around the resonance, whereas the velocity NFR

shows this trend overall. Moreover, the 2:1 resonance manifests itself through a small phase trough, which is responsible for a small jump of the PLL method. With attention to detail, this jump can be observed in Figure 44, approximately between 26.5 and 27 Hz, wherein the PLL curve has less resolution.

To assess the source of discrepancies, ten repeatability tests were performed for the ACBC and PLL. Statistics on the point cloud formed by

Fig. 44 : Shaker voltage (a) and laser velocity (b) NFRs around the second mode: SWS-up ( ), SWS-down ( ), ACBC ( ), PLL ( ) and SCBC ( ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-6788c875-8227-4f26-8694-73506e4e4a48.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-adf1db71-5590-4bfb-86ad-60f914b612ab.jpg)

the measurement were then obtained by fitting it with a spline with eleven control points and Gaussian process regressions [153], as explained in the supplementary materials. Figure 47 presents the results for the voltage NFR. Except for an outlier for the ACBC, both types of test yield relatively repeatable results, although the differences between test partially explain the discrepancies observed in Figure 44. The remaining reasons may come from the repeated tests at high amplitude

Fig. 45 : Harmonic amplitudes of the voltage signal: zeroth ( ), first ( ), second ( ), third ( ), fourth ( ) and fifth ( ) harmonic amplitudes, and total amplitude ( ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-5d122bb9-f736-4ae2-b340-7f30badf092a.jpg)

Fig. 46 : Phase of the first harmonic for the voltage ( ) and velocity ( ) NFRs.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-455a673e-efd3-4d68-bde2-5d6b09c66b83.jpg)

(it can be deduced from Figure 44b that the plate undergoes a displacement amplitude of more than twice its thickness at the laser measurement point at resonance) over a relatively long time span. Indeed, there was at least 30 min of heavy testing in between each curve in Figure 44. It is also interesting to note that the averaging process makes the small peak associated with the superharmonic resonance appear more clearly in both cases.

Fig. 47 : Average NFR (dark curve), and 95% confidence interval (light area) obtained from 10 repeatability tests: ACBC (a) and PLL (b).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-1d88d959-dd09-4a8e-9b25-0a318fbf6a47.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00138v1-figure-b369a5ac-ad79-4612-8749-1d8884d4a549.jpg)

## 9 Conclusion

With the growing need to exhaustively test nonlinear systems, experimental continuation approaches are emerging as robust and efficient means to obtain bifurcation diagrams experimentally, encompassing the open-loop unstable responses of the system under test. By reviewing their essential ingredients, namely continuation and control, this article presented the main variants in experimental continuation in a unified manner. In particular, a common framework for the derivative-free experimental continuation approaches which are usually presented in an individual way, namely SCBC, PLL and RCT, was given. It also objectively compared these different methods with an electronic Duffing oscillator and a thin clamped plate. The main conclusions coming from this work and the literature are outlined below, summarizing the strengths and weaknesses of the different methods in their current form.

SWS is typically the first type of test that one tries on a system. In spite of its inherent deficiencies, the fact that this method is open-loop makes it vastly simpler than the other methods presented in this work. It is an excellent approach to obtain preliminary results and to have an idea of the orders of magnitude involved with the system under test.

CBC offers the most mature and comprehensive suite for experimental continuation. With its numerous variants, it is possible to trade between simplicity and advanced nonlinear testing capabilities. The controller being separated from the continuation algorithm, it is possible to have an intuitive understanding of its effects. Yet, CBC can be the most complex of the approaches outlined herein.

PLL offers a quick and easy testing solution with a simple architecture. It is particularly effective and fast for testing around a specific resonance, obtaining the NFR close to this resonance, or tracking its backbone. Yet, tuning its parameters can sometimes prove difficult for the unstable parts of the NFR, and the fact that the phase folds as soon as there are multiple resonances (including secondary ones) inherently limits this approach for broadband testing.

The RCT is an attractive method for the possibility to implement it with commercial software, and the simple postprocessing associated with tests made at constant vibration amplitude. There remains yet to be formal proofs of its stabilization properties (or the conditions for them), as examples herein showed that this is not guaranteed.

Finally, we look back at the proposed ACBC procedure. The fact that it is based on arclength continuation and derivative-free makes it a powerful and general tool, allowing the experimenter to obtain complex bifurcation diagrams, as featured

in the plate example. Nevertheless, limitations were also highlighted, e.g., when the bifurcation diagram features a loop, or when adaptive filters fail to secure non-invasiveness. Future works will address these issues.

Since the pioneering work of Sieber and Krauskopf [48], substantial progress has been made over the last twenty years in experimental continuation, but great challenges also lie ahead. Some of them are listed below:

- 1. Gain selection and adaptivity: In almost all works, the control parameters are selected manually, often on a trial and error basis. A few theoretical studies have been made on a singledegree-of-freedom Duffing oscillator [28, 68, 115], and a general but time-consuming parameter space exploration procedure was proposed in [70]. More systematic procedures should be developed to enable the use of control-based methods for non-expert users. We note that recent works proposed to use adaptive procedures to control the system with solid theoretical properties [154-156]. These methods have stringent requirements for an experimental implementation, such as full-state feedback, but are very promising and represent a significant step forward. Alternatively, practical ways of setting the parameters of a PLL with limited information from the plant were proposed in [145]. Finally, the data-driven approach based on parameter sweeps recently proposed in [157] could form another interesting solution.
- 2. Stability and bifurcations: When a test is successfully carried out, the closed-loop stability of every equilibria is ascertained. However, since the system of interest is the open-loop one, it is interesting to know whether the measured equilibria are open-loop stable or not. Procedures based on temporarily shutting the control down [81] or on real-time linearization [97, 98] were proposed for the CBC, and similar procedures could be developed for the other control-based approaches. Going beyond that, different bifurcations could also be detected [158, 159] and even tracked in an experimental setting, as in the pioneering works of Renson et al [86] and Shen et al [160] (in the field of statics).
- 3. Isolated solutions: From a numerical standpoint, isolated solutions can be challenging to

find with traditional continuation tools, but may be reachable in experiments. As a few works demonstrated [81, 87, 90, 148], controlbased methods could offer robust tools to trace out a complete isola, but reaching or even detecting them in a systematic way still represents a challenge.

- 4. Aperiodic responses: None of the controlbased method is currently able to deal with aperiodic (quasiperiodic or chaotic) responses, even though they may generally appear in nonlinear systems. As stated in Section 3.2, it is possible to stabilize unstable periodic orbits embedded in a chaotic attractor [60, 64], but, to the authors' knowledge, little work has been done to stabilize an aperiodic attractor itself.
- 5. Actuator-system interaction: In most cases, a structure will not be perfectly actuated. This problem pertains to the breadth of nonlinear vibration testing problems [45] and not just to control-based methods. Although some attempts have been made [161, 162], no ideal solution currently exists in the literature.
- 6. Large-scale systems: Except a few notable exceptions, control-based methods were almost always tested on systems of relatively low dimensionality, being either single-degree-offreedom-like, or focusing on a single resonance widely spaced from the others. To achieve their full potential and provide convincing results for a future industrial use, control-based methods should turn toward more challenging, largescale systems. The use of multiple sensors and/or actuators and the study of underactuated problems also remain largely uncovered.

## References

- [1] Nayfeh, A.H., Mook, D.T.: Nonlinear Oscillations, p. 720. Wiley, Weinheim (1995). https://doi.org/10.1002/9783527617586 . https://onlinelibrary.wiley.com/doi/book/ 10.1002/9783527617586
- [2] Guckenheimer, J., Holmes, P.: Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Applied Mathematical Sciences, vol. 42. Springer, New York, NY (1983). https://doi.org/10.1007/ 978-1-4612-1140-2 . http://link.springer. com/10.1007/978-1-4612-1140-2

- [3] Kuznetsov, Y.A.: Elements of Applied Bifurcation Theory. Applied Mathematical Sciences, vol. 112. Springer, New York, NY (2004). https://doi.org/10.1007/ 978-1-4757-3978-7 . http://link.springer. com/10.1007/978-1-4757-3978-7
- [4] Wriggers, P.: Nonlinear Finite Element Methods. Springer, Berlin Heidelberg (2008)
- [5] Keller, H.B.: Numerical solution of bifurcation and nonlinear eigenvalue problems. In: Rabinowitz, P.H. (ed.) Applications of Bifurcation Theory: Proceedings of an Advanced Seminar, pp. 359-384. Academic Press, New York, NY (1977)
- [6] Allgower, E.L., Georg, K.: Introduction to Numerical Continuation Methods. SIAM, Philadelphia (2003)
- [7] Dhooge, A., Govaerts, W., Kuznetsov, Y.A.: Matcont: A matlab package for numerical bifurcation analysis of odes. ACM Transactions on Mathematical Software 29 (2), 141164 (2003) https://doi.org/10.1145/779359. 779362
- [8] Govaerts, W.J.: Numerical Methods for Bifurcations of Dynamical Equilibria. SIAM, Philadelphia (2000)
- [9] Nayfeh, A.H., Balachandran, B.: Applied Nonlinear Dynamics. Wiley, Weinheim (1995). https://doi.org/10.1002/ 9783527617548 . https://onlinelibrary.wiley. com/doi/book/10.1002/9783527617548
- [10] Knowles, J., Krauskopf, B., Lowenberg, M.: Numerical continuation applied to landing gear mechanism analysis. Journal of Aircraft 48 (4), 1254-1262 (2011) https://doi.org/10. 2514/1.C031247
- [11] Detroux, T., Renson, L., Masset, L., Kerschen, G.: The harmonic balance method for bifurcation analysis of large-scale nonlinear mechanical systems. Computer Methods in Applied Mechanics and Engineering 296 , 18-38 (2015) https://doi.org/10.1016/ j.cma.2015.07.017 arXiv:1604.05621
- [12] Dankowicz, H., Schilder, F.: An extended continuation problem for bifurcation analysis in the presence of constraints. Journal of Computational and Nonlinear Dynamics 6 (3) (2011) https://doi.org/10.1115/1. 4002684
- [13] Freour, V., Guillot, L., Masuda, H., Usa, S., Tominaga, E., Tohgi, Y., Vergez, C., Cochelin, B.: Numerical continuation of a physical model of brass instruments: Application to trumpet comparisons. Journal of the Acoustical SOciety of America 148 (2), 748 (2020) https://doi.org/10.1121/10.0001603
- [14] Haller, G., Ponsioen, S.: Nonlinear normal modes and spectral submanifolds: existence, uniqueness and use in model reduction. Nonlinear Dynamics 86 (3), 1493-1534 (2016) https://doi.org/10.1007/s11071-016-2974-z arXiv:1602.00560
- [15] Vizzaccaro, A., Shen, Y., Salles, L., Blahoˇ, s J., Touz´ e, C.: Direct computation of nonlinear mapping via normal form for reduced-order models of finite element nonlinear structures. Computer Methods in Applied Mechanics and Engineering 384 (2021) https://doi.org/10.1016/j.cma.2021. 113957 arXiv:2009.12145
- [16] Jain, S., Haller, G.: How to compute invariant manifolds and their reduced dynamics in high-dimensional finite element models. Nonlinear Dynamics 107 (2), 1417-1450 (2022) https://doi.org/10.1007/ s11071-021-06957-4 arXiv:2103.10264
- [17] Touz´ e, C., Vizzaccaro, A., Thomas, O.: Model order reduction methods for geometrically nonlinear structures: a review of nonlinear techniques. Nonlinear Dynamics 105 (2), 1141-1190 (2021) https: //doi.org/10.1007/s11071-021-06693-9 arXiv:2107.05077
- [18] Cenedese, M., Ax˚s, J., B¨uerlein, B., Avila, a a K., Haller, G.: Data-driven modeling and prediction of non-linearizable dynamics via spectral submanifolds. Nature Communications 13 (1), 872 (2022) https://doi.org/10. 1038/s41467-022-28518-y

- [19] Virgin, L.: Introduction to Experimental Nonlinear Dynamics. Cambridge University Press, Cambridge (2000)
- [20] Kerschen, G., Worden, K., Vakakis, A.F., Golinval, J.-C.: Past, present and future of nonlinear system identification in structural dynamics. Mechanical Systems and Signal Processing 20 (3), 505-592 (2006) https:// doi.org/10.1016/j.ymssp.2005.04.008
- [21] No¨ el, J.P., Kerschen, G.: Nonlinear system identification in structural dynamics: 10 more years of progress. Mechanical Systems and Signal Processing 83 , 2-35 (2017) https: //doi.org/10.1016/j.ymssp.2016.07.020
- [22] Anastasio, D., Marchesiello, S.: Nonlinear frequency response curves estimation and stability analysis of randomly excited systems in the subspace framework. Nonlinear Dynamics 111 (9), 8115-8133 (2023) https: //doi.org/10.1007/s11071-023-08280-6
- [23] Sieber, J., Krauskopf, B.: Control based bifurcation analysis for experiments. Nonlinear Dynamics 51 (3), 365-377 (2008) https: //doi.org/10.1007/s11071-007-9217-2
- [24] Sieber, J., Gonzalez-Buelga, A., Neild, S.A., Wagg, D.J., Krauskopf, B.: Experimental continuation of periodic orbits through a fold. Physical Review Letters 100 (24), 1-4 (2008) https: //doi.org/10.1103/PhysRevLett.100.244101 arXiv:0804.0320
- [25] Barton, D.A.W., Sieber, J.: Systematic experimental exploration of bifurcations with noninvasive control. Physical Review E 87 (5), 052916 (2013) https://doi.org/10. 1103/PhysRevE.87.052916
- [26] Peter, S., Leine, R.I.: Excitation power quantities in phase resonance testing of nonlinear systems with phase-locked-loop excitation. Mechanical Systems and Signal Processing 96 , 139-158 (2017) https://doi.org/ 10.1016/j.ymssp.2017.04.011
- [27] Karaa˘ ga¸lı, c T., ¨ zg¨ven, O u H.N.: Experimental identification of backbone curves
- of strongly nonlinear systems by using response-controlled stepped-sine testing (rct). Vibration 3 (3), 266-280 (2020) https://doi.org/10.3390/vibration3030019
- [28] Abeloos, G.: Control-based methods for the identification of nonlinear structures. PhD thesis, University of Li`ge (2022). https:// e hdl.handle.net/2268/295414
- [29] Kovacic, I., Brennan, M.J.: The Duffing Equation. John Wiley &amp; Sons, Ltd, Chichester, UK (2011). https://doi.org/10.1002/ 9780470977859 . http://doi.wiley.com/10. 1002/9780470977859
- [30] Seydel, R.: Practical Bifurcation and Stability Analysis vol. 5. Springer, New York Dordrecht Heidelberg London (2009)
- [31] Krack, M., Gross, J.: Harmonic Balance for Nonlinear Vibration Problems. Mathematical Engineering, p. 159. Springer, Cham (2019). https://doi.org/10.1007/978-3-030-14023-6 . http://link.springer.com/10.1007/ 978-3-030-14023-6
- [32] Dankowicz, H., Schilder, F.: Recipes for Continuation. SIAM, Philadelphia (2013)
- [33] Doedel, E., Keller, H.B., Kernevez, J.P.: Numerical analysis and control of bifurcation problems (i): Bifurcation in finite dimensions. International Journal of Bifurcation and Chaos 01 (03), 493-520 (1991) https: //doi.org/10.1142/S0218127491000397
- [34] Cochelin, B., Vergez, C.: A high order purely frequency-based harmonic balance formulation for continuation of periodic solutions. Journal of Sound and Vibration 324 (1-2), 243-262 (2009) https://doi.org/10.1016/j. jsv.2009.01.054
- [35] Krauskopf, B., Osinga, H.M., Gal´n-Vioque, a J.: Numerical Continuation Methods for Dynamical Systems vol. 2. Springer, Dordrecht (2007)
- [36] Haselgrove, C.B.: The solution of non-linear

equations and of differential equations with two-point boundary conditions. The Computer Journal 4 (3), 255-259 (1961) https: //doi.org/10.1093/comjnl/4.3.255

- [37] Riks, E.: An incremental approach to the solution of snapping and buckling problems. International Journal of Solids and Structures 15 (7), 529-551 (1979) https:// doi.org/10.1016/0020-7683(79)90081-7
- [38] Crisfield, M.A.: A fast incremental/iterative solution procedure that handles 'snapthrough'. Computers &amp; Structures 13 (13), 55-62 (1981) https://doi.org/10.1016/ 0045-7949(81)90108-5
- [39] Lewandowski, R.: Non-linear, steady-state vibration of structures by harmonic balance/finite element method. Computers &amp; Structures 44 (1-2), 287-296 (1992) https: //doi.org/10.1016/0045-7949(92)90248-X
- [40] Sundararajan, P., Noah, S.T.: Dynamics of forced nonlinear systems using shooting/are-length continuation methodapplication to rotor systems. Journal of Vibration and Acoustics, Transactions of the ASME 119 (1), 9-20 (1997) https://doi. org/10.1115/1.2889694
- [41] Slotine, J.-J.E.: Applied Nonlinear Control. Prentice-Hall, Englewood Cliffs (1991)
- [42] Parlitz, U., Lauterborn, W.: Superstructure in the bifurcation set of the Duffing equation. Physics Letters A 107 (8), 351-355 (1985) https://doi.org/10.1016/ 0375-9601(85)90687-5
- [43] Marchionne, A., Ditlevsen, P., Wieczorek, S.: Synchronisation vs. resonance: Isolated resonances in damped nonlinear oscillators. Physica D: Nonlinear Phenomena 380381 , 8-16 (2018) https://doi.org/10.1016/j. physd.2018.05.004
- [44] Stoker, J.J.: Nonlinear Vibrations in Mechanical and Electrical Systems vol. 2. Interscience publishers, New York (1950)
- [45] Pacini, B.R., Kuether, R.J., Roettgen,
- D.R.: Shaker-structure interaction modeling and analysis for nonlinear force appropriation testing. Mechanical Systems and Signal Processing 162 (May 2021), 108000 (2022) https://doi.org/10.1016/j. ymssp.2021.108000
- [46] Virgin, L.N., Todd, M.D., Begley, C.J., Trickey, S.T., Dowell, E.H.: Basins of attraction in experimental nonlinear oscillators. International Journal of Bifurcation and Chaos 08 (03), 521-533 (1998) https://doi. org/10.1142/S0218127498000334
- [47] Ludeke, C.A.: Predominantly subharmonic oscillations. Journal of Applied Physics 22 (11), 1321-1326 (1951) https://doi.org/ 10.1063/1.1699858
- [48] Sieber, J., Krauskopf, B.: Control-based continuation of periodic orbits with a timedelayed difference scheme. International Journal of Bifurcation and Chaos 17 (08), 2579-2593 (2007) https://doi.org/10.1142/ S0218127407018646
- [49] Mayr, O.: The origins of feedback control. Scientific American 223 (4), 110-118 (1970) https://doi.org/10.1038/ scientificamerican1070-110
- [50] Fuller, A.T.: Directions of research in control. Automatica 1 (4), 289-296 (1963) https: //doi.org/10.1016/0005-1098(63)90013-X
- [51] Bennett, S.: A brief history of automatic control. IEEE Control Systems 16 (3), 17-25 (1996) https://doi.org/10.1109/37.506394
- [52] Fuller, A.T.: The early development of control theory. Journal of Dynamic Systems, Measurement, and Control 98 (2), 109-118 (1976) https://doi.org/10.1115/1.3426994
- [53] Fuller, A.T.: The early development of control theory. ii. Journal of Dynamic Systems, Measurement, and Control 98 (3), 224-235 (1976) https://doi.org/10.1115/1.3427026
- [54] Atherton, D.P.: Early developments in nonlinear control. IEEE Control Systems 16 (3), 34-43 (1996) https://doi.org/10.1109/37.

## 506396

- [55] Kokotovi´ c, P., Arcak, M.: Constructive nonlinear control: a historical perspective. Automatica 37 (5), 637-662 (2001) https:// doi.org/10.1016/S0005-1098(01)00002-4
- [56] Iqbal, J., Ullah, M., Khan, S.G., Khelifa, B., ´ ukovi´, C c S.: Nonlinear control systems -a brief overview of historical and recent advances. Nonlinear Engineering 6 (4), 301-312 (2017) https://doi.org/10. 1515/nleng-2016-0077
- [57] Khalil, H.K.: Nonlinear Control. Pearson Education, Harlow (2002)
- [58] Sepulchre, R., Jankovic, M., Kokotovic, P.V.: Constructive Nonlinear Control. Springer, London (2012)
- [59] Wagg, D., Neild, S.: Nonlinear Vibration with Control. Solid Mechanics and Its Applications, vol. 218. Springer, Cham (2015). https://doi.org/10.1007/978-3-319-10644-1 . https://link.springer.com/10.1007/ 978-3-319-10644-1
- [60] Ott, E., Grebogi, C., Yorke, J.A.: Controlling chaos. Physical Review Letters 64 (11), 1196-1199 (1990) https://doi.org/10.1103/ PhysRevLett.64.1196
- [61] Shinbrot, T., Grebogi, C., Yorke, J.A., Ott, E.: Using small perturbations to control chaos. Nature 363 (6428), 411-417 (1993) https://doi.org/10.1038/363411a0
- [62] Boccaletti, S., Grebogi, C., Lai, Y.C., Mancini, H., Maza, D.: The control of chaos: Theory and applications. Physics Report 329 (3), 103-197 (2000) https://doi.org/10. 1016/S0370-1573(99)00096-4
- [63] Fradkov, A.L., Evans, R.J.: Control of chaos: Methods and applications in engineering. Annual Reviews in Control 29 (1), 33-56 (2005) https://doi.org/10.1016/j.arcontrol. 2005.01.001
- [64] Pyragas, K.: Continuous control of chaos by self-controlling feedback. Physics Letters A
- 170 (6), 421-428 (1992) https://doi.org/10. 1016/0375-9601(92)90745-8
- [65] Carroll, T.L., Triandaf, I., Schwartz, I., Pecora, L.: Tracking unstable orbits in an experiment. Physical Review A 46 (10), 6189-6192 (1992) https://doi.org/10.1103/ PhysRevA.46.6189
- [66] Gills, Z., Iwata, C., Roy, R., Schwartz, I.B., Triandaf, I.: Tracking unstable steady states: Extending the stability regime of a multimode laser system. Physical Review Letters 69 (22), 3169-3172 (1992) https:// doi.org/10.1103/PhysRevLett.69.3169
- [67] Neville, R.M., Groh, R.M.J., Pirrera, A., Schenk, M.: Beyond the fold: experimentally traversing limit points in nonlinear structures. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 476 (2233), 20190576 (2020) https: //doi.org/10.1098/rspa.2019.0576
- [68] Tatzko, S., Kleyman, G., Wallaschek, J.: Continuation methods for lab experiments of nonlinear vibrations. GAMM Mitteilungen 46 (2), 1-13 (2023) https://doi.org/10. 1002/gamm.202300009
- [69] Hayashi, S., Gutschmidt, S., Murray, R., Krauskopf, B.: Experimental bifurcation analysis of a clamped beam with designed mechanical nonlinearity. Nonlinear Dynamics (2024) https: //doi.org/10.1007/s11071-024-09873-5
- [70] Bureau, E., Schilder, F., Ferreira Santos, I., Juel Thomsen, J., Starke, J.: Experimental bifurcation analysis of an impact oscillator-tuning a non-invasive control scheme. Journal of Sound and Vibration 332 (22), 5883-5897 (2013) https://doi.org/10.1016/ j.jsv.2013.05.033
- [71] Bureau, E.: Experimental bifurcation analysis using control-based continuation. PhD thesis, Technical University of Denmark (2014). https://findit.dtu.dk/en/catalog/ 53be7cc5f7d9a21b490001de
- [72] Siettos, C.I., Kevrekidis, I.G., Maroudas,

- D.: Coarse bifurcation diagrams via microscopic simulators: A state-feedback control-based approach. International Journal of Bifurcation and Chaos 14 (01), 207-220 (2004) https: //doi.org/10.1142/S0218127404009193
- [73] Barton, D.A.W., Mann, B.P., Burrow, S.G.: Control-based continuation for investigating nonlinear experiments. Journal of Vibration and Control 18 (4), 509-520 (2012) https:// doi.org/10.1177/1077546310384004
- [74] Sieber, J., Krauskopf, B., Wagg, D., Neild, S., Gonzalez-Buelga, A.: Control-based continuation of unstable periodic orbits. Journal of Computational and Nonlinear Dynamics 6 (1), 1-9 (2011) https://doi.org/ 10.1115/1.4002101
- [75] Barton, D.A.W., Burrow, S.G.: Numerical continuation in a physical experiment: Investigation of a nonlinear energy harvester. Journal of Computational and Nonlinear Dynamics 6 (1), 1-6 (2011) https:// doi.org/10.1115/1.4002380
- [76] Kleyman, G., Jahn, M., Tatzko, S., Scheidt, L.P.-v.: A combined numericalexperimental approach for the damping evaluation of non-linear dissipative vibration systems. In: Lecture Notes in Applied and Computational Mechanics vol. 102, pp. 285-303. Springer, Cham (2024). https: //doi.org/10.1007/978-3-031-36143-2 15 . https://link.springer.com/10.1007/ 978-3-031-36143-2 { } 15
- [77] Iderstein, T., Wiebe, R.: Experimental path following of unstable static equilibria for snap-through buckling. Conference Proceedings of the Society for Experimental Mechanics Series 1 , 17-22 (2019) https: //doi.org/10.1007/978-3-319-74280-9 2
- [78] Shen, J., Groh, R.M.J., Schenk, M., Pirrera, A.: Experimental path-following of equilibria using newton's method. part i: Theory, modelling, experiments. International Journal of Solids and Structures 210-211 , 203-223 (2021) https://doi.org/10.1016/j. ijsolstr.2020.11.037
- [79] Beregi, S., Barton, D.A.W., Rezgui, D., Neild, S.A.: Real-time hybrid testing using iterative control for periodic oscillations (2023) arXiv:2312.06362
- [80] Melville, R.C., Suarez, A.: Experimental investigation of bifurcation behavior in nonlinear microwave circuits. IEEE Transactions on Microwave Theory and Techniques 65 (5), 1545-1559 (2017) https://doi.org/10. 1109/TMTT.2016.2640955
- [81] Bureau, E., Schilder, F., Elmeg˚rd, M., Sana tos, I.F., Thomsen, J.J., Starke, J.: Experimental bifurcation analysis of an impact oscillator-determining stability. Journal of Sound and Vibration 333 (21), 54645474 (2014) https://doi.org/10.1016/j.jsv. 2014.05.032
- [82] Schilder, F., Bureau, E., Santos, I.F., Thomsen, J.J., Starke, J.: Experimental bifurcation analysis-continuation for noisecontaminated zero problems. Journal of Sound and Vibration 358 , 251-266 (2015) https://doi.org/10.1016/j.jsv.2015.08.008
- [83] Blyth, M., Tsaneva-Atanasova, K., Marucci, L., Renson, L.: Numerical methods for control-based continuation of relaxation oscillations. Nonlinear Dynamics 111 (9), 7975-7992 (2023) https://doi.org/10.1007/ s11071-023-08288-y
- [84] Schilder, F., Dankowicz, H.: Continuation Core and Toolboxes (COCO). https:// sourceforge.net/projects/cocotools/
- [85] Cenedese, M., Haller, G.: How do conservative backbone curves perturb into forced responses? a Melnikov function analysis. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 476 (2234), 20190494 (2020) https://doi.org/10.1098/rspa.2019.0494 arXiv:1908.00721
- [86] Renson, L., Barton, D.A.W., Neild, S.A.: Experimental tracking of limitpoint bifurcations and backbone curves using control-based continuation. International Journal of Bifurcation and

- Chaos 27 (1), 1-19 (2017) https: //doi.org/10.1142/S0218127417300026
- [87] Woiwode, L., Krack, M.: Experimentally uncovering isolas via backbone tracking. Journal of Structural Dynamics, 122-143 (2024) https://doi.org/10.25518/2684-6500. 180
- [88] Kleyman, G., Paehr, M., Tatzko, S.: Application of control-based-continuation for characterization of dynamic systems with stiffness and friction nonlinearities. Mechanics Research Communications 106 , 103520 (2020) https://doi.org/10.1016/j. mechrescom.2020.103520
- [89] Renson, L., Gonzalez-Buelga, A., Barton, D.A.W., Neild, S.A.: Robust identification of backbone curves using control-based continuation. Journal of Sound and Vibration 367 , 145-158 (2016) https://doi.org/ 10.1016/j.jsv.2015.12.035
- [90] Renson, L., Shaw, A.D., Barton, D.A.W., Neild, S.A.: Application of control-based continuation to a nonlinear structure with harmonically coupled modes. Mechanical Systems and Signal Processing 120 , 449-464 (2019) https://doi.org/10.1016/j. ymssp.2018.10.008 arXiv:1808.01865
- [91] Abeloos, G., Renson, L., Collette, C., Kerschen, G.: Stepped and swept control-based continuation using adaptive filtering. Nonlinear Dynamics 104 (4), 3793-3808 (2021) https: //doi.org/10.1007/s11071-021-06506-z
- [92] Beregi, S.: Nonlinear analysis of the delayed tyre model with control-based continuation. Nonlinear Dynamics 110 (4), 3151-3165 (2022) https://doi.org/10.1007/ s11071-022-07796-7
- [93] Renson, L., Sieber, J., Barton, D.A.W., Shaw, A.D., Neild, S.A.: Numerical continuation in nonlinear experiments using local gaussian process regression. Nonlinear Dynamics 98 (4), 2811-2826 (2019) https: //doi.org/10.1007/s11071-019-05118-y
- [94] Beregi, S., Barton, D.A.W., Rezgui, D., Neild, S.A.: Robustness of nonlinear parameter identification in the presence of process noise using control-based continuation. Nonlinear Dynamics 104 (2), 885-900 (2021) https://doi.org/10.1007/ s11071-021-06347-w arXiv:2001.11008
- [95] Beregi, S., Barton, D.A.W., Rezgui, D., Neild, S.: Using scientific machine learning for experimental bifurcation analysis of dynamic systems. Mechanical Systems and Signal Processing 184 (June 2022), 109649 (2023) https://doi.org/10.1016/j. ymssp.2022.109649 arXiv:2110.11854
- [96] Lee, K.H., Barton, D.A.W., Renson, L.: Modelling of physical systems with a hopf bifurcation using mechanistic models and machine learning. Mechanical Systems and Signal Processing 191 (October 2022), 110173 (2023) https://doi.org/10. 1016/j.ymssp.2023.110173
- [97] Barton, D.A.W.: Control-based continuation: Bifurcation and stability analysis for physical experiments. Mechanical Systems and Signal Processing 84 , 5464 (2017) https://doi.org/10.1016/j.ymssp. 2015.12.039 arXiv:1506.04052
- [98] Dittus, A., Kruse, N., Barke, I., Speller, S., Starke, J.: Detecting stability and bifurcation points in control-based continuation for a physical experiment of the zeeman catastrophe machine. SIAM Journal on Applied Dynamical Systems 22 (2), 1275-1299 (2023) https://doi.org/10.1137/ 22M1503245
- [99] Tartaruga, I., Barton, D.A.W., Rezgui, D., Neild, S.A.: Experimental bifurcation analysis of a wing profile. In: International Forum on Aeroelasticity and Structural Dynamics : IFASD 2019, Savannah (2019). https://hdl.handle.net/1983/ 6774dc99-50b6-4e51-ba30-f06207581973
- [100] De Cesare, I., Salzano, D., Di Bernardo, M.,

Renson, L., Marucci, L.: Control-based continuation: A new approach to prototype synthetic gene networks. ACS Synthetic Biology 11 (7), 2300-2313 (2022) https://doi. org/10.1021/acssynbio.1c00632

- [101] Hayashi, S., Gutschmidt, S., Murray, R., Krauskopf, B.: Control-based continuation of an externally excited mems self-oscillator. In: ENOC 2024, Delft (2024)
- [102] Panagiotopoulos, I., Starke, J., Just, W.: Control of collective human behavior: Social dynamics beyond modeling. Physical Review Research 4 (4), 043190 (2022) https://doi.org/10.1103/PhysRevResearch. 4.043190
- [103] Panagiotopoulos, I., Starke, J., Sieber, J., Just, W.: Continuation with noninvasive control schemes: Revealing unstable states in a pedestrian evacuation scenario. SIAM Journal on Applied Dynamical Systems 22 (1), 1-36 (2023) https://doi.org/10.1137/ 22M1482032 arXiv:2203.02484
- [104] Babitsky, V.I.: Autoresonant mechatronic systems. Mechatronics 5 (5), 483-495 (1995) https://doi.org/10.1016/0957-4158(95) 00026-2
- [105] Sokolov, I.J., Babitsky, V.I.: Phase control of self-sustained vibration. Journal of Sound and Vibration 248 (4), 725-744 (2001) https: //doi.org/10.1006/jsvi.2001.3810
- [106] Mojrzisch, S., Wallaschek, J., Bremer, J.: An experimental method for the phase controlled frequency response measurement of nonlinear vibration systems. PAMM 12 (1), 253-254 (2012) https://doi.org/10. 1002/pamm.201210117
- [107] Davis, S., Bucher, I.: Automatic vibration mode selection and excitation; combining modal filtering with autoresonance. Mechanical Systems and Signal Processing 101 , 140-155 (2018) https://doi.org/10. 1016/j.ymssp.2017.08.009
- [108] Scheel, M.: Nonlinear modal testing of damped structures: Velocity
- feedback vs. phase resonance. Mechanical Systems and Signal Processing 165 (June 2021), 108305 (2022) https: //doi.org/10.1016/j.ymssp.2021.108305 arXiv:2108.06189
- [109] Abramovitch, D.: Phase-locked loops: a control centric tutorial. In: Proceedings of the 2002 American Control Conference (IEEE Cat. No.CH37301), vol. 1, pp. 1-15. American Automatic Control Council, Anchorage (2002). https://doi.org/10.1109/ACC. 2002.1024769 . http://ieeexplore.ieee.org/ document/1024769/
- [110] Balas, M.J.: Active control of flexible systems. Journal of Optimization Theory and Applications 25 (3), 415-436 (1978) https: //doi.org/10.1007/BF00932903
- [111] Niezrecki, C., Cudney, H.H.: Structural control using analog phase-locked loops. Journal of Vibration and Acoustics 119 (1), 104-109 (1997) https://doi.org/10.1115/1. 2889677
- [112] Connally, J.A., Brown, S.B.: Micromechanical fatigue testing. Experimental Mechanics 33 (2), 81-90 (1993) https://doi.org/10. 1007/BF02322482
- [113] Kern, D., Seemann, W.: Tracking of mechanical system parameters by phaselocked loops. PAMM 10 (1), 613-614 (2010) https://doi.org/10.1002/pamm.201010299
- [114] Kern, D., Brack, T., Seemann, W.: Resonance tracking of continua using self-sensing actuators. Journal of Dynamic Systems, Measurement, and Control 134 (5), 1-9 (2012) https://doi.org/10.1115/1.4006224
- [115] Denis, V., Jossic, M., Giraud-Audine, C., Chomette, B., Renault, A., Thomas, O.: Identification of nonlinear modes using phase-locked-loop experimental continuation and normal form. Mechanical Systems and Signal Processing 106 , 430-452 (2018) https: //doi.org/10.1016/j.ymssp.2018.01.014
- [116] Peter, S., Scheel, M., Krack, M., Leine, R.I.:

Synthesis of nonlinear frequency responses with experimentally extracted nonlinear modes. Mechanical Systems and Signal Processing 101 , 498-515 (2018) https://doi. org/10.1016/j.ymssp.2017.09.014

- [117] Scheel, M., Peter, S., Leine, R.I., Krack, M.: A phase resonance approach for modal testing of structures with nonlinear dissipation. Journal of Sound and Vibration 435 , 56-73 (2018) https://doi.org/10.1016/j.jsv. 2018.07.010 arXiv:2011.08500
- [118] Scheel, M., Weigele, T., Krack, M.: Challenging an experimental nonlinear modal analysis method with a new strongly friction-damped structure. Journal of Sound and Vibration 485 , 115580 (2020) https://doi.org/10.1016/j.jsv.2020.115580 arXiv:2011.08527
- [119] Peter, S., Schreyer, F., Leine, R.I.: A method for numerical and experimental nonlinear modal analysis of nonsmooth systems. Mechanical Systems and Signal Processing 120 , 793-807 (2019) https://doi.org/10. 1016/j.ymssp.2018.11.009
- [120] Scheel, M., Kleyman, G., Tatar, A., Brake, M.R.W., Peter, S., No¨ el, J.-P., Allen, M.S., Krack, M.: Experimental assessment of polynomial nonlinear state-space and nonlinear-mode models for near-resonant vibrations. Mechanical Systems and Signal Processing 143 , 106796 (2020) https://doi. org/10.1016/j.ymssp.2020.106796
- [121] Debeurre, M., Benacchio, S., Grolet, A., Grenat, C., Giraud-Audine, C., Thomas, O.: Phase resonance testing of highly flexible structures: Measurement of conservative nonlinear modes and nonlinear damping identification. Mechanical Systems and Signal Processing 215 (February), 111423 (2024) https://doi.org/10.1016/j. ymssp.2024.111423
- [122] Nagesh, M., Allemang, R.J., Phillips, A.W.: Characterization of nonlinearities in a structure using nonlinear modal testing methods. In: Conference Proceedings of the Society for Experimental Mechanics

Series vol. 1, pp. 167-178 (2022). https: //doi.org/10.1007/978-3-030-77135-5 19 . https://link.springer.com/10.1007/ 978-3-030-77135-5 { } 19

- [123] Givois, A., Giraud-Audine, C., De¨, J.-F., u Thomas, O.: Experimental analysis of nonlinear resonances in piezoelectric plates with geometric nonlinearities. Nonlinear Dynamics 102 (3), 1451-1462 (2020) https://doi. org/10.1007/s11071-020-05997-6
- [124] Givois, A., Tan, J.-J., Touz´ e, C., Thomas, O.: Backbone curves of coupled cubic oscillators in one-to-one internal resonance: bifurcation scenario, measurements and parameter identification. Meccanica 55 (3), 481-503 (2020) https://doi.org/10. 1007/s11012-020-01132-2
- [125] Schwarz, S., Kohlmann, L., Hartung, A., Gross, J., Scheel, M., Krack, M.: Validation of a turbine blade component test with frictional contacts by phase-locked-loop and force-controlled measurements. Journal of Engineering for Gas Turbines and Power 142 (5), 1-8 (2020) https://doi.org/10.1115/ 1.4044772
- [126] Tang, M., Stephan, C., B¨swald, M.: Phase o resonance method for nonlinear mechanical structures with phase locked loop control. Proceedings of ISMA 2020 - International Conference on Noise and Vibration Engineering and USD 2020 - International Conference on Uncertainty in Structural Dynamics, 1805-1818 (2020)
- [127] M¨ uller, F., Woiwode, L., Gross, J., Scheel, M., Krack, M.: Nonlinear damping quantification from phase-resonant tests under base excitation. Mechanical Systems and Signal Processing 177 (January), 109170 (2022) https: //doi.org/10.1016/j.ymssp.2022.109170
- [128] Abeloos, G., M¨ uller, F., Ferhatoglu, E., Scheel, M., Collette, C., Kerschen, G., Brake, M.R.W., Tiso, P., Renson, L., Krack, M.: A consistency analysis of phase-locked-loop testing and control-based

continuation for a geometrically nonlinear frictional system. Mechanical Systems and Signal Processing 170 (October 2021), 108820 (2022) https://doi.org/10.1016/j. ymssp.2022.108820

- [129] Peter, S.: Nonlinear experimental modal analysis and its application to the identification of nonlinear structures. PhD thesis, University of Stuttgart (2018). https://doi. org/10.2370/9783844061369 . https://www. researchgate.net/publication/329756057
- [130] Scheel, M.: Experimental nonlinear modal analysis -method development with particular focus on nonlinear damping. PhD thesis, University of Stuttgart (2022). https://www.dr.hut-verlag.de/ 978-3-8439-5187-6.html
- [131] Link, M., Boeswald, M., Laborde, S., Weiland, M., Calvi, A.: Non-linear experimental modal analysis and application to satellite vibration test data. ECCOMAS Thematic Conference -COMPDYN 2011: 3rd International Conference on Computational Methods in Structural Dynamics and Earthquake Engineering: An IACM Special Interest Conference, Programme (May 2014) (2011)
- [132] Karaa˘ ga¸lı, T., Ozg¨ven, H.N.: Experimenc ¨ u tal modal analysis of nonlinear systems by using response-controlled stepped-sine testing. Mechanical Systems and Signal Processing 146 , 107023 (2021) https://doi.org/10. 1016/j.ymssp.2020.107023
- [133] B´ eliveau, J.-G., Vigneron, F.R., Soucy, Y., Draisey, S.: Modal parameter estimation from base excitation. Journal of Sound and Vibration 107 (3), 435-449 (1986) https:// doi.org/10.1016/S0022-460X(86)80117-1
- [134] Carrella, A.: Introduction to Environmental Testing. Siemens (2020). https://community.sw.siemens.com/s/ article/simcenter-testlab-vibration-control
- [135] Szempli´ska-Stupnicka, n W.: The modified single mode method in the investigations of the resonant vibrations of non-linear
- systems. Journal of Sound and Vibration 63 (4), 475-489 (1979) https://doi.org/10. 1016/0022-460X(79)90823-X
- [136] Karaa˘ ga¸lı, c T., ¨ zg¨ven, O u H.N.: Experimental quantification and validation of modal properties of geometrically nonlinear structures by using response-controlled stepped-sine testing. Experimental Mechanics 62 (2), 199-211 (2022) https://doi.org/ 10.1007/s11340-021-00784-9
- [137] G¨ urb¨z, u M.F., Karaa˘ ga¸lı, c T., ¨ zer, O M.B., ¨ zg¨ven, O u H.N.: Bypassing the repeatability issue in nonlinear experimental modal analysis of jointed structures by using the RCT-HFS framework. In: Society for Experimental Mechanics Annual Conference and Exposition, pp. 75-80. Springer, Cham (2024). https: //doi.org/10.1007/978-3-031-36999-5 10 . https://link.springer.com/10.1007/ 978-3-031-36999-5 { } 10
- [138] Koyuncu, A., Karaa˘ ga¸lı, c T., S ¸ ahin, M., ¨ zg¨ven, H.N.: Experimental modal analyO u sis of nonlinear amplified piezoelectric actuators by using response-controlled steppedsine testing. Experimental Mechanics 62 (9), 1579-1594 (2022) https://doi.org/10.1007/ s11340-022-00878-y
- [139] Karaa˘ ga¸lı, T., Nevzat Ozg¨ven, H.: Experc ¨ u imental modal analysis of structures with high nonlinear damping by using responsecontrolled stepped-sine testing. In: Brake, M.R., Renson, L., Kuether, R.J., Tiso, P. (eds.) Nonlinear Structures &amp; Systems, Volume 1. SEM 2023. Conference Proceedings of the Society for Experimental Mechanics Series., pp. 125-132. Springer, Cham (2024). Chap. 19. https: //doi.org/10.1007/978-3-031-36999-5 17 . https://link.springer.com/10.1007/ 978-3-031-36999-5 { } 17
- [140] Karaaga¸ clı, T.: Nonlinear system identification and nonlinear experimental modal analysis by using response controlled stepped sine testing. PhD thesis, Middle East Technical University (2020). https://hdl.handle.

## net/11511/89603

- [141] Jones, B.K., Trefan, G.: The Duffing oscillator: A precise electronic analog chaos demonstrator for the undergraduate laboratory. American Journal of Physics 69 (4), 464-469 (2001) https://doi.org/10.1119/1. 1336838
- [142] Srinivasan, K., Thamilmaran, K., Venkatesan, A.: Effect of nonsinusoidal periodic forces in Duffing oscillator: Numerical and analog simulation studies. Chaos, Solitons &amp; Fractals 40 (1), 319-330 (2009) https://doi. org/10.1016/j.chaos.2007.07.090
- [143] Raze, G.: An electronic Duffing oscillator (2024). https://github.com/GhislainRaze/ Electronic-Duffing
- [144] No¨ el, J.P., Kerschen, G.: Frequency-domain subspace identification for nonlinear mechanical systems. Mechanical Systems and Signal Processing 40 (2), 701-717 (2013) https://doi.org/10.1016/j.ymssp. 2013.06.034
- [145] Hippold, P., Scheel, M., Renson, L., Krack, M.: Robust and fast backbone tracking via phase-locked loops. Mechanical Systems and Signal Processing 220 (June), 111670 (2024) https://doi.org/10.1016/j. ymssp.2024.111670 arXiv:2403.06639
- [146] Jain, S., Breunung, T., Haller, G.: Fast computation of steady-state response for high-degree-of-freedom nonlinear systems. Nonlinear Dynamics 97 (1), 313-341 (2019) https://doi.org/10.1007/ s11071-019-04971-1
- [147] Robbins, E., Kuether, R.J., Pacini, B.R., Moreu, F.: Stabilizing a strongly nonlinear structure through shaker dynamics in fixed frequency voltage control tests. Mechanical Systems and Signal Processing 190 (December 2022), 110118 (2023) https: //doi.org/10.1016/j.ymssp.2023.110118
- [148] Zhou, T., Kerschen, G.: Identification of secondary resonances of nonlinear systems using phase-locked loop testing. Journal of

Sound and Vibration, 118549 (2024) https: //doi.org/10.1016/j.jsv.2024.118549

- [149] Peeters, B., Van der Auweraer, H., Guillaume, P., Leuridan, J.: The polymax frequency-domain method: A new standard for modal parameter estimation? Shock and Vibration 11 (3-4), 395-409 (2004) https:// doi.org/10.1155/2004/523692
- [150] Ewins, D.J.: Modal Testing: Theory, Practice and Application, (2009)
- [151] Preumont, A.: Vibration Control of Active Structures, 3rd edn. Solid Mechanics and Its Applications, vol. 179. Springer, Dordrecht (2011). https://doi.org/10.1007/978-94-007-2033-6 . http://link.springer.com/10.1007/ 978-94-007-2033-6
- [152] Quaegebeur, S., Raze, G., Cheng, L., Kerschen, G.: A virtual acoustic black hole on a cantilever beam. Journal of Sound and Vibration 554 (March), 117697 (2023) https://doi.org/10.1016/j.jsv.2023.117697 arXiv:2212.05939
- [153] Goldberg, P.W., Williams, C.K.I., Bishop, C.M.: Regression with input-dependent noise a gaussian process treatment. Advances in Neural Information Processing Systems, 493-499 (1998)
- [154] Li, Y., Dankowicz, H.: Adaptive control designs for control-based continuation of periodic orbits in a class of uncertain linear systems. Nonlinear Dynamics 103 (3), 2563-2579 (2021) https://doi.org/10.1007/ s11071-021-06216-6
- [155] Li, Y., Dankowicz, H.: Model-free continuation of periodic orbits in certain nonlinear systems using continuous-time adaptive control. Nonlinear Dynamics 111 (6), 4945-4957 (2023) https://doi.org/10.1007/ s11071-022-08059-1 arXiv:2203.10306
- [156] Rezaee, H., Renson, L.: Noninvasive adaptive control of a class of nonlinear systems with unknown parameters, 1-21 (2023) arXiv:2307.09806

- [157] Kruse, N., Wallner, H., Dittus, A., Lukas, B.: Large basins of attraction for controlbased continuation of unstable periodic states, 1-24 (2024)
- [158] Lim, J., Epureanu, B.I.: Forecasting a class of bifurcations: Theory and experiment. Physical Review E 83 (1), 016203 (2011) https://doi.org/10.1103/PhysRevE. 83.016203
- [159] Habib, G.: Predicting saddle-node bifurcations using transient dynamics: a model-free approach. Nonlinear Dynamics 111 (22), 20579-20596 (2023) https: //doi.org/10.1007/s11071-023-08941-6
- [160] Shen, J., Groh, R.M.J., Schenk, M., Pirrera, A.: Experimental path-following of equilibria using newton's method. part ii: Applications and outlook. International Journal of Solids and Structures 213 , 2540 (2021) https://doi.org/10.1016/j.ijsolstr. 2020.11.038
- [161] Magnevall, M., Josefsson, A., Ahlin, K.: Experimental verification of a control algorithm for nonlinear systems. Conference Proceedings of the Society for Experimental Mechanics Series (2006)
- [162] Novak, A., Simon, L., Lotton, P.: A simple predistortion technique for suppression of nonlinear effects in periodic signals generated by nonlinear transducers. Journal of Sound and Vibration 420 , 104-113 (2018) https://doi.org/10.1016/j.jsv.2018.01.038
- [163] Raze, G., Abeloos, G., Kerschen, G.: Experimental data from 'experimental continuation in nonlinear dynamics: recent advances and future challenges' (2024) https://doi. org/10.6084/m9.figshare.26412187

## Statements and Declarations

- · Funding: Ga¨ etan Abeloos and Ghislain Raze are a FRIA Grantee and a Postdoctoral Researcher, respectively, of the Fonds de la Recherche Scientifique - FNRS, which is gratefully acknowledged.
- · Competing interests: The authors have no relevant financial or non-financial interests to disclose.
- · Data Availability: The data collected and used in this research is available in [163].
- · Author Contributions: Ghislain Raze: literature search, formal analysis, data analysis, investigation and writing -original draft preparation. Ga¨ etan Abeloos: conceptualization, methodology and writing -review and editing. Ga¨ etan Kerschen: literature search, funding acquisition, supervision and writing -review and editing.
Ghislain Raze Gaëtan Abeloos Gaëtan Kerschen
Experimental continuation in nonlinear dynamics: recent advances and future challenges
Experimental continuation encompasses a set of methods that combine control
and continuation to obtain the full bifurcation diagram of a nonlinear system
experimentally, including responses that would be unstable in the system
without feedback control. Such control-based methods thus allow the
experimenter to directly and exhaustively explore the dynamics of the system
without the need for a good mathematical model. The objective of this paper is
twofold, namely (i) to review and present the state-of-the-art methods in a
unified manner and (ii) to introduce a novel experimental derivative-free
arclength continuation procedure, termed arclength control-based continuation.
These methods are also demonstrated and compared on an electronic Duffing
oscillator and a clamped thin plate featuring geometrical nonlinearity.
Finally, the current state of the art is reflected upon, and the challenges
lying ahead for this growing field are discussed.