Salient-pole synchronous machine with fundamental or standard parameterization
Simscape / Electrical / Electromechanical / Synchronous
The Synchronous Machine Salient Pole block models a salient-pole synchronous machine using fundamental or standard parameters.
If the block is in a network that is compatible with the frequency-time simulation mode, you can perform a load-flow analysis on the network. A load-flow analysis provides steady-state values that you can use to initialize the machine.
For more information, see Perform a Load-Flow Analysis Using Simscape Electrical and Frequency and Time Simulation Mode. For an example that shows how initialize an synchronous machine using data from a load flow analysis, see Synchronous Machine Initialization with Loadflow.
The synchronous machine equations are expressed with respect to a rotating reference frame, defined by
where:
θe is the electrical angle.
N is the number of pole pairs.
θr is the rotor angle.
The Park transformation maps the synchronous machine equations to the rotating reference frame with respect to the electrical angle. The Park transformation is defined by
The Park transformation is used to define the per-unit synchronous machine equations. The stator voltage equations are defined by
and
where:
ed, eq, and e0 are the d-axis, q-axis, and zero-sequence stator voltages, defined by
va, vb, and vc are the stator voltages measured from port ~ to neutral port n.
ωbase is the per-unit base electrical speed.
ψd, ψq, and ψ0 are the d-axis, q-axis, and zero-sequence stator flux linkages.
ωr is the per-unit rotor rotational speed.
Ra is the stator resistance.
id, iq, and i0 are the d-axis, q-axis, and zero-sequence stator currents, defined by
ia, ib, and ic are the stator currents flowing from port ~ to port n.
The rotor voltage equations are defined by
and
where:
efd is the field voltage.
e1d, and e1q are the voltages across the d-axis damper winding 1 and q-axis damper winding 1. They are equal to 0.
ψfd, ψ1d, and ψ1q are the magnetic fluxes linking the field circuit, d-axis damper winding 1, and q-axis damper winding 1.
Rfd, R1d, and R1q are the resistances of rotor field circuit, d-axis damper winding 1, and q-axis damper winding 1.
ifd, i1d, and i1q are the currents flowing in the field circuit, d-axis damper winding 1, and q-axis damper winding 1.
The saturation equations are defined by
(if saturation is disabled),
(if saturation is enabled),
and
where:
ψad is the d-axis air-gap or mutual flux linkage.
ψaq is the q-axis air-gap or mutual flux linkage.
ψat is the air-gap flux linkage.
Ks is the saturation factor.
Ladu is the unsaturated mutual inductance of the stator d-axis.
Lad is the mutual inductance of the stator d-axis.
The saturation factor function, f, is calculated from the per-unit open-circuit lookup table as:
and
where:
Vag is the per-unit air-gap voltage.
In per-unit,
and
can be rearranged to
The stator flux linkage equations are defined by
and
where:
Ll is the stator leakage inductance.
Lad and Laq are the mutual inductances of the stator d-axis and q-axis.
The rotor flux linkage equations are defined by
and
where:
Lffd is the self-inductance of the rotor field circuit.
L11d is the self-inductance of the d-axis damper winding 1.
L11q is the self-inductance of the q-axis damper winding 1.
Lf1d is the roto field circuit and d-axis damper winding 1 mutual inductance.
The inductances are defined by these equations:
The inductance equations assume that per-unit mutual inductance L12q = Laq, that is, the stator and rotor currents in the q-axis all link a single mutual flux represented by Laq.
The rotor torque is defined by
You can perform plotting and display actions using the Electrical menu on the block context menu.
Right-click the block and, from the Electrical menu, select one of these options:
Display Base Values — Displays the machine per-unit base values in the MATLAB® Command Window.
Display Associated Base Values — Displays associated per-unit base values in the MATLAB Command Window.
Display Associated Initial Conditions — Displays associated initial conditions in the MATLAB Command Window.
Plot Open-Circuit Saturation (pu) — Plots air-gap voltage, Vag, versus field current, ifd, both measured in per-unit, in a MATLAB figure window. The plot contains three traces:
Unsaturated — Stator d-axis mutual inductance (unsaturated), Ladu you specify
Saturated — Per-unit open-circuit lookup table (Vag versus ifd) you specify
Derived — Open-circuit lookup table (per-unit) derived from the Per-unit open-circuit lookup table (Vag versus ifd) you specify. This data is used to calculate the saturation factor,Ks, versus magnetic flux linkage, ψat, characteristic.
Plot Saturation Factor (pu) — Plots saturation factor,Ks, versus magnetic flux linkage, ψat, both measured in per-unit, in a MATLAB figure window using the present machine parameters. This parameter is derived from other parameters that you specify:
Stator d-axis mutual inductance (unsaturated), Ladu
Per-unit field current saturation data, ifd
Per-unit air-gap voltage saturation data, Vag
The Variables settings allow you to specify the priority and initial target values for block variables before simulation. For more information, see Set Priority and Initial Target for Block Variables.
For this block, the Variables settings are visible only if, in the
Initial Conditions settings, the Initialization
option parameter is set to Set targets for rotor angle and
Park's transform variables
.
[1] Kundur, P. Power System Stability and Control. New York, NY: McGraw Hill, 1993.
[2] Lyshevski, S. E. Electromechanical Systems, Electric Machines and Applied Mechatronics. Boca Raton, FL: CRC Press, 1999.
Synchronous Machine Field Circuit | Synchronous Machine Measurement | Synchronous Machine Model 1.0 | Synchronous Machine Model 2.1 | Synchronous Machine Round Rotor