Fault Analysis Using Z Bus
Fault Analysis Using Z Bus
1.0 Introduction
(4)
2
(5)
(6)
3
(9)
4
Now eq. (12) is useful for computing fault
currents in the circuits. Consider Fig. 1.
Fig. 1
(16)
5
We can use eq. (16) to get the fault current
in the circuits. These values provide us with
the appropriate information for selecting the
circuit breakers in the lines.
6
The last observation can be utilized in an
effective fashion when performing fault
analysis. Let’s assume that we want to
compute the short circuit currents for a fault
at only one bus k. So we just want to get the
kth column of Z-bus, but we do not need the
entire Z-bus.
7
1. Column of Z: We can just consider a
single column of Z, instead of the entire
matrix. Call it Zk; it is a column vector.
The right hand side of (19) will just be a
column of I. Call it Ik. It will also be just
a vector and will contain zeros in every
row except for row k. The resulting
relation is:
(21)
2. LU-Decomposition: If you have taken
EE 456 or a linear algebra course in
math, then you are familiar with LU
decomposition. LU decomposition
provides a way to solve for the vector x
in the matrix relation
(22)
where A is a n×n square matrix, x is an
unknown n×1 vector, and b is an known
n×1 vector. The advantage to LU
decomposition is it does not require
inverting the matrix A. The basis of LU-
decomposition is that we may factor A
into a matrix product LU, i.e.,
8
(23)
where L is a lower diagonal n×n matrix
of the form
(24)
and U is an upper diagonal n×n matrix of
the form
(25)
Substitution of eq. (23) in (22) yields
(26)
Defining
(27)
provides that eq. (25) becomes:
(28)
If we have L and U, then (28) is easily
solved for w (without inverting L) using
forward substitution, and then (27) can
be easily solved for x (without inverting
U) using backwards substitution. More
details on LU decomposition may be
found in the notes called
“LUDecomposition.doc.”
9
Now observe that eq. (21) and (22) are in the
same form. Therefore we want to solve the
following equations in the order they are
given:
(29)
(30)
Factorization of a matrix Y can be done efficiently and easily using the matlab
command:
[L,U]=lu(Y)
Then it is easy to find w by hand using forward substitution from:
Alternatively, the manual steps of LU decomposition can be performed per the notes
from “LU Decomposition.”
10
a. Construct the Y-bus for this network (should be
a 4×4 matrix).
b. Consider that there is a three-phase
(symmetrical) fault at bus 2.
i. Use LU decomposition to obtain the 2nd
column of the Z-bus.
ii. Compute the subtransient fault current.
iii. Use eq. (12) to find the voltages during the
fault.
iv. Use eq. (17) to find the subtransient currents
in lines 3-2, 1-2, and 4-2.
3. See notes called “Examples.”
4. See notes called “Examples.”
11