Theoretical Study of Synchronous Behavior in a Dual-Pendulum-Rotor System

A dual-pendulum-rotor system widely appears in aero-power plant, mining screening machines, parallel robots, and the like of the other rotation equipment. Unfortunately, the synchronous behavior related to the dual-pendulum-rotor system is less reported. Based on the special backgrounds, a simplifiedmechanical model of the dual-pendulum-rotor system is proposed in the paper, and the intrinsic mechanisms of synchronous phenomenon in the system are further revealed with employing the Poincaré method. The research results show that the spring stiffness, the installation angular of the motor, and rotation direction of the rotors have a large influence on the existence and stability of the synchronization state in the coupling system, and the mass ratios of the system are irrelevant to the synchronous state of the system. It should be noted that to ensure the implementation of the synchronization of the system, the values of the parameters of the system should be far away to the two “critical points”.


Introduction
"Synchronization phenomenon" widely exists in nature and daily life.The so-called synchronization phenomenon refers to self-adjusting different frequencies of oscillating objects to a unified frequency relying on their internal weak couplings [1,2].In recent years, the study of synchronization phenomenon is involved in the fields of physics, chemistry, and biology, such as the practical application on complex dynamic network systems [3][4][5], nonlinear coupled chaotic systems [6,7], coupled pendulum system [8][9][10][11][12], and rotor system [1,[13][14][15][16][17][18][19].Nevertheless, the latter two systems can be categorized as synchronous problem of mechanics system in detail.In the synchronous study of the pendulum system, antiphase selfsynchronization of clock pendulum is firstly discovered by Huygens in 1665.In the modern computation and analysis, the difference of spring stiffness influenced on the synchronous state is concerned, and the synchronization in inphase and anti-phase state for Huygens model is analyzed.Meanwhile, the derivatizations of Huygens model, including two coupled double pendula, coupling pendula under elastic forces, and pendula connected with linear springs, have attracted many scholars' attention.In the self-synchronous study of the rotor system, I.I.Blekhman proposed Poincaré method to explore the synchronous mechanism in multirotor systems.Nowadays, the Poincaré method is still widely used in modern engineering design.Based on the foundation of I.I.Blekhman's method, many scientists proposed the other approaches to study self-synchronization of the rotor system.Wen Bangchun introduced an averaging method to calculate the synchronization and stability between multiple coupled rotors [13].Zhao Chunyu [14,15] gave a revised small parameter method to describe the synchronous process of rotor system, which greatly simplifies the solving process of synchronization problem in rotor system.Sperling presented a two-plane automatic balancing device for equilibration of rigid-rotor unbalance, on which the synchronization of the rigid rotors is determined with numerical method.Similarly, Balthazar [16] examined self-synchronization of four nonideal exciters in nonlinear vibration system via numerical simulations.Djanan.A. A. N [17] explored the system, three motors working on the same plate; the synchronous state depends on the physical characteristics of the motors and the plate.The objects of the above-mentioned researches are mainly synchronization of the coupled pendula or the coupled rotors; however, the synchronization of the pendulumrotor system is less concerned.The extensive application of the pendulum-rotor system in aero-power plant, mining screening machines, parallel robots, and the like of the other rotation equipment is important for our industrial production and daily life [20][21][22].The aim of this paper is to investigate a system composed of two unbalanced rotors coupled with two pendula in a horizontal plane.But mass of the motor and the rotor, which can be relatively heavy or light, may affect the system motion.Therefore, the proposed problem is different from the dynamics of classical pendulum systems presented in papers [8][9][10][11] and classical rotor systems presented in papers [12][13][14][15][16][17][18][19][20][23][24][25][26][27].We consider an ideal source of energy supply to asynchronous motors, which drag the unbalanced rotors.We also investigate selected aspects of the synchronization phenomenon.There are different synchronization states possible in mechanical systems.We will search for a complete synchronization when the rotors perform rotation with the difference direction or for the synchronous state when the phase differences between the rotors are stabilized.

Dynamic Model
A simplified dual-pendulum-rotor system is shown in Therefore, the oscillating characteristics of the motors are just like a pendulum or pendula with length  and lumped mass (equal to the mass of the motors)   [Kg] (for  = 3, 4).The angular installation of the motors is denoted by  ∈  1 in [rad];  represents oscillation of centroid of the system in xdirection, respectively. is the rotation direction of the rotors;  [m] is distance between pivot of the motor and centroid  of the system.
In rotating coordinate       , the rotating centers of the rotors Φ  1 , Φ  2 can be expressed as And the rotating center of the motors Φ  3 and Φ  4 is represented as In fixed coordinate , coordinates Φ 1 , Φ 2 , Φ 3 , Φ 4 can be obtained by where Φ 0 is displacement vector of the vibrating base, Φ 0 = [, 0] T .According to kinetic theory, kinetic energy T of the dualpendulum-rotor system should be written as follows: 1 and  2 are the rotational inertia of the rotors, and symbols ( • ) and ( • ) denote (•)/ and  2 (•)/ 2 , respectively.In addition, potential energy  of the dual-pendulumrotor system can be expressed by Moreover, dissipated energy  of the system can be obtained by According to the Lagrange equation the dynamic equation of the dual-pendulum-rotor system can be derived.In the system, coordinate matrix q = [,  1 ,  2 ,  3 ,  4 ] T is considered as the generalized coordinate matrix; meanwhile, the generalized force matrix of the system is assumed as Substituting ( 4), ( 5), (6), and ( 8) into (7) and, meanwhile, considering  1 ,  2 ≪  0 , ∑ 4 =0   = .As  1 and  2 are sums of rotational inertia of motor's rotor about spin axis and that of eccentric lump, compared with  1  2 and  1  2 , the rotational inertia of motor's rotor about spin axis is too small to be neglected; i.e.,  1 ≈  1  2 ,  2 ≈  2  2 , and the dynamic equation of the dual-pendulum-rotor system is written by where 1, the two rotor rotate in the same direction, −1, the two rotor rotate in the opposite direction.(10) It can be seen that the terms of φ 2 3 and φ 2 4 , nonlinear terms owing to the small periodic vibration in synchronous state, are included in (9c).The terms related to φ 2 3 and φ 2 4 can be rewritten as Poincaré style with small parameter .Meanwhile, damping coefficient of this system is very small, and so the terms with damping coefficient can be rearranged as Poincaré style.In synchronous state, velocity of the motors is fluctuated slightly and periodically; hence, the acceleration of the rotors is also considered as a small term.Therefore, introducing small parameter , the Poincaré style of (9a), (9b),(9c), (9d), (9e) is given as follows [27]: where Based on Zhao's method [15], the electromagnetic torque of the asynchronous motor is linearized at stable point φ  =   when the two motors are synchronously operated; i.e., where   is coefficient of mutual asynchronous in the ith motor;   is the ith coefficient of stator inductance in the ith motor;   is pole-pairs of the motors;   is synchronous velocity of the motors;   is rotor resistance in the ith motor;  0 is amplitude of stator voltage.

Approximate Solution
From formulas (9a) to (9e), it can be seen that the vibration characteristics of the system are coupled related to DOFs  1 ,  2 , ,  3 , and  4 .Neglecting the terms of small parameters, introduce the dimensionless parameter as follows: In addition, velocities φ 1 and φ 2 of the two motors approached to a constant are represented by   in the steady state; therefore, the accelerations of the two motors should be equal to zero; i.e., φ 1 = 0 and φ 2 = 0.In this case, (11c) -(11e) can be simplified as the dimensionless style: For simplifying the solution procedure, consider that the natural frequencies of the torsion spring are quasi-identical; thus, parameter  ̃ should be replaced by   .Applying Laplace transformation in (16a), (16b), (16c), the coupling terms can be decoupled.Therefore, we have From formulas (16a)-(16c), the approximate solutions of the system in complex domain can be given by where .Then, the numerators and denominators of (17a), (17b), (17c) in the right are divided by divisor   4   2 ; moreover, introduce frequency ratios into (17a)-(17c).In this case, the coefficients of spring stiffness   and   are converted into the frequency ratios   and   in light of ( 14), respectively.Finally, through applying the inverse Laplace transformation to (18), the approximate solutions of the system in time domain can be obtained by where In light of functions  1 (), 2 (), and  3 () in ( 14), the approximate solutions of the system can be rearranged as the following:

Synchronization
Based on the approximate solutions above, the synchronization and stability of the system can be determined with Poincaré method [20,[25][26][27].In addition, the phase angle of the rotors can be assumed by Substituting ( 22) into (21a), (21b), (21c), the approximate solutions of the system can be written by

(23c)
As the values of the phase difference of the rotors may be different when the asynchronous motors are operated in the steady state, the phase difference between the rotors is noted by The second derivatives of ,  3 , and  4 with respect to time t can be obtained.Substituting ẍ , φ 3 , and φ 4 into (12a) and (12b), and then integrating and averaging these equations related to t over period , respectively, we can obtain 1 and  2 represent the averaging residual torque of the first and second motor in single period T, respectively.As we all know, when the motors are operated in synchronous state,  is a small parameter approached to zero.Therefore, the residual torque in the motors is approximated to zero; i.e., Substituting ( 27) into ( 25) and ( 26), the balanced torques of the motor are written by ,   ( = 1, 2) represent average electromagnetic and friction torque in the motors over period , respectively.When the two rotors are operated synchronously, the vibrating system transmits electromagnetic torque between the two motors to overcome the differences of output torque for the two rotors by adjusting the phase difference between two rotors.The differences of output torque Δ between the motors can be expressed as It can be known that denominator  is a function related to variables , 12 , 21 , 23 , and , but parameter  indicates the rotation directions of the motors independent of nonzero solution of denominator .Therefore, we will discuss zero solution of parameter  related to installation angle and coupling coefficients.There are two cases causing the zero solution of denominator .
(1) Parameter  represents installation angle of the motors.Thus, parameter  with a particular value leads to D being zero.The first "critical point" is  =  − , ( = 1, 2, 3 . ..) . (33) When the installation angles of the motors are approximated or equal to this point, the absent-synchronization of the system will be implemented.In this situation, the vibrating characteristics of the oscillating body are unascertainable.
(2) According to (20), coupling coefficients  12 ,  21 , and  23 are fractions; thus parameters   and   with particular values lead to  being zero.The second "critical point" is The second "critical point" means that the excitation frequency of the motors is identical to the natural frequency of the springs.In this case, strong resonance exists in such value of the frequency ratios.However, frequency ratios   and   are the function related to stiffness coefficients   and   of the springs, which indicates that stiffness coefficients are the key parameters to determine phase difference .
Clearly, to ensure the synchronization and stability of the system, the value of frequency ratios   and   should be far away from 1.In light of different value of   and   , the coupling type of the system can be defined as follows: Type 1: system of before-resonance coupled beforeresonance (0.1 <   < 1 and 0.1 <   < 1).Type 2: system of after-resonance coupled beforeresonance (0.1 <   < 1 and 1.0 <   < 7, or 1.0 <   < 7 and 0.1 <   < 1).Type 3: system of after-resonance coupled after-resonance (1.0 <   < 7 and 1.0 <   < 7).
Eliminating the two kinds of "critical point", the phase difference between the two rotors can be calculated:

Synchronous Stability
From (35), some values of phase difference can be obtained; however, the stability of these values should be determined as I.I.Blekhman's method [1,2].If a certain constant  =  * satisfies equation phase difference  * is stable under the condition that the value of  is negative.
Thus, the residual torque difference can be obtained by As   ,   ( = 1, 2) are irrelevant to phase difference , in light of (33), the criterion of synchronous stability is expressed Rearranging the equation above, the criterion of synchronous stability is simplified by  11 cos  +  12 sin ( + ) +  21 sin ( + ) The formula above shows that synchronous stability of the system is also determined by rotation direction of the rotors, coupling coefficients, and installation location of the motors.Only should parameters of the system satisfy the balance equation (30) and the synchronous stability criterion (39), the synchronous operation of the rotors can be implemented.In this case, the phase difference between the rotors is called the stable phase difference.

Numerical Analysis
The above-mentioned sections have given some theoretical discussions in the simplified form on synchronization problem for the vibration system that the unbalanced rotors are coupled with pendulum.In this section, we will employ some numerical analysis to discuss the stable phase difference, which can be calculated according to (30) and (39).The phase difference is determined by rotation direction of the rotors, coupling coefficients, and installation location of the motors.However, the coupling coefficients are the function of frequency ratios (    ), mass ratios (  ,   ), and motor installation angle ().Therefore, it can be included that the stable phase difference may be influenced by the parameters above.The parameter values corresponding to general engineering application are as given in Table 1.
6.1.The Rotors Operated in Opposite Direction.The stable phase difference is determined by considering the different values of parameters ,   ( = 1, 2),   , and   .Moreover, in light of ( 35) and (39), the rotation direction of the two rotors is opposite, and the stable phase difference between the two rotors is shown in Figure 2. According to ( 20) and (30), we have (− 11 + 23 ) sin  = 0 when  = −1, and the synchronous balance condition as in (17a), (17b), (17c) is simply expressed as sin  = 0. Similarly, the synchronous stability criterion can be rewritten by (− 11 +  23 ) sin  > 0. It then follows that the value of the phase difference is stabilized at 0 [rad] in the blue area, which describes that the synphase motion is stable and the antiphase motion is unstable ( [1] describes that the motion, as the existence of  ∈ (−/2, /2), is called synphase synchronization; and the motion, as the existence of  ∈ (/2, 3/2), is called antiphase synchronization).On the contrary, the brown area represents that the synphase synchronization is unstable and the antiphase synchronization is stable.Firstly, installation angle  and mass ratio   are fixed at /3 and 0.02, respectively, and the mass ratios ( 1 ,  2 ) are equal to 1.0.In such values of the parameters, the stable phase difference is shown in Figure 2(a).It is indicated that the stable phase difference is related to frequency ratios   and   ; in other words, the phase difference is dependent on the stiffness of the torque spring in the pendulum and the pressure spring in the vibrating body.Secondly, mass ratio   is fixed at 0.02, and the frequency ratios   and   are considered identical.Figure 2(b) describes that the stable phase difference is changed with the installation angle .It can be seen that the synchronous state of the system is locked in the synphase synchronization when   and   > 1.5, and the variation of installation angle  has little influence on the stable phase difference.However, the synchronous state is changed in region of   and   < 1.5.Finally, the effect of mass ratios   ,  1 , and  2 on the phase difference needs to be discussed.Thus, the parameter values in Figure 2 respectively, it can be found that the values of the stable phase difference are independent of the mass ratios.Therefore, the mass variations of rotors and vibration body would not affect the synchronous motion when the two rotors operate in the opposite direction.
6.2.The Rotors Operated in the Same Direction.In the following calculations, the two rotors rotating in the same direction are considered according to (35) and (39).Firstly, the parameter values in Figure 3(a) are the same as Figure 2(a) except for rotation direction .From Figure 3(a), it follows that two rotors may synchronize either synphase or antiphase, also depending on the values of   and   .The obtained results reveal that, in the interval of 0.1 <   < 1 and 0.1 <   < 1, synphase synchronization is implemented between the rotors; in the interval of 0.1 <   < 1 and 1 <   < 7, antiphase synchronization is implemented between the rotors; in the interval of 1 <   < 7, 0.1 <   < 1 and 1 <   < 7, 1 <   < 7, both antiphase and synphase synchronization collectively remained; in the interval of 1 <   < 7 and for 1 <   < 7, the two rotors will be synchronously operated in synphase state., respectively, it can be found that the stable phase difference is also independent of the mass ratios.In other words, synchronous state is less influenced by mass of the rotors when the two rotors operate in the same direction.results of this sample can be calculated through (9a), (9b), (9c), (9d), (9e).When the two motors are supplied by the electric source at the same time, the velocities of the two rotors are compatible (in Figure 4(a)).When the angular velocities of the motors reach the nominal velocity, then the motors oscillate steadily with the identical responses.At this moment, the coupling torques (in Figure 4(c)), making the phase difference  stabilized nearby 0 [rad], are fluctuated.In this case, the two motors rotate stably in synphase synchronization.From Figure 4(c), it follows the displacements of the vibration body and the two motors.It is clear that when the rotation velocity of the two rotors passes through the resonant region of the coupling system, the resonant responses of the system in the  3 -and  4 -directions appear in the starting process.In the synchronous state, the displacements of the vibrating body and the motors are stable, and the amplitudes of them are 0[m], 0.01[rad], and 0.01[rad], respectively.Comparing simulation results with Figures 2(a) and 2(b), the stable phase difference obtained by the computer simulation is consistent with the theoretical computation (i.e., the stable phase difference in Figures 3(a

Conclusions
In this paper, the dual-pendulum-rotor system is concerned.The research results show that the spring stiffness, the installation angular of the motor, and rotation direction of the rotors have a large influence on the existence and stability of the synchronization state in the coupling system, and the mass ratios of the system are irrelevant to the synchronous state of the system.It should be noted that to ensure the implementation of the synchronization of the system, the values of the parameters of the system should be far away to the two "critical points".To verify the correctness of the theoretical computations, some example simulations are preformed, and the results of theoretical computation are in accordance with that of example simulations.
The dual-pendulum-rotor system widely appears in aeropower plant, mining screening machines, parallel robots, and the like of the other rotation equipment.However, the dynamic characteristics and performance accuracy of the dual-pendulum-rotor system are mainly influenced by synchronous behavior between the rotors.In the early stage, for the developing and understanding the internal characteristics of the system, we only consider the vibrating body under the assumption of horizontal displacement.What is synchronization state of the system that the vibrating body simultaneously vibrates in the pitch and vertical directions?We believe that finding the answer to this question is the next step in challenging task of getting a complete understanding of synchronization in such system.

Fig- ure 1 .
This system consists of a rigid oscillating body of mass  0 [Kg] elastically supported via a linear damping spring with stiffness   [N/m] and damping   [Ns/m].Unbalanced rotor  actuated by an asynchronous motor is modelled by a point mass   [Kg] (for  = 1, 2) and attached at the end of a massless rod of length  [m]; the rotation angle of the rotors is defined by   ∈  1 (for  = 1, 2) in [rad].It should be noted that all the motors are installed rigid vibrating body by linear torsion springs with stiffness   [N/rad] and damping   [N ⋅ m/(rad/s)].Thus, small oscillating should exist in the two motors when the motors are operated in steady state; the oscillating angle is denoted by   ∈  1 (for  = 3, 4) in [rad].
(c) are identical to Figure 2(a) except for   .And the parameter values in Figure 2(d) are identical to Figure 2(b) except for  1 .Comparing Figures 2(c) and 2(d) with Figures 2(a) and 2(b),

Figure 3 (
b) shows the variation of the stable phase difference with the installation angle.It is indicated that the synchronous state of the system is obviously influenced by .Finally, the effect of mass ratios   ,  1 , and  2 on the phase difference is discussed.Thus, the parameters in Figure3(c) are identical with Figure 3(a) except for   , and the parameters in Figure 3(d) are identical with Figure 3(b) except for  1 .Comparing Figures 3(c) and 3(d) with Figures 3(a) and 3(b)

Figure 4 :
Figure 4: The dynamic characteristics when the rotors rotate in the opposite direction.
2    { 11 sin  −  12 cos ( + ) −  21 cos ( + ) 2  2  2    {− 11 sin  +  12 cos ( + ) +  21 cos ( + ) −  23 sin [ + ( + 1) ]} Obviously, the value of the phase difference is related to the coupling coefficients ( 11 ,  12 ,  21 , and  23 ) and the installation angle.Coupling coefficients represent the mutual coupling ability among the oscillating body, rotors, and pendula through the springs.The larger the coupling coefficients are, the stronger the coupling ability of the system is.Obviously, coupling coefficients  12 ,  21 , and  23 contain sin ; thus the absence of the coupling ability appears when  = 0.In this situation, the rotors cannot implement synchronous operation.It should be noted that coupling coefficients are the functions of parameters   ,   , , and   , and these parameters may influence the value of phase difference .