# Forced Responses of the Parametric Vibration System for the Electromechanical Integrated Magnetic Gear

Authors: Xiu-hong Hao; Xue-jun Zhu
Journal: Shock and Vibration
 (2015)

---

## Abstract

Considering the magnetic fields modulating in the electromechanical integrated magnetic gear (EIMG), the electromagnetic coupling stiffnesses vary periodically and the expressions are given by the finite element method. The parametric vibration model and the dynamic differential equations are founded. The expressions of forced responses of EIMG system are deduced when the main resonances and the combination resonances occur. And then, the time and frequency responses are figured out. The dynamic characteristics of EIMG system are discussed. The results show that the dominant frequencies in the resonances are always the natural frequency of EIMG system. The relative amplitudes of the components have great difference and the components amplitudes of the main resonances are much bigger than the components amplitudes of the combination resonances. The time-varying meshing stiffness wave between the inner stator and the inner ferromagnetic pole-pieces has little influence on EIMG system.


---

## Body

## 1. Introduction

Since the high performance rare-earth permanent magnet materials came out in the 1980s, magnetic gears attracted much attention and were widely applied in industries because of their advantages, such as no mechanical fatigue, low noises, little maintenance, no lubrication, and inherent overload protection, while, because they adopt the parallel shaft topology, the traditional magnetic gears have small output torques and small torque density. So, the uses of the traditional magnetic gears are limited greatly [1].Field modulated magnetic gear (FMMG) adopts the coaxial topology and significantly improves the utilization of the permanent magnets (PMs) [2]. Except for the advantages of the traditional magnetic gears, FMMG can provide bigger torque and higher torque density [1–3]. So, FMMG can be widely used in the medicine, chemical, vehicle, navigation, and other fields [4, 5]. Based on the FMMG, the axial and linear field modulated magnetic gears were proposed [3, 6]. Many kinds of the direct current and alternating current permanent magnet motors were presented [7, 8]. Transmission mechanism [2, 3], torque characteristics [3, 4, 9], structural optimization [10], rotor eccentricity [11], and dynamics [12] have been discussed widely. The transmission performance of the FMMG system is improved greatly.The coaxial topology makes FMMG readily be integrated with electric machinery. In this paper, a new kind of the electromechanical integrated magnetic gear (EIMG) is proposed by authors, in which FMMG, drive, and control are integrated. Compared with FMMG, EIMG has two stage transmissions. So, EIMG has a compact structure and makes the electromechanical drive system be simplified drastically. Meanwhile, the output torque and speed can be controlled. EIMG can offer a 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.During running of EIMG system, the ferromagnetic pole-pieces take charge of modulating the magnetic fields in two air-gaps beside them and make the electromagnetic coupling stiffnesses among components vary periodically. In mechanical gear systems, the meshing stiffnesses vary with rectangular waves because of the periodical change of the meshing tooth pair number [13]. However, the meshing number of pole pairs is always a constant and the increment components of the electromagnetic coupling stiffnesses vary with a cosine function. Considering the time-varying electromagnetic coupling stiffnesses, EIMG system belongs to a typical parametric vibration system. When there are some outer exciting sources, the main resonances and the combination resonances will occur. The resonances will worsen the dynamic characteristics of EIMG system and should be avoided.In this paper, the finite element model of EIMG system is built and the expressions of the electromagnetic coupling stiffnesses are given. Considering the time-varying electromagnetic coupling stiffnesses, the parametric vibration dynamic model of EIMG system and the dynamic differential equations are founded. By the multiscale method, the approximate analytical solutions are deduced when the main resonances and the combination resonances occur. Then, the resonance responses are discussed and the results can provide the theoretical foundations for the parameter optimizations and the dynamic characteristics improvement of EIMG system.


## 2. Parametric Vibration Model of EIMG System

### 2.1. Electromagnetic Coupling Stiffnesses of EIMG System

EIMG shown in Figure1 is mainly composed of the inner stator, the inner ferromagnetic pole-pieces (FPs), the inner rotor, the outer FP, and the outer stator. Three-phase coils are arranged in the inner stator and can generate a rotating magnetic field in order to provide power. PMs of the N and S poles 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 take charge of modulating the magnetic fields in two air-gaps beside them in order to make the number of pole pairs of the PMs agree with the number of pole pairs of the space harmonic flux density of the air-gaps.Figure 1
Topology and prototype of the electromechanical integrated magnetic gear.The finite element model of the example EIMG system shown in Table1 can be founded. Here, the material of the inner stator is silicon steel sheets. The back iron materials of the inner rotor and the outer stator are Q235 steel. The materials of the magnetic components and nonmagnetic components in FPs are 23TW250 and epikote, respectively. The properties of the epikote basically agree with air. Meshing is carried out manually and the rapid calculations can be realized by controlling the grid refinement level.Table 1
Parameters of example EIMG system.
Number of pole pairs on the inner surface of the inner rotor
5
Number of pole pairs on the inner stator
2
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 FP/mm
100
Outer radius of the inner FP/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 FP/mm
204
Outer radius of the outer FP/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
Current frequency/Hz
50The meshes between the back iron and PMs of the inner rotor have common boundary nodes. During running of EIMG system, all unit nodes on PMs and the back iron of the inner rotor rotate together as a rigid body. Similarly, the back iron and PMs of the outer stator have common boundary nodes. The meshes on the magnetic components and nonmagnetic components of the inner and outer FPs have common boundary nodes, respectively. The middle positions of the air-gaps are not the common areas of the nodes. The magnetic field distribution is calculated by interpolations among the boundary nodes.The reluctances of the air-gaps are bigger than PMs and the back irons. Meanwhile, the thicknesses of the air-gaps are smaller. So, the grid density in the air-gaps is bigger. Although there are some errors, the errors are negligible if the meshes of the air-gaps are very small. The initial grid model of the example EIMG system and the magnetic field distribution can be figured out in Figure2.Figure 2
Finite element model of the example EIMG system.
(a)
Initial grid model
(b)
Magnetic field distributionThe magnetic flux density in the middle of the air-gaps isBj, where j is the number of four air-gaps, j=1,2,3,4. Bθj and Brj are the tangential and radial components of the magnetic flux density Bj, respectively. The tangential electromagnetic coupling forces among components can be calculated by Maxwell stress tensor method. Then, the electromagnetic coupling stiffnesses among components can be obtained as follows:(1)kit=dFitduit=dFitRjdθ=L/μ0∫02πBθjBrjdθdθ,where Fit and kit, i=Is, II, Io, oo, are the tangential components of the electromagnetic coupling forces and the electromagnetic coupling stiffnesses between the inner stator and the inner FP, the inner FP and the inner rotor, the inner rotor and the outer FP, and the outer FP and the outer stator, respectively; L is the axial effective length of EIMG system; Rj is the radius of the middle position of the air-gap among components; μ0 is the space permeability; θ is the mechanical rotation angle.When the coils in the inner stator are connected to the three-phase alternating currents, the number match of the pole pairs of the PMs is realized by the inner and outer FPs modulating the magnetic fields, while, because of the magnetic fields modulating, the electromagnetic coupling stiffnesses among components vary periodically. The dynamic electromagnetic coupling stiffnesses can be worked out by the finite element method and shown in Figure3.Figure 3
Dynamic electromagnetic coupling stiffnesses among components.
(a)
k
Is
(
t
) and kII(t)
(b)
k
Io
(
t
) and koo(t)
(c)
Frequency component ofkIs(t)
(d)
Frequency component ofkIo(t)Figure3 shows that the stiffnesses fluctuations between the inner stator and the inner FP, the inner rotor and the outer FP, are bigger than others and cannot be ignored. The dynamic electromagnetic coupling stiffnesses contain multiple harmonics. The main frequency of the stiffness fluctuation between the inner stator and the inner FP is the product of the current frequency and the number of pole pairs of the three-phase alternating currents on the inner stator. The main frequency of the stiffness fluctuation between the inner rotor and the outer FP is the product of the angular frequency and the number of pole pairs of the PMs on the outer surface of the inner rotor. Considering the amplitudes of other frequency components are much smaller, the time-varying electromagnetic coupling stiffnesses kIst and kIIt can be simplified as follows:(2)kIst=k-Is+2ΔkIscos⁡ωst=k-Is1+εIseωst+cc,kIot=k-Io+2ΔkIocos⁡ωot=k-Io1+εIoeωot+cc,where k-Is and k-Io are the average values of the electromagnetic coupling stiffnesses kIst and kIot, respectively; ΔkIs and ΔkIo are the increments of the stiffnesses kIst and kIot, respectively; ωs and ωo are the fluctuation frequencies of the stiffness increments ΔkIs and ΔkIo, respectively; namely, meshing frequencies ωs=2πnsps/60, ωo=2πnIpI/60; ns and ps are the frequency and the pole pairs of the three-phase alternating currents on the inner stator, respectively; nI and pI are the revolutions per minute and the pole pairs of PMs of the inner rotor, respectively; εIs and εIo are the small parameters of the electromagnetic coupling fluctuation, respectively, εIs=2ΔkIs/k-Is and εIo=2ΔkIo/k-Io; cc is the conjugate complex of the right expression in (2).The values of small parametersεIs and εIo are 0.041 and 0.0355, respectively. They are rather small and the numerical difference is smaller. Here, εIs and εIo are uniformly simplified with the same parameter ε and ε=0.04.


### 2.2. Parametric Vibration Model and Dynamic Differential Equations

When the outer FP is fixed, the parametric vibration model of EIMG system shown in Figure4 contains four subsystems, namely, the inner stator/inner FP subsystem, the inner FP/inner rotor subsystem, the inner rotor/outer FP subsystem, and the outer FP/outer stator subsystem. The model allows each part to rotate about their translation axes. The electromagnetic coupling dynamic model employs the following assumptions:(1)
Main components of EIMG system are considered to be rigid, assuming that the elastic deformations of all the components are negligible.(2)
Constraints between the inner stator, the inner FP, the outer stator, and the foundation are equivalent to the tangential linear spring, respectively.(3)
Although EIMG system has the maximum output torque, overload does not occur. Manufacturing and installation errors of all components are not included in this paper.(4)
All PMs on the inner and outer surfaces of the inner rotor are assumed to be identical with the same size and performance parameters, respectively. All PMs on the outer stator have the same size and performance parameters too. Magnetic components and nonmagnetic components in the inner and outer FPs have the same sizes and performance parameters, respectively.Figure 4
Dynamic model of the electromechanical integrated magnetic gear.
(a)
The inner stator/inner FP subsystem
(b)
The inner FP/inner rotor subsystem
(c)
The inner rotor/outer FP subsystem
(d)
The outer FP/outer rotor subsystemAfter defining the differential equations of individual subsystem, the differential equations of five degrees of freedom (DOF) of the overall EIMG system can be given in matrix form as(3)Mx¨+Cx˙+Kx=F.The mass matrixM, the damping matrix C, the stiffness matrix K, the load vector F, and the displacement vector x are derived as follows:(4)M=diag⁡MsMIfMIMofMo,C=diag⁡cscIfcIcofco,F=000-ΔTofcos⁡ωetRof0T,x=diag⁡usuIfuIuofuo,K=k-Iscos2⁡αIs+ks-k-Iscos2⁡αIs000-k-Iscos2⁡αIsk-Iscos2⁡αIs+kIIcos2⁡αII+kIf-kIIcos2⁡αII000-kIIcos2⁡αIIkIIcos2⁡αII+k-Iocos2⁡αIo-k-Iocos2⁡αIo000-k-Iocos2⁡αIok-Iocos2⁡αIo+koocos2⁡αoo-koocos2⁡αoo000-koocos2⁡αookoocos2⁡αoo+ko,where Ms, MIf, MI, Mof, and Mo are the masses of the inner stator, the inner FP, the inner rotor, the outer FP, and the outer stator, respectively; us, uIf, uI, uof, and uo are the torsional displacements of the inner stator, the inner FP, the inner rotor, the outer FP, and the outer stator, respectively; cs, cIf, cI, cof, and co are the damping coefficients on the inner stator, the inner FP, the inner rotor, the outer FP, and the outer stator, respectively; αIs is the meshing angle between the inner stator and the inner FP; αII is the meshing angle between the inner FP and the inner rotor; αIo is the meshing angle between the inner rotor and the outer FP; αoo is the meshing angle between the outer FP and the outer stator; ks, kIf, and ko are the supporting stiffnesses among the inner stator, the inner FP, the outer stator, and the foundation, respectively; ωe is the exciting frequency of the torque fluctuation on the outer FP; Rof is the turning radius of the outer FP; ΔTof is the torque increment on the outer FP.Considering the time-varying periodic stiffnesses, EIMG system is a typical parametric vibration system. The dynamic differential equations of the parametric vibration system can be got by substituting (2) into (3):(5)Mx¨+cx˙+Kx+ΔKtx=F,where(6)ΔKt=kIsεeiωstcos2⁡αIs-kIsεeiωstcos2⁡αIs000-kIsεeiωstcos2⁡αIskIsεeiωstcos2⁡αIs00000kIoεeiωotcos2⁡αIo-kIoεeiωotcos2⁡αIo000-kIoεeiωotcos2⁡αIokIoεeiωotcos2⁡αIo000000,in which ΔKt is the increment stiffness matrix.Based on the linear system in (3), (5) can be normalized into the following form:(7)x¨N+cNx˙N+kN+ΔkNxN=ΔFN.In (7), xN, kN, ΔkNt, and ΔFN are the normal displacement vector, the normal stiffness matrix, the normal increment stiffness matrix, and the normal load vector: (8)xN=uN1uN2uN3uN4uN5TΔFN=FN1FN2FN3FN4FN5TkN=diag⁡kN1kN2kN3kN4kN5=diag⁡ω12ω22ω32ω42ω52ΔkNt=⋯ΔkN1,n⋯⋮ΔkNk,n⋮⋯ΔkN5,n⋯,where(9)ΔkNk,n=Ak,neiωst+Bk,neiωot;Ak,n=k-IsεIscos2⁡αIsAN1n-AN2nAN1k-AN2k;Bk,n=k-IoεIocos2⁡αIoAN3n-AN4nAN3k-AN4k;FNk=-AN4kΔTofeiωet2Rof,in which ANk,n, k(n)=1,2,3,4,5, is the element on line k, column n of the normal shape matrix AN.Because the damping matrixc is a diagonal matrix, the elements on the primary diagonal of the normal damping matrix cN are much bigger than other elements. So, the normal damping matrix cN in (7) can be simplified into a diagonal matrix as follows:(10)cN=diag⁡cN1cN2cN3cN4cN5.


## 2.1. Electromagnetic Coupling Stiffnesses of EIMG System

EIMG shown in Figure1 is mainly composed of the inner stator, the inner ferromagnetic pole-pieces (FPs), the inner rotor, the outer FP, and the outer stator. Three-phase coils are arranged in the inner stator and can generate a rotating magnetic field in order to provide power. PMs of the N and S poles 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 take charge of modulating the magnetic fields in two air-gaps beside them in order to make the number of pole pairs of the PMs agree with the number of pole pairs of the space harmonic flux density of the air-gaps.Figure 1
Topology and prototype of the electromechanical integrated magnetic gear.The finite element model of the example EIMG system shown in Table1 can be founded. Here, the material of the inner stator is silicon steel sheets. The back iron materials of the inner rotor and the outer stator are Q235 steel. The materials of the magnetic components and nonmagnetic components in FPs are 23TW250 and epikote, respectively. The properties of the epikote basically agree with air. Meshing is carried out manually and the rapid calculations can be realized by controlling the grid refinement level.Table 1
Parameters of example EIMG system.
Number of pole pairs on the inner surface of the inner rotor
5
Number of pole pairs on the inner stator
2
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 FP/mm
100
Outer radius of the inner FP/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 FP/mm
204
Outer radius of the outer FP/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
Current frequency/Hz
50The meshes between the back iron and PMs of the inner rotor have common boundary nodes. During running of EIMG system, all unit nodes on PMs and the back iron of the inner rotor rotate together as a rigid body. Similarly, the back iron and PMs of the outer stator have common boundary nodes. The meshes on the magnetic components and nonmagnetic components of the inner and outer FPs have common boundary nodes, respectively. The middle positions of the air-gaps are not the common areas of the nodes. The magnetic field distribution is calculated by interpolations among the boundary nodes.The reluctances of the air-gaps are bigger than PMs and the back irons. Meanwhile, the thicknesses of the air-gaps are smaller. So, the grid density in the air-gaps is bigger. Although there are some errors, the errors are negligible if the meshes of the air-gaps are very small. The initial grid model of the example EIMG system and the magnetic field distribution can be figured out in Figure2.Figure 2
Finite element model of the example EIMG system.
(a)
Initial grid model
(b)
Magnetic field distributionThe magnetic flux density in the middle of the air-gaps isBj, where j is the number of four air-gaps, j=1,2,3,4. Bθj and Brj are the tangential and radial components of the magnetic flux density Bj, respectively. The tangential electromagnetic coupling forces among components can be calculated by Maxwell stress tensor method. Then, the electromagnetic coupling stiffnesses among components can be obtained as follows:(1)kit=dFitduit=dFitRjdθ=L/μ0∫02πBθjBrjdθdθ,where Fit and kit, i=Is, II, Io, oo, are the tangential components of the electromagnetic coupling forces and the electromagnetic coupling stiffnesses between the inner stator and the inner FP, the inner FP and the inner rotor, the inner rotor and the outer FP, and the outer FP and the outer stator, respectively; L is the axial effective length of EIMG system; Rj is the radius of the middle position of the air-gap among components; μ0 is the space permeability; θ is the mechanical rotation angle.When the coils in the inner stator are connected to the three-phase alternating currents, the number match of the pole pairs of the PMs is realized by the inner and outer FPs modulating the magnetic fields, while, because of the magnetic fields modulating, the electromagnetic coupling stiffnesses among components vary periodically. The dynamic electromagnetic coupling stiffnesses can be worked out by the finite element method and shown in Figure3.Figure 3
Dynamic electromagnetic coupling stiffnesses among components.
(a)
k
Is
(
t
) and kII(t)
(b)
k
Io
(
t
) and koo(t)
(c)
Frequency component ofkIs(t)
(d)
Frequency component ofkIo(t)Figure3 shows that the stiffnesses fluctuations between the inner stator and the inner FP, the inner rotor and the outer FP, are bigger than others and cannot be ignored. The dynamic electromagnetic coupling stiffnesses contain multiple harmonics. The main frequency of the stiffness fluctuation between the inner stator and the inner FP is the product of the current frequency and the number of pole pairs of the three-phase alternating currents on the inner stator. The main frequency of the stiffness fluctuation between the inner rotor and the outer FP is the product of the angular frequency and the number of pole pairs of the PMs on the outer surface of the inner rotor. Considering the amplitudes of other frequency components are much smaller, the time-varying electromagnetic coupling stiffnesses kIst and kIIt can be simplified as follows:(2)kIst=k-Is+2ΔkIscos⁡ωst=k-Is1+εIseωst+cc,kIot=k-Io+2ΔkIocos⁡ωot=k-Io1+εIoeωot+cc,where k-Is and k-Io are the average values of the electromagnetic coupling stiffnesses kIst and kIot, respectively; ΔkIs and ΔkIo are the increments of the stiffnesses kIst and kIot, respectively; ωs and ωo are the fluctuation frequencies of the stiffness increments ΔkIs and ΔkIo, respectively; namely, meshing frequencies ωs=2πnsps/60, ωo=2πnIpI/60; ns and ps are the frequency and the pole pairs of the three-phase alternating currents on the inner stator, respectively; nI and pI are the revolutions per minute and the pole pairs of PMs of the inner rotor, respectively; εIs and εIo are the small parameters of the electromagnetic coupling fluctuation, respectively, εIs=2ΔkIs/k-Is and εIo=2ΔkIo/k-Io; cc is the conjugate complex of the right expression in (2).The values of small parametersεIs and εIo are 0.041 and 0.0355, respectively. They are rather small and the numerical difference is smaller. Here, εIs and εIo are uniformly simplified with the same parameter ε and ε=0.04.


## 2.2. Parametric Vibration Model and Dynamic Differential Equations

When the outer FP is fixed, the parametric vibration model of EIMG system shown in Figure4 contains four subsystems, namely, the inner stator/inner FP subsystem, the inner FP/inner rotor subsystem, the inner rotor/outer FP subsystem, and the outer FP/outer stator subsystem. The model allows each part to rotate about their translation axes. The electromagnetic coupling dynamic model employs the following assumptions:(1)
Main components of EIMG system are considered to be rigid, assuming that the elastic deformations of all the components are negligible.(2)
Constraints between the inner stator, the inner FP, the outer stator, and the foundation are equivalent to the tangential linear spring, respectively.(3)
Although EIMG system has the maximum output torque, overload does not occur. Manufacturing and installation errors of all components are not included in this paper.(4)
All PMs on the inner and outer surfaces of the inner rotor are assumed to be identical with the same size and performance parameters, respectively. All PMs on the outer stator have the same size and performance parameters too. Magnetic components and nonmagnetic components in the inner and outer FPs have the same sizes and performance parameters, respectively.Figure 4
Dynamic model of the electromechanical integrated magnetic gear.
(a)
The inner stator/inner FP subsystem
(b)
The inner FP/inner rotor subsystem
(c)
The inner rotor/outer FP subsystem
(d)
The outer FP/outer rotor subsystemAfter defining the differential equations of individual subsystem, the differential equations of five degrees of freedom (DOF) of the overall EIMG system can be given in matrix form as(3)Mx¨+Cx˙+Kx=F.The mass matrixM, the damping matrix C, the stiffness matrix K, the load vector F, and the displacement vector x are derived as follows:(4)M=diag⁡MsMIfMIMofMo,C=diag⁡cscIfcIcofco,F=000-ΔTofcos⁡ωetRof0T,x=diag⁡usuIfuIuofuo,K=k-Iscos2⁡αIs+ks-k-Iscos2⁡αIs000-k-Iscos2⁡αIsk-Iscos2⁡αIs+kIIcos2⁡αII+kIf-kIIcos2⁡αII000-kIIcos2⁡αIIkIIcos2⁡αII+k-Iocos2⁡αIo-k-Iocos2⁡αIo000-k-Iocos2⁡αIok-Iocos2⁡αIo+koocos2⁡αoo-koocos2⁡αoo000-koocos2⁡αookoocos2⁡αoo+ko,where Ms, MIf, MI, Mof, and Mo are the masses of the inner stator, the inner FP, the inner rotor, the outer FP, and the outer stator, respectively; us, uIf, uI, uof, and uo are the torsional displacements of the inner stator, the inner FP, the inner rotor, the outer FP, and the outer stator, respectively; cs, cIf, cI, cof, and co are the damping coefficients on the inner stator, the inner FP, the inner rotor, the outer FP, and the outer stator, respectively; αIs is the meshing angle between the inner stator and the inner FP; αII is the meshing angle between the inner FP and the inner rotor; αIo is the meshing angle between the inner rotor and the outer FP; αoo is the meshing angle between the outer FP and the outer stator; ks, kIf, and ko are the supporting stiffnesses among the inner stator, the inner FP, the outer stator, and the foundation, respectively; ωe is the exciting frequency of the torque fluctuation on the outer FP; Rof is the turning radius of the outer FP; ΔTof is the torque increment on the outer FP.Considering the time-varying periodic stiffnesses, EIMG system is a typical parametric vibration system. The dynamic differential equations of the parametric vibration system can be got by substituting (2) into (3):(5)Mx¨+cx˙+Kx+ΔKtx=F,where(6)ΔKt=kIsεeiωstcos2⁡αIs-kIsεeiωstcos2⁡αIs000-kIsεeiωstcos2⁡αIskIsεeiωstcos2⁡αIs00000kIoεeiωotcos2⁡αIo-kIoεeiωotcos2⁡αIo000-kIoεeiωotcos2⁡αIokIoεeiωotcos2⁡αIo000000,in which ΔKt is the increment stiffness matrix.Based on the linear system in (3), (5) can be normalized into the following form:(7)x¨N+cNx˙N+kN+ΔkNxN=ΔFN.In (7), xN, kN, ΔkNt, and ΔFN are the normal displacement vector, the normal stiffness matrix, the normal increment stiffness matrix, and the normal load vector: (8)xN=uN1uN2uN3uN4uN5TΔFN=FN1FN2FN3FN4FN5TkN=diag⁡kN1kN2kN3kN4kN5=diag⁡ω12ω22ω32ω42ω52ΔkNt=⋯ΔkN1,n⋯⋮ΔkNk,n⋮⋯ΔkN5,n⋯,where(9)ΔkNk,n=Ak,neiωst+Bk,neiωot;Ak,n=k-IsεIscos2⁡αIsAN1n-AN2nAN1k-AN2k;Bk,n=k-IoεIocos2⁡αIoAN3n-AN4nAN3k-AN4k;FNk=-AN4kΔTofeiωet2Rof,in which ANk,n, k(n)=1,2,3,4,5, is the element on line k, column n of the normal shape matrix AN.Because the damping matrixc is a diagonal matrix, the elements on the primary diagonal of the normal damping matrix cN are much bigger than other elements. So, the normal damping matrix cN in (7) can be simplified into a diagonal matrix as follows:(10)cN=diag⁡cN1cN2cN3cN4cN5.


## 3. Resonances of EIMG System

### 3.1. Main Resonance of EIMG System

The multiscale method is adopted to solve (7). Meanwhile, in order to balance the effects of the damping forces and the increment stiffnesses and make them into the same perturbation equation, the following assumptions are adopted:(11)xNk=xNk0T0,T1+εxNk1T0,T1,cNk=εcNk′,ΔTof=εΔTof′,where Tn=εnt.By substituting (11) into (7) and making the same power coefficients of the small parameters ε on both sides of the equations equal, the following differential equations can be obtained:Zero power is(12)D02xNk0+ωk2xNk0=0.First power is(13)D02xNk1+ωk2xNk1=-2D0D1xNk0-cNk′D0xNk0+∑n=15Ak,neiωst+Bk,neiωotxNk0-AN4iΔTo′eiωet2Rof+cc.The general solution of (12) can be expressed:(14)xNk0=Ekeiωkt+cc,where Ek is the initial amplitude, which can be calculated by initial energy.By substituting (14) into (13), (13) can be rewritten as (15)D02xNk1+ωk2xNk1=-2iωkE˙keiωkt-cNk′iωkEkeiωkt-AN4kΔTo′eiωet2Rof+∑n=15Ak,neiωst+Bk,neiωotEkeiωkt+cc+cc.When the exciting frequency of the torque fluctuation on the outer FP is close to the natural frequencies of the derived EIMG system, the main resonances will occur. Here, the following assumption is introduced:(16)ωe=ω3+εσ1.By substituting (16) into (15) and eliminating the secular terms, the following differential equations can be obtained:(17)-2iω3E˙3-cN3′iω3E3-AN43ΔTo′eiσ1T12Rof=0,-2iωmE˙m-cNm′iωmEm=0,m=1,2,4,5.The solution of (17) can be worked out by the method of constant variation and can be expressed as follows:(18)E3=G3e-cN3′T1/2+iAN43ΔTo′ejσ1T12Rofω3cN3′+2iσ1+cc,Em=Gme-cNm′T1/2,where Gm is a constant.Equation (18) can be rewritten as(19)E3=G3e-cN3′T1/2+AN43ΔTo′eiθ1+σ1T12ω3RofcN3′2+4σ12+cc,where sin⁡θ1=cN3′/cN3′2+4σ12, cos⁡θ1=2σ1/cN3′2+4σ12.In (19), cNk′e-cNk′T1 will be close to zero with time increasing. The stable responses of the zero order approximate analytic solution of EIMG system can be given:(20)xN30=AN43ΔTo′eiθ1+σ1T1eiω3t2ω3RofcN3′2+4σ12+cc=AN43ΔTo′cos⁡θ1+σ1T1+ω3tω3RofcN3′2+4σ12,xNm0=0.By substituting (19) and (20) into (15), the stable responses of the first order approximate analytic solution of EIMG system can be calculated: (21)xN31=∑k=15EkA3,keiωs±ωktωs+ωk2-ω32+EkB3,keiωo±ωktωo+ωk2-ω32+cc=AN43ΔTo′A33eiωs±ω3t+θ1+σ1T12ω3Rofωs2+2ωsω3cN3′2+4σ12+AN43ΔTo′B33eiωo±ω3t+θ1+σ1T12ω3Rofωo2+2ωoω3cN3′2+4σ12+cc,xNm1=∑k=15EnA3,keiωs±ωktωs+ωk2-ωk2+EnB3,keiωo±ωktωo+ωk2-ωk2+AN4,mΔTo′eiωet2Rofωe2-ωm2+cc=AN43ΔTo′A33eiωs±ω3t+θ1+σ1T12ω3Rofωs2+2ωsω3cN3′2+4σ12+AN43ΔTo′B33eiωo±ω3t+θ1+σ1T12ω3Rofωo2+2ωoω3cN3′2+4σ12+AN4,mΔTo′eiωet2Rofωe2-ωm2+cc.Equation (21) can be simplified into(22)xN31=AN43ΔTo′A33cos⁡ωst±ω3t+θ1+σ1T12ω3Rofωs2+2ωsω3cN3′2+4σ12+AN43ΔTo′B33cos⁡ωot±ω3t+θ1+σ1T12ω3Rofωo2+2ωoω3cN3′2+4σ12+cc,xNm1=AN43ΔTo′A33cos⁡ωst±ω3t+θ1+σ1T12ω3Rofωs2+2ωsω3cN3′2+4σ12+AN43ΔTo′B33cos⁡ωot±ω3t+θ1+σ1T12ω3Rofωo2+2ωoω3cN3′2+4σ12+AN4,mΔTo′cos⁡ωet2Rofωe2-ωm2+cc.Similarly, when the exciting frequency is close to a certain frequencyωk, the resonance will occur and the stable responses of EIMG system in normal coordinate system can be deduced:(23)xNk0=AN4kΔTo′cos⁡θ1+σ1T1+ωktωkRofcNk′2+4σ12+cc,xNm0=0,xNk1=AN4kΔTo′Akkcos⁡ωst±ωkt+θ1+σ1T12ωkRofωs2+2ωsωkcNk′2+4σ12+AN4kΔTo′Bkkcos⁡ωot±ωkt+θ1+σ1T12ωkRofωo2+2ωoωkcNk′2+4σ12+ccxNm1=AN4kΔTo′Akkcos⁡ωst±ωkt+θ1+σ1T12ωkRofωs2+2ωsωkcNk′2+4σ12·AN4mΔTo′cos⁡ωet2Rofωe2-ωm2+AN4kΔTo′Bkkcos⁡ωot±ωkt+θ1+σ1T12ωkRofωo2+2ωoωkcNk′2+4σ12+cc,where m is any natural number except for k, k=1,2,3,4,5.By substituting (23) into (11), the approximate analytical solution of EIMG system can be got by(24)x=ANxN0+εxN1.


### 3.2. Combination Resonance of EIMG System

By the multiscale method, the following assumptions are adopted:(25)xNk=xNk0T0,T1+εxNk1T0,T1,cNk=εcNk′,ωe=ω3+ωo+εσ2.By substituting (25) into (7) and making the same power coefficients of the small parameters ε on both sides of the equations equal, the following differential equations can be obtained:Zero power is(26)D02xNk0+ωk2xNk0=AN4kΔToeiωet2Rof+cc.The first power is(27)D02xNk1+ωk2xNk1=-2D0D1xNk0-cNk′D0xNk0+∑n=15Ak,neiωst+Bk,neiωotxNk0+cc,k=1,2,…,5,n=1,2,…,5.The general solution of (26) can be expressed:(28)xNk0=Ekeiωkt+AN4kΔToeiωet2Rofωe2-ωk2+cc.By substituting (28) into (27) and eliminating the secular terms, the following differential equations can be obtained:(29)-2iω3dE3dT1-cN3′iω3E3+AN43ΔToB33eiσ2T12Rofωe2-ω32=0,-2iωmdEmdT1-cNm′iωmEm=0.The solution of (29) can be given:(30)E3=F3e-cN3′T1+AN4kΔToB33ejθ2+σ2T12ω3Rofωe2-ω32cN3′2+4σ22+cc,Em=Fme-cNm′T1,where sin⁡θ2=cN3′/cN3′2+4σ22, cos⁡θ2=2σ2/cN3′2+4σ22.In (30), Fke-cNk′T1 will gradually be close to zero with time going on. The stable responses of the zero order approximate analytic solutions of EIMG system can be given:(31)xN30=AN43ΔTo2Rofωe2-ω32B33eiω3t+θ2+σT1ω3cN3′2+4σ22+eiωet+cc,xNm0=AN4mΔToeiωet2Rofωe2-ωm2+cc.By substituting (28) and (31) into (27), the stable responses of the first order approximate analytic solution of (27) can be calculated:(32)
u
N
31
=
∑
n
=
1
5
A
3
n
E
n
e
i
ω
s
±
ω
j
t
ω
s
±
ω
n
2
-
ω
3
2
+
B
3
n
E
n
e
i
ω
o
±
ω
j
t
ω
o
±
ω
n
2
-
ω
3
2
+
A
43
Δ
T
o
A
3
n
e
i
ω
s
±
ω
e
t
2
R
o
f
ω
e
2
-
ω
n
2
ω
s
±
ω
e
2
-
ω
3
2
+
∑
m
=
1,2
,
4,5
A
43
Δ
T
o
B
3
m
e
i
ω
o
±
ω
e
t
2
R
o
f
ω
e
2
-
ω
m
2
ω
o
±
ω
e
2
-
ω
3
2
+
cc
=
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
3
2
A
33
B
33
e
i
ω
s
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
s
±
ω
3
2
-
ω
3
2
+
B
33
2
e
i
ω
o
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
o
±
ω
3
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
A
3
n
e
i
ω
s
±
ω
e
t
2
R
o
f
ω
e
2
-
ω
n
2
ω
s
±
ω
e
2
-
ω
3
2
+
∑
m
=
1,2
,
4,5
A
N
43
Δ
T
o
B
3
m
e
i
ω
o
±
ω
j
t
2
R
o
f
ω
e
2
-
ω
m
2
ω
o
±
ω
m
2
-
ω
3
2
+
cc
u
N
m
1
=
∑
n
=
1
5
A
3
n
E
n
e
i
ω
s
±
ω
j
t
ω
s
±
ω
n
2
-
ω
3
2
+
B
3
n
E
n
e
i
ω
o
±
ω
j
t
ω
o
±
ω
n
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
e
i
ω
s
±
ω
e
t
ω
s
±
ω
e
2
-
ω
3
2
+
B
3
n
e
i
ω
o
±
ω
e
t
ω
o
±
ω
e
2
-
ω
3
2
+
cc
=
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
3
2
A
33
B
33
e
i
ω
s
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
s
±
ω
3
2
-
ω
3
2
+
B
33
2
e
i
ω
o
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
o
±
ω
3
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
e
i
ω
s
±
ω
e
t
ω
s
±
ω
e
2
-
ω
3
2
+
B
3
n
e
i
ω
o
±
ω
e
t
ω
o
±
ω
e
2
-
ω
3
2
+
cc
.Equation (32) can be simplified into(33)
u
N
31
=
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
3
2
A
33
B
33
cos
⁡
ω
s
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
s
±
ω
3
2
-
ω
3
2
+
B
33
2
cos
⁡
ω
o
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
o
±
ω
3
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
A
3
n
cos
⁡
ω
s
±
ω
j
t
2
R
o
f
ω
e
2
-
ω
n
2
ω
s
±
ω
n
2
-
ω
3
2
+
∑
m
=
1,2
,
4,5
A
N
43
Δ
T
o
B
3
m
cos
⁡
ω
o
±
ω
m
t
2
R
o
f
ω
e
2
-
ω
m
2
ω
o
±
ω
m
2
-
ω
3
2
+
cc
,
u
N
m
1
=
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
3
2
A
33
B
33
cos
⁡
ω
s
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
s
±
ω
3
2
-
ω
3
2
+
B
33
2
cos
⁡
ω
o
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
o
±
ω
3
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
cos
⁡
ω
s
±
ω
n
t
ω
s
±
ω
n
2
-
ω
3
2
+
B
3
n
cos
⁡
ω
o
±
ω
n
t
ω
o
±
ω
n
2
-
ω
3
2
+
cc
.Similarly, when the exciting frequency is close to a certain frequencyωk, k=1,2,3,4,5, the stable responses of EIMG system in normal coordinate system can be deduced:(34)
x
N
k
0
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
B
k
k
e
i
ω
k
t
+
θ
2
+
σ
T
1
ω
k
c
N
k
′
2
+
4
σ
2
2
+
e
i
ω
e
t
+
cc
,
x
N
m
0
=
A
N
4
m
Δ
T
o
e
i
ω
e
t
2
R
o
f
ω
e
2
-
ω
m
2
+
cc
,
u
N
k
1
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
A
k
k
B
k
k
cos
⁡
ω
s
±
ω
k
t
+
θ
2
+
ε
σ
2
t
+
ω
k
t
ω
k
c
N
k
′
2
+
4
σ
2
2
ω
s
±
ω
k
2
-
ω
k
2
+
B
k
k
2
cos
⁡
ω
o
±
ω
k
t
+
θ
2
+
ε
σ
2
t
+
ω
k
t
ω
k
c
N
k
′
2
+
4
σ
2
2
ω
o
±
ω
k
2
-
ω
k
2
+
∑
n
=
1
5
A
4
k
Δ
T
o
A
k
k
cos
⁡
ω
s
±
ω
n
t
2
R
o
f
ω
e
2
-
ω
n
2
ω
s
±
ω
n
2
-
ω
k
2
+
∑
m
≠
k
A
4
k
Δ
T
o
B
m
m
cos
⁡
ω
o
±
ω
m
t
2
R
o
ω
e
2
-
ω
m
2
ω
o
±
ω
m
2
-
ω
k
2
+
cc
,
u
N
m
1
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
A
m
m
B
k
k
cos
⁡
ω
s
±
ω
k
t
+
θ
2
+
ε
σ
2
t
+
ω
k
t
ω
m
c
N
m
′
2
+
4
σ
2
2
ω
s
±
ω
n
2
-
ω
k
2
+
B
m
m
B
k
k
cos
⁡
ω
o
±
ω
k
t
+
θ
2
+
ε
σ
2
t
+
ω
k
t
ω
m
c
N
m
′
2
+
4
σ
2
2
ω
o
±
ω
m
2
-
ω
k
2
+
∑
n
=
1
5
A
4
n
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
cos
⁡
ω
s
±
ω
n
t
ω
s
±
ω
n
2
-
ω
k
2
+
B
3
n
cos
⁡
ω
o
±
ω
m
t
ω
o
±
ω
n
2
-
ω
k
2
+
cc
,where m is any natural number except for k.Meanwhile, the combination resonance responses can be worked out when the exciting frequency is close toωk+ωs:(35)xNk0=AN4kΔTo2Rofωe2-ωk2Akkeiωkt+θ3+σT1ωkcNk′2+4σ32+eiωet+cc,xNm0=AN4nΔToeiωet2Rofωe2-ωn2+cc,where ωe=ωk+ωs+εσ3, sin⁡θ3=cNk′/cNk′2+4σ32, cos⁡θ3=2σ3/cNk′2+4σ32. Consider(36)
u
N
k
1
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
A
k
k
2
cos
⁡
ω
s
±
ω
k
t
+
θ
3
+
ε
σ
2
t
+
ω
k
t
ω
k
c
N
k
′
2
+
4
σ
2
2
ω
s
±
ω
k
2
-
ω
k
2
+
A
k
k
B
k
k
cos
⁡
ω
o
±
ω
k
t
+
θ
3
+
ε
σ
2
t
+
ω
k
t
ω
k
c
N
k
′
2
+
4
σ
2
2
ω
o
±
ω
k
2
-
ω
k
2
+
∑
m
≠
k
A
4
k
Δ
T
o
A
m
m
cos
⁡
ω
s
±
ω
m
t
2
R
o
f
ω
e
2
-
ω
m
2
ω
s
±
ω
m
2
-
ω
k
2
+
∑
n
=
1
5
A
4
k
Δ
T
o
B
k
k
cos
⁡
ω
o
±
ω
n
t
2
R
o
ω
e
2
-
ω
n
2
ω
o
±
ω
n
2
-
ω
k
2
+
cc
,
u
N
m
1
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
A
k
k
2
cos
⁡
ω
s
±
ω
k
t
+
θ
3
+
ε
σ
2
t
+
ω
k
t
ω
m
c
N
m
′
2
+
4
σ
2
2
ω
s
±
ω
m
2
-
ω
k
2
+
A
k
k
B
k
k
cos
⁡
ω
o
±
ω
k
t
+
θ
3
+
ε
σ
2
t
+
ω
k
t
ω
m
c
N
m
′
2
+
4
σ
2
2
ω
o
±
ω
m
2
-
ω
k
2
+
∑
n
=
1
5
A
4
n
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
cos
⁡
ω
s
±
ω
n
t
ω
s
±
ω
n
2
-
ω
k
2
+
B
3
n
cos
⁡
ω
o
±
ω
n
t
ω
o
±
ω
n
2
-
ω
k
2
+
cc
.


## 3.1. Main Resonance of EIMG System

The multiscale method is adopted to solve (7). Meanwhile, in order to balance the effects of the damping forces and the increment stiffnesses and make them into the same perturbation equation, the following assumptions are adopted:(11)xNk=xNk0T0,T1+εxNk1T0,T1,cNk=εcNk′,ΔTof=εΔTof′,where Tn=εnt.By substituting (11) into (7) and making the same power coefficients of the small parameters ε on both sides of the equations equal, the following differential equations can be obtained:Zero power is(12)D02xNk0+ωk2xNk0=0.First power is(13)D02xNk1+ωk2xNk1=-2D0D1xNk0-cNk′D0xNk0+∑n=15Ak,neiωst+Bk,neiωotxNk0-AN4iΔTo′eiωet2Rof+cc.The general solution of (12) can be expressed:(14)xNk0=Ekeiωkt+cc,where Ek is the initial amplitude, which can be calculated by initial energy.By substituting (14) into (13), (13) can be rewritten as (15)D02xNk1+ωk2xNk1=-2iωkE˙keiωkt-cNk′iωkEkeiωkt-AN4kΔTo′eiωet2Rof+∑n=15Ak,neiωst+Bk,neiωotEkeiωkt+cc+cc.When the exciting frequency of the torque fluctuation on the outer FP is close to the natural frequencies of the derived EIMG system, the main resonances will occur. Here, the following assumption is introduced:(16)ωe=ω3+εσ1.By substituting (16) into (15) and eliminating the secular terms, the following differential equations can be obtained:(17)-2iω3E˙3-cN3′iω3E3-AN43ΔTo′eiσ1T12Rof=0,-2iωmE˙m-cNm′iωmEm=0,m=1,2,4,5.The solution of (17) can be worked out by the method of constant variation and can be expressed as follows:(18)E3=G3e-cN3′T1/2+iAN43ΔTo′ejσ1T12Rofω3cN3′+2iσ1+cc,Em=Gme-cNm′T1/2,where Gm is a constant.Equation (18) can be rewritten as(19)E3=G3e-cN3′T1/2+AN43ΔTo′eiθ1+σ1T12ω3RofcN3′2+4σ12+cc,where sin⁡θ1=cN3′/cN3′2+4σ12, cos⁡θ1=2σ1/cN3′2+4σ12.In (19), cNk′e-cNk′T1 will be close to zero with time increasing. The stable responses of the zero order approximate analytic solution of EIMG system can be given:(20)xN30=AN43ΔTo′eiθ1+σ1T1eiω3t2ω3RofcN3′2+4σ12+cc=AN43ΔTo′cos⁡θ1+σ1T1+ω3tω3RofcN3′2+4σ12,xNm0=0.By substituting (19) and (20) into (15), the stable responses of the first order approximate analytic solution of EIMG system can be calculated: (21)xN31=∑k=15EkA3,keiωs±ωktωs+ωk2-ω32+EkB3,keiωo±ωktωo+ωk2-ω32+cc=AN43ΔTo′A33eiωs±ω3t+θ1+σ1T12ω3Rofωs2+2ωsω3cN3′2+4σ12+AN43ΔTo′B33eiωo±ω3t+θ1+σ1T12ω3Rofωo2+2ωoω3cN3′2+4σ12+cc,xNm1=∑k=15EnA3,keiωs±ωktωs+ωk2-ωk2+EnB3,keiωo±ωktωo+ωk2-ωk2+AN4,mΔTo′eiωet2Rofωe2-ωm2+cc=AN43ΔTo′A33eiωs±ω3t+θ1+σ1T12ω3Rofωs2+2ωsω3cN3′2+4σ12+AN43ΔTo′B33eiωo±ω3t+θ1+σ1T12ω3Rofωo2+2ωoω3cN3′2+4σ12+AN4,mΔTo′eiωet2Rofωe2-ωm2+cc.Equation (21) can be simplified into(22)xN31=AN43ΔTo′A33cos⁡ωst±ω3t+θ1+σ1T12ω3Rofωs2+2ωsω3cN3′2+4σ12+AN43ΔTo′B33cos⁡ωot±ω3t+θ1+σ1T12ω3Rofωo2+2ωoω3cN3′2+4σ12+cc,xNm1=AN43ΔTo′A33cos⁡ωst±ω3t+θ1+σ1T12ω3Rofωs2+2ωsω3cN3′2+4σ12+AN43ΔTo′B33cos⁡ωot±ω3t+θ1+σ1T12ω3Rofωo2+2ωoω3cN3′2+4σ12+AN4,mΔTo′cos⁡ωet2Rofωe2-ωm2+cc.Similarly, when the exciting frequency is close to a certain frequencyωk, the resonance will occur and the stable responses of EIMG system in normal coordinate system can be deduced:(23)xNk0=AN4kΔTo′cos⁡θ1+σ1T1+ωktωkRofcNk′2+4σ12+cc,xNm0=0,xNk1=AN4kΔTo′Akkcos⁡ωst±ωkt+θ1+σ1T12ωkRofωs2+2ωsωkcNk′2+4σ12+AN4kΔTo′Bkkcos⁡ωot±ωkt+θ1+σ1T12ωkRofωo2+2ωoωkcNk′2+4σ12+ccxNm1=AN4kΔTo′Akkcos⁡ωst±ωkt+θ1+σ1T12ωkRofωs2+2ωsωkcNk′2+4σ12·AN4mΔTo′cos⁡ωet2Rofωe2-ωm2+AN4kΔTo′Bkkcos⁡ωot±ωkt+θ1+σ1T12ωkRofωo2+2ωoωkcNk′2+4σ12+cc,where m is any natural number except for k, k=1,2,3,4,5.By substituting (23) into (11), the approximate analytical solution of EIMG system can be got by(24)x=ANxN0+εxN1.


## 3.2. Combination Resonance of EIMG System

By the multiscale method, the following assumptions are adopted:(25)xNk=xNk0T0,T1+εxNk1T0,T1,cNk=εcNk′,ωe=ω3+ωo+εσ2.By substituting (25) into (7) and making the same power coefficients of the small parameters ε on both sides of the equations equal, the following differential equations can be obtained:Zero power is(26)D02xNk0+ωk2xNk0=AN4kΔToeiωet2Rof+cc.The first power is(27)D02xNk1+ωk2xNk1=-2D0D1xNk0-cNk′D0xNk0+∑n=15Ak,neiωst+Bk,neiωotxNk0+cc,k=1,2,…,5,n=1,2,…,5.The general solution of (26) can be expressed:(28)xNk0=Ekeiωkt+AN4kΔToeiωet2Rofωe2-ωk2+cc.By substituting (28) into (27) and eliminating the secular terms, the following differential equations can be obtained:(29)-2iω3dE3dT1-cN3′iω3E3+AN43ΔToB33eiσ2T12Rofωe2-ω32=0,-2iωmdEmdT1-cNm′iωmEm=0.The solution of (29) can be given:(30)E3=F3e-cN3′T1+AN4kΔToB33ejθ2+σ2T12ω3Rofωe2-ω32cN3′2+4σ22+cc,Em=Fme-cNm′T1,where sin⁡θ2=cN3′/cN3′2+4σ22, cos⁡θ2=2σ2/cN3′2+4σ22.In (30), Fke-cNk′T1 will gradually be close to zero with time going on. The stable responses of the zero order approximate analytic solutions of EIMG system can be given:(31)xN30=AN43ΔTo2Rofωe2-ω32B33eiω3t+θ2+σT1ω3cN3′2+4σ22+eiωet+cc,xNm0=AN4mΔToeiωet2Rofωe2-ωm2+cc.By substituting (28) and (31) into (27), the stable responses of the first order approximate analytic solution of (27) can be calculated:(32)
u
N
31
=
∑
n
=
1
5
A
3
n
E
n
e
i
ω
s
±
ω
j
t
ω
s
±
ω
n
2
-
ω
3
2
+
B
3
n
E
n
e
i
ω
o
±
ω
j
t
ω
o
±
ω
n
2
-
ω
3
2
+
A
43
Δ
T
o
A
3
n
e
i
ω
s
±
ω
e
t
2
R
o
f
ω
e
2
-
ω
n
2
ω
s
±
ω
e
2
-
ω
3
2
+
∑
m
=
1,2
,
4,5
A
43
Δ
T
o
B
3
m
e
i
ω
o
±
ω
e
t
2
R
o
f
ω
e
2
-
ω
m
2
ω
o
±
ω
e
2
-
ω
3
2
+
cc
=
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
3
2
A
33
B
33
e
i
ω
s
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
s
±
ω
3
2
-
ω
3
2
+
B
33
2
e
i
ω
o
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
o
±
ω
3
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
A
3
n
e
i
ω
s
±
ω
e
t
2
R
o
f
ω
e
2
-
ω
n
2
ω
s
±
ω
e
2
-
ω
3
2
+
∑
m
=
1,2
,
4,5
A
N
43
Δ
T
o
B
3
m
e
i
ω
o
±
ω
j
t
2
R
o
f
ω
e
2
-
ω
m
2
ω
o
±
ω
m
2
-
ω
3
2
+
cc
u
N
m
1
=
∑
n
=
1
5
A
3
n
E
n
e
i
ω
s
±
ω
j
t
ω
s
±
ω
n
2
-
ω
3
2
+
B
3
n
E
n
e
i
ω
o
±
ω
j
t
ω
o
±
ω
n
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
e
i
ω
s
±
ω
e
t
ω
s
±
ω
e
2
-
ω
3
2
+
B
3
n
e
i
ω
o
±
ω
e
t
ω
o
±
ω
e
2
-
ω
3
2
+
cc
=
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
3
2
A
33
B
33
e
i
ω
s
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
s
±
ω
3
2
-
ω
3
2
+
B
33
2
e
i
ω
o
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
o
±
ω
3
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
e
i
ω
s
±
ω
e
t
ω
s
±
ω
e
2
-
ω
3
2
+
B
3
n
e
i
ω
o
±
ω
e
t
ω
o
±
ω
e
2
-
ω
3
2
+
cc
.Equation (32) can be simplified into(33)
u
N
31
=
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
3
2
A
33
B
33
cos
⁡
ω
s
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
s
±
ω
3
2
-
ω
3
2
+
B
33
2
cos
⁡
ω
o
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
o
±
ω
3
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
A
3
n
cos
⁡
ω
s
±
ω
j
t
2
R
o
f
ω
e
2
-
ω
n
2
ω
s
±
ω
n
2
-
ω
3
2
+
∑
m
=
1,2
,
4,5
A
N
43
Δ
T
o
B
3
m
cos
⁡
ω
o
±
ω
m
t
2
R
o
f
ω
e
2
-
ω
m
2
ω
o
±
ω
m
2
-
ω
3
2
+
cc
,
u
N
m
1
=
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
3
2
A
33
B
33
cos
⁡
ω
s
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
s
±
ω
3
2
-
ω
3
2
+
B
33
2
cos
⁡
ω
o
±
ω
3
t
+
θ
2
+
ε
σ
2
t
+
ω
3
t
ω
3
c
N
3
′
2
+
4
σ
2
2
ω
o
±
ω
3
2
-
ω
3
2
+
∑
n
=
1
5
A
N
43
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
cos
⁡
ω
s
±
ω
n
t
ω
s
±
ω
n
2
-
ω
3
2
+
B
3
n
cos
⁡
ω
o
±
ω
n
t
ω
o
±
ω
n
2
-
ω
3
2
+
cc
.Similarly, when the exciting frequency is close to a certain frequencyωk, k=1,2,3,4,5, the stable responses of EIMG system in normal coordinate system can be deduced:(34)
x
N
k
0
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
B
k
k
e
i
ω
k
t
+
θ
2
+
σ
T
1
ω
k
c
N
k
′
2
+
4
σ
2
2
+
e
i
ω
e
t
+
cc
,
x
N
m
0
=
A
N
4
m
Δ
T
o
e
i
ω
e
t
2
R
o
f
ω
e
2
-
ω
m
2
+
cc
,
u
N
k
1
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
A
k
k
B
k
k
cos
⁡
ω
s
±
ω
k
t
+
θ
2
+
ε
σ
2
t
+
ω
k
t
ω
k
c
N
k
′
2
+
4
σ
2
2
ω
s
±
ω
k
2
-
ω
k
2
+
B
k
k
2
cos
⁡
ω
o
±
ω
k
t
+
θ
2
+
ε
σ
2
t
+
ω
k
t
ω
k
c
N
k
′
2
+
4
σ
2
2
ω
o
±
ω
k
2
-
ω
k
2
+
∑
n
=
1
5
A
4
k
Δ
T
o
A
k
k
cos
⁡
ω
s
±
ω
n
t
2
R
o
f
ω
e
2
-
ω
n
2
ω
s
±
ω
n
2
-
ω
k
2
+
∑
m
≠
k
A
4
k
Δ
T
o
B
m
m
cos
⁡
ω
o
±
ω
m
t
2
R
o
ω
e
2
-
ω
m
2
ω
o
±
ω
m
2
-
ω
k
2
+
cc
,
u
N
m
1
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
A
m
m
B
k
k
cos
⁡
ω
s
±
ω
k
t
+
θ
2
+
ε
σ
2
t
+
ω
k
t
ω
m
c
N
m
′
2
+
4
σ
2
2
ω
s
±
ω
n
2
-
ω
k
2
+
B
m
m
B
k
k
cos
⁡
ω
o
±
ω
k
t
+
θ
2
+
ε
σ
2
t
+
ω
k
t
ω
m
c
N
m
′
2
+
4
σ
2
2
ω
o
±
ω
m
2
-
ω
k
2
+
∑
n
=
1
5
A
4
n
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
cos
⁡
ω
s
±
ω
n
t
ω
s
±
ω
n
2
-
ω
k
2
+
B
3
n
cos
⁡
ω
o
±
ω
m
t
ω
o
±
ω
n
2
-
ω
k
2
+
cc
,where m is any natural number except for k.Meanwhile, the combination resonance responses can be worked out when the exciting frequency is close toωk+ωs:(35)xNk0=AN4kΔTo2Rofωe2-ωk2Akkeiωkt+θ3+σT1ωkcNk′2+4σ32+eiωet+cc,xNm0=AN4nΔToeiωet2Rofωe2-ωn2+cc,where ωe=ωk+ωs+εσ3, sin⁡θ3=cNk′/cNk′2+4σ32, cos⁡θ3=2σ3/cNk′2+4σ32. Consider(36)
u
N
k
1
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
A
k
k
2
cos
⁡
ω
s
±
ω
k
t
+
θ
3
+
ε
σ
2
t
+
ω
k
t
ω
k
c
N
k
′
2
+
4
σ
2
2
ω
s
±
ω
k
2
-
ω
k
2
+
A
k
k
B
k
k
cos
⁡
ω
o
±
ω
k
t
+
θ
3
+
ε
σ
2
t
+
ω
k
t
ω
k
c
N
k
′
2
+
4
σ
2
2
ω
o
±
ω
k
2
-
ω
k
2
+
∑
m
≠
k
A
4
k
Δ
T
o
A
m
m
cos
⁡
ω
s
±
ω
m
t
2
R
o
f
ω
e
2
-
ω
m
2
ω
s
±
ω
m
2
-
ω
k
2
+
∑
n
=
1
5
A
4
k
Δ
T
o
B
k
k
cos
⁡
ω
o
±
ω
n
t
2
R
o
ω
e
2
-
ω
n
2
ω
o
±
ω
n
2
-
ω
k
2
+
cc
,
u
N
m
1
=
A
N
4
k
Δ
T
o
2
R
o
f
ω
e
2
-
ω
k
2
A
k
k
2
cos
⁡
ω
s
±
ω
k
t
+
θ
3
+
ε
σ
2
t
+
ω
k
t
ω
m
c
N
m
′
2
+
4
σ
2
2
ω
s
±
ω
m
2
-
ω
k
2
+
A
k
k
B
k
k
cos
⁡
ω
o
±
ω
k
t
+
θ
3
+
ε
σ
2
t
+
ω
k
t
ω
m
c
N
m
′
2
+
4
σ
2
2
ω
o
±
ω
m
2
-
ω
k
2
+
∑
n
=
1
5
A
4
n
Δ
T
o
2
R
o
f
ω
e
2
-
ω
n
2
A
3
n
cos
⁡
ω
s
±
ω
n
t
ω
s
±
ω
n
2
-
ω
k
2
+
B
3
n
cos
⁡
ω
o
±
ω
n
t
ω
o
±
ω
n
2
-
ω
k
2
+
cc
.


## 4. Results and Discussion

### 4.1. Modal Characteristics of EIMG System

The dynamic parameters of the example EIMG system shown in Table1 can be given in Table 2. By substituting the dynamic parameters into (3) and neglecting the external load, the modal frequencies and corresponding modal shapes can be calculated, which are shown in Table 3.Table 2
Characteristic parameters of the example EIMG system.
k
oo (N/m)
k
Io (N/m)
k
II (N/m)
α
Io (°)
α
oo (°)
8.6444 × 105
3.1446 × 105
4.5467 × 105
0.1818
0.1014
k
Is (N/m)
k
If (N/m)
k
o (N/m)
α
II (°)
α
Is (°)
2.0416 × 105
3 × 106
3 × 106
0.1435
0.2022
c
s (N/(m/s))
c
If (N/(m/s))
c
of (N/(m/s))
M
s (kg)
M
If (kg)
0.1
0.1
0.05
1.3
1.3
c
I (N/(m/s))
c
o (N/(m/s))
Δ
T
of (N·m)
M
of (kg)
M
o (kg)
0.05
0.1
5
6.2
7.7
M
I (kg)
R
o
11
219Table 3
Modal frequencies and modal shapes of EIMG system.
Mode type
IRRM
OFRM
OSRM
ISRM
IFRM
Natural frequencies (rad/s)
277.8
456.4
827.7
1322.5
1624.3
Mode shapes
0.0489
−0.0333
0.0090
1.0000
0.3952
0.2860
−0.1806
0.0348
0.3804
−1.0000
1.0000
−0.5953
0.0892
−0.0195
0.0330
0.6747
1.0000
−0.9161
0.0014
−0.0014
0.3705
0.8125
1.0000
−0.0003
0.0001Table3 illustrates that there are five different modes, in which the natural frequencies and the modal shapes are completely different. When the natural frequencies increase, the relative displacements of the inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP come to the maximum in sequence. So, the corresponding modal shapes of the natural frequencies from low to high are named the rotational modes of the inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP (IRRM, OFRM, ORSM, IRSM, and IFRM), respectively. For example, when the natural frequency is 456.4 rad/s, the relative displacement of the outer FP is maximal and the corresponding mode shape is named rotational mode of the outer FP.The inner stator, the inner FP, and the outer stator are fixed on the frame of EIMG system. The inner rotor and the outer FP are coupled with other components by the electromagnetic coupling forces. Because the electromagnetic coupling stiffnesses are much smaller than the mechanical support stiffnesses, the modal frequencies of ROSM, RISM, and RIFM are bigger than the modal frequencies of RIRM and ROFM.


### 4.2. Main Resonances of EIMG System

By substituting the parameters of the example EIMG system into (24), the main resonances can be drawn and shown in Figure 5.Figure 5
Main resonance of the EIMG system.
(a)
(b)
(c)
(d)
(e)The resonance amplitudes vary widely in the main resonance of the EIMG system. Whenωe is close to ω1, except the inner rotor, the resonance displacements of the inner and the outer FP are bigger than others. When ωe is close to ω2, except the outer FP, the resonance displacements of the inner rotor and the outer stator are bigger than others. When ωe is close to ω3, the resonance displacements of the outer FP and the outer stator are much bigger than others. When ωe is close to ω4 or ω5, the resonance displacements of the inner stator and the inner FP are much bigger than others.When the exciting frequency is close to the rotational modal frequencies of the inner rotor, the outer FP, and the outer stator, the resonance amplitudes are much bigger than other resonances of EIMG system.


### 4.3. Combination Resonances of EIMG System

By substituting the parameters of the example EIMG system into (34) to (36), the time and frequency responses of the combination resonances can be drawn and shown in Figures 6 and 7.Figure 6
Combination resonance of the EIMG system withωe≈ωk+ωo.Figure 7
Combination resonance of the EIMG system withωe≈ωk+ωs.Figure6 shows that the dynamic displacements of the inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP will reach maximum, when the exciting frequency increases from small to big and the combination resonances occur. The dominant frequency of the combination resonances is the natural frequency ωk of EIMG system rather than the exciting frequency ωe.When the combination occurs, the mode energies of IRRM, OFRM, ORSM, IRSM, and IFRM will be motivated, respectively. The amplitudes of inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP will reach maximum, respectively. Considering the electromagnetic coupling interaction, the amplitudes of components will have a big difference, because the adjacent components are affected greatly and the separated components are affected less. For example, whenωe≈ω1+ωo, the mode energy of IRRM will be motivated and the amplitude of the inner rotor will be maximal. Meanwhile, the amplitudes of the inner and outer FP are bigger and the resonances of the inner and outer stator are relatively smaller. Compared with the main resonances, the amplitudes in combination resonances are much smaller.When the exciting frequency is close to the combination frequency between a certain frequency and the meshing frequencyωs, namely, ωe≈ωk+ωs, the combination resonances have similar characteristics with ωe≈ωk+ωo. Meanwhile, the resonance amplitudes are smaller. Because the inner stator is not adjacent with the outer FP and the outer exciting resource is on the outer FP, the influence of the meshing stiffness wave ΔkIs on the outer FP is very small. So the combination resonance will occur when the exciting frequency ωe must be very close to the combination frequency ωk+ωs and the resonance amplitude is very small.


### 4.4. Discussion on Resonances of EIMG System

All the components in the EIMG system are coupled by the electromagnetic fields and the electromagnetic coupling stiffnesses are much smaller than the supporting stiffnesses of the components. Because the interactions among all the components are soft-contact, the adjacent components are affected greatly and the separated components are affected less. So, the relative amplitudes of components in mode shapes have great difference. Meanwhile, the relative amplitudes of components in main resonances and the combination resonances have large difference too.Considering the mode shape characteristic of EIMG system, the equivalent loads on the inner rotor, the outer FP, and the outer stator are much bigger than the equivalent loads on the inner stator and the inner FP. Meanwhile, because the outer FP is the output component, the mode energies of OFRM, IRRM, and OFRM are motivated easily and the corresponding main resonances are relatively bigger.When the excited frequency is close to the combination frequencies and the natural frequencies, the mode energies will be motivated and the corresponding modal frequency will stand out. So the mode frequency, namely, the natural frequency, is the dominant frequency.Combination resonances must be excited by the inner and outer exciting, and the value of the combination frequency is bigger. So, the combination resonances are not relatively easy to take place and the amplitudes of components are smaller. The resonance amplitudes caused by the combination exciting of the outer torque wave and the electromagnetic coupling stiffness between the inner stator and inner FP are smaller and are much smaller.Considering other multiple harmonics components in the electromagnetic coupling stiffnesses, there are more resonance sources. In order to improve the dynamic characteristics, structure parameter optimizing and the dynamic parameter optimizing must be carried out.


## 4.1. Modal Characteristics of EIMG System

The dynamic parameters of the example EIMG system shown in Table1 can be given in Table 2. By substituting the dynamic parameters into (3) and neglecting the external load, the modal frequencies and corresponding modal shapes can be calculated, which are shown in Table 3.Table 2
Characteristic parameters of the example EIMG system.
k
oo (N/m)
k
Io (N/m)
k
II (N/m)
α
Io (°)
α
oo (°)
8.6444 × 105
3.1446 × 105
4.5467 × 105
0.1818
0.1014
k
Is (N/m)
k
If (N/m)
k
o (N/m)
α
II (°)
α
Is (°)
2.0416 × 105
3 × 106
3 × 106
0.1435
0.2022
c
s (N/(m/s))
c
If (N/(m/s))
c
of (N/(m/s))
M
s (kg)
M
If (kg)
0.1
0.1
0.05
1.3
1.3
c
I (N/(m/s))
c
o (N/(m/s))
Δ
T
of (N·m)
M
of (kg)
M
o (kg)
0.05
0.1
5
6.2
7.7
M
I (kg)
R
o
11
219Table 3
Modal frequencies and modal shapes of EIMG system.
Mode type
IRRM
OFRM
OSRM
ISRM
IFRM
Natural frequencies (rad/s)
277.8
456.4
827.7
1322.5
1624.3
Mode shapes
0.0489
−0.0333
0.0090
1.0000
0.3952
0.2860
−0.1806
0.0348
0.3804
−1.0000
1.0000
−0.5953
0.0892
−0.0195
0.0330
0.6747
1.0000
−0.9161
0.0014
−0.0014
0.3705
0.8125
1.0000
−0.0003
0.0001Table3 illustrates that there are five different modes, in which the natural frequencies and the modal shapes are completely different. When the natural frequencies increase, the relative displacements of the inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP come to the maximum in sequence. So, the corresponding modal shapes of the natural frequencies from low to high are named the rotational modes of the inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP (IRRM, OFRM, ORSM, IRSM, and IFRM), respectively. For example, when the natural frequency is 456.4 rad/s, the relative displacement of the outer FP is maximal and the corresponding mode shape is named rotational mode of the outer FP.The inner stator, the inner FP, and the outer stator are fixed on the frame of EIMG system. The inner rotor and the outer FP are coupled with other components by the electromagnetic coupling forces. Because the electromagnetic coupling stiffnesses are much smaller than the mechanical support stiffnesses, the modal frequencies of ROSM, RISM, and RIFM are bigger than the modal frequencies of RIRM and ROFM.


## 4.2. Main Resonances of EIMG System

By substituting the parameters of the example EIMG system into (24), the main resonances can be drawn and shown in Figure 5.Figure 5
Main resonance of the EIMG system.
(a)
(b)
(c)
(d)
(e)The resonance amplitudes vary widely in the main resonance of the EIMG system. Whenωe is close to ω1, except the inner rotor, the resonance displacements of the inner and the outer FP are bigger than others. When ωe is close to ω2, except the outer FP, the resonance displacements of the inner rotor and the outer stator are bigger than others. When ωe is close to ω3, the resonance displacements of the outer FP and the outer stator are much bigger than others. When ωe is close to ω4 or ω5, the resonance displacements of the inner stator and the inner FP are much bigger than others.When the exciting frequency is close to the rotational modal frequencies of the inner rotor, the outer FP, and the outer stator, the resonance amplitudes are much bigger than other resonances of EIMG system.


## 4.3. Combination Resonances of EIMG System

By substituting the parameters of the example EIMG system into (34) to (36), the time and frequency responses of the combination resonances can be drawn and shown in Figures 6 and 7.Figure 6
Combination resonance of the EIMG system withωe≈ωk+ωo.Figure 7
Combination resonance of the EIMG system withωe≈ωk+ωs.Figure6 shows that the dynamic displacements of the inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP will reach maximum, when the exciting frequency increases from small to big and the combination resonances occur. The dominant frequency of the combination resonances is the natural frequency ωk of EIMG system rather than the exciting frequency ωe.When the combination occurs, the mode energies of IRRM, OFRM, ORSM, IRSM, and IFRM will be motivated, respectively. The amplitudes of inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP will reach maximum, respectively. Considering the electromagnetic coupling interaction, the amplitudes of components will have a big difference, because the adjacent components are affected greatly and the separated components are affected less. For example, whenωe≈ω1+ωo, the mode energy of IRRM will be motivated and the amplitude of the inner rotor will be maximal. Meanwhile, the amplitudes of the inner and outer FP are bigger and the resonances of the inner and outer stator are relatively smaller. Compared with the main resonances, the amplitudes in combination resonances are much smaller.When the exciting frequency is close to the combination frequency between a certain frequency and the meshing frequencyωs, namely, ωe≈ωk+ωs, the combination resonances have similar characteristics with ωe≈ωk+ωo. Meanwhile, the resonance amplitudes are smaller. Because the inner stator is not adjacent with the outer FP and the outer exciting resource is on the outer FP, the influence of the meshing stiffness wave ΔkIs on the outer FP is very small. So the combination resonance will occur when the exciting frequency ωe must be very close to the combination frequency ωk+ωs and the resonance amplitude is very small.


## 4.4. Discussion on Resonances of EIMG System

All the components in the EIMG system are coupled by the electromagnetic fields and the electromagnetic coupling stiffnesses are much smaller than the supporting stiffnesses of the components. Because the interactions among all the components are soft-contact, the adjacent components are affected greatly and the separated components are affected less. So, the relative amplitudes of components in mode shapes have great difference. Meanwhile, the relative amplitudes of components in main resonances and the combination resonances have large difference too.Considering the mode shape characteristic of EIMG system, the equivalent loads on the inner rotor, the outer FP, and the outer stator are much bigger than the equivalent loads on the inner stator and the inner FP. Meanwhile, because the outer FP is the output component, the mode energies of OFRM, IRRM, and OFRM are motivated easily and the corresponding main resonances are relatively bigger.When the excited frequency is close to the combination frequencies and the natural frequencies, the mode energies will be motivated and the corresponding modal frequency will stand out. So the mode frequency, namely, the natural frequency, is the dominant frequency.Combination resonances must be excited by the inner and outer exciting, and the value of the combination frequency is bigger. So, the combination resonances are not relatively easy to take place and the amplitudes of components are smaller. The resonance amplitudes caused by the combination exciting of the outer torque wave and the electromagnetic coupling stiffness between the inner stator and inner FP are smaller and are much smaller.Considering other multiple harmonics components in the electromagnetic coupling stiffnesses, there are more resonance sources. In order to improve the dynamic characteristics, structure parameter optimizing and the dynamic parameter optimizing must be carried out.


## 5. Conclusion

Considering the magnetic field modulating, the electromagnetic coupling stiffnesses in EIMG system vary periodically and contain multiple harmonics, in which the main harmonic is the product of the angular frequency and the number of pole pairs. Considering the time-varying periodic stiffnesses, EIMG system is a typical parametric vibration system.Considering the exciting on the outer FP and the time-varying electromagnetic coupling stiffnesses, there are the main resonances and the combination resonances in EIMG system. When the main resonances and the combination resonances occur, the dominant frequencies are always the natural frequency of EIMG system. When the exciting frequency increases from small to big, the rotational displacements of the inner rotor, the outer FP, the outer stator, the inner stator, and the inner FP will reach maxima.The electromagnetic coupling stiffnesses are much smaller than the mechanical supporting stiffness. The relative resonance amplitudes of components have great difference. When the exciting frequency is close to the natural frequency of OFRM, IRRM, and OFRM, the resonances of EIMG system are bigger. The component amplitudes of the main resonances are much bigger than the components amplitudes of the combination resonances. Particularly, the meshing stiffness wave between the inner stator and the inner FP has little influence on EIMG system.The time-varying electromagnetic coupling stiffnesses make EIMG system be unstable. So the stability charts are the important research and must be discussed in next work.



---

*Source: 572937-2015-07-30.xml*