# Simplified Synchronous Machine

Model the dynamics of simplified three-phase synchronous machine

**Libraries:**

Simscape /
Electrical /
Electromechanical /
Synchronous

## Description

The Simplified Synchronous Machine block models both the electrical and mechanical characteristics of a simple synchronous machine.

The electrical system for each phase consists of a voltage source in series with an RL impedance, which implements the internal impedance of the machine. The value of R can be zero but the value of L must be positive.

The Simplified Synchronous Machine block implements the mechanical system described by

$$\begin{array}{c}\Delta \omega (t)=\frac{1}{2H}{\displaystyle \underset{0}{\overset{t}{\int}}\left(Tm-Te\right)\text{\hspace{0.17em}}}-Kd\Delta \omega (t)dt\\ \omega (t)=\Delta \omega (t)+{\omega}_{0},\end{array}$$

where

Δω = Speed variation with respect to speed of
operationH = constant of
inertiaTm = mechanical
torqueTe = electromagnetic
torqueKd = damping factor representing the effect of
damper windingsω(t) = mechanical
speed of the rotorω_{0} = speed of
operation (1 p.u.) | (1) |

Although the parameters can be entered in either SI units or per unit in the dialog box, the internal calculations are done in per unit. The following block diagram illustrates how the mechanical part of the model is implemented. The model computes a deviation with respect to the speed of operation; not the absolute speed itself.

The Kd damping coefficient simulates the effect of damper windings normally used in
synchronous machines. When the machine is connected to an infinite network (zero impedance),
the variation of machine power angle delta (δ) resulting from a change of mechanical power
(*P*_{m}) can be approximated by the following
second-order transfer function:

$$\frac{\delta}{{P}_{m}}=\frac{{\omega}_{s}/(2H)}{{s}^{2}+2\zeta {\omega}_{n}s+{\omega}_{n}^{2}},$$

where

δ | Power angle delta: angle of internal voltage |

P | Mechanical power in pu |

ω | Frequency of electromechanical oscillations = $$\sqrt{{\omega}_{s}{P}_{\text{max}}/(2H)}$$ in rad/s |

ζ | Damping ratio = $$\left({K}_{d}/4\right)\sqrt{2/\left({\omega}_{s}H{P}_{\text{max}}\right)}$$ |

ω | Electrical frequency in rad/s |

| Maximum power in pu transmitted through reactance X, where V_{t},
E, and X are in pu |

H | Inertia constant(s) |

K | Damping factor (pu_of_torque / pu_of_speed) |

This approximate transfer function, which has been derived by assuming sin(δ) = δ, is valid for small power angles (δ < 30 degrees). It follows from the preceding ζ expression that the Kd value required to obtain a given ζ damping ratio:

$${K}_{d}=4\zeta \sqrt{{\omega}_{s}H{P}_{\text{max}}/2}.$$

### Examples

The `power_simplealt`

example uses the Simplified Synchronous
Machine block to represent a 1000 MVA, 315 kV, 60 Hz equivalent source connected
to an infinite bus (Three-Phase Programmable Voltage Source block). The Simplified
Synchronous Machine (SI Units) block is used as a synchronous generator. The internal
resistance and reactance are set respectively to 0.02 pu (1.9845 Ω) and 0.2 pu (X = 19.845
Ω; L = 0.0526 H). The inertia of the machine is J = 168,870
kg.m^{2}, corresponding to an inertia constant H = 3 s. The
electrical frequency is ω_{s} = 2*π*60 = 377 rad/s. The machine has two
pairs of poles such that its synchronous speed is 2*π*60/2 = 188.5 rad/s or 1800 rpm.

The Load Flow Analyzer
app has been used to initialize the machine to start simulation in steady state with the
machine generating 500 MW. The required internal voltage computed by the load flow is 1.0149
pu. Therefore, an internal voltage E = 315e3*1.0149 = 319,690 Vrms phase-to-phase is
specified in the Constant block connected to the E input. The maximum power that can be
delivered by the machine with a terminal voltage V_{t} = 1.0 pu and an
internal voltage E = 1.0149 pu is P_{max} =
V_{t}*E/X = 1.0149/0.2 = 5.0745 pu.

The damping factor Kd is adjusted to obtain a damping ratio ζ = 0.3. The required Kd value is:

$${K}_{d}=4\zeta \sqrt{{\omega}_{s}H{P}_{\text{max}}/2}=64.3$$

Two Fourier blocks measure the power angle δ. This angle is computed as the difference between the phase angle of phase A internal voltage and the phase angle of phase A terminal voltage.

In this example, a step is performed on the mechanical power applied to the shaft. The machine is initially running in steady state with a mechanical power of 505 MW (mechanical power required for an output electrical power of 500 MW, considering the resistive losses). At t = 0.5 s the mechanical power is suddenly increased to 1000 MW.

Run the example and observe the electromechanical transient on the Scope block displaying the power angle δ in degrees, the machine speed in rpm, and the electrical power in MW.

For an initial electrical power Pe = 500 MW (0.5 pu), the load angle δ is 5.65 degrees, which corresponds to the expected value:

$$Pe=\frac{{V}_{t}E\mathrm{sin}\delta}{X}=\frac{1.0\cdot 1.0149\cdot \mathrm{sin}({5.65}^{\circ})}{0.2}=0.5\text{p}\text{.u}\text{.}$$

As the mechanical power is stepped from 0.5 pu to 1.0 pu, the load angle increases and goes through a series of under damped oscillations (damping ratio ζ = 0.3) before stabilizing to its new value of 11.3 degrees. The frequency of the oscillations is given by:

$${f}_{n}=\frac{1}{2\pi}\sqrt{\frac{{\omega}_{s}{P}_{\text{max}}}{2H}}=2.84\text{Hz}\text{.}$$

## Assumptions and Limitations

The electrical system of the Simplified Synchronous Machine block consists solely of a voltage source behind a synchronous reactance and resistance. All the other self- and magnetizing inductances of the armature, field, and damping windings are neglected. The effect of damper windings is approximated by the damping factor Kd. The three voltage sources and RL impedance branches are Y-connected (three wires or four wires). The load might or might not be balanced.

When you use Simplified Synchronous Machine blocks in discrete systems, you might have to use a small parasitic resistive load, connected at the machine terminals, to avoid numerical oscillations. Large sample times require larger loads. The minimum resistive load is proportional to the sample time. Remember that with a 25 μs time step on a 60 Hz system, the minimum load is approximately 2.5% of the machine nominal power. For example, a 200 MVA simplified synchronous machine in a power system discretized with a 50 μs sample time requires approximately 5% of resistive load or 10 MW. If the sample time is reduced to 20 μs, a resistive load of 4 MW should be sufficient.

## Ports

### Input

### Output

### Conserving

## Parameters

## Extended Capabilities

## Version History

**Introduced before R2006a**