# Influence of End Structure on Electromagnetic Forces on End Winding of a 1550 MW Nuclear Generator

Authors: Chong Zeng; Song Huang; Yongming Yang; Guanghou Zhou
Journal: International Journal of Rotating Machinery
 (2017)

---

## Abstract

A 3D electromagnetic model of the end region of a 1550 MW nuclear generator is set up. The electromagnetic forces on the involute and nose parts of the end winding under a rated operation are obtained through the 3D time-step finite element method. The electromagnetic forces on different coils in the same phase are analyzed. By changing the rotor’s relative length and stator coil’s linear length in the 3D electromagnetic model, the electromagnetic force distributions on the end winding are obtained. The influence of each structure change on the electromagnetic force in different directions is studied in detail. Conclusions that can be helpful in decreasing the electromagnetic forces on the end winding through optimizing the end region design are presented.


---

## Body

## 1. Introduction

Vibrations caused by electromagnetic forces can severely damage end-winding insulations [1, 2]. Electromagnetic force increases significantly with the increase of generator capacity, thereby making large generators more vulnerable to vibration. Vibration should be decreased to ensure the security of end winding. Simply enhancing and increasing the support structures can decrease the damage of electromagnetic force but will lead to other problems, such as low utilization of materials and high manufacturing cost [3]. Optimizing the electromagnetic design for the end region can decrease electromagnetic force and vibration without extra cost. Research on the influence of end structure on the electromagnetic force on end winding can be a reference for optimization.The accurate calculation of electromagnetic field in the end region and electromagnetic force on end winding is the basis of studies on electromagnetic force. Many studies focus on calculation.In the early years, analytical methods were used to calculate the electromagnetic field in the end region. Methods based on the Biot-Savart law are the most popular among them [4–6]. However, the errors in the results obtained through analytical methods are large because the end region’s structure is complex. To improve the accuracy, numerical methods are utilized in the calculation. Recently, 3D finite element (FE) method has become the most popular method. Wang et al. [7] proposed full 3D models and methods for coupled electromagnetic and temperature fields in the end region of a large hydrogenerator. They discussed the effect of mesh size on the result of eddy current loss and obtained several meshing principles. Waldhart et al. [8] analyzed the influence of several common simplifications on the accuracy of the 3D FE method’s results and found that the geometric details on the inner diameter, the accurate modeling of the stator winding, and the correct BH-curve of the eddy current material are essential to the electromagnetic field. Liang et al. [9] built a fine 3D geometrical model and calculated the magnetic field and eddy current loss in the end region. Wang et al. [10–14] calculated the electromagnetic field and the eddy current losses in the end region of a 330 MW turbogenerator. The influence of underexcitation [10], copper shield [11–13], and clamping plate [14] on the electromagnetic field and eddy current losses was calculated and discussed. Using the 3D time-step FE method, Huang et al. calculated the magnetic field and the eddy current loss on the end components of a 1550 MW nuclear generator and analyzed the influences of different structural parameters on the magnetic field and the eddy current loss [15].Relatively few studies focus on the electromagnetic force on end winding. Simplifying the end coil as straight line inductors, Richard et al. [16, 17] calculated the electromagnetic force on the end winding of a 600 MW turbogenerator using the 3D electromagnetic model of the stator end region. Senske et al. [18], Drubel et al. [19], and Grüning and Kulig [20] calculated the electromagnetic force on end winding based on the Biot-Savart law to line circuit segments. Albanese et al. [21] proposed a numerical approach based on an integral formulation to calculate the magnetic field and forces on the end winding of a large generator. Zhao et al. [22] built a 3D electromagnetic model of a 600 MW turbogenerator and calculated the electromagnetic force on the end winding using the 3D time-step FE method. Furthermore, the vibrations of the end winding were calculated and subjected to 3D FE structural analysis. None of the studies above analyzed the influence of the end structure on the electromagnetic force on end winding.In this paper, a 3D FE model of the end region of a 1550 MW nuclear generator is set up. The electromagnetic force on the end winding at rated operation is obtained by the FE method. The rotor relative length and the stator coil linear length are set to different values in the FE model. The radial, tangential, and axial force density distributions in different parts of the end winding are presented. The influence of the two parameters on the electromagnetic force on end winding is analyzed. Conclusions that can help in decreasing the electromagnetic forces on end winding through optimal design are presented.


## 2. 3D Finite Element Model for Electromagnetic Forces Calculation

### 2.1. Physical Model

The end region of a half-speed nuclear generator is quite complex. The end part of cores and windings is 3D in structure. The end region contains three end components: copper shield, press plate, and press finger. The structure of the end region is shown in Figure1.Figure 1
End region of a half-speed nuclear generator.


### 2.2. Mathematical Model for Electromagnetic Forces Calculation

To improve computational efficiency, the following are assumed:(i)
The distribution of current in the stator and rotor windings is uniform.(ii)
The high-order harmonics of current and displacement current are ignored.(iii)
The hysteresis effect is ignored. The core material is isotropic [11].(iv)
The cross section of the solution domain of the magnetic field calculation is shown in Figure2.Figure 2
Cross section of the solution domain.The solution domain is shown in Figure2. S1, S2, and S3 are the boundaries of the solution domain. By adopting A, φ-A method, we can express the mathematical model of the end region as follows [23, 24]:(1)∇×1μ∇×A+σ∂A∂t+σ∇φ=Js,∇·-σ∂A∂t-σ∇φ=0,where A denotes the magnetic vector potential, φ is the electric scalar potential, σ is the electrical conductivity, ω is the angular frequency, and Js is the source current density.The axial flux is very small inS1. Therefore, the symmetric boundary condition, which means no flux passes across the face, is set on S1. S2 and S3 are far from the stator and rotor cores, and their magnetic fields are weak. Therefore, S2 and S3 satisfy the first boundary condition A=0.The stator current cannot be calculated from the above model because the solution domain only contains the end region. Thus, an additional 2D model of the generator is required. The solution domain of the 2D model is shown in Figure3, where A1 is the inner area and E1 is the outer edge of the solution domain.Figure 3
Solution domain of the 2D model.The mathematical model can be described as follows:(2)∂∂xv∂A∂x+∂∂yv∂A∂y=-Jsσ∇φ=JsinA1A=0onE1,J
s is determined by the stator current and rotor current. When a generator operates at rated condition, the rotor current is constant and the 3-phase stator currents are given as follows:(3)iA=Imsin⁡ωt+θ,iB=Imsin⁡ωt+θ-120°,iC=Imsin⁡ωt+θ-240°,where θ is the angle between iA and d-axis, ω is the synchronous angular velocity, t is time, and Im is the magnitude of stator phase current.I
m can be calculated with the following formula:(4)Im=P3/2Ucos⁡φ,where P is the rated power, U is the rated voltage, and φ is the power factor angle at rated condition.Besides,θ and the rotor current can be calculated by iteration.With the magnetic field in the end region, the electromagnetic force density can be calculated as(5)f=IS×B,where f is the electromagnetic force density, I is the current in the coil, S is the cross-sectional area of the coil, and B is the magnetic flux density.


### 2.3. Modeling of the Proposed Generator

A 1550 MW nuclear generator is adopted because the electromagnetic force in a large generator is large. The basic parameters of the generator are shown in Table1.Table 1
Parameters of the 1550 MW nuclear generator.
Parameters
Values
Power
1550 MW
Voltage
27 kV
Frequency
50 Hz
Number of poles
4
Current
36926.9 A
Power factor
0.9 (lagging)JMAG 12.1 is used to build the model and perform the calculation. A 3D FE model of the end region is established per the actual design of the proposed generator. The meshing of the FE model is shown in Figure4. Using a sliding mesh, the rotation of the rotor is considered. The whole model contains 713386 elements and 129854 nodes.Figure 4
FE model mesh of the proposed generator end region. (a) Stator. (b) Rotor.
(a)
(b)A 2D FE model of the generator is set up with ANSYS Maxwell 16.1. The stator and rotor currents are obtained by iteration at a rated condition. The result is shown in Table2.Table 2
Iteration result of stator and rotor currents.
Parameters
Values
Rotor current
6668 A
θ
−246°


### 2.4. Validation of the Model

The experiment is performed on an 1150 MW nuclear generator because the 1550 MW nuclear generator is not available for the experiment. The FE model of the 1150 MW generator is built on the same principles, and the structure of the 1150 MW nuclear generator is almost the same as but smaller than that of the 1550 MW generator. Therefore, the relative errors of the two FE models have a similar order. Three probes are set as shown in Figure5. The comparison of the magnetic flux densities at no-load operation between the computational and experimental results is shown in Table 3. Table 3 shows that the calculation method is reasonable and the relative error of the calculated result is within 10%.Table 3
Comparison between the calculated and measured value of magnetic density.
Probes
Calculated result
Measured result
Relative error
A
0.25 T
0.27 T
−7.41%
B
0.70 T
0.76 T
−7.89%
C
0.77 T
0.79 T
−2.53%Figure 5
Distribution of the test points.


## 2.1. Physical Model

The end region of a half-speed nuclear generator is quite complex. The end part of cores and windings is 3D in structure. The end region contains three end components: copper shield, press plate, and press finger. The structure of the end region is shown in Figure1.Figure 1
End region of a half-speed nuclear generator.


## 2.2. Mathematical Model for Electromagnetic Forces Calculation

To improve computational efficiency, the following are assumed:(i)
The distribution of current in the stator and rotor windings is uniform.(ii)
The high-order harmonics of current and displacement current are ignored.(iii)
The hysteresis effect is ignored. The core material is isotropic [11].(iv)
The cross section of the solution domain of the magnetic field calculation is shown in Figure2.Figure 2
Cross section of the solution domain.The solution domain is shown in Figure2. S1, S2, and S3 are the boundaries of the solution domain. By adopting A, φ-A method, we can express the mathematical model of the end region as follows [23, 24]:(1)∇×1μ∇×A+σ∂A∂t+σ∇φ=Js,∇·-σ∂A∂t-σ∇φ=0,where A denotes the magnetic vector potential, φ is the electric scalar potential, σ is the electrical conductivity, ω is the angular frequency, and Js is the source current density.The axial flux is very small inS1. Therefore, the symmetric boundary condition, which means no flux passes across the face, is set on S1. S2 and S3 are far from the stator and rotor cores, and their magnetic fields are weak. Therefore, S2 and S3 satisfy the first boundary condition A=0.The stator current cannot be calculated from the above model because the solution domain only contains the end region. Thus, an additional 2D model of the generator is required. The solution domain of the 2D model is shown in Figure3, where A1 is the inner area and E1 is the outer edge of the solution domain.Figure 3
Solution domain of the 2D model.The mathematical model can be described as follows:(2)∂∂xv∂A∂x+∂∂yv∂A∂y=-Jsσ∇φ=JsinA1A=0onE1,J
s is determined by the stator current and rotor current. When a generator operates at rated condition, the rotor current is constant and the 3-phase stator currents are given as follows:(3)iA=Imsin⁡ωt+θ,iB=Imsin⁡ωt+θ-120°,iC=Imsin⁡ωt+θ-240°,where θ is the angle between iA and d-axis, ω is the synchronous angular velocity, t is time, and Im is the magnitude of stator phase current.I
m can be calculated with the following formula:(4)Im=P3/2Ucos⁡φ,where P is the rated power, U is the rated voltage, and φ is the power factor angle at rated condition.Besides,θ and the rotor current can be calculated by iteration.With the magnetic field in the end region, the electromagnetic force density can be calculated as(5)f=IS×B,where f is the electromagnetic force density, I is the current in the coil, S is the cross-sectional area of the coil, and B is the magnetic flux density.


## 2.3. Modeling of the Proposed Generator

A 1550 MW nuclear generator is adopted because the electromagnetic force in a large generator is large. The basic parameters of the generator are shown in Table1.Table 1
Parameters of the 1550 MW nuclear generator.
Parameters
Values
Power
1550 MW
Voltage
27 kV
Frequency
50 Hz
Number of poles
4
Current
36926.9 A
Power factor
0.9 (lagging)JMAG 12.1 is used to build the model and perform the calculation. A 3D FE model of the end region is established per the actual design of the proposed generator. The meshing of the FE model is shown in Figure4. Using a sliding mesh, the rotation of the rotor is considered. The whole model contains 713386 elements and 129854 nodes.Figure 4
FE model mesh of the proposed generator end region. (a) Stator. (b) Rotor.
(a)
(b)A 2D FE model of the generator is set up with ANSYS Maxwell 16.1. The stator and rotor currents are obtained by iteration at a rated condition. The result is shown in Table2.Table 2
Iteration result of stator and rotor currents.
Parameters
Values
Rotor current
6668 A
θ
−246°


## 2.4. Validation of the Model

The experiment is performed on an 1150 MW nuclear generator because the 1550 MW nuclear generator is not available for the experiment. The FE model of the 1150 MW generator is built on the same principles, and the structure of the 1150 MW nuclear generator is almost the same as but smaller than that of the 1550 MW generator. Therefore, the relative errors of the two FE models have a similar order. Three probes are set as shown in Figure5. The comparison of the magnetic flux densities at no-load operation between the computational and experimental results is shown in Table 3. Table 3 shows that the calculation method is reasonable and the relative error of the calculated result is within 10%.Table 3
Comparison between the calculated and measured value of magnetic density.
Probes
Calculated result
Measured result
Relative error
A
0.25 T
0.27 T
−7.41%
B
0.70 T
0.76 T
−7.89%
C
0.77 T
0.79 T
−2.53%Figure 5
Distribution of the test points.


## 3. Electromagnetic Forces on End Winding

We should analyze only the electromagnetic force on one winding phase because the generator and the electromagnetic force on each winding phase are symmetric. The electromagnetic forces on phase A are analyzed in this paper. In the proposed generator, a winding phase has eight coils. The eight coils are numbered as shown in Figure6. An end coil can be divided into several parts and the definition of these parts is shown in Figure 7.Figure 6
Numbering of phase A coils.Figure 7
Definition of the parts of the end coil.The magnetic fields on different coils in the same phase are different. Therefore, the electromagnetic forces on different coils are also different, as shown in Figure8. The forces on the involute parts change significantly. Generally, the forces on the upper involute parts increase with the coil number. Conversely, the forces on the lower involute parts decrease with the coil number. The upper involute part suffers the largest force among all the involute parts. The electromagnetic forces on the nose parts are smaller than the forces on the involute parts, and they change slightly.Figure 8
Electromagnetic forces on different coils of phase A.


## 4. Parametric Study on Electromagnetic Forces on End Winding

Although the force on the nose part is smaller than that on the involute part, the vibration on the nose part is even larger than that on the involute part because of the weak constraint [22]. The influences of the end structure parameters on the electromagnetic force on the involute part and the force on the nose part are analyzed in this section. To simplify the analysis, only the forces on coil 8 are studied in detail. Only the upper involute part and the nose are analyzed because the constraints of the upper and lower involute parts are almost the same.


### 4.1. Influence of the Relative Length of Rotor

Rotor relative length (lr) is defined as shown in Figure 9. The lengths of the rotor core and winding change with lr. When the rotor core is longer than the stator core, lr is positive. Conversely, when the rotor core is shorter than stator core, lr is negative. To study the influence of the rotor relative length on electromagnetic force, lr is set to −75, 0, and 75 mm in the FE model.Figure 9
Definition of rotor relative length.The electromagnetic force distribution on the involute part at 0.039 s is shown in Figure10. Most of the forces reach their maximums at the same time because the current in the coil reaches the maximum at 0.039 s. Figure 10 shows that radial force is the largest and increases with lr. The tangential and axial forces are smaller than the radial force and decrease with the increase of lr. For the three components, the changes are more significant at the side close to the linear part (the right side of each involute part in Figure 10). When lr varies from −75 mm to +75 mm, the maximum of the electromagnetic force density varies from 2.3 × 106 N/m3 to 2.6 × 106 N/m3, 1.2 × 106 N/m3 to 1.0 × 106 N/m3, and 1.8 × 106 N/m3 to 1.4 × 106 N/m3 in the radial, tangential, and axial directions, respectively.Figure 10
Electromagnetic force density distribution on the involute part at 0.039 s, whenlr = −75 mm (bottom), lr = 0 mm (middle), and lr = 75 mm (top). (a) Radial component. (b) Tangential component. (c) Axial component.
(a)
(b)
(c)Figure11 shows the maximum electromagnetic forces on the involute part with different lr. Figure 11 shows that the radial force is relatively larger than the tangential and axial forces. When lr increases from −75 mm to +75 mm, the radial force increases by 13% while the tangential and axial forces decrease by 9% and 12%, respectively.Figure 11
Electromagnetic forces on the involute part with different rotor relative lengths.The electromagnetic force distribution on the nose part is shown in Figure12. The magnitude of the force density is used to show the change in force strength because the direction of the force density changes to the opposite direction of the axis on some areas of the nose part. Figure 12 shows that the radial and axial force densities increase with lr. The radial force change mainly appears at the inner side of the nose part, and the axial force change mainly appears at the inner side of the top section of the nose part. When lr varies from −75 mm to +75 mm, the maximum of the electromagnetic force density varies from 1.9 × 106 N/m3 to 2.0 × 106 N/m3, 1.0 × 106 N/m3 to 1.0 × 106 N/m3, and 2.1 × 106 N/m3 to 2.2 × 106 N/m3 in the radial, tangential, and axial directions, respectively.Figure 12
Electromagnetic force density magnitude distribution on the nose part at 0.039 s, whenlr = −75 mm (bottom), lr = 0 mm (middle), and lr = 75 mm (top). (a) Radial component. (b) Tangential component. (c) Axial component.
(a)
(b)
(c)The maximum electromagnetic forces on the nose part with different rotor relative lengths are shown in Figure13. The radial force is smaller than the tangential and axial forces because the direction of the radial force on both sides of the nose part are opposite and almost offset each other. When lr varies from −75 mm to +75 mm, the radial and axial forces increase by 8% and 23%, respectively, while the tangential force almost stays constant.Figure 13
Electromagnetic forces on the nose part with different rotor relative lengths.For large generators, the radial vibration of end winding is the most serious [22]. A relatively short rotor can be adopted to reduce the radial vibration from the source. Both the electromagnetic force on the involute and nose parts are smaller with a relatively short rotor. The axial force on the nose part (especially on the top of the inner side) also decreases. However, the tangential and axial forces on the involute part (especially on the region close to the linear part of the end coil) increase at the same time.


### 4.2. Influence of Stator Coil Linear Length

The stator coil linear length (lsc) is defined in Figure 14 as the length of the lower linear part beyond the stator core. To study the influence of lsc on electromagnetic force, lsc is set to 221, 271, 321, 371, and 421 mm in the FE model.Figure 14
Definition of stator coil linear length.The electromagnetic force density distributions on the involute part at 0.039 s with different stator coil linear lengths are shown in Figure15. Figure 15 shows that the radial force decreases slightly with the increasing lsc, and the maximum of the radial electromagnetic force is reduced to 2.3 × 106 N/m3 when lsc is 421 mm from 2.7 × 106 N/m3 when lsc is 221 mm. The tangential and axial forces increase significantly with lsc. The maximums of the tangential and axial forces increase to 1.1 × 106 N/m3 and 1.8 × 106 N/m3, respectively, when lsc is 421 mm from 0.9 × 106 and 1.3 × 106 N/m3, respectively, when lsc is 221 mm.Figure 15
Electromagnetic force density distribution on the involute part at 0.039 s, whenlsc = 221 mm (bottom), lsc = 321 mm (middle), and lsc = 421 mm (top). (a) Radial component. (b) Tangential component. (c) Axial component.
(a)
(b)
(c)Figure16 shows the maximum electromagnetic forces on the involute part with different lsc. Figure 16 shows that the radial force is relatively larger than the tangential and axial forces. When lsc increases from 221 mm to 421 mm, the radial force decreases by 14% while the tangential and axial forces increase by 17% and 24%, respectively.Figure 16
Electromagnetic forces on the involute part with different stator coil linear length.The distribution of the electromagnetic force density magnitude on the nose part with differentlsc is shown in Figure 17. The influence of the stator coil linear length on electromagnetic force on the nose part is small. Generally, the radial force decreases with the increase of the stator coil’s linear length. Conversely, the tangential force increases with the stator coil’s linear length. The axial force is constant against the changes in stator coil linear length.Figure 17
Electromagnetic force density magnitude distribution on the nose part at 0.039 s, whenlsc = 221 mm (bottom), lsc = 321 mm (middle), and lsc = 421 mm (top). (a) Radial component. (b) Tangential component. (c) Axial component.
(a)
(b)
(c)Figure18 shows the maximum electromagnetic forces on the nose part with different lsc. When lsc increases from 221 mm to 421 mm, the radial and axial forces decrease by 9% and 17%, respectively, while the tangential force remains almost constant.Figure 18
Electromagnetic forces on the nose part with different stator coil linear lengths.To reduce the radial vibration on end winding, a relatively large stator coil linear length can be adopted. Consequently, the radial force on the involute and nose parts and the axial force on the nose part will decrease. In contrast, the tangential and axial forces on the involute part increase.


## 4.1. Influence of the Relative Length of Rotor

Rotor relative length (lr) is defined as shown in Figure 9. The lengths of the rotor core and winding change with lr. When the rotor core is longer than the stator core, lr is positive. Conversely, when the rotor core is shorter than stator core, lr is negative. To study the influence of the rotor relative length on electromagnetic force, lr is set to −75, 0, and 75 mm in the FE model.Figure 9
Definition of rotor relative length.The electromagnetic force distribution on the involute part at 0.039 s is shown in Figure10. Most of the forces reach their maximums at the same time because the current in the coil reaches the maximum at 0.039 s. Figure 10 shows that radial force is the largest and increases with lr. The tangential and axial forces are smaller than the radial force and decrease with the increase of lr. For the three components, the changes are more significant at the side close to the linear part (the right side of each involute part in Figure 10). When lr varies from −75 mm to +75 mm, the maximum of the electromagnetic force density varies from 2.3 × 106 N/m3 to 2.6 × 106 N/m3, 1.2 × 106 N/m3 to 1.0 × 106 N/m3, and 1.8 × 106 N/m3 to 1.4 × 106 N/m3 in the radial, tangential, and axial directions, respectively.Figure 10
Electromagnetic force density distribution on the involute part at 0.039 s, whenlr = −75 mm (bottom), lr = 0 mm (middle), and lr = 75 mm (top). (a) Radial component. (b) Tangential component. (c) Axial component.
(a)
(b)
(c)Figure11 shows the maximum electromagnetic forces on the involute part with different lr. Figure 11 shows that the radial force is relatively larger than the tangential and axial forces. When lr increases from −75 mm to +75 mm, the radial force increases by 13% while the tangential and axial forces decrease by 9% and 12%, respectively.Figure 11
Electromagnetic forces on the involute part with different rotor relative lengths.The electromagnetic force distribution on the nose part is shown in Figure12. The magnitude of the force density is used to show the change in force strength because the direction of the force density changes to the opposite direction of the axis on some areas of the nose part. Figure 12 shows that the radial and axial force densities increase with lr. The radial force change mainly appears at the inner side of the nose part, and the axial force change mainly appears at the inner side of the top section of the nose part. When lr varies from −75 mm to +75 mm, the maximum of the electromagnetic force density varies from 1.9 × 106 N/m3 to 2.0 × 106 N/m3, 1.0 × 106 N/m3 to 1.0 × 106 N/m3, and 2.1 × 106 N/m3 to 2.2 × 106 N/m3 in the radial, tangential, and axial directions, respectively.Figure 12
Electromagnetic force density magnitude distribution on the nose part at 0.039 s, whenlr = −75 mm (bottom), lr = 0 mm (middle), and lr = 75 mm (top). (a) Radial component. (b) Tangential component. (c) Axial component.
(a)
(b)
(c)The maximum electromagnetic forces on the nose part with different rotor relative lengths are shown in Figure13. The radial force is smaller than the tangential and axial forces because the direction of the radial force on both sides of the nose part are opposite and almost offset each other. When lr varies from −75 mm to +75 mm, the radial and axial forces increase by 8% and 23%, respectively, while the tangential force almost stays constant.Figure 13
Electromagnetic forces on the nose part with different rotor relative lengths.For large generators, the radial vibration of end winding is the most serious [22]. A relatively short rotor can be adopted to reduce the radial vibration from the source. Both the electromagnetic force on the involute and nose parts are smaller with a relatively short rotor. The axial force on the nose part (especially on the top of the inner side) also decreases. However, the tangential and axial forces on the involute part (especially on the region close to the linear part of the end coil) increase at the same time.


## 4.2. Influence of Stator Coil Linear Length

The stator coil linear length (lsc) is defined in Figure 14 as the length of the lower linear part beyond the stator core. To study the influence of lsc on electromagnetic force, lsc is set to 221, 271, 321, 371, and 421 mm in the FE model.Figure 14
Definition of stator coil linear length.The electromagnetic force density distributions on the involute part at 0.039 s with different stator coil linear lengths are shown in Figure15. Figure 15 shows that the radial force decreases slightly with the increasing lsc, and the maximum of the radial electromagnetic force is reduced to 2.3 × 106 N/m3 when lsc is 421 mm from 2.7 × 106 N/m3 when lsc is 221 mm. The tangential and axial forces increase significantly with lsc. The maximums of the tangential and axial forces increase to 1.1 × 106 N/m3 and 1.8 × 106 N/m3, respectively, when lsc is 421 mm from 0.9 × 106 and 1.3 × 106 N/m3, respectively, when lsc is 221 mm.Figure 15
Electromagnetic force density distribution on the involute part at 0.039 s, whenlsc = 221 mm (bottom), lsc = 321 mm (middle), and lsc = 421 mm (top). (a) Radial component. (b) Tangential component. (c) Axial component.
(a)
(b)
(c)Figure16 shows the maximum electromagnetic forces on the involute part with different lsc. Figure 16 shows that the radial force is relatively larger than the tangential and axial forces. When lsc increases from 221 mm to 421 mm, the radial force decreases by 14% while the tangential and axial forces increase by 17% and 24%, respectively.Figure 16
Electromagnetic forces on the involute part with different stator coil linear length.The distribution of the electromagnetic force density magnitude on the nose part with differentlsc is shown in Figure 17. The influence of the stator coil linear length on electromagnetic force on the nose part is small. Generally, the radial force decreases with the increase of the stator coil’s linear length. Conversely, the tangential force increases with the stator coil’s linear length. The axial force is constant against the changes in stator coil linear length.Figure 17
Electromagnetic force density magnitude distribution on the nose part at 0.039 s, whenlsc = 221 mm (bottom), lsc = 321 mm (middle), and lsc = 421 mm (top). (a) Radial component. (b) Tangential component. (c) Axial component.
(a)
(b)
(c)Figure18 shows the maximum electromagnetic forces on the nose part with different lsc. When lsc increases from 221 mm to 421 mm, the radial and axial forces decrease by 9% and 17%, respectively, while the tangential force remains almost constant.Figure 18
Electromagnetic forces on the nose part with different stator coil linear lengths.To reduce the radial vibration on end winding, a relatively large stator coil linear length can be adopted. Consequently, the radial force on the involute and nose parts and the axial force on the nose part will decrease. In contrast, the tangential and axial forces on the involute part increase.


## 5. Conclusions

In this paper, a 3D FE model of the end region of a 1550 MW nuclear generator is set up. The electromagnetic force on the stator end winding at a rated operation is calculated. The electromagnetic forces on different coils in the same phase are presented. By changing the rotor relative length and stator coil linear length in the FE model, the influence of these parameters on electromagnetic force is analyzed in detail. The following conclusions are drawn:(1)
The electromagnetic force on the involute part varies in different coils at the same phase and the upper involute part of the last coil in the rotation direction suffers the largest force. The electromagnetic force on the nose part is smaller than that on the involute part and changes slightly in different coils.(2)
Decreasing the rotor relative length decreases the radial forces on the involute and nose parts. However, decreasing the rotor relative length increases the tangential and axial forces on the involute parts.(3)
The stator coil linear length can be increased to decrease the radial forces on the end winding. The radial force on the involute and nose parts decreases significantly but the tangential and axial forces on the involute parts simultaneously increase.



---

*Source: 9545238-2017-10-29.xml*