- Division of Dynamics, Lodz University of Technology, Lodz, Poland
In this paper we discuss and explain the phenomenon of synchronization in lightly supported mechanical systems. The investigations are focused on the models of self–excited pendula hanged on the horizontally oscillating beam, which is lightly connected with the external support. Our results are based on the Centre-of-Mass (CoM) Theorem, which can be applied to the considered systems and allows to analytically confirm the observed behaviours. We present typical dynamical solutions, including periodic and quasiperiodic oscillations, within which the oscillators synchronize. The possible synchronous configurations are analyzed and examined, depending on the number of the pendula creating the system, their parameters and the initial conditions. We discuss bifurcations between different types of solutions, determining the regions and the conditions supporting the synchronization. Our investigations exhibit, that with the increase of the size of the network, the number of co–existing attractors also increases, leading to possible multistability and new types of behaviours (e.g., the traveling phase one). The results obtained numerically match with the analytical ones obtained from the CoM Theorem, which explains the existence of particular types of dynamical configurations. The study presented in this paper involves classical lightly supported pendula systems and due to their basic character, one can expect to observe similar behaviours in other types of mechanical models.
1 Introduction
Synchronization is one of the most fundamental types of behaviours found in nature [1]. The appearance of the coherent motion of dynamical systems has been reported in nonlinear vibrations [2], robotics [3], complex networks [4–6] or small–world systems [7], just to mention a few. In [8] one can find the study on the synchronization stability in mechanical networks with unstable local dynamics, while in [9, 10] the Authors investigate the phenomenon in the class of chaotic systems. Analysis of the noise effect on the synchronized oscillators and possible applications in Lagrangian systems have been discussed in [11, 12], respectively. The phenomenon is naturally related to control problems [13] and has been widely applied in signal processing, e.g. for secure communication [14]. A large part of synchronization problems refers to fundamental mechanical systems based on coupled pendula. The studies on synchronous dynamics have been performed for various pendula–type models, e.g. rotor-pendula [15], chaotic pendula [16] or Huygens’ coupling schemes [17], just to mention a few. In [18] Authors investigate the problem of controlled synchronization in coupled nodes, while in [19] one can find the discussion on pendula systems over digital communication channels.
The synchronous dynamics of nonlinear systems usually arises in various configurations, leading to the phenomenon of multistability [20, 21], when various attractors co–exist, depending on the system’s parameters and properties. The co–existence of solutions has been reported in many types of complex models, e.g. chaotic flows [22], mechanical oscillators [23] or brain dynamics [24]. Depending on the character of the considered system, such behaviour can be desired or not, which leads to the problem of control [25, 26] and basin stability [27]. In particular cases, one can even observe an expansion of co–existing stable solutions, i.e. extreme multistability, which has been found in both single [28] and coupled [29, 30] scenarios. The phenomenon is closely related to hidden oscillations [31], when the localization of possible attractors in the system’s phase space becomes not straightforward.
In this paper we investigate the dynamics and possible synchronous configurations for coupled pendula arranged in a lightly supported system. The considered model is schematically shown in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g001.jpg)
FIGURE 1. The model of n oscillating pendula with masses mi and lengths li suspended on the horizontally oscillating beam. The beam of mass M is connected with the support by the spring of stiffness k and the damper of damping coefficient c. Variable x stands for displacement of the beam, while φi is angular displacement of i–th pendulum. The paper uncovers various types of behaviours than can be found for different sizes of the network.
The system shown in Figure 1 consists of the beam of mass M [kg], which is connected with the support by the spring of stiffness k [N/m] and the damper of damping coefficient c [Ns/m]. The beam can oscillate in the horizontal direction and its position is denoted by variable x [m]. The considered structure supports the network of n pendula (nodes) with masses mi [kg] and lengths li [m], as shown in Figure 1 (i = 1, … , n). The pendula are equipped with the van der Pol type drives [32–34] (not shown in the figure) and their self–excited oscillations induce the motion of the beam, which allows to transfer the energy between the pendula. The angular displacement of the i–th node is given by variable φi ∈ (−π, π].
The dynamics of the system presented in Figure 1 can be investigated using the following equations of motion:
where i = 1, … , n.
During the analysis we have fixed the following parameters of system (1): M = 10 [kg] (the mass of the beam), k = 4 [N/m] (the stiffness of the spring), c = 1.53 [Ns/m] (the damping coefficient of the damper), g = 9.81 [m/s2] (the standard gravity term) and μ = 32.88 (the van der Pol type drive parameter). The damping of the i–th pendulum
The research presented in this paper continues the studies published in previous works [35–39]. In Refs. [35, 36] we have investigated model (1) with two nodes (n = 2), while in [37] the influence of the direction of the beam’s oscillations (horizontal, vertical and mixed) on the dynamics of two pendula has been studied. The clustering and experimental synchronization have been discussed in [38, 39], respectively.
In this paper we focus on a possible synchronization schemes for lightly supported systems, underlying the analytical solutions, which explain the observed behaviours and the results. We use the Centre-of-Mass (CoM) Theorem to show, that the performed analysis can be applied to any type of coupled oscillators arranged in a lightly supported model. The results uncovered in this paper expand the ones presented previously [35–39] and include new types of behaviours (e.g., the traveling phase solution), the appearance of the transient dynamics problems and possible multistability.
The analytical and the numerical investigations on system (1) are discussed in Section 2.
2 The results
When the beam in Figure 1 is lightly supported, i.e. k ≈ 0 [N/m] and c ≈ 0 [Ns/m], we can approximate system (1) as an isolated one. In such a case, the centre of mass of the model (the beam and the suspended pendula) is not moving, according to the Centre-of-Mass (CoM) Theorem. Consequently, the horizontal position of the centre of mass of the system shown in Figure 1 equals:
where
Excluding fixed coefficient
Assuming similar pendula lengths li ≈ l = const, i = 1, … , n, one gets:
The solutions of Eq. 4 depend on the number of pendula n and possible synchronous configurations. To compare the analytical results obtained using Eq. 4 with the numerical simulations, we have studied networks of different size, which is presented in the following Subsections.
To physically connect the beam with the external support (see Section 1. for details), we have analyzed model (1) with parameters k = 4 [N/m] and c = 1.53 [Ns/m]. Our analysis exhibits, that such scenario can be approximated by the CoM Theorem and Eq. 4 for isolated systems, as the numerical and the analytical results match with a good precision.
2.1 The case n = 2 pendula
In this Subsection we investigate system (1) with two pendula of identical masses, i.e. m1 = m2 = 1.0 [kg]. The lengths of the oscillators are fixed at l1 = 0.24848 [m] and l2 = 0.24849 [m]. A slight difference between the lengths is used to exclude coincidental synchronization of the nodes induced by the lack of numerical differentiability. According to formula
Two classical solutions of model (1) with n = 2 nodes are shown in Figure 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g002.jpg)
FIGURE 2. (Colour online). The in–phase synchronization of two pendula (variables φ1, φ2) oscillating in the anti–phase to the beam (variable x) (A) and the anti–phase synchronization of the nodes with the unmoving beam (C). The oscillations of the system’s centre of mass σ for case (A) are shown in panel (B). The results relate to classical synchronous behaviours found for coupled oscillators. System (1) with n = 2 (m1 = m2 = 1.0 [kg]).
Depending on the initial conditions, the pendula can stabilize on one of two possible synchronous configurations, i.e. the in–phase or the anti–phase solution.
In the first case, φ1 = φ2 = A sin(αt), where A [rad] and α [rad/s] denote the amplitude and the frequency of the oscillations, respectively, while t [s] is the dynamical time. For the in–phase scenario, the beam is moving, i.e.
Integrating Eq. 5 we obtain harmonic oscillations of the beam, with the same frequency α [rad/s] as the pendula and the anti–phase motion with respect to them. This scenario is presented in Figure 2A, where the displacements of the nodes are marked in red (φ1) and blue (φ2), while the movement of the beam is shown in black (x). The numerical time plots exhibit good match between the oscillations (φ1 = φ2) and allow to determine parameters A = 0.345 [rad] and α = 6.81 [rad/s]. The amplitude of the beam obtained from the simulations equals Xsim = 0.01412 [m], while the one calculated theoretically from Eq. 5 is given by
The case of two pendula with equal masses includes also the anti–phase synchronous configuration, when φ1 = A sin(αt), φ2 = A sin(αt + π) = −A sin(αt) and the beam is not moving (
Indeed, when the beam is not moving (
Assuming that the pendula have synchronized with phase β ∈ [0°, 180°], i.e. φ1 = A sin(αt) and φ2 = A sin(αt + β), one gets from Eq. 6:
System (7) has two solutions: (i) β = 0°, when m1 = −m2 < 0, which is physically contradictory and (ii) β = 180°, when m1 = m2, which is exactly the scenario of the anti–phase synchronization shown in Figure 2C. In this case A = 0.349 [rad] and α = 6.24 [rad/s].
When the masses of the pendula in model (1) are different (m1 ≠ m2), one can observe also the in–phase synchronization patterns with different amplitudes and quasiperiodic motion, depending on the parameters of the system. The discussion of such scenarios can be found in [40].
2.2 The case n = 3 pendula
In this Subsection we investigate model (1) with three pendula suspended on the beam (n = 3). The lengths of the units are fixed at l1 = 0.24848 [m], l2 = 0.24849 [m] and l3 = 0.2485 [m], while the masses m1, m2, m3 have been varied, depending on the considered synchronous scenario.
Due to the fact, that the lengths of the pendula are almost equal, the natural periods of the oscillators remain similar and it becomes natural to assume, that the pendula can synchronize with equal amplitudes and two phases locked, i.e. φ1 = A sin(αt), φ2 = A sin(αt + β2), and φ3 = A sin(αt + β3), where β2, β3 ∈ [0°, 360°) are the phases between the nodes (β2–the phase between the 1st and the 2nd node, β3–the phase between the 1st and the 3rd node). If
System (8) consists of two equations with two variables β2, β3 and allows to determine the synchronous configurations (the values of the phases) when the distribution of masses m1, m2, m3 is given. The examples of the behaviours observed in the considered case are shown in Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g003.jpg)
FIGURE 3. (Colour online). The periodic synchronous configurations for the identical masses (A) m1 = m2 = m3 = 1.0 [kg] and non–identical masses (B) m1 = 1.0, m2 = 0.7, m3 = 0.5 [kg]. The scenario of quasiperiodic synchronization between the 2nd (blue) and the 3rd (green) node is shown in panel (C) for m1 = 1.0, m2 = 0.5, m3 = 0.2 [kg]. In the case of three oscillators, both regular and irregular dynamics can arise, leading to more complex synchronous scenarios. System (1) with n = 3 pendula.
The scenario of identical masses m1 = m2 = m3 = 1.0 [kg] can be found in Figure 3A. The motion of the pendula has been marked in red (1st), blue (2nd) and green (3rd), while the small oscillations of the beam are shown in black (the latter time plot is multiplied 100 times–100x). Assuming that the velocity of the beam’s oscillations is small compared to the pendula, we can apply system (8) for finding the attractors. Indeed, Eq. 8 for m1 = m2 = m3 = 1.0 [kg] have two solutions: (i) β2 = 120°, β3 = 240°and (ii) β2 = 240°, β3 = 120°, which indicate the same state, since the masses of the pendula are equal. The example of the phase–locked synchronization pattern shown in Figure 3A corresponds to analytical result (i), with the phases calculated numerically at β2 = 120.21°and β3 = 240.68°. The centre of mass of the system is almost stationary, with the amplitude of its oscillations equal to Aσ = 3.17 ⋅ 10–6 [m].
When the pendula are non–identical, the values of phases β2, β3 strictly depend on the distribution of the masses and possible solutions of system (8). The example of the phase–locked state for m1 = 1.0, m2 = 0.7 and m3 = 0.5 [kg] is presented in Figure 3B, where the numerical values of the phases equal β2 = 207.63°, β3 = 139.08°. The analytical solution that can be obtained using Eq. 8 matches with the numerical one, exhibiting phases β2 = 207.72°and β3 = 139.32°. In this case, the amplitude of the centre of mass does not exceed Aσ = 2.62 ⋅ 10–6 [m]. It should be noted, that one can observe also the second synchronous configuration, for which β2 = 152.37°and β3 = 220.92°. The appearance of the pair of phase–locked states is associated with the following property of system (8): if
The pendula of model (1) can also synchronize in the in–phase, similar to n = 2 scheme (see Section 2.1. for details). In this case, the applicability of system (8) is excluded, since the beam is oscillating against the pendula (to keep the centre of mass stationary) and consequently
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g004.jpg)
FIGURE 4. (Colour online). The basins of attraction for different distribution of the pendula masses. In (A) the scenario of three co–existing periodic synchronous states is shown (IN: in–phase state, PL: the first phase–locked state and PL*: the second phase–locked state), while in (B) the case of bi–stability of the in–phase (IN) and the quasiperiodic (Q) configurations is presented. The values of the masses in panels (A) and (B) correspond to the ones considered in Figure 3B and Figure 3C, respectively. The results show, that depending on the initial conditions, various types of behaviours co–exist.
The results presented in Figure 4A correspond to the case of m1 = 1.0, m2 = 0.7, m3 = 0.5 [kg] (see the time plots in Figure 3B) and uncover the basins of attraction of three possible solutions: 1) the in–phase synchronization (IN–yellow), 2) the first phase–locked configuration with β2 = 152.37°, β3 = 220.92°(marked as PL–dark blue) and 3) the second phase–locked configuration with β2 = 207.72°, β3 = 139.32°(marked as PL*–light blue; to distinguish that the observed phase–locked solutions are different in the sense of β2, β3 phases distribution, we denote the second state by asterisk). The basins in Figure 4A are shown in subspace (φ2, φ3) ∈ (−π,π]2 (the initial positions of the 2nd and the 3rd pendulum), with the remaining conditions fixed at zeros, i.e.
As we have shown, the analytical solutions of system (8) for the chosen distribution of the pendula masses indicate possible attractors of model (1) with a good precision (differences between the analytical and the numerical phases β2, β3). Moreover, it can be easily shown, that for particular masses mi, i = 1, 2, 3, system (8) may not have any solution, which suggests the disappearance of phase–locked configurations. The example of such scenario is discussed in Figure 3C and Figure 4B for m1 = 1.0, m2 = 0.5 and m3 = 0.2 [kg]. Due to the fact, that Eq. 8 have no solution for the chosen masses, the oscillations of the pendula shown in Figure 3C become quasiperiodic and grouped into two clusters: 1) the 1st pendulum (red; φ1) and 2) the quasiperiodic, synchronous motion of the 2nd and the 3rd node (blue and green; φ2 = φ3). The quasiperiodic character of the vibrations has been confirmed using Poincare sections, which are included in the inbox in Figure 3C. The 1st (red) pendulum has been chosen as the reference node for creating the maps, i.e. the points have been collected every time it reached the local maximum. The observed quasiperiodic state co–exists with the in–phase synchronization and the basins of these two possible behaviours are shown in Figure 4B (the in–phase state–I in yellow; the quasiperiodic state–Q in brown). It should be noted, that the almost stationary position of the centre of mass in the case of quasiperiodic oscillations is preserved, with the amplitude Aσ = 1.2 ⋅ 10–4 [m].
To investigate possible transitions between different synchronous configurations of model (1) with n = 3 pendula, we have performed the bifurcation analysis in relation to system (8). The example of the bifurcation scenario for fixed masses m1 = 1.0, m2 = 0.5 [kg] and varied m3 is shown in Figure 5.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g005.jpg)
FIGURE 5. (Colour online). The bifurcation diagrams of model (1) for fixed m1 = 1.0, m2 = 0.5 [kg] and varied m3 ∈ [0.16, 2.92] [kg] (A). The interval is split into seven cases (I–VII), depending on the behaviour of the pendula and possible solutions of system (8). The bifurcations of phases β2 and β3 corresponding to the phase–locked synchronous configurations during stage IV are shown in panel (B). Using bifurcation analysis, one can trace transitions occurring in the system and the arise of different types of states.
The bifurcation diagrams shown in Figure 5A are split into seven regions: I–VII, depending on the dynamics of the pendula and possible synchronous configurations. The procedure applied for collecting the points for the diagrams is the same as the one used for the Poincare section shown in the inbox in Figure 3C (the 1st pendulum is chosen as the reference one).
The analysis begins at m3 = 0.2 [kg], when the behaviour of the system is quasiperiodic as shown in Figure 3C. This type of dynamics can be observed in a narrow region labelled as II for 0.1885 < m3 < 0.2158 [kg], which has been enlarged in the inbox in Figure 5A. Decreasing the mass of the 3rd pendulum, at m3 = 0.1885 [kg] the system bifurcates into region I, where all three pendula synchronize in the in–phase. Further decrease of parameter m3 does not change the dynamics of the model, since the synchronized pendula act as one clustered node, which oscillates in the anti–phase to the beam of a fixed mass. In both regions I and II, system (8) has no solutions.
Bifurcating model (1) in the opposite direction, i.e. starting from case II for m3 = 0.2 [kg] and increasing parameter m3, one can observe the vanishing of the quasiperiodic solution, which finally bifurcates into the periodic one for m3 = 0.2158 [kg] (region III). In this case, the 2nd (blue) and the 3rd (green) pendula are still synchronized (forming a cluster) but the amplitude of their oscillations is greater than the amplitude of the 1st (red) node. This scenario corresponds to the case of two pendula with different masses, which synchronize with various amplitudes [35, 40]. The reported dynamics resides in region III for 0.2158 ≤ m3 < 0.5 [kg], where the relation m1 > m2 + m3 is satisfied. With the increase of mass m3, the amplitude of the clustered pendula (blue and green) also increases and finally reaches the amplitude of the red node at m3 = 0.5 [kg], when m1 = m2 + m3. The system bifurcates into stage IV, where all of the pendula are oscillating with the same amplitude and system (8) has two analytical solutions leading to the phase–locked synchronous states. The phases between the pendula (β2, β3) observed in this scenario are shown in Figure 5B for 0.5 ≤ m3 ≤ 1.5 [kg]. In the beginning of stage IV, for m3 = 0.5 [kg] the phases are equal: β2 = β3 = 180°; this corresponds to the anti–phase pattern between the clustered nodes (blue and green) and the solitary one (red), since m1 = m2 + m3. With the increase of mass m3, the cluster breaks as the 2nd pendulum begins to converge to the 1st one (β2 decreases as shown in Figure 5B). The system resides on the phase–locked attractor with the phases depending on the value of m3. The convergence of the blue node to the red one is continued in a wide range of parameter m3, until the critical point at m3 = 1.5 [kg], for which β2 = 0°and β3 = 180°. In this scenario the 1st (red) and the 2nd (blue) pendula clustered, leaving the 3rd (green) one solitary. The synchronization becomes anti–phase, according to equality m1 + m2 = m3.
With further increase of mass m3, the system bifurcates into stage V, where Eq. 8 have no solutions and the pendula oscillate with different amplitudes. For 1.5 < m3 < 2.732 [kg] (region V), the amplitude of the clustered nodes (red and blue) is higher than the one observed for the 3rd (green) pendulum, as m1 + m2 < m3. The scenario changes from periodic to quasiperiodic motion for m3 = 2.732 [kg], when the dynamics reaches region VI. In this case, the clustered character of the pendula is preserved, as shown in the diagram in Figure 5A. With further increase of mass m3 [kg], the model finally stabilizes on the in–phase synchronization pattern for m3 ≥ 2.887 [kg], which has been labelled as region VII.
Scenarios I and VII are qualitatively the same. In scenarios (II and VI) and (III and V) one can observe the clustering of the pendula and their quasiperiodic or periodic dynamics, respectively. The clustering configuration depends on the distribution of masses m1, m2, m3 and one can obtain different patterns varying these parameters. The remaining case IV allows to observe phase–locked synchronous motion, which can be described by the analysis of the centre of mass, i.e. system (8).
To examine possible phase–locking patterns between the pendula, we have investigated Eq. 8 for fixed mass m1 = 1.0 [kg] and varied m2, m3. Our results are presented in Figure 6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g006.jpg)
FIGURE 6. (Colour online). The values of the phases within the phase–locked synchronization patterns for fixed m1 = 1.0 [kg] and varied m2, m3: values of β2 and β3 are shown in panels (A) and (B), respectively. The possible phases β2, β3, which can be realized by model (1) due to the solvability of system (8) are shown in (C). The results determine regions, where particular types of behaviours can be found, as well as their properties.
In Figures 6A,B one can see the values of phases β2 (panel a) and β3 (panel b) obtained analytically using system (8) for m2, m3 ∈ [0, 3] [kg] (mass m1 = 1.0 [kg] fixed). The results have been limited to the case, when β2 ∈ [0, 180°] (as described above, system (8) has a pair of solutions, i.e. (β2, β3) and (360° − β2, 360° − β3)). Using the maps presented in Figures 6A,B one can identify the regions of the pendula masses, for which phase–locked synchronization is possible. We have also calculated analogous maps by the numerical examination of model (1) for m2, m3 ∈ [0, 3] [kg] and the obtained results (phases) match with the ones observed for the analytical case. It should be noted, that the regions left blank in Figures 6A,B correspond to the in–phase synchronization and/or quasiperiodic motion, depending on the parameters.
The relation between the pendula masses mi, i = 1, 2, 3 and phases β2, β3 can be also inverted, i.e. one can investigate what values of the masses should be selected for the desired phase–locked configuration (particular values of the phases). Indeed, fixing the mass of a selected pendulum, e.g. m1 and considering β2, β3 as the parameters, we can calculate the values of the remaining masses m2 and m3 from Eq. 8 in the following way:
Using Eq. 9 one can calculate the values of m2 and m3, which allow to obtain the desired (β2, β3) phase–locked synchronous configuration. The solutions of (9) have the physical sense only when both formulas are greater than zero; the regions of (β2, β3) values corresponding to such scenario are shown in Figure 6C in black. The presented map allows to determine which values of the phases (the ones corresponding to the black domain) can be realized by system (1) within the phase–locked state. It should be noted, that close to the borders of the region shown in Figure 6C, masses m2, m3 in Eq. 9 converge to zero or infinity, which may influence the dynamics of system (1) and lead to unexpected behaviours.
2.3 The case n = 4 pendula
In this Subsection we study the case of model (1) with four pendula suspended on the beam (n = 4). The lengths of the oscillators are chosen as l1 = 0.24848 [m], l2 = 0.24849 [m], l3 = 0.2485 [m] and l4 = 0.24851 [m]. If the pendula synchronize with equal amplitudes and phases locked, i.e. φi = A sin(αt + βi), i = 1, 2, 3, 4, where β1 = 0°and β2, β3, β4 ∈ [0°, 360°), then for
Depending on the values of masses mi, i = 1, 2, 3, 4, system (10) can have no solutions or infinitely many solutions (two equations with three variables β2, β3, β4). The considered scenario is more complex than the one observed for n = 3 pendula, where the number of solutions of Eq. 8 has been limited (see Section 2.2. for details).
To investigate model (1) with four nodes, we have varied pendula masses in different scenarios. The examples of possible dynamical patterns are discussed in Figure 7.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g007.jpg)
FIGURE 7. (Colour online). Possible behaviours of the pendula for (A) m1 = m2 = m3 = m4 = 1.0 [kg] (identical masses and the traveling phase state), (B) m1 = 1.0, m2 = 0.99, m3 = 1.01, m4 = 0.985 [kg] (slightly different masses and the phase–locked state) and (C) m1 = 1.0, m2 = 0.7, m3 = 0.5, m4 = 0.45 [kg] (different masses and the clustering of the pendula). System (1) with n = 4 nodes. In the case of four coupled pendula, one can identify new types of states and even more complex synchronous scenarios.
We have investigated three main schemes of the masses distribution, i.e. identical masses, slightly different masses and different masses, which are discussed below.
If the pendula are identical, i.e. mi = const, i = 1, 2, 3, 4, then system (10) can be simplified by dividing the equations by the masses. It can be easily shown, that in such scenario, the solutions of (10) are given as follows: β2 = 180°,
The example of the described scenario is presented in Figure 7A for m1 = m2 = m3 = m4 = 1.0 [kg], where the time plots of the nodes are marked in red (1st), blue (2nd), green (3rd), and grey (4th), while the beam is shown in black. Since the oscillations within the clusters are anti–phase (red vs. blue and green vs. grey), the beam is not moving (as have been described previously for n = 2 pendula in Figure 2C; see Section 2.1. for details). As we have observed, the phase between the clusters (parameter
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g008.jpg)
FIGURE 8. (Colour online). In (A) the traveling of phases β2, β3 and β4 for identical pendula mi = 1.0 [kg], i = 1, 2, 3, 4 is shown. The basins of attraction of possible dynamical patterns (IN and TP) for varied (φ3, φ4) ∈ (−π,π]2 and (φ2, φ4) ∈ (−π,π]2 are presented in panels (B) and (C), respectively. The arise of the traveling phase state is related to the similarities in the pendula properties (the masses).
The diagrams included in Figure 8A present the values of β2, β3 and β4 phases in a long time interval of model (1) simulations. As one can see, the numerical solution matches with the analytical one (β2 = 180°,
As we have observed during various simulations, the configuration of the clusters can be correlated with the lengths of the pendula, i.e. the slight differences in lengths li and their distribution within the nodes determine the pairs, which will cluster. For the considered l1 = 0.24848 [m], l2 = 0.24849 [m], l3 = 0.2485 [m] and l4 = 0.24851 [m], system (1) has always finally converged to the scenario shown in Figure 8A and the observed clustering has been independent of the initial conditions. The described property has not been observed when varying the conditions from and between the clusters, which has been shown in Figures 8B,C for the basins of attraction. The scenarios of varied (φ3, φ4) ∈ (−π,π]2 and (φ2, φ4) ∈ (−π,π]2 (the remaining conditions equal to zeros) are presented in Figure 8B and Figure 8C, respectively and correspond to the initial conditions of the pendula forming the cluster (b) and the initial conditions of the pendula from different clusters (c). As one can see, the “traveling phase” state (TP, brown) co–exists with the in–phase synchronization pattern (IN, yellow) but the regions of convergence for both scenarios (b) and (c) are qualitatively indifferent.
The traveling character of the phases between the pendula is possible due to the identity of the masses. When the masses become slightly different, system (10) cannot be simplified as previously but its solutions are close to the ones described above (two clusters), since m1 ≈ m2 ≈ m3 ≈ m4. The example of such scenario is shown in Figure 7B for m1 = 1.0, m2 = 0.99, m3 = 1.01, m4 = 0.985 [kg]. As one can see, the 1st (red) and the 4th (grey) pendulum are synchronized almost in the anti–phase (the first cluster) and the 2nd (blue) and the 3rd (green) pendulum are also synchronized close to the anti–phase state (the second cluster). The slight shifts in the oscillations of the pendula within the clusters are shown in the enlargement in Figure 7B. Since the motion is not precisely anti–phased, the beam exhibits slight oscillations, which magnitude has been increased 1000 times for better clarity (1000x in the figure).
The slight differences in the masses of the pendula lead to the disappearance of the TP state and stabilization of the nodes on particular phase–locked synchronous configuration. The example of the evolution of phases β2, β3 and β4 for such scenario is presented in Figure 9A.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-g009.jpg)
FIGURE 9. (Colour online). The evolution of phases β2, β3, β4 (left panel) and possible synchronous configurations in (β2, β3) plane (right panel). Results (A,B) correspond to the case of slightly different masses: m1 = 1.0, m2 = 0.99, m3 = 1.01, m4 = 0.985 [kg], while results (C,D) have been obtained for different masses m1 = 1.0, m2 = 0.7, m3 = 0.5, m4 = 0.45 [kg]. The results exhibit, that the distribution of the pendula masses plays a major role in possible behaviours of the system.
The points for the diagrams included in Figure 9A have been collected after 1000 T [s], which has been considered as possible transient time for the stabilization of the system. As one can see, the time necessary for the full stabilization has reached around 500 000 T [s], as the phases are changing for a long time interval (time below 250 000 T [s] in Figure 9A). Our analysis for various initial conditions exhibits, that for slightly different masses of the pendula, system (1) finally converges to particular phase–locked state but the transient time of such stabilization can become very long (even exceeding 500 000 T [s]). These results uncover possible transient dynamics problems [41–43], which can arise when studying complex networks and become misleading for the researcher trying to determine the attractor(s) of the considered model.
As we have described above, system (10) can have infinitely many solutions for the fixed masses, leading to the possibility of the infinitely many phase–locked configurations. To investigate this property for the considered slightly different masses, we have performed a series of trials, simulating model (1) from random initial conditions and determining the final synchronous configuration. The results of our analysis are discussed in Figure 9B and are based on 1000 trials in total.
For each numerical trial we have simulated system (1) for 500 000 T [s], after which the observed pattern has been examined for stability (if the phases β2, β3 and β4 are still changing or not). In the case of the stable phases, the obtained phase–locked configuration has been included in Figure 9B in the form of the red dot, which coordinates indicate values β2 (horizontal axis) and β3 (vertical axis), while the dot’s label (in grey text) equals the value of phase β4. As one can see, the points marked in Figure 9B exhibit linear relation (the dashed red line along the points) and the values of the obtained phases are close to the solution, which can be calculated from system (10) for the identical masses, i.e. β4 ≈ 180°,
To prevent the overlapping of the states, when the phases β2, β3 and β4 for two particular trials are slightly different (which may be caused by numerical inaccuracies), we have applied the following criterion for distinguishing the states: “Configuration 1” with phases
In the last considered scenario we have studied model (1) with different masses, fixing m1 = 1.0, m2 = 0.7, m3 = 0.5 and m4 = 0.45 [kg]. The phase–locked pattern obtained for the given distribution is shown in Figure 7C, where the 2nd (blue) and the 3rd (green) pendulum formed the cluster (φ2 ≈ φ3) and the scenario is similar to the one observed for n = 3 coupled nodes (see Section 2.2. for details). Indeed, when the masses of the pendula are different, the system (apart from the in–phase synchronization) tends to the clustering of the nodes [44–47]. This property has been also observed for larger networks, when the oscillators group into three clusters [38].
In the case of clustering, possible phases between the pendula can be easily obtained using CoM Theorem and system (10). Indeed, when two of the pendula form a cluster, e.g. the 3rd and the 4th, we obtain φ3 = φ4, β3 = β4 and the simplification of Eq. 10 in the following form:
where β34 = β3 = β4. The above system is analogous to system (8) for n = 3 pendula and possesses solutions depending on the values of m1, m2 and (m3 + m4).
For four pendula one can obtain six different configurations of clustering (m1 + m2, m1 + m3, m1 + m4, m2 + m3, m2 + m4 and m3 + m4), of which three lead to solutions of system (10) for the considered masses m1 = 1.0, m2 = 0.7, m3 = 0.5 and m4 = 0.45 [kg]. Each clustering configuration involves a pair of phase–locked synchronous solutions (see Section 2.2. for details), which results in six different synchronization patterns in total. The latter ones have been listed as “phase–locking (1)–(6)” in Table 1, where in the first column the labels of the states are shown (attractor type), while in the second one the values of the phases obtained analytically from (10) are presented (analytical phases). The table includes also the in–phase synchronization pattern and the “phase–locking (*)” state, which is discussed further.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021836/fphy-10-1021836-HTML/image_m/fphy-10-1021836-t001.jpg)
TABLE 1. The comparison of the analytical and the numerical phases for the phase–locked synchronous configurations (1)–(6). The probabilities of the occurrence of possible states are given in the last column. System (1) with different masses: m1 = 1.0, m2 = 0.7, m3 = 0.5, m4 = 0.45 [kg].
The example of the phases convergence for “phase–locking (2)” state is included in Figure 9C and shows the evolution of phases β2, β3 and β4 for 50 000 T [s] (the points collected after 1000 T [s]). In this case, the phases stabilize faster than for the slightly different masses scenario (compare Figure 9A and Figure 9C) and the transient motion becomes easier to omit. The phases of the 2nd (blue, β2) and the 3rd (green, β3) node are very similar, which exhibits a good match between the analytical results and the ones obtained from the numerical simulations. The comparison of both cases can be found in Table 1 in the second (analytical phases) and the third (numerical phases) columns.
As one can see in Table 1, two patterns “phase–locking (1)” and “phase–locking (3)” are labelled as unstable. Even though system (10) can have particular solution (e.g., the “phase–locking (1)”state), its stability has to be examined numerically by simulating the model. To investigate the attractors of system (1) for different masses, we have performed a series of 40 000 trials, determining the synchronous configurations after 10 000 T [s]. The obtained results are discussed in Figure 9D in (β2, β3) plane, with the values of β4 labeling the points (see also the description of Figure 9B above). Four black dots correspond to solutions “phase–locking (2), (4), (5) and (6)” described in Table 1, while the single red dot indicates the solution “phase–locking (*)” (the last row in Table 1), for which the pendula have not clustered. During the simulations we have not observed the stability of patterns “phase–locking (1)” and “phase–locking (3)”, which suggest that these solutions are unstable. It should be noted though, that the phases obtained for the “phase–locking (*)” state are close to the analytical ones for “phase–locking (3)”, which may indicate that these configurations are closely related and the difference between the analytics and the numerics results from the assumptions of the considered model (not completely isolated).
The probability analysis performed for the solutions included in Table 1 shows, that the in–phase synchronous configuration is the most probable state (probability around 43%), while the phase–locked patterns appear with similar probabilities ranging from 5% to 15%. The presented results are based on a series of 40 000 trials as mentioned above; for each trial, the initial positions of the pendula have been sampled from the uniform distribution U[ − π, π), with the remaining initial conditions of the system considered zeros.
It should be noted, that the considered scenario of four coupled pendula includes also other possible behaviours, such as quasiperiodic dynamics or synchronous patterns with different amplitudes. When the pendula group into three clusters (as shown in Figure 7C), the scenario is analogous to the synchronization of three oscillators (n = 3 in Section 2.2.) and we can perform similar bifurcation analysis as the one presented in Figure 5.
Generalizing the results, when model (1) consists of n > 4 pendula synchronized within the phase–locked state, i.e. φi = A sin(αt + βi), i = 1, … , n, β1 = 0°and β2, … , βn ∈ [0°, 360°), system (4) transforms into the following equations:
System (12) can possess infinitely many solutions and generalizes system (10) considered for n = 4 nodes. In this case, one can expect to obtain similar synchronous configurations and dynamics, depending on the parameters of the pendula and the beam, as well as the initial conditions. The increased number of nodes can lead to the appearance of more complex behaviours, including chimeras [48, 49] and chimera–like states [50, 51], which occur naturally in the networks of coupled oscillators.
3 Conclusion
In this paper we have investigated the dynamics of coupled self–excited pendula, arranged in a lightly supported system. Depending on the network’s size (the number of the nodes suspended on the beam), one can observe different types of synchronous configurations and behaviours. Our results are based on the Centre-of-Mass (CoM) Theorem, which can be applied to the considered, almost isolated system. We have described typical dynamical structures, including the in–phase and the anti–phase synchronization, as well as more complex patterns like the clustering of the pendula or the phase–locked solutions. The basins of attraction included in the study uncover possible co–existence configurations of the system’s states, while the bifurcation analysis allows to trace the transitions between different synchronous solutions. The analysis can be further applied for determining the relations between the pendula masses and their phases within the phase–locking scenarios, which has been included in the results. During the research we have uncovered the “traveling phase” state, which is characterized by the continuous change of the phases between the synchronized oscillators and transient dynamics problems, leading to the uncertainty about the stabilization of the observed synchronous states. The analysis of the considered models has exhibited possible high multistability with many co–existing attractors, that can be observed for slightly different nodes.
The analytical results obtained from the CoM Theorem match with the numerical ones calculated during simulations. The former results explain and allow to understand the dynamics that is observed within the models, especially the synchronous configurations between the pendula. As we have shown, the theory approximates lightly supported systems with a good precision and can be applied for studying similar networks of the mechanical type. The latter can be especially beneficial in the case of experimental investigations, since the structure of presented model allows to confirm the results practically with the use of simple mechanical setup.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This work has been supported by the National Science Centre, Poland: SONATA Programme (Project No 2019/35/D/ST8/00412) and OPUS Programmes (Project No 2018/29/B/ST8/00457 and Project No 2021/43/B/ST8/00641).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Pikovsky A, Rosenblum M, Kurths J. Synchronization: A universal concept in nonlinear sciences. Cambridge University Press (2003).
2. Blekhman II. Synchronization in science and technology. American Society of Mechanical Engineers (ASME) (1988).
3. Nijmeijer H, Rodriguez-Angeles A. Synchronization of mechanical systems. World Scientific (2003).
4. Awerbuch B. Complexity of network synchronization. J ACM (1985) 32(4):804–23. doi:10.1145/4221.4227
5. Arenas A, Diaz-Guilera A, Kurths J, Moreno Y, Zhou C. Synchronization in complex networks. Phys Rep (2008) 469(3):93–153. doi:10.1016/j.physrep.2008.09.002
6. Ghosh D, Frasca M, Rizzo A, Majhi S, Rakshit S, Alfaro-Bittner K, et al. The synchronized dynamics of time-varying networks. Phys Rep (2022) 949:1–63. doi:10.1016/j.physrep.2021.10.006
7. Barahona M, Pecora LM. Synchronization in small-world systems. Phys Rev Lett (2002) 89:054101. doi:10.1103/physrevlett.89.054101
8. Nair S, Leonard NE. Stable synchronization of mechanical system networks. SIAM J Control Optim (2008) 47(2):661–83. doi:10.1137/050646639
9. Mahmoud GM, Aly SA, Al-Kashif MA. Dynamical properties and chaos synchronization of a new chaotic complex nonlinear system. Nonlinear Dyn (2008) 51:171–81. doi:10.1007/s11071-007-9200-y
10. Boccaletti S, Kurths J, Osipov G, Valladares DL, Zhou CS. The synchronization of chaotic systems. Phys Rep (2002) 366(1):1–101. doi:10.1016/s0370-1573(02)00137-0
11. Ver Hoeye S, Suarez A, Sancho S. Analysis of noise effects on the nonlinear dynamics of synchronized oscillators. IEEE Microw Wirel Compon Lett (2001) 11(9):376–8. doi:10.1109/7260.950766
12. Spong MW, Chopra N. Synchronization of networked Lagrangian systems. In: Lagrangian and Hamiltonian methods for nonlinear control 2006. Springer Berlin Heidelberg (2007). p. 47.
13. Borhaug E. Nonlinear control and synchronization of mechanical systems (Doctoral thesis). Norway: Norwegian university of science and technology (2008).
14. De Sousa Vieira M, Lichtenberg AJ, Lieberman MA. Nonlinear dynamics of self-synchronizing systems. Int J Bifurcation Chaos (1991) 1(3):691–9. doi:10.1142/s0218127491000506
15. Fang P, Hou Y. Synchronization characteristics of a rotor-pendula system in multiple coupling resonant systems. Proc Inst Mech Eng C: J Mech Eng Sci (2018) 232(10):1802–22. doi:10.1177/0954406217711468
16. Baker GL, Blackburn JA, Smith HJT. Intermittent synchronization in a pair of coupled chaotic pendula. Phys Rev Lett (1998) 81:554–7. doi:10.1103/physrevlett.81.554
17. Ramirez JP, Alvarez J. Rotating waves in oscillators with Huygens’ coupling**This work was partly supported by the CONACyT under Grant CB2012-180011-Y. IFAC-PapersOnLine (2015) 48(18):71–6. doi:10.1016/j.ifacol.2015.11.013
18. Pogromsky AY, Belykh VN, Nijmeijer H. Controlled synchronization of pendula. In: 42nd IEEE International Conference on Decision and Control (IEEE Cat. No.03CH37475), 5 (2003). p. 4381.
19. Fradkov AL, Andrievsky B, Ananyevskiy M. State estimation and synchronization of pendula systems over digital communication channels. Eur Phys J Spec Top (2014) 223(4):773–93. doi:10.1140/epjst/e2014-02140-0
21. Feudel U. Complex dynamics in multistable systems. Int J Bifurcation Chaos (2008) 18(06):1607–26. doi:10.1142/s0218127408021233
22. Li C, Hu W, Sprott JC, Wang X. Multistability in symmetric chaotic systems. Eur Phys J Spec Top (2015) 224:1493–506. doi:10.1140/epjst/e2015-02475-x
23. Brzeski P, Perlikowski P. Sample-based methods of analysis for multistable dynamical systems. Arch Comput Methods Eng (2019) 26:1515–45. doi:10.1007/s11831-018-9280-5
24. Kelso JAS. Multistability and metastability: Understanding dynamic coordination in the brain. Phil Trans R Soc B (2012) 367(1591):906–18. doi:10.1098/rstb.2011.0351
25. Pisarchik AN, Feudel U. Control of multistability. Phys Rep (2014) 540(4):167–218. doi:10.1016/j.physrep.2014.02.007
26. Feudel U, Grebogi C. Multistability and the control of complexity. Chaos (1997) 7(4):597–604. doi:10.1063/1.166259
27. Rakshit S, Bera BK, Majhi S, Hens C, Ghosh D. Basin stability measure of different steady states in coupled oscillators. Sci Rep (2017) 7:45909. doi:10.1038/srep45909
28. Jafari S, Ahmadi A, Khalaf AJM, Abdolmohammadi HR, Pham V-T, Alsaadi FE. A new hidden chaotic attractor with extreme multi-stability. AEU - Int J Elect Commun (2018) 89:131–5. doi:10.1016/j.aeue.2018.03.037
29. Chakraborty P, Poria S. Extreme multistable synchronisation in coupled dynamical systems. Pramana - J Phys (2019) 93:19. doi:10.1007/s12043-019-1789-0
30. Hens CR, Banerjee R, Feudel U, Dana SK. How to obtain extreme multistability in coupled dynamical systems. Phys Rev E (2012) 85:035202. doi:10.1103/physreve.85.035202
31. Kuznetsov NV, Leonov GA. Hidden attractors in dynamical systems: Systems with no equilibria, multistability and coexisting attractors. IFAC Proc Volumes (2014) 47(3):5445–54. doi:10.3182/20140824-6-za-1003.02501
32. Holmes PJ, Rand DA. Bifurcations of the forced van der Pol oscillator. Q Appl Math (1978) 35:495–509. doi:10.1090/qam/492551
33. Kennedy M, Chua L. Van der Pol and chaos. IEEE Trans Circuits Syst (1986) 33(10):974–80. doi:10.1109/tcs.1986.1085855
34. Guckenheimer J. Dynamics of the van der Pol equation. IEEE Trans Circuits Syst (1980) 27(11):983–9. doi:10.1109/tcs.1980.1084738
35. Czołczyński K, Perlikowski P, Stefański A, Kapitaniak T. Why two clocks synchronize: Energy balance of the synchronized clocks. Chaos (2011) 21(2):023129. doi:10.1063/1.3602225
36. Dudkowski D, Czołczyński K, Kapitaniak T. Multistability and basin stability in coupled pendulum clocks. Chaos (2019) 29(10):103140. doi:10.1063/1.5118726
37. Dudkowski D, Czołczyński K, Kapitaniak T. Multistability and synchronization: The co-existence of synchronous patterns in coupled pendula. Mech Syst Signal Process (2022) 166:108446. doi:10.1016/j.ymssp.2021.108446
38. Czołczyński K, Perlikowski P, Stefanski A, Kapitaniak T. Clustering of non-identical clocks. Prog Theor Phys (2011) 125(3):473–90. doi:10.1143/ptp.125.473
39. Kapitaniak M, Czołczyński K, Perlikowski P, Stefanski A, Kapitaniak T. Synchronization of clocks. Phys Rep (2012) 517(1):1–69. doi:10.1016/j.physrep.2012.03.002
40. Dudkowski D, Czołczyński K, Kapitaniak T. Synchronization of two self-excited pendula: Influence of coupling structure’s parameters. Mech Syst Signal Process (2018) 112:1–9. doi:10.1016/j.ymssp.2018.04.025
41. Grebogi C, Ott E, Yorke JA. Critical exponent of chaotic transients in nonlinear dynamical systems. Phys Rev Lett (1986) 57:1284–7. doi:10.1103/physrevlett.57.1284
42. Frank KT, Petrie B, Fisher JAD, Leggett WC. Transient dynamics of an altered large marine ecosystem. Nature (2011) 477:86–9. doi:10.1038/nature10285
43. Tarnowski W, Neri I, Vivo P. Universal transient behavior in large dynamical systems on networks. Phys Rev Res (2020) 2:023333. doi:10.1103/physrevresearch.2.023333
44. Golomb D, Hansel D, Shraiman B, Sompolinsky H. Clustering in globally coupled phase oscillators. Phys Rev A (Coll Park) (1992) 45:3516–30. doi:10.1103/physreva.45.3516
45. Okuda K. Variety and generality of clustering in globally coupled oscillators. Physica D: Nonlinear Phenomena (1993) 63(3):424–36. doi:10.1016/0167-2789(93)90121-g
46. Hansel D, Mato G, Meunier C. Clustering and slow switching in globally coupled phase oscillators. Phys Rev E (1993) 48:3470–7. doi:10.1103/physreve.48.3470
47. Hou H, Zhang Q, Zheng M. Cluster synchronization in nonlinear complex networks under sliding mode control. Nonlinear Dyn (2016) 83:739–49. doi:10.1007/s11071-015-2363-z
48. Kuramoto Y, Battogtokh D. Coexistence of coherence and incoherence in nonlocally coupled phase oscillators. Nonlinear Phenom Complex Syst (2002) 380(5).
49. Abrams DM, Strogatz SH. Chimera states for coupled oscillators. Phys Rev Lett (2004) 93:174102. doi:10.1103/physrevlett.93.174102
50. Dudkowski D, Wojewoda J, Czołczyński K, Kapitaniak T. Transient chimera-like states for forced oscillators. Chaos (2020) 30(1):011102. doi:10.1063/1.5141929
Keywords: pendula systems, synchronization, multistability, bifurcation analysis, classical mechanics
Citation: Dudkowski D, Jaros P and Kapitaniak T (2022) Different coherent states for lightly supported coupled pendula. Front. Phys. 10:1021836. doi: 10.3389/fphy.2022.1021836
Received: 17 August 2022; Accepted: 12 September 2022;
Published: 26 September 2022.
Edited by:
Daniel Adrián Stariolo, Fluminense Federal University, BrazilReviewed by:
Manish Dev Shrimali, Central University of Rajasthan, IndiaDibakar Ghosh, Indian Statistical Institute, India
Edward Hellen, University of North Carolina at Greensboro, United States
Copyright © 2022 Dudkowski, Jaros and Kapitaniak. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dawid Dudkowski, dawid.dudkowski@p.lodz.pl