Steady-State and Transient Performance Analysis of Permanent-Magnet Machines Using Time-Stepping Finite Element Technique

© 2012 Kurihara, licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Steady-State and Transient Performance Analysis of Permanent-Magnet Machines Using Time-Stepping Finite Element Technique


Introduction
Recently, global warming has become an important problem.High-efficiency machines have been needed in a large variety of industrial products in order to save electrical energy.For many applications, permanent-magnet (PM) synchronous machines can be designed which is smaller in size but more efficient as compared to induction machines [1][2][3].Besides, PMs have been employed as an alternative to current carrying coils for magnetic field excitation in synchronous machines for over 50 years.The lack of slip rings, brushes and field winding losses have always been viewed as distinct advantages over that of conventional wound field machines.However, when the machine size becomes small, the efficiency becomes low.This is mainly due to the reason that the iron loss and the copper loss are large, because the iron core of the stator in the small machine generally does not have annealing and the resistance of the stator windings is comparatively large.This chapter presents a successful design of the high-efficiency small but novel Interior permanent-magnet (IPM) machines using Neodymium-Boron-Iron (NdBFe) magnets.It is designed to operate with both high-efficiency line-start IPM motors [3] and generators with damper bars [4].Time-stepping finite element analysis has been used to successfully predict the dynamic and transient performances of the prototype machines.Time-stepping finite element analysis [3][4][5][6] has been used to successfully predict the dynamic and transient performances of the prototype IPM machines.The computed performance has been validated by tests in the prototype machine.

IPM machine configuration
The photograph of an IPM rotor, the cross section of a quarter of the high-efficiency motor and the demagnetization curve of the NdBFe magnet used for finite-element analysis are shown respectively in Figures 1-3, respectively [3].
A frame size of a 600 W, 3-phase, 4-pole, Y-connected, 50 Hz, 200 V squirrel-cage induction machine was used for testing the IPM rotor shown in Figure 1.The four-pole magnets arrangement in the rotor is oriented for a high-field type IPM synchronous machine.The experimentally developed rotor has the following distinctive design features [3]: 1.The fluxes from both sides of the magnet are concentrated effectively in the middle of the magnetic poles of the rotor.2. The reluctance of the d axis is larger than that of the q axis, because the d-axis flux passes across the magnet with high reluctance.Large reluctance torque can be obtained.3. The conducting material between the magnet and the rotor core is made from aluminum and has both functions of the flux barrier and cage bar.Furthermore, the number and configuration of rotor slots have been successfully designed by using the finite-element method so that the waveform of the electromotive force (EMF) due to the PMs was close to the sine waveform and the cogging torque was low.

Method for analysis
The analysis for taking the eddy currents into account, in general becomes essential to solve the three-dimensional problem.In this paper, it is assumed that the eddy currents flow approximately in the axial direction, because the rotor shown in Figure 1 is equipped with end rings.This reduces the analysis to a two-dimensional problem.The fundamental equations for the magnetic field are represented in the two-dimensional rectangular co-ordinates as where A is the z component of magnetic vector potential A, J0 is the stator-winding current density, Je is the eddy current density, Jm is the equivalent magnetizing current density, Mx, My are x and y components of the magnetization M, respectively. is the conductivity, and  is the reluctivity.The value of  in the PM is assumed the same as the reluctivity of free space 0.Jm is assumed zero, outside the PM.
The effect of the eddy current for the rotor ends is taken into account by multiplying by the coefficient kc as described below.It is done to reduce the analysis to two-dimensional.The equivalent resistance R2 for rotor bars including the rotor end rings can be given below if the bars are distributed at equal intervals in the rotor [7].
where Rb is the resistance of a bar, Re is the resistance of the end rings, Z2 is the number of rotor slots and p is the pole pair number.
Therefore, kc is given by This coefficient kc is found effective to take into account the rotor-bar current for the fundamental space harmonic.Moreover, it has been found that the agreement between computed and measured results of the starting performance characteristics in the IPM motor is good [3].Therefore, it is considered that design use of the kc is acceptable, even if the higher space harmonics exists [5].The value of coefficient kc is 0.55 in this paper.

Voltage, current and dynamic equations and calculation steps for IPM synchronous motor
where va, vb, and vc are the phase voltages, subscripts a, b, and c represent stator quantities in lines a, b, and c, respectively.vn is the potential of the neutral n, when the potential of the neutral of the supply source is zero, ia, ib, and ic are the line currents, r1 and L1 are the resistance and end-winding leakage inductance of the stator winding per phase, respectively.ea, eb, and ec are the induced phase voltages; and ea is given by the line integral of the vector potential round ca which is along the stator windings of phase a where A t is A at time t.t is the time step.eb, and ec can be obtained similarly, as in [5].vn can be obtained by adding each side of ( 6)-( 8) and then applying ( 9) and ( 11) One obtains the following equation by substituting (12) in ( 6)-( 8) [5]: The dynamic equation is given as [3] 0 where T is the instantaneous electromagnetic torque, J is the rotational inertia, r is the rotor angular speed, B0 is the friction coefficient, and Tl is the load torque.The torque T is calculated by using the Bil rule [8].The angular speed, r is given by where  is the rotational angle of the rotor.
One obtains the following equation by substituting (17) in (16): In this paper, the forward difference method is used to obtain the rotational angle at time t because the vector potential, currents and rotational angle at time t-t are all known One obtains the following equation by substituting ( 19) and ( 20) in (18) [6]: ].
In the case when the effect of the friction is negligibly small, the above equation can be represented simply as follows: 2.
Next, the calculation steps for this analysis are shown in Figure 5.
1. First, the terminal voltage Vl, its initial phase angle 0, Tl, and t are set, respectively.
Each voltage for the three stator phase windings can be represented by 2. The vector potential A at t = 0 is set, where the static field caused by only PMs is given as the initial value.
3. At t = t + t, the value of  t at new t is set.
4. At t = t + t, each voltage at new t is set. 5.The initial values for A t , ia t , ib t , and ic t are set.6.The matrix equation constructed by the time-stepping finite element technique is solved [5].7. The convergence of A t is tested.Unless A t converges, the process returns to step 6).
8. After the convergence of A t , ia t , ib t , and ic t , T t can be calculated.Then, the  t is determined from (22). 9.The calculation process from step 3) to step7) continues till the steady-state currents are obtained.

Voltage and current equations and calculation steps for IPM synchronous generator
Figure 6 shows the circuit of the three-phase IPM synchronous generator.The voltage and current equations of the IPM generator are given as where RL is a load resistance per phase.
  One can obtain the vector potential, currents by assuming a constant speed and then solving (1), ( 31)-(33) using the time-stepping finite-element technique [5].
Next, the calculation steps for this analysis are shown in Figure 7. 2. The vector potential A at time t = 0 is set, where the static field caused by only PMs is given as the initial value.3.At t = t + Δt, the initial values for A t , ia t , ib t and ic t at new t are set.4. The matrix equation constructed by the time-stepping finite-element technique is solved [5]. 5.The convergence of A t is tested.Unless A t converges, the process returns to step 4).6.After the convergence of A t , ia t , ib t and ic t are obtained.The calculation process from step 3) to 5) continues till the steady-state currents are obtained.

Steady-state synchronous and transient performance
This paper contains the steady-state synchronous and transient performance characteristics of the IPM synchronous machine shown in Figure 2. The good agreement between computed and measured results validates the proposed method for the finite-element analysis to predict the machine performance exactly.

Steady-state synchronous and transient performance of Line-start IPM synchronous motor
Figure 9 shows the load performance characteristics at 140V.It is clear from Figure 9 that the power factor is almost unity at all loads.The efficiency and power factor of the IPM motor were 86.2% and 0.986, respectively for the output of 600 W. The efficiency-power-factor product is 85.0%.It is about 35% higher than that for the induction motor.These values of the IPM motor are very high when compared to those of the induction motor for the same 600 W nameplate rating [3]. Figure 10 shows the computed and measured results of the input current versus the output power at 140V.It is shown that the agreement between the measured and computed results is excellent.In the figure 10, two kinds of computed curves are given, and the agreement is also good.It is used to determine the suitable value of the t.This value must be determined by taking into account the effects due to the space harmonics [5].The space harmonics effect is also the source of the cogging and ripple torques in the IPM motor.It can be compensated by skewing the stator by one slot pitch.Therefore, t.should be smaller than to ts to include the influence of the ripple harmonics on the starting with ts, which is to move by one stator slot pitch at synchronous speed of the motor where f is the line frequency, Ns is the number of stator slots.Herein, the following four values for the t are chosen: 208, 104, 52 and 26 μs are an eighth, a sixteenth, a thirty-second, and a sixty-fourth of ts, respectively.It is evident from Figure 10 that the choice of 208 μs is suitable at synchronous speed.However, this value is not sufficient in starting the IPM motor with large load inertia.
Figure 11 shows the computed speed-time responses at no load condition with the eddycurrent brake disc coupled to the shaft, when the stator of the motor was supplied with balanced three-phase voltages at rated frequency of 50 Hz and rated voltage 140 V.The inertia of the disc is about 18 times the experimental rotor inertia, and the initial phase angle 0 of (23-25) is /2 in the figure.It is seen that the agreement between the curves of 52 μs and 26 μs is good and that those are superposed.The choice of a time step of 52 μs is suitable when the starting of the IPM motor.
Figure 12 shows the computed and measured speed-time responses with time during runup and synchronizing period when t and 0 are 52 μs and /2, respectively.It can be seen that the good agreement between the measured and computed results exits.Figure 22 shows the measured results of the efficiency versus output.The efficiency was 85.8% at 600 W and 90% at 100W of light load.It is found that the efficiency is very high.

Conclusion
A successful design of a high-efficiency small but novel IPM machine with cage bars was developed and tested.It is designed to operate with both high-efficiency line-start IPM

Figure 4
Figure4shows the circuit of the three-phase line-start IPM synchronous motor.It has three stator phase windings, which are star connected with neutral.The voltage and current equations of the IPM motor are given as and Transient Performance Analysis of Permanent-Magnet Machines Using Time-Stepping Finite Element Technique 187

Figure 4 .
Figure 4. Circuit of three-phase line-start IPM synchronous motor

Figure 5 .
Figure 5. Flowchart of three-phase line-start IPM synchronous motor

Figure 6 .
Figure 6. Circuit of three-phase IPM synchronous generator

Figure 7 .
Figure 7. Flowchart of three-phase IPM synchronous generator 1.First, N, Δt and the corresponding rotational step s are set, respectively.

Figure 8
Figure 8 shows the terminal voltage waveform generated by PMs in driving the IPM synchronous machine at 1500 r/min by the external motor.It is shown that the agreement between the computed and measured values of the generated voltage is excellent.

Figure 9 .Figure 10 .
Figure 9. Measured results of load performance characteristics of IPM motor

Figure 13 .
Figure 13.Experimental setup for IPM synchronous generator

Figure 14 .Figure 15 .Figure 16 .
Figure 14.Measured results of terminal voltage versus time in IPM generator

Figure 17 .Figure 18 .
Figure 17.Computed results of line current versus time in IPM generator

Figure 20
Figure20shows the measured and computed results of the terminal voltage versus the output.It can be seen that the good agreement between the measured and computed values exists except near maximum output.

Figure 19 .
Figure 19.Steady-state line current versus time in IPM generator

Figure 20 .Figure 21 .Figure 22 .
Figure 20.Measured and computed results of terminal voltage versus output in IPM generator

Finite
Element Analysis -New Trends and Developments 200 motor and generator with damper bars.The IPM motor can start and synchronize fully with large load inertia.Beside, the effects of the damper bars on stability during load change and efficiency were investigated.Time-stepping finite element analysis has been used to successfully predict the steady-state and transient performances of the prototype IPM machines.It is clear that cage bars are used effectively to start up in a line-start IPM motor, and to operate stably in the IPM generator with damper bars.It has been found that the proposed design has yielded successful simulation and experimental results.