Summary

In the magnetic field analysis of a squirrel-cage induction machine (induction motor), a three-dimensional model is considered necessary to properly model the rotor bars and end rings in the rotor (conceptual diagram Fig. 1).　However, since 3D analysis requires a large amount of computation time and induction motors generally reach a steady state very slowly, 2D analysis that does not take into account the effects of the end rings is often used.　An example of induction motor analysis by using EMSolution has already been shown in "Induction Motor Analysis". There, the end rings are not included in the analysis and their effects (mainly due to resistance) are approximated by reducing the conductivity of the rotor bars. In this case, the area-integrated axial current in the rotor bar should be zero, so SUFCUR (Surface Inflow Current Source) is used to set the total axial current area-integrated value to zero. In the case of a two-dimensional analysis, the axial current integral value will generally not be zero unless such a constraint condition is put in place.

Fig.1 Rotor bars and end rings of squirrel-cage induction machine

Explanation

A more realistic approximation method for two-dimensional analysis is to simulate an end ring with a resistor and incorporate it into an electric circuit, coupling it with the finite element method.　Fig. 2 shows the equivalent circuit. The rectangles in the finite element analysis area surrounded by dotted lines represents a rotor bar.
The resistors represent end rings, which connect the rotor bars adjacent to each other at the top and bottom. Assuming that the rotor bars are equally spaced, the resistance values are all equal. In addition, it is possible to add resistance at the ends of the rotor bars or to insert inductance to include inductive effects, but for simplicity, we will only deal with resistance values here. Also, since the same current flows in the upper and lower sections, the lower section can be connected with zero resistance between the nodes (black circles in Fig. 2) and the upper section’s resistance can be doubled. An electrical circuit like the one in Fig. 2 can be easily created using the NETWORK module and can be coupled with the finite element method.　 For a 360-degree model without rotational period symmetry, it is sufficient to consider all rotor bars and end rings as shown in Fig. 2, but in many cases there is symmetry. The analysis can be performed by modeling only the periodic part (the mesh part).

As an example, let’s assume an N-fold symmetrical model with N series circuits on the primary side and model it for 1/N. If REGION_FACTOR=N, the primary input voltage and coil resistance can be set as a full cycle model. In this case, set the end ring resistance as N times so that the secondary currents flowing in the rotor bar and end ring are the same as in the full circumference model、so the resulting currents are equal to the primary and secondary currents flowing in the full circumference model. If REGION_FACTOR=1, the end ring resistance can be set as the full circumference model, but the primary input voltage and coil resistance must be set to 1/ N to make the resulting current equal to the primary and secondary current values flowing in the full circumference model. In this case, CURRENT for SUFCUR can be set to 1 in both cases.

Fig. 4 shows an example of a finite element mesh in which the end rings are approximated as resistance in the model used in "Induction Motor Analysis". The rotor bars of the rotor must all have different property numbers. The rotor bars on the periodically symmetric surface are rotationally symmetric, so the same property numbers are set for them. In "Induction Motor Analysis", the power circuit is entered using CIRCUIT and the conductivity of the rotor bars is set to take into account the effect of end rings, as described above. In this example, NETWORK is used here instead of CIRCUIT, and NETWORK allows the resistance of the end rings to be set separately.

To set up a three-phase AC power supply, you can also set up the circuits for the three phases independently, but here we consider a Y-shaped connection. When using CIRCUIT to set up, CONNECTION_MATRIX must be entered, which is not easy to set up intuitively. NETWORK can be used to set up the actual power circuit as it is. In this example, the Y-connection is made so that the voltage between phases is $100 \sqrt 3V$ (RMS value). If a current power input is used, the neutral point must be connected with a large resistor. In "Induction Motor Analysis", Delta-connection is used so that the voltage between phases is 100 V (RMS value). See "Delta and Y-connection" for an example of the wiring method and input, as well as some precautions for wiring a three-phase AC power supply.

Fig.2 Equivalent circuit of rotor bars and end rings of rotor

Fig.3 Equivalent circuit in the case with periodic symmetry

(a) Overall viewル

(b) Enlarged view of rotor section

Fig.4 Analytical mesh considering end rings by using NETWORK

The current and torque waveforms of the stator primary winding obtained from the analysis at slip S = 0.25 are shown in Figs. 5 and 6. It can be seen that this model is in good agreement with the results obtained by correcting the conductivity of the rotor bars in "Induction Motor Analysis". Fig. 7 shows the current waveform of the rotor bar for the property number shown in Fig. 4, and Fig. 8 shows the end ring current waveform connected to the rotor bar in Fig. 4.

In this way, the currents flowing in each rotor bar and the current waveforms of the end rings can be seen even in the two-dimensional analysis.

In reality, it is difficult to determine the resistance of the end ring in general, since it is considered to be frequency dependent due to skin effects and also depends on the shape of the end ring. Rather, it may be realistic to adjust the resistance value to better match the experimental results.

Fig.5 Current variation in stator primary winding at slip S=0.25

Fig.6 Torque acting on the rotor (half) at slip S=0.25

Fig.7 Rotor bar current variation at slip S=0.25

Fig.8 End ring current variation at slip S=0.25, speed =1125 rpm

How to use

1. In Handbook “16.1 Volume Element Properties," set the conductivity to the rotor bar property. Here, the conductivity of aluminum is set as it is.

* MAT_ID * POTENTIAL * B_H_CURVE_ID * SIGMA * MU * PACKING * ANISOTROPY * IRON_LOSS 51 0 0 2.9841e+007 1 1 0 0 * MAT_ID * POTENTIAL * B_H_CURVE_ID * SIGMA * MU * PACKING * ANISOTROPY * IRON_LOSS 52 0 0 2.9841e+007 1 1 0 0 ・・・ * MAT_ID * POTENTIAL * B_H_CURVE_ID * SIGMA * MU * PACKING * ANISOTROPY * IRON_LOSS 67 0 0 2.9841e+007 1 1 0 0

2. Set the rotor bar in Handbook "17.5 Surface Inflow Current Source (SUFCUR)".

   * SUFCUR * SERIES_NO * TIME_ID * SURFACE_MAT_ID * CURRENT * IN_ROTOR * SUFCUR 11 0 51 1 1 * SUFCUR * SERIES_NO * TIME_ID * SURFACE_MAT_ID * CURRENT * IN_ROTOR * SUFCUR 12 0 52 1 1 ・・・ * SUFCUR * SERIES_NO * TIME_ID * SURFACE_MAT_ID * CURRENT * IN_ROTOR * SUFCUR 27 0 67 1 1

3. In Handbook "17.8 Power Supply and Wiring (NETWORK)", connect the rotor bars set in 2 with resistors that simulate end rings. Also set up a three-phase AC power supply; since this is a 180-degree model, set REGION_FACTOR=2, which doubles the resistance of the end rings. The resistance of the end ring is $10.50 \mu \Omega$, calculated from the dimensions with reference to Reference [1].

******************************************************** * NETWORK * REGION_FACTOR * NETWORK 2 ******************************************************** * Rotor Bars ****************************************** * FEM * ID * NODE1 * NODE2 * SERIED_ID * FEM 1 1 201 11 FEM 2 2 202 12 ・・・ FEM 17 17 217 27 ********************************************************* * Rotor bars connected at the top and bottom with resistance that simulates end rings ** * Upper side ************************************************** * R * ID * NODE1 * NODE2 * RESISTANCE * R 21 1 2 2.1e-005 R 22 2 3 2.1e-005 ・・・ R 37 17 1 2.1e-005 * Bottom side ************************************************** R 41 201 202 2.1e-005 R 42 202 203 2.1e-005 ・・・ R 57 217 201 2.1e-005 * Three-phase AC power supply ****************************************** * FEM * ID * NODE1 * NODE2 * SERIED_ID * FEM 101 106 110 1 FEM 102 107 111 2 FEM 103 108 112 3 * R * ID * NODE1 * NODE2 * RESISTANCE * R 115 105 106 2.92 R 116 105 107 2.92 R 117 105 108 2.92 * VPS * ID * NODE1 * NODE2 * TIME_ID * VPS 130 120 110 1 VPS 131 120 111 2 VPS 132 120 112 3 ********************************************************* * END NETWORK * END *********************************************************

References

[1] "Electromagnetic Field Analysis Techniques for Virtual Engineering of Rotating Machines," IEEJ Technical Report No. 776 (2000/3)

Download

Analysis Model

・ inputNETWORK_Y_AC.ems : AC steady-state analysis
・ inputNETWORK_Y_1125rpm.ems : Transient analysis
・ pre_geom2D.neu : Stator mesh file
・ rotor_mesh2D : Rotor mesh file