Shock and Vibration
Volume 2016 (2016), Article ID 6973597, 14 pages
http://dx.doi.org/10.1155/2016/6973597

Research Article

On the Synchronization of Two Eccentric Rotors with Common Rotational Axis: Theory and Experiment

Xiaozhe Chen, Xiangxi Kong, Xueliang Zhang, Lingxuan Li, and Bangchun Wen

School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110819, China

Received 28 March 2016; Revised 22 May 2016; Accepted 22 June 2016

Academic Editor: Hassan Haddadpour

Copyright © 2016 Xiaozhe Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We study synchronization of two eccentric rotors (ERs) with common rotational axis in the vibration system of the far-resonant spatial motion. We deduce the dimensionless coupling equation of two ERs with applying the average method of small parameters. We convert the synchronization problem into the existence and stability of solving the zero solutions for the dimensionless coupling equations. By introducing the synchronization torque and the difference between the residual torques of two motors, we obtain the synchronization condition that two ERs achieve the synchronized motion. We derive the stability condition of the synchronized motion, which satisfies Routh-Hurwitz criterion. We discuss numerically the choosing motion feature of the vibration system, which indicates that the vibration system has two steady motion modes. The synchronization torque forces the phase difference to approach when the structural parameters of the vibration system satisfy the condition of the spatial cone motion, and the synchronization torque drives the phase difference to approach zero when those satisfy the condition of the spatial circle motion. Finally, through the comparison and analysis of experimental data, the structural parameters of the vibration system satisfying the above two conditions can guarantee the synchronization stability for two ERs.

1. Introduction

In vibration utilization engineering, synchronization of ERs (also called exciter) is a special phenomenon of the nonlinear vibration system, which is used to design the vibration conveyer, the vibration screen, the vibration centrifuge, and so forth [1]. The vibration mill is a kind of grinding equipment with high efficiency and energy saving, which has been widely used in many fields, such as metallurgy, chemical, mining, medical, ceramics, cement. In order to satisfy the increasing of powder processing, the vibration mill is large-scale development. However, for the large-scale vibration mill, its longer barrel is not driven by only one motor. There are many questions about the design of this vibration mill excited by multiple ERs, for example, two motors drive separately two ERs running in synchronization state to replace forced synchronization. Hence, applying the synchronization theory to design the large-scale vibration mill can simplify structure of the vibration mill, decrease production cost of that, and enhance work efficiency of that.

For a long time, utilizing the advantage of synchronization in engineering and technology fields has been the main aims of scholars, who have paid attention to synchronization of two or more ERs. When it comes to synchronization, we may date back to that of pendulum clocks found firstly by Huygens [2]. In the 1960s, Blekhman first clarified synchronization principle of two ERs in a vibration system and has successfully solved some problem about synchronization [3–8]. Since then, there are many excellent results which pushed the researches of the synchronization theory. Wen et al. selected the phase difference of ERs as object and applied the method of integral average to simplify the solution of synchronization [1]. Based on this concept, Zhao et al. applied the average method of small parameters to obtain the synchronization condition and the synchronization stability for two ERs [9–11]. This method considers the effect of damping and dynamic characteristics of induced motor. By applying the average method of small parameters, Zhang deduced the condition that vibration synchronization transmission for two ERs with reverse rotational direction and gave the result of the select motion characteristic for the vibration system driven by two ERs with same rotational direction [12, 13]. He also gave the result that using three ERs instead of two on a single base cannot improve the effective power in the far-resonant vibration system [14]. There are also scholars who applied numerical analysis to study the synchronization problems. Balthazar analyzed the nonlinear phenomenon of synchronization of two nonideal sources by means of numerical simulations [15, 16]. Bonkobara et al. applied shooting method to solve the displacement response curves of the vibration system to judge whether the vibration system achieves synchronization [17]. Rumyantsev and Tarasov addressed some comments on the nonlinear dynamics of vibration transport machines by the numerical simulation [18].

The above proposed result all selected two ERs or multiple ERs in a vibration system of plane motion as the object. However, there are few researches for spatial motion. Zhao et al. investigated the coupling dynamic characteristic for a vertical conveyer [19]. Miklós presented a dual-rotors exciter to produce vibrations by independently adjusting frequency and amplitude [20–22].

Hence, the research of synchronization of two ERs rotating with the common rotational axis is of great significance.

In this paper, synchronization of two ERs with the common rotational axis in the vibrating system of the far-resonant spatial motion is investigated with applying the average method of small parameters. We addressed that this type of vibration system has the choosing motion feature, and the theory approach is verified to be useful and feasible by experiments. This paper is summarized as follows: first, the motion equations of the vibrating system are derived. Second, the synchronization condition and the synchronization stability are deduced. Third, the quantitative numeric discussions are provided. Fourth, experiments are shown. Finally, conclusions are given.

2. Dynamic Model and Motion Equation of the Vibration System

As shown in Figure 1(a), the dynamic model of the vibration system consists of a rigid body and two ERs. Four symmetrical isolation springs support the rigid body. Two induction motors drive two ERs with the same rotational direction, respectively, and their rotational axes of two ERs are common. and are the rotational centers of ER 1 and ER 2, respectively.

(a)

(b)

Figure 1: Dynamical model of the vibration system: (a) schematic diagram of structure; (b) the coordinate frames: , the earth fixed inertial frame, , the body fixed frame.

In Figure 1(b), let denote the earth fixed inertial frame and let denote the body fixed frame rigidly attached to the vibration system with origin in the center of mass.

Generally, because of the elastic foundation, the rigid body has six degrees of freedom. In engineering, the motion of this type of machine is designed for the spatial circle motion. According to the design, rubber spring is chosen. The spring of ROSTA only has the plane motion because of the bigger rigidity in the axis direction. In order to simplify analysis, the -direction motion is considered to be fixed. Due to the symmetrical structure of the rigid body, the -axis is assumed to coincide with the principle axis of the rigid body. Since points , , and are on -axis, so the resultant force of the exciting forces passes the center of mass, and thus the rotation around -axis is not existing.

To sum up, the vibration system has six degrees of freedom as follows: the rotational directions of ER 1 and ER 2 are denoted, respectively, by and . The center of mass translations is and , and two angular rotations are and .

The coordinates of ER 1 and ER 2 are and in , respectively, and their expressions are described as

By applying the method of the rotation matrix , the coordinates of ER 1 and ER 2 are gained in (see (A.1)). Assuming that coordinates of the center of mass are in , the coordinates of ER 1 and ER 2 in can be expressed as follows:

The kinetic energy of system is

where is the mass of the rigid body; are the masses of ER , respectively; are moments of inertia of the rigid body with rotating about -axis and -axis, respectively; are moment of inertia of two motors’ rotors, respectively; denotes , .

In fact, the translational and rotational displacement of the rigid body are very small in the vibration condition. So the rotation matrix (see (A.2)) can be simplified into , and the kinetic energy of the vibration system can also be simplified.

As illustrated in Figure 1(a), four springs have the same physical characteristics. If the vibration system does not run, the coordinates of the connection points between the springs and the rigid body overlap for and :

So the coordinates of the connected points () in can be expressed in the form

Then the potential energy of the vibration system is

where are stiffness matrix of spring .

The viscous dissipation function of the vibration system is

where are damping matrix of spring .

By using Lagrange equation, the motion equations are

where are the generalized coordinates and are the generalized forces of the vibration system.

The generalized coordinates are , and the generalized forces are , , and , in which and are the electromagnetic torques of motor 1 and motor 2, respectively.

Use (3), (6), and (7) and apply Lagrange equation (8) to obtain the equations of motion of the vibration system as follows [9–14]:

where and , and , and , , and are moment of inertia of ERs, and are damping coefficient of the axis of motor, and denotes .

3. Synchronization of Two ERs

Assuming that the average phase of two ERs is , the phase difference between two ERs is :

Since the response of the vibration system is periodic, the changes of angular velocities of two motors are periodic as well. Assuming the lowest common multiple of period between two ERs is , the average value of over must be constant; that is,

We assume that the synchronization velocity between two motors is , and the instantaneous change coefficients of and are and , respectively; then we obtain

Substituting (10) into the first four formulas of (9), we obtain the responses of the steady state of the vibration system as follows [9–11]:

where

, , , and are the natural frequency of the vibration system in -, -, -, and -directions; , , , and are the corresponding damping ratios of the springs in -,-, -, and -directions; , , , and are phase angle in -, -, -, and -directions, respectively; is mass ratio between the mass of ER 1 and the mass of the vibration system; is mass ratio between the mass of ER 1 and that of ER 2; () is defined as the ratio between and (); and are the equivalent rotating radius of the vibration system about the -axis and -axis of the rigid body.

3.1. The Dimensionless Coupling Equation

Inserting , , , and (with respect to time ) into the last two formulas of (9) and integrating them over , respectively, we have

with

The coefficients of , , , and in (16) are listed in Appendix B. Based on the direct motions separation [9–11], , , , , and are considered to be their integral midvalves , , , , and , respectively.

According to the electromagnetic torques formula of induction motor, we have [23]

where

is the number of pole pairs; is the stator inductance, is the rotor inductance, and is the mutual inductance; is the stator resistance and is the rotor resistance; is the synchronization electric angular velocity and is the rotor electric angular velocity; is the phase voltage.

Substituting (16) and (17) into (15), we obtain the dimensionless coupling equations of two ERs as follows:

where , and are matrices, , , ,

3.2. The Synchronization Condition

Substituting and into (19), we obtain and . According to , we have

where is called the synchronization torque, ; is called the difference of the residual electromagnetic torques between the two motors, ; and .

Since , the synchronization condition that the vibration system carries out the synchronized motion of two ERs is

Hence, the synchronization condition of two ERs achieving the synchronized motion is that the synchronization torque is greater than or equal to the absolute value of the difference between the residual electromagnetic torques of two motors .

We define ratio value between the synchronization torque and the load torque as the coefficient of synchronization ability, , as follows:

where is the load torque that the vibration system acts on two motors, . The greater the coefficient of synchronization ability, the stronger the synchronization ability of the vibration system and the easier to carry out synchronized motion of two motors.

3.3. Stability of the Synchronized Motion

If the vibration system satisfies the synchronization condition, we solve and for and , which are denoted by and , respectively. Linearize (19) around , consider , use the notation , and write the dimensionless coupling equation as follows:

where

Assuming in the form of exponential time-dependent, we deduce the characteristic equation for the eigenvalue by solving the determinant equation , as follows:

where , , , and

In engineering, because the maximum of damping constant of the vibration systems is very small [9–11], can be neglected with comparing in , , and . Then, we simplify as follows:

Applying Routh-Hurwitz criterion [9–11], we deduce the stability condition of the synchronized motion of two ERs when the trivial solution is stable, as follows:

Based on the sign of , we rewrite inequality (29) as inequalities (30) and (31) as follows:

Since and , we can obtain from and ,

According to , we obtain

Substituting the expressions of into ,

When , inequality (34) is true. Hence, inequalities (32) and (33) satisfy inequality (34).

If , needs and needs . However, the left side of inequality (34) is false. Hence, , , and cannot satisfy the condition of .

In a word, inequality (30) satisfies Routh-Hurwitz criterion. So we obtain the stability condition of two ERs achieving the synchronized motion, which is . If , satisfies the stability condition, and if , satisfies the stability condition. So the vibration system has two kinds of stable synchronized motion state when the synchronization condition is satisfied.

4. Numerical Results and Discussions

Some simplified forms of theoretical analysis for the synchronization problem have been given in previous section. This section will present numerical results to compare with theoretical analysis for the synchronization condition and the stability condition.

The dimensionless parameters of the vibration system, , , , and , have an effect on the synchronization condition and the stability condition and are also functions of , , , , , , , and . Since the motion trajectory of this type of vibration machine is assumed as circle in -plane, such as the vibration mill that achieves the spatial circular motion, we assume , , and . However, , , , and change little (24/25~99/100) that can be treated as constant [12]. Herein, we will study the function of dimensionless parameters , , and on the synchronization ability and the stability of synchronization.

The actual parameter values of the vibration system are  kg,  kg,  kg·m²,  kN/m,  kN/rad,  kNs/m, and  kNms/rad. The parameters of the two motors are the same:  Ω,  Ω,  mH,  mH, and  mH.

We give the expression of versus and to reflect effects on the synchronization ability. Based on the expression of , its maximum can be simplified in the form

The synchronized motion of the vibration system is stable when it satisfies (33) and (34). Figure 2 shows that for ; we will arrange the value of from 0 to 7 in next discussions.

Figure 2: Maximum of .

According to the theoretical analysis in previous section, , to satisfy the synchronization condition. When two identical motors drive the two different masses of ERs, we can obtain the simplified form of (22):

Figure 3 shows the regions of implementing the synchronized motion in -plane. We divide the -plane into three regions. The phase difference that the vibration system carries out during synchronized motion is in Region I, and it is in Region II. The vibration system cannot carry out the synchronized motion in Regions I and II. Regions I, II, and III converge into a point ( and ), and at this point. We know that this point is the optimal parameters of two ERs to enhance the synchronization ability.

Figure 3: Regions of implementing synchronization.

4.1. Implementing the Synchronized Motion

Since the dimensionless parameter is the function of , its value changes with changing and has zero point. The synchronization torque is the function of , so it has similar changes with . As mentioned in theoretical analysis, the greater the values of the synchronization torque, the easier the two ERs to implement synchronization. Figure 4 shows the changes of with different mass ratio between ERs. has zero point in every case, and the point at about depends on structure parameters of the vibration system. According to the regions divided in Figure 3, the vibration system cannot implement the synchronized motion when . obtains the maximal values at , and the vibration system operates at its best synchronization state. In addition, the greater , the greater . Comparing with , the vibration system can easier carry out the synchronized motion for . When , increases with increasing . By comparing the synchronization torque with the different mass ratio, the better the symmetry of the structural parameter, the greater the synchronization torque.

(a)

(b)

Figure 4: Synchronization torque: (a) ; (b) .

Figure 5 shows the changes of with different mass ratio between ERs. Since is the function of , , and and they are the function of , has no effect on . In addition, has zero point for every . From Figure 5, it can be seen that all first decrease with increasing and then decrease. To sum up, plays a decisive role in the adjustment of the synchronization ability.

Figure 5: Coefficients of synchronization ability.

4.2. Stability of the Synchronized Motion

The stability condition of two ERs achieving the synchronized motion is . When the stability condition is satisfied, determines the change of the phase difference . In addition, the phase difference determines the motion type of the vibration system. In other words, if the structural parameters of the vibration system satisfy the synchronization condition, the motion type is determined by . This is the choosing motion characteristic that the motion type of the vibration system depends on its structural parameters.

Figure 6 shows the regions where the phase difference changes with the different mass ratio between ERs. From (21), plane is divided into two regions of and . changes at with every . If the structure parameters of the vibration system are completely symmetrical, such as the mass-radius product of two ERs, achieves the change of to with increasing . By comparing with , the better the symmetry of the structural parameter of the vibration system is, the easier approach and .

Figure 6: Approximate values of versus .

For further study of the effect of the structural parameters of the vibration system on the synchronization stability, we apply bisection method to obtain the stability domains of the synchronized motion. Figure 7 shows the changes of the coefficients of the synchronization stability , , and with the different mass ratio between ERs. From the stability condition, we know that the vibration system runs in the steady state when , , and . Since is the function of , it has zero value at . By comparing the changes of , , and , it can be seen that is always greater than zero. In the region of the phase difference tending to , increases with increasing , and and decrease with increasing . It is noteworthy that goes earlier through zero than . According to the stability condition, the vibration system is unstable at the moment. Comparing comprehensively , , and , we know that the vibration system can operate in steady state when . Therefore, the stability condition of two ERs achieving the synchronized motion is . In addition, the smaller , the bigger the synchronization stability of the vibration system.

(a)

(b)

Figure 7: Coefficients of stability: (a) ; (b) .

5. Experiments

In this section, we will verify the validity of the dynamic model, theoretical analysis, and numerical results of the above sections, by comparing to experimental results for a laboratory model. Figure 8 shows the mechanical composition of the laboratory model, which consists of two AC vibration motors, the rigid vibration base, the springs, and the support base. Number 1 represents AC motor, number 2 represents ER, number 3 represents the rigid vibration base, number 4 represents the rubber spring, number 5 and number 6 represent the positions measured by the acceleration sensors, and number 7 represents the angular velocity measured by Hall sensors. The angular velocities of two motors and the signals of the vibrations in the vertical and horizontal directions are collected by NI data acquisition card, while two motors operate in 45 Hz power frequency by Siemens inverters.

Figure 8: Experiment equipment ().

As shown in Figure 9, , motor 1 first runs in the starting process, and motor 2 begins to run after 3 s. Because of the identical motors and the same inertial moments of two ERs (), the synchronization torque begins to adjust the phase difference between two ERs to achieve the synchronized motion, with approaching the rotational velocities of two motors. Since the adjustment of the synchronization torque, the synchronization velocity and the phase difference are stabilizing fast, as shown in Figures 9(a) and 9(b), the synchronization velocity is about  r/min, and the phase difference is about . Due to the fact that two ERs operate in far-resonant, they excite the resonant responses through the resonant region of the vibration system in - and -directions, as shown in Figures 9(c)–9(f). Figures 9(c) and 9(d) show the responses of the position for number 5, and Figures 9(e) and 9(f) show the responses of the position for number 6. From comparison between the response of number 5 and that of number 6, it can be seen that the main motion type of the vibration system is the spatial cone motion.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 9: Experiment results with , , and : (a) rotational velocities; (b) the phase difference; (c) response in -direction for number 5; (d) response in -direction for number 5; (e) response in -direction for number 6; (f) response in -direction for number 6; (g) current of motor 1.

At about 30 s, the power of motor 2 is switched off, which breaks the original synchronization state between two ERs. Under the function of the synchronization torque, the vibration system can achieve another synchronization state by adjusting the phase difference. This phenomenon is called the vibration synchronization transmission [1]. Figure 9 shows that, after switching off the power of motor 2, the synchronization velocity is about  r/min, and the phase difference is about . The responses of number 5 are slightly larger, and those of number 6 are slightly smaller.

Figure 9(g) shows the root mean square of stator current of motor 1. Since the load torque of motor 1 is the biggest at rest, there is a peak value on the current curve at the beginning of power supply, and the current begins to decrease after motor 1 rotating. There is another peak value of current before the moment when synchronized motion is carried out, and the current curve is fluctuating because the phase of motor 2 lags behind that of motor 1. The synchronization torque acts the load torque on motor 1, which increases the current of motor 1. After the moment synchronized motion is achieved, the current curve becomes smooth, which illustrates that the vibration system operates in the steady synchronization state because the load torque of motor 1 is close to constant. At about 30 s, the power of motor 2 is switched off. The current value of motor 1 increases with increasing phase difference, but the current value is still smooth under the synchronized motion. According to the change of current, we know the energy distribution of the vibration system. In other words, we can see how the power consumption changes after cutting off motor 2 by observing the change of current of motor 1.

As shown in Figure 10, , motor 1 first runs in the starting process, and motor 2 begins to run after 3 s. Figures 10(a) and 10(b) show that the synchronization velocity is about  r/min, and the phase difference is about . In the starting process, the phase difference oscillates between and because of the difference in rotational velocities of two motors. Similar to scheme in Figures 9(c)–9(f), Figure 10(c) and Figure 10(d) show the responses of the position for number 5, and Figure 10(e) and Figure 10(f) show the responses of the position for number 6. From comparison of the response of number 5 and number 6, it can be seen that the main motion type of the vibration system is the spatial circle motion. Due to the different type of the vibration system, the responses of number 6 in Figure 10 are greater than that in Figure 9.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 10: Experiment results with , , and : (a) rotational velocities; (b) the phase difference; (c) response in -direction for number 5; (d) response in -direction for number 5; (e) response in -direction for number 6; (f) response in -direction for number 6.

At about 27 s, the power of motor 2 is switched off, which breaks the original synchronization state between two ERs. Since the condition of the vibration synchronization transmission cannot be satisfied, the rotational velocity of motor 2 rapidly decreases, oscillation occurs in the phase difference, and the responses of number 5 and number 6 also show irregularity oscillation, as shown in Figure 9. After 3 s, the power of motor 2 is switched on. The vibration system again implements the synchronized motion. The new synchronization state is same with the old synchronization state.

From the above two groups of experiments, we know that the initial state of two motors has no effect on the synchronization state. When the synchronization torque is large enough, the vibration synchronization transmission can be achieved. As shown in Figure 4, the synchronization torque of is greater than that of . So the vibration system with can achieve the vibration synchronization transmission and that with cannot. When , the motion type of the vibration system is the spatial cone motion, while when , the motion type of the vibration system is the spatial circle motion, which comply with Figure 6. From Figure 7, we know that when , the synchronized motion of the vibration system is stable, which conforms to the responses of number 5 and number 6. The comparison between numeric results and experiments shows they are very good consistency. Based on experiments, the validity of theoretical analysis is proved, which helps the design of this type of vibration machine.

6. Conclusions

Based on the theoretical analysis and experiment result, the following remarks should be explained.

With introducing the average method of small parameters, synchronization of two ERs with the common rotational axis in the vibration system of the far-resonant spatial motion is investigated. Based on the dimensionless coupling equation, the synchronization condition of two ERs is obtained, and the stability condition of two ERs achieving the synchronized motion is derived, which satisfies Routh-Hurwitz criterion. The two conditions can be used to evaluate and discriminate whether the vibration system is able to carry out the synchronized motion.

The coupling dynamic characteristic, such as the choosing motion feature of the vibration system, is discussed by numerical method. Synchronization of two ERs depends on the coupling dynamic characteristic. The better the symmetry for two ERs, the greater the synchronization torque, the easier the achievement of the synchronized motion, and the stronger the synchronization ability. In addition, should be as small as possible for strengthening the stability of the synchronized motion.

From experiments, it can be seen that if , the spatial cone motion of the vibration system forces the phase difference to approach , and if , the spatial circle motion of the vibration system forces the phase difference to approach zero. In engineering, should be selected as far as possible. Only in this way we can obtain the circle motion of spatial for this type of vibration machine.

Comparing theoretical analysis, numerical results, and experiments, the feasibility of the theory method in this paper is proved.

Appendix

A. Rotation Matrix

According to Cartesian coordinates, the rotation matrix is [19]

Under the hypothesis of small angel and , the matrix can be simplified as

B. Dimensionless Parameters

Consider the following:

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

This work is supported by National Natural Science Foundation of China, Grants nos. 51305070 and 51375080.

References

B. C. Wen, J. Fan, C. Y. Zhao, and W. L. Xiong, Vibratory Synchronization and Controlled Synchronization in Engineering, Science Press, Beijing, China, 2009.

C. Huygens, Horologium Oscillatorium, F. Muguet, Paris, France, 1673.

I. I. Blekhman, “Self-synchronization of vibrators in some types of vibrational machines,” Inzhenerny Sbornik, vol. 16, pp. 49–72, 1953 (Russian).

I. I. Blekhman, A. L. Fradkov, H. Nijmeijer, and A. Y. Pogromsky, “On self-synchronization and controlled synchronization,” Systems & Control Letters, vol. 31, no. 5, pp. 299–305, 1997.

I. I. Blekhman, Synchronization in Science and Technology, ASME Press, New York, NY, USA, 1988.

I. I. Blekhman, A. L. Fradkov, O. P. Tomchina, and D. E. Bogdanov, “Self-synchronization and controlled synchronization: general definition and example design,” Mathematics and Computers in Simulation, vol. 58, no. 4–6, pp. 367–384, 2002.

I. I. Blekhman and N. P. Yaroshevich, “Extension of the domain of applicability of the integral stability criterion (extremum property) in synchronization problems,” Journal of Applied Mathematics and Mechanics, vol. 68, no. 6, pp. 839–846, 2004.

I. I. Blekhman and V. S. Sorokin, “On the separation of fast and slow motions in mechanical systems with high-frequency modulation of the dissipation coefficient,” Journal of Sound and Vibration, vol. 329, no. 23, pp. 4936–4949, 2010.

C. Y. Zhao, H. T. Zhu, R. Z. Wang, and B. C. Wen, “Synchronization of two non-identical coupled exciters in a non-resonant vibrating system of linear motion. Part I: theoretical analysis,” Shock and Vibration, vol. 16, no. 5, pp. 505–515, 2009.

C. Y. Zhao, H. T. Zhu, T. J. Bai, and B. C. Wen, “Synchronization of two non-identical coupled exciters in a non-resonant vibrating system of linear motion. Part II: numeric analysis,” Shock and Vibration, vol. 16, no. 5, pp. 517–528, 2009.

C. Y. Zhao, Q. H. Zhao, Y. M. Zhang, and B. C. Wen, “Synchronization of two non-identical coupled exciters in a non-resonant vibrating system of plane motion,” Journal of Mechanical Science and Technology, vol. 25, no. 1, pp. 49–60, 2011.

X.-L. Zhang, C.-Y. Zhao, and B.-C. Wen, “Theoretical and experimental study on synchronization of the two homodromy exciters in a non-resonant vibrating system,” Shock and Vibration, vol. 20, no. 2, pp. 327–340, 2013.

X. L. Zhang, B. C. Wen, and C. Y. Zhao, “Vibratory synchronization transmission of two exciters in a super-resonant vibrating system,” Journal of Mechanical Science and Technology, vol. 28, no. 6, pp. 2049–2058, 2014.

X. L. Zhang, B. C. Wen, and C. Y. Zhao, “Synchronization of three non-identical coupled exciters with the same rotating directions in a far-resonant vibrating system,” Journal of Sound and Vibration, vol. 332, no. 9, pp. 2300–2317, 2013.

J. M. Balthazar, J. L. P. Felix, and R. M. L. R. F. Brasil, “Some comments on the numerical simulation of self-synchronization of four non-ideal exciters,” Applied Mathematics and Computation, vol. 164, no. 2, pp. 615–625, 2005.

J. M. Balthazar, J. L. P. Felix, and R. M. L. R. F. Brasil, “Short comments on self-synchronization of two non-ideal sources supported by a flexible portal frame structure,” JVC/Journal of Vibration and Control, vol. 10, no. 12, pp. 1739–1748, 2004.

Y. Bonkobara, T. Ono, and T. Kondou, “Development of a generation mechanism of synchronous vibration suitable for hand-held vibrating tools: investigation of an impact model with two oscillators,” Journal of System Design and Dynamics, vol. 5, no. 6, pp. 1361–1375, 2011.

S. Rumyantsev and D. Tarasov, “Numerical simulation of non-linear dynamics of vibration transport machines in case of three independently rotating vibration exciters: recent advances in applied mathematics,” in Proceedings of the American Conference on Applied Mathematics (AMERICAN-MATH ’10), pp. 191–194, Cambridge, Mass, USA, January 2010.

C. Y. Zhao, H. T. Zhu, Y. M. Zhang, and B. C. Wen, “Synchronization of two coupled exciters in a vibrating system of spatial motion,” Acta Mechanica Sinica, vol. 26, no. 3, pp. 477–493, 2010.

Á. Miklós and Z. Szabó, “Vibrator with DC motor driven eccentric rotors,” Periodica Polytechnica, Mechanical Engineering, vol. 56, no. 1, pp. 49–53, 2012.

Á. Miklós and Z. Szabó, “Mechanical synchronization in dual-rotor vibroactuator,” Proceedings in Applied Mathematics and Mechanics, vol. 13, no. 1, pp. 41–42, 2013.

A. Miklós and Z. Szabó, “Simulation and experimental validation of the dynamical model of a dual-rotor vibrotactor,” Journal of Sound and Vibration, vol. 334, pp. 98–107, 2015.

Y. Ruan and B. S. Chen, Control Systems of Electric Derives-Motion Control Systems, China Machine Press, Beijing, China, 2010 (Chinese).