# Saturable Transformer

Implement two- or three-winding saturable transformer

**Libraries:**

Simscape /
Electrical /
Specialized Power Systems /
Power Grid Elements

## Description

The Saturable Transformer block model shown consists of three coupled windings wound on the same core.

The model takes into account the winding resistances (R1 R2 R3) and the leakage inductances (L1 L2 L3) as well as the magnetizing characteristics of the core, which is modeled by a resistance Rm simulating the core active losses and a saturable inductance Lsat.

You can choose one of the following two options for the modeling of the nonlinear flux-current characteristic

Model saturation without hysteresis. The total iron losses (eddy current + hysteresis) are modeled by a linear resistance, Rm.

Model hysteresis and saturation. Specification of the hysteresis is done by means of the

**Hysteresis Design Tool**of the Powergui block. The eddy current losses in the core are modeled by a linear resistance, Rm.**Note**Modeling the hysteresis requires additional computation load and therefore slows down the simulation. The hysteresis model should be reserved for specific applications where this phenomenon is important.

### Saturation Characteristic Without Hysteresis

When the hysteresis is not modeled, the saturation characteristic of the Saturable Transformer block is defined by a piecewise linear relationship between the flux and the magnetization current.

Therefore, if you want to specify a residual flux, phi0, the second point of the saturation characteristic should correspond to a null current, as shown in the figure (b).

The saturation characteristic is entered as (i, phi) pair values in per units, starting with pair (0, 0). The software converts the vector of fluxes Φpu and the vector of currents Ipu into standard units to be used in the saturation model of the Saturable Transformer block:

Φ = Φ_{pu}Φ_{base}I =
I_{pu}I_{base}, | (1) |

where the base flux linkage (Φ_{base}) and base current
(*I*_{base}) are the peak values obtained at nominal
voltage power and frequency:

$$\begin{array}{c}{I}_{\text{base}}=\frac{Pn}{{V}_{1}}\sqrt{2}\\ {\Phi}_{\text{base}}=\frac{{V}_{1}}{2\pi {f}_{n}}\sqrt{2}.\end{array}$$

The base flux is defined as the peak value of the sinusoidal flux (in webers) when
winding 1 is connected to a 1 pu sinusoidal voltage source (nominal voltage). The
Φ_{base} value defined above represents the base flux linkage (in
volt-seconds). It is related to the base flux by the following equation:

Φ_{base} = Base flux × number of turns of winding
1. | (2) |

When they are expressed in pu, the flux and the flux linkage have the same value.

### Saturation Characteristic with Hysteresis

The magnetizing current I is computed from the flux Φ obtained by integrating voltage across the magnetizing branch. The static model of hysteresis defines the relation between flux and the magnetization current evaluated in DC, when the eddy current losses are not present.

The hysteresis model is based on a semi-empirical characteristic, using an arctangent
analytical expression Φ(I) and its inverse I(Φ) to represent the operating point
trajectories. The analytical expression parameters are obtained by curve fitting empirical
data defining the major loop and the single-valued saturation characteristic. The **Hysteresis design tool** of the Powergui block is used
to fit the hysteresis major loop of a particular core type to basic parameters. These
parameters are defined by the remanent flux (Φr), the coercive current (Ic), and the slope
(dΦ/dI) at (0, Ic) point as shown in the next figure.

The major loop half cycle is defined by a series of N equidistant points connected by
line segments. The value of N is defined in the **Hysteresis design
tool** of the Powergui block. Using N = 256 yields a smooth curve
and usually gives satisfactory results.

The single-valued saturation characteristic is defined by a set of current-flux pairs defining a saturation curve which should be asymptotic to the air core inductance Ls.

The main characteristics of the hysteresis model are summarized below:

A symmetrical variation of the flux produces a symmetrical current variation between -Imax and +Imax, resulting in a symmetrical hysteresis loop whose shape and area depend on the value of Φmax. The major loop is produced when Φmax is equal to the saturation flux (Φs). Beyond that point the characteristic reduces to a single-valued saturation characteristic.

In transient conditions, an oscillating magnetizing current produces minor asymmetrical loops, as shown in the next figure, and all points of operation are assumed to be within the major loop. Loops once closed have no more influence on the subsequent evolution.

The trajectory starts from the initial (or residual) flux point, which must lie on the vertical axis inside the major loop. You can specify this initial flux value phi0, or it is automatically adjusted so that the simulation starts in steady state.

### The Per Unit Conversion

In order to comply with industry practice, the block allows you to specify the resistance and inductance of the windings in per unit (pu). The values are based on the transformer rated power Pn in VA, nominal frequency fn in Hz, and nominal voltage Vn, in Vrms, of the corresponding winding. For each winding the per unit resistance and inductance are defined as

$$\begin{array}{c}R(\text{p}\text{.u}\text{.})=\frac{R(\Omega )}{{R}_{\text{base}}}\\ L(\text{p}\text{.u}\text{.})=\frac{L(H)}{{L}_{\text{base}}}.\end{array}$$

The base resistance and base inductance used for each winding are

$$\begin{array}{c}{R}_{\text{base}}=\frac{{V}_{n}^{2}}{Pn}\\ {L}_{\text{base}}=\frac{{R}_{\text{base}}}{2\pi {f}_{n}}.\end{array}$$

For the magnetization resistance Rm, the pu values are based on the transformer rated power and on the nominal voltage of winding 1.

The default parameters of winding 1 specified in the dialog box section give the following base values:

$$\begin{array}{c}{R}_{\text{base}}=\frac{{\left(735\cdot {10}^{3}/\sqrt{3}\right)}^{2}}{250\cdot {10}^{6}}=720.3\Omega \\ {L}_{\text{base}}=\frac{720.3}{2\pi \cdot 60}=1.91H.\end{array}$$

For example, if winding 1 parameters are R1 = 1.44 Ω and L1 = 0.1528 H, the corresponding values to enter in the dialog box are

$$\begin{array}{c}{R}_{1}=\frac{1.44\Omega}{720.3\Omega}=0.002\text{p}\text{.u}\text{.}\\ {L}_{1}=\frac{0.1528H}{1.91H}=0.08\text{p}\text{.u}\text{.}\end{array}$$

### Examples

The `power_xfosaturable`

example illustrates the energization of one
phase of a three-phase 450 MVA, 500/230 kV transformer on a 3000 MVA source. The transformer
parameters are

Nominal power and frequency | Pn = 150e6 VA | fn = 60 Hz | |

Winding 1 parameters (primary) | V1 = 500e3 Vrms/sqrt(3) | R1 = 0.002 pu | L1 = 0.08 pu |

Winding 2 parameters (secondary) | V2 = 230e3 Vrms/sqrt(3) | R2 = 0.002 pu | L2 = 0.08 pu |

Saturation characteristic | [0 0; 0.0 1.2; 1.0 1.52] | ||

Core loss resistance and initial flux | Rm = 500 pu | phi0 = 0.8 pu |

Simulation of this circuit illustrates the saturation effect on the transformer current and voltage.

As the source is resonant at the fourth harmonic, you can observe a high fourth- harmonic content in the secondary voltage. In this circuit, the flux is calculated in two ways:

By integrating the secondary voltage

By using the Multimeter block

### Assumptions and Limitations

Windings can be left floating (that is, not connected by an impedance to the rest of the circuit). However, the floating winding is connected internally to the main circuit through a resistor. This invisible connection does not affect voltage and current measurements.

## Ports

### Conserving

## Parameters

## References

[1] Casoria, S., P. Brunelle, and G. Sybille, “Hysteresis Modeling in the MATLAB/Power System Blockset,” Electrimacs 2002, École de technologie supérieure, Montreal, 2002.

[2] Frame, J.G., N. Mohan, and Tsu-huei Liu, “Hysteresis modeling in an Electro-Magnetic Transients Program,” presented at the IEEE PES winter meeting, New York, January 31 to February 5, 1982.

## Extended Capabilities

## Version History

**Introduced before R2006a**