Fast Calculation Model and Theoretical Analysis of Rotor Unbalanced Magnetic Pull for Inter-Turn Short Circuit of Field Windings of Non-Salient Pole Generators

Inter-turn short circuit of field windings (ISCFW) may cause the field current of a generator to increase, output reactive power to decrease, and unit vibration to intensify, seriously affecting its safe and stable operation. Full integration of mechanical and electrical characteristics can improve the sensitivity of online monitoring, and detect the early embryonic period fault of small turns. This paper studies the calculations and variations of unbalanced magnetic pull (UMP), of which the excitation source of rotor vibration is the basis and key to online fault monitoring. In grid load operation, ISCFW are first calculated with the multi-loop method, so as to obtain the numerical solutions of the stator and the rotor currents during the fault. Next, the air-gap magnetic field of the ISCFW is analyzed according to the actual composition modes of the motor loops in the fault, so as to obtain the analytic expressions of the air-gap magnetic motive force (MMF) and magnetic density. The UMP of the rotor is obtained by solving the integral of the Maxwell stress. The correctness of the electric quantity calculation is verified by the ISCFW experiment, conducted in a one pair-pole non-salient pole model machine. On this basis, comparing the simulation analysis with the calculation results of the model in this paper not only verifies the accuracy of the electromagnetic force calculation, but also proves that the latter has the advantages of a short time consumption and high efficiency. Finally, the influencing factors and variation law of UMP are analyzed by means of an analytic model. This develops a base for the online monitoring of ISCFW with the integration of mechanical and electrical information.


Introduction
The requirements for the safe and stable operation of synchronous generators continually increase with the continuous development of power systems.Inter-turn short circuit of field windings (ISCFW) is a common electrical fault in generators.A small turns short circuit fault may not affect the normal operation of generators [1,2].A spreading fault will cause the field current of the generator to increase, output reactive power to reduce, and unit vibration to intensify.This causes burning of shaft necks and bearings, seriously threatening the safe and stable operation of the generator unit and the power system [3,4].
The traditional method of monitoring the ISCFW is to use the electrical quantity in the fault to achieve fault monitoring [5].Although most of these methods are mature, there are still many problems, such as the need to rely on branch current transformers (the vast majority of turbines cannot Energies 2017, 10, 732 2 of 19 be installed), lower sensitivity, and larger monitoring dead zones.As a result, domestic and foreign scholars are looking for new monitoring methods.
In the motor, owing to mechanical and electrical coupling, the ISCFW will cause partial magnetic poles to short circuit, thereby causing air gap magnetic field distortion.This produces the air gap electromagnetic force wave and unbalanced magnetic pull (UMP) differing from normal operation, causing the radial vibration of the motor rotors.Many domestic and foreign scholars have tried to use the vibration characteristics of the generator to monitor the fault.Trutt proved that unit bearing vibration and shaft vibration would exceed the normal standard when the ISCFW occurred [6].The increase of vibration amplitude and field current, and the decrease of output reactive power were tending toward uniformity [7,8].The calculation model for fault electromagnetic torque was established based on the multi-loop and finite element methods.The correctness of the model was verified experimentally, and the harmonic characteristics of the steady-state electromagnetic torque in the fault were analyzed in the literature [4,9].The electromagnetic forces on the motor eccentric fault were calculated using a numerical method from the literature [10,11], and the obtained results confirm the experimental results.
A great deal of research has been achieved in the literature [12,13].The vibration equations of stator and rotor were derived, and the vibration characteristics of the stator and rotor after faults were determined [14,15].From this, we can see that it is feasible to use the vibration characteristics of the generator rotor as an auxiliary fault criterion.However, most of the existing studies focused on the vibration response of the rotor in the fault, and studies on vibration excitation sources remain insufficient and are not detailed.UMP is the excitation source of rotor vibration in the fault, and its numerical size is closely related to the severity of the fault and the running state of the generator.The accurate calculation of UMP in a fault, the variation law, and the influencing factors are the basis and key to detecting ISCFW using rotor vibration characteristics.
Aiming at the common one-pair pole turbo generator, this paper puts forward a method to obtain the UMP on a rotor by analytic calculation of the air gap magnetic field in ISCFW.First of all, the ISCFW of the turbo generator in the grid connection operation are calculated using the multi-loop method, and the numerical solutions of the stator and the rotor currents on the faults are obtained.A fault experiment on the model machine is conducted, and the experimental results are compared with theoretical calculation results.On this basis, the analytic calculations of air gap magnetic motive force (MMF) and magnetic density on inter-turn short circuit faults are conducted; the size and direction of the UMP on the rotor in faults are calculated using the Maxwell stress tensor method, and the analytic calculation models of the UMP are obtained.Then, in order to solve the problem that UMP cannot be measured directly by experiment, finite element calculations of UMP are conducted by means of the finite element physical generator model established with ANSYS simulation software.The correctness of the rapid calculation model is verified by comparing the finite element calculation and analytic model calculation results.Finally, the influences of the short position, number of turns, and the running state on the rotor UMP in the generator rotor inter-turn short circuit faults are studied.The variation law of UMP is obtained by means of the rapid calculation model.The conclusion not only provides the theoretical basis for the further study of rotor vibration response characteristics, but also lays the foundation for the online monitoring of rotor inter-turn short circuits with mechanical and electrical information integration.

Fast Calculation Model of UMP
The accurate calculation of the UMP in the fault is the basis of the study on the rotor vibration excitation source.Because the UMP of the rotor cannot be directly measured experimentally, obtaining of the UMP is achieved by finite element simulation.Although this method can fully consider the saturation of the magnetic field, the internal magnetic field distribution of the motor is set closer to the actual situation so that the final simulation results error is small.The biggest drawback of this method is that its calculation process is too time-consuming.For a real generator, the higher the Energies 2017, 10, 732 3 of 19 number of winding coils, the slower the calculation will be.This is difficult to accept in engineering.The rapid calculation model of UMP based on the multi-loop method proposed in this paper can take into account both accuracy and calculation speed.

Loop Currents in Rotor Inter-Turn Short Circuit Faults
The basic thinking of calculating ISCFW with the multi-loop method is as follows: First, the inductive parameters of coils are obtained by taking a single coil as a research object.Next, in accordance with the actual winding connection and composition mode in the faults, the winding loop parameters are calculated through the related coil parameters [16,17].This method involves the stator phase and rotor winding, and calculates current and voltage of all branches of the stator and field windings, along with the damp circuit [18].
Under the premise of field winding, the circuit topology structure changes.Unbalanced current inside the stator phase winding from the fault is considered, and the equation is stated with all loop currents (with the stator and the rotor as variables) and is obtained based on the multi-loop theory [19]: where M and R are a loop inductance matrix and a loop resistance matrix, respectively.The value of M as a time-varying matrix varies with the relative motion between the stator and the rotor.The values of M T and R T as constant matrices are related to the grid line parameters, the transformer resistance, leakage inductance, the internal resistance, and the internal inductance of the field system.Column vector I is the current in each circuit of the stator and rotor, including the normal circuit and the fault additional circuit of the rotor field winding, the damp circuit, and the stator armature winding circuit.
Column vector E, as a known quantity, is composed of the grid voltage and the power supply voltage of the field system.D is a differential operator, d/dt.Equation ( 1) is the multi-loop mathematical model of the synchronous generator on field winding inter-turn short circuit faults.The numerical solutions of the current and the voltage of the stator and the rotor on the faults can be obtained by means of an appropriate numerical solving method, so as to realize the transient simulation of the faults.The analyses and calculations below are based on the results of this model.

Air Gap MMF and Magnetic Density in Rotor
Inter-Turn Short Circuit Fault

Air Gap MMF Generated by Rotor Field Current
The field winding of each pole of the non-salient pole synchronous generator is composed of distributed concentric coils with different pitches in series.Let the starting point of the rotor spatial coordinate θ be the d-axis of the rotor.The MMF alone produced by each concentric coil is a rectangular wave symmetrical with the d-axis center line, and distributed around the whole motor air gap, as shown in Figure 1.
Energies 2017, 10, 732 3 of 19 in engineering.The rapid calculation model of UMP based on the multi-loop method proposed in this paper can take into account both accuracy and calculation speed.

Loop Currents in Rotor Inter-Turn Short Circuit Faults
The basic thinking of calculating ISCFW with the multi-loop method is as follows: First, the inductive parameters of coils are obtained by taking a single coil as a research object.Next, in accordance with the actual winding connection and composition mode in the faults, the winding loop parameters are calculated through the related coil parameters [16,17].This method involves the stator phase and rotor winding, and calculates current and voltage of all branches of the stator and field windings, along with the damp circuit [18].
Under the premise of field winding, the circuit topology structure changes.Unbalanced current inside the stator phase winding from the fault is considered, and the equation is stated with all loop currents (with the stator and the rotor as variables) and is obtained based on the multi-loop theory [19]: where ' M and ' R are a loop inductance matrix and a loop resistance matrix, respectively.The value of ' M as a time-varying matrix varies with the relative motion between the stator and the rotor.The values of T M and T R as constant matrices are related to the grid line parameters, the transformer resistance, leakage inductance, the internal resistance, and the internal inductance of the field system.Column vector '  I is the current in each circuit of the stator and rotor, including the normal circuit and the fault additional circuit of the rotor field winding, the damp circuit, and the stator armature winding circuit.Column vector E, as a known quantity, is composed of the grid voltage and the power supply voltage of the field system.D is a differential operator, d / dt .Equation (1) is the multi- loop mathematical model of the synchronous generator on field winding inter-turn short circuit faults.The numerical solutions of the current and the voltage of the stator and the rotor on the faults can be obtained by means of an appropriate numerical solving method, so as to realize the transient simulation of the faults.The analyses and calculations below are based on the results of this model.The expression of the air gap MMF F fd,h generated by the hth (h = 1, 2, 3, . . . ) concentric coil in the field winding is:

Air Gap MMF and Magnetic Density in Rotor Inter-Turn Short Circuit Fault
where y h and w fd,h are, respectively, the pitch and turn number of the hth concentric coil.α h,1 and α h,2 are the rotor coordinate electric angles corresponding to its two edges.On field winding inter-turn short circuits with a synchronous generator, the field winding is divided into a normal field circuit (with current i f ) and the field fault additional circuit (with current i fkL ), as shown in Figure 2.
Energies 2017, 10, 732 4 of 19 The expression of the air gap MMF fd,h F generated by the hth (h = 1, 2, 3, …) concentric coil in the field winding is: where h y and fd,h w are, respectively, the pitch and turn number of the hth concentric coil.
,1 h  and ,2 h  are the rotor coordinate electric angles corresponding to its two edges.
On field winding inter-turn short circuits with a synchronous generator, the field winding is divided into a normal field circuit (with current if) and the field fault additional circuit (with current ifkL), as shown in Figure 2.
Loops of field windings with inter-turn short circuit.
In Figure 2, EZF is electromotive force, RfkL is the short-circuit transition resistance, and RZF is internal resistance of the field power supply.In the rotor inter-turn short circuit faults, in addition to the DC component, there is an AC component in the field current.The AC component is small, so its influences on the air gap magnetic field are not considered here.MMF where f fd,n,h F and f fd,f,h F are, respectively, the MMF generated by the normal field circuit and the fault additional circuit in the hfth fault coil.

Air Gap MMF Generated by the Rotor Damp Current
The damp winding of a synchronous motor is generally a cage circuit composed of damping bars [20], as shown in Figure 3.

N S S
Rotation direction of the rotor ...   In Figure 2, E ZF is electromotive force, R fkL is the short-circuit transition resistance, and R ZF is internal resistance of the field power supply.In the rotor inter-turn short circuit faults, in addition to the DC component, there is an AC component in the field current.The AC component is small, so its influences on the air gap magnetic field are not considered here.MMF F fd,h f in the air gap-generated by the fault winding-is equivalent to the superposition of the MMF generated by the normal field circuit and reverse MMF generated by the fault additional circuit, as shown in Equation (3): where F fd,n,h f and F fd,f,h f are, respectively, the MMF generated by the normal field circuit and the fault additional circuit in the h f th fault coil.

Air Gap MMF Generated by the Rotor Damp Current
The damp winding of a synchronous motor is generally a cage circuit composed of damping bars [20], as shown in Figure 3.
Energies 2017, 10, 732 4 of 19 The expression of the air gap MMF fd,h F generated by the hth (h = 1, 2, 3, …) concentric coil in the field winding is: where h y and fd,h w are, respectively, the pitch and turn number of the hth concentric coil.
,1 h  and ,2 h  are the rotor coordinate electric angles corresponding to its two edges.
On field winding inter-turn short circuits with a synchronous generator, the field winding is divided into a normal field circuit (with current if) and the field fault additional circuit (with current ifkL), as shown in Figure 2.
Loops of field windings with inter-turn short circuit.
In Figure 2, EZF is electromotive force, RfkL is the short-circuit transition resistance, and RZF is internal resistance of the field power supply.In the rotor inter-turn short circuit faults, in addition to the DC component, there is an AC component in the field current.The AC component is small, so its influences on the air gap magnetic field are not considered here.MMF where F are, respectively, the MMF generated by the normal field circuit and the fault additional circuit in the hfth fault coil.

Air Gap MMF Generated by the Rotor Damp Current
The damp winding of a synchronous motor is generally a cage circuit composed of damping bars [20], as shown in Figure 3.

N S S
Rotation direction of the rotor ...   Unlike field winding, there is only an AC harmonic component in the current of the damp winding without a DC component.The expression for the µth harmonic current of the jth (j = 1, 2, . . . ) damp circuit: where ω is the synchronous angular frequency, t is time, I d,j,µ and ϕ d,j,µ are, respectively, the root mean square (RMS) and the phase of the µth harmonic current (µ is the round number for common one-pair pole generators) of the jth damp circuit.The MMF generated in the air gap is: where y d,j is the pitch of the jth damp circuit, α d,j,1 and α d,j,2 are the electric angles of the rotor coordinate corresponding to its two edges, and w d is the number of turns of the damp circuit coil.The MMF in the air gap generated by all damp circuits of the whole rotor can then be obtained: When a synchronous generator is running normally, the current in the stator branch is a fundamental wave.After the ISCFW, the stator winding is still normal, but the space harmonic magnetic field generated by the fault field winding would induce an even-number harmonic unbalanced current inside the stator winding [21].
The zero point of the stator spatial coordinates x is placed in the center line of stator slot number one.A counterclockwise direction is taken as the forward direction of x, as shown in Figure 4.The expression for the µth harmonic current of the jth coil (j = 1, 2, . . . ) of a stator armature A phase ai (i = 1, 2, ...) branch is: i s,ai,j,µ = √ 2I s,ai,j,µ cos(µωt + ϕ s,ai,j,µ ) Energies 2017, 10, 732 5 of 19 Unlike field winding, there is only an AC harmonic component in the current of the damp winding without a DC component.The expression for the  th harmonic current of the jth (j = 1, 2, …) where  is the synchronous angular frequency, t is time, d, , j I  and d, , j   are, respectively, the root mean square (RMS) and the phase of the  th harmonic current (  is the round number for common one-pair pole generators) of the jth damp circuit.The MMF generated in the air gap is: ,0 or 2 2 where d,j y is the pitch of the jth damp circuit, are the electric angles of the rotor coordinate corresponding to its two edges, and d w is the number of turns of the damp circuit coil.The MMF in the air gap generated by all damp circuits of the whole rotor can then be obtained: When a synchronous generator is running normally, the current in the stator branch is a fundamental wave.After the ISCFW, the stator winding is still normal, but the space harmonic magnetic field generated by the fault field winding would induce an even-number harmonic unbalanced current inside the stator winding [21].
The zero point of the stator spatial coordinates x is placed in the center line of stator slot number one.A counterclockwise direction is taken as the forward direction of x, as shown in Figure 4.
The expression for the  th harmonic current of the jth coil (j = 1, 2, …) of a stator armature A phase Its MMF generated within the air gap is: Its MMF generated within the air gap is: 2π w s i s,ai,j,µ , α s,ai,j,1 ≤ x < α s,ai,j,2 − y s 2π w s i s,ai,j,µ , 0 ≤ x < α s,ai,j,1 or α s,ai,j,2 ≤ x < 2π (8) Energies 2017, 10, 732 6 of 19 where y s is the coil pitch of the stator winding, w s is the number of coil turns, and α s,ai,j,1 and α s,ai,j,2 are the stator coordinate electric angles corresponding to the two edges of the coil.The MMF generated by the stator A phase current and the MMF generated by ABC three phase currents can be obtained:

Air Gap Magnetic Density in Fault
In order to facilitate this calculation, the calculation results of 2.2.1-2.2.3 are converted to the stator coordinate system.Let γ 0 represent the electric angle at t = 0 of the rotor d-axis.The stator coordinate axis, γ is the rotor position angle (i.e., the electric angle of the rotor d-axis counterclockwise rotation leading over the stator coordinate axis): where the rotor rotates at a constant speed, γ = ωt + γ 0 .Replacing θ with x − γ obtains the expression of the air gap MMF in the fault: For a uniform magnetic circuit, the formula for air gap magnetic density [22] can be obtained according to Equation (12): where λ = µ 0 /l represents the permanent per unit area, µ 0 is air gap magnetic permeability, and l is equivalent air gap length.

Analysis of Rotor UMP on Inter-Turn Short Circuit Fault
When the ISCFW occurs, the magnetic field in the air gap is no longer symmetrical, and the rotor would be subjected to UMP.The rotor UMP can be divided into radial magnetic pull and tangential magnetic pull.However, the value of radial magnetic pull is much higher than that of tangential magnetic pull.Only the radial magnetic pull of the rotor is studied, and the tangential magnetic pull is no longer considered [23,24].Based on the distribution of air gap magnetic density on faults calculated in Section 2.2.3, electromagnetic force density (N/m 2 )-radially distributed at an arbitrary position along the circumference of the generator-can be obtained by means of the Maxwell stress tensor method: In order to calculate the UMP on the rotor, the above equation shall be integrated along a closed path on the rotor surface: here the numerical method is used to solve the above integral.Discretization of continuous data is required.N δ units may be evenly taken on this path.An electric angle of the stator coordinates corresponding to each discrete unit is x m (1 ≤ m ≤ N δ ).The radial force withstood by each unit can be obtained by calculation: where D rotor is the rotor diameter, L rotor is the length of rotor core, and ∆x is the difference of the stator coordinate electric angles corresponding to the two adjacent discrete units.The X-axial and Y-axial components can be obtained in the global coordinate system.Their superposition will be the X-axis component F X and Y-axis component F Y of the rotor UMP.
On the basis of the obtained F X and F Y , the modulus value F UMP and the direction α UMP of the UMP (α UMP is the angle between the UMP and X-axis of the global coordinate system) can be obtained:

Experiment and Simulation Verification
The UMP of the rotor cannot be measured in the experiment.Verification in this paper for a fast calculation model of the UMP is as follows.First, the correctness of the model electric quantity calculation is verified by comparing the experimental results with the simulation results of the stator and rotor currents on faults.The UMP of the rotor is calculated by means of the established generator finite element physical model, with the calculated stator and rotor currents.Then, the correctness of the fast calculation model can be indirectly verified by comparing the results of finite element analysis with the calculation results of the fast model proposed by this paper.

Experiment on ISCFW
It is important to verify the accuracy of the multi-loop mathematical model.This is accomplished by calculating a one-pair pole A1552 model machine grid loaded on the inter-turn short circuit faults.A generator rotor inter-turn short circuit experiment by means of a machine model is necessary.Figure 5 is the experimental wiring diagram.The A1552 model machine is driven by a DC motor which is controlled by a SIMENS 6RA70 digital speed regulating system, with a rated speed of 3000 RPM.The machine model connects to the power grid through a three-phase voltage regulator, and the field winding is supplied by a DC stabilized voltage source with minimal ripple.The parameters of the machine model A1552 and the structure of the rotor are shown in Table 1 and Figure 6, respectively.The rotor is a laminated iron core with two poles, comprising 24 uniformly-spaced slots, of which the number of the actual field slots is 20.There are two dead slots under each pole.Each pole has five field concentric coils, placed symmetrically in 10 field slots.The total turn number of field windings in series is 660.Seven taps are led from inside of the field winding The parameters of the machine model A1552 and the structure of the rotor are shown in Table 1 and Figure 6, respectively.The rotor is a laminated iron core with two poles, comprising 24 uniformly-spaced slots, of which the number of the actual field slots is 20.There are two dead slots under each pole.Each pole has five field concentric coils, placed symmetrically in 10 field slots.The total turn number of field windings in series is 660.Seven taps are led from inside of the field winding of the machine model.In addition to two taps of the initial and terminal ends, the other five taps are used for the field winding inter-turn short circuit experiment.The parameters of the machine model A1552 and the structure of the rotor are shown in Table 1 and Figure 6, respectively.The rotor is a laminated iron core with two poles, comprising 24 uniformly-spaced slots, of which the number of the actual field slots is 20.There are two dead slots under each pole.Each pole has five field concentric coils, placed symmetrically in 10 field slots.The total turn number of field windings in series is 660.Seven taps are led from inside of the field winding of the machine model.In addition to two taps of the initial and terminal ends, the other five taps are used for the field winding inter-turn short circuit experiment.An experimental comparison with such a large number of turns is designed to eliminate deviation in the simulation and calculation errors caused by motor manufacturing factors.The experimental current waveforms of the a2 branch of the grid-connected loaded machine model agree well with the simulation waveforms using the multi-loop method, as shown in Figure 7.An experimental comparison with such a large number of turns is designed to eliminate deviation in the simulation and calculation errors caused by motor manufacturing factors.The experimental current waveforms of the a2 branch of the grid-connected loaded machine model agree well with the simulation waveforms using the multi-loop method, as shown in Figure 7.In order to further verify the accuracy of the calculation with the multi-loop method, the RMS of each harmonic component of the stator and rotor steady-state currents is obtained.Fast Fourier transform (FFT) analysis is compared with the calculated steady state results, and each harmonic component is calculated with the multi-loop method, as shown in Table 2:  In order to further verify the accuracy of the calculation with the multi-loop method, the RMS of each harmonic component of the stator and rotor steady-state currents is obtained.Fast Fourier transform (FFT) analysis is compared with the calculated steady state results, and each harmonic component is calculated with the multi-loop method, as shown in Table 2 The experimental and simulation results with the multi-loop method are in good agreement.The mathematical model is verified as correct for calculating a grid-connected one-pair pole synchronous generator on field winding inter-turn short circuit faults.On this basis, we continue to use the multi-loop method to calculate the RMS.The phase of the rotor damp circuit and stator armature steady-state current is also calculated using Equation (1), as shown in Tables 3 and 4:

Finite Element Model of a Generator
In this paper, the finite element analysis ANSYS software is used to establish the two-dimensional finite element physical model of the entire cross-section of machine model A1552, as shown in Figure 8.

Finite Element Model of a Generator
In this paper, the finite element analysis ANSYS software is used to establish the twodimensional finite element physical model of the entire cross-section of machine model A1552, as shown in Figure 8.The intelligent partition for the model of triangle mesh with six nodes is made.Magnetic field distribution in the air gap is the key to the coil coupling parameters analysis.In order to ensure partition density, the starting position of the partition is set at the air gap.The model after the partition is shown in Figure 9.In each unit of the partition, the kept material magnetic permeability is considered to be constant.The intelligent partition for the model of triangle mesh with six nodes is made.Magnetic field distribution in the air gap is the key to the coil coupling parameters analysis.In order to ensure partition density, the starting position of the partition is set at the air gap.The model after the partition is shown in Figure 9.In each unit of the partition, the kept material magnetic permeability is considered to be constant.

Finite Element Model of a Generator
In this paper, the finite element analysis ANSYS software is used to establish the twodimensional finite element physical model of the entire cross-section of machine model A1552, as shown in Figure 8.The intelligent partition for the model of triangle mesh with six nodes is made.Magnetic field distribution in the air gap is the key to the coil coupling parameters analysis.In order to ensure partition density, the starting position of the partition is set at the air gap.The model after the partition is shown in Figure 9.In each unit of the partition, the kept material magnetic permeability is considered to be constant.

UMP Analysis Procedure with the Finite Element Method
The calculation of the UMP under a short circuit fault with the magnetic field method of finite element analysis was as follows.First, the rotor field current, damp current, stator armature current, and other related running data (using an A1552 model machine grid-connected and loaded on a rotor inter-turn short circuit fault) are calculated with the multi-loop method.The two-dimensional finite element physical model of a generator is built based on the actual size.The running data obtained in the first step are incorporated into the finite element physical model established in the second step.Next, the generator internal electromagnetic field on inter-turn short circuit faults is calculated.The magnetic density of each unit along the air gap circumference is obtained and kept.The electromagnetic force of each unit along the rotor surface, and the UMP on the whole rotor in fault are subsequently calculated.
Steps three and four are repeated for the changed rotor position until the rotor rotates one revolution.In addition, an air gap line set at the center position of the air gap by the model divides the air gap into stator and rotor air gaps.The automatic rotation of the rotor is realized by the nodes on a dividing line between the stator and rotor air gap composing different sets of coupling points.

Comparative Analysis of Analytical Model Calculation and Finite Element Calculation Results
After establishing mathematical and physical models of a generator, the air gap magnetic density distribution and the UMP on the rotor in ISCFW are obtained using the analytical model and finite element calculations.The comparison of the obtained results is as follows: The air gap magnetic density distributions of a synchronous generator in ISCFW obtained with the analytical and finite element models are shown in Figure 10.Here the rotor position angle γ = 0.The X-axis and Y-axis components of the UMP on the rotor are shown in Figures 11 and 12. Figures 13  and 14 are the direction and the modulus value of the UMP.Results show that the UMP and the rotor synchronously rotate.The amplitude periodically changes six times per rotor revolution.As seen in Figures 11-14, the results obtained with the two methods are in good agreement, verifying the correctness of the analytical calculation model of the UMP.In addition, the finite element calculation took up to two hours, while the entire analytic model calculation took only one minute.Compared to the finite element method, the analytical model calculation has the advantages of a shorter time consumption and a higher efficiency.
the first step are incorporated into the finite element physical model established in the second step.Next, the generator internal electromagnetic field on inter-turn short circuit faults is calculated.The magnetic density of each unit along the air gap circumference is obtained and kept.The electromagnetic force of each unit along the rotor surface, and the UMP on the whole rotor in fault are subsequently calculated.Steps three and four are repeated for the changed rotor position until the rotor rotates one revolution.
In addition, an air gap line set at the center position of the air gap by the model divides the air gap into stator and rotor air gaps.The automatic rotation of the rotor is realized by the nodes on a dividing line between the stator and rotor air gap composing different sets of coupling points.

Comparative Analysis of Analytical Model Calculation and Finite Element Calculation Results
After establishing mathematical and physical models of a generator, the air gap magnetic density distribution and the UMP on the rotor in ISCFW are obtained using the analytical model and finite element calculations.The comparison of the obtained results is as follows: The air gap magnetic density distributions of a synchronous generator in ISCFW obtained with the analytical and finite element models are shown in Figure 10.Here the rotor position angle 0   .The X-axis and Y-axis components of the UMP on the rotor are shown in Figures 11 and 12. Figures 13 and 14 are the direction and the modulus value of the UMP.Results show that the UMP and the rotor synchronously rotate.The amplitude periodically changes six times per rotor revolution.As seen in Figures 11-14, the results obtained with the two methods are in good agreement, verifying the correctness of the analytical calculation model of the UMP.In addition, the finite element calculation took up to two hours, while the entire analytic model calculation took only one minute.Compared to the finite element method, the analytical model calculation has the advantages of a shorter time consumption and a higher efficiency.

Effects of Field Winding Short Circuit Position on UMP
Post-fault MMF is the sum of the MMF produced by the normal field circuit and the MMF produced by the fault additional circuit.When the additional harmonic current of field and damp winding (caused by post-fault stator armature reaction MMF) is not considered, the MMF properties produced by the normal field circuit are the same as those before the faults.Only the effects of the MMF produced by the fault field circuit on the UMP must be analyzed.Fourier decomposition for MMF produced by the fault field concentric coil is carried out: The rectangular wave MMF is symmetrical around the rotor d-axis, so there is only the cosine term in Equation ( 20).This calculation is factored into: where I f is the DC component of field current, w is the number of short circuit turns, and k y,v is pitch shortening factor of the fault field concentric coils to the vth harmonic.Its computational formula is as follows: where β is pitch shortening ratio of the coil.
From Equation ( 21) we can see that the vth harmonic MMF produced by the fault coil is proportional to the pitch shortening factor k y,v of the coil with regard to the vth harmonic for definite field current and short turns.Therefore, the RMS of the vth harmonic current induced in the stator branch and the amplitude of the space harmonic MMF induced in the air gap are both proportional to k y,v .Each harmonic of fault UMP is generated by magnetic field interactions in the air gap.The amplitude of each harmonic is only affected by k y,v .In the field test, the fundamental component of UMP is measured and observed.This paper takes the fundamental component of UMP as the main research object.UMP thereafter refers to the fundamental component of the UMP.
In order to study the influence of short circuit position on UMP, 11-55 foot field concentric coils on 30 turn short circuit faults are respectively calculated.Conditions of operation are grid-connected.The simulation results are shown in Table 5.When the two field concentric coils with the same pitch shortening factor k y,2 respectively undergo faults, the values of UMP are similar.The UMP with larger pitch shortening factor k y,1 is slightly less than that with smaller k y,1 .
When the field concentric coils with different pitches undergo short circuit faults, UMP is generated by an interaction between the fundamental MMF and the second MMF in the air gap.
The fundamental MMF includes the fundamental field MMF and the MMF produced by the stator harmonic current.The fundamental field MMF can be regarded as the sum of the fundamental MMF produced by the normal field circuit and the fundamental MMF produced by the fault additional circuit.The fundamental MMF produced by the fault circuit increased with increasing k y,1 , but the fundamental MMF produced by the normal field circuit remained unchanged.The superposition of the two will show a decreasing trend, but the change is minimal overall.
The fundamental current induced in the stator winding and the fundamental MMF produced by the fundamental current will also have the same variation law.The combined fundamental MMF in the air gap will decrease with increasing k y,1 , but the overall change is also minimal.The size of the second MMF in the air gap is approximately proportional to k y,2 .Therefore, UMP is approximately proportional to pitch shortening factor k y,2 of the fault coil.For the same k y,2 , the size of UMP is negatively correlated with k y,1 .

Effects of Distributed Field Winding Short Circuit Turns on UMP
In order to study the effects of short circuit turns on UMP, 55 feet of the field winding concentric coil undergoing 10 to 80 turns were calculated.The calculation results of harmonic currents and UMP of branch a2 are shown in Table 6 and Figure 15. on 30 turn short circuit faults are respectively calculated.Conditions of operation are grid-connected.
The simulation results are shown in Table 5.When the two field concentric coils with the same pitch shortening factor ky,2 respectively undergo faults, the values of UMP are similar.The UMP with larger pitch shortening factor ky,1 is slightly less than that with smaller ky,1.
When the field concentric coils with different pitches undergo short circuit faults, UMP is generated by an interaction between the fundamental MMF and the second MMF in the air gap.The fundamental MMF includes the fundamental field MMF and the MMF produced by the stator harmonic current.The fundamental field MMF can be regarded as the sum of the fundamental MMF produced by the normal field circuit and the fundamental MMF produced by the fault additional circuit.The fundamental MMF produced by the fault circuit increased with increasing ky,1, but the fundamental MMF produced by the normal field circuit remained unchanged.The superposition of the two will show a decreasing trend, but the change is minimal overall.
The fundamental current induced in the stator winding and the fundamental MMF produced by the fundamental current will also have the same variation law.The combined fundamental MMF in the air gap will decrease with increasing ky,1, but the overall change is also minimal.The size of the second MMF in the air gap is approximately proportional to ky,2.Therefore, UMP is approximately proportional to pitch shortening factor ky,2 of the fault coil.For the same ky,2, the size of UMP is negatively correlated with ky,1.

Effects of Distributed Field Winding Short Circuit Turns on UMP
In order to study the effects of short circuit turns on UMP, 55 feet of the field winding concentric coil undergoing 10 to 80 turns were calculated.The calculation results of harmonic currents and UMP of branch a2 are shown in Table 6 and Figure 15.UMP shows the approximate linear increase variation law with the increasing short circuit turns.From Table 6, the fundamental current of the stator branch gradually reduces with increasing short circuit turns, but the change is small.The second harmonic and field currents increase with increasing short circuit turns.The growth rate of the former increases with an increasing number of short circuit turns.The latter is proportional to the number of short circuit turns.Therefore, the second harmonic MMF in the air gap will also increase, so UMP increases with increasing short circuited turns.The fundamental MMF changes little, and overall shows a decreasing trend.The increased rate of UMP slows, and the growth rate will not increase with increasing short circuited turns.This shows a variation that is approximately proportional to the short circuit turns.

Influence of Generator Running State on UMP
When a field winding inter-turn short circuit fault occurs in a generator, the position of the short circuit and the turn number is fixed in a certain time.With the change of the generator's running state, the UMP on the rotor will change.In order to determine whether or not the field winding inter-turn short circuit fault occurs, it is very important to study the relationship between UMP and the generator's running state.

Effect of Output Active Power on UMP
In order to study the effect of active power on UMP, its power factor angle (PFA) is changed.The output active power is also changed, and the reactive power is unchanged according to taps 2-4 short-circuit in the field winding of the machine model.The short circuit faults under these states are calculated.
As shown in Table 7, the output active power will decrease and UMP will increase with increasing PFA.This is due to the fundamental current of the stator branch reducing with increasing PFA ϕ under the premise of constant output reactive power.The size of the secondary harmonic current is essentially invariant.This is because there are no changes in the field current.The fundamental MMF and the secondary MMF produced by the field current, as well as the secondary MMF produced by the stator branch second harmonic current, are all unchanged.Only the fundamental MMF produced by the stator branch fundamental current becomes smaller.Because ϕ > 0, the fundamental MMF of the stator armature plays a demagnetizing role in field MMF.Therefore, the combined MMF in the air gap reduces and UMP decreases.As seen in Figure 16, UMP decreases faster with increasing active power.The fundamental MMF produced by the stator branch fundamental current will also increase.Demagnetizing armature reaction will also be stronger, and the saturation degree in the air gap magnetic field will decrease.The higher the air gap saturation, the lower the UMP change rate.The lower the air gap saturation, the greater the UMP change rate.
As seen in Figure 16, UMP decreases faster with increasing active power.The fundamental MMF produced by the stator branch fundamental current will also increase.Demagnetizing armature reaction will also be stronger, and the saturation degree in the air gap magnetic field will decrease.The higher the air gap saturation, the lower the UMP change rate.The lower the air gap saturation, the greater the UMP change rate.

Effect of Field Current on UMP
In order to study the effect of field current on UMP, its PFA is changed so that the field current and output reactive power are also changed.The active power is unchanged under the premise that taps 2-4 short circuit in the field winding of the machine model.The short circuit faults under these states are calculated.
As seen in Table 8 and Figure 17

Effect of Field Current on UMP
In order to study the effect of field current on UMP, its PFA is changed so that the field current and output reactive power are also changed.The active power is unchanged under the premise that taps 2-4 short circuit in the field winding of the machine model.The short circuit faults under these states are calculated.
As seen in Table 8 and Figure 17 the output reactive power and UMP are increased with increasing PFA.When the PFA increases (ϕ > ϕ), due to unchanged output active power (I N cos ϕ = I N cos ϕ ), the fundamental current of the stator branch and the output reactive power are increased (I N > I N , I N sin ϕ > I N sin ϕ).As the generator voltage .U remains unchanged, the field electromotive force becomes larger (E 0 > E 0 ), and the field current also increases (I f > I f ), as shown in Figure 18.The fundamental MMF and the second MMF produced by the rotor and stator will increase with increasing field current.Although the fundamental MMF produced by the stator branch current plays a demagnetizing role compared with the stator fundamental MMF, the influences of the fundamental field MMF on the combined MMF in the air gap are more obvious.UMP will also increase with increasing field current.The fundamental MMF and the second MMF produced by the rotor and stator will increase with increasing field current.Although the fundamental MMF produced by the stator branch current plays a demagnetizing role compared with the stator fundamental MMF, the influences of the fundamental field MMF on the combined MMF in the air gap are more obvious.UMP will also increase with increasing field current.

Conclusions
In this paper, a pair pole non-salient synchronous generator in ISCFW is calculated with the rapid calculation model of UMP and the finite element model.The comparison of the two results verifies the correctness and efficiency of the rapid calculation model of UMP.On this basis, the factors affecting UMP in ISCFW are analyzed.The following conclusions are drawn: (1) Under the premise of the same short circuit turns, UMP is approximately proportional to pitch shortening factor k y,2 of the fault coil.For the same k y,2 , the size of UMP is negatively correlated with k y,1 .Under the premise of the determined fault location, UMP is approximately proportional to the number of short circuit turns.(2) When the field current is fixed, UMP is negatively correlated with the output active power.
When the output active power is unchanged, UMP is positively correlated with the field current.

Figure 1 .
Figure 1.Magnetic motive force (MMF) produced by concentric field coil of field windings.

F
in the air gapgenerated by the fault winding-is equivalent to the superposition of the MMF generated by the normal field circuit and reverse MMF generated by the fault additional circuit, as shown in Equation (3):

Figure 2 .
Figure 2. Loops of field windings with inter-turn short circuit.

F
in the air gapgenerated by the fault winding-is equivalent to the superposition of the MMF generated by the normal field circuit and reverse MMF generated by the fault additional circuit, as shown in Equation (3):

Figure 4 .
Figure 4.The stator coordinate and stator slot position.

Figure 6 .
Figure 6.Rotor lamination and field winding's taps of A1552 model machine.Table 1. Parameters of A1552 non-salient pole model machine.Variable Parameter Variable Parameter Rated capacity (kW) 12 Number of series coils per branch of stator 6 Rated voltage (V) 400 Pitch of stator coil 15 Rated field current (A) 8.5 Number of conductors per slot of stator 12 No-load field current (A) 6.0 Number of rotor slots 24 Rated power factor 0.8 (hysteresis) Series turns per pole of field winding 330 Rated frequency (Hz) 50 Field winding parallel branch 1 Number of pole-pair (P) 1 Actual rotor slots 20 Gap length (mm) 1.5 Number of damp bars per pole of rotor 7 Number of parallel branches per phase of stator 2

Figure 6 .
Figure 6.Rotor lamination and field winding's taps of A1552 model machine.

Figure 7 .
Figure 7. Steady-state waveform in branch a2 under inter-turn circuit of taps 2-4 under load condition connected to power grid.(a) Simulation waveform; (b) Experimental waveform.

Figure 7 .
Figure 7. Steady-state waveform in branch a2 under inter-turn circuit of taps 2-4 under load condition connected to power grid.(a) Simulation waveform; (b) Experimental waveform.

Figure 8 .
Figure 8. Finite element model of the generator.

Figure 9 .
Figure 9. Finite element model meshed by triangle elements.

Figure 8 .
Figure 8. Finite element model of the generator.

Figure 8 .
Figure 8. Finite element model of the generator.

Figure 9 .
Figure 9. Finite element model meshed by triangle elements.Figure 9. Finite element model meshed by triangle elements.

Figure 9 .
Figure 9. Finite element model meshed by triangle elements.Figure 9. Finite element model meshed by triangle elements.

Figure 10 .
Figure 10.The flux density distribution of air gap in fault.

Figure 15 .
Figure 15.Relationship between short circuit turns and UMP.Figure 15.Relationship between short circuit turns and UMP.

Figure 15 .
Figure 15.Relationship between short circuit turns and UMP.Figure 15.Relationship between short circuit turns and UMP.

Figure 16 .
Figure 16.Relationship between UMP and output active power.
the output reactive power and UMP are increased with increasing PFA.When the PFA increases (    ), due to unchanged output active power (  ), the fundamental current of the stator branch and the output reactive power are increased ( the generator voltage U remains unchanged, the field electromotive force becomes larger ( 0 0 EE   ), and the field current also increases ( f f II   ), as shown in Figure 18.

Figure 16 .
Figure 16.Relationship between UMP and output active power.

Figure 17 .
Figure 17.Relationship between UMP and field current.

Figure 17 .
Figure 17.Relationship between UMP and field current.

Figure 17 .
Figure 17.Relationship between UMP and field current.

Table 1 .
Parameters of A1552 non-salient pole model machine.

Table 2 .
Different harmonics root mean square (RMS) of steady state current under inter-turn circuit

Table 2 .
: Different harmonics root mean square (RMS) of steady state current under inter-turn circuit of taps 2-4 under load condition connected to power grid.

Table 3 .
Rotor damping loop steady-state current RMS and phase of fault experiment.

Table 4 .
Stator armature steady state current RMS and phase of fault experiment.

Table 5 .
UMP and the RMS of stator harmonic current when field concentric coils with different pitches undergo the same turn short circuit.

Table 6 .
The RMS of the harmonic current of stator and rotor when field concentric coil undergoes different short circuit turns.

Table 5 .
UMP and the RMS of stator harmonic current when field concentric coils with different pitches undergo the same turn short circuit.

Table 7 .
Operation data of A1552 model machine field winding on taps 2-4 short circuit.PFA: power factor angle.

Table 8 .
Operation data of A1552 model machine field winding on taps 2-4 short circuit.

Table 8 .
Operation data of A1552 model machine field winding on taps 2-4 short circuit.