Concept:
Bus Admittance Matrix (YBus):
Consider a small power system network consisting of two generating stations, three transmission lines, one load, and a static capacitor connected to load bus 3.
Assumed that the network is symmetrical and operating under the balance conditions.
The node voltage equation of the system can be written as,
\({I_1} = \left( {{y_{12}} + {y_{31}}} \right){V_1} - {y_{12}}{V_2} - {y_{31}}{V_3}\)
\({I_2} = - {y_{12}}{V_1} + \left( {{y_{12}} + {y_{23}}} \right){V_2} - {y_{23}}{V_3}\)
\( - {I_3} = - {y_{31}}{V_1} - {y_{23}}{V_2} + \left( {{y_{31}} + {y_{23}} + {y_{30}}} \right){V_3}\)
Where,
\({y_{12}} = \frac{1}{{{z_{12}}}},;{y_{23}} = \frac{1}{{{z_{23}}}},;{y_{31}} = \frac{1}{{{z_{31}}}}\)
The above equation can be written in form of a matrix,
\(\left[ {\begin{array}{*{20}{c}} {{I_1}}\ {{I_2}}\ { - {I_3}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{y_{12}} + {y_{31}}}&{ - {y_{12}}}&{ - {y_{31}}}\ { - {y_{12}}}&{{y_{12}} + {y_{23}}}&{ - {y_{23}}}\ { - {y_{31}}}&{ - {y_{23}}}&{{y_{31}} + {y_{23}} + {y_{30}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{V_1}}\ {{V_2}}\ {{V_3}} \end{array}} \right]\)
\(\left[ {\begin{array}{*{20}{c}} {{I_1}}\ {{I_2}}\ { - {I_3}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{Y_{11}}}&{{Y_{12}}}&{{Y_{13}}}\ {{Y_{21}}}&{{Y_{22}}}&{{Y_{23}}}\ {{Y_{31}}}&{{Y_{32}}}&{{Y_{33}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{V_1}}\ {{V_2}}\ {{V_3}} \end{array}} \right]\)
Where,
Y11 = y12 + y31, Y22 = y21 + y23, Y33 = y30 + y31 + y32
Y12 = Y21 = - y12, Y23 = Y32 = -y23, Y31 = Y13 = -y31
The element Y11, Y22, Y33 is termed as self admittances.
The elements Y12, Y13, Y21, Y23, Y31, Y32 termed as mutual admittance.
Calculation:
Given circuit can be drawn as,
Applying KCL at node V,
I1 + I2 + I3 + I4 = 0
\(\frac{{{E_1} - V}}{{j1}} + \frac{{{E_2} - V}}{{j1}} + \frac{{{E_3} - V}}{{j1}} + \frac{{ - V}}{{j1}} = 0\)
4V = E1 + E2 + E3
\({I_1} = \frac{{{E_1} - V}}{{j1}} = \frac{{{E_1} - \left( {\frac{{{E_1} + {E_2} + {E_3}}}{4}} \right)}}{{j1}}\)
\({I_1} = j\frac{{ - 3}}{4}{E_1} + j\frac{1}{4}{E_2} + j\frac{1}{4}{E_3}\)
\({I_2} = \frac{{{E_2} - V}}{{j1}} = \frac{{{E_2} - \left( {\frac{{{E_1} + {E_2} + {E_3}}}{4}} \right)}}{{j1}}\)
\({I_2} = j\frac{1}{4}{E_1} + j\frac{{ - 3}}{4}{E_2} + j\frac{1}{4}{E_3}\)
\({I_3} = \frac{{{E_3} - V}}{{j1}} = \frac{{{E_3} - \left( {\frac{{{E_1} + {E_2} + {E_3}}}{4}} \right)}}{{j1}}\)
\({I_1} = j\frac{1}{4}{E_1} + j\frac{1}{4}{E_2} + j\frac{{ - 3}}{4}{E_3}\)
The equation of current can be written in form of the matrix,
\(\left[ {\begin{array}{*{20}{c}} {{I_1}}\ {{I_2}}\ {{I_3}} \end{array}} \right] = j\left[ {\begin{array}{*{20}{c}} {\frac{{ - 3}}{4}}&{\frac{1}{4}}&{\frac{1}{4}}\ {\frac{1}{4}}&{\frac{{ - 3}}{4}}&{\frac{1}{4}}\ {\frac{1}{4}}&{\frac{1}{4}}&{\frac{{ - 3}}{4}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{E_1}}\ {{E_2}}\ {{E_3}} \end{array}} \right]\)
\({Y_{bus}} = j\left[ {\begin{array}{*{20}{c}} {\frac{{ - 3}}{4}}&{\frac{1}{4}}&{\frac{1}{4}}\ {\frac{1}{4}}&{\frac{{ - 3}}{4}}&{\frac{1}{4}}\ {\frac{1}{4}}&{\frac{1}{4}}&{\frac{{ - 3}}{4}} \end{array}} \right]\)
\(Y_{bus}=\left[ \begin{array}{20 {c*}} {-\frac 3 4 j}&{\frac 1 4 j}&{\frac 1 4 j}\\ {\frac 1 4 j}&{-\frac 3 4 j}&{\frac 1 4j}\\ {\frac 1 4 j}&{\frac 1 4 j}&{-\frac 3 4 j} \end{array} \right]\)