PLECS 2.2 Online Help

Saturable Induction Machine

Purpose

Induction machine with slip-ring rotor and main flux saturation.

Description

pict

The Saturable Induction Machine models main flux saturation by means of a continuous function. The model is implemented with machine equations in the stationary reference frame (Clarke transformation). Since the machine terminals have current-source characteristic, no external inductors may be connected. They must be included in the leakage inductances.

This machine model can be used with both the continuous and the discrete state-space method. To facilitate the discretization of the model the Integrator block is used which allows non-linear feedback in the circuit.

The machine operates as a motor or generator; the sign of the mechanical torque determines the mode of operation (positive for motoring, negative for generating). All electrical variables and parameters are viewed from the stator side. In the component icon, phase a of the stator and rotor winding is marked with a dot.

Electrical System:

pict

pict

The rotor flux is defined as

Yr,d = L ′lri′r,d + Lm (is,d + i′r,d)

                 (       )
Yr,q = L′lri′r,q +Lm   is,q + i′r,q  .

The electrical system of the machine model is implemented with state-variable equations that are derived from the equivalent circuit in the stationary reference frame. The value of the main flux inductance Lm   is not constant but depends on the main flux linkage Ym   as illustrated in the Ym ∕im   diagram.

pict

For flux linkages far below the transition flux YT  , the relationship between flux and current is almost linear and determined by the unsaturated magnetizing inductance Lm,0  . For large flux linkages the relationship is governed by the saturated magnetizing inductance Lm,sat  . YT   defines the knee of the transition between unsaturated and saturated main flux inductance. The tightness of the transition is defined with the form factor fT  . If you do not have detailed information about the saturation characteristic of your machine, fT = 1   is a good starting value. The function

plsaturation(Lm0, Lmsat, PsiT, fT)

plots the main flux vs. current curve and the magnetizing inductance vs. current curve for the parameters specified.

The model accounts for steady-state cross-saturation, i.e. the steady-state magnetizing inductances along the d-axis and q-axis are functions of the currents in both axes. In the implementation, the stator currents and the main flux linkage are chosen as state variables. With this type of model, the representation of dynamic cross-saturation can be neglected without affecting the machine's performance. The computation of the time derivative of the main flux inductance is not required.

In order to inspect the implementation, please select the component in your circuit and choose Look under mask from the Edit menu. If you want to make changes, you must first choose Break library link and then Unprotect, both from the Edit menu.

Electro-Mechanical System

Electromagnetic torque:

      3
Te =  2 p (is,qYs,d - is,dYs,q)

Mechanical System

Mechanical rotor speed ωm  :

ωm  =  1-(Te - Fωm - Tm)
       J

ω  = pωm

Mechanical rotor angle θm  :

θm  = ωm

θ = pθm

Parameters and Dialog Box

[Picture]

Stator resistance
Stator winding resistance Rs   in ohms (_O_  ).
Stator leakage inductance
Stator leakage inductance Lls   in henries (H).
Rotor resistance
Rotor winding resistance  ′
Rr    in ohms (_O_  ), referred to the stator side.
Rotor leakage inductance
Rotor leakage inductance L′lr   in henries (H), referred to the stator side.
Unsaturated magnetizing inductance
Unsaturated main flux inductance Lm,0  , in henries (H), referred to the stator side.
Saturated magnetizing inductance
Saturated main flux inductance Lm,sat    in henries (H), referred to the stator side. If you do not want to model saturation, set Lm,sat = Lm,0  .
Magnetizing flux at saturation transition
Transition flux linkage YT  , in Vs, defining the knee between unsaturated and saturated main flux inductance.
Tightness of saturation transition
Form factor fT   defining the tightness of the transition between unsaturated and saturated main flux inductance. The default is 1.
Inertia
Combined rotor and load inertia J   in Nms2  .
Friction coefficient
Viscous friction F   in Nms.
Number of pole pairs
Number of pole pairs p  .
Initial rotor speed
Initial mechanical rotor speed ωm,0   in s- 1  .
Initial rotor position
Initial mechanical rotor angle θm,0   in radians. If θm,0   is an integer multiple of 2π∕p   the stator windings are aligned with the rotor windings at simulation start.
Initial stator currents
A two-element vector containing the initial stator currents is,a,0   and is,b,0   of phases a and b in amperes (A).
Initial stator flux
A two-element vector containing the initial stator flux Ys,d,0   and Ys,q,0   in the stationary reference frame in Vs.

Inputs and Outputs

Mechanical torque
The input signal Tm   represents the mechanical torque at the rotor shaft, in Nm.

The output vector "m" contains the following 3 signals:

(1) Rotational speed
The rotational speed ωm   of the rotor in radians per second (s-1  ).
(2) Rotor position
The mechanical rotor angle θm   in radians.
(3) Electrical torque
The electrical torque Te   of the machine in Nm.

Probe Signals

Stator phase currents
The three-phase stator winding currents is,a  , is,b    and is,c  , in A. Currents flowing into the machine are considered positive.
Rotor phase currents
The three-phase rotor winding currents i′r,a  , i′r,b    and i′s,c   in A, referred to the stator side. Currents flowing into the machine are considered positive.
Stator flux (dq)
The stator flux linkages Y
  s,d    and Y
  s,q    in the stationary reference frame in Vs.
Magnetizing flux (dq)
The magnetizing flux linkages Ym,d    and Ym,q    in the stationary reference frame in Vs.
Rotor flux (dq)
The rotor flux linkages Y′r,d    and Y ′r,q    in the stationary reference frame in Vs, referred to the stator side.
Rotational speed
The rotational speed ωm   of the rotor in radians per second (s-1  ).
Rotor position
The mechanical rotor angle θm   in radians.
Electrical torque
The electrical torque Te   of the machine in Nm.

References

K. A. Corzine, B. T. Kuhn, S. D. Sudhoff, H. J. Hegner, "An improved method for incorporating magnetic saturation in the Q-D synchronous machine model", IEEE Transactions on Energy Conversion, Vol. 13, No. 3, Sept. 1998.
E. Levi, "A unified approach to main flux saturation modelling in D-Q axis models of induction machines", IEEE Transactions on Energy Conversion, Vol. 10, No. 3, Sept. 1995.
E. Levi, "Impact of cross-saturation on accuracy of saturated induction machine models", IEEE Transactions on Energy Conversion, Vol. 12, No. 3, Sept. 1997.