# The Combined Internal and Principal Parametric Resonances on Continuum Stator System of Asynchronous Machine

Authors: Baizhou Li; Qichang Zhang
Journal: Shock and Vibration
 (2014)

---

## Abstract

With the increasing requirement of quiet electrical machines in the civil and defense industry, it is very significant and necessary to predict the vibration and noise characteristics of stator and rotor in the early conceptual phase. Therefore, the combined internal and principal parametric resonances of a stator system excited by radial electromagnetic force are presented in this paper. The stator structure is modeled as a continuum double-shell system which is loaded by a varying distributed electromagnetic load. The nonlinear dynamic equations are derived and solved by the method of multiple scales. The influences of mechanical and electromagnetic parameters on resonance characteristics are illustrated by the frequency-response curves. Furthermore, the Runge-Kutta method is adopted to numerically analyze steady-state response for the further understanding of the resonance characteristics with different parameters.


---

## Body

## 1. Introduction

As the electromagnetic conversion device between mechanical and electrical energy, electrical machines are widely used in the civil and defense industry. The electromagnetic vibration of the stator system, which is one of the main noise sources, can directly affect the noise level of electrical machines. Therefore, with the development of quiet electrical machines, it is very significant to study electromagnetic vibration mechanism of the stator system for noise reduction.Consideration of the levels of vibration and noise at the design stage has been of importance for many years, and it is still a major consideration for machine designers and operators. It is well known that vibration and noise can be reduced to a very large extent when the electromagnetic force is not allowed to excite any kind of resonances of electrical machines. In fact, the electromagnetic load still contains many frequency components and it will excite radial vibrations of the stator system and generate strong noise. Moreover, the resonance responses of the stator system make the electromagnetic load change by affecting the air-gap magnetic field. Therefore, the interaction mechanism between mechanical and electromagnetic behavior is of great concern in this paper.Many researches have been developed in electromagnetic vibration on the stator system of multifarious electrical machines over the years in the analytical and the numerical methods. Some of them pay attention to establishing the models of stator and rotor and calculating the natural characteristics by finite element method (FEM) which could provide accurate solutions for dynamic system with complicated stator and rotor shapes [1–5]. The purpose is to achieve lower levels of vibrations by designing the electromagnetic force to have excitation frequencies away from natural frequencies of stator and rotor. And the models have been developed from rigid bodies to elastic bodies. However, some of the researchers studied the origins of electromagnetic excitation which are the Maxwell forces, the Lorentz forces, and the magnetostrictive forces [6, 7]. The Maxwell forces are the predominant forces that act on the boundary regions from the air-gap to the stator core [8–13]. In particular, Pennacchi and Frosini presented a model to calculate the magnetic pull for the dynamical behavior analysis of large-size generators [9]. And the method proposed in [9] has been validated using experimental results in [10]. Furthermore, nonlinear effects due to electromechanical interaction have been analyzed based on the proposed method in [11]. Therefore, the analysis of the vibration characteristics of the stator system under Maxwell forces is very essential.On the other hand, for the continuum vibration, Qing et al. [14] analyzed the natural frequencies and mode shapes by a semianalytical method based on thick double-shell systems. Luo and Wang [12] established a continuum beam model for the nonlinear vibration of rotor-bearing system. Nguyen [15] and Hegazy [16] presented the nonlinear resonances involving two mode shapes of thin plates. It has been shown that when a parametric resonance is excited in the presence of an internal resonance, the coincidence of these two types of resonances will give rise to the simultaneous resonances.In spite of publications on the electromagnetic vibration of the electrical machines, very few of them focused on the continuum stator system under the coupling electromagnetic forces between the electromagnetic and displacement fields. In this paper, the complicated stator system of a three-phase asynchronous machine is simplified to a continuous elastic structure with two coaxial cylindrical shells. The electromagnetic force expression with the effect of resonance response of the stator system is given. Based on the energy principle the nonlinear dynamic equation is derived. And the method of multiple scales is employed to develop the nonlinear dynamic analysis of the continuous system. The jump phenomenon is pointed out, and the influences of detuning parameters, damping coefficients, and electromagnetic parameters are illustrated by the frequency-response curves. The steady-state responses of the combined internal and principal parametric resonances are obtained with the Runge-Kutta method. Diagrams of time history, phase chart, the fast Fourier transformation (FFT), and Poincaré section are presented to confirm both the effects of different parameters and the response of the system.


## 2. Distributed Electromagnetic Forces

The electromagnetic force contains radial and tangential components. The radial component can excite the radial vibration of the stator core, and it is the main source of the electromagnetic noise. The tangential component, which is the reaction force corresponding to the electromagnetic torque, bends the tooth root, and it is the minor source of the electromagnetic noise. Therefore, it is generally accepted that the tangential component is negligible in the analysis of radial vibration. The radial electromagnetic force acts on the stator core and generates the radial displacementws of the stator core. Meanwhile, the radial displacement changes the thickness of the air-gap and then changes the air-gap permeance Λ as
(1)Λx,θ=Λ01+W=Λ01-W+W2-W3+W4+⋯,
where Λ0 is the average air-gap permeance. W=ws/δ0 is the dimensionless radial displacement of the stator core, and δ0 is the average eccentricity. It is well known that the dynamic eccentricity and mechanical imbalance can produce the rotational frequency vibrations. But the eccentricity is neglected because it merely modifies the expression of (1) and will not change the mechanism of producing vibrations. Then, the radial electromagnetic force [8–13] is proportional to the square of the magnetic flux-density which is the product of the magnetomotive force and air-gap permeance, and it can be derived as
(2)qr=Z0+Z1ws+Z2ws2+Z3ws3+Z4cos⁡⁡pθ-Ωt+Z5wscos⁡⁡pθ-Ωt+Z6ws2cos⁡⁡pθ-Ωt+Z7ws3cos⁡(pθ-Ωt)+⋯,
where
(3)Z0=14μ0Bδ2+Bν2+Bμ2,Z1=-2Z0δ0,Z2=3Z0δ02,Z3=-4Z0δ03,Z4=12μ0BνBμ,Z5=-2Z4δ0,Z6=3Z4δ02,Z7=-4Z4δ03.
And Ω=ων±ωμ are the frequencies of electromagnetic force corresponding to p=ν±μ. Bδ, Bν, and Bμ are fundamental components, the harmonic components of stator windings and rotor windings, respectively. The phase angles are neglected in vibration analysis. Equation (2) shows that there is mutual interaction between the electromagnetic and mechanical behaviors. That is the major problem focused on in this paper.


## 3. Double-Shell Stator Model

Several stator models are presented in previous works for kinds of stator structures, such as single thin circular ring model, two concentric circular rings model, single thick cylindrical shell model [3], and two coaxial cylindrical shells named double-shell system here. A three-phase asynchronous machine is considered in this paper; its stator core and frame are interconnected by key bars. The length of the stator core is shorter than that of the frame. Thus, the double-shell stator model can be established as in Figure 1.Double-shell stator system.
(a)
(b)Based on Hamilton’s Principle, the displacement variational principle for nonlinear vibration of thin elastic shell can be expressed as(4)δ∫0tKE-PEdt+∫0t∬sqxδu+qθδv+qrδwRdxdθdt=0,
where KE is kinetic energy, and PE is potential energy of thin shell. qx, qθ, and qr are components of electromagnetic force, and x, θ, and r are cylindrical coordinates, respectively. The coincident part of the double-shell stator system can be considered as being fixed together by a supposed force P(x,θ)sinωit. Then, the force conditions of the stator core and the frame can be written as
(5)Stator core:qrsx,θ,t=-RαR1Px,θsinωit+RβR1qrx,θ,t,iiiiiiFrame:qrrx,θ,t=RαR1Px,θsinωit,
where Rβ is the centroid radius of the stator core. The superscript s is for the stator core and r is for the frame. The displacement solutions can be written as the products of periodic function Ti(t) and modal shape functions Ui,Vi,Wi determined by the boundary condition as
(6)ux,θ,t=∑iTitUix,θ,vx,θ,t=∑iTitVix,θ,wx,θ,t=∑iTitWix,θ,
where i is orders of mode shapes. Equation (6) can be used to determine the displacement components of the stator core and the frame, respectively. Considering the orthogonality of mode shapes and the displacement boundary condition of ws≡wr in the coincident region s1, the vibration equation with modal damping can be obtained as
(7)T¨it+2μiT˙it+ωi2Tit=RβMis+Mir∬s1qrx,θ,tWisdxdθ,
where μi is modal damping coefficient, and ωi is the natural frequency of double-shell system.


## 4. Nonlinear Dynamic Equations

Substituting (2) into (7), the nonlinear dynamic equations of the stator system with the interaction between the radial displacement of the stator core and the air-gap magnetic field are obtained as(8a)T¨it+2μiT˙it+ω-i2Tit=ηi1∑j∑kΓijkTjtTkt+ηi4∑j∑kΓijkTj(t)Tk(t)cos⁡Ωt+⋯(i≠p),(8b)T¨p(t)+2μpT˙p(t)+ω-p2Tp(t)=ηp1∑j∑kΓpjkTjtTkt+ηp4∑j∑kΓppjkTj(t)Tk(t)cos⁡Ωt+2Fpcos⁡Ωt+⋯0000000000000000000000000000000000000000(i=p),where j,k are orders of mode shapes distinguishing from i, and
(9)ω-i2=ωi2-RβMis+Mir∬s1Z1wsWisdxdθ
are the frequencies of the stator system in working state. Γijk are the coupling coefficients which present the interaction of different mode shapes. And ηin and Γijk are expressed as
(10)ηi,p,1=1δ023Z0l1RβMi,p,s+Mi,p,r,ηi,p,4=1δ023Z4l1RβMi,p,s+Mi,p,r,2Fp=πZ4l1RβMps+Mpr,Γijk=∫02πcos⁡iθ·cos⁡jθ·cos⁡kθdθ.
Equations (8a) and (8b) show that the forced excitation term exists in the dynamic equations if and only if i=p. Therefore, both of the two equations are necessary for the combined internal and principal parametric resonances. All the other nonlinear terms are not important for the steady-state analysis and they will decay with time.The relationship of internal resonance can be expressed by introducing detuning parameterδ1, and the nearness of excitation frequency Ω to one natural frequency ω-p can be expressed by introducing the external detuning parameter δ2 according to
(11)ω-4=ω-1+ω-3+εδ1,Ω=ω-4+εδ2.We introduce small parameterε≪1 to transform the scale μi→εμi′ and ηin→εηin′, and the method of multiple scales [17] is applied to solve (8a) and (8b) as
(12)Tit,ε=Ti0T0,T1+εTi1T0,T1i=1,3,4,
where T0=t is a fast time scale and T1=εt is a slow time scale. Inserting (12) into (8a) and (8b), and separating at each order of ε, one obtains(13a)O(ε0):D02Ti0+ω-i2Ti0=0,O(ε1):D02Ti1+ω-i2Ti1=-2D0D1Ti0-2μi′D0Ti0+ηi1′∑j∑kΓijkTj0Tk0+ηi4′∑j∑kΓpijkTj0Tk0cos⁡ΩT0+⋯i≠p,(13b)O(ε0):D02Tp0+ω-p2Tp0=0,O(ε1):D02Tp1+ω-p2Tp1=-2D0D1Tp0-2μp′D0Tp0+ηp1′∑j∑kΓpjkTj0Tk0+ηp4′∑j∑kΓppjkTj0Tk0cos⁡ΩT0+2Fp′cos⁡ΩT0+⋯000000000000000000000000000000000000000(i=p),where Dn=∂/∂Tn. At O(ε0) the solution can be written in the following convenient form:
(14)Ti0T0,T1=aiT12exp⁡jωiT0+ϕiT1+c.c.,
where ai and ϕi are real. ai are the amplitudes; ϕi are the phases. c.c. stands for complex conjugate and j2=-1.Substituting (14) into the ε1-order equations in (13a) and (13b) leads to the differential equations, which can be written as a six-dimensional dynamical system by separating real and imaginary parts as
(15)a1′=μ1a1-η114ω-1Γ143a3a4sinγ1,a1ϕ1′=-η114ω-1Γ143a3a4cos⁡γ1,a3′=μ3a3-η314ω-3Γ341a1a4sinγ1,a3ϕ3′=-η314ω-3Γ341a1a4cos⁡γ1,a4′=μ4a4+η414ω-4Γ413a1a3sinγ1-η444ω-4Γ4444a42sinγ2-F4ω-4sinγ2,a4ϕ4′=-η414ω-4Γ413a1a3cos⁡γ1-η444ω-4Γ4444a42cos⁡γ2-F4ω-4cos⁡γ2,
where ·′ stands for the derivation with respect to T1, γ1=δ1T1-ϕ1-ϕ3+ϕ4, and γ2=δ2T1-ϕ4.For the steady-state solution, by canceling the time derivatives (a1′=a3′=a4′=γ1′=γ2′=0), one can obtain that only two kinds of solutions are possible. In the absence of the internal resonance, the forced excitation can excite only one mode at a time. In this case, the amplitude of the stator system can be only expressed by the fourth mode as
(16)a1=0,a3=0,Λ442a44+2ΛΛ44-μ42-εδ22a42+Λ2=0.But in the presence of the internal resonance, the combined internal and principal parametric resonances occur, and the amplitude of the stator system can be obtained as(17)a12=1Λ3Λ4μ3μ1P1±P22-P32,a32=μ1Λ3μ3Λ1a12,a42=μ1μ3Λ1Λ31+εδ1+εδ22μ1+μ32,
where
(18)P1=μ1μ3μ1+μ3εδ2εδ1+εδ2-μ4μ1+μ3,P2=Λ1Λ3Λ+Λ44a42,P3=μ1μ3μ1+μ3εδ2μ1+μ3+μ4εδ1+εδ2,(19)Λ1=η114ω-1Γ143,Λ3=η314ω-3Γ341,Λ4=η414ω-4Γ413,Λ44=η444ω-4Γ4444,Λ=F4ω-4.The stability analysis can be derived by computing the Jacobian matrix of (15) and the existence condition of real solutions. If all the real parts of the eigenvalues from the Jacobian matrix are negative, the corresponding motion is stable; otherwise it is unstable. Finally, the following stability condition of the continuous system is obtained as
(20)La=P1±P22-P32>0,Lb=P22-P32≥0,Lc=μ1μ3Λ1Λ31+(εδ1+εδ2)2(μ1+μ3)21/2≥a4.


## 5. Numerical Simulation and Discussion

In this section, the fourth-order Runge-Kutta method is used to solve (15) for more insight into the combined internal and principal parametric resonances by assuming μ1=μ3=μ4=μ. The stationary frequency-response curves are illustrated in Figures 2–4. In these figures, the ordinate a represents the steady-state amplitude and the abscissa εδ denotes the frequency detuning parameter. Solid and broken lines in these figures represent the stable and unstable solutions, respectively.The frequency-response curves. (a)δ1=0, (b) δ1>0, and (c) δ1<0.
(a)
(b)
(c)The frequency-response curves. (a)δ2=0, (b) δ2>0, and (c) δ2<0.
(a)
(b)
(c)The frequency-response curves.
(a)
(b)
(c)
(d)


### 5.1. Influence of Detuning Parameters

Figure2 shows the stationary frequency-response curves with three cases of internal resonance detuning parameter δ1. All the modal damping coefficients are assumed to be 2. When the internal resonance detuning parameter δ1 is zero, the response curves are symmetric and present M-shape. If δ1>0, the symmetry axis moves to the left based on (17) and inclines towards right. And if δ1<0, the symmetry axis moves to the right and inclines towards left. This means the fourth mode is forcedly and parametrically excited, and the presences of the first and third modes are possible only by transferring energy from the fourth mode to them through internal mechanism. The peak value of a4, which originates from the combined forced and parametric resonance, is independent of the internal resonance detuning parameter δ1. Once the amplitude a4 is larger than Lc, the solution becomes unstable. The internal resonance detuning parameter δ1 has a direct influence on the energy transfer when the simultaneous resonances are excited.Figure3 shows the stationary frequency-response curves with three cases of forced resonance detuning parameter δ2. Similarly, only when δ2=0, the response curves are symmetric. When a1=a3=0, there is only the fourth mode of the stator system left under parametrically forced resonance with high amplitude. If δ2>0, the right of a1 and a3 curves rises, and all the curves move towards left based on (17). On the contrary, if δ2<0, the left of a1 and a3 curves rises, and all the curves move towards right. This means that the amplitudes of a1 and a3 increase since the simultaneous resonances occur and then drop to zero when the energy transfer vanishes in Figure 3(b). But the amplitudes of a1 and a3 jump to a certain value when simultaneous resonances occur and then decrease to zero when the energy transfer vanishes in Figure 3(c).


### 5.2. Influence of Damping Coefficients

Take the frequency-response curves of forced resonance detuning parameter as an example to discuss the influence of damping coefficients. Both the resonance region and the amplitudes of three modes are seen reduced due to the increasing of the damping coefficients in Figure4. If μ<2, a4 presents an unstable double-peak curve. And the fourth mode is on combined forced and parametric resonances. When the combined internal and principal parametric resonances are excited, the stable resonance amplitude of a4 increases slowly, the unstable one decreases, and the resonance amplitudes of a1 and a3 reduce gradually with the increase of damping coefficients. When the damping coefficients increase to 6, the resonance amplitude of a4 has been larger than the other modes. Until the damping coefficients increase to 8, there is only the fourth mode left for the parametrically forced resonance with the absence of the first and third modes. The reason of this phenomenon is that the energy transferred from the fourth mode to the other modes is reducing with the damping coefficients increasing.


### 5.3. Influence of Electromagnetic Parameters

Figure5 illustrates the relationship between the amplitude and the magnetic flux-density. The simultaneous resonances can be excited at a certain value. The amplitudes are decreasing with the fundamental component of the magnetic flux-density increasing, which presents the fundamental component as a damping for the resonance of the stator system in Figure 5(a). According to (2), the influences of the stator and rotor windings harmonic components are similar. And Figure 5(b) shows that the stator windings harmonic components make the amplitudes of the first and third modes increase in the resonance region. That means the harmonic magnetic field of the stator and rotor is the main source of electromagnetic excitation.The amplitude-magnetic flux-density response curves.
(a)
(b)


### 5.4. Dynamic Characteristics

In the following figures, diagrams of time history, phase chart, FFT, and Poincaré section are presented for further understanding of the combined internal and principal parametric resonances of the stator system. Based on the parameters in Figures2–5, by assuming εδ1=0.12, εδ2=0.48, μ=8, Bδ=0.75T, and Bν=Bμ=0.1T, Figures 6(a) through 6(d) confirm that the system is period-1 motion according to Figure 4(d). The phase chart is a closed circle, and the energy of the fourth mode is much larger than that of the first and third modes. There is only one nonzero point, and the other two zero points coincide with each other in the Poincaré section. That means the stator system is in stable single-frequency and periodic motion. With the damping coefficients reducing from 8 to 6, there are three points in Figure 7(d) and three frequencies in Figure 7(c). The energy of the fourth mode transfers to the first and third modes, but still higher than the other modes. The stator system is in quasiperiod motion. When the damp coefficients reduce from 6 to 2, the energy of the first mode is the highest in Figure 8. And the phase chart presents reciprocating motion in a bounded domain. With the amplitude of electromagnetic excitation increasing to Bδ=0.5T and Bν=Bμ=0.2T, the energy of all three modes is higher than the former, and the amplitude of the stator system increases. Figures 9(a) through 9(d) confirm the conclusion of Figure 5.Dynamic characteristic withεδ1=0.12, εδ2=0.48, μ=8, Bδ=0.75T, and Bμ=0.1T.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré sectionDynamic characteristic withεδ1=0.12, εδ2=0.48, μ=6, Bδ=0.75T, and Bμ=0.1T.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré sectionDynamic characteristic withεδ1=0.12, εδ2=0.48, μ=2, Bδ=0.75T, and Bμ=0.1T.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré sectionDynamic characteristic withεδ1=0.12, εδ2=0.48, μ=2, Bδ=0.5T, and Bμ=0.2T.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré sectionIt is essential to point out that whenεδ1=0.52, there is an additional relation between the first and fourth natural frequencies as ω-4≈3ω-1. But the coupling coefficient Γ4111=0 which means there is only frequency relationship, and no interaction between these two modes. The expression in (11) is still effective. Figures 10(a) through 10(d) show that the stator system presents the period-2 motion.Dynamic characteristic withεδ1=0.52, εδ2=0.48, and ω-4≈3ω-1.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré section


## 5.1. Influence of Detuning Parameters

Figure2 shows the stationary frequency-response curves with three cases of internal resonance detuning parameter δ1. All the modal damping coefficients are assumed to be 2. When the internal resonance detuning parameter δ1 is zero, the response curves are symmetric and present M-shape. If δ1>0, the symmetry axis moves to the left based on (17) and inclines towards right. And if δ1<0, the symmetry axis moves to the right and inclines towards left. This means the fourth mode is forcedly and parametrically excited, and the presences of the first and third modes are possible only by transferring energy from the fourth mode to them through internal mechanism. The peak value of a4, which originates from the combined forced and parametric resonance, is independent of the internal resonance detuning parameter δ1. Once the amplitude a4 is larger than Lc, the solution becomes unstable. The internal resonance detuning parameter δ1 has a direct influence on the energy transfer when the simultaneous resonances are excited.Figure3 shows the stationary frequency-response curves with three cases of forced resonance detuning parameter δ2. Similarly, only when δ2=0, the response curves are symmetric. When a1=a3=0, there is only the fourth mode of the stator system left under parametrically forced resonance with high amplitude. If δ2>0, the right of a1 and a3 curves rises, and all the curves move towards left based on (17). On the contrary, if δ2<0, the left of a1 and a3 curves rises, and all the curves move towards right. This means that the amplitudes of a1 and a3 increase since the simultaneous resonances occur and then drop to zero when the energy transfer vanishes in Figure 3(b). But the amplitudes of a1 and a3 jump to a certain value when simultaneous resonances occur and then decrease to zero when the energy transfer vanishes in Figure 3(c).


## 5.2. Influence of Damping Coefficients

Take the frequency-response curves of forced resonance detuning parameter as an example to discuss the influence of damping coefficients. Both the resonance region and the amplitudes of three modes are seen reduced due to the increasing of the damping coefficients in Figure4. If μ<2, a4 presents an unstable double-peak curve. And the fourth mode is on combined forced and parametric resonances. When the combined internal and principal parametric resonances are excited, the stable resonance amplitude of a4 increases slowly, the unstable one decreases, and the resonance amplitudes of a1 and a3 reduce gradually with the increase of damping coefficients. When the damping coefficients increase to 6, the resonance amplitude of a4 has been larger than the other modes. Until the damping coefficients increase to 8, there is only the fourth mode left for the parametrically forced resonance with the absence of the first and third modes. The reason of this phenomenon is that the energy transferred from the fourth mode to the other modes is reducing with the damping coefficients increasing.


## 5.3. Influence of Electromagnetic Parameters

Figure5 illustrates the relationship between the amplitude and the magnetic flux-density. The simultaneous resonances can be excited at a certain value. The amplitudes are decreasing with the fundamental component of the magnetic flux-density increasing, which presents the fundamental component as a damping for the resonance of the stator system in Figure 5(a). According to (2), the influences of the stator and rotor windings harmonic components are similar. And Figure 5(b) shows that the stator windings harmonic components make the amplitudes of the first and third modes increase in the resonance region. That means the harmonic magnetic field of the stator and rotor is the main source of electromagnetic excitation.The amplitude-magnetic flux-density response curves.
(a)
(b)


## 5.4. Dynamic Characteristics

In the following figures, diagrams of time history, phase chart, FFT, and Poincaré section are presented for further understanding of the combined internal and principal parametric resonances of the stator system. Based on the parameters in Figures2–5, by assuming εδ1=0.12, εδ2=0.48, μ=8, Bδ=0.75T, and Bν=Bμ=0.1T, Figures 6(a) through 6(d) confirm that the system is period-1 motion according to Figure 4(d). The phase chart is a closed circle, and the energy of the fourth mode is much larger than that of the first and third modes. There is only one nonzero point, and the other two zero points coincide with each other in the Poincaré section. That means the stator system is in stable single-frequency and periodic motion. With the damping coefficients reducing from 8 to 6, there are three points in Figure 7(d) and three frequencies in Figure 7(c). The energy of the fourth mode transfers to the first and third modes, but still higher than the other modes. The stator system is in quasiperiod motion. When the damp coefficients reduce from 6 to 2, the energy of the first mode is the highest in Figure 8. And the phase chart presents reciprocating motion in a bounded domain. With the amplitude of electromagnetic excitation increasing to Bδ=0.5T and Bν=Bμ=0.2T, the energy of all three modes is higher than the former, and the amplitude of the stator system increases. Figures 9(a) through 9(d) confirm the conclusion of Figure 5.Dynamic characteristic withεδ1=0.12, εδ2=0.48, μ=8, Bδ=0.75T, and Bμ=0.1T.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré sectionDynamic characteristic withεδ1=0.12, εδ2=0.48, μ=6, Bδ=0.75T, and Bμ=0.1T.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré sectionDynamic characteristic withεδ1=0.12, εδ2=0.48, μ=2, Bδ=0.75T, and Bμ=0.1T.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré sectionDynamic characteristic withεδ1=0.12, εδ2=0.48, μ=2, Bδ=0.5T, and Bμ=0.2T.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré sectionIt is essential to point out that whenεδ1=0.52, there is an additional relation between the first and fourth natural frequencies as ω-4≈3ω-1. But the coupling coefficient Γ4111=0 which means there is only frequency relationship, and no interaction between these two modes. The expression in (11) is still effective. Figures 10(a) through 10(d) show that the stator system presents the period-2 motion.Dynamic characteristic withεδ1=0.52, εδ2=0.48, and ω-4≈3ω-1.
(a)
Time history
(b)
Phase chart
(c)
FFT
(d)
Poincaré section


## 6. Conclusion

The combined internal and principal parametric resonances on the stator system of an asynchronous machine excited by the varying distributed electromagnetic force are investigated in this paper. The stator system is modeled as a continuum coupled by two cylindrical shells. And the distributed electromagnetic force expression with the interaction between mechanical and electromagnetic behavior is given. The numerical simulation results show the influences of mechanical and electromagnetic parameters on the resonances.When there is a resonance relationship between natural and excitation frequencies, and the corresponding coupling coefficients of modes are nonzero, the forced resonances occur. Simultaneously, in the presence of the resonance relationship of natural frequencies, the combined internal and principal parametric resonances occur and excite three modes. The jump phenomenon is pointed out and the amplitudes of the first and third modes are much higher than the fourth mode. Energy can be transferred directly from the excited fourth mode to the other modes through the internal mechanism. The harmonic magnetic field of the stator and rotor is the main source of electromagnetic excitation. Furthermore, diagrams of time history, phase chart, FFT, and Poincaré section are presented to show the complicated dynamic characteristic of the stator system with the response of period-1, quasiperiod, and particularly double period motions which should be avoided in engineering practice.



---

*Source: 835104-2014-12-22.xml*