Shock and Vibration
Volume 2014 (2014), Article ID 835104, 13 pages
http://dx.doi.org/10.1155/2014/835104

Research Article

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

Baizhou Li1,2 and Qichang Zhang¹

¹School of Mechanical Engineering, Tianjin University, Tianjin 300072, China
²Tianjin Key Laboratory of High Cutting and Precision Machining, Tianjin University of Technology and Education, Tianjin 300222, China

Received 31 March 2014; Accepted 30 July 2014; Published 22 December 2014

Academic Editor: Dumitru I. Caruntu

Copyright © 2014 Baizhou Li and Qichang Zhang. 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.

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.

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–12]. 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. [13] analyzed the natural frequencies and mode shapes by a semianalytical method based on thick double-shell systems. Established a continuum beam model for the nonlinear vibration of rotor-bearing system. Nguyen [14] and Hegazy [15] 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 displacement of the stator core. Meanwhile, the radial displacement changes the thickness of the air-gap and then changes the air-gap permeance as

where is the average air-gap permeance. is the dimensionless radial displacement of the stator core, and 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–12] 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

where

And are the frequencies of electromagnetic force corresponding to . , , and 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.

(a)

(b)

Figure 1: Double-shell stator system.

Based on Hamilton’s Principle, the displacement variational principle for nonlinear vibration of thin elastic shell can be expressed as

where is kinetic energy, and is potential energy of thin shell. , , and are components of electromagnetic force, and , , and are cylindrical coordinates, respectively. The coincident part of the double-shell stator system can be considered as being fixed together by a supposed force . Then, the force conditions of the stator core and the frame can be written as

where is the centroid radius of the stator core. The superscript is for the stator core and is for the frame. The displacement solutions can be written as the products of periodic function and modal shape functions determined by the boundary condition as

where 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 in the coincident region , the vibration equation with modal damping can be obtained as

where is modal damping coefficient, and 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

where are orders of mode shapes distinguishing from , and

are the frequencies of the stator system in working state. are the coupling coefficients which present the interaction of different mode shapes. And and are expressed as

Equations (8a) and (8b) show that the forced excitation term exists in the dynamic equations if and only if . 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 , and the nearness of excitation frequency to one natural frequency can be expressed by introducing the external detuning parameter according to

We introduce small parameter to transform the scale and , and the method of multiple scales [16] is applied to solve (8a) and (8b) as

where is a fast time scale and is a slow time scale. Inserting (12) into (8a) and (8b), and separating at each order of , one obtains

where . At the solution can be written in the following convenient form:

where and are real. are the amplitudes; are the phases. . stands for complex conjugate and .

Substituting (14) into the -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

where stands for the derivation with respect to , , and .

For the steady-state solution, by canceling the time derivatives (), 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

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

where

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

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 . The stationary frequency-response curves are illustrated in Figures 2–4. In these figures, the ordinate 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.

(a)

(b)

(c)

Figure 2: The frequency-response curves. (a) , (b) , and (c) .

(a)

(b)

(c)

Figure 3: The frequency-response curves. (a) , (b) , and (c) .

(a)

(b)

(c)

(d)

Figure 4: The frequency-response curves.

5.1. Influence of Detuning Parameters

Figure 2 shows the stationary frequency-response curves with three cases of internal resonance detuning parameter . All the modal damping coefficients are assumed to be 2. When the internal resonance detuning parameter is zero, the response curves are symmetric and present M-shape. If , the symmetry axis moves to the left based on (17) and inclines towards right. And if , 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 , which originates from the combined forced and parametric resonance, is independent of the internal resonance detuning parameter . Once the amplitude is larger than , the solution becomes unstable. The internal resonance detuning parameter has a direct influence on the energy transfer when the simultaneous resonances are excited.

Figure 3 shows the stationary frequency-response curves with three cases of forced resonance detuning parameter . Similarly, only when , the response curves are symmetric. When , there is only the fourth mode of the stator system left under parametrically forced resonance with high amplitude. If , the right of and curves rises, and all the curves move towards left based on (17). On the contrary, if , the left of and curves rises, and all the curves move towards right. This means that the amplitudes of and increase since the simultaneous resonances occur and then drop to zero when the energy transfer vanishes in Figure 3(b). But the amplitudes of and 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 Figure 4. If , 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 increases slowly, the unstable one decreases, and the resonance amplitudes of and reduce gradually with the increase of damping coefficients. When the damping coefficients increase to 6, the resonance amplitude of 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

Figure 5 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.

(a)

(b)

Figure 5: The amplitude-magnetic flux-density response curves.

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 Figures 2–5, by assuming , , , , and , 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 and , 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.

(a) Time history

(b) Phase chart

(c) FFT

(d) Poincaré section

Figure 6: Dynamic characteristic with , , , , and .

(a) Time history

(b) Phase chart

(c) FFT

(d) Poincaré section

Figure 7: Dynamic characteristic with , , , , and .

(a) Time history

(b) Phase chart

(c) FFT

(d) Poincaré section

Figure 8: Dynamic characteristic with , , , , and .

(a) Time history

(b) Phase chart

(c) FFT

(d) Poincaré section

Figure 9: Dynamic characteristic with , , , , and .

It is essential to point out that when , there is an additional relation between the first and fourth natural frequencies as . But the coupling coefficient 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.

(a) Time history

(b) Phase chart

(c) FFT

(d) Poincaré section

Figure 10: Dynamic characteristic with , , and .

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.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by a grant from the innovation team training plan of Tianjin colleges (TD12-5043).

References

C. Wang and J. C. S. Lai, “Vibration analysis of an induction motor,” Journal of Sound and Vibration, vol. 224, no. 4, pp. 733–756, 1999.

F. Ishibashi, M. Matsushita, S. Noda, and K. Tonoki, “Change of mechanical natural frequencies of induction motor,” IEEE Transactions on Industry Applications, vol. 46, no. 3, pp. 922–927, 2010.

S. P. Verma and A. Balan, “Experimental investigations on the stators of electrical machines in relation to vibration and noise problem,” IEE Proceedings-Electric Power Applications, vol. 145, no. 5, pp. 455–461, 1998.

M. Al Nahlaoui, D. Braunisch, B. Eichinger, S. Kulig, B. Ponick, and U. Werner, “Calculation methods for electromagnetically excited noise in induction motors,” in Proceedings of the 1st International Electric Drives Production Conference (EDPC ’11), pp. 124–131, Nuremberg, Germany, September 2011.

P. Pellerey, V. Lanfranchi, and G. Friedrich, “Coupled numerical simulation between electromagnetic and structural models. influence of the supply harmonics for synchronous machine vibrations,” IEEE Transactions on Magnetics, vol. 48, no. 2, pp. 983–986, 2012.

A. Shahaj and S. D. Garvey, “A possible method for magnetostrictive reduction of vibration in large electrical machines,” IEEE Transactions on Magnetics, vol. 47, no. 2, pp. 374–385, 2011.

A. Belahcen, “Magnetoelastic coupling in rotating electrical machines,” IEEE Transactions on Magnetics, vol. 41, no. 5, pp. 1624–1627, 2005.

D. Guo, F. Chu, and D. Chen, “The unbalanced magnetic pull and its effects on vibration in a three-phase generator with eccentric rotor,” Journal of Sound and Vibration, vol. 254, no. 2, pp. 297–312, 2002.

P. Pennacchi and L. Frosini, “Dynamical behaviour of a three-phase generator due to unbalanced magnetic pull,” IEE Proceedings-Electric Power Applications, vol. 152, no. 6, pp. 1389–1400, 2005.

P. Pennacchi, “Computational model for calculating the dynamical behaviour of generators caused by unbalanced magnetic pull and experimental validation,” Journal of Sound and Vibration, vol. 312, no. 1, pp. 332–353, 2008.

P. Pennacchi, “Nonlinear effects due to electromechanical interaction in generators with smooth poles,” Nonlinear Dynamics, vol. 57, no. 4, pp. 607–622, 2009.

H. Im, H. H. Yoo, and J. Chung, “Dynamic analysis of a BLDC motor with mechanical and electromagnetic interaction due to air gap variation,” Journal of Sound and Vibration, vol. 330, no. 8, pp. 1680–1691, 2011.

G. H. Qing, Y. H. Liu, J. J. Qiu, and X. J. Meng, “A semi-analytical method for the free vibration analysis of thick double-shell systems,” Finite Elements in Analysis and Design, vol. 42, no. 10, pp. 837–845, 2006.

H. Nguyen, “Simultaneous resonances involving two mode shapes of parametrically-excited rectangular plates,” Journal of Sound and Vibration, vol. 332, no. 20, pp. 5103–5114, 2013.

U. H. Hegazy, “Nonlinear vibrations of a thin plate under simultaneous internal and external resonances,” Journal of Vibration and Acoustics, Transactions of the ASME, vol. 132, no. 5, Article ID 051004, 2010.

A. H. Nayfeh and D. T. Mook, Nonlinear Oscillations, John Wiley & Sons, New York, NY, USA, 1979.