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

Authors: Fang Pan; Hou Yongjun; Dai Liming; Du Mingjun
Journal: Shock and Vibration
 (2018)

---

## Abstract

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 simplified mechanical 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”.


---

## Body

## 1. 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–5], nonlinear coupled chaotic systems [6, 7], coupled pendulum system [8–12], and rotor system [1, 13–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 self-synchronization 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 in-phase 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 pendulum-rotor 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–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–11] and classical rotor systems presented in papers [12–20, 23–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.


## 2. Dynamic Model

A simplified dual-pendulum-rotor system is shown in Figure1. This system consists of a rigid oscillating body of mass m0 [Kg] elastically supported via a linear damping spring with stiffness kx [N/m] and damping fx [Ns/m]. Unbalanced rotor i actuated by an asynchronous motor is modelled by a point mass mi [Kg] (for i=1,2) and attached at the end of a massless rod of length r [m]; the rotation angle of the rotors is defined by φi∈S1 (for i=1,2) in [rad]. It should be noted that all the motors are installed rigid vibrating body by linear torsion springs with stiffness kφN/rad and damping fφ[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 φj∈S1 (for j=3,4) in [rad]. Therefore, the oscillating characteristics of the motors are just like a pendulum or pendula with length l and lumped mass (equal to the mass of the motors) mj [Kg] (for j=3,4). The angular installation of the motors is denoted by β∈S1 in [rad]; x represents oscillation of centroid of the system in x-direction, respectively. σ is the rotation direction of the rotors; l [m] is distance between pivot of the motor and centroid O of the system.Figure 1
The dual-pendulum-rotor system.In rotating coordinateo′x″y″, the rotating centers of the rotors Φ1″, Φ2″ can be expressed as(1)Φ1″=σrcos⁡φ1-lcos⁡φ3+βrsin⁡φ1+lsin⁡φ3+β,Φ2″=rcos⁡φ2+lcos⁡φ4+βrsin⁡φ2+lsin⁡φ4+β.And the rotating center of the motorsΦ3″ and Φ4″ is represented as(2)Φ3″=-lcos⁡β+φ3lsin⁡β+φ3,Φ4″=lcos⁡β+φ4lsin⁡β+φ4.In fixed coordinateoxy, coordinates Φ1,Φ2,Φ3,Φ4 can be obtained by(3)Φi=Φ0+Φi″.i=1,2,3,4whereΦ0 is displacement vector of the vibrating base, Φ0=x,0T.According to kinetic theory, kinetic energy T of the dual-pendulum-rotor system should be written as follows:(4)T=12m0x˙+12∑i=14miΦ˙iTΦ˙i+12∑i=12Jiφ˙i2.J 1 and J2 are the rotational inertia of the rotors, and symbols •˙ and •¨ denote d•/dt and d2•/dt2, respectively.In addition, potential energyV of the dual-pendulum-rotor system can be expressed by(5)V=12kxx2+12kφφ32+12kφφ42Moreover, dissipated energyD of the system can be obtained by(6)D=12fxx˙2+12f1φ˙12+12f2φ˙22+12f3φ˙32+12f4φ˙42According to the Lagrange equation(7)ddt∂T∂q˙-∂T-V∂q+∂D∂q=Qi,the dynamic equation of the dual-pendulum-rotor system can be derived. In the system, coordinate matrixq=x,φ1,φ2,φ3,φ4T is considered as the generalized coordinate matrix; meanwhile, the generalized force matrix of the system is assumed as(8)QxQφ1Qφ2Qφ3Qφ4T=0Me1-Re1Me2-Re200TSubstituting (4), (5), (6), and (8) into (7) and, meanwhile, considering m1,m2≪m0,∑k=04mk=M. As J1 and J2 are sums of rotational inertia of motor’s rotor about spin axis and that of eccentric lump, compared with m1r2 and m1r2, the rotational inertia of motor’s rotor about spin axis is too small to be neglected; i.e., J1≈m1r2,J2≈m2r2, and the dynamic equation of the dual-pendulum-rotor system is written by (9a)J1φ¨1=Me1-Re1+σm1rsinφ1x¨-m1rlφ¨3cos⁡β+σφ1+m1rlφ˙32sin⁡β+σφ1(9b)J2φ¨2+J2φ¨2=Me2-Re2+m2rsin⁡φ2x¨-m2rlφ¨4cos⁡β-φ2+m2rlφ˙42sin⁡β-φ2(9c)Mx¨+fxx˙+kxx=σm1rφ¨1sin⁡φ1+φ˙12cos⁡φ1+m2rφ¨2sin⁡φ2+φ˙22cos⁡φ2-m1+m3lφ¨3sin⁡β+φ˙32cos⁡β+m2+m4lφ¨4sin⁡β+φ˙42cos⁡β(9d)m1l2+m3l2φ¨3+fφφ˙3+kφφ3=-m1+m3lsin⁡βx¨-m1lrφ¨1cos⁡β+σφ1-σφ˙12sin⁡β+σφ1(9e)m2l2+m4l2φ¨4+fφφ˙4+kφφ4=m2+m4lsin⁡βx¨-m2lrφ¨2cos⁡β-φ2+φ˙22sin⁡β-φ2where(10)σ=1,thetworotorrotateinthesamedirection,-1,thetworotorrotateintheoppositedirection.It can be seen that the terms ofφ˙32 and φ˙42, nonlinear terms owing to the small periodic vibration in synchronous state, are included in (9c). The terms related to φ˙32 and φ˙42 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]:(11a)J1φ¨1=μΦ1(11b)J2φ¨2=μΦ2(11c)Mx¨+kxx=σm1rφ¨1sin⁡φ1+φ˙12cos⁡φ1+m2rφ¨2sin⁡φ2+φ˙22cos⁡φ2-m1+m3lφ¨3sin⁡β+m2+m4lφ¨4sin⁡β-μfxx˙+m1+m3lφ˙32cos⁡β-m2+m4lφ˙42cos⁡β(11d)m1l2+m3l2φ¨3+kφφ3=-m1+m3lsin⁡βx¨-m1lrφ¨1cos⁡β+σφ1-μfφφ˙3-σm1lrφ˙12sin⁡β+σφ1(11e)m2l2+m4l2φ¨4+kφφ4=m2+m4lsin⁡βx¨-m2lrφ¨2cos⁡φ2-β-μfφφ˙4+m2lrφ˙22sin⁡φ2-βwhere(12a)μΦ1=Me1-Re1+σm1rsin⁡φ1x¨-m1rlφ¨3cos⁡β+σφ1+m1rlφ˙32sin⁡β+σφ1(12b)μΦ2=Me2-Re2+m2rsin⁡φ2x¨-m2rlφ¨4cos⁡β-φ2+m2rlφ˙42sin⁡β-φ2Based on Zhao’s method [15], the electromagnetic torque of the asynchronous motor is linearized at stable point φ˙i=ωs when the two motors are synchronously operated; i.e.,(13)Mei=npLmi2US02Lsi2ωmRriωs-npωmi=1,2whereLmi is coefficient of mutual asynchronous in theith motor; Lsi is theith coefficient of stator inductance in theith motor; np is pole-pairs of the motors; ωm is synchronous velocity of the motors; Rri is rotor resistance in theith motor; US0 is amplitude of stator voltage.


## 3. 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,x,φ3, and φ4. Neglecting the terms of small parameters, introduce the dimensionless parameter as follows:(14)rl=rl,mM=rm,mim=ηii=1,2,3,…,4,ωx=kxM,ωφ=kφm1l2+m3l2,ω~φ=kφ/m1l2+m3l2,f1t=rmrlση1cos⁡φ1+η2cos⁡φ2,f2t=η1rlη1+η3sin⁡φ1+σβ,f3t=η2rlη2+η4sin⁡φ2-β.In addition, velocitiesφ˙1 and φ˙2 of the two motors approached to a constant are represented by ωm 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:(15a)x¨+ωx2x=ωm2lf1t-η1+η3rmlsin⁡βφ¨3+η2+η4rmlsin⁡βφ¨4,(15b)φ¨3+ωφ2φ3=-sin⁡βlx¨+ωm2f2t,(15c)φ¨4+ω~φ2φ4=sin⁡βlx¨+ωm2f3t.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 (16a)s2X+ωx2X=ωm2lF1s-η1+η3rmls2sin⁡βϕ3+η2+η4rmls2sin⁡βϕ4(16b)s2ϕ3+ωφ2ϕ3=-s2sin⁡βlX+ωm2F2s(16c)s2ϕ4+ωφ2ϕ4=s2sin⁡βlX+ωm2F3sFrom formulas (16a)-(16c), the approximate solutions of the system in complex domain can be given by(17a)X=ωm2lG1sF1s-η1+η3ωm2rmls2sin⁡βG1sF2s+η2+η4ωm2rmls2sin⁡βG1sF3s(17b)ϕ3=-s2sin⁡βωm2G1sF1s+ωm2G2s+ωm2η1+η3rms4sin2βG1sG2sF2s-ωm2η2+η4rms4sin2βG1sG2sF3s(17c)ϕ4=ωm2s2sin⁡βG1sF1s-ωm2η1+η3rms4sin2βG1sG2sF2s+ωm2G2s+ωm2η2+η4rms4sin2βG1sG2sF3swhereG1s=s2+ωx2s2+ωφ2-η1+η3+η2+η4rms4sin2β,G2s=s2+ωφ2.Then, the numerators and denominators of (17a), (17b), (17c) in the right are divided by divisor ωφ4ωx2; moreover, introduce frequency ratios(18)nx=ωmωx,nφ=ωmωφinto (17a)-(17c). In this case, the coefficients of spring stiffness kx and kφ are converted into the frequency ratios nx and nφ 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(19a)x=μ11lf1t+μ12η1+η3lrmf2t+μ13η2+η4lrmf3t,(19b)φ3=μ21f1t+μ22η1+η3rmf2t+μ23η2+η4rmf3t,(19c)φ4=μ31f1t+μ32η1+η3rmf2t+μ33η2+η4rmf3t,where(20)μ11=nx21-nφ2g,μ12=μ21=nx2nφ2sin⁡βg,μ13=μ31=-nx2nφ2sin⁡βg,μ22=μ33=nx2nφ4sin2β1-nφ2g+nφ21-nφ2η1+η3rm,μ23=μ32=-nx2nφ4sin2β1-nφ2g,g=1-nx21-nφ2-η1+η3+η2+η4rmnx2nφ2sin2βIn light of functionsf1t,f2t, and f3t in (14), the approximate solutions of the system can be rearranged as the following:(21a)x=σμ11η1rrmcos⁡φ1+μ11η2rrmcos⁡φ2+μ12η1rrmsin⁡φ1+σβ+μ13η2rrmsin⁡φ2-β,(21b)φ3=σμ21η1rmrlcos⁡φ1+μ21η2rmrlcos⁡φ2+μ22η1rlrmsin⁡φ1+σβ+μ23η2rlrmsin⁡φ2-β,(21c)φ4=μ31rmrlση1cos⁡φ1+η2cos⁡φ2+μ32η1rlrmsin⁡φ1+σβ+μ33η2rlrmsin⁡φ2-β.


## 4. Synchronization

Based on the approximate solutions above, the synchronization and stability of the system can be determined with Poincaré method [20, 25–27]. In addition, the phase angle of the rotors can be assumed by (22)φ1=ωmt+α1φ2=ωmt+α2Substituting (22) into (21a), (21b), (21c), the approximate solutions of the system can be written by (23a)x=μ11rrmση1cos⁡ωmt+α1+η2cos⁡ωmt+α2+μ12rrmη1sin⁡ωmt+α1+σβ+μ13rrmη2sin⁡ωmt+α2-β,(23b)φ3=μ21rmrlση1cos⁡ωmt+α1+η2cos⁡ωmt+α2+μ22rlrmη1sin⁡ωmt+α1+σβ+μ23rlrmη2sin⁡ωmt+α2-β,(23c)φ4=μ31rmrlση1cos⁡ωmt+α1+η2cos⁡ωmt+α2+μ32η1rlrmsin⁡ωmt+α1+σβ+μ33η2rlrmsin⁡ωmt+α2-β.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(24)α=α1-α2The second derivatives ofx,φ3, and φ4 with respect to timet can be obtained. Substituting x¨,φ¨3, and φ¨4 into (12a) and (12b), and then integrating and averaging these equations related tot over period T, respectively, we can obtain (25)P1=1T∫0TμΦ1dt=M-e1-R-e1+12η1η2mr2ωm2rm-σμ11sin⁡α+σμ12cos⁡α+β+μ21cosα+σβ-μ23sinα+σ+1β(26)P2=1T∫0TμΦ2dt=M-e2-R-e2+12η1η2mr2ωm2rm-σμ11sin⁡α+σμ12cosα+β+μ21cosα+σβ-μ23sinα+σ+1βP 1 and P2 represent the averaging residual torque of the first and second motor in single periodT, 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., (27)P1≈P2≈0Substituting (27) into (25) and (26), the balanced torques of the motor are written by (28a)M-e1-R-e1=12η1η2mr2ωm2rmσμ11sin⁡α-σμ12cos⁡α+β-μ21cos⁡α+σβ+μ23sin⁡α+σ+1β(28b)M-e2-R-e2=12η1η2mr2ωm2rm-σμ11sin⁡α+σμ12cos⁡α+β+μ21cos⁡α+σβ-μ23sin⁡α+σ+1βM - e i , R - e i   i = 1,2 represent average electromagnetic and friction torque in the motors over period T,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 ΔM between the motors can be expressed as(29)ΔM=M-e1-R-e1-M-e2-R-e2=η1η2mr2ωm2rmσμ11sin⁡α-σμ12cos⁡α+β-μ21cos⁡α+σβσμ12cos⁡α+β+μ21cos⁡α+σβAs the parameter of the two asynchronous motors is identical, the differences of output torque should be zero. Therefore, to implement synchronous operation, the values of the parameters in the system must satisfy the following equation:(30)-σμ11sin⁡α+σμ12cos⁡α+β+μ21cos⁡α+σβ-μ23sin⁡α+σ+1β=0which is called balance equation of synchronous state. Obviously, the phase difference is influenced by the rotation direction of the motors, the coupling coefficients between rotor and pendula, and the installation location of the motors. According to (30), we can determine phase difference α with numerical computations. With trigonometric function, the equation above can be expanded as (31)cot⁡α=σμ11+σμ12sin⁡β+σμ21sin⁡β+μ23cos⁡σ+1βσμ12cos⁡β+μ21cos⁡β-μ23sin⁡σ+1βIf some solutions of phase differenceα exist in (31), the denominator of this equality in the right should be nonzero. Define the denominator of this equality with a symbol D; i.e.,(32)D=σμ12cos⁡β+μ21cos⁡β-μ23sin⁡σ+1βIt can be known that denominatorD 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 D. Therefore, we will discuss zero solution of parameter D related to installation angle and coupling coefficients. There are two cases causing the zero solution of denominator D.( 1 ) Parameter β represents installation angle of the motors. Thus, parameter β with a particular value leads toD being zero. The first "critical point" is (33)β=kπ-π,k=1,2,3….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 nx and nφ with particular values lead to D being zero. The second "critical point" is (34)nx=1ornφ=1.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 ratiosnx and nφ are the function related to stiffness coefficients kx and kφ 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 ratiosnx and nφ should be far away from 1. In light of different value of nφ and nφ, the coupling type of the system can be defined as follows:Type 1: system of before-resonance coupled before-resonance (0.1<nx<1 and 0.1<nφ<1).Type 2: system of after-resonance coupled before-resonance (0.1<nx<1 and 1.0<nφ<7, or 1.0<nx<7 and 0.1<nφ<1).Type 3: system of after-resonance coupled after-resonance (1.0<nφ<7 and 1.0<nx<7).Eliminating the two kinds of "critical point", the phase difference between the two rotors can be calculated:(35)α=arccot⁡σμ11+σμ12sin⁡β+σμ21sin⁡β+μ23cos⁡σ+1βσμ12cos⁡β+μ21cos⁡β-μ23sin⁡σ+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 nx, nφ, β, and rm, and these parameters may influence the value of phase difference α.


## 5. 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 (36)∂P1-P2∂α-χ=0,phase differenceα∗ is stable under the condition that the value of χ is negative.Thus, the residual torque difference can be obtained by(37)P1-P2=μΦ1=M-e1-R-e1-M-e2-R-e2+η1η2mr2ωm2rm-σμ11sin⁡α+σμ12cos⁡α+β+μ21cos⁡α+σβ-μ23sin⁡α+σ+1βAsM-ei,R-eii=1,2 are irrelevant to phase difference α, in light of (33), the criterion of synchronous stability is expressed(38)χ=∂P1-P2∂α=η1η2mr2ωm2rm-σμ11cos⁡α-σμ12sin⁡α+β-μ21sin⁡α+σβ-μ23cos⁡α+σ+1β<0.Rearranging the equation above, the criterion of synchronous stability is simplified by(39)σμ11cos⁡α+σμ12sin⁡α+β+μ21sin⁡α+σβ+μ23cos⁡α+σ+1β>0The 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.


## 6. 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 (nxnφ), mass ratios (rm, ηi), 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.Table 1
Parameter values.
(a)
Parameter values for system equation (6)
Unbalanced rotor fori=1,2 Vibroplatform Pendulum rod Asynchronous motor m i = 2  [ k g ] m 0 = 100  [ k g ] l = 0.3  [ m ] m 3 = m 4 = 14  [ k g ] r = 0.05  [ m ] k x = 246490000 - 50307  [ N / m ] k φ = 4980 - 24402510  [ N m / r a d ] L m i = 0.13  [ H ] ω m = 152 - 157  [ r a d / s ] f x = 1064 [ N s / m ] f φ = 15  [ N m / ( r a d / s ) ] L s i = 0.1  [ H ] - - β = 0 - 5 π / 12  [ r a d ] n p = 2 - - - R r i = 0.54  [ Ω ] - - - U S 0 = 220  [ V ]  (b)
Parameter values according to dimensionless equations (11a), (11b), (11c), (11d), (11e)
η 1 = 0.2,1 η 2 = 1 η 3 = η 4 = 12 σ = - 1,1 r m = 0.02 n x = 0.1 ~ 7 n φ = 0.1 ~ 7


### 6.1. The Rotors Operated in Opposite Direction

The stable phase difference is determined by considering the different values of parametersβ, ηii=1,2, nx, and nφ. 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+μ23sin⁡α=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+μ23sin⁡α>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 rm 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 nx and nφ; 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 rm is fixed at 0.02, and the frequency ratios nx and nφ 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 nx and nφ>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 nx and nφ<1.5. Finally, the effect of mass ratios rm, η1, and η2 on the phase difference needs to be discussed. Thus, the parameter values in Figure 2(c) are identical to Figure 2(a) except for rm. 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), 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.Figure 2
Stable phase difference whenσ=-1.
(a)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   β = π / 3 (b)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   n x ≈ n φ (c)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.04 ,   β = π / 3 (d)
η 1 = 0.2 ,   η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   n x ≈ n φ


### 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 nx and nφ. The obtained results reveal that, in the interval of 0.1<nx<1 and 0.1<nφ<1, synphase synchronization is implemented between the rotors; in the interval of 0.1<nx<1 and 1<nφ<7, antiphase synchronization is implemented between the rotors; in the interval of 1<nx<7, 0.1<nφ<1 and 1<nx<7, 1<nφ<7, both antiphase and synphase synchronization collectively remained; in the interval of 1<nx<7 and for 1<nφ<7, the two rotors will be synchronously operated in synphase state. 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 rm, η1, and η2 on the phase difference is discussed. Thus, the parameters in Figure 3(c) are identical with Figure 3(a) except for rm, 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), 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.Figure 3
Stable phase difference whenσ=1.
(a)
η 1 = η 2 = 1 , η 3 = η 4 = 12 ,   r m = 0.02 ,   β = π / 3 (b)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   n x ≈ n φ (c)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.04 ,   β = π / 3 (d)
η 1 = 0.2 ,   η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   n x ≈ n φ


## 6.1. The Rotors Operated in Opposite Direction

The stable phase difference is determined by considering the different values of parametersβ, ηii=1,2, nx, and nφ. 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+μ23sin⁡α=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+μ23sin⁡α>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 rm 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 nx and nφ; 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 rm is fixed at 0.02, and the frequency ratios nx and nφ 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 nx and nφ>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 nx and nφ<1.5. Finally, the effect of mass ratios rm, η1, and η2 on the phase difference needs to be discussed. Thus, the parameter values in Figure 2(c) are identical to Figure 2(a) except for rm. 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), 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.Figure 2
Stable phase difference whenσ=-1.
(a)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   β = π / 3 (b)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   n x ≈ n φ (c)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.04 ,   β = π / 3 (d)
η 1 = 0.2 ,   η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   n x ≈ n φ


## 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 nx and nφ. The obtained results reveal that, in the interval of 0.1<nx<1 and 0.1<nφ<1, synphase synchronization is implemented between the rotors; in the interval of 0.1<nx<1 and 1<nφ<7, antiphase synchronization is implemented between the rotors; in the interval of 1<nx<7, 0.1<nφ<1 and 1<nx<7, 1<nφ<7, both antiphase and synphase synchronization collectively remained; in the interval of 1<nx<7 and for 1<nφ<7, the two rotors will be synchronously operated in synphase state. 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 rm, η1, and η2 on the phase difference is discussed. Thus, the parameters in Figure 3(c) are identical with Figure 3(a) except for rm, 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), 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.Figure 3
Stable phase difference whenσ=1.
(a)
η 1 = η 2 = 1 , η 3 = η 4 = 12 ,   r m = 0.02 ,   β = π / 3 (b)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   n x ≈ n φ (c)
η 1 = η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.04 ,   β = π / 3 (d)
η 1 = 0.2 ,   η 2 = 1 ,   η 3 = η 4 = 12 ,   r m = 0.02 ,   n x ≈ n φ


## 7. Sample Verifications

### 7.1. Forσ=-1, nx=6, nφ=6, η1=η2=1, η3=η4=12, rm=0.02, and β=π/3

This subsection refers to the case that the system of after-resonance is coupled with the after-resonance; i.e.,nx=6, nφ=6. The mass of the two rotors operated in opposite direction is identical; i.e., η1=η2=1. In the simulation model, the coefficients of the spring stiffness are kx=65969[N/m] and kφ=1053[N·m/rad], and the other parameter is identical with Table 1(a). The numerical 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) and 3(b) is equal to 0 [rad]; here, the stable phase difference is fluctuated near 0 [rad]).Figure 4
The dynamic characteristics when the rotors rotate in the opposite direction.
(a)
Velocity (b)
Phase difference (c)
Electromagnetic torque (d)
Vibration response


### 7.2. Forσ=1, nx=6, nφ=0.3, η1=η2=1, η3=η4=12, rm=0.04, and β=π/3

Simulation results obtained by (9a), (9b), (9c), (9d), (9e) are shown in Figure 5 whenσ=1,nx=6,nφ=0.3,η1=η2=1,η3=η4=12,rm=0.04, and β=π/3. Here, the coefficients of spring stiffness are kx=65969[N/m] and kφ=976100[N·m/rad], and the rotors are operated in the same direction. The coupling type of the system belongs to type 2. When the two motors are supplied by the electric source at the same time, the angular accelerations of the two rotors are compatible (in Figure 5(a)). The reason is that the inertia moments of the rotors are identical and the spring stiffness is extremely stronger. During the staring process of the system, the velocity difference exists between the two rotors, which leads to the phase difference instability, shown as in Figure 5(b). However, when the angular velocities of the motors reach the rated speed and the motor oscillate steadily, the synchronization phenomenon occurs. It should be noted that the velocity fluctuation of the motors is smaller than that of Section 7.1, because the vibration of the motors is feebler. At this moment, the average coupling torques (in Figure 5(c)), making the phase difference α stabilized at -0.8984 [rad], are approximated to 4.08[N·m]. In this case, the two motors rotate stably in the synphase synchronization, and the synchronous velocity is 157 [rad/s]. From Figure 5(d) it follows the displacements of the vibrating body and the two motors. It can be seen that the displacement responses of the vibrating body and the two motors are stable, and the amplitudes of them are 1.5×10-3[m], 4.1×10-4[rad] and 7.3×10-4[rad], respectively. Comparing simulation results with Figure 3(b), it should be noted that the value of the stable phase difference is in agreement with the results obtained for the case of the theoretical solutions (i.e., the stable phase difference of such parameters in Figure 3(a) is equal to -1.08 [rad]; here, the stable phase difference is equal to -0.8984 [rad]).Figure 5
The dynamic characteristics when the rotors rotate in the same direction.
(a)
Velocity (b)
Phase difference (c)
Electromagnetic torque (d)
Vibration response


## 7.1. Forσ=-1, nx=6, nφ=6, η1=η2=1, η3=η4=12, rm=0.02, and β=π/3

This subsection refers to the case that the system of after-resonance is coupled with the after-resonance; i.e.,nx=6, nφ=6. The mass of the two rotors operated in opposite direction is identical; i.e., η1=η2=1. In the simulation model, the coefficients of the spring stiffness are kx=65969[N/m] and kφ=1053[N·m/rad], and the other parameter is identical with Table 1(a). The numerical 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) and 3(b) is equal to 0 [rad]; here, the stable phase difference is fluctuated near 0 [rad]).Figure 4
The dynamic characteristics when the rotors rotate in the opposite direction.
(a)
Velocity (b)
Phase difference (c)
Electromagnetic torque (d)
Vibration response


## 7.2. Forσ=1, nx=6, nφ=0.3, η1=η2=1, η3=η4=12, rm=0.04, and β=π/3

Simulation results obtained by (9a), (9b), (9c), (9d), (9e) are shown in Figure 5 whenσ=1,nx=6,nφ=0.3,η1=η2=1,η3=η4=12,rm=0.04, and β=π/3. Here, the coefficients of spring stiffness are kx=65969[N/m] and kφ=976100[N·m/rad], and the rotors are operated in the same direction. The coupling type of the system belongs to type 2. When the two motors are supplied by the electric source at the same time, the angular accelerations of the two rotors are compatible (in Figure 5(a)). The reason is that the inertia moments of the rotors are identical and the spring stiffness is extremely stronger. During the staring process of the system, the velocity difference exists between the two rotors, which leads to the phase difference instability, shown as in Figure 5(b). However, when the angular velocities of the motors reach the rated speed and the motor oscillate steadily, the synchronization phenomenon occurs. It should be noted that the velocity fluctuation of the motors is smaller than that of Section 7.1, because the vibration of the motors is feebler. At this moment, the average coupling torques (in Figure 5(c)), making the phase difference α stabilized at -0.8984 [rad], are approximated to 4.08[N·m]. In this case, the two motors rotate stably in the synphase synchronization, and the synchronous velocity is 157 [rad/s]. From Figure 5(d) it follows the displacements of the vibrating body and the two motors. It can be seen that the displacement responses of the vibrating body and the two motors are stable, and the amplitudes of them are 1.5×10-3[m], 4.1×10-4[rad] and 7.3×10-4[rad], respectively. Comparing simulation results with Figure 3(b), it should be noted that the value of the stable phase difference is in agreement with the results obtained for the case of the theoretical solutions (i.e., the stable phase difference of such parameters in Figure 3(a) is equal to -1.08 [rad]; here, the stable phase difference is equal to -0.8984 [rad]).Figure 5
The dynamic characteristics when the rotors rotate in the same direction.
(a)
Velocity (b)
Phase difference (c)
Electromagnetic torque (d)
Vibration response


## 8. 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 aero-power 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.



---

*Source: 9824631-2018-07-09.xml*