This block implements a doubly-fed induction machine (DFIM) with magnetizing inductance saturation in the stationary (stator) reference frame along with the temperature effects on the stator and rotor resistances.

The DFIM block implements a three-phase induction machine (asynchronous machine) with an accessible wound rotor model with resolvers and encoders. The machine can operate in both motoring mode, when the mechanical torque is positive, and generating mode when the mechanical torque is negative.

Model Formulation

The DFIM model formulation is decribed in the sections below. Recall that the DFIM rotoric tension are null.

Unable to render {include} The included page could not be found.

Resolver Model

A resolver is a rotary transformer where the magnitude of the signal through the resolver windings varies sinusoidally as the shaft rotates.

Single excitation resover

The equations of the resolver can be expressed as follows:

where is the resolver angle, is the mechanical angle of the machine, is the angle offset, is the number of pole pairs of the resolver, , , , and are the resolver sine cosine gains and is the direction of the sensor rotation (0 = clockwise, 1 = counterclockwise).

Note: does not have an effect to the resolver angle, but change the resolver sine sign to match the counterclockwise convention.

When the gains are set to its default values (, , , and ), the sensor rotation is clockwise, and the excitation is a sinusoidal wave (, where is the excitation frequency), the outputs are given by the following expressions:

which means the excitaition signals are modulated by the sinus and the cosinus of the rotor position angle.

Excitation source

The excitation signal can be selected among one of the following options:

DC: a continuos signal is used as excitation, so the outputs are and ;
AC: a constant-frequency sinusoidal, unitary amplitude waveform is used as excitation , so the outputs are those expressed in the previous section.
External: a external signal, coming from one of the analog inputs, is used as excitation. This incoming signal may be rescaled by the Analog Input Differential Rescaling module (AIR).

Resolver gains

The resolver gains can be used to simulate non-ideal conditions of the sensor, occuring due to manufacturing imperfections.

Amplitude imbalance

Output phases have unequal inductances. Assuming the coupling between the rotor excitation winding and the cosine winding as the reference:

where .

Imperfect quadrature

Output phases are not in perfect spatial quadrature. Assuming the angle when the rotor excitation winding is aligned with the cosine winding:

where is the angle between the cosine winding and the sine winding (for there is a perfect quadrature, so one retrives the default gain values).

Faults

One of the windings is not connected due to a fault. In this case, the corresponding gain is equal to zero.

Excitation time shift

When the excitation signal comes from an external source (one of the analog input channels) the signal can be shifted by a given amount of time.

Paraemeters and Measurements

The parameters and measurements of the DFIM with saturation are separated in 4 different categories: Electrical, Mechanical, Resolver and Encoder.

Electrical Parameters and Measurements

Symbol	Name	Description	Unit	Type
R_snom	Stator resistance	Stator winding resistance of phase a, b, and c	Ω	Input
R_rnom'	Rotor resistance	Equivalent rotor winding resistance referred to the stator of phase a, b, and c	Ω	Input
Δ_θ	Temperature difference	Temperature difference w.r.t the initial temperature of stator and rotor windings	°C	Input
α_condstator	Stator temperature coefficient of resistance	Temperature coefficient of resistance of stator winding	°C^-1	Input
α_condrotor	Rotor temperature coefficient of resistance	Temperature coefficient of resistance of rotor winding	°C^-1	Input
L_ls	Stator leakage inductance	Stator winding leakage inductance of phase a, b, and c	H	Input
L_lr'	Rotor leakage inductance	Equivalent rotor winding leakage inductance referred to the stator of phase a, b, and c	H	Input
E_sat	Electrical saturation profile input	Choose your saturation profile between (more details) Stator Voltage Line to Line vs. Stator Current Magnetizing Inductance vs. Magnetizing Flux	N/A	Dropdown
V_sll	Stator voltage line to line table	No Load saturation curve parameters	V	File (more details)
L_m	Magnetizing inductance table	Stator-rotor mutual (magnetizing) inductance of phase a, b, and c	H	File (more details)
F_cLm	Cut-off frequency magnetizing inductance	Cut-off frequency associated with the frequency response of the magnetic coil, which acts as a filtering element. The default value is 3 kHz.	Hz	Input
N_sr	Turns ratio	Stator to rotor windings turns ratio	N/A	Input
pp	Number of pole pairs	Number of pole pairs	N/A	Input
i_s	Stator phase currents	Stator currents measured at phases a, b and c	A	Measurement
i_r	Rotor phase currents	Rotor equivalent phase a, b, and c currents, referred to the stator	A	Measurement
i_sα_β	Stator αβ currents	Stator currents in αβ reference frame	A	Measurement
i_sdq	Stator dq currents	Stator currents in dq reference frame	A	Measurement
i_rα_β	Rotor αβ currents	Rotor currents in αβ reference frame, referred to the stator	A	Measurement
i_rdq	Rotor dq currents	Rotor currents in dq reference frame, referred to the stator	A	Measurement
Frame_sel	Synchronous reference frame orientation	Selection of the synchronous reference frame depending on the use case: Stator flux-oriented reference frame (default value); Rotor flux-oriented reference frame; Air gap flux-oriented reference frame.	N/A	Dropdown
Frame_offset	Synchronous reference offset	Offset angle for synchronous reference frame	degree	Input
θ_est	Synchronous Angle	Estimated synchronous angle displayed depending of the orientation	degree	Measurement
V_s_α_β	Stator αβ voltages	Stator voltages in αβ reference frame	V	Measurement
Φ_s_α_β	Stator αβ fluxes	Stator fluxes in αβ reference frame	Wb	Measurement
Φ_r_α_β	Rotor αβ fluxes	Rotor fluxes in αβ reference frame, referred to the stator	Wb	Measurement
R_sn	Snubber resistance	Resistances of the snubber on phase a, b and c	Ω	Input
C_sn	Snubber capacitance	Capacitance of the snubber on phase a, b and c	F	Input

Mechanical Parameters and Measurements

Symbol	Name	Description	Unit	Type
J	Rotor inertia	Moment of inertia of the rotor	kg*m²	Input
F_v	Viscous friction coefficient	Viscous friction	Nms/rad	Input
F_s	Static friction torque	Static friction	N*m	Input
ctrl	Mechanical control elmode	Control mode of the mechanical model. Has two possible values: speed or torque. In speed mode, the mechanical model is bypassed and the speed command is sent directly. In torque mode, the torque command is used to measure the speed using the mechanical parameters of the machine.		Input
T	Torque command	Torque command sent to the mechanical model	N*m	Input
ω_rc	Rotor speed command	Speed command sent to the mechanical model	rpm	Input
ω_r	Rotor speed	Speed of the rotor	rpm	Measurement
T_e	Electromagnetic torque	Torque measured at the rotor	N*m	Measurement
θ₀	Initial rotor angle	Rotor position at time t = 0	°	Input
θ	Rotor angle	Rotor position from 0 to 360 degrees	°	Measurement

Resolver Parameters and Measurements

Symbol	Name	Description	Unit	Type
R_en	Enable resolver	Whether or not to enable the resolver	N/A	Input
R_sc	Resolver feedback signals	The two two-phase windings producing a sine and cosine feedback current proportional to the sine and cosine of the angle of the motor	N/A	Measurement
R_pp	Number of resolver pole pairs	Number of pole pairs of the resolver	N/A	Input
R_dir	Direction of the sensor rotation	Direction in which the sensor is turning, either clockwise or counterclockwise	N/A	Input
R_θ	Angle offset Δθ ( Sensor- Rotor )	Angle offset between the resolver and the rotor position from 0 to 360 degrees	°	Input
R_k	Resolver sine cosine gains	The sine/cosine modulation output sine/cosine component amplitude. Default value are 1, 0, 0 and 1	N/A	Input
E_type	Excitation source type	The source from which the excitation of the resolver is generated. Can either be AC, which is generated inside the FPGA with the specified frequency, DC, which is generated with a 90° from the rotor and External, which is generated from outside the model	N/A	Input
E_f	Excitation frequency	Frequency of the excitation when in AC mode	Hz	Input
E_src	Excitation source	Source of the external excitation source when in External mode	N/A	Input
E_ts	Excitation time shift	This parameter is used to compensate the time offset between the carrier generation's input in the system and modulated signals' output	s	Input

Encoder Parameters and Measurements

Symbol	Name	Description	Unit	Type
Enc_en	Enable encoder	Whether or not to enable the encoder	N/A	Input
Enc_type	Encoder type	Encoder type, either Quadrature or Hall Effect	N/A	Input
Q_ABZ	A B Z encoder signals	A B and Z signals of the encoder	N/A	Measurement
Q_ppr	Number of pulses per revolution	Number of pulses in one full revolution of the encoder	N/A	Input
Q_dir	Direction of the sensor rotation	Direction in which the sensor is turning, either A leads B or B leads A	N/A	Input
Q_θ	Angle offset Δθ ( Sensor - Rotor )	Angle offset between the encoder and the rotor position from 0 to 360 degrees	°	Input
Q_rat	Encoder speed ratio ( sensor to mechanical position )	Mechanical to encoder ratio. Angle of Encoder = Qrat * machine mechanical angle.	N/A	Input
H_θ	Hall effect sensor position	Position of sensor phases A, B and C in Hall effect mode	°	Input

Visualization of Resolver Encoder Parameters Effects

The number of resolver pole pairs affects the number of electrical turns per mechanical turns. On the left figure, the number of resolver pole pairs is 2, on the right figure, the number of resolver pole pairs is 4.

Resolver sine cosine gains affect the sine ( first axe ) and cosine ( second axe ) modulation output. Default values set to 1, 0, 0, 1 make it so the sine modulation has a sine form and the cosine modulation has a cosine form. If set to 0, 1, 1, 0, the sine modulation would have a cosine and the cosine modulation would have a sine form.

Excitation frequency, in AC excitation source type, affects the frequency of the carrier signal. We can see the time step highlighted in red. ( Figure 3 is a zoomed in view of Figure 2 )

Number of pulses per revolution (Q_ppr) defines how many times signals A and B pulse between two Z pulses ( one full rotation ).

Direction of the sensor rotation describes if A leads B ( Clockwise ) or if B leads A ( Counterclockwise )

Doubly-Fed Induction Machine with Saturation