Free vibration of the electromechanical integrated magnetic gear system

Hao, Xiuhong; Zhu, Xuejun; Zhang, Hongfei

Journal of Vibroengineering

Browse Journal

Submit article

Published: 15 May 2015

Check for updates

Free vibration of the electromechanical integrated magnetic gear system

Xiuhong Hao¹

Xuejun Zhu²

Hongfei Zhang³

^{1, 3}School of Mechanical Engineering, Yanshan University, Qinhuangdao, China

²Key Laboratory of Advanced Forging & Stamping Technology and Science, Yanshan University, Ministry of Education of China, Qinhuangdao, China

²Hebei Provincial Key Laboratory of Parallel Robot and Mechatronic System, Yanshan University, Qinhuangdao, China

Corresponding Author:

Xiuhong Hao

Cite the article Download PDF

Downloads 1103

Abstract

The electromechanical integrated magnetic gear (EIMG), in which the field modulated magnetic gear, drive and control are integrated, is proposed in this paper. The dynamic model of the EIMG system with four subsystems is founded and the model assumptions are given. Then, the electromagnetic coupling stiffnesses are calculated by the finite element method and the dynamic differential equations are deduced. On the basis of the modal analyses of the EIMG system, the changes of the natural frequencies with the system parameters are discussed. The results show that the electromagnetic coupling sitffnesses change periodically with the relative rotation angles. The EIMG system has five torsional modes and five transverse modes, which have entirely different modal characteristics. The natural frequencies of the EIMG system are affected greatly by the system parameters.

1. Introduction

Magnetic gear is a kind of the magnetic transmission and offers multiple advantages, such as non-contact, non-wear, low vibration and noises, and so on. Magnetic gears overcome the mechanical fatigues and other disadvantages of the mechanical gears, and have significant advantages, such as little maintenance, improved reliability, no lubrication, inherent overload protection [1]. Field modulated magnetic gear (FMMG) proposed by K.Atallah adopts the coaxial topology and has a high utilization of the permanent magnets (PMs) [2]. So, FMMG can provide larger torque and higher torque density than the traditional magnetic gears, which adopt the parallel shaft topology. FMMG can be widely used in the medicine, chemical industry, vehicle, navigation and other fields [3].

FMMG has attracted the attentions of many scholars and extensive researches are carried out. Transmission mechanism [5], torque characteristics [6], structural optimization [7], transmission efficiency [8], rotor eccentricity [9] and dynamics [10] have been discussed. A variety of new transmissions [11, 12] have been proposed. FMMGs have been rapidly applied to many fields. However, FMMGs must be driven by the excellent-performance motors. The whole drive system will take up a bigger space and the transmission performance will be affected by the motor performance.

In this paper, a new kind of an electromechanical integrated magnetic gear is presented by authors, in which FMMG, drive and control are integrated. It is named the electromechanical integrated magnetic gear (EIMG). Except for the advantages of FMMG, EIMG has a compact structure and makes the electromechanical system simplify drastically. It can provide larger torque at a lower speed and can be applied in robot control, aerospace, navigation, vehicle and other drive fields, in which high control precisions are required.

In order to optimize the design parameters and improve dynamic characteristics, the dynamics of the EIMG system must be studied. Meanwhile, the dynamic model and the free vibration are the bases of the forced vibration, parametric vibration and nonlinear vibration. In this paper, an electromagnetic coupling dynamic model of the EIMG system with four subsystems is presented. The modal characteristics and the influences of the main design parameters on the natural frequencies are discussed.

2. Dynamic model and the dynamic differential equations of the EIMG system

2.1. Model assumptions

EIMG shown in Fig.1 is composed of the inner stator, the inner and outer ferromagnetic pole-pieces (FPs), the inner rotor, the outer stator and four air-gaps. Three-phase coils are arranged in the inner stator and generate a rotating magnetic field in order to provide power. PMs are arranged uniformly on the inner surface of the outer stator, the inner and outer surfaces of the inner rotor. The inner and outer FPs are composed of the permeability and non-permeability materials at regular intervals. The FPs takes charge of modulating the magnetic fields in two air-gaps beside them in order to make the number of pole pairs of the PMs on the inner and outer rotors agree with the number of pole pairs of the space harmonic flux density of the air-gaps.

Fig. 1Topology and prototype of the electromechanical integrated magnetic gear

When the outer FPs is fixed, the dynamic model of the EIMG system shown in Fig. 2 contains four subsystems, namely, the inner stator/inner FPs subsystem, the inner FPs/inner rotor subsystem, the inner rotor/outer FPs subsystem and the outer FPs/outer stator subsystem. The electromagnetic coupling dynamic model employs following assumptions [13]:

(1) Main components of the EIMG system are considered to be rigid, assuming that the elastic deformations of PMs, FPs, the rotor and the stator are negligible.

(2) Magnetic interactions among all parts are modeled as the linear spring along tangential direction and normal to their axes. Supports of the inner rotor and the outer FPs can be equivalent to the transverse (i.e. radial) linear spring normal to their axes. Constraints between the inner stator, the inner FPs, the outer stator and the foundation are equivalent to the tangential linear spring and the transverse linear spring, respectively.

(3) Frictional forces among parts are considered to be negligible. Time-varying components of the electromagnetic coupling stiffnesses due to the magnetic field modulating are negligible.

(4) Although the EIMG system has the maximum output torque, out-of-step because of the overload doesn’t occur. Manufacturing and installation errors of all components are not included in this paper.

(5) The effects of the inertias and load waves connected to the EIMG system at the input and output sides are excluded here, considering that these can only cause resonances at very low frequencies.

(6) Amplitudes and frequencies of the electric currents in the inner stator are constants always.

(7) All PMs on the inner and outer surfaces of the inner rotor are assumed to be identical with the same size, performance parameters, respectively. All PMs on the outer stator have the same size, performance parameters too. Magnetic components and non-magnetic components in the inner and outer FPs have the same sizes and performance parameters, respectively.

Fig. 2Dynamic model of the electromechanical integrated magnetic gear

a) The inner stator/inner FPs subsystem

b) The inner FPs/inner rotor subsystem

c) The inner rotor/outer FPs subsystem

d) The outer FPs/outer rotor subsystem

The dynamic model of the EIMG system allows each part translate in $y$ directions and rotate about their translation axes. Torsional displacements of the inner stator, the inner FPs, the inner rotor, the outer FPs and the outer stator, are $θ_{1}$ , $θ_{2}$ , $θ_{3}$ , $θ_{4}$ , $θ_{5}$ , respectively. For convenience, the torsional angular displacements are all replaced by their corresponding translational displacements as:

u_{i} = R_{i} θ_{i}, i = 1, 2, 3, 4, 5,

where $R_{i}$ , are the equivalent radius of gyration of the inner stator, the inner FPs, the inner rotor, the outer FPs and the outer stator, respectively.

Transverse vibration displacements of the inner stator, the inner FPs, the inner rotor, the outer FPs and the outer stator are $y_{1}$ , $y_{2}$ , $y_{3}$ , $y_{4}$ , $y_{5}$ , respectively. Then, the generalized displacement vector of the EIMG system can be written as:

x = {[\begin{array}{l} u_{1} & y_{1} & u_{2} & y_{2} & u_{3} & y_{3} & u_{4} & y_{4} & u_{5} & y_{5} \end{array}]}^{T} .

2.2. Inner stator/inner FPs subsystem

Fig. 2(a) illustrates the dynamic model of the inner stator/inner FPs subsystem. The undamped differential equations of 4 degrees of freedom (DOF) model of the inner stator/inner FPs subsystem are given as follows:

1

\{\begin{array}{l} M_{1} {\ddot{u}}_{1} + k_{I s} x_{12} c o s α_{12} + k_{1} u_{1} = 0, \\ m_{1} {\ddot{y}}_{1} + k_{I s} x_{12} s i n α_{12} + k_{y 1} y_{1} = 0, \\ M_{2} {\ddot{u}}_{2} - k_{I s} x_{12} c o s α_{12} + k_{2} u_{2} = 0, \\ m_{2} {\ddot{y}}_{2} - k_{I s} x_{12} s i n α_{12} + k_{y 2} y_{2} = 0, \end{array}

where $m_{1}$ and $m_{2}$ are the masses of the inner stator and the inner FPs, respectively; $M_{1}$ and $M_{2}$ are the equivalent masses of the inner stator and the inner FPs along their torsional vibration direction, respectively, $M_{1} = J_{1} / R_{1}^{2}$ , $M_{2} = J_{2} / R_{2}^{2}$ ; $J_{1}$ and $J_{2}$ the mass moments of inertia of the inner stator and the inner FPs, respectively, $J_{1} = m_{1} \cdot R_{1}^{2} / 2$ , $J_{2} = m_{2} \cdot R_{2}^{2} / 2$ ; $k_{y 1}$ and $k_{y 2}$ are the transverse supporting stiffnesses of the inner stator and the inner FPs, respectively; $k_{1}$ and $k_{2}$ are the torsional supporting stiffnesses of the inner stator and the inner FPs around their axes, respectively; $k_{I s}$ is the electromagnetic coupling stiffness between the inner stator and the inner FPs, $k_{I s} = \sqrt{k_{I s r}^{2} + k_{I s t}^{2}}$ ; $k_{I s r}$ and $k_{I s t}$ are the radial and tangential components of $k_{I s}$ , respectively; $x_{12}$ and $α_{12}$ are the relative displacement and the meshing angle between the inner stator and the inner FPs, respectively, $α_{12} = a t a n (k_{I s r} / k_{I s t})$ .

The relative displacement between the inner stator and the inner FPs $x_{12}$ can be calculated as:

2

x_{12} = (u_{1} - u_{2}) c o s α_{12} + (y_{1} - y_{2}) s i n α_{12} .

By substituting Eq. (2) into Eq. (1), the undamped differential equations of the inner stator/inner FPs subsystem are gotten as follows:

3

\{\begin{array}{l} M_{1} {\ddot{u}}_{1} + k_{I s} (u_{1} - u_{2}) c o s^{2} α_{12} + k_{12} (y_{1} - y_{2}) s i n α_{12} c o s α_{12} + k_{1} u_{1} = 0, \\ m_{1} {\ddot{y}}_{1} + k_{I s} (u_{1} - u_{2}) s i n α_{12} c o s α_{12} + k_{I s} (y_{1} - y_{2}) s i n^{2} α_{12} + k_{y 1} y_{1} = 0, \\ M_{2} {\ddot{u}}_{2} - k_{I s} (u_{1} - u_{2}) c o s^{2} α_{12} - k_{12} (y_{1} - y_{2}) s i n α_{12} c o s α_{12} + k_{2} u_{2} = 0, \\ m_{2} {\ddot{y}}_{2} - k_{I s} (u_{1} - u_{2}) s i n α_{12} c o s α_{12} - k_{I s} (y_{1} - y_{2}) s i n^{2} α_{12} + k_{y 2} y_{2} = 0 . \end{array}

2.3. Inner FPs/inner rotor subsystem

Fig. 2(b) shows the 4 DOF dynamic model of the inner FPs/inner rotor subsystem. The undamped differential equations of the subsystem can be expressed as follows:

4

\{\begin{array}{l} M_{2} {\ddot{u}}_{2} - k_{I I} x_{23} c o s α_{23} + k_{2} u_{2} = 0, \\ m_{2} {\ddot{y}}_{2} - k_{I I} x_{23} s i n α_{23} + k_{y 2} y_{2} = 0, \\ M_{3} {\ddot{u}}_{3} + k_{I I} x_{23} c o s α_{23} = 0, \\ m_{3} {\ddot{y}}_{3} + k_{I I} x_{23} s i n α_{23} + k_{y 3} y_{3} = 0, \end{array}

where $m_{3}$ and $M_{3}$ are the mass and the equivalent mass of the inner rotor, respectively, $M_{3} = J_{3} / R_{3}^{2}$ ; $J_{3}$ is the mass moment of inertia of the inner rotor, $J_{3} = m_{3} \cdot R_{3}^{2} / 2$ ; $k_{I I}$ is the electromagnetic coupling stiffness between the inner FPs and the inner rotor, $k_{I I} = \sqrt{k_{I I r}^{2} + k_{I I t}^{2}}$ ; $k_{I I r}$ and $k_{I I t}$ are the radial and tangential components of $k_{I I}$ , respectively; $k_{y 3}$ is the transverse supporting stiffness of the inner rotor; $x_{23}$ and $α_{23}$ are the relative displacement and the meshing angle between the inner FPs and the inner rotor, respectively, $α_{23} = a t a n (k_{I I r} / k_{I I t})$ .

The relative displacement between the inner rotor and the inner FPs $x_{23}$ can be calculated as:

5

x_{23} = (u_{3} - u_{2}) c o s α_{23} + (y_{3} - y_{2}) s i n α_{23} .

By substituting Eq. (5) into Eq. (4), the dynamic differential equations of the inner FPs/inner rotor subsystem are derived:

6

\{\begin{array}{l} M_{2} {\ddot{u}}_{2} - k_{I I} (u_{3} - u_{2}) c o s^{2} α_{23} - k_{I I} (y_{3} - y_{2}) s i n α_{23} c o s α_{23} + k_{2} u_{2} = 0, \\ m_{2} {\ddot{y}}_{2} - k_{I I} (u_{3} - u_{2}) s i n α_{23} c o s α_{23} - k_{I I} (y_{3} - y_{2}) s i n^{2} α_{23} + k_{y 2} y_{2} = 0, \\ M_{3} {\ddot{u}}_{3} + k_{I I} (u_{3} - u_{2}) c o s^{2} α_{23} + k_{I I} (y_{3} - y_{2}) s i n α_{23} c o s α_{23} = 0, \\ m_{3} {\ddot{y}}_{3} + k_{I I} (u_{3} - u_{2}) s i n α_{23} c o s α_{23} + k_{I I} (y_{3} - y_{2}) s i n^{2} α_{23} + k_{y 3} y_{3} = 0 . \end{array}

2.4. Inner rotor/outer FPs subsystem

The 4 DOF dynamic model shown in Fig.2(c) is used to describe the coupling between the inner rotor and the outer FPs. The undamped differential equations of the subsystem can be got:

7

\{\begin{array}{l} M_{3} {\ddot{u}}_{3} + k_{I o} x_{34} c o s α_{34} = 0, \\ m_{3} {\ddot{y}}_{3} + k_{I o} x_{34} s i n α_{34} + k_{y 3} y_{3} = 0, \\ M_{4} {\ddot{u}}_{4} - k_{I o} x_{34} c o s α_{34} = T_{o f} / R_{4}, \\ m_{4} {\ddot{y}}_{4} - k_{I o} x_{34} s i n α_{34} + k_{y 4} y_{4} = 0, \end{array}

where $m_{4}$ and $M_{4}$ are the mass and the equivalent mass of the outer FPs, respectively, $M_{4} = J_{4} / R_{4}^{2}$ ; $J_{4}$ is the mass moment of inertia of the outer FPs, $J_{4} = m_{4} \cdot R_{4}^{2} / 2$ ; $k_{I o}$ is the electromagnetic coupling stiffness between the outer FPs and the inner rotor, $k_{I o} = \sqrt{k_{I o r}^{2} + k_{I o t}^{2}}$ ; $k_{I o r}$ and $k_{I o t}$ are the transverse and tangential components of $k_{I o}$ , respectively; $k_{y 4}$ is the transverse supporting stiffness of the outer FPs; $x_{34}$ and $α_{34}$ are the relative displacement and the meshing angle between the outer FPs and the inner rotor, respectively, $α_{34} = a t a n (k_{I o r} / k_{I o t})$ ; $T_{o}$ and $T_{o f}$ are the torque on the inner rotor and the outer FPs, respectively.

The relative displacement between the inner rotor and the outer FPs $x_{34}$ can be given as:

8

x_{34} = (u_{3} - u_{4}) c o s α_{34} + (y_{3} - y_{4}) s i n α_{34} .

By substituting Eq. (8) into Eq. (7), the dynamic differential equations of the inner rotor/outer FPs subsystem can be obtained:

9

\{\begin{array}{l} M_{3} {\ddot{u}}_{3} + k_{I o} (u_{3} - u_{4}) c o s^{2} α_{34} + k_{I o} (y_{3} - y_{4}) s i n α_{34} c o s α_{34} = 0, \\ m_{3} {\ddot{y}}_{3} + k_{I o} (u_{3} - u_{4}) s i n α_{34} c o s α_{34} + k_{I o} (y_{3} - y_{4}) s i n^{2} α_{34} + k_{y 3} y_{3} = 0, \\ M_{4} {\ddot{u}}_{4} - k_{I o} (u_{3} - u_{4}) c o s^{2} α_{34} - k_{I o} (y_{3} - y_{4}) s i n α_{34} c o s α_{34} = T_{o f} / R_{4}, \\ m_{4} {\ddot{y}}_{4} - k_{I o} (u_{3} - u_{4}) s i n α_{34} c o s α_{34} - k_{I o} (y_{3} - y_{4}) s i n^{2} α_{34} + k_{y 4} y_{3} = 0 . \end{array}

2.5. Outer FPs/outer stator subsystem

The 4 DOF dynamic model shown in Fig. 2(d) illustrates the subsystem of the outer FPs and the outer stator. The undamped differential equations of the subsystem can be expressed:

10

\{\begin{array}{l} M_{4} {\ddot{u}}_{4} - k_{o o} x_{45} c o s α_{45} = T_{o f} / R_{4}, \\ m_{4} {\ddot{y}}_{4} - k_{o o} x_{45} s i n α_{45} + k_{y 5} y_{5} = 0, \\ M_{5} {\ddot{u}}_{5} + k_{o o} x_{45} c o s α_{45} + k_{5} u_{5} = 0, \\ m_{5} {\ddot{y}}_{5} + k_{o o} x_{45} s i n α_{45} + k_{y 5} y_{5} = 0, \end{array}

where $m_{5}$ and $M_{5}$ are the mass and the equivalent mass of the outer stator, respectively, $M_{5} = J_{5} / R_{5}^{2}$ ; $J_{5}$ is the mass moment of inertia of the outer stator, $J_{5} = m_{5} \cdot R_{5}^{2} / 2$ ; $k_{o o}$ is the electromagnetic coupling stiffness between the outer FPs and the outer stator, $k_{o o} = \sqrt{k_{o o r}^{2} + k_{o o t}^{2}}$ ; $k_{o o r}$ and $k_{o o t}$ are the radial and tangential components of $k_{o o}$ , respectively; $k_{y 5}$ is the transverse supporting stiffness of the outer stator; $x_{45}$ and $α_{45}$ are the relative displacement and the meshing angle between the outer FPs and the outer stator, respectively, $α_{45} = a t a n (k_{o o r} / k_{o o t})$ ; $k_{5}$ are the torsional supporting stiffnesses of the outer stator around its axis.

The relative displacement between the outer FPs and the outer stator $x_{45}$ can be calculated as follows:

11

x_{45} = (u_{5} - u_{4}) c o s α_{45} + (y_{5} - y_{4}) s i n α_{45} .

By substituting Eq. (11) into Eq. (10), the dynamic differential equations of the outer FPs/outer stator subsystem are:

12

\{\begin{array}{l} M_{4} {\ddot{u}}_{4} - k_{o o} (u_{5} - u_{4}) c o s^{2} α_{45} - k_{o o} (y_{5} - y_{4}) s i n α_{45} c o s α_{45} = 0, \\ m_{4} {\ddot{y}}_{4} - k_{o o} (u_{5} - u_{4}) s i n α_{45} c o s α_{45} - k_{o o} (y_{5} - y_{4}) s i n^{2} α_{45} + k_{y 4} y_{4} = 0, \\ M_{5} {\ddot{u}}_{5} + k_{o o} (u_{5} - u_{4}) c o s^{2} α_{45} + k_{o o} (y_{5} - y_{4}) s i n α_{45} c o s α_{45} + k_{5} u_{5} = 0, \\ m_{5} {\ddot{y}}_{5} + k_{o o} (u_{5} - u_{4}) s i n α_{45} c o s α_{45} + k_{o o} (y_{5} - y_{4}) s i n^{2} α_{45} + k_{y 5} y_{5} = 0 . \end{array}

2.6. The overall system

Eq. (3), Eq. (6), Eq. (9) and Eq. (12), which define the differential equations of individual subsystem, are combined systematically to obtain the differential equations of the overall EIMG system. The undamped differential equations of the 10 DOF of the overall EIMG system can be given in matrix form as:

13

m \ddot{x} + k x = F .

The mass matrix $m$ and the load vector $F$ are given respectively as follows:

m = d i a g ([\begin{array}{l} M_{1} & m_{1} & M_{2} & m_{2} & M_{3} & m_{3} & M_{4} & m_{4} & M_{5} & m_{5} \end{array}]),

F = {[\begin{array}{l} 0 & 0 & 0 & 0 & 0 & 0 & T_{o f} / R_{4} & 0 & 0 & 0 \end{array}]}^{T} .

Because $k$ is a 10×10 matrix and bigger, the elements in the matrix $k$ can be expressed as follows, respectively:

k_{11} = k_{I I} c o s^{2} α_{12} + k_{1}, k_{12} = k_{I s} s i n α_{12} c o s α_{12}, k_{13} = - k_{I s} c o s^{2} α_{12},

k_{14} = - k_{I s} s i n α_{12} c o s α_{12}, k_{15} = \dots = k_{1 i} = \dots = k_{110} = 0, k_{21} = k_{12},

k_{22} = k_{I s} s i n^{2} α_{12} + k_{y 1}, k_{23} = - k_{I s} s i n α_{12} c o s α_{12}, k_{24} = - k_{I s} s i n^{2} α_{12},

k_{25} = \dots = k_{2 i} = \dots = k_{210} = 0, k_{31} = k_{13}, k_{32} = k_{23},

k_{33} = k_{I s} c o s^{2} α_{12} + k_{I I} c o s^{2} α_{23} + k_{2}, k_{34} = k_{I s} s i n α_{12} c o s α_{12} + k_{I I} s i n α_{23} c o s α_{23},

k_{35} = - k_{I I} c o s^{2} α_{23}, k_{36} = - k_{I I} s i n α_{23} c o s α_{23}, k_{37} = k_{38} = k_{39} = k_{310},

k_{41} = k_{14}, k_{42} = k_{14}, k_{43} = k_{34}, k_{44} = k_{I s} s i n^{2} α_{12} + k_{I I} s i n^{2} α_{23} + k_{y 2},

k_{45} = - k_{I I} s i n α_{23} c o s α_{23}, k_{46} = - k_{I I} s i n^{2} α_{23}, k_{47} = k_{48} = k_{49} = k_{410} = 0,

k_{51} = \dots = k_{61} = \dots = k_{101} = 0, k_{52} = \dots = k_{62} = \dots = k_{102} = 0, k_{53} = k_{35},

k_{54} = k_{45}, k_{55} = k_{I I} c o s^{2} α_{23} + k_{I o} c o s^{2} α_{34},

k_{56} = k_{I I} s i n α_{23} c o s α_{23} + k_{I o} s i n α_{34} c o s α_{34}, k_{57} = - k_{I o} c o s^{2} α_{34},

k_{58} = - k_{I o} s i n α_{34} c o s α_{34}, k_{59} = k_{510} = 0, k_{63} = k_{36}, k_{64} = k_{46}, k_{65} = k_{56},

k_{66} = k_{I I} s i n^{2} α_{23} + k_{I o} s i n^{2} α_{34} + k_{y 3}, k_{67} = - k_{I o} s i n α_{34} c o s α_{34}, k_{68} = - k_{I o} s i n^{2} α_{34},

k_{69} = k_{610} = 0, k_{73} = k_{74} = 0, k_{75} = k_{57}, k_{76} = k_{67},

k_{77} = k_{I o} s i n α_{34} c o s α_{34} + k_{o o} s i n α_{45} c o s α_{45}, k_{78} = k_{I o} s i n α_{34} c o s α_{34} + k_{o o} s i n α_{45} c o s α_{45},

k_{79} = - k_{o o} c o s^{2} α_{45}, k_{710} = - k_{o o} s i n α_{45} c o s α_{45}, k_{83} = k_{84} = 0, k_{85} = k_{58},

k_{86} = k_{68}, k_{87} = k_{78}, k_{88} = k_{I o} s i n^{2} α_{34} + k_{o o} s i n^{2} α_{45} + k_{y 4},

k_{89} = - k_{o o} s i n α_{45} c o s α_{45}, k_{810} = - k_{o o} s i n^{2} α_{45}, k_{93} = k_{94} = k_{95} = k_{96} = 0,

k_{97} = k_{79}, k_{98} = k_{89}, k_{99} = k_{o o} c o s^{2} α_{45} + k_{5}, k_{910} = k_{o o} s i n α_{45} c o s α_{45},

k_{103} = k_{104} = k_{105} = k_{106} = 0, k_{107} = k_{710}, k_{108} = k_{810}, k_{109} = k_{910},

k_{1010} = k_{o o} s i n^{2} α_{45} + k_{y 5} .

Free vibration differential equations of the EIMG system are written in matrix form as:

14

m \ddot{x} + k x = 0,

where $0$ is a $10 \times 1$ zero vector.

3. Electromagnetic coupling stiffness

All components in the EIMG system are coupled by the electromagnetic field. Here, magnetic circuits are considered to be linear and unsaturated. The two-degree finite element model of the example EIMG system shown in Table 1 can be founded in Ansys. The initial meshes and the magnetic flux density can be calculated and is shown in Fig. 3.

The starting points of the electromagnetic coupling FEM computations of the EIMG system are four basic laws of Maxwell electromagnetic field, namely Ampere circuit law, Gauss law, Faraday’s law of induction and law of continuity of magnetic flux. The whole computational domain is divided into thousands of grid cells. By the variational principle, the partial differential equations are converted to multiple algebraic equations and can be solved by numerical method. The reluctances of the air-gaps are bigger than permanent magnets and the back iron. Meanwhile, the thicknesses of the air-gaps are smaller. So, the grids of the air-gap domains are small and dense.

The meshes on the back iron of the inner rotor and the meshes on the permanent magnets of the inner rotor have common boundary nodes. During running of the EIMG system, all unit nodes on the permanent magnets and the back iron of the inner rotor rotate together as a rigid body. Similarly, the meshes on the back iron of the outer stator and the meshes on the permanent magnets of the outer stator have common boundary nodes. The meshes on the magnetic components and non-magnetic components of the inner and outer FPs have common boundary nodes. The meshes on the magnetic components and non-magnetic components of the outer FPs rotate together as a rigid body. The middle positions of the air-gaps aren’t the common areas of the nodes. During running of the EIMG system, the magnetic field distribution is calculated by interpolations among the boundary nodes. Although there are some errors, the errors are negligible if the meshes of the air-gaps are very small.

Fig. 3The finite element model of the EIMG system

a) Initial meshes of the EIMG system

b) Magnetic field distribution

c) Arrangements of PMs segments

The magnetic flux densities of the middle positions of the air-gaps $B_{j}$ , $j$ = 1, 2, 3, 4, 5, can be obtained. $B_{r j}$ and $B_{θ j}$ are the radial and tangential components of the magnetic flux densities $B_{j}$ , respectively. The average torque on each component can be achieved by integrating the Maxwell stress tensor in Ansys. Meanwhile, the electromagnetic coupling forces on all parts can be calculated as follows:

15

F_{r i} = \frac{{R'}_{j} L}{μ_{0}} \int_{0}^{2 π} \frac{B_{j} r^{2}}{2 \times μ} d θ,

16

F_{t i} = \frac{{R'}_{j} L}{μ_{0}} \int_{0}^{2 π} B_{θ j} B_{r j} d θ,

where $F_{r i}$ and $F_{t i}$ , $i$ = 1, 2, 3, 4, 5, are the radial and tangential components of the electromagnetic coupling forces, respectively; $L$ is the axial length of the EIMG system; ${R'}_{j}$ , $j$ = 1, 2, 3, 4, is the radius of the middle position of the each air-gap; $μ_{0}$ is the space permeability.

The electromagnetic coupling stiffnesses among components can be obtained in Ansys:

17

k_{i r} = d F_{i r} / d r = \frac{R_{j} L}{μ_{0}} \int_{0}^{2 π} \frac{B_{j} r^{2}}{2 \times μ} d θ / d r,

18

k_{i t} = d F_{i t} / d u_{i t} = d F_{i t} / (R_{j} d θ) = (\frac{L}{μ_{0}} \int_{0}^{2 π} B_{θ j} B_{r j} d θ) / d θ,

where $k_{i r}$ and $k_{i t}$ , $i$ = 1, 2, 3, 4, 5, are the radial and tangential components of the electromagnetic coupling forces among components.

When the three-phase alternating currents in the inner stator are switched on, the electromagnetic stiffnesses among components will vary with the relative position angles among components in Ansys and can be worked out in Fig. 4.

Fig. 4Electromagnetic coupling stiffnesses among parts of the EIMG system

Fig. 4 indicates that the tangential electromagnetic coupling stiffnesses, $k_{I s t}$ , $k_{I I t}$ , $k_{I o t}$ and $k_{o o t}$ , vary sinusoidally with the relative position angles among the components. When there isn’t relative rotation between the inner stator and the inner rotor, or between the inner rotor and the outer stator, the EIMG system will stabilize at zero point, also called no-load point or balance point. Now, there isn’t output torque on the outer FPs, namely, the tangential magnetic coupling stiffnesses are zero. The tangential electromagnetic coupling stiffnesses increase with the relative rotation angle among parts increasing. When the relative rotation angle comes to a certain angle, $θ = {2 π / p}_{i}$ , the tangential electromagnetic coupling stiffnesses will achieve the maxima. Where, $p_{i}$ is the number of pole pairs of PMs on the inner rotor or the number of pole pairs of the three-phase alternating currents on the inner stator.

Fig. 4 shows that the radial electromagnetic stiffnesses $k_{I s r}$ , $k_{I I r}$ , $k_{I o r}$ and $k_{o o r}$ vary sinusoidally with the relative position angles. When the relative rotation angle comes to a certain angle, $θ = {2 π / p}_{i}$ , the radial electromagnetic coupling stiffnesses will achieve the maxima. But theaverage stiffnesses aren’t zero.

Meanwhile, Fig. 4 indicates that the radial electromagnetic stiffnesses are much smaller than the tangential electromagnetic stiffnesses. This is because magnetic fields in Fig. 3(b) and the radial electromagnetic forces on all parts are nearly evenly distributed along the entire circumference. So, the radial resultant forces are less than tangential resultant forces. Accordingly, the radial electromagnetic stiffnesses are smaller.

Table 1Parameters of example EIMG system

Number of pole pairs on the inner stator	2	Number of pole pairs on the inner surface of the inner rotor	5
Number of pole pairs on the outer surface of the inner rotor	4	Number of pole pairs on the outer stator	17
Number of the inner ferromagnetic pole pieces	7	Number of the outer ferromagnetic pole pieces	21
Outer radius of the inner stator / mm	98	Inner radius of the inner FPs / mm	100
Outer radius of the inner FPs / mm	130	Thickness of PMs on the inner rotor / mm	10
Inner radius of the inner rotor yoke / mm	152	Outer radius of the inner rotor yoke / mm	182
Inner radius of the outer FPs / mm	204	Outer radius of the outer FPs / mm	234
Thickness of PMs on the outer stator / mm	10	Inner radius of the outer stator yoke / mm	256
Outer radius of the outer stator yoke / mm	276	Axial length / mm	40
Remanence of PMs / T	1.3	Coercive force of PMs / KOe	11.6
Magnitude of currents on the inner stator / A	20

4. Free vibration of the EIMG system

When the output torque on the outer FPs is the maximum, the electromagnetic coupling stiffnesses among all components will be close to maxima. Now, the characteristic parameters of the EIMG system are shown in Table 2. Substituting the design parameters into the Eq. (14), the natural frequencies and the modal shapes can be worked out in Matlab and shown in Table 3.

Table 2Characteristic parameters of the example EIMG system

$k_{oo}$ (N/m)	$k_{Io}$ (N/m)	$k_{II}$ (N/m)	$k_{Is}$ (N/m)	$α_{45}$ (rad)	$α_{34}$ (rad)	$α_{23}$ (rad)
8.6444×10⁵	3.1446×10⁵	4.5467×10⁵	2.0416×10⁵	0.1014	0.1818	0.1435
$α_{12}$ (rad)	$m_{1}$ (kg)	$m_{2}$ (kg)	$m_{3}$ (kg)	$m_{4}$ (kg)	$m_{5}$ (kg)	$k_{y 1}$ (N/m)
0.2022	1.3	1.3	11	6.2	7.7	3×10⁶
$k_{y 2}$ (N/m)	$k_{2}$ (N/m)	$k_{y 3}$ (N/m)	$k_{y 4}$ (N/m)	$k_{y 5}$ (N/m)	$k_{5}$ (N/m)
3×10⁶	3×10⁶	3×10⁶	3×10⁶	3×10⁶	3×10⁶

Table 3 indicates that there are ten different modes, in which the natural frequencies and the modal shapes are completely different. In the mode shapes corresponding to the natural frequencies 1322.1 rad/s, 1639.2 rad/s 277.2 rad/s, 453.0 rad/s and 832.6 rad/s, the relative deflections of the torsional degree of freedoms (DOFs) are much bigger than the relative deflections of the transverse DOFs. Also, the torsional relative deflection of the inner stator, the inner FPs, the inner rotor, the outer FPs and the outer stator reaches maximum, respectively. Considering the shape characteristics, the above five modes are named the torsional modes of the inner stator, the inner FPs, the inner rotor, the outer FPs and the outer stator (ISRM, IFRM, IRRM, OFRM and OSRM). Other five modes, in which the transverse relative deflections of all parts are much bigger than the torsional relative deflections, are named the transverse modes of the inner stator, the inner FPs, the inner rotor, the outer FPs and the outer stator (ISTM, IFTM, IRTM, OFTM and OSTM).

Table 3Mode frequencies and mode shapes of the EIMG system

Mode type	ISRM	ISTM	IFRM	IFTM	IRRM
Natural frequencies (rad/s)	1322.1	1529.5	1639.2	1501.9	277.2
Mode shapes	-1.0000 0.0326 -0.3773 0.0011 0.0194 0.0017 -0.0015 -0.0002 0.0003 0.0000	-0.0015 1.0000 0.1705 0.1796 -0.0074 -0.0006 0.0004 0.0000 -0.0000 -0.0000	0.4003 0.1331 -1.0000 -0.2413 0.0335 0.0027 -0.0015 -0.0002 0.0001 0.0000	-0.1691 0.1442 0.4506 -1.0000 -0.0120 -0.0010 0.0006 0.0001 -0.0001 -0.0000	0.0489 0.0033 0.2841 0.0125 1.0000 -0.0294 0.6704 -0.0032 0.3670 0.0109
Mode type	IRTM	OFRM	OFTM	ORRM	ORTM
Natural frequencies (rad/s)	530.6	453.0	694.2	832.6	623.5
Mode shapes	-0.0186 -0.0012 -0.0956 -0.0039 -0.1608 -1.0000 0.1569 -0.0210 0.1719 -0.0018	0.0313 0.0020 0.1688 0.0071 0.5793 -0.1422 -1.0000 0.0598 -0.7975 -0.0120	0.0010 0.0001 0.0046 0.0002 0.0129 0.0027 0.0995 -1.0000 -0.2081 -0.0254	0.0097 0.0005 0.0370 0.0012 0.0925 0.0140 -0.9268 -0.1766 1.0000 0.0726	-0.0009 -0.0001 -0.0041 -0.0002 -0.0119 -0.0037 0.0602 -0.0155 -0.1031 1.0000

Meanwhile, Table 3 shows that the relative influences among the adjacent components are bigger and the relative vibrations among the intervals components are smaller. Namely, when the resonance of the EIMG system occurs and the displacement of a certain DOF reaches maximum, such as the torsional displacement of the inner rotor, the torsional and transverse displacements of the adjacent components, such as the inner and outer FPs, will be bigger and the displacements of the inner and outer stators will be smaller. This is because all components in the EIMG system are coupled by the magnetic field forces rather than the direct contacts. The weak magnetic coupling forces make the influence among intervals components be smaller.

5. Influence of the design parameters on the natural frequencies

When the design parameters change, the natural frequencies will change too. Fig. 5 shows the influence of the design parameters on the natural frequencies.

(1) Fig. 5(a) shows that all the natural frequencies increase with the remanence of PMs $B_{r}$ increasing. The natural frequencies of the torsional modes increase obviously, and others are nearly invariable. That is because that the electromagnetic coupling forces and the stiffnesses among components will increase proportionally with $B_{r}$ increasing. But the radial electromagnetic stiffnesses are smaller than the tangential electromagnetic stiffnesses. When the transverse supporting stiffnesses are almost unchanged, the natural frequencies will change as above.

(2) Fig. 5(b) shows that all the natural frequencies decrease when the axial length $L$ increases. Especially, the natural frequencies of the ISTM and IFTM decrease greatly. This is because the electromagnetic torques on all the components and the electromagnetic coupling stiffnesses among the components will increase, when the axial length $L$ increases. Meanwhile, the masses of all the components increase greatly too. When the torsional and transverse supporting stiffnesses are unchanged, all the natural frequencies decrease.

(3) The torsional and transverse supporting stiffnesses will change a lot, when the supporting systems of the components are different. Fig. 5(c) and Fig. 5(d) indicate that the corresponding mode frequency will increase when a certain supporting stiffnesses increases, and other frequencies are similarly constant.

(4) When the structure parameters or materials are different, the masses of the components will also change. Fig. 5(e) and Fig. 5(f) show that only the natural frequencies of the torsional and transverse modes of this component are affected significantly, when the mass of a certain component increases. That is, the mass change of a certain component hardly affects the other natural frequencies of other components.

Fig. 5Changes of the natural frequencies along with the main parameters

a)

b)

c)

d)

e)

f)

In conclusion, the remanence of PMs and the axial length have great influences on the natural frequencies. The supporting systems only affect the vibration frequencies of a certain component. When the space of the EIMG system is certain, the masses of parts change little and has smaller influences on the natural frequencies.

To optimize the dynamic characteristics of the EIMG system and reduce the vibration, the natural frequencies always need be adjusted. On the premise of the enough output capacity, the better options are to change the remanence of PMs and the effective axial length of the EIMG system. When a certain frequency of a component is adjusted, the supporting system can be changed.

6. Conclusion

The electromagnetic coupling stiffnesses in the EIMG system are much smaller than the meshing stiffness of the mechanical gear drive system. The tangential electromagnetic coupling stiffnesses are much bigger than the radial electromagnetic coupling stiffnesses. Meanwhile, the electromagnetic coupling stiffnesses change periodically with the relative rotation angles among components. The EIMG system with four subsystems includes five torsional modes and five transverse modes, which have entirely different modal characteristics. The natural frequencies of the EIMG system are affected greatly by the remanence of PMs and structure size. The component support of a certain part affects the mode frequency of the corresponding component. All these can provide some theory basis for the parameter designs and optimizations of the EIMG system in order to work steadily in a certain working condition.

References

Rasmussen P. O., Andersen T. O., Joergense F. T., Nielsen O. Development of a high performance magnetic gear. IEEE Transactions on Industry Applications, Vol. 41, Issue 3, 2005, p. 764-770.

Search CrossRef
Atallah K., Howe D. A novel high-performance magnetic gear. IEEE Transactions on Magnetics, Vol. 37, Issue 4, 2001, p. 2844-2846.

Search CrossRef
Baghli L., Gouda E., Mezni S., Rezzoug A. Hybrid vehicle with a magnetic planetary gear. EFEEA’10 International Symposium on Environment Friendly Energies in Electrical Applications, Ghardaïa, Algeria, 2010.

Search CrossRef
Frank N. W., Toliyat H. A. Gearing ratios of a magnetic gear for marine applications. Electric Ship Technologies Symposium, Maryland, USA, 2009, p. 477-481.

Search CrossRef
Atallah K., Wang J., Howe D. A high-performance linear magnetic gear. Journal of Applied Physics, Vol. 97, Issue 10, 2005, p. N516-N516-3.

Search CrossRef
Niguchi N., Hirata K., Hayakawa Y. Study on transmission torque characteristics of a surface-permanent-magnet-type magnetic gear. IEEJ Transactions on Industry Applications, Vol. 131, Issue 3, 2011, p. 396-402.

Search CrossRef
Jian L. N., Xu G., Song J., Xue H., et al. Optimum design for improving modulating effect of coaxial magnetic gear using response surface methodology and genetic algorithm. Progress In Electromagnetics Research, Vol. 116, 2011, p. 297-312.

Search CrossRef
Du Shiqin, Jiang Jianzhong, Zhang Yuejin, Gong Yu A magnetic gearing. Transactions of China Electrotechnical Society, Vol. 25, Issue 9, 2010, p. 41-46.

Search CrossRef
Percebon L. A., Ferraz R., Ferreira D. L., Mauricio V. Modelling of a magnetic gear considering rotor eccentricity. IEEE International Electric Machines and Drives Conference, 2011, p. 1237-1241.

Search CrossRef
Hao Xiuhong, Yuan Xiaoming, Zhang Hongfei, Zhang Lijie Combination resonances of parametric vibration system of the field modulated magnetic gear. Journal of Vibroengineering, Vol. 16, Issue 3, 2014, p. 1590-1601.

Search CrossRef
Jian Linni, Chau K. T., Jiang J. Z. A magnetic-geared outer-rotor permanent-magnet brushless machine for wind power generation. IEEE Transactions on Industry Applications, Vol. 45, Issue 3, 2009, p. 954-962.

Search CrossRef
Rasmussen P. O., Mortensen H. H., Matzen T. N., et al. Motor integrated permanent magnet gear with a wide torque-speed range. IEEE Energy Conversion Congress and Exposition, 2009, p. 1510-1518.

Search CrossRef
Ahmet Kahraman Free torsional vibration characteristics of compound planetary gear sets. Mechanism and Machine Theory, Vol. 39, 2004, p. 953-971.

Search CrossRef

About this article

Received

28 September 2014

Accepted

29 April 2015

Published

15 May 2015

SUBJECTS

Mechanical vibrations and applications

Keywords

magnetic gear

electromechanical integrated

field modulated

dynamic model

free vibration

Acknowledgements

This project is supported by Natural Science Foundation of China (51205341), Research Program of Natural Science at Universities of Hebei province (Q2012032), the Joint Fund of Specialized Research Fund for the Doctoral Program of Higher Education and Hebei Provincial Education Office (20121333120004).

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.