Vectorized Synchronous Machine GENROU

R2026b

Vectorized GENROU synchronous generator model

Since R2026b

• 

Vectorized Synchronous Machine GENROU block

To add a block to a model, double-click the canvas and start typing the block name. Then, select the block from the list.

Libraries:
Simscape / Electrical / Electromechanical / Synchronous / Vectorized Synchronous Machines

Description

The Vectorized Synchronous Machine GENROU block models a GENROU synchronous generator that supports vectorized inputs and parameters and allows you to simulate multiple machines by using a single block.

To model a single GENROU synchronous machine, consider using the Synchronous Machine GENROU block instead.

GENROU synchronous generator models are widely used in power plant model verification and power system transient stability studies.

This diagram shows the equivalent circuit for the GENROU synchronous generator model.

Equations

The block expresses the synchronous machine equations with respect to a rotating reference frame,

$${\theta }_e(t)=N{\theta }_r(t),$$

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. This equation defines the Park transformation:

$$P_s=\frac{2}{3}\left[\begin{array}{ccc}\cos {\theta }_e& \cos ({\theta }_e-\frac{2\pi }{3})& \cos ({\theta }_e+\frac{2\pi }{3})\\ -\sin {\theta }_e& -\sin ({\theta }_e-\frac{2\pi }{3})& -\sin ({\theta }_e+\frac{2\pi }{3})\\ \frac{1}{2}& \frac{1}{2}& \frac{1}{2}\end{array}\right].$$

The block uses the Park transformation to define the per-unit synchronous machine equations. These equations define the d-axis and q-axis stator voltages:

$$\begin{array}{l}{v}_d=-{\omega }_r{{\psi }^{″}}_q-R_a{i}_d+{\omega }_r{X^{″}}_d{i}_q-\frac{{X^{″}}_d}{{\omega }_{base}}\frac{d{i}_d}{dt}\\ v_q={\omega }_r{{\psi }^{″}}_q-R_a{i}_q-{\omega }_r{X^{″}}_d{i}_d-\frac{{X^{″}}_d}{{\omega }_{base}}\frac{d{i}_q}{dt}\end{array}$$

In these equations:

• ω_r is the per-unit rotor rotational speed.

• ω_base is the electrical base speed.

• $${{\psi }^{″}}_d=\frac{{X^{\prime }}_d-{X^{″}}_d}{{X^{\prime }}_d-X_l}{{\psi }^{\prime }}_d+\frac{{X^{″}}_d-X_l}{{X^{\prime }}_d-X_l}{e^{\prime }}_q$$ is the d-axis flux linkage behind the subtransient reactance X^''_d.

• $${{\psi }^{″}}_q=-\frac{{X^{\prime }}_q-{X^{″}}_q}{{X^{\prime }}_q-X_l}{{\psi }^{\prime }}_q+\frac{{X^{″}}_q-X_l}{{X^{\prime }}_q-X_l}{e^{\prime }}_d$$ is the q-axis flux linkage behind the subtransient reactance X^''_q.

• ψ^'_d is the d-axis flux linkage behind the transient reactance X^'_d.

• ψ^'_q is the q-axis flux linkage behind the transient reactance X^'_q.

• R_a is the stator resistance.

• i_d is the d-axis stator current.

• i_q is the q-axis stator current.

• X^'_d is the d-axis transient reactance.

• X^''_d is the d-axis subtransient reactance.

• X^'_q is the q-axis transient reactance.

• X^''_q is the q-axis subtransient reactance.

• X_l is the stator leakage reactance.

• e^'_d is the d-axis voltage behind the transient reactance.

• e^'_q is the q-axis voltage behind the transient reactance.

These equations define the voltages behind the transient reactances:

$$\begin{array}{l}\frac{d{{\psi }^{\prime }}_d}{dt}=\frac{1}{{T^{″}}_{d0}}\left[{e^{\prime }}_q-{{\psi }^{\prime }}_d-\left({X^{\prime }}_d-X_l\right)i_d\right]\\ \frac{d{{\psi }^{\prime }}_q}{dt}=\frac{1}{{T^{″}}_{q0}}\left[{e^{\prime }}_d-{{\psi }^{\prime }}_q-\left({X^{\prime }}_q-X_l\right)i_q\right]\\ \frac{d{e^{\prime }}_d}{dt}=\frac{1}{{T^{\prime }}_{q0}}\left\{\left[i_q-\frac{{X^{\prime }}_q-{X^{″}}_q}{{\left({X^{\prime }}_q-X_l\right)}^2}\left({e^{\prime }}_d-{{\psi }^{\prime }}_q+\left({X^{\prime }}_q-X_l\right)i_q\right)\right]\left(X_q-{X^{\prime }}_q\right)+S\frac{X_q-X_l}{X_d-X_l}{{\psi }^{″}}_q-{e^{\prime }}_d\right\}\\ \frac{d{e^{\prime }}_q}{dt}=\frac{1}{{T^{\prime }}_{d0}}\left\{E_{fd}-\left[i_d+\frac{{X^{\prime }}_d-{X^{″}}_d}{{\left({X^{\prime }}_d-X_l\right)}^2}\left({e^{\prime }}_q-{{\psi }^{\prime }}_d+\left({X^{\prime }}_d-X_l\right)i_d\right)\right]\left(X_d-{X^{\prime }}_d\right)-S{{\psi }^{″}}_d-{e^{\prime }}_q\right\}\end{array}$$

In these equations:

• T^'_d0 is the d-axis transient open-circuit time constant.

• T^'_q0 is the q-axis transient open-circuit time constant.

• T^''_d0 is the d-axis subtransient open-circuit time constant.

• T^''_q0 is the q-axis subtransient open-circuit time constant.

• $$S=\frac{f\left({\psi }_{at}\right)}{{\psi }_{at}}$$ is the saturation factor:

  • $${\psi }_{at}=\sqrt{{{\psi }^{″}}_d{}^2+{{\psi }^{″}}_q{}^2}$$ is the air-gap flux linkage.

  • If you set the Magnetic saturation representation parameter to None, then the saturation factor is equal to zero.

  • If you set the Magnetic saturation representation parameter to Quadratic, Scaled quadratic, or Exponential, then the block calculates the saturation factor function, f, using the values of the Saturation factor S10 and Saturation factor S12 parameters.

  • If you set the Magnetic saturation representation parameter to Open-circuit lookup table, then the block calculates the saturation factor function using the values of the Per-unit air-gap voltage saturation data, Vag and Non-reciprocal per-unit field current saturation data, Ifd parameters.

This equation defines the per-unit field current in a non-reciprocal per-unit system:

$$I_{fd}=\left\{i_d+\frac{{X^{\prime }}_d-{X^{″}}_d}{{\left({X^{\prime }}_d-X_l\right)}^2}\left[{e^{\prime }}_q-{{\psi }^{\prime }}_d-\left({X^{\prime }}_d-X_l\right)i_d\right]\right\}\left(X_d-{X^{\prime }}_d\right)+S{{\psi }^{″}}_d+{e^{\prime }}_q.$$

This equation defines the rotor torque:

$$T{}_e={{\psi }^{″}}_d{i}_q-{{\psi }^{″}}_q{i}_d.$$

Plotting and Display Options

You can perform these plotting and display actions by clicking the button next to the associated parameter in the Utilities section:

• Base values — Display the machine per-unit base values in the MATLAB^® Command Window.

• Associated initial conditions — Display associated initial conditions in the MATLAB Command Window.

• Open-circuit saturation (pu) — Plot air-gap voltage, V_ag, versus field current, i_fd, both of which are per-unit measurements, in a MATLAB figure window. The plot contains an unsaturated and a saturated trace.

• Saturation factor (pu) — Plot the saturation factor, f, versus magnetic flux linkage, ψ_at, both of which are per-unit measurements, in a MATLAB figure window using the machine parameters. If you set the Magnetic saturation representation parameter to Quadratic, Scaled quadratic, or Exponential, then the block derives the saturation factor function from the value of the Saturation factor, S10 and Saturation factor, S12 parameters.

  If you set the Magnetic saturation representation to Open-circuit lookup table, then the block derives the saturation factor function from the value of the Non-reciprocal per-unit field current saturation data, Ifd and Per-unit air-gap voltage saturation data, Vag parameters.

• Power capability curves — Plot active power P versus reactive power Q, both of which are per-unit measurements, in a MATLAB figure window. The plot can contain multiple traces. Each trace corresponds to a maximum field circuit voltage measured per unit (non-reciprocal per-unit system). The plot shows the maximum reactive power that the generator produces when operating with a lagging power factor and the minimum reactive power that the generator absorbs when operating with a leading power factor.

Ports

Input

Efd — Per-unit field voltage

vector

Per-unit field voltage, specified as a vector.

Pm — Per-unit mechanical power

vector

Per-unit mechanical power, specified as a vector.

Dependencies

To enable this port, in the Mechanical Parameters tab, set the Input type parameter to Power Pm (per unit).

Omega — Per-unit machine speed

vector

Per-unit machine speed, specified as a vector.

Dependencies

To enable this port, in the Mechanical Parameters tab, set the Input type parameter to Speed Omega (per unit).

vABC — Three-phase voltages at stator terminals, V

vector

Three-phase voltages at the stator terminals, in volts, specified as a vector.

Output

pu — Machine per-unit measurements

vector

Machine per-unit measurements, returned as a vector with these elements:

• Stator terminal voltage, Vt

• Stator terminal current, It

• Rotor velocity, wr

• Field current (field circuit base), Ifd

• Electrical torque, Te

• Real power, P

• Reactive power, Q

• Rotor electrical angle, thetae

To connect to this port, use the Synchronous Machine Measurement block.

iABC — Three-phase currents at stator terminals, A

vector

Three-phase currents at the stator terminals, in ampere, returned as a vector.

Parameters

To edit block parameters interactively, use the Property Inspector. From the Simulink^® Toolstrip, on the Simulation tab, in the Prepare gallery, select Property Inspector.

Main

Number of machines — Number of machines to model

3 (default) | positive scalar

Number of machines to model using a single block.

Rated apparent power (VA) — Rated apparent power

[247.5e6 192e6 128e6] (default) | vector of positive scalars

Rated apparent powers of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Rated voltage (V) — RMS rated line-line voltage

[16.5e3 18e3 13.8e3] (default) | vector of positive scalars

Root mean square (RMS) of the rated line-line voltages of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Rated electrical frequency (Hz) — Nominal electrical frequency

[60 60 60] (default) | vector of positive scalars

Nominal electrical frequencies at which the block quotes the rated apparent power of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Number of pole pairs — Number of machine pole pairs

[20 1 1] (default) | vector of positive scalars

Number of pole pairs in each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Rotor angle definition — Reference point for rotor angle measurement

Angle between the a-phase magnetic axis and the d-axis (default) | Angle between the a-phase magnetic axis and the q-axis

Reference point for the rotor angle measurement.

• Angle between the a-phase magnetic axis and the d-axis — The rotor d-axis and stator a-phase magnetic axis are aligned when the rotor angle is zero.

• Angle between the a-phase magnetic axis and the q-axis — The rotor q-axis and stator a-phase magnetic axis are aligned when the rotor angle is zero.

Sample time (s) — Time between consecutive block executions

5e-5 (default) | nonnegative scalar

Time between consecutive block executions. During execution, the block produces outputs and, if appropriate, updates its internal state. For more information, see What Is Sample Time? and Specify Sample Time.

For discrete-time operation, specify a positive scalar. For continuous-time operation, specify 0.

If this block is in a masked subsystem, or other variant subsystem that allows you to switch between continuous operation and discrete operation, promote the sample time parameter. Promoting the sample time parameter ensures correct switching between the continuous and discrete implementations of the block. For more information, see Promote Block Parameters to a Mask.

Impedances (per unit)

Stator resistance, Ra — Resistances of stator

[0.01 0.01 0.01] (default) | vector of positive scalars

Resistances of the stator of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Stator leakage reactance, Xl — Leakage reactances of stator

[0.0832 0.1 0.095] (default) | vector of positive scalars

Leakage reactance of the stator of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

d-axis synchronous reactance, Xd — Direct-axis synchronous reactances

[0.3614 1.7199 1.68] (default) | vector of positive scalars

Direct-axis synchronous reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

q-axis synchronous reactance, Xq — Quadrature-axis synchronous reactances

[0.2398 1.6598 1.61] (default) | vector of positive scalars

Quadrature-axis synchronous reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

d-axis transient reactance, Xd' — Direct-axis transient reactances

[0.1505 0.23 0.2321] (default) | vector of positive scalars

Direct-axis transient reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

q-axis transient reactance, Xq' — Quadrature-axis transient reactances

[0.1854 0.378 0.32] (default) | vector of positive scalars

Quadrature-axis transient reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

d-axis subtransient reactance, Xd'' — Direct-axis subtransient reactances

[0.0990 0.1728 0.128] (default) | vector of positive scalars

Direct-axis subtransient reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

q-axis subtransient reactance, Xq'' — Quadrature-axis subtransient reactances

[0.1114 0.192 0.192] (default) | vector of positive scalars

Quadrature-axis subtransient reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Time Constants (s)

d-axis transient open-circuit, Td0' — Direct-axis transient open-circuit time constants

[8.96 6 5.89] (default) | vector of positive scalars

Direct-axis transient open-circuit time constant of each machine, in seconds. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

d-axis subtransient open-circuit, Td0'' — Direct-axis subtransient open-circuit time constants

[0.05 0.05 0.05] (default) | vector of positive scalars

Direct-axis subtransient open-circuit time constant of each machine, in seconds. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

q-axis transient open-circuit, Tq0' — Quadrature-axis transient open-circuit time constants

[0.48 0.535 0.6] (default) | vector of positive scalars

Quadrature-axis transient open-circuit time constant of each machine, in seconds. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

q-axis subtransient open-circuit, Tq0'' — Quadrature-axis subtransient open-circuit time constants

[0.09 0.04 0.04] (default) | vector of positive scalars

Quadrature-axis subtransient open-circuit time constant of each machine, in seconds. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Saturation

Magnetic saturation representation — Magnetic saturation representation

None (default) | Quadratic | Scaled quadratic | Exponential | Open-circuit lookup table

Option to represent the magnetic saturation of the block:

• None — The block does not model the magnetic saturation.

• Quadratic — The block models the magnetic saturation by using a quadratic function.

• Scaled quadratic — The block models the magnetic saturation by using a scaled function.

• Exponential — The block models the magnetic saturation by using an exponential function.

• Open-circuit lookup table — The block models the magnetic saturation by using per-unit non-reciprocal field current and per-unit air-gap voltage saturation data.

Saturation factor, S10 — Saturation factor, S10

0.005 (default) | positive scalar

Saturation factor that corresponds to 1.0 per-unit terminal voltage. The value of this parameter must be less than the value of the Saturation factor, S12 parameter.

Dependencies

To enable this parameter, set Magnetic saturation representation to Quadratic, Scaled quadratic, or Exponential.

Saturation factor, S12 — Saturation factor, S12

0.08 (default) | positive scalar

Saturation factor that corresponds to 1.2 per-unit terminal voltage. The value of this parameter must be greater than the value of the Saturation factor, S10 parameter.

Dependencies

To enable this parameter, set Magnetic saturation representation to Quadratic, Scaled quadratic, or Exponential.

Non-reciprocal per-unit field current saturation data, Ifd — Non-reciprocal per-unit field current saturation data

[ 0.00 0.80 1.26 2.29 2.97 ] (default) | vector of nonnegative values

Nonreciprocal field current, I_fd, data that populates the lookup table for air-gap voltage, V_ag, versus field current. The value of this parameter must be a vector of at least five elements.

Dependencies

To enable this parameter, set Magnetic saturation representation to Open-circuit lookup table.

Per-unit air-gap voltage saturation data, Vag — Per-unit air-gap voltage saturation data

[ 0.00 0.80 1.08 1.31 1.40 ] (default) | vector of nonnegative values

Air-gap voltage, V_ag, data that populates the lookup table for air-gap voltage versus field current, i_fd. The value of this parameter must be a vector of at least five elements.

Dependencies

To enable this parameter, set Magnetic saturation representation to Open-circuit lookup table.

Mechanical Parameters

Input type — Type of mechanical input

Power Pm (per unit) (default) | Speed Omega (per unit)

Option to specify the type of the mechanical input:

• Power Pm (per unit) — The block accepts the per-unit mechanical power as input. This value exposes the Pm port.

• Speed Omega (per unit) — The block accepts the per-unit machine speed as input. This value exposes the Omega port.

Inertia constant, H (s) — Inertia constants

[9.5515 3.3333 2.3516] (default) | vector of positive scalars

Inertia constant of each machine, in seconds. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Input type to Power Pm (per unit).

Damping coefficient, D (per unit) — Damping coefficients

[0 0 0] (default) | vector of nonnegative scalars

Per-unit damping coefficient of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Input type to Power Pm (per unit).

Initial Conditions

Specify initialization by — Option to specify values at start of simulation

Electrical power and voltage (default) | Mechanical and magnetic states

Option to specify values for certain parameters and variables at the start of simulation:

• Electrical power and voltage — Specify the magnitude and angle of the terminal voltages, and the active and reactive generated powers.

• Mechanical and magnetic states — Specify the mechanical and magnetic states of each machine, such as the rotor electrical angle and frequency, the initial dq-axes voltages behind the transient reactances, the initial direct-quadrature (dq) axes currents, and the initial flux linkages behind the dq-axes transient reactances.

Terminal voltage magnitude (V) — Magnitudes of terminal voltages

[16.5e3 18e3 13.8e3] (default) | vector of positive scalars

Magnitude of the terminal voltage of each machine, in volts. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Electrical power and voltage.

Terminal voltage angle (deg) — Angles of terminal voltages

[0 0 0] (default) | vector of scalars

Angle of the terminal voltage of each machine, in degrees. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Electrical power and voltage.

Active power generated — Generated active powers

[50e6 50e6 50e6] (default) | vector of positive scalars

Active power that each machine generates. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Electrical power and voltage.

Reactive power generated — Generated reactive powers

[0 0 0] (default) | vector of scalars

Reactive power that each machine generates. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Electrical power and voltage.

Rotor electrical angle (rad) — Rotor electrical angles

[0 0 0] (default) | vector of scalars

Rotor electrical angle of each machine, in rad. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Mechanical and magnetic states.

Rotor electrical frequency (Hz) — Rotor electrical frequencies

[60 60 60] (default) | vector of positive scalars

Rotor electrical frequency of each machine, in Hz. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Mechanical and magnetic states.

Psid' — Initial d-axis flux linkages behind transient reactance

[0 0 0] (default) | vector

Initial direct-axis flux linkage behind the transient reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Mechanical and magnetic states.

Psiq' — Initial q-axis flux linkages behind transient reactance

[0 0 0] (default) | vector

Initial quadrature-axis flux linkage behind the transient reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Mechanical and magnetic states.

Ed' — Initial d-axis voltages behind transient reactance

[0 0 0] (default) | vector

Initial direct-axis voltage behind the transient reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Mechanical and magnetic states.

Eq' — Initial q-axis voltages behind transient reactance

[0 0 0] (default) | vector

Initial quadrature-axis voltage behind the transient reactance of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Mechanical and magnetic states.

Id — Initial d-axis stator currents

[0 0 0] (default) | vector

Initial direct-axis stator current of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Mechanical and magnetic states.

Iq — Initial q-axis stator currents

[0 0 0] (default) | vector

Initial quadrature-axis stator current of each machine. The value of this parameter must be a vector of the same size as the value of the Number of machines parameter.

Dependencies

To enable this parameter, set Specify initialization by to Mechanical and magnetic states.

Utilities

Choose a machine — Choose machine for which to display and plot values and conditions

1 (default) | positive scalar

Choose the machine for which you want to display and plot specific values and conditions. The value of this parameter must be equal to or greater than 1 and equal to or less than the Number of machines parameter value.

Display base values — Display machine per-unit base values

button

Display the per-unit base values of the machine you specified in the Choose a machine parameter in the MATLAB Command Window.

Display associated initial conditions — Display associated initial conditions

button

Display the associated initial conditions of the machine you specified in the Choose a machine parameter in the MATLAB Command Window.

Plot open-circuit saturation (pu) — Plot air-gap voltage against field current

button

Plot the per-unit air-gap voltage, V_ag, against the per-unit field current, i_fd, of the machine you specified in the Choose a machine parameter in a MATLAB figure window.

Plot saturation factor (pu) — Plot saturation factor against magnetic flux linkage

button

Plot the per-unit saturation factor, f, against the per-unit magnetic flux linkage, ψ_at, of the machine you specified in the Choose a machine parameter in a MATLAB figure window.

Plot power capability curves — Plot active power against reactive power

button

Plot the per-unit active power P against the per-unit reactive power Q of the machine you specified in the Choose a machine parameter in a MATLAB figure window.

References

[1] Sauer, Peter W., and M. A. Pai. "Power System Dynamics and Stability". Prentice Hall, 1998.

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using Simulink® Coder™.

Version History

Introduced in R2026b