MATLAB Examples

Transient Analysis of a Linear Circuit

This example shows steady-state and transient simulation of a linear circuit.

H. Le-Huy (Universite Laval, Quebec) and G. Sybille (Hydro-Quebec)

Contents

Description

This circuit is a simplified model of a 230 kV three-phase power system. Only one phase of the transmission system is represented. The equivalent source is modeled by a voltage source (230 kV rms/sqrt(3) or 187.8 kV peak, 60 Hz) in series with its internal impedance (Rs Ls) corresponding to a 3-phase 2000 MVA short circuit level and X/R = 10. (X = 230e3^2/2000e6 = 26.45 ohms or L = 0.0702 H, R = X/10 = 2.645 ohms). The source feeds a RL load through a 150 km transmission line. The line distributed parameters (R = 0.035ohm/km, L = 0.92 mH/km, C = 12.9 nF/km) are modeled by a single pi section (RL1 branch 5.2 ohm; 138 mH and two shunt capacitances C1 and C2 of 0.967 uF). The load (75 MW - 20 Mvar per phase) is modeled by a parallel RLC load block.

A circuit breaker is used to switch the load at the receiving end of the transmission line. The breaker which is initially closed is opened at t = 2 cycles, then it is reclosed at t = 7 cycles. Current and Voltage Measurement blocks provide signals for visualization purpose.

Simulation

1. Simulation using a continuous solver (ode23tb)

Start the simulation and observe line voltage and load current transients during load switching and note that the simulation starts in steady-state. Use the zoom buttons of the oscilloscope to observe the transient voltage at breaker reclosing.

2. Using the Powergui to obtain steady-state phasors and set initial states

Open the Powergui block and select "Steady State Voltage and Currents" to measure the steady-state voltage and current phasors. Using the Powergui select now "Initial States Setting" to obtain the initial state values (voltage across capacitors and current in inductances). Now, reset all the initial states to zero by clicking the "to zero" button and then "Apply" to confirm changes. Restart the simulation and observe transients at simulation starting. Using the same Powergui window, you can also set selected states to specific values.

3. Discretizing your circuit and simulating at fixed steps

The Powergui block can also be used to discretize your circuit and simulate it at fixed steps. Open the Powergui. Select "Discretize electrical model" and specify a sample time of 50e-6 s. The state-space model will now be discretized using trapezoidal fixed step integration. The precision of results is now imposed by the sample time. Restart the simulation and compare simulation results with the continuous integration method. Vary the sample time of the discrete system and note the impact on precision of fast transients.

4. Using the phasor simulation method

You will now use a third simulation technique. The "phasor simulation" method consists to replace the circuit state-space model by a set of algebraic equations evaluated at a fixed frequency and to replace sinusoidal voltage and current sources by phasors (complex numbers). This method allows a fast computation of voltage and current phasors at a selected frequency, disregarding fast transients. It is particularly efficient to study electromechanical transients of generators and motors involving low frequency oscillation modes. Open the Powergui block and select "Phasor simulation". Restart the simulation. Observe that the magnitude of 60 Hz voltage and current is now displayed on the scope. If you double click on the voltage or current measurement block you can choose to output phasor signals in four different formats: Complex, Real/Imag, Magnitude/Angle (in degrees), or just Magnitude (default value). Notice that you cannot send a complex signal to an oscilloscope.