Synchronization and coupling dynamic characteristics of a dual-rotors exciter

In this work, some theoretical analyses, numerical simulations and experimental results on synchronization of a dual-rotors exciter are given. The exciter is made up of two rotors with eccentric masses (REMs) respectively driven by two DC motors with common axis. By adjusting the phase difference between two REMs to change the response amplitude, the decoupling between response amplitude and exciting frequency can be realized. The motion equations of the vibration system are established by using Lagrange equation, and the dimensionless coupling equations of that are obtained by applying the average method of small parameter. According to the existence condition of the zero solution of the dimensionless coupling equations, the synchronization condition of the vibration system is obtained. The stability condition of the vibration system implementing synchronization motion is acquired based on the principle of Hamilton. Through the comparison between numerical simulations and experimental results, the validity of theoretical analyses is proved, which helps the design of the dual-rotors exciter.


Introduction
Synchronization is a motion form which objectively exists in the natural world and the human society, as well as the fields of engineering and technology [1].
Huygens [2] is the first person who found the synchronization phenomenon in pendulum clock system.Blekhman [3][4][5][6][7] studied motion stability of double exciters by using the method of Poincare-Lyapunov small parameter, and first proposed the synchronization theory of exciters.
Wen [1,[8][9][10][11] considered the effect of system damping on the basis of the former, and converted the motion equations of REMs into the equations of phase differences of REMs after average integral.Balthazar [12,13] gave some comments on self-synchronization mechanisms of two non-identical exciters and four non-ideal exciters by numerical simulations.Quinn [14] studied the dynamics and synchronization under resonance case.Nijmeijer [15] studied dynamic control on synchronization from a control theory perspective.Yamapi [16,17] studied the dynamics and synchronization of two coupled self-excited devices.Perlikowski [18] described the relation between the complete, phase and generalized synchronization of the mechanical oscillators driven by the chaotic signal generated by the driven system.Acebron [19,20] applied the Kuramoto model to explain synchronization phenomenon in large population of phase oscillators.Tanaka [21] analyzed a large system of nonlinear phase oscillators with sinusoidal nonlinearity, uniformly distributed natural frequencies and global all-to-all coupling, which is an extension of Kuramotor's model to second-order system.
Further research on vibration synchronization theory and broaden its application domain, in order to satisfy the demands of practical engineering, is one of the current research hotspots in the field of synchronization.
When the phase difference between two REMs in the vibration system is stable in a certain value, it will make the system work in a particular motion form [8][9][10][11].We can apply this principle to design vibration machines, which satisfy various demands such as motion trajectory and amplitude.The rotation axes of two REMs in many articles are parallel, there is not the case that the rotation axes of two REMs are common.A novel design for small vibrotactors called the dualrotors exciter is first presented by Miklos [22,23], which makes it possible to produce vibrations with independently adjustable frequency and amplitude.However, for vibration equipment, synchronization is the result that vibration system automatically selects energy distribution by itself.In engineering, whatever kind of synchronization is adopted, only understand the coupling dynamic characteristic of the vibration system can reasonable limit the phase difference between two REMs.
This paper discusses the coupling dynamic characteristic of the dual-rotors exciter based on synchronization theory, and it contains the following elements: First, the motion equations of the dual-rotors exciter are established.Second, the synchronization condition and the stability condition of synchronization state are deduced.Third, numerical discussions are provided.Fourth, experiments are given.Finally, concludes this paper.

Dynamic model of the dual-rotors exciter
Fig. 1 illustrates the dynamic model of the dual-rotors exciter, in which springs are connected to a rigid frame.The two DC motors are installed in the rigid frame with the common axis, which drive two REMs rotating in the same direction to excite the vibration system.The frame is a fixed frame, and its origin is the equilibrium point of centroid of the rigid frame, meanwhile the frame ′ ′ is a non-rotation moving frame.The motions of the rigid frame are vibrations inand -directions, denoted by and .Each REM rotates about its rotation axis, denoted by and . is the mass of the rigid frame.( where = + + is the vibration mass; ≈ is the moment of inertia; is the damping coefficients; and are the electromagnetic torques of two motors; (⦁ ) and (⦁ ) denote ⦁/ and ⦁/ , respectively; Above mentioned, = 1, 2.

Synchronization condition and stability condition of synchronization state
We assume the average phase between two REMs and their phase difference are and 2 , respectively [8][9][10][11], then we obtain, = + , = − .Since the periodical vibration of the vibration system, the changes of mechanical angular velocities of two motors are periodic as well.The average value of their average velocities over the least common multiple period between two REMs must be constant, = constant.Assuming the instantaneous change coefficients of and to be and , we obtain = (1 + + ) , = (1 + − ) .According to the previous works [8][9][10][11], and selecting = = , = = , the responses of Eq. ( 1) can be approximately expressed in the forms: where: and − is the phase angle, ( = , ).
Differentiating and in Eq. ( 2) with respect to time , we obtain , and .Substituting them into the last two equations in Eq. ( 1), and integrating them over = 0~2 , and neglecting the high-order terms of and , the motion equations of two REMs are expressed in the forms: where: Compared with the change of ( = ) with respect to time , those of , , , and are very small.Therefore, they are considered to be slow-changing parameters, while the change of is considered as fast-changing parameter [8][9][10][11].During the aforementioned integration over = 0~2 , , , , and can be assumed to be the middle valves of their integrations , ̅ , ̅ , ̅ and ̅ , respectively.
If the variation coefficient of mechanical angular velocity of the REM is during the steady state operation of the vibration system, i.e., = (1 + ) .Applying the first-order Taylor expression around , the electromagnetic torque of the DC motor in the steady state can be obtained as follow [24] where: and is voltage, is armature resistance, Φ is exciting flux, is electromotive force constant, is torque constant, is rotational speed of the motor.The electromagnetic torques of two motors in Eq. ( 4) can be expressed as follows: We assume the parameters of two motors are the same, = = , = = .Substituting Eq. ( 5) into Eq.( 3), and writing Eq. ( 3) in matrix form that adding two formulas in equations to obtain the first row, and subtracting second formula from the first one to get the second row, then we obtain: where: In Eq. ( 6), ̅ and ̅ represent the perturbation parameters that the average angular velocities of two motors and the phase difference between two REMs.Eq. ( 6) describes the coupling relation of two REMs and is referred to as the dimensionless coupling equation of those.

Synchronization condition
If the vibration system can implement synchronization motion of two REMs, the average values and during the interval of the single period = 2 / must be zero, i.e., ̅ = 0 and ̅ = 0.Then, we can obtain = 0 and = 0 from Eq. ( 6) as follows: where the term + is the sum of the electromagnetic torques of two motors; the term 2 is the sum of the damping torques of two REMs.= 2 (1 + cos2 ) is the sum of the load torques that the vibration system acts on two motors.So Eq. ( 7) is the torque balance equation that the vibration system operates in the steady state.
According to Eq. ( 8), the electromagnetic torque difference between two motors is defined as = − , the torque of frequency capture is defined as = 2 .Rearranging Eq. ( 8) as follow: Since |sin2 | ≤ 1, the synchronization condition is that the torque of frequency capture must be greater than or equal to the absolute value of the electromagnetic torque difference between two motors, ≥ | |.From Eq. ( 3), when the vibration system operates in steady synchronization state, the load torques that the vibration system acts on the two motors can be expressed as follows: The adjusting torque is defined as sin2 .Assuming the phase of motor 1 leads that of motor 2, the adjusting torque acts on the motor of leading phase as load torque to decrease its angular velocity, and acts that of lagging phase as driving torque to increase its angular velocity, which ensure two REMs operate in synchronization state.If the voltages of two motors are different, their electromagnetic torques are also different.The load torques of two motor are different because the adjusting torque, which ensure there is a phase difference between two REMs.Applying this principle, this phase difference can change the response amplitude with satisfying the synchronization condition, which will be explained in the following chapter.
Under the condition of satisfying the synchronization condition, we define absolute ratio between the torque of frequency capture and the sum of the load torques as synchronization ability coefficient , which represents the ability that the vibration system adjusts the load torque of each motor to implement synchronization motion of two REMs.The greater value of synchronization ability coefficient, the easier implementing synchronization motion: where is the nonlinear function of .We will explain its nonlinear characteristic in the next section.

Stability condition of synchronization state
The kinetic energy and potential energy of the vibration system: In one vibration periodic, Hamilton's average action can be expressed as: According to Eq. ( 9), the solution of 2 has two values, which correspond two different states, and using 2 * to represent the steady solution.Based on the principle of Hamilton, when the second-order derivative of Hamilton's average action with respect to α is greater than zero, the 2125.SYNCHRONIZATION AND COUPLING DYNAMIC CHARACTERISTICS OF A DUAL-ROTORS EXCITER.
From the Eq. ( 15), we can know that only the product of stability ability coefficient and cosine function of the phase difference is greater than zero to ensure the stability of synchronization state.Namely when > 0, the stable domain of 2 * is (− 2 ⁄ , 2 ⁄ ), and when < 0, the stable domain of 2 * is ( 2 ⁄ , 3 2 ⁄ ).

Numeric discussions
The previous chapter has given some theoretical analyses in the simplified form on synchronization issue.This section will analyze quantitatively the numerical results of dynamical characteristics of the vibration system to verify the validity of the previous theoretical results.The parameters of the vibration system in this paper are: = 744 g, = 23 g, = 10 mm, = = 4930 N/m, = = 30N•s/m, = 0.0001, Φ = 0.0043 Nm/A, Φ = 0.0409 Vs/rad, = 5.5 Ω.According to Eq. ( 2), when two REMs operate in synchronization state, we can deduce the max amplitude of response in simply form with using the method of mathematical extreme solution as follow: where: According to Eq. ( 16), we can obtain the relation between the phase difference and the max amplitude of response as shown in Fig. 2. The dimensionless parameter is decreasing with increasing the phase difference.This is show that the response amplitude of the vibration system changes with changing the phase difference, which can achieve the decoupling between the response amplitude and exciting frequency.XIAOZHE CHEN, XIANGXI KONG, YUNSHAN LIU, BANGCHUN WEN should be noteworthy that, the vertical line = 1, which splits the domain in Fig. 3 into two parts.
One is sub-resonant domain ( < ) on the left hand side, and another is super-resonant domain ( > ) on the other side.As shown in Fig. 3, we can see that the closer resonance point, the smaller synchronization ability coefficient.
Fig. 4 shows the stable ability coefficient of synchronization state with different mass ratio.Its tendency is similar to Fig. 3 that there is a break point at = 1.When the vibration system operates in sub-resonant state, > 0 and the stable domain of 2 * is (− 2 ⁄ , 2 ⁄ ), and when it operates in super-resonant state, < 0 and the stable domain of 2 * is ( 2 ⁄ , 3 2 ⁄ ).

Numerical simulation
In this section, there are two simulation results are carried out, one case is two REMs operate in sub-resonant domain and another is two REMs operate in super-resonant domain.For the convenience of observing the initial state, we put the record time of data delay two seconds in post-processing, and below no longer explain it.
Fig. 5 shows the case that two REMs operate in a sub-resonant domain.During the starting process in Fig. 5(a), the angular velocities of two motors are equal to each other because their physical parameters are same.When time at about 10 s, the phase of motor 2 adds a disturbance of /6 in Fig. 5(b), after that, the state of the vibration system gradually returns to the previous one, which indicate that the stability of synchronization motion is very strong.In this stage, the synchronization angular velocity nears to 758 r/min, the phase difference is equal to zero.Later, the voltage of motor 2 is subtracted 0.1 V, 0.1 V, 0.15 V, 0.05 V, after those two REMs are de-synchronization.Before de-synchronization of two REMs, the phase difference reaches 70.3°, which satisfies synchronization condition.From the above process of changing voltage of motor 1, it is seen that the phase difference increases gradually, and the amplitude of response in vertical direction decreases gradually in Fig. 5(c), which accord with Fig. 2. In Fig. 6, two REMs operate in a super-resonant domain.The same as Fig. 5, at about 10 s, a disturbance of /6 is also added to the phase of motor 1, after about 5 s, the state of the vibration system gradually return to the previous one.In this stage, the synchronization angular velocity XIAOZHE CHEN, XIANGXI KONG, YUNSHAN LIU, BANGCHUN WEN nears to 1654 r/min in Fig. 6(a), the phase difference is equal to 180° in Fig. 6(b).Later, the voltage of motor 1 is added 0.2 V every 10 s until de-synchronization of two REMs.Comparing the results in Fig. 5, two REMs stably operate in the super-resonant with opposite phase.As shown in Fig. 6(c), the response amplitude of -direction is very small because the phase difference leads to decrease the sum exciting forces of two REMs.From these analyses, these dynamic characteristic accord with the theoretical results in previous sections.We can know that the work point of the vibration system has better select at sub-near-resonant to acquire big enough amplitude.The experiment schemes according to simulation ones are achieved in Fig. 5 and Fig. 6.Fig. 8 shows the case that two REMs operate in a sub-resonant state.When time at about 10 s, the power of motor 2 is cut during 1 s, the angular velocity of motor 2 quickly descend, the phase difference exhibits non-periodicity behavior, and the response amplitude of -direction also changes.

Experiments
In Fig. 9, two REMs operate in a super-resonant domain.Since the physical parameters of two motors and the structure symmetry of the vibration system, etc., the response amplitude is not equal to zero when the phase difference closes to 180°.
The results of matching between numerical simulations and experimental results show they are very good consistency.Through the comparison between numerical simulations and experimental results, the validity of theoretical analyses is proved, which helps the design of the dual-rotors exciter.

Conclusions
From the theoretical analyses, numerical simulations and experimental results in the present paper, the following remarks should be stressed: This paper proposes an analytical approach to investigate synchronization and coupling dynamic characteristic of a dual-rotors exciter.By introducing the average method of small parameters, the dimensionless coupling equation of the vibration system is deduced, which describes the coupling dynamic characteristics of the REMs.
The synchronization condition of two REMs operated in synchronous state is derived, which is that the torque of frequency capture is equal to greater than the absolute value of the difference between the residual torques of two motors.With the adjustment of the phase difference between two REMs to change response amplitude, the decoupling between response amplitude and exciting frequency can be realized under the synchronization condition.
According to the principle of Hamilton, the stability condition of the vibration system implementing synchronization is acquired.The coupling dynamic characteristic that the vibration system has selecting motion is discussed.For a dual-rotors exciter system with small damping, the phase difference between two REMs is stabilized at (− /2, /2) with a sub-resonant state, while it is stabilized at ( 2 ⁄ , 3 /2 ) with a super-resonant state.By the comparison of numerical simulations and experimental results, the feasibility of theory method is proved.The synchronization condition and stability condition of synchronous state can help to the design of a dual-rotors exciter.

Fig. 1 .
Fig.1illustrates the dynamic model of the dual-rotors exciter, in which springs are connected to a rigid frame.The two DC motors are installed in the rigid frame with the common axis, which drive two REMs rotating in the same direction to excite the vibration system.The frame is a fixed frame, and its origin is the equilibrium point of centroid of the rigid frame, meanwhile the frame ′ ′ is a non-rotation moving frame.The motions of the rigid frame are vibrations inand -directions, denoted by and .Each REM rotates about its rotation axis, denoted by and . is the mass of the rigid frame.and are the masses of two REMs.and are the eccentric radiuses of two REMs.and are the angular velocities of two REMs.and are the constants of springs, and are the damping constants in -and -directions, respectively.

Fig. 5 .Fig. 6 .
Fig. 5. Dynamic characteristic of the vibration system with the sub-resonant state: a) angular velocities, b) the phase difference, c) the response amplitude in -direction AND COUPLING DYNAMIC CHARACTERISTICS OF A DUAL-ROTORS EXCITER.

Fig. 7
Fig. 7 shows the mechanical composition of the dual-rotors exciter, which consists of two DC vibration motors, the vibration rigid, the springs and the support base.No. 1 represents REMs, No. 2 represents DC motors, No. 3 represents the suspension system, No. 4 represents the acceleration sensors, and No. 5 represents photoelectric rotation velocity sensors.The signals of the vibrations in the horizontal and vertical directions, and the rotation velocities of two REMs are collected by 3650D measurement frontends of B&K.While the powers of the two motors are supplied by GPD-3303S of Good Will Instrument.

Fig. 8 .Fig. 9 .
Fig. 8. Experiment results of the vibration system with the sub-resonant state: a) angular velocities, b) the phase difference, c) the response amplitude in -direction

Yunshan
Liu is currently a Ph.D. candidate in the School of Mechanical Engineering and Automation, Northeastern University, and a lecturer of Gudaojiaotong Polytechnic Institute.His research interests are mechanical dynamics, nonlinear vibration, vibration utilization engineering.Bangchun Wen is currently an Academician of Chinese Academy of Science and a professor at Northeastern University, China.He has systematically studied and developed the new course of "Vibration Utilization Engineering", his research interests include nonlinear vibration in engineering, vibration utilization engineering, and vibration synchronization theory.