Nodal admittance matrix

Electric power transmission needs optimization in order to determine the necessary real and reactive power flows in a system for a given set of loads, as well as the voltages and currents in the system. Power flow studies are used not only to analyze current power flow situations, but also to plan ahead for anticipated disturbances to the system, such as the loss of a transmission line to maintenance and repairs. The power flow study would determine whether or not the system could continue functioning properly without the transmission line. Only computer simulation allows the complex handling required in power flow analysis because in most realistic situations the system is very complex and extensive and would be impractical to solve by hand. The Y Matrix is a tool in that domain. It provides a method of systematically reducing a complex system to a matrix that can be solved by a computer program. The equations used to construct the Y matrix come from the application of Kirchhoff's current law and Kirchhoff's voltage law to a circuit with steady-state sinusoidal operation. These laws give us that the sum of currents entering a node in the circuit is zero, and the sum of voltages around a closed loop starting and ending at a node is also zero. These principles are applied to all the nodes in a power flow system and thereby determine the elements of the admittance matrix, which represents the admittance relationships between nodes, which then determine the voltages, currents and power flows in the system.

Starting from the single line diagram of a power system, there are three main steps before writing the equations that form the ${\displaystyle Y}$ Matrix. First, the single line diagram is converted to an impedance diagram. Next, all voltage sources are converted to their equivalent current source representations. From here, the impedance diagram is then converted to an admittance diagram. Following these three steps, the admittance matrix can be created in a straightforward manner: For an admittance diagram with ${\displaystyle N}$ buses, the admittance between the bus in consideration, k, and another bus, i, connected to k, can be described by ${\displaystyle y_{ik}=g_{ik}+jb_{ik}}$. The term ${\displaystyle y_{k}}$ should be introduced here; this term accounts for the admittance of linear loads connected to bus ${\displaystyle k}$ as well as the admittance-to-ground at bus ${\displaystyle k}$. The general mathematical expression follows:

${\displaystyle Y_{ij}={\begin{cases}y_{i}+\sum _{k=1,2,\ldots ,N \atop k\neq i}{y_{ik}},&{\mbox{if}}\quad i=j\\-y_{ij},&{\mbox{if}}\quad i\neq j\end{cases}}}$

It is important to note that ${\displaystyle y_{ki}}$ is non-zero only where a physical connection exists between two buses.[1] This consideration is not seen in the following example because each node is connected to both of the other nodes. Each ${\displaystyle y_{ki}}$ defines one element of the ${\displaystyle N\times {N}}$ matrix. From the general case where N is greater than 2, it is desirable to solve these equations as a system, namely through matrix algebra. The general matrix appears as follows: The nodal admittance matrix form: ${\displaystyle Y={\begin{bmatrix}Y_{11}&Y_{12}&\cdots &Y_{1n}\\Y_{21}&Y_{22}&\cdots &Y_{2n}\\\cdots &\cdots &\cdots &\cdots \\Y_{n1}&Y_{n2}&\cdots &Y_{nn}\end{bmatrix}}}$

The admittance diagram of a three bus network.

Once the admittance matrix has been formed, the admittance matrix can be input to solve the matrix form of Ohm's law—the equation ${\displaystyle YV=I}$. In this instance ${\displaystyle V}$ is an ${\displaystyle N\times {1}}$ vector of the voltage at each node and ${\displaystyle I}$ is the ${\displaystyle N\times {1}}$ vector of corresponding currents. In matrix form, Ohm's law is as follows:

${\displaystyle {\begin{pmatrix}Y_{11}&Y_{12}&\cdots &Y_{1N}\\Y_{21}&Y_{22}&\cdots &Y_{2N}\\\vdots &\ddots &\cdots &\vdots \\Y_{NN}&Y_{N2}&\cdots &Y_{NN}\\\end{pmatrix}}{\begin{pmatrix}V_{1}\\V_{2}\\\vdots \\V_{N}\\\end{pmatrix}}={\begin{pmatrix}I_{1}\\I_{2}\\\vdots \\I_{N}\\\end{pmatrix}}}$

To illustrate this process with the admittance matrix of the three bus network in the figure would be:

${\displaystyle Y={\begin{pmatrix}y_{1}+y_{12}+y_{13}&-y_{12}&-y_{13}\\-y_{12}&y_{2}+y_{12}+y_{23}&-y_{23}\\-y_{13}&-y_{23}&y_{3}+y_{13}+y_{23}\\\end{pmatrix}}}$

The Y Matrix diagonal elements ${\displaystyle Y_{11},Y_{22},...,Y_{nn}}$ are called the self-admittances at the nodes, and each equals the sum of all the admittances terminating on the node identified by the repeated subscripts. The other admittances are the mutual admittances of the nodes, and each equals the negative of the sum of all admittances connected directly between the nodes identified by the double subscripts. The admittance matrix ${\displaystyle Y}$ is typically a symmetric matrix as ${\displaystyle y_{mn}=y_{nm}}$. However, extensions of the line model and models of other components may make ${\displaystyle Y}$ asymmetrical.[2] An example being a phase-shifting transformer, which will cause ${\displaystyle Y}$ to become asymmetrical.

For small transmission systems of about less than 10 nodes or buses, the Y matrix can be calculated manually. But for a realistic system with relatively large number of nodes or buses, say 1000 nodes, a computer program for computing Y is more practical to use.

Node Network

To help motivate the importance of using a system of equations in matrix form, see the adjacent figure. Not only does it become impractical to calculate the current vector ${\displaystyle V}$ by hand, it becomes necessary to use computational power to form the admittance matrix itself.

Example:[3]

Two node power flow diagram courtesy of Grainger and Stevenson

To take a look at a generalizable ${\displaystyle 2\times 2}$ matrix, consider the figure of the two-node network. By Kirchhoff's Current Law, it can be shown that :${\displaystyle {\begin{pmatrix}I_{m}\\I_{n}\\\end{pmatrix}}={\begin{pmatrix}1\\-1\\\end{pmatrix}}I_{a}}$ since there are no other currents entering or exiting nodes ${\displaystyle m}$ or ${\displaystyle n}$. The voltage drop across the line can be expressed as :${\displaystyle V_{a}={\begin{pmatrix}1,-1\\\end{pmatrix}}{\begin{pmatrix}V_{m}\\V_{n}\\\end{pmatrix}}}$. Next, use Ohm's law with admittance instead of impedance. Using substitution to get :${\displaystyle {\begin{pmatrix}Y_{a}&-Y_{a}\\-Y_{a}&Y_{a}\end{pmatrix}}{\begin{pmatrix}V_{m}\\V_{n}\\\end{pmatrix}}={\begin{pmatrix}I_{m}\\I_{n}\\\end{pmatrix}}}$. To reintroduce some generality, ${\displaystyle Y_{11}=Y_{mm},Y_{21}=Y_{nm},Y_{12}=Y_{mn},}$ and ${\displaystyle Y_{22}=Y_{nn}}$. Thus, this example can be taken as a first step in understanding how to generally construct an ${\displaystyle N\times N}$ matrix by hand.

• Admittance parameters
• Nodal analysis
• Zbus

• A C/C++ Program and Source Code for Computing Ybus and Zbus Matrices