Synchronization behavior of a weakly damped far-resonance vibrating system

Chen, Bang; Xia, Xiao’ou; Wang, Xiaobo

doi:10.21595/vp.2018.20199

Vibroengineering Procedia

Browse Procedia

Published: 19 October 2018

Check for updates

Synchronization behavior of a weakly damped far-resonance vibrating system

Bang Chen¹

Xiao’ou Xia²

Xiaobo Wang³

^{1, 2}School of Mechanical Engineering, University of Science and Technology Beijing, Beijing, China

³BGRIMM Technology Group, Beijing, China

Corresponding Author:

Bang Chen

Cite the article Download PDF

Downloads 1395

Abstract

In order to reveal synchronization characteristics of a weakly damped system with two rotors mounted on different vibrating bodies, we propose a simplified physical model. Vibration of the system is discussed by the average method, which can separate fast motions (high frequency) from slow motions (low frequency). Theoretical research shows that vibration torque is the key factor to balance the energy distribution between rotors. For the system with rotational frequencies larger than the nature frequencies, the coupling characteristic frequency or characteristic frequency curve should be considered. As the coupling frequency is close to the characteristic frequency, or the vibration state is close to the characteristic frequency curve, self-synchronization of two rotors can be obtained easily.

1. Introduction

The so-called self-synchronization phenomenon corresponds to the consistency or certain relationship of systems’ parameters caused by their internal couplings and has been widely involved in non-linear vibration, hydraulic [1-3], electromechanical coupling, automatic control theory and other fields [4-8].

Huygens was the first person who observed the synchronization of pendulum clocks in the 17th Century. In [4, 5], Czolczynski et al. presented different synchronous behaves of two or $n$ pendula installed on a frame. The self-synchronization theory of rotors was developed by Bleckman [1, 2] with averaging method in the middle of the 20th century. Wen and Zhao et al. [7] modified the averaging method and proposed the average method with two small parameters. Zhang [8] deduced the synchronization condition and the synchronization stability for the vibrating system with three rotors. Hou and Fang [9, 10] investigated a vibrating screen based on the model of a rotor-pendulum system.

The above researches are mostly focused on the synchronization of pendula or rotors installed on the same vibrating frame. In this paper, we propose a vibrating system with two rotors mounted on two different vibrating bodies.

2. Dynamical equations of the vibrating system

As shown in Fig. 1, two rotors are mounted on different vibrating bodies. The vibrating body ( $M_{i}$ ) ( $i =$ 1, 2) can move in horizontal direction ( $x_{i}$ ) and is installed on the foundation by the spring. The two bodies are connected by a coupling spring. Counterclockwise direction is taken as positive. Inertia moment and eccentricity of the rotor on its mass center are given by $j_{i}$ and $r_{i}$ . Other variables are show in Fig. 1. In this paper, synchronization of rotors is analyzed in a non-resonant vibrating system, in which rotation frequencies of rotors are larger than nature frequencies of vibrating bodies. The system is denoted as after-resonance system. We assume that ${\dot{φ}}_{1}$ , ${\dot{φ}}_{2} > 2 ω_{1}$ , $2 ω_{2}$ , where $ω_{1} = \sqrt{k_{1} / M_{1}}$ , $ω_{2} = \sqrt{k_{2} / M_{2}}$ .

Fig. 1Simplified model of the system

The electromagnetic torque and resistance moment of the driving motor $i$ are $T_{e i}$ and $T_{f i}$ ( $T_{f i} = f_{r i} ∙ {\dot{φ}}_{i}$ ) respectively. When ${\dot{φ}}_{i}$ ( $i =$ 1, 2) fluctuates near the frequency $ω_{n}$ , the influence of electromagnetic leakage can be neglected, and the driving force of induction motor can be linearized as $T_{e i} = n L_{m i}^{2} U_{0}^{2} (ω_{s i} - n ω_{n}) / ({L_{s i}^{2} ω_{s i} R_{r i} ω}_{s i})$ , where $n_{i}$ , $L_{m i}$ , $ω_{s i}$ , $L_{s i}$ and $R_{r i}$ ( $i =$ 1, 2) correspond to the pole number, mutual inductance, synchronous speed, stator inductance and rotor resistance of the motor; $U_{0}$ is the voltage amplitude. As self-synchronization of rotors is achieved, speed fluctuations of rotors are small [1, 2]. Therefore, small variables can be neglected. Introducing $ω_{p} = \sqrt{k_{p} / M_{1}}$ , $η = M_{2} / M_{1}$ , $J_{1} = j_{1} + m_{1} r_{1}^{2}$ and $J_{2} = j_{2} + m_{2} r_{2}^{2}$ , we have:

1

{\ddot{x}}_{1} + ω_{1}^{2} x_{1} + η ω_{p}^{2} (x_{1} - x_{2}) = \frac{m_{1}}{M_{1}} r_{1} {\dot{φ}}_{1}^{2} \cos φ_{1},

{\ddot{x}}_{2} + ω_{2}^{2} x_{2} + η ω_{p}^{2} (x_{2} - x_{1}) = \frac{m_{2}}{M_{2}} δ r_{2} {\dot{φ}}_{2}^{2} \cos φ_{2},

J_{1} {\ddot{φ}}_{1} = T_{e 1} - T_{f 1} + m_{1} r_{1} {\ddot{x}}_{1} \sin φ_{1},

J_{2} {\ddot{φ}}_{2} = T_{e 2} - T_{f 2} + δ m_{2} r_{2} {\ddot{x}}_{2} \sin φ_{2} .

The synchronous speed of two rotors is denoted by $ω_{n}$ . When self-synchronization of rotors is achieved [1, 3], phases of the two rotors can be denoted as $φ_{1} = ω_{n} t + α_{1}$ , $φ_{2} = ω_{n} t + α_{2}$ , where $α_{1}$ and $α_{2}$ are slowly-varying parameters. From Eq. (1), we obtain:

2

x_{1} = μ_{11} \cos φ_{1} + μ_{12} \cos φ_{2},

x_{2} = μ_{21} \cos φ_{1} + μ_{22} \cos φ_{2},

where:

μ_{11} = \frac{m_{1} r_{1} ω_{n}^{2} (ω_{2}^{2} + ω_{p}^{2} - ω_{n}^{2})}{N}, μ_{12} = \frac{δ m_{2} r_{2} ω_{n}^{2} ω_{p}^{2}}{N},

μ_{21} = \frac{m_{1} r_{1} ω_{n}^{2} ω_{p}^{2}}{N}, μ_{22} = \frac{δ m_{2} r_{2} ω_{n}^{2} (ω_{1}^{2} + η ω_{p}^{2} - ω_{n}^{2})}{(η N)},

N = M_{1} [(ω_{1}^{2} + η ω_{p}^{2} - ω_{n}^{2}) (ω_{2}^{2} + ω_{p}^{2} - ω_{n}^{2}) - η ω_{p}^{4}] .

$μ_{11}$ , $μ_{12}$ , $μ_{21}$ and $μ_{22}$ show the coupling effects in the system. Rotors are driven by motors, and the resistance is approximately proportional to its speed. The average values of resultant torques of rotors are denoted by $P_{1}$ , $P_{2}$ . As the system is stable, we have:

3

P_{1} = \int_{t}^{t + T} T_{e 1} - T_{f 1} + m_{1} r_{1} {\ddot{x}}_{1} \sin φ_{1} d t = T_{e 1} (ω_{n}) - T_{f 1} (ω_{n}) - \frac{1}{2} T_{v},

P_{2} = \int_{t}^{t + T} T_{e 2} - T_{f 2} + δ m_{2} r_{2} {\ddot{x}}_{2} \sin φ_{2} d t = T_{e 2} (ω_{n}) - T_{f 2} (ω_{n}) + \frac{1}{2} T_{v},

where $T_{v} = δ m_{2} r_{2} μ_{21} ω_{n}^{2} \sin (α_{1} - α_{2}) = m_{1} r_{1} μ_{12} ω_{n}^{2} \sin (α_{1} - α_{2})$ is the vibration torque (VT). Introducing the variable substitutions:

∆ α = α_{1} - α_{2}, ∆ T_{e} = T_{e 1} (ω_{n}) - T_{e 2} (ω_{n}) ∆ T_{f} (ω_{n}) = T_{f 1} (ω_{n}) - T_{f 2} (ω_{n}),

\sin ∆ α = \frac{[∆ T_{e} (ω_{n}) - ∆ T_{f} (ω_{n})]}{(m_{1} r_{1} μ_{12} ω_{n}^{2})},

can be obtained. Thus, the synchronization condition can be expressed as follows [1]:

4

|[∆ T_{e} (ω_{n}) - ∆ T_{f} (ω_{n})] ∙ \frac{1}{δ m_{2} r_{2} μ_{21} ω_{n}^{2}}| \leq 1 .

The stability criterion of the synchronous state can be discussed based on Lyapunov stability theory, it can be deduced as:

5

\frac{δ}{(ω_{1}^{2} + η ω_{p}^{2} - ω_{n}^{2}) (ω_{2}^{2} + ω_{p}^{2} - ω_{n}^{2}) - η ω_{p}^{4}} \cos ∆ α > 0 .

3. Discussions of theoretical results

In this paper, two rotors rotate in the same direction, that is, $δ =$ 1. In this paper, $L_{m i} =$ 0.14 H, $L_{s i} =$ 0.12 H, $R_{r i} =$ 0.6 $Ω$ , $ω_{s i} =$ 314 rad/s, $U_{0} =$ 220 V, $n =$ 2 and other parameters of the system are shown in Table 1.

Table 1Parameters of the system

Parameters	$M_{i}$ [kg]	$m_{i}$ [kg]	$j_{i}$ [kg·m²]	$r_{i}$ [m]	$f_{i}$ [N·s/m]	$k_{i}$ [N/m]	$f_{r i}$ [N·m·s/rad]
Rotor 1	300	3.5	0.3	0.15	200	7.5×10⁵	3×10^-2
Rotor 2	200	2.5	0.3	0.1	200	7.4×10⁵	1.47×10^-1

We take $T_{v m a x} (ω_{p}^{2}) = δ m_{1} r_{1} m_{2} r_{2} ω_{n}^{4} ω_{p}^{2} / N$ into consideration. $T_{v} = \sin ∆ α T_{v m a x} (ω_{p}^{2})$ can be obtained easily. Thus, $T_{v m a x} (ω_{p}^{2})$ is the maximum vibration torque of the system (MVT). For the after-resonance system, denominator of $N$ may be zero when $ω_{p}$ come to be a specific value $χ$ . $χ^{2}$ is deduced as:

6

χ = \sqrt{\frac{(ω_{1}^{2} - ω_{n}^{2}) (ω_{2}^{2} - ω_{n}^{2})}{ω_{n}^{2} - ω_{1}^{2} + η (ω_{n}^{2} {- ω}_{2}^{2})}} .

$χ$ is called the characteristic frequency (CF) of the system, as shown in Fig. 2(a) and Fig. 3. In the coordinates of $ω_{n}$ and $ω_{p}^{2}$ , the characteristic curve composed of characteristic frequencies at different synchronous speeds is shown in Fig. 2(b). The curve of $T_{v m a x} (ω_{p}^{2})$ accompanied with $ω_{p}^{2}$ is shown in Fig. 3 when $ω_{n}$ is a certain value (for example, $ω_{n} =$ 155 rad/s). The synchronous speed varies with the stiffness of the coupling spring and it can be obtained by numerical simulation.

When $ω_{p}^{2}$ approach $χ^{2}$ from the left side, $T_{v m a x}$ tends to infinity; when $ω_{p}^{2}$ continues to increase over $χ^{2}$ to infinity, $T_{v m a x}$ gradually decreases and tends to a constant value. In Fig. 3, the curve $T_{v m a x}$ is divided into four parts by the curves of $\pm (∆ T_{e} (ω_{n}) - ∆ T_{f} (ω_{n}))$ . The four parts are denoted by LA (passing through point A), LB (passing through point B), LC (passing through point C) and LD (passing through point D) respectively. According to Eq. (4), self-synchronization of two rotors cannot be obtained when system state occurs on LA and LD; On the contrary, rotations of two rotors can be self-synchronizing when system state occurs on LB and LC.

As coupling frequency $ω_{p}$ is close to the characteristic frequency $χ$ , or the system state is near the characteristic frequency curve, the system coupling performance is strong, and the self synchronization of two rotors can be obtained easily. And it is convenient to control the synchronization performance by adjusting the coupling spring stiffness $k_{p}$ .

Fig. 2Surface of the maximum vibration momentvarying with the coupling frequency and synchronous speed in the after-resonance system

a) The three-dimensional diagram

b) The two-dimensional diagram

Fig. 3Relationship between the maximum vibration moment and coupling frequency

4. Simulations for synchronization of two rotors

Simulations are carried out with $ω_{p}^{2}$ set to be 6800 (rad/s)², 11600 (rad/s)², 16000 (rad/s)², 35600 (rad/s)² and infinity respectively, corresponding to points A, B, C and D.

Fig. 4 shows the simulation results with $ω_{p}^{2} =$ 6800 (rad/s)². From Figs. 2 and 3, we can suggest that rotations of two rotors would not be self-synchronizing in this case. As shown in Fig. 4, speeds of the two rotors are not consistent. Fig. 5 shows the simulation results with $ω_{p}^{2} =$ 11600 (rad/s)². In this case, $|T_{v m a x}| > |{∆ T}_{e} (ω_{n}) - {∆ T}_{f} (ω_{n})|$ . $∆ α$ is calculated to be 3.28 rad. As shown in Fig. 5, speeds of the two rotors reach the same value around 1.5 s, and $∆ α$ is stable at 15.83 rad (15.83 – 4 × $π$ = 3.26 (rad)). Fig. 6 shows the simulation results with $ω_{p}^{2} =$ 16000 (rad/s)². Parameters of system satisfy the self-synchronization conditions in this case too. $∆ α$ is calculated to be 3.65 rad. As shown in Fig. 6, speeds of the two rotors reach the same value around 1.5 s, and $∆ α$ is stable at 16.21 rad (16.21 – 4 × $π$ = 3.64 (rad)). The numerical results in Figs. 5 and 6 are consistent with the theoretical analysis.

Fig. 4Simulation results of the after-resonance system when the coupling stiffness is 1.36×106 N/m