3516iama 2 .pdf


Nom original: 3516iama_2.pdfTitre: untitled

Ce document au format PDF 1.3 a été généré par / Acrobat Distiller 8.0.0 (Windows), et a été envoyé sur fichier-pdf.fr le 02/02/2014 à 16:17, depuis l'adresse IP 197.205.x.x. La présente page de téléchargement du fichier a été vue 645 fois.
Taille du document: 1 Mo (4 pages).
Confidentialité: fichier public


Aperçu du document


3516

IEEE TRANSACTIONS ON MAGNETICS, VOL. 46, NO. 8, AUGUST 2010

Study of Static and Dynamic Eccentricities of a
Synchronous Generator Using 3-D FEM
B. A. T. Iamamura1 , Y. Le Menach2 , A. Tounzi2 , N. Sadowski3 , and E. Guillot1
EDF-Lamel 1 Av. Général de Gaulle Clamart 92141 Cedex, France
L2EP-Lamel University of Lille Nord de France, Cité Scientifique Villeneuve d’Ascq 59655, France
GRUCAD/UFSC Florianópolis-SC CP. 476, 88040-900, Brazil
This paper deals with the study of static and dynamic rotor eccentricities of a synchronous generator using 3-D FEM. First, both
eccentricity cases are introduced, as well as the used approach. Then the models of the structure, based on the magnetic scalar potential
formulation and on the magnetic vector potential formulation, are presented. Two eccentricity modeling methods are introduced, the
2-D eccentricity and the 3-D eccentricity. Results obtained for magnetostatic and magnetodynamic cases related to both defects and their
combination are presented and discussed.
Index Terms—Fault diagnosis, finite element methods, magnetic flux measurement, turbogenerators.

I. INTRODUCTION
YNCHRONOUS generators are widely used in high power
conversion plants and their predictive maintenance constitutes a very important task in order to avoid an interruption of
production. Therefore, the forecast and/or the detection of eventual defects are essential.
A very common defect in the rotor of such generators is eccentricity. Many techniques for diagnosis of this type of failure
have been proposed, including analysis of the current in parallel
windings [1], analysis of the rotor and stator vibrations [2], analysis of the shaft voltage, use of capacitive sensors [3] and use of
flux probes in the stator iron stack [4]. A 3-D model of rotor eccentricity with force analysis has already been presented [5], but
only static eccentricity was analyzed. Another approach is based
on the measure of the magnetic field density in the air gap of the
machine [6]. In this paper, we focus on this last method using
numerical modeling to find signatures of eccentricity defects
to detect them. Furthermore, to avoid numerical errors which
can occur using different meshes, we introduce a simple procedure which leads to study the operating of the machine in both
healthy and defected cases using the same mesh.
In this paper, we present the study of the rotor static and
dynamic eccentricities of a synchronous generator using 3-D
finite element analysis. Once both defects and the detection
method presented, we introduce the used 3-D-FEM model
which is based on the scalar and on the magnetic vector
potential formulations. Both 2-D and 3-D eccentricity cases
are investigated. Results for a specific synchronous generator
are given for both defects and their combination. In the case
of 2-D eccentricity, the calculations are achieved using both
potential formulations to compare the results obtained by both
approaches. Then, the effects of the induced currents in the
damping bars on the defect signatures are highlighted.
In the case of 3-D eccentricity, calculations have been carried
out using the scalar potential formulation and the results are
compared to the 2-D case signatures.

S

Manuscript received December 23, 2009; revised January 28, 2010; accepted
February 07, 2010. Current version published July 21, 2010. Corresponding author: B. A. T. Iamamura (e-mail: ba.tanno-iamamura@ed.univ-lille1.fr).
Color versions of one or more of the figures in this paper are available online
at http://ieeexplore.ieee.org.
Digital Object Identifier 10.1109/TMAG.2010.2043347

Fig. 1. Eccentricities: (a) static; (b) dynamic.

Fig. 2. Measurement method.

II. ECCENTRICITY CASES
Generally, we can distinguish two eccentricity cases. In the
static case, the rotor turns around its own rotation axis but this
one is displaced with regard to the stator axis [Fig. 1(a)].
In the case of dynamic eccentricity, the center of the rotor
axis is the same as the center of the stator one, but it does not
correspond to the rotation axis of the rotor [Fig. 1(b)]. These
cases can be found alone or combined. When they are combined,
the center of the rotor is not the same as the center of the stator
and the rotor rotation axis is displaced.
III. MEASUREMENT METHOD
It must be noticed that, in the experimental detection of these
defects, when measuring the magnetic flux density signal on
only one point of the machine air-gap, it is impossible to verify
whether there is a static eccentricity or not, because there is no
variation of the radial flux in only one point. To overcome this
problem, two 90 spatially shifted signals, in the case of a 2 pole
pair machine, are measured. These signals are then added and
treated with the Fast Fourier Transform (FFT) (Fig. 2).

0018-9464/$26.00 © 2010 IEEE

IAMAMURA et al.: STUDY OF STATIC AND DYNAMIC ECCENTRICITIES OF A SYNCHRONOUS GENERATOR USING 3-D FEM

3517

In the case of a healthy machine, the measured magnetic flux
densities are in phase opposition and their sum is normally nil.
On the other hand, in the case of a defect, the sum of both signals
is no more nil. Then, a FFT study of the resultant signal leads
to specific harmonics which could be used as signatures of the
different eccentricity cases.
IV. FINITE ELEMENT METHOD APPROACH
To study the structure, we use 3-D finite element method. In
the magnetostatic case, the computation is carried out with the
magnetic scalar potential formulation

Fig. 3. Eccentricities techniques: (a) healthy machine; (b) static; (c) dynamic.

(1)
the magnetic field due
with the magnetic permeability and
the current density, and the magnetic vector potential formulation
(2)
with the magnetic permeability and the current density.
For the conductive region (damper bars), - formulation is
used to compute the eddy current
(3)
with the electrical conductivity and the electrical scalar potential.
In the 2-D case, formulation needs fewer unknowns than
formulation while it is the opposite in 3-D. Consequently,
to spare computations time, it is more interesting to use the
formulation in 2-D cases and the formulation in 3-D ones.
Furthermore, the magnetic materials are modeled taking into
account the nonlinear
curve. At last, the movement is
implemented through the locked step method.
V. ECCENTRICITY APPROACH
A. General Approach
In order to simulate the machine in the healthy and the deficiently cases with the same mesh, i.e., without re-meshing, we
developed some methods. In our calculation code, the movement is implemented through the locked step method with a slip
surface which must be regular and centered in the origin. Furthermore, this surface shall not be moved or deformed.
In the case of an eccentricity, the rotor and the stator do not deform, only the air gap changes. Thus, to simulate an eccentricity
with the same mesh, only the elements which compose the air
gap can be deformed. So, to take into account all the eccentricity
cases while respecting the constraints due to the movement implementation method, we divide the air gap in two layers of elements [Fig. 3(a)]. The nodes that form this border between the
layers in the air gap will build the slip surface, and shall not be
moved in relation to the geometric center, represented by the
dashed lines [Fig. 3(a)].
Then, for the healthy machine, the stator, rotor and slip surface are centered in the geometric center. To take into account
the static eccentricity, all the nodes outside the slip surface (i.e.,
the nodes of the stator) must be moved [Fig. 3(b)] with the same
distance. In the case of dynamic eccentricity, only the nodes of
the rotor have to be moved while the slip surface remains, always, at its initial position [Fig. 3(c)].

Fig. 4. Three-dimensional eccentricity modeling method: (a) healthy machine
(b) static; (c) dynamic.

To simulate the combined eccentricity (static and dynamic)
the nodes of the rotor and of the stator must be moved, but not
the slip surface.
B. Eccentricity Modeling Methods
When the intensity of the eccentricity is similar on both opposite sides of the machine, a two dimension finite analysis is
sufficient. We call this case 2-D eccentricity. Really, the eccentricities are often due to the fact that one side of the shaft is not
aligned whit the other one. Thus, to study this case, 3-D modeling and analysis is necessary. This case is called 3-D eccentricity.
The studied structure is meshed with prism elements with
only one layer of elements in the 2-D eccentricity case. Thus,
to use the nodes displacement procedure, the nodes in the two
machine ends are moved of the same distance.
To study the 3-D eccentricity, we kept the same cross section
mesh as the one used in the 2-D case but 3 layers of prism elements are built along the machine rotation axis. Then, to take
into account the defect, each slice of nodes is moved with the
adapted value to get closer to a real case of a 3-D eccentricity. In
Fig. 4 we can see the transversal section of the simulated cases.
We can notice that, for all cases, the slip surface remains parallel
in the same position. Only the air-gap is deformed to simulate
the defect cases.
VI. APPLICATION
A. Description of the Machine
The proposed approach is then used to study a turbo-generator of small size, similar to a French nuclear plant generator.

3518

Fig. 5. Cross section view of the turbo-generator mesh.

IEEE TRANSACTIONS ON MAGNETICS, VOL. 46, NO. 8, AUGUST 2010

Fig. 7. Two-dimensional eccentricity magnetostatic case, A formulation.

Fig. 6. Magnetic field distribution in the machine.

This machine, of about 26.7 kVA at 50 Hz, has four poles, 48
slots in the stator and 36 slots in the rotor. Besides, the rotor has
36 short-circuited damping bars and the constant air-gap thickness is 1.5 mm.
The same mesh is used for all the eccentricity cases. One
layer of elements is used for 2-D eccentricities (36226 nodes
and 36128 elements) while three layers of elements are built for
3-D eccentricities (72452 nodes and 108384 elements). A cross
section view of the mesh is given in Fig. 5 with a zoom of the
air gap region.
Fig. 6 shows the magnetic field distribution, at no load, in a
cross section of the healthy machine with an excitation current
of 10 A. We can clearly see the symmetry of the 4 poles.
B. Results
Two-Dimensional Eccentricities: Calculations were carried
out for 2-D eccentricity cases, at no load and an excitation current
of 10 A, with a constant magnitude of 0.3 mm along the
rotation axis.
Figs. 7 and 8 show the results obtained with the magnetostatic
formulation (18017 unknowns). In the case of the healthy machine, as expected, the sum of the two magnetic flux density signals is almost nil, and then harmonic free. On the other hand, the
static eccentricity induces a significant 50 Hz-harmonic while
dynamic eccentricity yields harmonics of 25 Hz and 75 Hz [6].
For the static-dynamic eccentricity case, there is a superposition
of both eccentricity signals.
Figs. 9 and 10 show the results obtained using the magnetostatic formulation (36226 unknowns). We can notice that they
are very close to the ones shown above.
In the case of eccentricities, the air-gap width is no more constant all around the rotor. Then, even at no load, the variation

Fig. 8. FFT of two-dimensional eccentricity magnetostatic case,
tion.

A formula-

Fig. 9. Two-dimensional eccentricity magnetostatic case,
formulation.

Fig. 10. FFT of two-dimensional eccentricity magnetostatic case,
formulation.

of the air-gap reluctance induces currents in the damping bars.
Calculations were thus carried out using the magnetodynamic

IAMAMURA et al.: STUDY OF STATIC AND DYNAMIC ECCENTRICITIES OF A SYNCHRONOUS GENERATOR USING 3-D FEM

Fig. 11. Two-dimensional eccentricity magnetodynamic case,
tion.

A-' formula-

Fig. 14. FFT of three-dimensional eccentricity magnetostatic case,
lation.

3519


formu-

The magnitude of the signals is about twice lesser than in the
2-D eccentricity case. However, the same magnetostatic characteristics are found and then, eccentricity may be detected even
on 3-D case.
VII. CONCLUSION

Fig. 12. FFT of two-dimensional eccentricity magnetodynamic case, A -' formulation.

In this paper, 3-D-FEM is used to study static and dynamic
eccentricities of a synchronous generator at no load. Two 90
shifted magnetic flux density signals in the air gap, may lead to
the detection of eccentricities and to differentiate their type by
studying their harmonic content.
A modeling technique for 3-D rotor eccentricity has been presented. This method does not require remeshing for fault case
simulations, which reduces simulation time.
For all studied cases (magnetostatic and , and magnetodynamic - formulations) it is possible to detect eccentricity
and specify whether it is static, dynamic or static-dynamic eccentricity. Faults may be detected even when there are eddy currents in damper bars.
ACKNOWLEDGMENT
This work has been achieved within the framework of
MEDEE2, supported by EDF and the Nord-Pas-de Calais
Council.

Fig. 13. Three-dimensional eccentricity magnetostatic case,
formulation.

- formulation (18017 unknowns). For the same conditions,
Figs. 11 and 12 show the results of 2-D eccentricity.
We can notice that, for the static eccentricity, the damping
currents induce a decrease of the magnetic flux density
50 Hz-harmonic with respect to the magnetostatic case signature while harmonics at 25 and 75 Hz appear. Nevertheless, it
is still possible to differentiate the eccentricity cases since only
static eccentricity leads to a 50 Hz harmonic.
Three-Dimensional Eccentricities: To simulate 3-D eccentricity, the nodes on one extremity are displaced of 0.3 mm. Calculations have then been performed using the magnetostatic
formulation (72452 unknowns). Figs. 13 and 14 show the obtained results. The magnetic flux densities have been calculated
in an air gap center element along the z-axis (i.e., when the nodes
have been moved of about 0.1 mm, 0.2 mm and 0.3 mm).

REFERENCES
[1] A. Foggia, J. Torlay, C. Corewinder, A. Audoli, and J. Hérigault, “Circulating current analysis in the parallel-connected windings of synchronous generators under abnormal operating conditions,” in Proc.
IEMD Int. Conf. Electric Machines and Drives, Seattle, WA, 1999, pp.
634–636.
[2] D. Canha, W. Cronje, A. Meyer, and S. J. Hoffe, “Methods for diagnosing static eccentricity in a synchronous 2 pole generator,” in Proc.
IEEE Power Technol., Lausanne, Switzerland, 2007, pp. 2162–2167.
[3] G. B. Pollock and J. F. Lyles, “Vertical hydraulic generators experience
with dynamic air gap monitoring,” IEEE Trans. Energy Conversion,
vol. 7, no. 4, pp. 660–668, 1992.
[4] J. J. Simond, M. Tu Xuan, and R. Wetter, “An innovative inductive
air-gap monitoring for large low-speed hydro-generators,” presented at
the ICEM Int. Conf. Electrical Machines, Vilamoura, 2008.
[5] Y. Kawase, N. Mimura, and K. Ida, “3-D electromagnetic force analysis of effects of off-center of rotor in interior permanent magnet synchronous motor,” IEEE Trans. Magn., vol. 36, no. 4, pp. 1858–1862,
Apr. 2000.
[6] R. L. Stoll and A. Hennache, “Method of detecting and modelling presence of shorted turns in DC field windings of cylindrical rotor synchronous machines using two airgaps search coils,” IEE Proc., Electric
Power Appl., vol. 135, pp. 281–294, 1988.


3516iama_2.pdf - page 1/4


3516iama_2.pdf - page 2/4


3516iama_2.pdf - page 3/4

3516iama_2.pdf - page 4/4


Télécharger le fichier (PDF)


3516iama_2.pdf (PDF, 1 Mo)

Télécharger
Formats alternatifs: ZIP



Documents similaires


3516iama 2
ijmet 04 03 021
04277904
chen2015
physreva 84 023818
tutorial magnetostatics